Double-well potential Created 2026-10-25 Updated 2026-11-05
A double-well potential has two local minima separated by a potential energy barrier. A symmetric confining example is
Its minima are at and the central barrier has height . Near either minimum, writing gives , a harmonic potential with local angular frequency for mass .
For this symmetric one-dimensional bound-state problem, the lowest two energy eigenstates can be chosen as an even function and an odd function. They are delocalized across both wells, with splitting produced by quantum tunnelling. Define orthonormal combinations and . For well-separated wells with a sufficiently high barrier, these approximate states localized on the left and right. An initial evolves within this doublet with
The first full transfer takes . A two-mode approximation for a Bose gas uses analogous left/right orbitals; its validity also requires suppression of higher energy levels and appropriate control of interactions. Two local minima alone do not guarantee that every arbitrary nonconfining potential has two discrete bound levels.
In the classical thermal model, a coordinate in a harmonic potential obeys , where is the drag coefficient and . Its stationary one-sided spectral density per hertz is the displayed Lorentzian line shape. The integral over positive frequencies equals , consistent with the equipartition theorem. On a log-log plot the spectrum is flat for and has slope for . The crossover is determined by the restoring stiffness and drag; an added detector response multiplies the mechanical spectrum by the squared magnitude of its transfer function.
The WKB approximation describes a wave whose local wavelength changes slowly compared with one wavelength. For the one-dimensional stationary Schrödinger equation, put . In a classically allowed region,
The phase is the classical action divided by . The amplitude varies as so that a single travelling component has constant probability current. A useful local validity condition is .
For , write . The wave instead has growing and decaying components proportional to . At a simple turning point, vanishes and neither local WKB expression is valid. Linearizing the potential gives an Airy-function solution that connects the two sides and supplies the turning-point phase shift. It is not legitimate to carry the divergent WKB amplitude straight through the turning point.[ref]
A useful application is tunnelling through a smooth opaque barrier. With turning points , the leading transmission probability is
The exponent is twice the forbidden-region action in units of , because transmission is a probability rather than a wave amplitude. For an ideal rectangular barrier of height and width , the same leading exponential is ; matching at the abrupt edges supplies an additional prefactor.
Figure 1.
Original smooth-barrier construction and leading WKB transmission exponent; the shaded interval is classically forbidden
. Details.
For a bound state between two smooth turning points, the related WKB quantization rule is . In a harmonic potential the integral equals , recovering . This exact agreement of energies for the harmonic oscillator does not make WKB wavefunctions uniformly exact near its turning points.
Irrotational motion in this simply connected fluid region has a velocity potential, . Incompressibility gives , hence .
The moving surface is a material surface: differentiating along a particle gives at . Linearizing about a flat stationary surface drops the product term and evaluates at zero, giving there. The Unsteady Bernoulli equation, after dropping quadratic velocity terms and absorbing the time-only constant into , is . The linear capillary pressure at the surface is relative to atmosphere, in the sign convention of the source. Therefore
For choose the harmonic potential decaying at infinite depth,
Its vertical derivative at zero equals . Substitution into the dynamic condition gives . Thus
The two positive terms are the gravitational and capillary contributions to the restoring force.
Put , where . The Taylor expansion at a circular orbit is
The constant has no effect on the force. The paper's displayed harmonic potential should be read as this expansion after subtracting that constant and neglecting higher orders; it is not an equality for the absolute effective potential in general. For example, when and , and , whereas the printed quadratic expression is zero at .
Differentiating again and then using gives
The linearized radial equation of motion is , so the radial frequency in a power-law central potential is
For the circular orbit is unstable and radial disturbances grow exponentially, with growth rate . At the linear restoring force vanishes, as described in part (c). Thus an attractive central force alone does not guarantee small oscillations.
Ehrenfest theorem relates the dynamics of quantum expectation values to operator commutators. For a normalized state obeying and a possibly time-dependent observable , differentiate . Using the adjoint Schrödinger equation for the bra gives
Thus
The identity assumes a common domain on which the required products and derivatives are defined.
For , the canonical commutation relation gives and . Therefore
These resemble Hamilton's equations and Newton's second law, but generally . If the state is localized about , a local expansion gives
The classical trajectory equation for the mean is accurate when the force changes little across the wave packet. For a quadratic harmonic potential, is linear and this closure is exact for every state with the relevant finite moments: . This exact classical evolution of the centroid does not eliminate its quantum position uncertainty or make every state a point particle.
The correspondence principle is broader: quantum predictions recover classical ones in an appropriate limit, such as large action compared with and observations that cannot resolve individual quantum levels. It is not the assertion that every large quantum number automatically makes every possible state classical. Delocalized superpositions and interference can remain significant unless the physical preparation and measurement justify a classical description.
A spectral example is useful. For a quantum harmonic oscillator, neighboring levels differ by , giving transition frequency , equal to its classical oscillation frequency. A coherent state makes the connection more concrete: its mean position and momentum oscillate classically while the finite-width Gaussian profile retains its quantum uncertainty. For hydrogenic Bohr levels , transitions between nearby large have , which equals the classical circular orbital angular frequency in that model. Ehrenfest's theorem supplies an exact expectation-value bridge; a classical trajectory follows only when the force can also be evaluated at the mean position.
A Gaussian distribution with mean and standard deviation has probability density
Its symmetry, finite variance and simple behavior under linear combinations make it useful for repeated measurements. Roughly , and of observations lie within one, two and three true standard deviations of the mean. These are coverage probabilities for a specified distribution, not guarantees that every experimental error is normal.
Why it appears. The central limit theorem explains why many sums of small independent finite-variance contributions approach a Gaussian law after centering and scaling, provided no single contribution dominates. Detector noise, small mechanical disturbances and multiple electronic fluctuations can combine approximately this way. The theorem does not remove a common calibration bias, and it does not justify assuming normal tails for rare large disturbances or arbitrarily correlated errors. A large-mean Poisson distribution also approaches a Gaussian with variance equal to its mean, making normal approximations useful for high-count photon or particle data; low counts require the discrete model.
A classical particle at thermal equilibrium in a one-dimensional harmonic potential has a Gaussian position distribution directly from the Boltzmann distribution:
Thus a Gaussian law can arise from the energy dependence of equilibrium probability, rather than only from an informal appeal to many independent errors.
Precision of repeated measurements. For independent observations with the same mean and variance, the sample mean has variance . For Gaussian observations this mean is exactly Gaussian:
Consequently quadrupling the number of independent readings halves the random standard error. When the true variance is unknown, a small-sample confidence interval for a normal mean uses a Student t-distribution with statistical degrees of freedom; replacing that by a normal quantile can overstate precision. Shared systematic errors and correlations do not decrease by the same averaging rule.
Combining results and fitting models. Independent Gaussian estimates of one quantity, with known variances , combine through the inverse-variance weighted mean
For a model with independent zero-mean Gaussian noise, maximizing the likelihood is equivalent to minimizing the weighted least-squares statistic
For a correct linear model with known independent Gaussian errors and fitted parameters of full rank, the minimized statistic has a chi-squared distribution with degrees of freedom. Nonlinear fits only inherit that familiar approximation under additional regularity and adequate-data assumptions. Known correlated Gaussian errors require the full inverse covariance matrix in the quadratic form.
These properties support confidence intervals, consistency checks, calibration and error propagation. They also give practical diagnostics: inspect residuals for trends, changing spread, correlations or outliers before trusting a Gaussian fit. Gaussian models connect physical fluctuations to quantitative uncertainty, but their assumptions must be checked rather than inferred from the shape of a convenient bell curve.
Stokes theorem converts circulation around an oriented boundary into curl flux through a spanning surface. For a sufficiently smooth vector field on an oriented piecewise smooth surface ,
The boundary orientation follows the right-hand rule relative to . The field must be smooth on the surface, not merely on its boundary; a singularity crossing the surface cannot be ignored. For a fixed loop, applying the theorem to Faraday's differential law gives the integral flux law of induction.
The planar form is Green theorem. For a region with positively oriented boundary and sufficiently smooth ,
The outer boundary is counterclockwise and any holes have clockwise boundaries. For instance makes the area of equal to half the line integral of ; the unit circle gives area .
A closely related use of “Green's theorems” in field theory is Green's identities, obtained by applying the divergence theorem to :
and, after subtracting the expression with exchanged,
These relate interior sources to boundary data in electrostatics and other boundary-value problems. For example, a harmonic potential difference vanishing on a Dirichlet boundary condition has zero integral of its squared gradient by the first identity, establishing uniqueness. Both the planar circulation theorem and the field identities are included here to make the convention clear.
Take an overdamped particle with viscous drag coefficient , so its deterministic drift in the harmonic potential is . The one-dimensional Fokker-Planck equation has probability current
At thermal equilibrium, and the Boltzmann distribution is , with . Inserting this into the zero-current condition gives
The result is the Einstein diffusion relation between diffusion and mobility. It assumes a bath at one temperature, linear drag and thermal equilibrium. The same Gaussian probability density gives , providing an independent calibration of trap stiffness.