Papers
Topics
Authors
Recent
Search
2000 character limit reached

Double-Morse Potential: Dynamics & Quantum Features

Updated 12 November 2025
  • Double-Morse potential is a versatile, exactly-solvable model defined by the superposition of two Morse wells with controllable symmetry and barrier height.
  • It underpins studies in quantum state engineering, roaming dynamics, and optimal parameter estimation by revealing intricate phase-space and topological structures.
  • Its analytical solutions, including coherent states and nonclassicality metrics, provide practical insights for chemical dynamics and advanced quantum metrology.

The double-Morse potential is a paradigmatic model in theoretical physics and chemical dynamics, defined by the superposition of two Morse-type wells. Its symmetry, controllable shape, and exact solvability properties underpin its widespread use in studies of quantum state engineering, nonlinear dynamics, quantum metrology, and roaming reaction mechanisms. The model exhibits a rich topological and phase-space structure, including multiple potential wells, an index-one saddle, and an unbounded flat region enabling intricate transport phenomena.

1. Mathematical Formulation and Parameter Landscape

The one-dimensional double-Morse potential is constructed by centering two identical Morse wells at positions x=±x0x = \pm x_0, each of depth DD and width parameter α\alpha: VDM(x)=D[Acosh(αx)1]2,V_{\rm DM}(x) = D\,\bigl[\,A\,\cosh(\alpha x) - 1\,\bigr]^2, where A2eαx0A \equiv 2e^{-\alpha x_0}. For large αx0>ln2\alpha x_0 > \ln 2 (i.e., A<1A<1), this form describes a symmetric double-well structure of spacing 2x02x_0. The parameters have the following roles:

  • DD: well depth (dissociation energy)
  • α\alpha: width (asymmetry) parameter, controlling well separation and barrier height
  • x0x_0: half the separation between centers

In two dimensions, the planar double-Morse model generalizes to

VDM(x,y)=De[1eα(r+re)]2+De[1eα(rre)]2De,V_{\rm DM}(x, y) = D_e \bigl[1 - e^{-\alpha(r_+ - r_e)}\bigr]^2 + D_e \bigl[1 - e^{-\alpha(r_- - r_e)}\bigr]^2 - D_e,

where r±(x,y)=(xb)2+y2r_\pm(x, y) = \sqrt{(x\mp b)^2 + y^2}, rer_e is the equilibrium bond length, and bb sets the well separation. The resulting surface supports two equivalent minima at (x,yx, y) = (±b,0)(\pm b, 0), an index-one saddle at (0,0)(0,0), and a flat plateau as rr\to \infty (Carpenter et al., 2017).

2. Classical Dynamics and Roaming Mechanisms

Classical dynamics on the double-Morse surface are governed by a two-degree-of-freedom Hamiltonian,

H(x,y,px,py)=px2+py22m+VDM(x,y),H(x, y, p_x, p_y) = \frac{p_x^2 + p_y^2}{2m} + V_{\rm DM}(x, y),

with phase-space structures determined by unstable periodic orbits (POs), which serve as Normally Hyperbolic Invariant Manifolds (NHIMs) (Montoya et al., 2019, Carpenter et al., 2017). Three key PO families,

  • Type-1: in the far flat region (controls escape/dissociation)
  • Type-2: near the saddle (mediates direct interwell transfer)
  • Type-3: encircling each well (controls well-to-flat "roaming" transitions)

dictate the transport between wells and through the roaming region.

Roaming processes are identified by sequences of dividing surface crossings: an orbit originating in well A exits via the Type-3 DS into the flat region, wanders, and may enter well B or escape. In contrast, direct A→B transport crosses the Type-2 DS at the saddle without traversing the flat region (Carpenter et al., 2017). The boundaries between roaming, dissociation, and direct transfer are codified by the tube manifolds Ws,u(Γk)W^{s,u}(\Gamma_k) associated with each PO (Montoya et al., 2019).

Lagrangian descriptors, defined as arc-length integrals of trajectory velocities,

M+(x0,t0,τ+)=t0t0+τ+i=1nx˙i(t)2dt,M_+(\mathbf{x}_0, t_0, \tau_+) = \int_{t_0}^{t_0+\tau_+} \sum_{i=1}^n |\dot{x}_i(t)|^2\, dt,

reveal these manifolds and their intricate phase-space tangles without recourse to Poincaré maps, providing a diagnostic for roaming-dissociation mechanisms (Montoya et al., 2019).

3. Quantum Structure and Exact Solutions

The Schrödinger equation for a particle of mass mm in the double-Morse potential,

22mψ(x)+VDM(x)ψ(x)=Eψ(x),-\frac{\hbar^2}{2m} \psi''(x) + V_{\rm DM}(x)\psi(x) = E\psi(x),

admits exact, closed-form ground-state and energy expressions. With the change of variables y=αx/2y = \alpha x / 2, μ2=8mD/(2α2)\mu^2 = 8mD/(\hbar^2\alpha^2), and A=2eαx0A=2e^{-\alpha x_0}, the ground-state is

ψ0(x;α)=1K0(A)exp ⁣[A2cosh(αx)],\psi_0(x;\alpha) = \frac{1}{\sqrt{K_0(A)}} \exp\!\left[-\frac{A}{2} \cosh(\alpha x)\right],

with corresponding energy E0(α)=2α28m[1+A2]E_0(\alpha)=\frac{\hbar^2\alpha^2}{8m}[1 + A^2] (Chogle et al., 10 Nov 2025). Quasi-exact solvability occurs for μ=n+1\mu = n+1 (nNn\in\mathbb{N}), yielding n+1n+1 analytic energy levels.

The two-dimensional separable case supports a product of Morse eigenfunctions, with degeneracies controlled by the arithmetic nature of the spectral parameter pp (irrational pp eliminates all accidental degeneracies except for the two-fold symmetry nmn\leftrightarrow m) (Moran, 2021). Ladder operators and generalized coherent states can be constructed, reducing to Barut–Girardello-type states in suitable limits.

4. Nonclassicality and Quantum Metrology

The double-Morse potential is an analytically tractable, tunable source of non-Gaussianity and quantum nonclassicality. Two primary quantifiers are:

  • Non-Gaussianity ηNG\eta_{\rm NG}, given for pure states by

ηNG=h(detσ),\eta_{\rm NG} = h\left(\sqrt{\det \sigma}\right),

where h(u)=(u+12)ln(u+12)(u12)ln(u12)h(u) = (u + \frac{1}{2})\ln(u + \frac{1}{2}) - (u - \frac{1}{2})\ln(u - \frac{1}{2}) and σ\sigma is the covariance matrix. For the double-Morse ground state, ηNG(α)\eta_{\rm NG}(\alpha) increases strictly monotonically with α\alpha, from 0.0615\approx 0.0615 at the bifurcation threshold (αx0=ln2\alpha x_0 = \ln 2) to unbounded values as α\alpha\to\infty (Chogle et al., 10 Nov 2025).

  • Wigner negativity ηNC\eta_{\rm NC}, defined via the volume of the negative part of the Wigner function,

ηNC=νν+1,ν=W0(x,p)dxdp1,\eta_{\rm NC} = \frac{\nu}{\nu+1}, \quad \nu = \iint |W_0(x,p)|\,dx\,dp - 1,

also increases monotonically with α\alpha, approaching unity for strong well-separation.

The entanglement potential—bipartite entanglement generated by a 50:50 beam splitter with vacuum—also grows with α\alpha. These features position the double-Morse as a resource for continuous-variable quantum technologies. The analytic ground state serves as a testbed for studies in nonclassical state generation, continuous-variable quantum computation, and quantum error correction (Chogle et al., 10 Nov 2025).

In quantum metrology, estimation of the width parameter α\alpha via position measurements on the ground state reaches the quantum Cramér-Rao bound: Δα1NF(α),\Delta\alpha \geq \frac{1}{\sqrt{N\,\mathcal{F}(\alpha)}}, with

F(α)=4dx[αψ0(x;α)]2.\mathcal{F}(\alpha) = 4\int dx\, [\partial_\alpha \psi_0(x;\alpha)]^2.

The optimal measurement strategy requires only projective measurement of xx, with no advantage from nonclassical measurement protocols (Chogle et al., 10 Nov 2025).

5. Topology, Spectral Degeneracies, and Coherent States

Topological aspects of the double-Morse model manifest in both the classical and quantum regime. In two dimensions, energy levels are indexed by quantum numbers n,mn, m with energies

En,m=2β22m[(pn)2+(pm)2],E_{n,m} = -\frac{\hbar^2\beta^2}{2m} [(p-n)^2 + (p-m)^2],

where p=(ν1)/2p = (\nu-1)/2. For irrational pp, all accidental degeneracies are excluded except for the symmetry nmn\leftrightarrow m; for rational pp, further degeneracies remain due to sum-of-squares relations (Moran, 2021).

Generalized coherent states are constructed as eigenstates of lowering operators built from the nondegenerate, mixed states of the spectrum. The spatial probability distribution of these states exhibits sharp localization for Ψ1|\Psi| \lesssim 1 and spreads across several islands for larger Ψ|\Psi|, reflecting superpositions over higher eigenstates.

Uncertainty products (ΔQs)2(ΔPs)2(\Delta Q_s)^2(\Delta P_s)^2 for s=x,ys = x, y approach the minimum $1/4$ in units =1\hbar=1 for small Ψ|\Psi|, increasing with delocalization (Moran, 2021). These coherent states generalize the single-well Morse minimum-uncertainty states to multiwell settings, under precise degeneracy control.

6. Applications and Physical Significance

The double-Morse potential provides a versatile, exactly-solvable reference for a range of phenomena:

  • Chemical reaction dynamics and roaming: Models the phase-space topology underlying roaming processes in chemical reactions, encoding a flat “no-man’s-land” region through which trajectories can relocate between wells without additional saddles (Carpenter et al., 2017, Montoya et al., 2019).
  • Quantum state engineering: Acts as a tunable resource for preparing non-Gaussian, Wigner-negative states, exploitable in continuous-variable quantum computation, metrology, and quantum error correction (Chogle et al., 10 Nov 2025).
  • Metrology and parameter estimation: Enables optimal, high-precision estimation of system parameters such as α\alpha; high quantum Fisher information values in shallow double wells offer new pathways for field or asymmetry sensing.
  • Experimental realizability: Implementation proposals include cold atom traps in engineered optical potentials and superconducting circuit-QED systems with tailored Kerr or Morse-like nonlinearities.

A plausible implication is that as double-well separations and barrier heights become easily tunable (e.g., via α\alpha in optics or nanomechanics), the double-Morse model will see further application as both a benchmark and a resource for quantum-enhanced technologies.

7. Summary Table: Core Double-Morse Potential Properties

Property Formula / Value References
Potential (1D symmetric) VDM(x)=D[Acosh(αx)1]2V_{\rm DM}(x) = D\,[A\cosh(\alpha\,x)-1]^2 (Chogle et al., 10 Nov 2025)
Double-well separation control A=2eαx0A = 2e^{-\alpha x_0}, αx0>ln2    \alpha x_0>\ln 2 \implies double-well (Chogle et al., 10 Nov 2025)
2D potential minima At (x,±y)=(±b,0)(x,\pm y) = (\pm b,0), V=DeV=-D_e (Montoya et al., 2019)
Index-one saddle At (0,0)(0,0), V(0,0)=2De[1eα(bre)]2DeV(0,0) = 2D_e[1-e^{-\alpha(b - r_e)}]^2 - D_e (Carpenter et al., 2017)
Exact ground state (1D) ψ0(x;α)=(1/K0(A))exp[Acosh(αx)/2]\psi_0(x;\alpha) = (1/\sqrt{K_0(A)})\exp[-A\cosh(\alpha x)/2] (Chogle et al., 10 Nov 2025)
Non-Gaussianity ηNG(α)\eta_{\rm NG}(\alpha)\uparrow with α\alpha (Chogle et al., 10 Nov 2025)
Wigner negativity ηNC(α)\eta_{\rm NC}(\alpha)\uparrow with α\alpha (Chogle et al., 10 Nov 2025)
Roaming mechanism Trajectories cross Type-3 DS into CC, then to BB or escape (Carpenter et al., 2017)

The double-Morse potential thus constitutes a technically accessible but topologically rich model system, crucial to the modern theory of roaming transitions, quantum state engineering, and parameter estimation in quantum technologies.

Topic to Video (Beta)

No one has generated a video about this topic yet.

Whiteboard

No one has generated a whiteboard explanation for this topic yet.

Follow Topic

Get notified by email when new papers are published related to Double-Morse Potential.