Aharonov–Anandan Phase and Non-Adiabatic Generalizations

The article Geometric (Berry) Phase attaches its phase to a slow round trip of external parameters. That framing invites a suspicion: is the geometry an artifact of adiabaticity? In 1987 Aharonov and Anandan answered with a clean no. Take any evolution — fast, driven, wildly non-adiabatic — that is cyclic: the physical state returns to itself, ∣ψ(T)⟩=eiϕ∣ψ(0)⟩|\psi(T)\rangle = e^{i\phi}|\psi(0)\rangle. Subtract from the total phase ϕ\phi the dynamical bookkeeping −1ℏ∫0T⟨H⟩ dt-\tfrac{1}{\hbar}\int_0^T\langle H\rangle\,dt, and what remains is a geometric phase β\beta that depends on exactly one thing: the closed curve the state's ray traces in projective Hilbert space. Not the speed, not the Hamiltonian, not any parameter loop — the curve. This article builds that statement and stress-tests it by machine. We derive the phase decomposition and prove β=∮i⟨ψ~∣dψ~⟩\beta = \oint i\langle\tilde\psi|d\tilde\psi\rangle is gauge- and reparametrization-invariant; we prove the Hamiltonian-independence theorem by driving the same Bloch-sphere circle with two entirely different Hamiltonians — a static field and a rotating field — and finding identical β\beta (4×10−164\times10^{-16} apart) while the dynamical phases differ by radians. For spin-½ we solve the rotating-field problem exactly and watch a cancellation delete every trace of the sweep rate, leaving β=−12Ω\beta = -\tfrac12\Omega with Ω\Omega the solid angle of the circle the spin actually traces — verified for a wobbly non-circular loop to 4×10−134\times10^{-13}. In the slow-drive limit the actual path snaps onto the parameter path and Berry's phase re-emerges as a special case, with the convergence measured.

Key takeaways
  • A cyclic evolution returns the ray to itself: ∣ψ(T)⟩=eiϕ∣ψ(0)⟩|\psi(T)\rangle = e^{i\phi}|\psi(0)\rangle. No adiabaticity is required — only closure of the loop in ray space.
  • The total phase splits uniquely: ϕ=ϕdyn+β\phi = \phi_{\text{dyn}} + \beta, with ϕdyn=−1ℏ∫0T⟨ψ∣H∣ψ⟩ dt\phi_{\text{dyn}} = -\tfrac{1}{\hbar}\int_0^T \langle\psi|H|\psi\rangle\,dt and β\beta the Aharonov–Anandan (AA) geometric phase.
  • β=∮i⟨ψ~∣dψ~⟩\beta = \oint i\langle\tilde\psi|d\tilde\psi\rangle over the closed curve in projective Hilbert space — invariant under gauge choices, reparametrization, and even changes of the Hamiltonian that preserve the curve.
  • Spin-½: β=−12Ω~\beta = -\tfrac12\tilde\Omega, half the solid angle of the actual Bloch-sphere trajectory — exact at any speed, not just adiabatically.
  • Berry's phase is the adiabatic limit: as the drive slows, the actual path converges to the parameter path and β→γBerry\beta \to \gamma_{\text{Berry}} (deviation measured linear in sweep rate).
  • Key idea: strip the unphysical overall phase and watch the ray — whenever it closes a loop, the leftover phase after removing the dynamical part is pure geometry of that loop.
The core insight

Strip the unphysical overall phase from the state and watch the ray. Whenever the ray returns to where it started — a cyclic evolution — the leftover phase, after removing the dynamical phase −1ℏ∫⟨H⟩ dt-\tfrac1\hbar\int\langle H\rangle\,dt, is a geometric quantity fixed entirely by the loop in projective Hilbert space: no adiabaticity, no parameters, no Hamiltonian details.

Where this lives today

The AA formulation is the working language wherever geometric phases are engineered rather than merely observed. Non-adiabatic holonomic quantum gates — realized in trapped ions, NV centers and superconducting circuits — deliberately drive qubits around fast closed Bloch-sphere loops so that the gate angle is a solid angle: rate-independent by Derivation B, hence insensitive to timing noise that would scramble a dynamical phase, and fast enough to beat decoherence precisely because adiabaticity is not required. The design rule is this article's Step 3 cancellation, run in reverse: pick the loop for the phase you want, then pick any convenient Hamiltonian that drives it. The same phase is routine in interferometry. NMR experiments read β\beta by spin-echo sequences that cancel ϕdyn\phi_{\text{dyn}}; polarization optics implements Pancharatnam's version in geometric-phase lenses and metasurface holograms, where a spatially varying loop on the Poincaré sphere prints a phase profile on a beam with no optical path difference at all. In all of these the operational content is the one proved above: the phase survives anything that preserves the loop.

The physical picture

Berry's construction needs a supporting cast: a Hamiltonian H(R)H(\mathbf{R}), a slowly closed parameter loop, an eigenstate faithfully tracking it. Aharonov and Anandan noticed that the phase never cared about the cast — only about the trajectory the state itself performs. The natural stage is ray space: the set of physical states, with overall phase discarded. For a two-level system, ray space is exactly the Bloch sphere. Any evolution whose Bloch vector comes back to its starting point is cyclic, whether it got there by gentle steering or violent driving. Picture a spin in a magnetic field whose direction sweeps a cone of half-angle θ\theta at rate ω\omega. Adiabatically (ω→0\omega \to 0) the spin hugs the field and traces that same cone. At finite ω\omega it cannot keep up: the exact cyclic state traces a different cone, of half-angle θ~\tilde\theta, tilted away from the field's. The AA phase reads the solid angle of the cone the spin actually traces — β=−12Ω(θ~)\beta = -\tfrac12\Omega(\tilde\theta) — and the Berry value −12Ω(θ)-\tfrac12\Omega(\theta) is recovered only in the limit θ~→θ\tilde\theta \to \theta.

Figure
The loop that matters: the actual Bloch path, not the parameter path
z field cone θ (parameter path) actual cyclic path θ̃ great circle (θ̃=90°): β = −π the AA reading β = −½ × (solid angle of the actual path) — at any speed slow drive: θ̃ → θ, β → Berry value −π(1−cos θ) fast drive: θ̃ → π, β → −2π ≡ 0 (trivial) θ=60°, ω=0.5ω₀: θ̃=90° exactly — the spin traces the equator while the field leans

The deepest property is what the phase ignores. Two Hamiltonians that push a state around the same ray-space loop — however different their energies, matrix elements, or time profiles — imprint the same β\beta. The punchline figure shows the numerical experiment: a static field and a rotating field drive the identical equatorial circle; their dynamical phases differ by 5.445.44 rad, their geometric phases agree to 4×10−164\times10^{-16}. On the right, the measured convergence of β\beta to Berry's value as the drive slows.

Figure
Punchline: same loop, different Hamiltonians, same phase — and the measured Berry limit
two drives, one circle (θ̃ = 90°) the shared equatorial loop drive A: static H = (ħω/2)σ_z φ_dyn = 0 drive B: rotating field, θ=60° φ_dyn = −5.44 rad β_A = β_B = −π (Δ = 4×10⁻¹⁶) measured Berry limit, θ = 60° sweep rate ω/ω₀ |β−γ_B| 0.5: 1.571 0.2: 0.543 0.1: 0.254 0.02: 0.048 deviation ∝ ω/ω₀ (linear) — Berry's phase is the slow endpoint

The phase anatomy of a cyclic evolution

Fix a solution of the time-dependent Schrödinger equation over one period. Cyclic means the physical state closes, which constrains the vector only up to phase. That single phase ϕ\phi is the raw observable-in-waiting:

∣ψ(T)⟩  =  eiϕ ∣ψ(0)⟩(1)|\psi(T)\rangle \;=\; e^{i\phi}\,|\psi(0)\rangle \qquad (1)

where: TT — the return time of the ray, not of any external dial; ϕ\phi — the total phase, defined mod 2π2\pi. A stationary state satisfies (1) trivially with ϕ=−ET/ℏ\phi = -E T/\hbar (the standing phase factor of Stationary States and Energy Eigenstates); the interesting cases are loops that genuinely move. Part of ϕ\phi is mere clockwork: every instant, the state accrues phase at rate −⟨H⟩/ℏ-\langle H\rangle/\hbar, exactly as an energy eigenstate would at rate −E/ℏ-E/\hbar. Integrating that rate defines the dynamical phase,

ϕdyn  =  − 1ℏ∫0T⟨ψ(t)∣ H(t) ∣ψ(t)⟩ dt(2)\phi_{\text{dyn}} \;=\; -\,\frac{1}{\hbar}\int_0^T \langle\psi(t)|\,H(t)\,|\psi(t)\rangle\,dt \qquad (2)

where: the expectation value is taken in the actual evolving state — for non-adiabatic motion this is not −EnT/ℏ-E_nT/\hbar, because ∣ψ⟩|\psi\rangle need not be an instantaneous eigenstate. Whatever (2) fails to explain is the Aharonov–Anandan phase, β=ϕ−ϕdyn\beta = \phi - \phi_{\text{dyn}}, and the derivation below shows β\beta is a property of the closed ray-space curve alone.

Derivation ASplitting the total phase into dynamical and geometric parts

We are deriving the AA decomposition: for any cyclic evolution, ϕ−ϕdyn\phi - \phi_{\text{dyn}} equals a loop integral ∮i⟨ψ~∣dψ~⟩\oint i\langle\tilde\psi|d\tilde\psi\rangle built from a single-valued lift ∣ψ~⟩|\tilde\psi\rangle of the ray curve, and that integral is unchanged by rephasing the lift or reparametrizing time. The strategy follows Aharonov–Anandan's 1987 paper: peel the offending phase off by hand, then let the Schrödinger equation tell us what the peeled state's phase velocity is. Watch the energy integral cancel — that cancellation is the whole theorem.

Step 1. Define the lift ∣ψ~(t)⟩=e−if(t) ∣ψ(t)⟩|\tilde\psi(t)\rangle = e^{-i f(t)}\,|\psi(t)\rangle with any smooth ff obeying f(T)−f(0)=ϕf(T) - f(0) = \phi. Then ∣ψ~(T)⟩=e−if(T)eiϕ∣ψ(0)⟩=∣ψ~(0)⟩|\tilde\psi(T)\rangle = e^{-if(T)}e^{i\phi}|\psi(0)\rangle = |\tilde\psi(0)\rangle: the lift is single-valued around the loop — a closed curve in Hilbert space projecting onto the closed ray curve.

Step 2. Differentiate the lift and use the Schrödinger equation iℏ ∂t∣ψ⟩=H∣ψ⟩i\hbar\,\partial_t|\psi\rangle = H|\psi\rangle:

i⟨ψ~∣∂tψ~⟩  =  f˙  +  1ℏ⟨ψ∣H∣ψ⟩.i\langle\tilde\psi|\partial_t\tilde\psi\rangle \;=\; \dot f \;+\; \frac{1}{\hbar}\langle\psi|H|\psi\rangle .

Why this is the pivot: the left side is a purely kinematic object (a connection evaluated along the curve, the same i⟨⋅∣d⋅⟩i\langle\cdot|d\cdot\rangle structure as the Berry connection); the right side splits it into the phase bookkeeping f˙\dot f and the energy rate. No adiabatic assumption was available to be used.

Step 3. Integrate over the period. ∫0Tf˙ dt=ϕ\int_0^T \dot f\,dt = \phi and the energy term integrates to −ϕdyn-\phi_{\text{dyn}} by definition (2):

β  ≡  ∮i⟨ψ~∣ dψ~⟩  =  ϕ−ϕdyn.\beta \;\equiv\; \oint i\langle\tilde\psi|\,d\tilde\psi\rangle \;=\; \phi - \phi_{\text{dyn}} .

Why the energy integral cancels out of the geometry: it has been subtracted explicitly — what remains is the integral of a one-form along the closed lift, with the Hamiltonian appearing nowhere.

Step 4. Invariance. Rephasing the lift, ∣ψ~⟩→eiχ(t)∣ψ~⟩|\tilde\psi\rangle \to e^{i\chi(t)}|\tilde\psi\rangle with χ\chi single-valued around the loop, shifts the integral by −∮dχ=0-\oint d\chi = 0. Reparametrizing time changes dtdt but not the differential dψ~d\tilde\psi — the final expression contains no clock. So β\beta is fixed by the closed curve in ray space together with its orientation, mod 2π2\pi. Verified: the decomposition is checked end-to-end in the article's verification script — an independent TDSE integration of the rotating-field problem reproduces the exact total phase to 2×10−132\times10^{-13} with return overlap ∣⟨χ∣ψ(T)⟩∣=1|\langle\chi|\psi(T)\rangle| = 1, and ϕ−ϕdyn\phi - \phi_{\text{dyn}} lands on the geometric prediction to 10−910^{-9}.

The phase belongs to the curve, not the Hamiltonian

Derivation A leaves a loose thread: many Hamiltonians drive a state around the same ray-space loop. If β\beta deserves the name geometric, all of them must produce the same value. They do, and the proof is two lines.

Go deeper · Derivation B · 1 minHamiltonian-independence of the geometric phase
Derivation BHamiltonian-independence of the geometric phase

We are proving that two different evolutions tracing the same closed curve in ray space, in the same orientation, carry identical β\beta — the theorem that divorces the phase from dynamics. The idea: two state histories with the same rays differ only by a time-dependent phase, and Derivation A already showed such phases drop out of the loop integral. We then confirm by machine with the least similar pair of drives we can build.

Step 1. Let ∣ψ1(t)⟩|\psi_1(t)\rangle and ∣ψ2(t)⟩|\psi_2(t)\rangle be solutions (of possibly different Hamiltonians H1H_1, H2H_2) whose rays coincide at every reparametrized instant. Equal rays means ∣ψ2⟩=eiα(t)∣ψ1⟩|\psi_2\rangle = e^{i\alpha(t)}|\psi_1\rangle for some smooth α(t)\alpha(t).

Step 2. Their single-valued lifts (Step 1 of Derivation A) then differ by a phase eiχ(t)e^{i\chi(t)} that is itself single-valued around the loop (both lifts close). By Derivation A Step 4, the loop integrals agree: β2=β1\beta_2 = \beta_1. The total and dynamical phases can differ arbitrarily — only their difference is pinned.

Step 3. Machine exhibit (the punchline figure). Drive A: static HA=ℏω2σzH_A = \tfrac{\hbar\omega}{2}\sigma_z precesses a state tilted at 90°90° around the equator; since ⟨σz⟩=0\langle\sigma_z\rangle = 0 on the equator, ϕdynA=0\phi_{\text{dyn}}^A = 0 — the total phase is the geometric phase. Drive B: the rotating field of the next section with θ=60°\theta = 60°, ω=0.5 ω0\omega = 0.5\,\omega_0, whose exact cyclic state traces the same equator with ϕdynB=−5.44\phi_{\text{dyn}}^B = -5.44 rad. Computed: βA=βB=−π\beta_A = \beta_B = -\pi with difference 4×10−164\times10^{-16}. Verified: the twin-drive test in the verification script confirms equal circles (θ~\tilde\theta equal to 10−1210^{-12}), equal geometric phases to 4×10−164\times10^{-16}, and dynamical phases differing by 5.445.44 rad.

The right arena for all of this is projective Hilbert space CPN−1\mathbb{CP}^{N-1}: states modulo phase. For a two-level system that space is the Bloch sphere — the half-angle dictionary of Spinors and Two-Component Wavefunctions maps every ray to a unit vector, and cyclic evolutions to closed circuits. The AA phase is a functional on closed circuits; the Berry connection and curvature of the neighboring article organize the same geometry field-theoretically, and Berry Connection and Berry Curvature takes that formal story up.

Spin-½ exactly: the solid angle of the path actually traced

Now the flagship solvable case. Let a field of fixed magnitude sweep a cone: H(t)=ℏω02 B^(t)⋅σH(t) = \tfrac{\hbar\omega_0}{2}\,\hat{\mathbf{B}}(t)\cdot\boldsymbol\sigma with B^=(sin⁡θcos⁡ωt, sin⁡θsin⁡ωt, cos⁡θ)\hat{\mathbf{B}} = (\sin\theta\cos\omega t,\ \sin\theta\sin\omega t,\ \cos\theta). The rotating-frame construction that solves the driven two-level problem exactly (the machinery of Transition Probability for a Periodic Two-Level System) turns this into a static problem: in the frame rotating at ω\omega about zz,

Heff  =  ℏ2[ ω0sin⁡θ σx  +  (ω0cos⁡θ−ω) σz](3)H_{\text{eff}} \;=\; \frac{\hbar}{2}\Big[\,\omega_0\sin\theta\,\sigma_x \;+\; (\omega_0\cos\theta - \omega)\,\sigma_z\Big] \qquad (3)

where: the −ωσz/2-\omega\sigma_z/2 piece is the frame's fictitious field. The eigenstates of (3) are the cyclic states: dragged back to the lab frame they trace perfect cones about zz and return exactly after T=2π/ωT = 2\pi/\omega. Their cone half-angle is

tan⁡θ~  =  ω0sin⁡θω0cos⁡θ−ω(4)\tan\tilde\theta \;=\; \frac{\omega_0\sin\theta}{\omega_0\cos\theta - \omega} \qquad (4)

where: θ~\tilde\theta — the tilt of the actual trajectory. Slow drive (ω≪ω0\omega \ll \omega_0): θ~→θ\tilde\theta \to \theta, the spin hugs the field. Fast drive: θ~→π\tilde\theta \to \pi, the state pins to the (inverted) rotation axis. And when ω=ω0cos⁡θ\omega = \omega_0\cos\theta the σz\sigma_z term vanishes: θ~=90°\tilde\theta = 90° — the spin rides the equator, a great circle, even though the field still leans at θ\theta.

Go deeper · Derivation C · 1 minExact AA phase for a spin in a rotating field
Derivation CExact AA phase for a spin in a rotating field

We are computing β\beta for the cyclic states of the rotating-field problem exactly — no adiabatic approximation — and watching every trace of the drive details cancel, leaving −12×-\tfrac12\times(solid angle of the actual cone). This is the calculation the article's claims stand on: the total phase from the exact propagator, the dynamical phase from a time-independent expectation value, and a cancellation of the effective splitting λ\lambda between them.

Step 1. Total phase. With a=ω0sin⁡θa = \omega_0\sin\theta, b=ω0cos⁡θ−ωb = \omega_0\cos\theta - \omega, λ=a2+b2\lambda = \sqrt{a^2 + b^2}, the lab propagator over one period factorizes as U(T)=e−iωTσz/2 e−iHeffT/ℏU(T) = e^{-i\omega T\sigma_z/2}\,e^{-iH_{\text{eff}}T/\hbar}. Acting on the upper cyclic state ∣χ⟩|\chi\rangle (eigenvalue +ℏλ/2+\hbar\lambda/2), and using ωT=2π\omega T = 2\pi so e−iπσz=−1e^{-i\pi\sigma_z} = -\mathbb{1}:

U(T)∣χ⟩=− e−iλT/2 ∣χ⟩    ⟹    ϕ=π−λT2(mod2π).U(T)|\chi\rangle = -\,e^{-i\lambda T/2}\,|\chi\rangle \;\;\Longrightarrow\;\; \phi = \pi - \frac{\lambda T}{2} \pmod{2\pi}.

Step 2. Dynamical phase. ⟨H⟩\langle H\rangle is time-independent on a cyclic state (checked to 10−1210^{-12} in the script). At t=0t=0 the Bloch vector sits at tilt θ~\tilde\theta and the field at θ\theta, so ⟨H⟩=ℏω02cos⁡(θ−θ~)\langle H\rangle = \tfrac{\hbar\omega_0}{2}\cos(\theta - \tilde\theta). Expanding with a=λsin⁡θ~a = \lambda\sin\tilde\theta, b=λcos⁡θ~b = \lambda\cos\tilde\theta:

ω0cos⁡(θ−θ~)=(b+ω)cos⁡θ~+asin⁡θ~=λ+ωcos⁡θ~    ⟹    ϕdyn=−T2(λ+ωcos⁡θ~).\omega_0\cos(\theta - \tilde\theta) = (b+\omega)\cos\tilde\theta + a\sin\tilde\theta = \lambda + \omega\cos\tilde\theta \;\;\Longrightarrow\;\; \phi_{\text{dyn}} = -\frac{T}{2}\big(\lambda + \omega\cos\tilde\theta\big).

Why this expansion: it re-expresses the field direction in terms of the trajectory's own geometry — the step that lets the dynamics cancel.

Step 3. Subtract. The λT/2\lambda T/2 terms — the only place the effective splitting, hence the drive strength and detuning, appear — cancel identically:

β=ϕ−ϕdyn=π+ωT2cos⁡θ~=π+πcos⁡θ~  ≡  − π (1−cos⁡θ~)(mod2π).\beta = \phi - \phi_{\text{dyn}} = \pi + \frac{\omega T}{2}\cos\tilde\theta = \pi + \pi\cos\tilde\theta \;\equiv\; -\,\pi\,(1 - \cos\tilde\theta) \pmod{2\pi}.

That is −12 Ω(θ~)-\tfrac12\,\Omega(\tilde\theta), half the solid angle Ω=2π(1−cos⁡θ~)\Omega = 2\pi(1-\cos\tilde\theta) of the cone actually traced. The sweep rate ω\omega survives only inside θ~\tilde\theta — only through the shape of the path, exactly as the general theorem demands.

Step 4. Beyond circles. The law is not about cones. Drive the spin with the parallel-transport Hamiltonian H=ℏ2(r×r˙)⋅σH = \tfrac{\hbar}{2}(\mathbf{r}\times\dot{\mathbf{r}})\cdot\boldsymbol\sigma along a wobbly closed path θ(t)=1.05+0.22sin⁡(3φ)\theta(t) = 1.05 + 0.22\sin(3\varphi): then ⟨H⟩=0\langle H\rangle = 0 identically (the field is perpendicular to the Bloch vector), so the total phase is the geometric phase. The machine returns β=−1.5972843\beta = -1.5972843, matching −12∮(1−cos⁡θ) dφ-\tfrac12\oint(1-\cos\theta)\,d\varphi to 4×10−134\times10^{-13}. Verified: solid-angle law confirmed for the cone family at three drive rates (10−910^{-9}), by independent TDSE integration (2×10−132\times10^{-13} on the total phase), and on the non-circular path above (4×10−134\times10^{-13}); the great-circle case gives β=−π\beta = -\pi exactly.

The adiabatic limit: recovering Berry

Send ω/ω0→0\omega/\omega_0 \to 0 in (4): θ~→θ\tilde\theta \to \theta, the actual cone collapses onto the field cone, and β→−π(1−cos⁡θ)\beta \to -\pi(1-\cos\theta) — the solid-angle result carried by Geometric (Berry) Phase for exactly this setup. Berry's phase is thus the adiabatic edge of the AA phase: the special case in which the ray-space loop is inherited from a parameter loop. The convergence is measurable and linear. At θ=60°\theta = 60°, where the Berry value is −π/2-\pi/2, the computed deviations ∣β−γBerry∣|\beta - \gamma_{\text{Berry}}| are

ω/ω0: 0.5→1.571,0.2→0.542,0.1→0.254,0.05→0.122,0.02→0.048 rad(5)\omega/\omega_0:\ 0.5 \to 1.571,\quad 0.2 \to 0.542,\quad 0.1 \to 0.254,\quad 0.05 \to 0.122,\quad 0.02 \to 0.048\ \text{rad} \qquad (5)

where: the first entry is the great-circle case — at half the Larmor rate the trajectory is already 30°30° away from the field cone, and the 'correction' is as large as the Berry phase itself. A five-fold slowdown (0.1→0.020.1 \to 0.02) shrinks the deviation five-fold: the leading error is first order in ω/ω0\omega/\omega_0. Adiabaticity here is a convergence statement about paths, not a precondition for geometry.

Limiting cases

caseresultmeaning
ω/ω0→0\omega/\omega_0 \to 0 (adiabatic)θ~→θ\tilde\theta \to \theta, β→−π(1−cos⁡θ)\beta \to -\pi(1-\cos\theta)actual path → parameter path; Berry's phase recovered as the slow limit
ω=ω0cos⁡θ\omega = \omega_0\cos\theta (great circle)θ~=90°\tilde\theta = 90°, β=−π\beta = -\pimaximal geometric phase; with a static drive this loop has ϕdyn=0\phi_{\text{dyn}} = 0 — the phase is pure geometry
ω/ω0→∞\omega/\omega_0 \to \infty (fast drive)θ~→π\tilde\theta \to \pi, β→−2π≡0\beta \to -2\pi \equiv 0the state pins to the rotation axis; the loop degenerates and the phase trivializes
stationary eigenstateray never moves; β=0\beta = 0, ϕ=−ET/ℏ\phi = -ET/\hbarthe trivial cyclic case — all phase is dynamical clockwork
⟨H⟩=0\langle H\rangle = 0 along the pathϕ=β\phi = \betaparallel-transport drives: the total phase is directly geometric (the wobbly-path exhibit)
non-cyclic evolutionβ\beta undefined as statedthe ray must close; the Samuel–Bhandari geodesic-closure extension handles open paths (one line, not owned here)

Worked examples

Example 1 — a great-circle spin, end to end. A two-level system has Larmor frequency ω0/2π=100\omega_0/2\pi = 100 MHz. Its field sweeps a cone of half-angle θ=60°\theta = 60° at the great-circle rate. Find the drive frequency, the period, and both phases of the upper cyclic state.

Step 1. Great-circle condition: ω=ω0cos⁡60°=ω0/2\omega = \omega_0\cos 60° = \omega_0/2, so ω/2π=50\omega/2\pi = 50 MHz and T=2π/ω=20 nsT = 2\pi/\omega = \mathbf{20\ ns}.

Step 2. Geometric phase: θ~=90°\tilde\theta = 90°, solid angle Ω=2π\Omega = 2\pi, so β=−12Ω=\beta = -\tfrac12\Omega = \boldsymbol{-\pi} — independent of ω0\omega_0, ω\omega, and TT individually (rate-independence checked in the script to 10−910^{-9}).

Step 3. Dynamical phase: from Derivation C Step 2, ϕdyn=−T2(λ+ωcos⁡θ~)=−λT/2\phi_{\text{dyn}} = -\tfrac{T}{2}(\lambda + \omega\cos\tilde\theta) = -\lambda T/2 on the great circle. Numerically ∣ϕdyn∣=5.44|\phi_{\text{dyn}}| = 5.44 rad — comparable to ∣β∣=π|\beta| = \pi and carrying all the Hamiltonian details that β\beta forgot.

Step 4. Reading: an interferometric comparison of this run against a reference that cancels the 5.445.44 rad of clockwork isolates a clean −π-\pi: the state comes back minus itself geometrically — the AA face of the spinor sign flip.

Example 2 — one field cone, three trajectories. Keep θ=60°\theta = 60°, ω0=1\omega_0 = 1 (units of the Larmor rate), and drive at ω=0.2, 0.5, 0.8\omega = 0.2,\ 0.5,\ 0.8. Find each actual cone and its AA phase.

Step 1. From (4), tan⁡θ~=sin⁡60°/(cos⁡60°−ω)\tan\tilde\theta = \sin 60°/(\cos 60° - \omega), giving θ~=70.9°, 90.0°, 109.1°\tilde\theta = 70.9°,\ 90.0°,\ 109.1° — the trajectory leans further from the axis as the drive speeds up, crossing the equator exactly at ω=0.5\omega = 0.5.

Step 2. Solid angles Ω=2π(1−cos⁡θ~)=4.227, 6.283, 8.340\Omega = 2\pi(1-\cos\tilde\theta) = 4.227,\ 6.283,\ 8.340 sr.

Step 3. Phases β=−Ω/2=−2.113, −3.142, −4.170\beta = -\Omega/2 = \mathbf{-2.113,\ -3.142,\ -4.170} rad (the last ≡+2.113\equiv +2.113 mod 2π2\pi — only eiβe^{i\beta} is physical). All three verified against the exact propagator to 10−910^{-9}; the Berry value −π/2=−1.571-\pi/2 = -1.571 is approached only far to the left of this table, per (5).

Origins · 1 minYakir Aharonov and how this was discovered
YA
Yakir Aharonov
1932– · Haifa, Bristol, Columbia (South Carolina), Tel Aviv

Yakir Aharonov has spent a career locating the physics that hides in phases. As a doctoral student with David Bohm in Bristol he showed in 1959 that a charged particle records the vector potential of a solenoid it never touches — the Aharonov–Bohm effect, the archetype every geometric-phase argument leans on (its story is told with Geometric (Berry) Phase). Twenty-eight years later, working with Jeeva Anandan at the University of South Carolina, he pulled Berry's discovery off its adiabatic scaffolding: the 1987 paper reframed the geometric phase as a property of any cyclic evolution in projective Hilbert space, with Berry's phase as its slow limit and the solid-angle rule promoted to the actual trajectory. The pair went further in 1988, showing ray space carries a natural metric that turns quantum evolution itself into geometry. Aharonov's other signature — weak values and weak measurement (1988, with Albert and Vaidman) — extends the same instinct: interrogate the quantum phase gently and it answers with structure no textbook expectation value contains.

History
1956
— Pancharatnam formulates a phase for cyclic changes of light polarization, the optical ancestor of the geometric phase
1959
— Aharonov and Bohm predict the solenoid phase shift; the potential acquires loop-level physical meaning
1984
— Berry publishes the adiabatic geometric phase for slowly cycled Hamiltonians
1987
— Aharonov and Anandan remove adiabaticity: the phase belongs to any closed loop in projective Hilbert space
1988
— Samuel and Bhandari extend the construction to non-cyclic evolutions via geodesic closure
1988
— Suter, Mueller and Pines measure the AA phase directly by NMR interferometry

Going further

References. Aharonov & Anandan, Phys. Rev. Lett. 58, 1593 (1987) · Anandan & Aharonov, Phys. Rev. D 38, 1863 (1988) · Berry, Proc. R. Soc. Lond. A 392, 45 (1984) · Sakurai & Napolitano, Modern Quantum Mechanics, geometric-phase section · Griffiths & Schroeter, 3rd ed., Ch. 11 · Cohen-Tannoudji, Diu & Laloë, geometric-phase complements · Bohm, Mostafazadeh, Koizumi, Niu & Zwanziger, The Geometric Phase in Quantum Systems (Springer, 2003)

Provenance. Phase-decomposition derivation reproduces Aharonov–Anandan, Phys. Rev. Lett. 58, 1593 (1987), Eqs. (3)–(5); spin-1/2 solid-angle result β=−π(1−cos⁡θ)\beta=-\pi(1-\cos\theta) checkable against Sakurai & Napolitano Modern QM geometric-phase section and Griffiths Ch. 11. Verifier should: (1) confirm the cancellation of the energy integral in the lift derivation; (2) check the rotating-frame spin solution and that ω\omega cancels from β\beta; (3) verify Ω=2π(1−cos⁡θ)\Omega=2\pi(1-\cos\theta) as the solid angle of a cone of half-angle θ\theta; (4) confirm the gauge term ∮dχ=0\oint\mathrm{d}\chi=0 for single-valued χ\chi. No nonstandard numerical constants are asserted.

Builds on: The Interaction Picture and Time-Dependent Hamiltonians · Dyson Series and the Time-Ordered Evolution Operator · Transition Probability for a Periodic Two-Level System Leads to: Semiclassical Theory of Atom–Radiation Interaction · Holonomic Quantum Computation · Pancharatnam Phase in Optics

Do I need a Hamiltonian at all to define the AA phase?

Only to produce the loop and to define the dynamical phase you subtract. The value of β\beta itself is a functional of the closed curve in ray space — Derivation B shows any Hamiltonian driving the same curve gives the same answer. In the extreme case of a parallel-transport drive with ⟨H⟩=0\langle H\rangle = 0, nothing is subtracted at all and the total phase is directly geometric. That is why the modern view treats β\beta as belonging to the geometry of projective Hilbert space, with the Hamiltonian demoted to a means of transportation.

If the overall phase of a state is unphysical, what exactly is measured?

A relative phase. Split a system into two arms (interferometry, or two spin components), let one undergo the cyclic evolution, recombine, and the fringe shift reads ϕ=ϕdyn+β\phi = \phi_{\text{dyn}} + \beta. Isolating β\beta is an engineering problem — spin-echo and phase-cycling sequences cancel the dynamical part (that is precisely how the 1988 NMR experiment measured the AA phase). The geometric piece survives such cancellation exactly because it does not care about the energies being echoed away.

The rotating field never stops leaning at 60° — why does the spin trace the equator instead?

Because the cyclic state is an eigenstate of the rotating-frame Hamiltonian (3), and that operator includes the fictitious −ωσz/2-\omega\sigma_z/2 field of the frame transformation. At ω=ω0cos⁡θ\omega = \omega_0\cos\theta the fictitious field exactly cancels the zz-component of the real one, leaving an effective field in the equatorial plane; the cyclic state aligns with it, not with the lab field. The lab-frame trajectory is that tilted state carried around zz once per period — the equator. The AA phase reads this actual circle, which is why it is −π-\pi and not the adiabatic −π/2-\pi/2.

Is the evolution's period T doing any work in the answer?

None beyond defining the loop. β=∮i⟨ψ~∣dψ~⟩\beta = \oint i\langle\tilde\psi|d\tilde\psi\rangle contains no dtdt — reparametrize the traversal arbitrarily and the integral is untouched (Derivation A, Step 4). The verification script exhibits this concretely: the great-circle phase is −π-\pi at any drive rate satisfying the great-circle condition, while the period and the dynamical phase change underneath it.

What happens if the state only approximately returns to its initial ray?

Then ϕ\phi in equation (1) is not defined and neither, strictly, is β\beta. Two standard responses: close the gap explicitly — follow the evolution with a short geodesic in ray space connecting endpoint to start (the Samuel–Bhandari construction, 1988), which reduces to the AA phase when the gap vanishes; or accept a small systematic error of order the closure defect. In experiments one usually engineers exact cyclicity instead — as with the rotating-frame eigenstates here, which return exactly by construction.

Is the minus sign in β = −Ω/2 a convention?

Half of it. The sign of β\beta flips with the orientation of the loop and with which spin state (upper or lower) you transport — the lower cyclic state acquires +Ω~/2+\tilde\Omega/2. What is not conventional: the relative sign between the two states, the magnitude 12Ω~\tfrac12\tilde\Omega, and the mod-2π2\pi identification (Example 2's −4.170≡+2.113-4.170 \equiv +2.113). The verification script pins the stated conventions: upper state, loop traversed with increasing φ\varphi, β=−π(1−cos⁡θ~)\beta = -\pi(1-\cos\tilde\theta).

Before this

Where this leads