The Adiabatic Theorem

Grab the string of a swinging pendulum and yank it half as long: chaos. Shorten it gently over many swings and the motion stays "the same swing," adapting energy and frequency in lockstep — classical mechanics knew this as the adiabatic invariant. Quantum mechanics sharpens it into a theorem with teeth: change a Hamiltonian slowly compared with the clock set by its energy gaps, and a system starting in the nn-th eigenstate stays in the nn-th eigenstate of the moving Hamiltonian — even if, by the end, that Hamiltonian has been deformed beyond recognition. This is the opposite regime from perturbation theory: there the change is small but may be fast; here it may be enormous in total but must be gentle in rate. We build the machinery honestly — expand in the instantaneous eigenbasis, derive the exact coupled equations, and extract the identity ⟨m∣n˙⟩=⟨m∣H˙∣n⟩/(En−Em)\langle m|\dot n\rangle = \langle m|\dot H|n\rangle/(E_n - E_m) that powers the entire adiabatic family. A boundary-term argument then shows why slow driving cannot pump amplitude across a healthy gap. Then we do what makes the theorem believable: solve a live system exactly. A spin-½ riding a field that sweeps around a cone is one of the few genuinely time-dependent problems with a closed-form answer, and it exhibits adiabatic following quantitatively — the leak collapses as (sweep rate/Larmor rate)², confirmed by direct numerical integration of the Schrödinger equation to machine precision. The price of following is two phases, one dynamical and one geometric; the geometric one opens the door Geometric (Berry) Phase walks through next.

Key takeaways
  • If H(t)H(t) changes slowly compared with ℏ/(Em−En)\hbar/(E_m-E_n) gaps, a system prepared in ∣n(0)⟩|n(0)\rangle tracks the instantaneous eigenstate ∣n(t)⟩|n(t)\rangle, acquiring only phases.
  • The engine is exact: c˙m=−cm⟨m∣m˙⟩−∑n≠mcn ⟨m∣H˙∣n⟩En−Em ei(θn−θm)\dot c_m = -c_m\langle m|\dot m\rangle - \sum_{n\neq m} c_n\,\frac{\langle m|\dot H|n\rangle}{E_n-E_m}\,e^{i(\theta_n-\theta_m)}, with the off-diagonal identity derived from differentiating the eigenproblem.
  • The surviving phase splits into the dynamical θn=−1ℏ∫En dt\theta_n = -\frac{1}{\hbar}\int E_n\,\mathrm dt and a geometric piece γn=i∫⟨n∣n˙⟩ dt\gamma_n = i\int\langle n|\dot n\rangle\,\mathrm dt — the seed of the Berry phase.
  • Exact laboratory: a spin-½ in a field sweeping a cone at rate ω\omega leaks with maximum probability (ωsin⁡α/λ)2(\omega\sin\alpha/\lambda)^2 — 6.4×10−86.4\times10^{-8} at a 1 MHz sweep against a 2.8 GHz Larmor frequency, 0.1020.102 at 1 GHz.
  • Classical shadow: the quantum number nn is the quantum adiabatic invariant, heir to the pendulum's E/ωE/\omega.
  • Key idea: slowness is measured against the gap clock — fast internal phases average away the coupling between levels, locking the system to its moving eigenstate.
The core insight

Expand the state in the instantaneous eigenbasis and the slowness of H˙\dot H relative to the energy gaps suppresses every off-diagonal coupling: each eigenstate evolves independently, deforming with the Hamiltonian and collecting only a dynamical and a geometric phase.

Where this lives today

The adiabatic theorem is quietly load-bearing across physics. The Born-Oppenheimer approximation — all of quantum chemistry — is the theorem with nuclear coordinates as the slow parameters: electrons ride instantaneous eigenstates while nuclei crawl, and potential-energy surfaces are just En(R)E_n(\mathbf R) landscapes. Adiabatic-passage techniques in NMR and atomic physics sweep fields to transfer population with near-unit fidelity, robust against pulse imperfections because following is gap-protected. And adiabatic quantum computation stakes a computing paradigm on the theorem: encode the answer in a complicated Hamiltonian's ground state, start from an easy one, interpolate slowly — runtime governed by the minimum gap, which is why gap engineering, not qubit count, is the field's hard currency. Where protocols must outrun the theorem, "shortcuts to adiabaticity" add counterdiabatic terms that cancel the leak — engineering aimed squarely at the off-diagonal identity derived above.

The physical picture

Think of a compass needle in a field you are rotating. Turn the field slowly compared with the needle's precession and the needle glides along, always pointing with the field — it follows. Jerk the field and the needle is left behind, precessing about the new axis. The quantum spin-½ does exactly this, and the criterion is a race between two clocks: the external sweep clock (ω\omega) and the internal Larmor clock of the gap (ω0\omega_0). Adiabaticity is ω≪ω0\omega \ll \omega_0 — nothing about the total angle turned, everything about the rate.

Figure
Two clocks decide everything: slow sweeps are followed, fast sweeps leave the spin behind
adiabatic: ω ≪ ω₀ violated: ω ~ ω₀ B(t) ⟨S⟩ hugs the field cone swept at ω B(t) ⟨S⟩ left behind, precessing internal clock ω₀ = gap/ħ wins the race → the state deforms with H(t) external clock competes → amplitude leaks across the gap

Why does slow win? Because the coupling between levels is not just small — it is fast-oscillating. In the moving basis, the term feeding level mm from level nn beats at the gap frequency (Em−En)/ℏ(E_m-E_n)/\hbar; a slowly varying prefactor times a fast oscillation integrates to almost nothing, leaving boundary terms of order ℏH˙/(ΔE)2\hbar\dot H/(\Delta E)^2. The gap does double duty — it divides the coupling and drives the protective oscillation — which is why the criterion involves the gap squared, and why everything dies at a level crossing. The theorem is quantitative, and one system lets us watch it with no approximations. For the spin on the cone the exact leak probability is known in closed form, and its ceiling collapses like the square of the sweep rate — the curve below spans six decades, with both worked-example points on it.

Figure
The adiabatic theorem measured: maximum leak probability vs sweep rate for a spin-½ on a 45° cone, Larmor 2.8 GHz
max P(leak) = (ω sinα/λ)², exact — log-log, slope 2 in the adiabatic corner 10⁶ 10⁷ 10⁸ 10⁹ 10¹⁰ sweep frequency ω/2π (Hz) 10⁻⁸ 10⁻² 1 1 MHz sweep: P = 6.4×10⁻⁸ — following is essentially perfect 1 GHz sweep: P = 0.102 — adiabaticity violated ω = ω₀: the two clocks tie — the curve bends, slope-2 law ends slope 2: halve the rate, quarter the leak curve and both markers verified by direct TDSE integration (agreement 5×10⁻¹⁵) in the verification script

The same physics runs in a box. Widen an infinite well slowly and a ground-state particle stays in the instantaneous ground state, its energy sliding down E1(t)=π2ℏ2/2mL(t)2E_1(t) = \pi^2\hbar^2/2mL(t)^2 as the walls recede — not lost, but paid out as work by quantum pressure on the retreating wall. Yank the wall instead and the wavefunction has no time to adapt: it stays frozen and must be re-expanded in the new box — the sudden regime, owned by The Sudden Approximation later in this chapter.

Figure
The quantum piston: a slowly widening box keeps its particle in the ground state while quantum pressure pays out work
slow wall: n = 1 is preserved, E₁(L) = π²ħ²/2mL² tracks down ψ₁ at L = 1 nm slow still ψ₁ — now of the 2 nm box wall does not outrun the internal clock ħ/(E₂−E₁) energy ledger: E₁(1 nm) = 0.376 eV → E₁(2 nm) = 0.094 eV; the missing 0.282 eV = work done on the wall by quantum pressure F = 2E₁/L

The moving eigenbasis: exact machinery

Freeze the Hamiltonian at each instant and diagonalize it. This defines the instantaneous eigenproblem — a family of snapshots, with no dynamics in it yet:

H(t) ∣n(t)⟩=En(t) ∣n(t)⟩(1)H(t)\,|n(t)\rangle = E_n(t)\,|n(t)\rangle \qquad (1)

where: ∣n(t)⟩|n(t)\rangle — the nn-th eigenstate of the frozen H(t)H(t), assumed non-degenerate; En(t)E_n(t) — its eigenvalue. These are not solutions of the Schrödinger equation; they are a moving coordinate system the true solution will cling to in the slow regime. Expand the true state over this moving basis, peeling off the phase each snapshot would accumulate if it were genuinely stationary — the choice that cancels the EnE_n term identically:

∣ψ(t)⟩=∑ncn(t) eiθn(t) ∣n(t)⟩,θn(t)=−1ℏ∫0tEn(t′) dt′(2)|\psi(t)\rangle = \sum_n c_n(t)\, e^{i\theta_n(t)}\, |n(t)\rangle, \qquad \theta_n(t) = -\frac{1}{\hbar}\int_0^t E_n(t')\,\mathrm dt' \qquad (2)

where: cn(t)c_n(t) — the amplitude actually carrying the dynamics; θn\theta_n — the dynamical phase, the running generalization of −Ent/ℏ-E_nt/\hbar. Were HH constant, every cnc_n would freeze — the statistics-frozen stationary states of Stationary States and Energy Eigenstates; everything beyond that lives in the motion of the basis itself.

Derivation AThe exact equations of motion in a moving eigenbasis

We want the exact — no approximation yet — equations for the amplitudes cm(t)c_m(t), and from them the two structural results powering the whole adiabatic family: the off-diagonal identity ⟨m∣n˙⟩=⟨m∣H˙∣n⟩/(En−Em)\langle m|\dot n\rangle = \langle m|\dot H|n\rangle/(E_n-E_m), and the boundary-term argument showing why a slow sweep across a healthy gap transfers almost nothing. The strategy is substitution of expansion (2) into the Schrödinger equation, then projection; the trap is forgetting that the basis kets themselves move, so d/dt\mathrm d/\mathrm dt hits three factors, not two.

Step 1. Insert (2) into iℏ ∂t∣ψ⟩=H∣ψ⟩i\hbar\,\partial_t|\psi\rangle = H|\psi\rangle. The derivative produces three terms per summand — c˙n\dot c_n, the phase derivative, and the moving ket ∣n˙⟩|\dot n\rangle. The phase was engineered so its derivative cancels the right-hand side exactly, leaving

∑neiθn(c˙n ∣n⟩+cn ∣n˙⟩)=0.\sum_n e^{i\theta_n}\left(\dot c_n\,|n\rangle + c_n\,|\dot n\rangle\right) = 0.

Why the cancellation is exact. θ˙n=−En/ℏ\dot\theta_n = -E_n/\hbar by construction — the dynamical phase absorbs all the "trivial" evolution, isolating what the moving basis does.

Step 2. Project with ⟨m∣\langle m| and use orthonormality:

c˙m=− cm ⟨m∣m˙⟩  −  ∑n≠mcn ⟨m∣n˙⟩ ei(θn−θm).\dot c_m = -\,c_m\,\langle m|\dot m\rangle \;-\; \sum_{n\neq m} c_n\,\langle m|\dot n\rangle\, e^{i(\theta_n - \theta_m)}.

Exact, for any speed. The diagonal term will become a phase; the off-diagonal terms are the leak channels.

Step 3. Convert ⟨m∣n˙⟩\langle m|\dot n\rangle into Hamiltonian language. Differentiate the eigenproblem H∣n⟩=En∣n⟩H|n\rangle = E_n|n\rangle in time and project on ⟨m∣\langle m| with m≠nm \neq n: the E˙n\dot E_n term drops (orthogonality), and ⟨m∣H∣n˙⟩=Em⟨m∣n˙⟩\langle m|H|\dot n\rangle = E_m\langle m|\dot n\rangle pulls out, giving

⟨m∣n˙⟩=⟨m∣H˙∣n⟩En−Em(m≠n).\langle m|\dot n\rangle = \frac{\langle m|\dot H|n\rangle}{E_n - E_m} \qquad (m\neq n).

Why this identity is the engine. It converts "how fast the basis turns" into "how fast the Hamiltonian changes, divided by the gap" — the only two quantities a validity criterion could ever involve. On the two-level family H=ℏω02(cos⁡θ σz+sin⁡θ σx)H = \frac{\hbar\omega_0}{2}(\cos\theta\,\sigma_z + \sin\theta\,\sigma_x) both sides evaluate to 12\frac{1}{2} per unit θ\theta — checked symbolically in the verification script.

Step 4. Let the drive be slow, starting in cn(0)=1c_n(0)=1. To first order,

cm(T)≈−∫0T⟨m∣H˙∣n⟩En−Em  eiℏ∫0t(Em−En) dt′ dt,c_m(T) \approx -\int_0^T \frac{\langle m|\dot H|n\rangle}{E_n-E_m}\; e^{\frac{i}{\hbar}\int_0^t (E_m-E_n)\,\mathrm dt'}\,\mathrm dt,

a slowly varying prefactor times a phase beating at the gap frequency. Integrate by parts: each pass pulls out ℏ/(Em−En)\hbar/(E_m-E_n) and leaves boundary terms, so

∣cm∣∼ℏ ∣⟨m∣H˙∣n⟩∣(Em−En)2≪1|c_m| \sim \frac{\hbar\,|\langle m|\dot H|n\rangle|}{(E_m-E_n)^2} \ll 1

for smooth, slow drives. The gap appears squared — once as the identity's denominator, once from the oscillation the integration by parts exploits. Why this is a sketch and not a proof. The boundary-term estimate assumes the prefactor varies slowly and smoothly; making "slowly" precise, quantifying the error, and cataloguing the ways the naive criterion can fail is the business of Adiabatic Approximation and Validity Conditions, which owns the dimensionless criterion built from exactly this ratio. Rigorous proofs are the Born-Fock (1928) and Kato (1950) results cited in the history band. Verified: the 2×2 symbolic check of the off-diagonal identity, the eigenvalue bookkeeping, and the vanishing of the diagonal connection for the real two-level family are in verify_aq-qm-adiabatic-theorem.py (22 checks green).

The phases that survive

Suppose the leak terms are negligible. The exact equation of Step 2 collapses to c˙n=−cn⟨n∣n˙⟩\dot c_n = -c_n\langle n|\dot n\rangle, which integrates to a pure phase — ⟨n∣n˙⟩\langle n|\dot n\rangle is purely imaginary since ∂t⟨n∣n⟩=0\partial_t\langle n|n\rangle = 0. The adiabatic solution is

∣ψ(t)⟩≈eiθn(t) eiγn(t) ∣n(t)⟩,γn(t)=i∫0t⟨n(t′)∣n˙(t′)⟩ dt′(3)|\psi(t)\rangle \approx e^{i\theta_n(t)}\, e^{i\gamma_n(t)}\, |n(t)\rangle, \qquad \gamma_n(t) = i\int_0^t \langle n(t')|\dot n(t')\rangle\,\mathrm dt' \qquad (3)

where: θn\theta_n — the dynamical phase of (2), clocking the running energy; γn\gamma_n — the geometric phase, fed not by energy but by the motion of the eigenvector itself. For real eigenvector families γn\gamma_n vanishes identically (checked symbolically), which is why it went unnoticed for decades. When the parameter path closes into a loop, γn\gamma_n becomes gauge-invariant and measurable — that loop physics, the spin-½ solid-angle result, and the curvature picture are owned by Geometric (Berry) Phase; here we establish only that the phase exists and where it comes from. "Stays in the eigenstate" thus has fine print: the system stays in the eigenstate ray, and the two phases riding on it are physics — interference experiments read them.

An exact laboratory: the spin on a cone

Claims about slowness deserve a system where nothing is approximated. Take an electron spin (Larmor frequency ω0=γB\omega_0 = \gamma B) in a field of fixed magnitude whose direction sweeps a cone of half-angle α\alpha at rate ω\omega:

H(t)=ℏω02 n^(t)⋅σ,n^=(sin⁡αcos⁡ωt,  sin⁡αsin⁡ωt,  cos⁡α)(4)H(t) = \frac{\hbar\omega_0}{2}\,\hat{\mathbf n}(t)\cdot\boldsymbol\sigma, \qquad \hat{\mathbf n} = (\sin\alpha\cos\omega t,\; \sin\alpha\sin\omega t,\; \cos\alpha) \qquad (4)

where: ω0\omega_0 — the internal gap clock (E+−E−=ℏω0E_+ - E_- = \hbar\omega_0); ω\omega — the external sweep clock. The instantaneous eigenstates are the spin-along-n^\hat{\mathbf n} spinors (cos⁡α2, eiωtsin⁡α2)(\cos\frac{\alpha}{2},\, e^{i\omega t}\sin\frac{\alpha}{2}) — verified as exact eigenvectors in the script. This is the compass needle of the opening picture, exactly solvable because the time dependence is a rigid rotation.

Go deeper · Derivation B · 2 minAdiabatic following in a rotating field, solved exactly
Derivation BAdiabatic following in a rotating field, solved exactly

We want the exact probability that the spin, prepared along the initial field direction, is later found against the instantaneous field — the leak the adiabatic theorem says should vanish for ω≪ω0\omega \ll \omega_0. The strategy is the rotating-frame transformation: co-rotating with the field, the Hamiltonian becomes time-independent, the problem reduces to precession about a fixed effective field, and the exact answer drops out. The move to watch: changing frame adds a term to the Hamiltonian — the frame's own rotation — which is precisely what tilts the effective field away from n^\hat{\mathbf n} and causes the leak.

Step 1. Pass to the frame rotating at ω\omega about zz: ∣ϕ⟩=eiωt σz/2∣ψ⟩|\phi\rangle = e^{i\omega t\,\sigma_z/2}|\psi\rangle. The transformed Hamiltonian is time-independent (computed symbolically in the script):

Hrot=ℏ2[ ω0sin⁡α σx+(ω0cos⁡α−ω) σz].H_{\text{rot}} = \frac{\hbar}{2}\Big[\,\omega_0\sin\alpha\,\sigma_x + \big(\omega_0\cos\alpha - \omega\big)\,\sigma_z\Big].

Why the −ω-\omega appears. The unitary is time-dependent, so H→UHU†+iℏ U˙U†H \to UHU^\dagger + i\hbar\,\dot U U^\dagger; the second piece is the fictitious field of the rotating frame — the quantum version of the centrifugal bookkeeping.

Step 2. Read off the effective precession vector: magnitude

λ=ω02+ω2−2 ω0 ωcos⁡α,\lambda = \sqrt{\omega_0^2 + \omega^2 - 2\,\omega_0\,\omega\cos\alpha},

pointing along m^∝(ω0sin⁡α, 0, ω0cos⁡α−ω)\hat{\mathbf m} \propto (\omega_0\sin\alpha,\, 0,\, \omega_0\cos\alpha - \omega). In this frame the field direction n^(0)\hat{\mathbf n}(0) is fixed, so "instantaneous eigenstate" means the static spinor along n^(0)\hat{\mathbf n}(0), while the state precesses about m^\hat{\mathbf m} at rate λ\lambda.

Step 3. The leak is now a precession problem: starting along n^\hat{\mathbf n}, the probability of being found along −n^-\hat{\mathbf n} at time tt is sin⁡2β sin⁡2(λt/2)\sin^2\beta\,\sin^2(\lambda t/2), with β\beta the angle between n^\hat{\mathbf n} and m^\hat{\mathbf m}. The cross product gives sin⁡β=ωsin⁡α/λ\sin\beta = \omega\sin\alpha/\lambda (verified symbolically), so

Pleak(t)=(ωsin⁡αλ)2sin⁡2 ⁣(λt2).P_{\text{leak}}(t) = \left(\frac{\omega\sin\alpha}{\lambda}\right)^{2}\sin^2\!\Big(\frac{\lambda t}{2}\Big).

Exact at every sweep rate — the general driven-two-level flopping structure, whose resonant anatomy belongs to Transition Probability for a Periodic Two-Level System; we need only its adiabatic corner.

Step 4. Read the corner. For ω≪ω0\omega \ll \omega_0: λ→ω0\lambda \to \omega_0 and the ceiling collapses quadratically,

Pmax=(ωsin⁡αλ)2≈(ωω0)2sin⁡2α  ⟶  0.P_{\text{max}} = \left(\frac{\omega\sin\alpha}{\lambda}\right)^2 \approx \left(\frac{\omega}{\omega_0}\right)^2\sin^2\alpha \;\longrightarrow\; 0.

The spin follows the field with fidelity 1−O(ω/ω0)21 - O(\omega/\omega_0)^2 — the adiabatic theorem, exact and quantitative. At the other extreme ω≫ω0\omega \gg \omega_0, λ→ω\lambda \to \omega and Pmax→sin⁡2αP_{\text{max}} \to \sin^2\alpha: the spin stands still while the field runs away — the sudden regime. Verified: the rotating-frame Hamiltonian, λ, the cross-product amplitude, and the full P(t) are checked symbolically; direct lab-frame TDSE integration reproduces the exact formula to 5×10⁻¹⁵ at both 1 GHz and 1 MHz sweeps, and the log-log slope of the ceiling in the adiabatic corner measures 2.002 (verification script).

The classical bridge and the quantum piston

The theorem has a classical ancestor, and honoring it sharpens the quantum statement. The mechanics track's Adiabatic Invariants shows a pendulum with a slowly shortened string conserving not its energy but J=E/ωJ = E/\omega, the energy growing as E∝ℓ−1/2E \propto \ell^{-1/2} and the amplitude as θmax⁡∝ℓ−3/4\theta_{\max} \propto \ell^{-3/4}; its Adiabatic Theorem for Slowly Varying Hamiltonians is honest about where invariance must fail — a pendulum pumped toward the top of its arc violates adiabaticity at the swing-to-spin crossover no matter how slow the pump. The quantum heir of the invariant action is the quantum number: since En/ω∝(n+12)ℏE_n/\omega \propto (n+\frac12)\hbar for an oscillator, "conserve E/ωE/\omega" and "stay in level nn" are one sentence in two languages — and the classical failure mode translates too: a diverging period is a collapsing gap. The infinite square well makes the energy bookkeeping vivid. Stay in n=1n=1 while the wall recedes from LL to L′L' and the energy must slide from E1(L)E_1(L) to E1(L′)E_1(L') — where does the difference go? Into mechanical work: the confined state pushes on the wall with the quantum pressure force F=−dE1/dL=2E1/LF = -\mathrm dE_1/\mathrm dL = 2E_1/L (derived symbolically in the script), and the wall, retreating slowly, is paid exactly the energy the level loses. For an electron in a 1 nm box, E1=0.376E_1 = 0.376 eV and F=1.2×10−10F = 1.2\times10^{-10} N; widening to 2 nm delivers 0.282 eV to the wall. Adiabatic following conserves the quantum number, never the energy.

Limiting cases

caseresultmeaning
slow, smooth drive; healthy gapstate tracks ∣n(t)⟩\vert n(t)\rangle, leak ∼[ℏH˙mn/ΔE2]2\sim [\hbar\dot H_{mn}/\Delta E^2]^2the adiabatic theorem; criterion quantified in Adiabatic Approximation and Validity Conditions
instantaneous changestate frozen; re-expand in new eigenbasisthe sudden regime — opposite limit, owned by The Sudden Approximation
ω≪ω0\omega \ll \omega_0 (spin on cone)Pleak≤(ω/ω0)2sin⁡2αP_{\text{leak}} \le (\omega/\omega_0)^2\sin^2\alphaquadratic collapse of the leak — the measured slope-2 law
ω≫ω0\omega \gg \omega_0 (spin on cone)Pmax→sin⁡2αP_{\text{max}} \to \sin^2\alphaspin cannot follow; sudden limit recovered inside the exact solution
gap closes mid-sweepcriterion fails at any finite rateavoided-crossing physics — Landau–Zener Transitions owns the sweep-through problem
degenerate levelphases become a matrixnon-Abelian (Wilczek-Zee) transport — one line here, graduate-track territory
HH returns to its start (γn\gamma_n over a loop)gauge-invariant geometric phaseowned by Geometric (Berry) Phase

Worked examples

Example 1 — how slow is slow: the cone numbers. An electron spin sits in a B=0.1B = 0.1 T field (ω0/2π=28.02 GHz/T×0.1 T=2.80\omega_0/2\pi = 28.02\ \text{GHz/T} \times 0.1\ \text{T} = 2.80 GHz) whose direction sweeps a cone of half-angle α=45°\alpha = 45°. Compute the worst-case leak probability for sweep rates ω/2π=1\omega/2\pi = 1 MHz and 1 GHz.

Step 1. Form the exact ceiling Pmax=(ωsin⁡α/λ)2P_{\text{max}} = (\omega\sin\alpha/\lambda)^2 with λ=ω02+ω2−2ω0ωcos⁡α\lambda = \sqrt{\omega_0^2 + \omega^2 - 2\omega_0\omega\cos\alpha}.

Step 2. Slow sweep, ω/2π=1\omega/2\pi = 1 MHz: ω/ω0=3.57×10−4\omega/\omega_0 = 3.57\times10^{-4}, λ≈ω0\lambda \approx \omega_0, so Pmax≈(3.57×10−4×0.707)2=6.4×10−8P_{\text{max}} \approx (3.57\times10^{-4}\times0.707)^2 = \textbf{6.4}\times\textbf{10}^{-8} — following is near-perfect through every turn.

Step 3. Fast sweep, ω/2π=1\omega/2\pi = 1 GHz: ω/ω0=0.357\omega/\omega_0 = 0.357, λ/ω0=0.789\lambda/\omega_0 = 0.789, so Pmax=(0.357×0.707/0.789)2=0.102P_{\text{max}} = (0.357\times0.707/0.789)^2 = \textbf{0.102} — a ten-percent leak; the intermediate ladder runs 6.4×10−66.4\times10^{-6} at 10 MHz, 6.7×10−46.7\times10^{-4} at 100 MHz, riding the slope-2 law.

Step 4. Cross-check: direct numerical integration of the lab-frame Schrödinger equation reproduces the exact P(t)P(t) at both rates to 5×10−155\times10^{-15} — the "theorem" here is not asymptotic hand-waving but a measured property of the exact solution.

Example 2 — the quantum piston's energy ledger. An electron occupies the ground state of a 1 nm infinite well whose right wall recedes slowly to 2 nm. Track the energy.

Step 1. Initial energy: E1(1 nm)=π2ℏ2/2mL2=0.376 eVE_1(1\ \text{nm}) = \pi^2\hbar^2/2mL^2 = \textbf{0.376 eV}.

Step 2. Adiabatic following keeps n=1n = 1, so the final energy is E1(2 nm)=0.376/4=0.094 eVE_1(2\ \text{nm}) = 0.376/4 = \textbf{0.094 eV}.

Step 3. The difference, 0.282 eV\textbf{0.282 eV}, is work done on the wall: the level exerts F=2E1/LF = 2E_1/L (1.2×10−101.2\times10^{-10} N at 1 nm — a tenth of a nanonewton of pure confinement pressure), and ∫F dL\int F\,\mathrm dL reproduces the energy drop exactly, since F=−dE1/dLF = -\mathrm dE_1/\mathrm dL.

Step 4. Contrast the yank: pull the wall instantly and the energy expectation is unchanged (no time to do work) — the old ψ1\psi_1 is now a superposition over the new box's levels. Which superposition is The Sudden Approximation's story.

Origins · 1 minAdiabatic Theorem and how this was discovered
AT
Adiabatic Theorem
1916–1950 · Leiden, Göttingen, Tokyo

The adiabatic theorem predates the mechanics it now lives in. In 1916 Paul Ehrenfest, hunting for classical quantities that could survive quantization, proposed the adiabatic hypothesis: slow deformations carry allowed motions into allowed motions, and the invariant actions E/ωE/\omega are exactly the quantities to quantize — a guiding principle of the old quantum theory. When wave mechanics arrived, Max Born and Vladimir Fock recast the hypothesis as a theorem about the Schrödinger equation (1928), proving eigenstate-following for non-degenerate spectra. The modern rigorous form is Tosio Kato's 1950 proof, which relaxed the hypotheses and introduced the gap condition in the operator language still used in quantum-annealing runtime bounds. One theorem, three eras: quantization heuristic, wave-mechanics result, functional-analytic workhorse.

History
1916
— Ehrenfest states the adiabatic hypothesis: slow change maps allowed motions to allowed motions, singling out adiabatic invariants for quantization
1928
— Born and Fock prove the adiabatic theorem for the Schrödinger equation with non-degenerate, non-crossing spectra
1932
— Landau, Zener, Stückelberg and Majorana independently solve the sweep through an avoided crossing, quantifying where adiabaticity fails
1950
— Kato publishes the rigorous adiabatic theorem with the gap condition, in the operator language of modern proofs
1984
— Berry isolates the geometric phase hiding in adiabatic transport, turning the theorem's fine print into a field
2011
— commercial quantum annealers ship, wagering computation on slow-sweep ground-state tracking

Going further

References. Griffiths & Schroeter, 3rd ed., Ch. 11 · Sakurai & Napolitano, 2nd ed., §5.6 · Messiah, Vol. II, Ch. XVII · Born & Fock, Z. Phys. 51, 165 (1928) · Kato, J. Phys. Soc. Japan 5, 435 (1950) · Griffiths Example 10.1 (rotating-field spin, 2nd ed. numbering).

Provenance. Equations derived from the time-dependent Schrödinger equation in the instantaneous eigenbasis. A verifier should: (1) reproduce the coupled amplitude equations and the identity ⟨m∣n˙⟩=⟨m∣H˙∣n⟩/(En−Em)\langle m|\dot n\rangle=\langle m|\dot H|n\rangle/(E_n-E_m) by differentiating the eigenproblem (Griffiths §11.1, Sakurai §5.6); (2) confirm the (Em−En)2(E_m-E_n)^2 scaling via integration by parts (Messiah Vol. II Ch. XVII); (3) recompute the spin-½ Berry connection and the ∓12Ω\mp\tfrac12\Omega result against Berry (1984) and Griffiths §11.4; (4) check the Landau–Zener exponent against standard references. All numerical factors (½ solid angle, 2π2\pi in Ω\Omega) are standard textbook values.

Builds on: The Interaction Picture and Time-Dependent Hamiltonians · Time-Dependent Perturbation Theory: General Formalism · First-Order Transition Amplitudes and Probabilities. Leads to: Berry Phase and Geometric Phases · Landau–Zener Transitions · The Born–Oppenheimer Approximation. (Scope note: the base export's Berry-connection/curvature formulas and spin-½ solid-angle derivation are owned by Geometric (Berry) Phase per the chapter ownership map; the Landau-Zener exponent by Landau–Zener Transitions.)

Are the instantaneous eigenstates solutions of the Schrödinger equation?

No — that is the crucial subtlety. They solve the frozen eigenproblem at each instant: a family of snapshots with no dynamics. The theorem's content is that the true solution of the time-dependent Schrödinger equation clings to one snapshot when the driving is slow. At finite speed the true state always contains small admixtures of other snapshots — that admixture is the leak we computed.

Does the system's energy stay constant during adiabatic evolution?

No. The quantum number is conserved, not the energy: E_n(t) tracks the moving level, and the difference is exchanged as work with the agent moving the Hamiltonian (the piston: 0.282 eV to the receding wall) — mirroring the classical statement that the pendulum's E/ω is invariant while E changes.

The theorem needs a gap. What happens with continuous spectra or exact degeneracies?

The non-degenerate gapped case is the theorem treated here (Born-Fock/Kato). An exactly degenerate level has no internal gap, and adiabatic transport within the subspace becomes a path-ordered matrix — the non-Abelian Wilczek-Zee holonomy, graduate-track territory. Continuous spectra need the operator-theoretic treatments descended from Kato's 1950 paper; the naive criterion does not parse there.

Is the geometric phase just a mathematical artifact of the basis choice?

On an open path you can reshuffle it by redefining the phase of |n(t)⟩ — gauge-dependent bookkeeping. But around a closed loop the total γₙ is gauge-invariant: no redefinition removes it, and interferometry measures it. That invariance, the spin-½ solid-angle formula, and the curvature reading are what Geometric (Berry) Phase establishes next.

How is the adiabatic theorem different from first-order time-dependent perturbation theory?

Complementary corners. Perturbation theory assumes the change is small in amplitude (but possibly fast) and expands around a fixed basis. The adiabatic theorem assumes the change is slow in rate (but possibly huge in total) and rides a moving basis. A field rotated through 1000 full turns is hopeless for perturbation theory and trivial for the adiabatic theorem — provided ω ≪ ω₀.

In the spin-cone laboratory, is the adiabatic result an approximation or exact?

The solution is exact at every sweep rate — that is the laboratory's point. "Adiabatic" enters only when reading the exact ceiling (ω sinα/λ)² in the corner ω ≪ ω₀, where it collapses as (ω/ω₀)²sin²α. The theorem is then a verified property, not an assumption: numerical integration of the full Schrödinger equation matches the closed form to 5×10⁻¹⁵.

Why do people say quantum annealing runtime scales with the inverse gap squared?

The leak estimate |c_m| ~ ħ|⟨m|Ḣ|n⟩|/ΔE² must stay small at the worst moment — the minimum gap. Holding the leak fixed while H sweeps a fixed range forces the runtime to grow like 1/Δ²_min (refinements and caveats are Adiabatic Approximation and Validity Conditions' territory). Hard optimization instances tend to develop exponentially small gaps, so the gap — not qubit count — is the honest cost metric.

Before this

Where this leads