Skip to content

Scattering, Bound States, and Tunneling

Scrollable formula. Use Left and Right arrow keys to reveal hidden content.

Scrollable table. Use Left and Right arrow keys to reveal additional columns.

Scrollable figure. Use Left and Right arrow keys to inspect hidden content.

Bound states, scattering states, virtual states, resonances, and tunneling are related but not interchangeable. Boundary conditions define the physical problem, conserved flux fixes probabilities, and analytic continuation reveals which poles belong to which class. This page develops those distinctions for short-range, single-particle potential scattering and then states the limited bridge to black-hole wave equations.

Required background. Calculus and Linear Algebra supplies complex amplitudes, inner products, and orthogonal mode expansions; Ordinary Differential Equations supplies second-order boundary-value problems, matching, and Wronskians; Quantum States, Observables, and Measurement supplies stationary states, self-adjoint Hamiltonians, and probability. Helpful background. Waves, Green Functions, and Geometric Optics supplies incoming and outgoing waves, Green functions, and barrier propagation.

Flux fixes reflection and transmission probabilities

Section titled “Flux fixes reflection and transmission probabilities”

Consider a particle of mass MM in a real one-dimensional potential,

−ℏ22Mψ′′+V(x)ψ=Eψ,-\frac{\hbar^2}{2M}\psi''+V(x)\psi=E\psi,

with time dependence e−iEt/ℏe^{-iEt/\hbar}. The probability current is

j=ℏMIm⁡(ψ∗ψ′).j=\frac{\hbar}{M}\operatorname{Im}(\psi^*\psi').

The stationary equation and its complex conjugate give j′=0j'=0. Suppose V(x)V(x) approaches constants VLV_L and VRV_R, and take EE above both asymptotic thresholds. Define

kL=2M(E−VL)ℏ,kR=2M(E−VR)ℏ.k_L=\frac{\sqrt{2M(E-V_L)}}{\hbar}, \qquad k_R=\frac{\sqrt{2M(E-V_R)}}{\hbar}.

A unit-amplitude wave incident from the left obeys

ψ(x)∼{eikLx+r(E)e−ikLx,x→−∞,t(E)eikRx,x→+∞.\psi(x)\sim \begin{cases} e^{ik_Lx}+r(E)e^{-ik_Lx},&x\to-\infty,\\ t(E)e^{ik_Rx},&x\to+\infty. \end{cases}

Since Ae±ikxAe^{\pm ikx} carries current j=±ℏk∣A∣2/Mj=\pm\hbar k|A|^2/M, the flux ratios are

R=∣r∣2,T=kRkL∣t∣2.R=|r|^2, \qquad T=\frac{k_R}{k_L}|t|^2.

Current conservation therefore gives

R+T=1.R+T=1.

The factor kR/kLk_R/k_L is essential when the asymptotic channels have different velocities. If instead E<VRE\lt V_R, the right-hand wavenumber is imaginary and that closed channel carries no asymptotic transmitted flux. These statements assume a real potential and no absorptive or additional inelastic channel. Tong derives the equal-threshold version from the conserved current in Tong 2017, §6.1.1, Eqs. (6.3)–(6.4), printed pp. 190–192.

Scattering states occur for a continuum of real energies. They are delta-normalized rather than square-integrable and contain a prescribed incoming component together with the allowed outgoing components.

Boundary conditions and sheets classify poles

Section titled “Boundary conditions and sheets classify poles”

For a short-range potential with continuum threshold E=0E=0, write E=ℏ2k2/(2M)E=\hbar^2k^2/(2M). The variable kk uniformizes the two-sheeted energy surface: the physical branch has k>0k\gt0 just above the continuum cut and Im⁡k>0\operatorname{Im}k\gt0 away from it, while continuation through the cut enters Im⁡k<0\operatorname{Im}k\lt0.

The pole-bearing object must also be named. For a self-adjoint Hamiltonian HH, the resolvent

R(E)=(H−E)−1\mathcal R(E)=(H-E)^{-1}

is initially defined away from the real spectrum. Meromorphic continuation requires decay and regularity hypotheses stronger than simply labeling a potential short range: compactly supported bounded potentials provide a standard sufficient class, while slower decay can restrict the continuation domain. In such a class, suitably localized resolvent matrix elements can be continued meromorphically in kk; the explicitly solvable point interaction below has the same continuation directly from its amplitudes. The one-channel scattering amplitudes exhibit the same relevant bound-state and resonance pole positions in these examples. The following table summarizes that setting. Threshold k=0k=0 is a special point and must be analyzed separately.

ClassComplex momentum or energyBoundary conditionHilbert-space status
Scattering statek>0k\gt0 realPrescribed incoming wave plus outgoing wavesDelta-normalized continuum state, not a pole
Bound statek=iκk=i\kappa, κ>0\kappa\gt0Exponential decay at both endsSquare-integrable discrete eigenstate; physical-sheet pole
Virtual statek=−iκk=-i\kappa, κ>0\kappa\gt0Exponential growth on the continued solutionNon-normalizable pole on the continued sheet
Decaying resonancek=k0−iγk=k_0-i\gamma, k0,γ>0k_0,\gamma\gt0Purely outgoing at both endsNon-normalizable continued-sheet pole

For a resonance it is often clearer to use

Epole=ER−iΓ2,Γ>0.E_{\rm pole}=E_R-\frac{i\Gamma}{2}, \qquad \Gamma\gt0.

Its time factor is

e−iEpolet/ℏ=e−iERt/ℏe−Γt/(2ℏ).e^{-iE_{\rm pole}t/\hbar} =e^{-iE_Rt/\hbar}e^{-\Gamma t/(2\hbar)}.

Thus the amplitude decays on scale 2ℏ/Γ2\hbar/\Gamma and the probability on the lifetime scale τ=ℏ/Γ\tau=\hbar/\Gamma in a regime dominated by this isolated pole. A normalizable wave packet also contains non-pole spectral contributions, so this exponential need not describe its very early or very late evolution. The purely outgoing spatial solution grows at infinity and is not an eigenvector in the self-adjoint Hamiltonian’s Hilbert space. It is useful because it controls the decay and the nearby real-energy line shape.

The sheet and boundary condition are part of the definition: a negative-energy pole is not automatically a bound state, and not every pole is a decaying resonance. As a rigorous model, Zworski treats PV=−Δ+VP_V=-\Delta+V with compactly supported VV on R3\mathbb R^3 and defines its resonances through meromorphic continuation of the resolvent in Zworski 2017, §2.2, Theorem 2 and Definition 1, printed p. 18. The associated outgoing resonant-state condition appears in Definition 2 and Eq. (2.24), printed p. 19. Other dimensions, long-range potentials, and black-hole radial operators require their own continuation theorems.

Partial waves make elastic unitarity visible

Section titled “Partial waves make elastic unitarity visible”

For a real, short-range central potential in three dimensions, choose an incident plane wave along the zz axis. The far field is

ψ(r)∼eikz+f(θ)eikrr,dσdΩ=∣f(θ)∣2.\psi(\mathbf r)\sim e^{ikz}+f(\theta)\frac{e^{ikr}}{r}, \qquad \frac{d\sigma}{d\Omega}=|f(\theta)|^2.

The angular-momentum page explains the rotation eigenfunctions. Here axial symmetry about the incident direction selects m=0m=0, for which Yℓ0Y_{\ell0} is proportional to the Legendre polynomial Pℓ(cos⁡θ)P_\ell(\cos\theta). Write the reduced radial solution as

uℓ(r)∼12i[Sℓei(kr−ℓπ/2)−e−i(kr−ℓπ/2)].u_\ell(r)\sim\frac{1}{2i} \left[ S_\ell e^{i(kr-\ell\pi/2)}-e^{-i(kr-\ell\pi/2)} \right].

The second term is incoming and the first is outgoing; this equation fixes their relative sign and defines the channel S-matrix SℓS_\ell. In a purely elastic channel, equal incoming and outgoing fluxes imply ∣Sℓ∣=1|S_\ell|=1, so one may write

Sℓ=e2iδℓS_\ell=e^{2i\delta_\ell}

with real phase shift δℓ\delta_\ell. The corresponding angular amplitude is

f(θ)=1k∑ℓ=0∞(2ℓ+1)eiδℓsin⁡δℓ Pℓ(cos⁡θ),f(\theta) =\frac1k\sum_{\ell=0}^{\infty}(2\ell+1) e^{i\delta_\ell}\sin\delta_\ell\, P_\ell(\cos\theta),

Orthogonality of the Legendre polynomials then gives

σtot=∑ℓ=0∞σℓ,σℓ=4πk2(2ℓ+1)sin⁡2δℓ.\sigma_{\rm tot} =\sum_{\ell=0}^{\infty}\sigma_\ell, \qquad \sigma_\ell =\frac{4\pi}{k^2}(2\ell+1)\sin^2\delta_\ell.

Each elastic channel obeys the unitarity bound

σℓ≤4πk2(2ℓ+1),\sigma_\ell\le\frac{4\pi}{k^2}(2\ell+1),

and the same expansion yields the optical theorem

σtot=4πkIm⁡f(0).\sigma_{\rm tot}=\frac{4\pi}{k}\operatorname{Im}f(0).

Tong derives the scattering amplitude, partial-wave expansion, optical theorem, and bound in Tong 2017, §§6.2.2–6.2.4, Eqs. (6.23)–(6.37), printed pp. 209–215. Long-range potentials such as the unscreened Coulomb potential require modified asymptotics and are outside this short-range formula.

Exact transfer check: one attractive delta well

Section titled “Exact transfer check: one attractive delta well”

Let

V(x)=−gδ(x),g>0,α=Mgℏ2.V(x)=-g\delta(x), \qquad g\gt0, \qquad \alpha=\frac{Mg}{\hbar^2}.

The wavefunction is continuous at the origin. Integrating the stationary equation across it gives

ψ′(0+)−ψ′(0−)=−2αψ(0).\psi'(0^+)-\psi'(0^-)=-2\alpha\psi(0).

For left incidence with the same wavenumber on both sides, continuity gives t=1+rt=1+r. Combining it with the derivative jump yields

t(k)=kk−iα,r(k)=iαk−iα.t(k)=\frac{k}{k-i\alpha}, \qquad r(k)=\frac{i\alpha}{k-i\alpha}.

For real kk,

∣t∣2=k2k2+α2,∣r∣2=α2k2+α2,|t|^2=\frac{k^2}{k^2+\alpha^2}, \qquad |r|^2=\frac{\alpha^2}{k^2+\alpha^2},

so the flux check gives R+T=1R+T=1. Both amplitudes have a pole at k=iαk=i\alpha. Directly imposing decay at both ends gives

ψb(x)=α e−α∣x∣,Eb=−ℏ2α22M=−Mg22ℏ2.\psi_b(x)=\sqrt{\alpha}\,e^{-\alpha|x|}, \qquad E_b=-\frac{\hbar^2\alpha^2}{2M} =-\frac{Mg^2}{2\hbar^2}.

The pole’s positive imaginary momentum and its decaying wavefunction identify a bound state, not a resonance. Tong derives the jump condition and bound state in Tong, Quantum Mechanics, undated, §2.3.2, Eqs. (2.39)–(2.40), printed pp. 48–49.

Tunneling is a mechanism, not a spectral class

Section titled “Tunneling is a mechanism, not a spectral class”

In a forbidden region, define

κ(x)=2M[V(x)−E]ℏ.\kappa(x)=\frac{\sqrt{2M[V(x)-E]}}{\hbar}.

For a smooth, semiclassically broad barrier between simple turning points x1x_1 and x2x_2, matching the WKB branches gives

T≈Cexp⁡ ⁣[−2∫x1x2κ(x) dx],T\approx C\exp\!\left[-2\int_{x_1}^{x_2}\kappa(x)\,dx\right],

where CC is a nonexponential connection prefactor in the controlled regime. The approximation-methods page derives the local WKB control and its turning-point failure. Tong obtains the tunneling exponent by connecting the decaying and growing branches in Tong 2017, §2.2.5, Eqs. (2.29)–(2.31), printed pp. 47–49.

An exact barrier calculation checks the exponent independently. Set the common exterior potential to zero and take a rectangular barrier of height V0V_0 and width LL, with 0<E<V00\lt E\lt V_0. Put

κ=2M(V0−E)ℏ.\kappa=\frac{\sqrt{2M(V_0-E)}}{\hbar}.

Direct matching gives

T=[1+V02sinh⁡2(κL)4E(V0−E)]−1.T=\left[ 1+\frac{V_0^2\sinh^2(\kappa L)}{4E(V_0-E)} \right]^{-1}.

In the broad-barrier limit κL≫1\kappa L\gg1,

T∼16E(V0−E)V02e−2κL.T\sim \frac{16E(V_0-E)}{V_0^2}e^{-2\kappa L}.

The prefactor depends on the sharp interfaces, but the exponential agrees with the WKB action because ∫0Lκ dx=κL\int_0^L\kappa\,dx=\kappa L. Tong derives the exact coefficient and this limit in Tong, Quantum Mechanics, undated, §2.4.2, Eqs. (2.51)–(2.54), printed pp. 60–62.

This estimate does not justify continuing an outer WKB branch through a turning point. Thin or discontinuous barriers require direct matching, while multiple barriers can produce interference and resonant transmission. A classically allowed well separated from infinity by a controlled barrier can produce a shape resonance whose small width is set by barrier leakage. Zworski analyzes this mechanism and its exponentially narrow resonances in Zworski 2017, §1.3, printed pp. 8–10.

Tunneling can therefore transmit a scattering wave, split bound levels, or give a quasibound state a finite width. It is not itself a fourth kind of state.

The bounded bridge to compact-object waves

Section titled “The bounded bridge to compact-object waves”

For a massless perturbation of a nonrotating, asymptotically flat black hole, separation often gives the schematic radial equation

d2udr∗2+[ω2−Veff(r)]u=0.\frac{d^2u}{dr_*^2}+[\omega^2-V_{\rm eff}(r)]u=0.

With time dependence e−iωte^{-i\omega t}, its quasinormal modes are ingoing at the future horizon, u∼e−iωr∗u\sim e^{-i\omega r_*} as r∗→−∞r_*\to-\infty, and outgoing at infinity, u∼e+iωr∗u\sim e^{+i\omega r_*} as r∗→+∞r_*\to+\infty. Its complex frequency is an outgoing resonance, not a normalizable bound-state energy. The terminology is geometric: both fluxes leave the exterior barrier region, one by falling through the future horizon and the other by escaping to infinity. These boundary conditions and the resulting lower-half-plane decay convention are stated in Berti, Cardoso, and Starinets 2009, §3.1, Eqs. (27) and (30), printed pp. 17–18. Rotating horizons replace the horizon wavenumber by a co-rotating combination, while de Sitter, anti-de Sitter, and massive-field problems impose different outer conditions.

A greybody problem instead prescribes incident data and compares Wronskian-normalized fluxes through the effective curvature barrier. At real frequency in a conservative scattering problem, flux conservation survives, but long-range tails modify the asymptotic states, rotation gives the horizon flux a signed co-rotating frequency factor, and coupled channels replace the scalar equation R+T=1R+T=1 by a summed or matrix relation. The canonical reflection, transmission, and absorption page derives those extensions. This page supplies the boundary-condition grammar, not a derivation of black-hole spectra.

Fixed-particle quantum mechanics also cannot create or annihilate quanta. Continue to Flat-Spacetime Quantum Field Theory for Gravitation for fields, Fock space, propagators, and interactions.

Using amplitudes as probabilities. Transmission is a ratio of fluxes. A wavenumber or group-velocity factor remains when the two open channels differ.

Classifying a pole without its sheet. The same real energy can sit beside different analytic continuations. Boundary conditions and sheet data separate bound, virtual, threshold, and resonance behavior.

Treating an outgoing resonance as a normalizable state. Its spatial growth is the analytic counterpart of temporal decay. It lies outside the self-adjoint Hamiltonian’s Hilbert-space spectrum.

Calling every small transmission a resonance. Generic tunneling can be weak without any nearby pole. Resonant transmission requires the additional phase and trapping conditions that create a quasibound mode.

Near an isolated elastic partial-wave resonance in an open channel with E0>0E_0\gt0, neglect the slowly varying background phase and take the Breit–Wigner form reviewed in Tong 2017, §6.2.7, Eq. (6.53), printed pp. 223–225:

Sℓ(E)=E−E0−iΓ/2E−E0+iΓ/2,Γ>0.S_\ell(E) =\frac{E-E_0-i\Gamma/2}{E-E_0+i\Gamma/2}, \qquad \Gamma\gt0.

Show that ∣Sℓ∣=1|S_\ell|=1 for real EE, locate its pole, determine the probability lifetime, and show that the partial cross-section saturates its unitarity bound at E=E0E=E_0.

Solution

For real EE, the numerator is the complex conjugate of the denominator, so ∣Sℓ∣=1|S_\ell|=1. The denominator vanishes at

Epole=E0−iΓ/2,E_{\rm pole}=E_0-i\Gamma/2,

which lies on the decaying continued sheet. Its time-dependent probability contains e−Γt/ℏe^{-\Gamma t/\hbar}, hence τ=ℏ/Γ\tau=\hbar/\Gamma.

At E=E0E=E_0,

Sℓ(E0)=−1=e2iδℓ,S_\ell(E_0)=-1=e^{2i\delta_\ell},

so δℓ=π/2\delta_\ell=\pi/2 modulo π\pi. Therefore

σℓ(E0)=4πk2(2ℓ+1),\sigma_\ell(E_0) =\frac{4\pi}{k^2}(2\ell+1),

where k=k(E0)k=k(E_0) is real. This is the elastic unitarity bound. A nonzero background phase or an inelastic channel changes the observed line shape and must be included in a less idealized model.

  • Berti, Emanuele, Vitor Cardoso, and Andrei O. Starinets. “Quasinormal Modes of Black Holes and Black Branes.” Classical and Quantum Gravity 26 (2009): 163001. doi:10.1088/0264-9381/26/16/163001. Open PDF.
  • Tong, David. Quantum Mechanics, §2.3.2, Eqs. (2.39)–(2.40), printed pp. 48–49; §2.4.2, Eqs. (2.51)–(2.54), printed pp. 60–62. Cambridge lecture notes, undated; official PDF inspected August 30, 2026. Official PDF.
  • Tong, David. Topics in Quantum Mechanics, §2.2.5, printed pp. 47–49; §6.1, printed pp. 189–205; §§6.2.2–6.2.4, printed pp. 209–215; §6.2.7, printed pp. 223–225. Cambridge lecture notes, 2017. Official PDF.
  • Zworski, Maciej. “Mathematical Study of Scattering Resonances.” Bulletin of Mathematical Sciences 7 (2017): 1–85. doi:10.1007/s13373-017-0099-4. Open PDF.