Skip to content

Ehrenfest theorem relates the time evolution of expectation values to classical-looking equations of motion. For

H=p22m+V(x),H=\frac{p^2}{2m}+V(x),

it gives

ddt⟨x⟩=⟨p⟩m,ddt⟨p⟩=−⟨V′(x)⟩.\frac{d}{dt}\langle x\rangle =\frac{\langle p\rangle}{m}, \qquad \frac{d}{dt}\langle p\rangle =-\langle V'(x)\rangle.

For an operator that may depend explicitly on time, the general identity is

ddt⟨A⟩=iℏ⟨[H,A]⟩+⟨∂A∂t⟩.\frac{d}{dt}\langle A\rangle =\frac{i}{\hbar}\langle[H,A]\rangle +\left\langle\frac{\partial A}{\partial t}\right\rangle.

It assumes that the state and the relevant operator products lie in domains for which the differentiations and integrations by parts are legitimate. This qualification matters for unbounded operators such as xx and pp.

From the Heisenberg equation,

dAHdt=iℏ[H,AH]\frac{dA_H}{dt} =\frac{i}{\hbar}[H,A_H]

for operators with no explicit time dependence. Taking expectation values gives

ddt⟨A⟩=iℏ⟨[H,A]⟩.\frac{d}{dt}\langle A\rangle = \frac{i}{\hbar}\langle[H,A]\rangle.

For A=xA=x and H=p2/(2m)+V(x)H=p^2/(2m)+V(x),

ddt⟨x⟩=⟨p⟩m.\frac{d}{dt}\langle x\rangle =\frac{\langle p\rangle}{m}.

For A=pA=p,

ddt⟨p⟩=−⟨V′(x)⟩.\frac{d}{dt}\langle p\rangle =-\langle V'(x)\rangle.

Classical Newtonian motion would use

md2dt2xcl=−V′(xcl).m\frac{d^2}{dt^2}x_{\rm cl} =-V'(x_{\rm cl}).

Ehrenfest theorem gives instead

md2dt2⟨x⟩=−⟨V′(x)⟩.m\frac{d^2}{dt^2}\langle x\rangle =-\langle V'(x)\rangle.

These agree with the classical equation for ⟨x⟩\langle x\rangle only when

⟨V′(x)⟩≈V′(⟨x⟩).\langle V'(x)\rangle \approx V'(\langle x\rangle).

This approximation is exact for potentials at most quadratic in xx, and often useful for narrow wave packets in slowly varying potentials.

To quantify the approximation, put

q=⟨x⟩,δx=x−q,σx2=⟨(δx)2⟩.q=\langle x\rangle, \qquad \delta x=x-q, \qquad \sigma_x^2=\langle(\delta x)^2\rangle.

Expanding the force about qq gives

⟨V′(x)⟩=V′(q)+12V′′′(q)σx2+16V(4)(q)⟨(δx)3⟩+⋯ .\langle V'(x)\rangle =V'(q) +\frac12V'''(q)\sigma_x^2 +\frac16V^{(4)}(q)\langle(\delta x)^3\rangle +\cdots.

This moment expansion requires a smooth force over the appreciable support of the packet and controlled higher moments; a small variance alone need not bound rare tails. The first displayed correction involves the packet width and the third derivative of the potential. For a quartic potential V(x)=λx4/4V(x)=\lambda x^4/4,

⟨V′(x)⟩=λ(q3+3qσx2+⟨(δx)3⟩).\langle V'(x)\rangle =\lambda\left(q^3+3q\sigma_x^2 +\langle(\delta x)^3\rangle\right).

Even a symmetric packet therefore acquires the finite-width correction 3λqσx23\lambda q\sigma_x^2.

  • For a free particle, d⟨p⟩/dt=0d\langle p\rangle/dt=0 and the center moves uniformly.
  • For a uniform force, V′(x)V'(x) is constant, so the mean obeys the exact classical acceleration law.
  • For a harmonic oscillator, V′(x)=mω2xV'(x)=m\omega^2x, so the first moments close exactly.
  • For generic anharmonic potentials, first moments couple to variances and higher central moments.

The width already has its own equation. If

Cxp=12⟨δx δp+δp δx⟩,C_{xp}=\frac12\langle\delta x\,\delta p+\delta p\,\delta x\rangle,

then

dσx2dt=2mCxp.\frac{d\sigma_x^2}{dt}=\frac{2}{m}C_{xp}.

The covariance in turn couples to momentum variance and force-position correlations. This hierarchy explains why correct center motion does not imply a rigid, classical packet.

For A=HA=H, the same general identity gives

ddt⟨H⟩=⟨∂H∂t⟩.\frac{d}{dt}\langle H\rangle =\left\langle\frac{\partial H}{\partial t}\right\rangle.

Energy is conserved for a time-independent Hamiltonian, provided the domain assumptions needed for the commutator argument hold.

Ehrenfest theorem does not by itself explain the full classical limit. It says something about expectation values, not about definite trajectories, measurement records, decoherence, or the suppression of interference between macroscopically different histories. For the broader conceptual bridge, see Correspondence Principle.

Even when ⟨x⟩\langle x\rangle follows a classical-looking trajectory, the wave packet can spread. A free Gaussian packet has uniform center motion while its width changes. The semiclassical discussion asks how long a localized packet remains a useful trajectory approximation.

  • Replacing ⟨V′(x)⟩\langle V'(x)\rangle by V′(⟨x⟩)V'(\langle x\rangle) without justification.
  • Saying Ehrenfest theorem proves that particles follow classical paths.
  • Forgetting that broad or split wave packets can have expectation values that are poor summaries of the state.
  • Ignoring wave-packet spreading.
  • Treating the theorem as a substitute for semiclassical analysis or decoherence.
  • P. Ehrenfest, “Bemerkung über die angenäherte Gültigkeit der klassischen Mechanik innerhalb der Quantenmechanik,” Zeitschrift für Physik 45, 455-457 (1927).
  • D. J. Griffiths and D. F. Schroeter, Introduction to Quantum Mechanics, 3rd ed., Cambridge University Press, 2018.
  • J. J. Sakurai and J. Napolitano, Modern Quantum Mechanics, 3rd ed., Cambridge University Press, 2020.
  • R. Shankar, Principles of Quantum Mechanics, 2nd ed., Springer, 1994.
  1. For the harmonic oscillator V(x)=mω2x2/2V(x)=m\omega^2x^2/2, with m>0m>0 and ω>0\omega>0, show that ⟨x⟩\langle x\rangle obeys the classical oscillator equation. Solve for the mean position given ⟨x⟩0=x0\langle x\rangle_0=x_0 and ⟨p⟩0=p0\langle p\rangle_0=p_0, assuming the moments and Ehrenfest derivatives exist. Does this determine the full quantum state?
Solution

Here V′(x)=mω2xV'(x)=m\omega^2x, so

⟨V′(x)⟩=mω2⟨x⟩.\langle V'(x)\rangle =m\omega^2\langle x\rangle.

Ehrenfest theorem gives

md2dt2⟨x⟩=−mω2⟨x⟩,m\frac{d^2}{dt^2}\langle x\rangle =-m\omega^2\langle x\rangle,

or

d2dt2⟨x⟩+ω2⟨x⟩=0.\frac{d^2}{dt^2}\langle x\rangle+\omega^2\langle x\rangle=0.

Since d⟨x⟩/dt=⟨p⟩/md\langle x\rangle/dt=\langle p\rangle/m, the initial data give

⟨x⟩(t)=x0cos⁡(ωt)+p0mωsin⁡(ωt).\langle x\rangle(t) =x_0\cos(\omega t)+\frac{p_0}{m\omega}\sin(\omega t).

This is exact for the mean under the stated assumptions, without a narrow-packet approximation. It does not determine the full state: different states can have the same x0,p0x_0,p_0 and different widths, higher moments, and interference structure.

  1. Show that expectation-value motion closes exactly for every state if V(x)V(x) is at most quadratic.
Solution

For V(x)=ax2+bx+cV(x)=ax^2+bx+c, the force is affine:

V′(x)=2ax+b.V'(x)=2ax+b.

Linearity of expectation values gives

⟨V′(x)⟩=2a⟨x⟩+b=V′(⟨x⟩).\langle V'(x)\rangle=2a\langle x\rangle+b =V'(\langle x\rangle).

Thus the equations for ⟨x⟩\langle x\rangle and ⟨p⟩\langle p\rangle form a closed classical system for every state, even though higher moments may still evolve.

  1. For V(x)=λx4/4V(x)=\lambda x^4/4 and a packet symmetric about its mean, compute the leading finite-width correction to the force on the mean.
Solution

Symmetry about the mean gives ⟨(δx)3⟩=0\langle(\delta x)^3\rangle=0. Since

⟨x3⟩=q3+3qσx2,\langle x^3\rangle=q^3+3q\sigma_x^2,

Ehrenfest’s equation becomes

mq¨=−λ(q3+3qσx2).m\ddot q=-\lambda(q^3+3q\sigma_x^2).

The correction to the point-particle force is −3λqσx2-3\lambda q\sigma_x^2.