Skip to content

Coherent-State Dynamics

A coherent state of the harmonic oscillator is a minimum-uncertainty Gaussian whose center follows the classical oscillator orbit without changing shape. This exact stability under harmonic evolution makes coherent states a model of classical-like quantum dynamics, but not classical states: they remain nonorthogonal, retain irreducible fluctuations, and can leave the coherent-state family under squeezing or nonlinear interactions.

Coherent States owns the Hilbert-space construction, number-basis expansion, and displacement-operator identities. Coherent States in Phase Space owns the Wigner-function geometry. This page owns the time evolution and the conditions under which coherent dynamics is preserved or lost. Coherent-State Semiclassics Preview compares this exact dynamics with projected one-packet motion, coherent-state saddle propagation, and frozen-Gaussian trajectory sums.

When the chosen coherent or Gaussian family is not exactly invariant, the Time-Dependent Variational Principle supplies the canonical tangent-space projection and residual diagnostics.

For an oscillator with

H0=ℏω(a†a+12),H_0 = \hbar\omega \left( a^\dagger a+\frac{1}{2} \right),

a coherent state satisfies

a∣α⟩=α∣α⟩.a\lvert\alpha\rangle = \alpha\lvert\alpha\rangle.

Equivalently,

∣α⟩=D(α)∣0⟩,\lvert\alpha\rangle = D(\alpha)\lvert0\rangle,

where

D(α)=exp⁡(αa†−α∗a).D(\alpha) = \exp\left( \alpha a^\dagger-\alpha^*a \right).

The coherent label is related to the first moments by

xc=2 ℓ Re⁡α,x_c = \sqrt{2}\,\ell\, \operatorname{Re}\alpha,

and

pc=2ℏℓIm⁡α,p_c = \frac{\sqrt{2}\hbar}{\ell} \operatorname{Im}\alpha,

with oscillator length

ℓ=ℏmω.\ell = \sqrt{\frac{\hbar}{m\omega}}.

The label α\alpha is dimensionless. It packages a position and momentum into one complex phase-space coordinate adapted to the oscillator.

Let

U0(t)=e−iH0t/ℏ.U_0(t) = e^{-iH_0t/\hbar}.

The annihilation operator evolves as

U0†(t)aU0(t)=ae−iωt.U_0^\dagger(t)aU_0(t) = ae^{-i\omega t}.

Act on the evolved state with aa:

aU0(t)∣α⟩=U0(t)[U0†(t)aU0(t)]∣α⟩=αe−iωtU0(t)∣α⟩.\begin{aligned} aU_0(t)\lvert\alpha\rangle &= U_0(t) \left[ U_0^\dagger(t)aU_0(t) \right] \lvert\alpha\rangle \\ &= \alpha e^{-i\omega t} U_0(t)\lvert\alpha\rangle. \end{aligned}

Therefore U0(t)∣α⟩U_0(t)\lvert\alpha\rangle is another coherent state with label αe−iωt\alpha e^{-i\omega t}. Including the zero-point phase,

U0(t)∣α⟩=e−iωt/2∣αe−iωt⟩.U_0(t)\lvert\alpha\rangle = e^{-i\omega t/2} \lvert\alpha e^{-i\omega t}\rangle.

The global phase does not affect ordinary expectation values, but it matters when this branch interferes with another amplitude.

The coherent-state family is therefore invariant under harmonic evolution:

α(t)=α0e−iωt.\alpha(t) = \alpha_0e^{-i\omega t}.

Writing

α0=x02ℓ+iℓp02ℏ,\alpha_0 = \frac{x_0}{\sqrt{2}\ell} + i\frac{\ell p_0}{\sqrt{2}\hbar},

the first moments evolve as

xc(t)=x0cos⁡ωt+p0mωsin⁡ωt,x_c(t) = x_0\cos\omega t + \frac{p_0}{m\omega}\sin\omega t,

and

pc(t)=p0cos⁡ωt−mωx0sin⁡ωt.p_c(t) = p_0\cos\omega t - m\omega x_0\sin\omega t.

These are exactly the classical harmonic-oscillator equations. In scaled phase-space coordinates, the label rotates on a circle; in dimensional (x,p)(x,p) coordinates, the center follows the classical energy ellipse.

The mean energy is

⟨H0⟩α=ℏω(∣α∣2+12).\langle H_0\rangle_\alpha = \hbar\omega \left( \lvert\alpha\rvert^2+\frac{1}{2} \right).

The classical contribution is associated with the action scale

J∼ℏ∣α∣2.J \sim \hbar\lvert\alpha\rvert^2.

Large ∣α∣\lvert\alpha\rvert therefore corresponds to large action compared with ℏ\hbar, while the zero-point contribution remains.

For every coherent state,

(Δx)2=ℏ2mω,(\Delta x)^2 = \frac{\hbar}{2m\omega},

and

(Δp)2=mℏω2.(\Delta p)^2 = \frac{m\hbar\omega}{2}.

The symmetrized covariance is zero in the natural oscillator axes:

Cxp=0.C_{xp}=0.

Hence

Δx Δp=ℏ2.\Delta x\,\Delta p = \frac{\hbar}{2}.

Under harmonic evolution these variances remain constant. The Wigner Gaussian translates around the classical orbit without shearing or changing area.

If the classical position amplitude is

Ax=2 ℓ∣α∣,A_x = \sqrt{2}\,\ell\lvert\alpha\rvert,

then

ΔxAx=12∣α∣.\frac{\Delta x}{A_x} = \frac{1}{2\lvert\alpha\rvert}.

Thus the packet becomes narrow relative to its orbit for ∣α∣≫1\lvert\alpha\rvert\gg1. This is a relative classicality statement. The absolute quantum width does not vanish at fixed mm and ω\omega.

Number Statistics and Relative Fluctuations

Section titled “Number Statistics and Relative Fluctuations”

The occupation-number distribution is Poissonian:

Pα(n)=e−∣α∣2∣α∣2nn!.P_\alpha(n) = e^{-\lvert\alpha\rvert^2} \frac{\lvert\alpha\rvert^{2n}}{n!}.

Therefore

⟨n⟩=∣α∣2,\langle n\rangle = \lvert\alpha\rvert^2,

and

Δn=∣α∣.\Delta n = \lvert\alpha\rvert.

The relative number fluctuation is

Δn⟨n⟩=1∣α∣.\frac{\Delta n}{\langle n\rangle} = \frac{1}{\lvert\alpha\rvert}.

Large-amplitude coherent states have small relative number fluctuations even though their absolute fluctuations grow. This scaling is one reason classical fields can emerge from highly occupied coherent modes without eliminating quantum shot noise.

Coherent states are not orthogonal:

⟨α∣β⟩=exp⁡[−∣α∣22−∣β∣22+α∗β].\langle\alpha\vert\beta\rangle = \exp\left[ -\frac{\lvert\alpha\rvert^2}{2} -\frac{\lvert\beta\rvert^2}{2} +\alpha^*\beta \right].

The survival probability under harmonic evolution is

∣⟨α∣U0(t)∣α⟩∣2=exp⁡[−4∣α∣2sin⁡2(ωt2)].\begin{aligned} \left\lvert \langle\alpha\vert U_0(t)\vert\alpha\rangle \right\rvert^2 &= \exp\left[ -4\lvert\alpha\rvert^2 \sin^2\left( \frac{\omega t}{2} \right) \right]. \end{aligned}

For large ∣α∣\lvert\alpha\rvert, the state rapidly becomes nearly orthogonal to its initial state as the packet center moves away, even though it remains perfectly coherent. After one oscillator period,

T=2πω,T=\frac{2\pi}{\omega},

the label returns and the survival probability is one. Coherent evolution combines classical orbital recurrence with quantum overlap geometry.

Add a prescribed linear drive:

H(t)=ℏω(a†a+12)+ℏ[η(t)a†+η∗(t)a],\begin{aligned} H(t) &= \hbar\omega \left( a^\dagger a+\frac{1}{2} \right) \\ &\quad+ \hbar \left[ \eta(t)a^\dagger +\eta^*(t)a \right], \end{aligned}

where η(t)\eta(t) has units of angular frequency. A coherent initial state remains coherent up to a global phase:

∣ψ(t)⟩=eiϕ(t)∣α(t)⟩.\lvert\psi(t)\rangle = e^{i\phi(t)} \lvert\alpha(t)\rangle.

The label obeys the classical linear equation

iα˙=ωα+η(t),i\dot\alpha = \omega\alpha+\eta(t),

or

α˙=−iωα−iη(t).\dot\alpha = -i\omega\alpha-i\eta(t).

With initial time t0t_0,

α(t)=e−iω(t−t0)α(t0)−i∫t0tds e−iω(t−s)η(s).\begin{aligned} \alpha(t) &= e^{-i\omega(t-t_0)} \alpha(t_0) \\ &\quad- i \int_{t_0}^{t} ds\, e^{-i\omega(t-s)} \eta(s). \end{aligned}

The first term is free rotation. The integral is the displacement produced by the drive. The phase ϕ(t)\phi(t) includes the classical action and ordering information but does not change the coherent label.

For a position force −F(t)x-F(t)x, one may identify

η(t)=−F(t)ℓ2ℏ\eta(t) = -\frac{F(t)\ell} {\sqrt{2}\hbar}

in the convention used above. The resulting first moments obey the forced classical oscillator equation

mx¨c+mω2xc=F(t).m\ddot x_c + m\omega^2x_c = F(t).

Linear forcing therefore moves the Gaussian center without changing its intrinsic covariance.

The coherent-state family is preserved by Hamiltonians that generate phase rotations and displacements:

H(t)=ℏω(t)a†a+ℏ[η(t)a†+η∗(t)a]+c(t)I,H(t) = \hbar\omega(t)a^\dagger a + \hbar \left[ \eta(t)a^\dagger +\eta^*(t)a \right] + c(t)I,

with suitable time ordering. The state remains coherent, although its label and global phase change.

A general quadratic Hamiltonian can also contain

ℏ2[ζ(t)a†2+ζ∗(t)a2].\frac{\hbar}{2} \left[ \zeta(t)a^{\dagger 2} +\zeta^*(t)a^2 \right].

These parametric terms mix aa and a†a^\dagger. Gaussianity is preserved, but a coherent state generally becomes a squeezed coherent state. The coherent family is not invariant.

An anharmonic term such as

HKerr∝(a†a)2H_{\rm Kerr} \propto (a^\dagger a)^2

assigns nonlinear phases to number states. An initial coherent state then shears in phase space and can develop collapse, revival, and cat-like superpositions. Coherence in the technical Glauber sense is not preserved by generic nonlinear dynamics.

Coherent states combine several classical-like features:

  • their first moments follow the classical oscillator equations;
  • their covariance remains minimal and fixed under harmonic motion;
  • large amplitudes make relative position and number fluctuations small;
  • a linear drive transports them according to the forced classical equation;
  • normally ordered field correlations factorize for ideal coherent radiation.

They remain quantum because

  • ΔxΔp=ℏ/2\Delta x\Delta p=\hbar/2, not zero;
  • different coherent states have nonzero overlap;
  • number fluctuations and shot noise remain;
  • coherent superpositions can interfere;
  • nonlinear dynamics can generate squeezing, entanglement, and non-Gaussian states.

The correct phrase is classical-like quantum state, not classical state.

Canonical coherent states resolve the identity:

I=∫Cd2απ∣α⟩⟨α∣.I = \int_{\mathbb C} \frac{d^2\alpha}{\pi} \lvert\alpha\rangle \langle\alpha\rvert.

Inserting this overcomplete identity between short-time evolution factors produces a path integral over complex labels. A common continuum action is

S[α∗,α]=∫titfdt[iℏ2(α∗α˙−α˙∗α)−Hcs(α∗,α)],\begin{aligned} S[\alpha^*,\alpha] &= \int_{t_i}^{t_f}dt \Bigg[ \frac{i\hbar}{2} \left( \alpha^*\dot\alpha -\dot\alpha^*\alpha \right) \\ &\qquad\qquad- H_{\rm cs}(\alpha^*,\alpha) \Bigg], \end{aligned}

together with endpoint terms determined by the boundary convention. The Hamiltonian symbol HcsH_{\rm cs} depends on the operator ordering and time-slicing prescription.

Stationary variation gives

iℏα˙=∂Hcs∂α∗,i\hbar\dot\alpha = \frac{\partial H_{\rm cs}} {\partial\alpha^*},

and

−iℏα˙∗=∂Hcs∂α.-i\hbar\dot\alpha^* = \frac{\partial H_{\rm cs}} {\partial\alpha}.

For the harmonic oscillator these are the classical complex phase-space equations. The preview is useful because it explains why coherent-state labels become field variables in many-body and QFT path integrals. It is not a license to ignore boundary terms, overcompleteness, or symbol-ordering corrections.

Why Path Integrals? supplies the general amplitude and action logic. Coherent-State Path Integrals Preview fixes the thermal trace, finite-slice discretization, Grassmann boundary sign, and operator symbol before assigning a continuum action.

A free bosonic field decomposes into oscillator modes. A multimode coherent state satisfies

ak∣{α}⟩=αk∣{α}⟩a_{\mathbf k} \lvert\{\alpha\}\rangle = \alpha_{\mathbf k} \lvert\{\alpha\}\rangle

for every retained mode. Under a free field Hamiltonian,

αk(t)=αk(0)e−iωkt.\alpha_{\mathbf k}(t) = \alpha_{\mathbf k}(0) e^{-i\omega_{\mathbf k}t}.

Field expectation values built from these amplitudes obey the corresponding classical linear field equations. This is the field analogue of the oscillator center following its classical orbit.

The state is still quantum. Vacuum and shot-noise fluctuations remain, coherent states are not number eigenstates, and interacting evolution generally entangles modes and drives the state away from a simple product of coherent states. Gauge fields also require physical-mode and constraint care.

The oscillator-to-field dictionary continues in Harmonic Oscillator to Fields.

  • Treating α\alpha as a dimensional position rather than a complex oscillator coordinate.
  • Dropping the zero-point phase when comparing interfering amplitudes.
  • Saying minimum uncertainty means zero uncertainty.
  • Assuming every minimum-uncertainty state is an unsqueezed coherent state.
  • Calling Poissonian number fluctuations noise-free.
  • Assuming every quadratic Hamiltonian preserves coherent states; parametric terms generally squeeze them.
  • Assuming a coherent state remains coherent under an anharmonic interaction.
  • Writing a coherent-state path integral without specifying endpoint and ordering conventions.
  • Treating a large field expectation value as proof that quantum fluctuations and entanglement are absent.
  • E. Schrödinger, “Der stetige Übergang von der Mikro- zur Makromechanik,” Die Naturwissenschaften 14, 664–666, 1926.
  • R. J. Glauber, “Coherent and incoherent states of the radiation field,” Physical Review 131, 2766–2788, 1963, doi:10.1103/PhysRev.131.2766.
  • J. R. Klauder and B.-S. Skagerstam, eds., Coherent States: Applications in Physics and Mathematical Physics, World Scientific, 1985.
  • A. Perelomov, Generalized Coherent States and Their Applications, Springer, 1986.
  • G. A. Hagedorn, “Semiclassical quantum mechanics. I. The ℏ→0\hbar\to0 limit for coherent states,” Communications in Mathematical Physics 71, 77–93, 1980.
  • W. P. Schleich, Quantum Optics in Phase Space, Wiley-VCH, 2001.
  1. Prove that harmonic evolution maps a coherent state to another coherent state.
Solution

Use

U0†(t)aU0(t)=ae−iωt.U_0^\dagger(t)aU_0(t) = ae^{-i\omega t}.

Then

aU0(t)∣α⟩=U0(t)[U0†(t)aU0(t)]∣α⟩=αe−iωtU0(t)∣α⟩.\begin{aligned} aU_0(t)\lvert\alpha\rangle &= U_0(t) \left[ U_0^\dagger(t)aU_0(t) \right] \lvert\alpha\rangle \\ &= \alpha e^{-i\omega t} U_0(t)\lvert\alpha\rangle. \end{aligned}

Thus the evolved state is an eigenstate of aa with eigenvalue αe−iωt\alpha e^{-i\omega t}. Comparing the vacuum component or using the number expansion gives

U0(t)∣α⟩=e−iωt/2∣αe−iωt⟩.U_0(t)\lvert\alpha\rangle = e^{-i\omega t/2} \lvert\alpha e^{-i\omega t}\rangle.
  1. Derive the classical first-moment trajectory from α(t)=α0e−iωt\alpha(t)=\alpha_0e^{-i\omega t}.
Solution

Write

α0=x02ℓ+iℓp02ℏ.\alpha_0 = \frac{x_0}{\sqrt{2}\ell} + i\frac{\ell p_0}{\sqrt{2}\hbar}.

Using

xc(t)=2ℓ Re⁡α(t)x_c(t) = \sqrt{2}\ell\, \operatorname{Re}\alpha(t)

and

pc(t)=2ℏℓIm⁡α(t),p_c(t) = \frac{\sqrt{2}\hbar}{\ell} \operatorname{Im}\alpha(t),

expand e−iωt=cos⁡ωt−isin⁡ωte^{-i\omega t}=\cos\omega t-i\sin\omega t. This gives

xc(t)=x0cos⁡ωt+p0mωsin⁡ωt,x_c(t) = x_0\cos\omega t + \frac{p_0}{m\omega}\sin\omega t,

and

pc(t)=p0cos⁡ωt−mωx0sin⁡ωt.p_c(t) = p_0\cos\omega t - m\omega x_0\sin\omega t.
  1. Compute the harmonic survival probability.
Solution

The global phase cancels from the modulus, so

∣⟨α∣U0(t)∣α⟩∣2=∣⟨α∣αe−iωt⟩∣2=exp⁡[−∣α−αe−iωt∣2].\begin{aligned} \left\lvert \langle\alpha\vert U_0(t)\vert\alpha\rangle \right\rvert^2 &= \left\lvert \langle\alpha \vert \alpha e^{-i\omega t} \rangle \right\rvert^2 \\ &= \exp\left[ -\left\lvert \alpha-\alpha e^{-i\omega t} \right\rvert^2 \right]. \end{aligned}

Now

∣1−e−iωt∣2=4sin⁡2(ωt2).\left\lvert 1-e^{-i\omega t} \right\rvert^2 = 4\sin^2\left( \frac{\omega t}{2} \right).

Therefore

∣⟨α∣U0(t)∣α⟩∣2=exp⁡[−4∣α∣2sin⁡2(ωt2)].\left\lvert \langle\alpha\vert U_0(t)\vert\alpha\rangle \right\rvert^2 = \exp\left[ -4\lvert\alpha\rvert^2 \sin^2\left( \frac{\omega t}{2} \right) \right].
  1. A resonant drive has η(t)=η0e−iωt\eta(t)=\eta_0e^{-i\omega t} for 0≤t≤T0\leq t\leq T. Find α(T)\alpha(T).
Solution

With t0=0t_0=0,

α(T)=e−iωTα0−i∫0Tds e−iω(T−s)η0e−iωs.\begin{aligned} \alpha(T) &= e^{-i\omega T}\alpha_0 \\ &\quad- i \int_0^Tds\, e^{-i\omega(T-s)} \eta_0e^{-i\omega s}. \end{aligned}

The integrand is η0e−iωT\eta_0e^{-i\omega T}, so

α(T)=e−iωT(α0−iη0T).\alpha(T) = e^{-i\omega T} \left( \alpha_0-i\eta_0T \right).

In the frame rotating at ω\omega, the coherent amplitude is displaced linearly with time.

  1. Show that the relative number fluctuation vanishes for large coherent amplitude.
Solution

For a Poisson distribution,

⟨n⟩=∣α∣2,\langle n\rangle = \lvert\alpha\rvert^2,

and

(Δn)2=∣α∣2.(\Delta n)^2 = \lvert\alpha\rvert^2.

Hence

Δn⟨n⟩=∣α∣∣α∣2=1∣α∣.\frac{\Delta n}{\langle n\rangle} = \frac{\lvert\alpha\rvert} {\lvert\alpha\rvert^2} = \frac{1}{\lvert\alpha\rvert}.

This tends to zero as ∣α∣→∞\lvert\alpha\rvert\to\infty, although the absolute fluctuation Δn=∣α∣\Delta n=\lvert\alpha\rvert grows.

  1. Why does a parametric term proportional to a†2+a2a^{\dagger2}+a^2 preserve Gaussianity but not the coherent-state family?
Solution

A quadratic Hamiltonian implements a linear symplectic transformation of the quadratures, so a Gaussian Wigner function remains Gaussian. A parametric term mixes aa and a†a^\dagger through a Bogoliubov transformation. That transformation changes the covariance matrix and generally creates squeezing. An ordinary coherent state is a displaced ground-state Gaussian with the unsqueezed ground covariance, so the evolved squeezed Gaussian is not an eigenstate of the original aa operator and is not an ordinary coherent state.