Skip to content

Double-Well Potential

A double-well potential is a one-dimensional potential with two separated minima. The symmetric double well is the canonical first encounter with tunneling splitting: low-energy states localized in the left and right wells mix through the barrier, producing even and odd energy eigenstates separated by a small energy difference.

A standard smooth example is

V(x)=λ(x2−a2)2,λ>0.V(x)=\lambda(x^2-a^2)^2, \qquad \lambda\gt0.

It has minima near x=±ax=\pm a and a barrier near x=0x=0. This page explains the qualitative and elementary two-state structure. Double-Well Tunneling develops the graduate physical picture, while Tunneling Splittings owns the quantitative Herring, WKB, and instanton comparison.

Double-Well Splitting is the graduate worked calculation for the quartic model, with a two-state reduction, WKB action, instanton saddle, and numerical benchmark.

The double well combines three ideas from earlier wave mechanics:

  • each individual well has approximately discrete bound states;
  • the barrier allows exponentially small leakage between wells;
  • parity symmetry organizes the exact eigenstates when V(x)=V(−x)V(x)=V(-x).

The result is not two independent copies of a single well. If the barrier is finite, the two wells communicate quantum mechanically.

The exactly solvable Double Delta Potential is the clean algebraic model of this same splitting mechanism.

Near either minimum, a smooth double well is approximately harmonic. For the quartic example,

V′′(a)=8λa2,V''(a)=8\lambda a^2,

so the local oscillator frequency is

ω0=8λa2m.\omega_0 =\sqrt{\frac{8\lambda a^2}{m}}.

If the barrier were impenetrable, the left and right local ground states would be separate states, denoted schematically by

∣L⟩,∣R⟩.\lvert L\rangle, \qquad \lvert R\rangle.

They would have nearly the same local energy E0E_0. For a finite barrier, these localized states are useful approximations, but they are not exact stationary states of the symmetric Hamiltonian.

When

V(x)=V(−x),V(x)=V(-x),

the Hamiltonian commutes with parity. Exact energy eigenstates can be chosen even or odd:

ψ+(x)=ψ+(−x),ψ−(x)=−ψ−(−x).\psi_+(x)=\psi_+(-x), \qquad \psi_-(x)=-\psi_-(-x).

For a high barrier, the lowest pair is approximately

∣+⟩=12(∣L⟩+∣R⟩),\lvert +\rangle =\frac{1}{\sqrt2} \left( \lvert L\rangle+\lvert R\rangle \right),

and

∣−⟩=12(∣L⟩−∣R⟩).\lvert -\rangle =\frac{1}{\sqrt2} \left( \lvert L\rangle-\lvert R\rangle \right).

The even state is usually lower in energy because it has no node between the wells, while the odd state has a node at the center.

In the localized basis, the lowest doublet is captured by

Heff=(E0−K−KE0),K>0.H_{\mathrm{eff}} = \begin{pmatrix} E_0 & -K \\ -K & E_0 \end{pmatrix}, \qquad K\gt0.

The eigenstates are the even and odd combinations above. The eigenvalues are

E+=E0−K,E−=E0+K.E_+ =E_0-K, \qquad E_- =E_0+K.

The splitting is

ΔE=E−−E+=2K.\Delta E =E_- - E_+ =2K.

The coupling KK is controlled by barrier penetration. It becomes smaller when the barrier is higher, wider, or when the mass is larger.

A localized state is a superposition of stationary states:

∣L⟩=12(∣+⟩+∣−⟩).\lvert L\rangle =\frac{1}{\sqrt2} \left( \lvert +\rangle+\lvert -\rangle \right).

Under time evolution,

∣ψ(t)⟩=e−i(E++E−)t/(2ℏ)[cos⁡(ΔE t2ℏ)∣L⟩+isin⁡(ΔE t2ℏ)∣R⟩],\lvert \psi(t)\rangle =e^{-i(E_++E_-)t/(2\hbar)} \left[ \cos\left(\frac{\Delta E\,t}{2\hbar}\right) \lvert L\rangle +i\sin\left(\frac{\Delta E\,t}{2\hbar}\right) \lvert R\rangle \right],

up to an overall phase convention for ∣R⟩\lvert R\rangle. Therefore

PR(t)=sin⁡2(ΔE t2ℏ).P_R(t) =\sin^2\left(\frac{\Delta E\,t}{2\hbar}\right).

The first complete transfer from left to right occurs at

ttransfer=πℏΔE.t_{\mathrm{transfer}} =\frac{\pi\hbar}{\Delta E}.

Small splitting means slow tunneling oscillation. This is why exponentially small energy differences can still have observable dynamical consequences.

At a rough semiclassical level, the splitting has the form

ΔE∝e−S/ℏ,\Delta E \propto e^{-S/\hbar},

where SS is an action through the classically forbidden barrier region. In a WKB estimate,

S∼∫x1x22m(V(x)−E0) dx,S \sim \int_{x_1}^{x_2} \sqrt{2m(V(x)-E_0)}\,dx,

with x1x_1 and x2x_2 the turning points inside the barrier. This formula is a scaling guide here, not a complete derivation. The derivation and prefactors are part of Barrier Penetration and Tunneling and the graduate double-well tunneling page.

Real double wells are often not exactly symmetric. A useful two-state Hamiltonian is

Heff=(EL−K−KER).H_{\mathrm{eff}} = \begin{pmatrix} E_L & -K \\ -K & E_R \end{pmatrix}.

The energy splitting becomes

Ehigh−Elow=2K2+(EL−ER2)2.E_{\mathrm{high}}-E_{\mathrm{low}} =2\sqrt{ K^2+\left(\frac{E_L-E_R}{2}\right)^2 }.

If ∣EL−ER∣≫K\lvert E_L-E_R\rvert\gg K, the eigenstates are mostly localized in one well. Small tunneling coupling cannot overcome a large energy bias. This is the basic reason why controllable two-state systems often tune a bias parameter through an avoided crossing.

Double-well language appears in many places:

  • ammonia inversion, where the nitrogen position relative to the hydrogen plane has two equivalent configurations;
  • tunneling splittings in molecular conformations;
  • flux and charge qubits, where two macroscopic circuit configurations can form an effective double well;
  • symmetric coupled wells in semiconductor and cold-atom settings.

These examples differ in microscopic detail. The shared structure is a pair of nearly degenerate localized states mixed by tunneling.

  • Treating ∣L⟩\lvert L\rangle and ∣R⟩\lvert R\rangle as exact stationary states in a symmetric finite double well.
  • Forgetting that exact symmetric-well eigenstates have definite parity.
  • Confusing the splitting ΔE\Delta E with the barrier height.
  • Assuming tunneling requires energy above the barrier.
  • Applying infinite-wall boundary conditions at the central barrier.
  • Ignoring asymmetry; even a small bias can suppress left-right oscillations if it is large compared with KK.
  • Using a WKB exponential estimate outside its semiclassical regime.
  • D. J. Griffiths and D. F. Schroeter, Introduction to Quantum Mechanics, 3rd ed., Cambridge University Press, 2018.
  • R. Shankar, Principles of Quantum Mechanics, 2nd ed., Springer, 1994.
  • C. Cohen-Tannoudji, B. Diu, and F. Laloë, Quantum Mechanics, Wiley, 1977.
  • L. D. Landau and E. M. Lifshitz, Quantum Mechanics: Non-Relativistic Theory, 3rd ed., Pergamon, 1977.
  1. Diagonalize the symmetric two-state Hamiltonian.
Solution

For

Heff=(E0−K−KE0),H_{\mathrm{eff}} = \begin{pmatrix} E_0 & -K \\ -K & E_0 \end{pmatrix},

the symmetric vector (1,1)/2(1,1)/\sqrt2 has eigenvalue

E+=E0−K,E_+=E_0-K,

and the antisymmetric vector (1,−1)/2(1,-1)/\sqrt2 has eigenvalue

E−=E0+K.E_-=E_0+K.

Thus the splitting is ΔE=2K\Delta E=2K.

  1. A doublet has splitting ΔE\Delta E. If the particle begins in ∣L⟩\lvert L\rangle, when is the first maximum probability to find it in ∣R⟩\lvert R\rangle?
Solution

The right-well probability is

PR(t)=sin⁡2(ΔE t2ℏ).P_R(t) =\sin^2\left(\frac{\Delta E\,t}{2\hbar}\right).

The first maximum occurs when

ΔE t2ℏ=π2.\frac{\Delta E\,t}{2\hbar}=\frac{\pi}{2}.

Therefore

t=πℏΔE.t=\frac{\pi\hbar}{\Delta E}.
  1. Explain qualitatively why increasing the particle mass reduces the tunneling splitting.
Solution

In the forbidden region, the WKB decay rate contains

2m(V(x)−E).\sqrt{2m(V(x)-E)}.

A larger mass increases the decay rate under the barrier, reducing the overlap between left and right localized states. Since the splitting is controlled by this overlap, ΔE\Delta E becomes smaller.