research34 min studyLesson 828 in Mathematics

Tunneling and Agmon estimates

Use weighted energy identities to control forbidden-region decay and distinguish decay estimates from a full tunneling-splitting asymptotic.

On this page

What you’ll understand

  • Derive the mass-one weighted Schrödinger identity and Agmon metric factor.
  • Compute exact attenuation through a constant forbidden region.
  • Explain the extra hypotheses needed to turn decay into a sharp double-well splitting formula.

A classically forbidden region does not force a quantum eigenfunction to vanish. Instead its amplitude can decay exponentially on the h scale. Agmon's method proves such decay through a weighted energy identity and an associated distance. Tunneling between wells uses that decay, but a matching lower bound and an asymptotic prefactor require additional analysis. An upper estimate on one tail is not by itself a complete formula for an energy splitting.

The weight competes with the potential barrier

For P_h=−h²Δ/2+V with real smooth V, take a compactly supported smooth test v and a real smooth weight φ. Integration by parts gives the identity below, with no boundary term because of compact support. The gradient of the weight subtracts |∇φ|²/2 from the effective potential. This factor one half follows from the mass-one kinetic normalization; using a distance formula from a source with kinetic term −h²Δ without rescaling would give the wrong decay exponent.

For a normalized eigenfunction, one applies the identity to a cutoff times that eigenfunction. The cutoff creates controlled source terms, so those terms must be retained rather than declaring that every weighted energy term vanishes. In a region where V−E−|∇φ|²/2 is positive, the identity controls the weighted amplitude and gradient there using contributions from the allowed region and the cutoff. This is an energy proof of decay, distinct from assuming a pointwise WKB ansatz is globally valid.

Re⁡⟨eϕ/h(Ph−E)v,eϕ/hv⟩=h22∥∇(eϕ/hv)∥22+∫(V−E−∣∇ϕ∣22)e2ϕ/h∣v∣2 dx\operatorname{Re}\langle e^{\phi/h}(P_h-E)v,e^{\phi/h}v\rangle=\frac{h^2}{2}\|\nabla(e^{\phi/h}v)\|_2^2+\int\left(V-E-\frac{|\nabla\phi|^2}{2}\right)e^{2\phi/h}|v|^2\,dx
The identity is for compactly supported smooth v and real φ. Cutoffs of eigenfunctions introduce nonzero commutator source terms.↔ Scroll wide equations horizontally.

Distance depends on forbidden action

For a fixed energy E, the Agmon length of a path integrates sqrt(2(V−E)_+) against ordinary arclength. Taking the infimum from the allowed set V≤E defines an Agmon distance. It satisfies the corresponding eikonal inequality almost everywhere, and a slightly reduced weight (1−ε) times that distance leaves a positive margin in the forbidden region. Smooth approximations and localization justify using this generally Lipschitz weight in the energy argument. Constants must be controlled over any energy interval being considered.

Under the standard confining self-adjoint Schrödinger hypotheses and an appropriate separated forbidden region, these weighted estimates yield exponential decay up to specified constants or polynomial h losses. The distance uses the relevant eigenenergy, which may depend on h. A bottom-energy action can instead describe a leading exponential scale for low-lying states, but it must not be mistaken for an exact zero-energy eigenfunction. The kinetic energy usually raises the true ground energy above the minimum of V.

dE(x,A)=inf⁡γ:A→x∫2(V(γ(s))−E)+ ∣γ˙(s)∣ dsd_E(x,A)=\inf_{\gamma:A\to x}\int\sqrt{2(V(\gamma(s))-E)_+}\,|\dot\gamma(s)|\,ds
A is the allowed set or a specified reference well set. The factor √2 matches −h²Δ/2; a reduced distance weight is used in rigorous estimates.↔ Scroll wide equations horizontally.

A splitting formula needs more than an upper tail

In a symmetric smooth double well with nondegenerate minima, low-lying states can be built near the individual wells and coupled through the barrier. Agmon estimates give small overlaps and useful upper bounds. A sharp leading splitting additionally needs a controlled interaction matrix, nonvanishing coupling and suitable barrier geometry; in higher dimensions multiple paths can complicate the coefficient. Symmetry, well equivalence and the chosen spectral pair must be stated before interpreting a tiny difference as tunneling.

It is therefore unsafe to write splitting equals exp(−distance/h) as a universal identity. The exponent may encode the leading action, while h powers and multiplicative constants depend on the well and barrier structure. A decay estimate alone supplies no matching lower bound and can miss cancellations. The elementary barrier calculation below deliberately concerns a local decaying solution, so its exact exponential attenuation is not advertised as a full transmission coefficient or a double-well spectral theorem.

Exact decay in a constant barrier

In an interval with V=2 and E=1/2, solve the decaying branch of (−h²u″/2+Vu)=Eu. Find the amplitude and squared-amplitude ratios across length one, including h=0.2.

  1. The equation reduces to u″=(3/h²)u. The right-decaying branch is u(x)=C exp(−√3 x/h); the growing branch is excluded by the chosen local matching condition.

  2. Across length one the amplitude ratio is exp(−√3/h), equal to exp(−5√3) when h=0.2. The Agmon action is exactly the constant integrand √(2(2−1/2)) times the length, hence √3.

  3. The squared-amplitude ratio is exp(−2√3/h), or exp(−10√3) at h=0.2. Halving h squares the previous amplitude ratio. Interface matching and incident flux have not been specified, so these are not asserted to be the full scattering transmission probability.

The result

Amplitude decay and squared-amplitude decay differ by a factor two in the exponent. This exact local model checks the kinetic normalization in the Agmon distance.

Explore the model

Interactive lab

Interactive lab · heat

Compare amplitude with its square

Pause and Reset, then set mode n = 1 and diffusivity κ = 1. After two Steps, t=1 and amplitude is exp(−1); its square is exp(−2). Reset and set both n = 1 and κ = 2, then take two Steps: amplitude is exp(−2), the square of the previous amplitude. This periodic heat model illustrates exponential attenuation arithmetic only. It has no barrier, eigenenergy, or h parameter and does not simulate quantum tunneling or determine a splitting prefactor.

Compare amplitude with its squareThe exact periodic heat-equation mode is u(x,t) = exp(−κn²t) sin(nx), on [0,2π]. Playback advances physical time. The initial sine is shown separately, and the unnormalized L² norm is √π exp(−κn²t). This analytic instrument does not run a finite-difference solver or test discretization stability. Current readings: Current amplitude: 1; Spatial L² norm: 1.772; Decay rate · κn²: 3.6 /s. All parameters are adjustable with the labeled controls below.Periodic heat equation mode 3 at time 0: the amplitude is 1.0π/2π3π/22π−1−0.500.51x · periodic intervalTemperature perturbation uExact mode · n = 3
Current temperature modeInitial sineCurrent amplitude
Current amplitude
1
Spatial L² norm
1.772
Decay rate · κn²
3.6 /s
t = 0 s

The exact periodic heat-equation mode is u(x,t) = exp(−κn²t) sin(nx), on [0,2π]. Playback advances physical time. The initial sine is shown separately, and the unnormalized L² norm is √π exp(−κn²t). This analytic instrument does not run a finite-difference solver or test discretization stability. Parameter changes pause playback and recompute the model at the current elapsed time. Reset restores both parameters and returns time to zero. Playback starts paused. Adjust either parameter, then choose Play or Step.

Look forAn attenuation factor for a wave amplitude is not automatically the same quantity as its squared intensity or a spectral splitting.

Check your understanding

What can an Agmon decay upper estimate alone establish about a double-well splitting?

Try it yourself

For V(x)=(x²−1)², compute the bottom-energy reference Agmon action between −1 and 1 at E=0 with the mass-one kinetic convention. Explain why this does not assert a zero-energy eigenstate or an exact splitting.

Take a moment to make a prediction or work it out on paper. Open a hint when you need one.

Give me a hint

Inside the interval, |x²−1|=1−x². Integrate √2(1−x²) and distinguish a reference action from the actual h-dependent eigenenergy.

Compare with the solution

The reference action is 2∫−11(1−x2) dx=42/3\sqrt2\int_{-1}^{1}(1-x^2)\,dx=4\sqrt2/3. It gives a natural leading forbidden-action scale for this symmetric quartic double well. The true low-lying eigenenergy is raised above zero by kinetic localization, so this calculation is not the distance at an exact E=0 eigenstate. Nor does it prove splitting equals exp(−4√2/(3h)): a spectral pair, nonvanishing coupling, h-dependent prefactor and error control must still be analyzed. The positive reference barrier action is one input to that analysis rather than the complete result.

Take it with you

  • Agmon weights turn forbidden-region decay into a controlled weighted energy estimate.
  • The metric factor follows the kinetic normalization and the relevant energy.
  • Exponential decay alone is not a sharp tunneling-splitting asymptotic or a transmission coefficient.

Read further

Original sources and deeper explanations to follow at your own pace.

Optional progress, saved in this browser when available. Revisit any lesson at any time.

Lesson not marked complete.