Quantum dynamics describes how a quantum state evolves in time. It is relevant whenever you need to predict the time-dependence of a system: the oscillation of a qubit in a magnetic field, the scattering of two particles, the tunneling of an electron through a barrier, or the propagation of a wave packet. The dynamics is governed by the Hamiltonian, which encodes the interactions (see Interactions).
Definition. The state ∣ψ(t)⟩ evolves according to iℏ∂t∣ψ(t)⟩=H^∣ψ(t)⟩, where H^ is the Hamiltonian operator. The formal solution is ∣ψ(t)⟩=U^(t,t0)∣ψ(t0)⟩ with U^(t,t0)=Texp(−ℏi∫t0tH^(t′)dt′). For a time-independent Hamiltonian: U^(t)=e−iH^t/ℏ.
Details
From the classical Hamiltonian to the quantum one §
In classical mechanics, the state evolves via Hamilton’s equations: q˙=∂H/∂p, p˙=−∂H/∂q. In quantum mechanics, the state is a vector in Hilbert space (see Quantum states) and evolves via the Schrödinger equation. The classical Hamiltonian function H(q,p) becomes an operator H^ by the correspondence q→q^, p→p^=−iℏ∂/∂x.
Step 1: Time-independent case.
If H^ does not depend on time, the evolution operator is U^(t)=e−iH^t/ℏ. If the initial state is an energy eigenstate H^∣n⟩=En∣n⟩, then ∣ψ(t)⟩=e−iEnt/ℏ∣n⟩ — the state picks up only a global phase, and all measurable probabilities are constant. Superpositions of energy eigenstates, however, oscillate: ∣ψ(t)⟩=c1e−iE1t/ℏ∣1⟩+c2e−iE2t/ℏ∣2⟩ has a time-dependent relative phase e−i(E2−E1)t/ℏ that produces oscillating probabilities.
Step 2: Time-dependent case.
When H^ depends on time (e.g., a time-varying electromagnetic field), the evolution operator involves a time-ordering T. There is no simple exponential form, and one must use perturbation theory or numerical methods.
Example: a two-level system in a resonant field.
A qubit with energy splitting ℏω0 is driven by an oscillating field of frequency ω≈ω0. In the rotating-wave approximation, the effective Hamiltonian in the interaction picture is H^eff=2ℏ(ΔΩΩ−Δ) where Δ=ω−ω0 is the detuning and Ω is the Rabi frequency. Starting in ∣0⟩, the probability of finding the system in ∣1⟩ oscillates as P1(t)=Ω2+Δ2Ω2sin2(2Ω2+Δ2t). At resonance (Δ=0), full population transfer occurs at t=π/Ω (a π-pulse).
Example problem.
Problem
A spin-½ particle in a uniform magnetic field B=B0z^ has Hamiltonian H^=−γB0S^z=−2ℏω0σz where ω0=γB0. If the initial state is ∣↑x⟩=(∣↑z⟩+∣↓z⟩)/2, find ∣ψ(t)⟩ and the probability of measuring Sz=+ℏ/2.
The probability is constant at 1/2 — the spin precesses around z at frequency ω0, but the z-component remains equally likely to be up or down. However, ⟨Sx⟩(t)=2ℏcos(ω0t) oscillates, showing the precession.
Definition. For H^=H^0+V^(t), the transition amplitude from state ∣i⟩ to state ∣f⟩ (first order) is cf(1)(t)=−ℏi∫t0t⟨f∣V^(t′)∣i⟩eiωfit′dt′ where ωfi=(Ef−Ei)/ℏ. Fermi’s golden rule gives the transition rate: Γi→f=ℏ2π∣⟨f∣V^∣i⟩∣2ρ(Ef).”
Many physical situations involve a small time-dependent perturbation: an atom in a weak electromagnetic wave, a particle encountering a short potential barrier, a qubit subjected to a control pulse. Exact solutions are rare; perturbation theory expands in powers of V^.
Step 1: First-order amplitude.
Assume the system starts in an eigenstate ∣i⟩ of H^0 at t=−∞. The first-order correction to the amplitude of finding it in state ∣f⟩ at time t is an integral of the matrix element of V^ weighted by the phase factor eiωfit′. The integral is large when ωfi≈0 (energy conservation) — this is the resonance condition.
Step 2: Fermi’s golden rule.
For a perturbation that is on for a long time T, the transition probability grows as T2sinc2(ωfiT/2). The rate (probability per unit time) in the limit T→∞ becomes proportional to 2πδ(Ef−Ei), which for a continuum of final states gives the density of states ρ(Ef). This is Fermi’s golden rule.
Example: spontaneous emission.
An excited atom coupled to the electromagnetic vacuum (a time-dependent perturbation from the quantized field) decays with rate Γ=3πϵ0ℏc3ω03∣d∣2 where d is the dipole matrix element. This is a first-order process in the field coupling.
Example problem.
Problem
A particle in a 1D infinite well of width a is in the ground state ψ1(x)=2/asin(πx/a). A time-dependent perturbation V^(t)=V0xe−t2/τ2 is applied. Write the first-order amplitude for transition to the first excited state ψ2(x)=2/asin(2πx/a).
Definition. Scattering describes the asymptotic behavior of a wave packet that approaches, interacts with, and recedes from a potential. The S-matrix S^ maps incoming to outgoing states: ∣ψout⟩=S^∣ψin⟩. The differential cross-section is dσ/dΩ=∣f(θ,ϕ)∣2 where f is the scattering amplitude.
A scattering experiment sends a particle toward a target and measures where it ends up. Quantum mechanically, the incoming state is a wave packet (or, in the idealized limit, a plane wave). The outgoing state consists of an undistorted forward component plus a spherical outgoing wave f(θ,ϕ)eikr/r.
Step 1: The scattering amplitude.
For a spherically symmetric potential V(r), the scattering amplitude in the partial-wave expansion is:
f(θ)=k1∑ℓ=0∞(2ℓ+1)eiδℓsinδℓPℓ(cosθ)
where δℓ are the phase shifts and Pℓ are Legendre polynomials. Each partial wave ℓ contributes a term weighted by (2ℓ+1) and the phase shift δℓ.
Step 2: Low-energy (s-wave) scattering.
At low energies, only ℓ=0 contributes. The cross-section becomes isotropic: σ=4πa2 where a=−tanδ0/k is the scattering length. This single parameter characterizes low-energy scattering and is central to ultracold atom physics and neutron scattering.
Connection to the Standard Model.
In relativistic quantum field theory (see Quantum field theory), scattering is described by Feynman diagrams and the S-matrix elements are computed from loop integrals. The non-relativistic scattering theory above is the low-energy limit of the full QFT calculation.
Example problem.
Problem
A neutron of wavelength λ=1 Å scatters from a small spherical potential well of depth V0=100 meV and radius R=1 Å. In the s-wave approximation, estimate the scattering length if kR≪1.
Solution
k=2π/λ=2π/10−10=6.28×1010 m⁻¹. kR=6.28×1010×10−10=6.28. This is not ≪1, so the s-wave approximation is marginal. For a true low-energy estimate, we would need λ≫R.
For a hard sphere of radius R (infinite potential), a=R and σ=4πR2=4π(10−10)2≈1.3×10−19 m² = 130 barns.
For a finite well, the scattering length depends on the bound-state condition: if the well supports a virtual state near zero energy, a can be much larger than R (resonant scattering).
Definition. The path integral (Feynman) formulation expresses the propagator as a sum over all paths: K(xf,tf;xi,ti)=∫D[x(t)]exp(ℏi∫titfL(x,x˙)dt). Each path contributes a phase proportional to the classical action. In the limit ℏ→0 (or large action), the path of stationary action (the classical trajectory) dominates — this is the correspondence principle.
Alongside matrix mechanics (Heisenberg) and wave mechanics (Schrödinger), the path integral (Feynman, 1948) provides a conceptually distinct formulation. Instead of evolving a state vector, one sums over all possible histories weighted by eiS/ℏ.
Step 1: The propagator as a sum over paths.
The amplitude for a particle to go from (xi,ti) to (xf,tf) is:
where ϵ=(tf−ti)/N. Each term in the sum is the Lagrangian evaluated along a discrete path.
Step 2: The classical limit.
When the action S≫ℏ, paths near the classical trajectory (where δS=0) contribute constructively (their phases align), while paths far from the classical one cancel by destructive interference. The result is that the particle follows the classical trajectory — the correspondence principle connecting to Mechanics.
Path integrals are not merely a calculational tool for QFT. They provide:
A natural framework for quantum gravity (sum over geometries).
A direct connection between quantum and statistical mechanics (Wick rotation t→−iℏβ turns the path integral into a partition function; see Thermodynamics).
A geometric picture of quantum interference (the double-slit experiment is a sum over paths through both slits).
Example: the free-particle propagator.
For V=0, the path integral can be evaluated exactly. The result is:
This is a Gaussian in (xf−xi) with a complex phase — the quantum analogue of the classical free-particle trajectory.
Example problem.
Problem
Using the path integral, show that for a free particle, the propagator at tf−ti=T satisfies the composition law K(xf,T1+T2;xi)=∫dxK(xf,T1;x)K(x,T2;xi).
Complete the square in the exponent. The coefficient of x2 is −im/(2ℏ)(1/T1+1/T2)=−im/(2ℏ)⋅(T1+T2)/(T1T2). The linear term gives the shift. After Gaussian integration, the result is:
Problem. A qubit with Hamiltonian H^=2ℏωσz is initially in state ∣0⟩=∣↑z⟩. At time t=0, a sudden perturbation V^=ℏΩσx is applied for a duration τ=π/(2Ω) (a π/2-pulse). Find the state after the pulse (to first order in Ω/ω, assuming Ω≪ω).
Solution
During the pulse, the total Hamiltonian is H^+V^=2ℏωσz+ℏΩσx. To first order in Ω/ω, the dominant evolution is from H^: ∣ψ(t)⟩≈e−iωt/2∣0⟩. The perturbation gives a first-order correction:
∣ψ(τ)⟩≈e−iωτ/2[∣0⟩−iΩ∫0τdteiωtσx∣0⟩]
σx∣0⟩=∣1⟩, so:
∣ψ(τ)⟩≈e−iωτ/2[∣0⟩−iΩeiωτ/2iωeiωτ−1∣1⟩]
For Ω≪ω and τ=π/(2Ω), the integral ∫0τeiωtdt=(eiωτ−1)/(iω)≈−1/(iω) (since ωτ≫1 and the oscillatory integral is small). So the correction is of order Ω/ω≪1.
To leading order, the state remains approximately ∣0⟩ (up to a phase). The π/2-pulse would fully transfer population only if ω≈0 (resonant drive) — here the large detuning ω≫Ω suppresses the transition. This illustrates the importance of the resonance condition in time-dependent perturbation theory.