The Time-Evolution Operator

The formalism chapter built the stage: states as vectors, observables as Hermitian operators, measurement as projection. This chapter sets the stage in motion — and opens with a change of strategy. Instead of solving the Schrödinger equation state by state, quantum dynamics packages all motion into one object: the time-evolution operator U^(t,t0)\hat U(t, t_0), the unitary map taking any state at t0t_0 to its value at tt. Solve for U^\hat U once and every initial condition is solved forever. The payoff is structural. Three physical requirements — evolving for zero time does nothing, evolving in two stages equals evolving in one, and total probability stays 1 — force the operator to be unitary and its generator to be Hermitian. The Schrödinger equation becomes the infinitesimal fingerprint of those requirements, with the Hamiltonian revealed as the generator of time translation, exactly as momentum generates displacement. That is the viewpoint this chapter runs on: spectral decompositions, the Dyson series, the Schrödinger and Heisenberg pictures are all one operator U^\hat U worn different ways. The completeness machinery is already paid for. Completeness Relation and Resolution of Identity assembled functions of operators from projectors and stopped, deliberately, at the eigenphases of e−iH^t/ℏe^{-i\hat H t/\hbar}. This article cashes that promise: the phases become precession, oscillation, and the two-level beat it left waiting.

Key takeaways
  • U^(t,t0)\hat U(t, t_0) propagates states: ∣ψ(t)⟩=U^(t,t0)∣ψ(t0)⟩|\psi(t)\rangle = \hat U(t, t_0)|\psi(t_0)\rangle — one operator, every initial condition.
  • Physics forces its shape: identity at t=t0t = t_0, composition U^(t2,t0)=U^(t2,t1)U^(t1,t0)\hat U(t_2, t_0) = \hat U(t_2, t_1)\hat U(t_1, t_0), unitarity U^†U^=1^\hat U^\dagger \hat U = \hat 1 (probability conservation).
  • The infinitesimal step U^(t0+dt,t0)=1^−iℏH^ dt\hat U(t_0 + dt, t_0) = \hat 1 - \tfrac{i}{\hbar}\hat H\,dt identifies H^\hat H as the generator; it obeys the operator Schrödinger equation iℏ ∂tU^=H^U^i\hbar\,\partial_t \hat U = \hat H\hat U.
  • Unitarity ⇔\Leftrightarrow Hermiticity of H^\hat H — conservation of probability and reality of energy are the same fact.
  • Time-independent H^\hat H: U^=e−iH^(t−t0)/ℏ\hat U = e^{-i\hat H(t - t_0)/\hbar}, acting on each energy eigenstate as a pure phase e−iEnt/ℏe^{-iE_n t/\hbar}.
  • Three regimes: constant H^\hat H (one exponential) · [H^(t′),H^(t′′)]=0[\hat H(t'), \hat H(t'')] = 0 (integrate the exponent) · non-commuting (time-ordered exponential — Time-Ordered Exponential and the Dyson Series).
  • Key idea: all dynamics of a closed system is one unitary operator generated by the Hamiltonian — its unitarity is probability conservation, its composition law is the consistency of time itself.
The core insight

All quantum dynamics of a closed system is the action of one unitary operator U(t,t0)U(t,t_0) generated by the Hamiltonian via iℏ ∂tU=HUi\hbar\,\partial_t U = H U; its unitarity is conservation of probability and its group law U(t2,t1)U(t1,t0)=U(t2,t0)U(t_2,t_1)U(t_1,t_0)=U(t_2,t_0) is the consistency of time evolution.

Where this lives today

Every gate in a quantum computer is a hand-built U^\hat U: a calibrated pulse making the qubit Hamiltonian act for exactly the right duration. Digital quantum simulation is the product-limit construction run in reverse — approximate e−iH^t/ℏe^{-i\hat Ht/\hbar} for an intractable H^\hat H by Trotterized slices of easy pieces, the 1/N1/N error scaling of equation (2) setting the gate budget. And the 28-GHz-per-tesla precession of Example 1 is working hardware: electron spin resonance reads molecular environments off exactly that clock, while its nuclear cousin runs every MRI scanner. When this chapter later derives pictures, propagators, and perturbation series, each will be this one operator re-dressed.

The physical picture

Think of U^(t,t0)\hat U(t, t_0) as a rigid rotation of the whole Hilbert space. A rotation preserves lengths and angles; U^\hat U preserves norms and inner products — which is to say, probabilities. The Hamiltonian is the axis and rate of the rotation: just as momentum generates translations (the Weyl-form card of Canonical Commutation Relation), energy generates motion through time. States with definite energy are the directions the rotation does not mix — each spins in place at its own phase rate En/ℏE_n/\hbar — and all visible dynamics is interference between phases spinning at different rates.

Figure
One operator moves every state — U as a rigid rotation of Hilbert space, H as its generator
|ψ(t₀)⟩ |ψ(t)⟩ = U(t,t₀)|ψ(t₀)⟩ U rotates norm ‖ψ‖ = 1 preserved: motion stays on the sphere generator: Ĥ infinitesimal step U(t₀+dt) = 1 − iĤdt/ℏ energy = rate of rotation one operator, every state: solve for U once, apply to any initial condition

The second picture is composition. Evolving from 9:00 to 9:02 and then 9:02 to 9:05 must equal evolving from 9:00 to 9:05 — time evolution has no seams: U^(t2,t0)=U^(t2,t1)U^(t1,t0)\hat U(t_2, t_0) = \hat U(t_2, t_1)\hat U(t_1, t_0), the earlier interval acting first (rightmost). This unassuming rule, plus unitarity, is nearly the whole theory: it forces an infinitesimal-generator structure, and for constant H^\hat H it upgrades the family to a one-parameter group — operators depending only on elapsed time, multiplying by adding arguments. Third: what makes the problem easy or hard is not whether H^\hat H depends on time, but whether H^\hat H at one time commutes with H^\hat H at another. Constant Hamiltonians give a clean exponential. Time-dependent but self-commuting ones just integrate the exponent. Non-commuting ones are like composing rotations about different axes — order matters at every instant, and the honest bookkeeping is the time-ordered exponential built two articles ahead.

Figure
Which U for which H — the three-regime decision tree of quantum dynamics
does H depend on t? no yes U = e^(−iH(t−t₀)/ℏ) one exponential · group law [H(t′), H(t″)] = 0 for all pairs? yes no U = e^(−(i/ℏ)∫H dt′) exponent integrates · still exact 𝒯-ordered exp Dyson series article the fork that matters is self-commutation, not time dependence itself

The defining requirements — and the equation they force

Quantum dynamics asks for a family of operators doing three things. The contract:

U^(t0,t0)=1^,U^(t2,t0)=U^(t2,t1) U^(t1,t0),U^†U^=1^(1)\hat U(t_0, t_0) = \hat 1, \qquad \hat U(t_2, t_0) = \hat U(t_2, t_1)\,\hat U(t_1, t_0), \qquad \hat U^\dagger \hat U = \hat 1 \qquad (1)

where: the first clause says zero elapsed time changes nothing; the second is seamlessness (read right to left: earlier segment first); the third is unitarity — the inner-product-preserving class built at Unitary Operators and Transformations — demanded so a normalized state stays normalized forever. Linearity is also assumed: evolution must respect superpositions, or interference physics falls apart. From this contract alone, the dynamical law follows. That construction — Sakurai's opening move of quantum dynamics — is the first derivation.

Derivation AFrom three physical requirements to the operator Schrödinger equation

We derive the equation of motion for the evolution operator without assuming the Schrödinger equation — the requirements of contract (1) plus smoothness in time do all the work. The strategy: look at an infinitesimally short evolution, expand to first order, and let unitarity and composition dictate the structure of the linear term. What emerges is a Hermitian generator whose physical identity — the Hamiltonian — is fixed by the classical correspondence of energy with time translation, and an operator differential equation that contains the state-level Schrödinger equation as a corollary. Watch how little input produces how much output.

Step 1. Expand a short step. For an infinitesimal interval dtdt, smoothness lets us write

U^(t0+dt,t0)=1^−iℏ Ω^ dt+O(dt2)\hat U(t_0 + dt, t_0) = \hat 1 - \frac{i}{\hbar}\,\hat\Omega\,dt + O(dt^2)

with some operator Ω^\hat\Omega to be constrained. The factor −i/ℏ-i/\hbar is a convention that makes Ω^\hat\Omega come out Hermitian with units of energy. Why start linear: the identity at dt=0dt = 0 is clause one of the contract; any smooth family departs from it linearly, and everything hinges on what the linear coefficient must be.

Step 2. Impose unitarity at first order:

U^†U^=(1^+iℏΩ^†dt)(1^−iℏΩ^ dt)=1^+iℏ(Ω^†−Ω^) dt+O(dt2)\hat U^\dagger \hat U = \Big(\hat 1 + \frac{i}{\hbar}\hat\Omega^\dagger dt\Big)\Big(\hat 1 - \frac{i}{\hbar}\hat\Omega\,dt\Big) = \hat 1 + \frac{i}{\hbar}\big(\hat\Omega^\dagger - \hat\Omega\big)\,dt + O(dt^2)

For this to equal 1^\hat 1 at order dtdt, we need Ω^†=Ω^\hat\Omega^\dagger = \hat\Omega: the generator must be Hermitian. An observable has appeared out of pure bookkeeping. Why this is remarkable: probability conservation alone manufactures a Hermitian operator — the mirror-symmetric class of Hermitian Operators and Observables — before any physics names it.

Step 3. Name the generator. Which observable generates time translation? Classical mechanics answers: the Hamiltonian generates time evolution, as momentum generates displacement — and the translation-operator analogy e−iap^/ℏe^{-ia\hat p/\hbar} (the Weyl card of Canonical Commutation Relation) makes ℏ\hbar's bookkeeping match. Set Ω^=H^\hat\Omega = \hat H, the energy operator. This identification is a physical postulate — the same one Time-Dependent Schrödinger Equation makes at the state level — judged by its scorecard, not derived.

Step 4. Build the differential equation from composition. Compose a finite evolution with one more infinitesimal step:

U^(t+dt,t0)=U^(t+dt,t) U^(t,t0)=(1^−iℏH^(t) dt)U^(t,t0)\hat U(t + dt, t_0) = \hat U(t + dt, t)\,\hat U(t, t_0) = \Big(\hat 1 - \frac{i}{\hbar}\hat H(t)\,dt\Big)\hat U(t, t_0)

Subtract U^(t,t0)\hat U(t, t_0), divide by dtdt, and take the limit:

iℏ ∂∂tU^(t,t0)=H^(t) U^(t,t0),U^(t0,t0)=1^i\hbar\,\frac{\partial}{\partial t}\hat U(t, t_0) = \hat H(t)\,\hat U(t, t_0), \qquad \hat U(t_0, t_0) = \hat 1

Why this contains the state equation: apply both sides to a fixed ∣ψ(t0)⟩|\psi(t_0)\rangle and the operator equation becomes iℏ ∂t∣ψ(t)⟩=H^∣ψ(t)⟩i\hbar\,\partial_t|\psi(t)\rangle = \hat H|\psi(t)\rangle — the ket-space law of the previous article, now recognized as the infinitesimal fingerprint of unitary, seamless evolution. *Verified: the exponential solution satisfies the operator equation symbolically (sympy, 2×2 H^\hat H: i dU^/dt−H^U^=0i\,d\hat U/dt - \hat H\hat U = 0 exactly), and the composition law U^(1.9)=U^(0.7)U^(1.2)\hat U(1.9) = \hat U(0.7)\hat U(1.2) holds to 10−1210^{-12} on a random 3×3 Hermitian Hamiltonian (verify script, 20 checks green).*

Unitarity and Hermiticity are the same fact

The construction above went from unitarity to a Hermitian generator. The traffic runs both ways, and the round trip is worth one compact derivation because it is the exam-grade statement of "probability is conserved."

Go deeper · Derivation B · 1 minHermiticity of H is exactly unitarity of U
Derivation BHermiticity of H is exactly unitarity of U

We prove that U^†U^\hat U^\dagger\hat U is constant in time if and only if H^=H^†\hat H = \hat H^\dagger, so the initial condition U^(t0,t0)=1^\hat U(t_0, t_0) = \hat 1 pins U^†U^=1^\hat U^\dagger\hat U = \hat 1 forever. The strategy is to differentiate the product U^†U^\hat U^\dagger\hat U and watch the two equations of motion — one for U^\hat U, one for its adjoint — cancel precisely when the Hamiltonian equals its own adjoint. This is the operator-level twin of the wave-mechanics norm-conservation proof, which lives at Conservation of Normalization and is recalled, not re-run, here.

Step 1. Take the adjoint of the operator Schrödinger equation. Using the product-reversal rule of Hermitian Conjugate and Adjoint Operators:

iℏ ∂tU^=H^ U^⟹−iℏ ∂tU^†=U^† H^†i\hbar\,\partial_t \hat U = \hat H\,\hat U \quad\Longrightarrow\quad -i\hbar\,\partial_t \hat U^\dagger = \hat U^\dagger\,\hat H^\dagger

Step 2. Differentiate the product and insert both equations:

∂t(U^†U^)=(∂tU^†) U^+U^† (∂tU^)=iℏ U^†(H^†−H^)U^\partial_t\big(\hat U^\dagger\hat U\big) = (\partial_t\hat U^\dagger)\,\hat U + \hat U^\dagger\,(\partial_t\hat U) = \frac{i}{\hbar}\,\hat U^\dagger\big(\hat H^\dagger - \hat H\big)\hat U

Why the structure is decisive: the derivative of the product is sandwiched between U^†\hat U^\dagger and U^\hat U, so it vanishes for every state exactly when H^†−H^=0\hat H^\dagger - \hat H = 0 — no cancellation is available from anything else.

Step 3. Read both directions. If H^=H^†\hat H = \hat H^\dagger: the derivative vanishes, U^†U^\hat U^\dagger\hat U keeps its initial value 1^\hat 1, and every norm ⟨ψ(t)∣ψ(t)⟩=⟨ψ(t0)∣U^†U^∣ψ(t0)⟩\langle\psi(t)|\psi(t)\rangle = \langle\psi(t_0)|\hat U^\dagger\hat U|\psi(t_0)\rangle is frozen — probability conserved. If instead H^\hat H has an anti-Hermitian part −iΓ/2-i\Gamma/2: the same algebra gives exponential norm decay — precisely the optical-potential model of Conservation of Normalization, whose τ=ℏ/Γ\tau = \hbar/\Gamma lifetime is the controlled way to break this theorem on purpose. Verified: numerically — ∥U^†U^−1^∥<10−12\|\hat U^\dagger\hat U - \hat 1\| < 10^{-12} and norm conservation to 10−1410^{-14} for a random 3×3 Hermitian H^\hat H; inverse = adjoint (U^(−t)=U^(t)†\hat U(-t) = \hat U(t)^\dagger) to 10−1210^{-12} (verify script).

Solving the contract: the exponential and its spectral form

For a time-independent Hamiltonian the operator equation is solved by compounding the infinitesimal step NN times over the interval and letting N→∞N \to \infty — legal because every factor involves the same H^\hat H and commuting factors multiply like numbers:

U^(t,t0)=lim⁡N→∞(1^−iℏH^ t−t0N) ⁣N=exp⁡ ⁣[−iℏH^ (t−t0)](2)\hat U(t, t_0) = \lim_{N\to\infty}\Big(\hat 1 - \frac{i}{\hbar}\hat H\,\frac{t - t_0}{N}\Big)^{\!N} = \exp\!\Big[-\frac{i}{\hbar}\hat H\,(t - t_0)\Big] \qquad (2)

where: the exponential means its power series, and the limit is the operator version of (1+x/N)N→ex(1 + x/N)^N \to e^x. The verify script's error ladder makes the convergence concrete: truncating at N=10,100,1000N = 10, 100, 1000 steps leaves errors 0.24,0.022,0.00220.24, 0.022, 0.0022 — shrinking like 1/N1/N, which is also why product-formula ("Trotterized") quantum simulation works. Differentiating (2) returns the operator Schrödinger equation, closing the loop. How does the exponential act? Completeness Relation and Resolution of Identity assembled functions of an operator from its spectral projectors — eigenvalues go through the function, eigenvectors stand still — and applied it to e−iH^t/ℏe^{-i\hat Ht/\hbar}, stopping at the eigenphases. Cash the promise:

U^(t,0)=∑ne−iEnt/ℏ ∣En⟩⟨En∣(3)\hat U(t, 0) = \sum_n e^{-iE_n t/\hbar}\,|E_n\rangle\langle E_n| \qquad (3)

where: each energy eigenstate is multiplied by a pure phase turning at rate En/ℏE_n/\hbar — the spectral form of time evolution. A single eigenstate goes nowhere observable (the stationary-state story, next articles); superpositions evolve only through phase differences. Equation (3) is verified against the direct matrix exponential to 10−1210^{-12} on the two-level system of Eigenvalues and Eigenstates — and Example 2 below turns its eigenphase drift into a real oscillation. The full craft of using (3) — the expand–rotate–resum algorithm, beats, revivals — is the next-but-two article, Spectral Decomposition of Time Evolution. One honest caveat with a famous name. In infinite dimensions, "Hermitian generator" needs the self-adjoint fine print flagged at Hermitian Operators and Observables — and the rigorous statement is Stone's theorem: every strongly continuous one-parameter unitary group is e−iA^t/ℏe^{-i\hat At/\hbar} for a unique self-adjoint A^\hat A, and conversely. Self-adjointness is exactly the credential that buys a Hamiltonian its evolution operator. Proving that real atoms hold this credential was a landmark — the portrait below.

Which U for which Hamiltonian

The taxonomy every dynamics problem starts with:

regimeevolution operatorwhy
H^\hat H constantexp⁡[−iℏH^(t−t0)]\exp[-\tfrac{i}{\hbar}\hat H(t - t_0)]one generator; the product limit (2) applies verbatim
H^(t)\hat H(t), but [H^(t′),H^(t′′)]=0[\hat H(t'), \hat H(t'')] = 0exp⁡[−iℏ∫t0tH^(t′) dt′]\exp[-\tfrac{i}{\hbar}\int_{t_0}^t \hat H(t')\,dt']commuting generators — exponents add like numbers
[H^(t′),H^(t′′)]≠0[\hat H(t'), \hat H(t'')] \neq 0Texp⁡[−iℏ∫t0tH^(t′) dt′]\mathcal T\exp[-\tfrac{i}{\hbar}\int_{t_0}^t \hat H(t')\,dt']order matters at every instant; see Time-Ordered Exponential and the Dyson Series
any regime, infinitesimal1^−iℏH^(t) dt\hat 1 - \tfrac{i}{\hbar}\hat H(t)\,dtthe generator statement itself — always true

Row three is not the naive integrated exponential — that formula is wrong whenever H^\hat H fails to commute with itself across times, and the failure is quantitative (a two-piece counterexample is worked to a number in the Dyson article). The composition law survives in every regime; only time-translation invariance — dependence on t−t0t - t_0 alone — is special to constant H^\hat H.

Limiting cases

caseresultmeaning
t→t0t \to t_0U^→1^−iℏH^ (t−t0)\hat U \to \hat 1 - \tfrac{i}{\hbar}\hat H\,(t - t_0)the generator is the short-time fingerprint
eigenstate inputU^∣En⟩=e−iEnt/ℏ∣En⟩\hat U|E_n\rangle = e^{-iE_n t/\hbar}|E_n\ranglepure phase — stationary states, next article
[A^,H^]=0[\hat A, \hat H] = 0U^†A^U^=A^\hat U^\dagger\hat A\hat U = \hat A, ⟨A^⟩\langle\hat A\rangle frozenconstants of the motion — the conservation story ahead
H^→H^+c\hat H \to \hat H + cU^→e−ict/ℏ U^\hat U \to e^{-ict/\hbar}\,\hat Uglobal phase only — the energy-shift theorem of the previous article
run time backwardU^(t0,t)=U^†(t,t0)\hat U(t_0, t) = \hat U^\dagger(t, t_0)unitary evolution is reversible; nothing is lost between measurements
ℏ→0\hbar \to 0 at fixed energiesphases spin infinitely fastclassical limit lives in stationary-phase cancellations (propagator articles ahead)

Worked examples

Example 1 — spin precession, the cleanest U in physics (confidence builder). An electron sits in a uniform field B=Bz^\mathbf B = B\hat z with B=1B = 1 T. Its Hamiltonian is H^=ωS^z\hat H = \omega \hat S_z with ω=geμBB/ℏ\omega = g_e \mu_B B/\hbar. Find ω\omega, write U^(t)\hat U(t), and evolve the state ∣+x⟩|{+x}\rangle.

Step 1. The rate. ω=geμBB/ℏ=2.0023×9.274×10−24×11.055×10−34=1.761×1011 rad/s\omega = g_e\mu_B B/\hbar = \dfrac{2.0023 \times 9.274\times10^{-24} \times 1}{1.055\times10^{-34}} = \boxed{1.761\times10^{11}\ \mathrm{rad/s}} — i.e. f=ω/2π=28.02 GHzf = \omega/2\pi = \boxed{28.02\ \mathrm{GHz}}, the electron-spin-resonance frequency at one tesla.

Step 2. The operator. Since S^z=(ℏ/2)σz\hat S_z = (\hbar/2)\sigma_z: U^(t)=e−iωt σz/2\hat U(t) = e^{-i\omega t\,\sigma_z/2} — a rotation about z^\hat z by angle ωt\omega t, diagonal in the ∣±z⟩|\pm z\rangle basis with eigenphases e∓iωt/2e^{\mp i\omega t/2}. The spectral form (3) with two terms.

Step 3. Evolve ∣+x⟩=(∣+z⟩+∣−z⟩)/2|{+x}\rangle = (|{+z}\rangle + |{-z}\rangle)/\sqrt2. Each branch takes its own phase; the relative phase e−iωte^{-i\omega t} is the physics: ⟨σx⟩(t)=cos⁡ωt\langle\sigma_x\rangle(t) = \cos\omega t, ⟨σy⟩(t)=+sin⁡ωt\langle\sigma_y\rangle(t) = +\sin\omega t (machine-checked at 41 times) — the spin expectation precesses in the xyxy-plane at ω\omega, one full turn every T=2π/ω=35.7 psT = 2\pi/\omega = \boxed{35.7\ \mathrm{ps}}.

Step 4. A quantum signature for free: at t=Tt = T the expectations return, but U^(T)=e−iπσz=−1^\hat U(T) = e^{-i\pi\sigma_z} = -\hat 1 — the spinor comes back with a minus sign (verified exactly). A 2π2\pi rotation is not the identity on spin-12\frac12; the sign is physical in interferometry, and its full story belongs to the angular-momentum chapter's SU(2) article.

Example 2 — cashing the eigenphase handshake: the two-level oscillation (exam level). Completeness Relation and Resolution of Identity built e−iH^t/ℏe^{-i\hat Ht/\hbar} for the coupled two-level system of Eigenvalues and Eigenstates — ε1=1.0\varepsilon_1 = 1.0, ε2=2.0\varepsilon_2 = 2.0, coupling g=0.3g = 0.3 eV, eigenvalues E∓=0.917/2.083E_\mp = 0.917 / 2.083 eV — and reported the eigenphase beat time 2πℏ/(E+−E−)=3.552\pi\hbar/(E_+ - E_-) = 3.55 fs, promising this article would turn it into oscillation. Prepare the system in the bare state ∣1⟩|1\rangle and find the probability of finding it in ∣2⟩|2\rangle at time tt.

Step 1. Expand the initial state over eigenstates (mixing angle tan⁡2θ=2g/(ε2−ε1)⇒θ=15.5°\tan 2\theta = 2g/(\varepsilon_2 - \varepsilon_1) \Rightarrow \theta = 15.5°), apply (3), resum. The off-diagonal element of U^\hat U gives, exactly,

P1→2(t)=sin⁡22θ sin⁡2 ⁣((E+−E−) t2ℏ)P_{1\to2}(t) = \sin^2 2\theta\,\sin^2\!\Big(\frac{(E_+ - E_-)\,t}{2\hbar}\Big)

verified against the direct matrix exponential to 10−1210^{-12} across two full periods.

Step 2. Read the two factors. The amplitude sin⁡22θ=0.265\sin^2 2\theta = \boxed{0.265} is set by geometry — how far the eigenbasis is rotated from the bare basis. The clock runs at the Bohr frequency: full oscillation period =2πℏ/(E+−E−)=3.55 fs= 2\pi\hbar/(E_+ - E_-) = \boxed{3.55\ \mathrm{fs}}, the promised beat.

Step 3. The peak transfer: Pmax⁡=26.5%P_{\max} = 26.5\% at t=1.77 fst = \boxed{1.77\ \mathrm{fs}} (half the beat period), then the population flows back — forever, since evolution is unitary and reversible. Complete transfer would need θ=45°\theta = 45°: degenerate bare levels (ε1=ε2\varepsilon_1 = \varepsilon_2), the resonance condition. This formula is the Rabi oscillation of a static coupling; the driven version rides on the same algebra in the time-dependent chapter.

Step 4. Conservation check, both flavors: P1→1+P1→2=1P_{1\to1} + P_{1\to2} = 1 at every instant (unitarity), and ⟨H^⟩\langle\hat H\rangle is constant while ⟨σx⟩\langle\sigma_x\rangle oscillates — [H^,H^]=0[\hat H, \hat H] = 0 but [σx,H^]≠0[\sigma_x, \hat H] \neq 0 (both machine-checked), previewing the constants-of-the-motion criterion.

The man who proved the hydrogen atom evolves

Origins · 1 minTosio Kato and how this was discovered
TK
Tosio Kato
1917–1999 · University of Tokyo, then Berkeley

Stone's theorem says a self-adjoint Hamiltonian generates a unitary group — but for twenty years after von Neumann sharpened "Hermitian" into "self-adjoint," nobody knew whether the Hamiltonian of a real atom, Coulomb singularities and all, actually held that credential. If not, e−iH^t/ℏe^{-i\hat Ht/\hbar} for hydrogen would be formalism without a referent. Kato settled it in 1951: the Coulomb Hamiltonians of atoms and molecules are essentially self-adjoint — the singular attraction is dominated by the kinetic term (the Kato–Rellich smallness argument), so the evolution operator of every atom in the universe rigorously exists. He had proved it years earlier in wartime Tokyo, working while evacuated to the countryside with tuberculosis, unable to publish; the paper appeared in an American journal after a colleague circulated it. Kato's Perturbation Theory for Linear Operators (1966) became the standard reference for how spectra and evolution respond when Hamiltonians change; the adiabatic theorem in modern operator form is his too (1950). A quiet figure who loved botany and precise estimates, he supplied quantum dynamics with the thing physicists rarely notice needing: the proof that its central object exists. (Von Neumann's self-adjointness fine print lives at Hermitian Operators and Observables — this band is the existence-of-dynamics story only.)

History
1926
— Schrödinger's equation gives dynamics as a differential law; solving it state by state is the only tool
1927
— Weyl's group-theoretic quantum mechanics recasts evolution as a one-parameter unitary family
1932
— Stone proves the theorem pairing unitary groups with self-adjoint generators; von Neumann's Grundlagen frames dynamics as U(t)
1949
— Dyson's series organizes the time-ordered exponential for QED's interaction picture
1951
— Kato proves atomic Hamiltonians self-adjoint: the hydrogen atom's U(t) rigorously exists
1996
— Lloyd's universal quantum simulation runs on Trotterized products of short-time evolution operators

Going further

References. Sakurai & Napolitano, Modern Quantum Mechanics, 2nd ed., Sec. 2.1–2.2 · Griffiths, Introduction to Quantum Mechanics, 3rd ed., Ch. 2 · Shankar, Principles of Quantum Mechanics, 2nd ed., Ch. 4 · Cohen-Tannoudji, Diu & Laloë, Quantum Mechanics, Vol. I, Ch. III · Townsend, A Modern Approach to Quantum Mechanics, 2nd ed., Ch. 4.

Provenance. Definitions and properties (operator Schrodinger equation, unitarity, composition law, exponential and time-ordered forms) follow Sakurai & Napolitano Modern QM 2nd ed. Sec. 2.1-2.2, Shankar Ch. 4, Cohen-Tannoudji Vol. I Ch. III, Griffiths 3rd ed. Ch. 2. Verifier should: (1) confirm ∂t(U†U)=0\partial_t(U^\dagger U)=0 algebra uses only H=H†H=H^\dagger; (2) check the (1+x/N)N→ex(1+x/N)^N\to e^x operator limit requires commuting factors; (3) verify the spectral-form resummation; (4) confirm the infinite-square-well revival time Trev=4mL2/(πℏ)T_{rev}=4mL^2/(\pi\hbar) from En=n2π2ℏ2/(2mL2)E_n=n^2\pi^2\hbar^2/(2mL^2) and the Rabi population sin⁡2(ΩRt/2)\sin^2(\Omega_R t/2). All numbers are standard textbook values. For this edition all quantitative claims were machine-verified (numpy/scipy/sympy, 20 checks green): the operator-ODE and composition/inverse/unitarity identities; the product-limit error ladder (0.24/0.022/0.00220.24/0.022/0.0022 at N=10/102/103N = 10/10^2/10^3); the spectral assembly against the matrix exponential (10−1210^{-12}); the precession numbers (ω=1.761×1011\omega = 1.761\times10^{11} rad/s, f=28.02f = 28.02 GHz, T=35.7T = 35.7 ps, ⟨σx,y⟩\langle\sigma_{x,y}\rangle traces, the 2π2\pi spinor sign); and the two-level oscillation (E∓=0.917/2.083E_\mp = 0.917/2.083 eV, θ=15.5°\theta = 15.5°, amplitude 0.2650.265, beat 3.553.55 fs, peak at 1.771.77 fs). Per the chapter ownership map: the state-level equation anatomy belongs to Time-Dependent Schrodinger Equation (previous serial); the expand–rotate–resum craft, beats, and the square-well revival (Trev=4mL2/πℏT_{rev} = 4mL^2/\pi\hbar, quoted from the base) to Spectral Decomposition of Time Evolution; the time-ordered exponential and its counterexample number to Time-Ordered Exponential and the Dyson Series; the rank-1 spectral assembly consumed in equation (3) to Completeness Relation and Resolution of Identity (its 3.55 fs eigenphase handshake is cashed in Example 2); the unitary class and generator doctrine to Unitary Operators and Transformations; conservation-law development to Constants of the Motion and Conservation Laws; the propagator kernel ⟨x∣U^∣x′⟩\langle x|\hat U|x'\rangle to The Quantum Propagator.

Builds on: Time-Dependent Schrodinger Equation · Stationary States and Energy Eigenstates · The Schrodinger Picture

Leads to: Spectral Decomposition of Time Evolution · Time-Ordered Exponential and the Dyson Series · The Heisenberg Picture · Equivalence of the Dynamical Pictures · Constants of the Motion and Conservation Laws

Is U(t,t₀) something I measure, or just bookkeeping?

Its matrix elements are measured constantly: ∣⟨f∣U^(t)∣i⟩∣2|\langle f|\hat U(t)|i\rangle|^2 is the transition probability from ∣i⟩|i\rangle to ∣f⟩|f\rangle — Example 2's 26.5% is exactly ∣U21∣2|U_{21}|^2 at its peak, and quantum-computing labs reconstruct entire gate unitaries element by element (process tomography). What you never measure is U^\hat U's overall phase — consistent with the energy-shift theorem: adding a constant to H^\hat H re-phases U^\hat U without touching any probability.

Why does U sit to the LEFT in the composition law, with the earliest interval rightmost?

Operators act on kets from the left: in U^(t2,t1)U^(t1,t0)∣ψ(t0)⟩\hat U(t_2,t_1)\hat U(t_1,t_0)|\psi(t_0)\rangle the state meets the t0→t1t_0 \to t_1 segment first, as chronology demands. Misordering is harmless for constant H^\hat H but wrong the moment the Hamiltonian's direction turns — the same order-sensitivity that forces time-ordering in the Dyson series. Reading products right-to-left is the most protective habit in this chapter.

Does U exist when the Hamiltonian depends on time? The exponential formula seems to die.

The formula (2) dies; the operator survives. The defining contract never referenced constancy of H^\hat H. For self-commuting time dependence the exponent just integrates; in general the solution is the time-ordered exponential, built two articles ahead. Losing the simple formula costs only time-translation invariance (U^\hat U now depends on both endpoints, not the difference).

Is the Hamiltonian the only operator that could generate time evolution?

Unitarity forces some Hermitian generator (Derivation A, Step 2); which observable it is, is physics. The identification with energy comes from correspondence: classically the Hamiltonian generates time translation through Poisson brackets, and the momentum-generates-displacement analogy fixes the pairing and the ℏ\hbar. The deep version — time-translation symmetry pairs with energy conservation — is Noether's logic, developed at this chapter's symmetry articles.

If unitary evolution is reversible, where does irreversibility in the real world come from?

Not from U^\hat U — a closed system's evolution can always be undone by U^†\hat U^\dagger (spin echoes do exactly this, rewinding dephased magnetization). Irreversibility enters two ways: measurement (the projection postulate is not unitary — the collapse story of Measurement and Wavefunction Collapse), and openness — a system entangling with an unmonitored environment has no unitary of its own. Inside this chapter, everything is closed and reversible by construction.

My textbook writes U(t) = exp(−iHt/ħ) with t₀ = 0 everywhere. When does the t₀ matter?

For constant H^\hat H, never — only elapsed time enters, and setting t0=0t_0 = 0 is free. For time-dependent Hamiltonians the family genuinely depends on both arguments: a pulse applied at 9:00 and the same pulse applied at 9:05 produce different evolutions relative to the lab clock, and U^(t2,t1)≠U^(t2−t1,0)\hat U(t_2, t_1) \neq \hat U(t_2 - t_1, 0) in general. Keeping both slots is the notational tell that you are ready for driven systems.

Before this

Where this leads