Ehrenfest theorem

The Ehrenfest theorem, named after Paul Ehrenfest, an Austrian theoretical physicist at Leiden University, relates the time derivative of the expectation values of the position and momentum operators x and p to the expectation value of the force F = −dV/dx on a massive particle moving in a scalar potential,

$m\frac{d}{dt}\langle x\rangle = \langle p\rangle,\;\; \frac{d}{dt}\langle p\rangle = -\left\langle \frac{\partial V(x)}{\partial x}\right\rangle ~.$

Loosely speaking, one can thus say that "quantum mechanical expectation values obey Newton’s classical equations of motion". (This loose statement needs some caveats, see.^[1])

The Ehrenfest theorem is a special case of a more general relation between the expectation of any quantum mechanical operator and the expectation of the commutator of that operator with the Hamiltonian of the system ^[2][3]

$\frac{d}{dt}\langle A\rangle = \frac{1}{i\hbar}\langle [A,H] \rangle+ \left\langle \frac{\partial A}{\partial t}\right\rangle ~,$

where A is some QM operator and ⟨A⟩ is its expectation value. This more general theorem was not actually derived by Ehrenfest (it is due to Werner Heisenberg).

It is most apparent in the Heisenberg picture of quantum mechanics, where it is just the expectation value of the Heisenberg equation of motion. It provides mathematical support to the correspondence principle.

The reason is that Ehrenfest's theorem is closely related to Liouville's theorem of Hamiltonian mechanics, which involves the Poisson bracket instead of a commutator. Dirac's rule of thumb suggests that statements in quantum mechanics which contain a commutator correspond to statements in classical mechanics where the commutator is supplanted by a Poisson bracket multiplied by iħ. This makes the operator expectation values obey corresponding classical equations of motion, provided the Hamiltonian is at most quadratic in the coordinates and momenta. Otherwise, the evolution equations still may hold approximately, provided fluctuations are small.

Derivation in the Schrödinger picture

Suppose some system is presently in a quantum state Φ. If we want to know the instantaneous time derivative of the expectation value of A, that is, by definition

$\begin{align} \frac{\mathrm{d}}{\mathrm{d}t}\langle A\rangle &= \frac{\mathrm{d}}{\mathrm{d}t}\int \Phi^* A \Phi\,\mathrm{d}^3x \\ &= \int \left( \frac{\partial \Phi^*}{\partial t} \right) A\Phi\,\mathrm{d}^3x + \int \Phi^* \left( \frac{\partial A}{\partial t}\right) \Phi\,\mathrm{d}^3x +\int \Phi^* A \left( \frac{\partial \Phi}{\partial t} \right)\,\mathrm{d}^3x \\ &= \int \left( \frac{\partial \Phi^*}{\partial t} \right) A\Phi\,\mathrm{d}^3x + \left\langle \frac{\partial A}{\partial t}\right\rangle + \int \Phi^* A \left( \frac{\partial \Phi}{\partial t} \right)\,\mathrm{d}^3x \end{align}$

where we are integrating over all space. If we apply the Schrödinger equation, we find that

$\frac{\partial \Phi}{\partial t} = \frac{1}{i\hbar}H\Phi$

By taking the complex conjugate we find

$\frac{\partial \Phi^*}{\partial t} = -\frac{1}{i\hbar}\Phi^*H^* = -\frac{1}{i\hbar}\Phi^*H.$ ^[4]

Note H = H ^∗, because the Hamiltonian is Hermitian. Placing this into the above equation we have

$\frac{d}{dt}\langle A\rangle = \frac{1}{i\hbar}\int \Phi^* (AH-HA) \Phi~dx^3 + \left\langle \frac{\partial A}{\partial t}\right\rangle = \frac{1}{i\hbar}\langle [A,H]\rangle + \left\langle \frac{\partial A}{\partial t}\right\rangle.$

Often (but not always) the operator A is time independent, so that its derivative is zero and we can ignore the last term.

Derivation in the Heisenberg Picture

In the Heisenberg picture, the derivation is trivial. The Heisenberg picture moves the time dependence of the system to operators instead of state vector. Starting with the Heisenberg equation of motion

$\frac{d}{dt}A(t) = \frac{\partial A(t)}{\partial t} + \frac{1}{i \hbar}[A(t),H],$

we can derive Ehrenfest's theorem simply by projecting the Heisenberg equation onto $|\Psi\rangle$ from the right and $\langle\Psi|$ from the left, or taking the expectation value, so

$\langle\Psi|\frac{d}{dt}A(t)|\Psi\rangle = \langle\Psi|\frac{\partial A(t)}{\partial t}|\Psi\rangle + \langle\Psi|\frac{1}{i \hbar}[A(t),H]|\Psi\rangle,$

We can pull the d/dt out of the first term since the state vectors are no longer time dependent in the Heisenberg Picture. Therefore,

$\frac{d}{dt}\langle A(t)\rangle = \left\langle\frac{\partial A(t)}{\partial t}\right\rangle + \frac{1}{i \hbar}\langle[A(t),H]\rangle$

General example

The expectation values of the theorem, however, are the very same in the Schrödinger picture as well. For the very general example of a massive particle moving in a potential, the Hamiltonian is simply

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

where x is the position of the particle.

Suppose we wanted to know the instantaneous change in momentum p. Using Ehrenfest's theorem, we have

$\frac{d}{dt}\langle p\rangle = \frac{1}{i\hbar}\langle [p,H]\rangle + \left\langle \frac{\partial p}{\partial t}\right\rangle = \frac{1}{i\hbar}\langle [p,V(x,t)]\rangle,$

since the operator p commutes with itself and has no time dependence.^[5] By expanding the right-hand-side, replacing p by −iħ∇, we get

$\frac{d}{dt}\langle p\rangle = \int \Phi^* V(x,t)\frac{d}{dx}\Phi~dx - \int \Phi^* \frac{d}{dx} (V(x,t)\Phi)~dx ~.$

After applying the product rule on the second term, we have

$\begin{align} \frac{d}{dt}\langle p\rangle &= \int \Phi^* V(x,t) \frac{d}{dx}\Phi~dx - \int \Phi^* (\frac{d}{dx} V(x,t))\Phi ~dx - \int \Phi^* V(x,t) \frac{d}{dx}\Phi~dx \\ &= - \int \Phi^* (\frac{d}{dx} V(x,t))\Phi ~dx \\ &= \langle - \frac{d}{dx} V(x,t)\rangle = \langle F \rangle, \end{align}$

but we recognize this as Newton's second law. This is an example of the correspondence principle: the result manifests as Newton's second law in the case of having so many excitations superposed in the wavefunction that the net motion is given by the expectation value simulating a classical particle.

Similarly we can obtain the instantaneous change in the position expectation value.

$\begin{align} \frac{d}{dt}\langle x\rangle &= \frac{1}{i\hbar}\langle [x,H]\rangle + \left\langle \frac{\partial x}{\partial t}\right\rangle \\ &= \frac{1}{i\hbar} \left \langle \left [x,\frac{p^2}{2m} + V(x,t) \right ] \right \rangle + 0 \\ &= \frac{1}{i\hbar} \left \langle \left [x,\frac{p^2}{2m} \right] \right \rangle \\ &= \frac{1}{i\hbar 2 m} \left \langle [x,p] \frac{d}{dp} p^2 \right\rangle \\ &= \frac{1}{i\hbar 2 m}\langle i \hbar 2 p\rangle \\ &= \frac{1}{m}\langle p\rangle \end{align}$

This result is again in accord with the classical equation.

Derivation of the Schrödinger equation from the Ehrenfest theorems

It was established above that the Ehrenfest theorems are consequences of the Schrödinger equation. However, the converse is also true: the Schrödinger equation can be inferred from the Ehrenfest theorems.^[6] We begin from

$\begin{align} m\frac{d}{dt} \left \langle \Psi(t) \right | \hat{x} \left | \Psi(t) \right \rangle &= \left \langle \Psi(t) \right | \hat{p} \left | \Psi(t) \right \rangle, \\ \frac{d}{dt} \left \langle \Psi(t) \right | \hat{p} \left | \Psi(t) \right \rangle &= \left \langle \Psi(t) \right | -V'(\hat{x}) \left | \Psi(t) \right \rangle. \end{align}$

Applications of the product rule leads to

$\begin{align} \left \langle \frac{d\Psi}{dt} \Big | \hat{x} \Big | \Psi \right \rangle + \left \langle \Psi \Big | \hat{x} \Big | \frac{d\Psi}{dt} \right \rangle &= \left \langle \Psi \Big | \frac{\hat{p}}{m} \Big | \Psi \right \rangle, \\ \left \langle \frac{d\Psi}{dt} \Big | \hat{p} \Big | \Psi \right \rangle + \left \langle \Psi \Big | \hat{p} \Big | \frac{d\Psi}{dt} \right \rangle &= \langle \Psi | -V'(\hat{x}) | \Psi \rangle, \end{align}$

into which we substitute a consequence of Stone's theorem

$i\hbar \left | \frac{d\Psi}{dt} \right \rangle = \hat{H} | \Psi(t) \rangle ~,$

where ħ was introduced as a normalization constant to the balance dimensionality. Since these identities must be valid for any initial state, the averaging can be dropped and the system of commutator equations for the unknown quantum generator of motion Ĥ are derived

$im [\hat{H}, \hat{x}] = \hbar \hat{p}, \qquad i [\hat{H}, \hat{p}] = -\hbar V'(\hat{x}).$

Assuming that observables of the coordinate and momentum obey the canonical commutation relation [x̂, p̂] = iħ. Setting $\hat{H} = H(\hat{x}, \hat{p})$, the commutator equations can be converted into the differential equations^[6][7]

$m \frac{\partial H (x,p)}{\partial p} = p, \qquad \frac{\partial H(x,p)}{\partial x} = V'(x),$

whose solution is the familiar quantum Hamiltonian

$\hat{H} = \frac{\hat{p}^2}{2m} + V(\hat{x}).$

Whence, the Schrödinger equation was derived from the Ehrenfest theorems by assuming the canonical commutation relation between the coordinate and momentum. If one assumes that the coordinate and momentum commute, the same computational method leads to the Koopman–von Neumann classical mechanics, which is the Hilbert space formulation of classical mechanics.^[6] Therefore, this derivation as well as the derivation of the Koopman–von Neumann mechanics shows that the essential difference between quantum and classical mechanics reduces to the value of the commutator [x̂, p̂].

Notes

1. Wheeler, Nicholas. "Remarks concerning the status & some ramifications of Ehrenfest's theorem" (PDF).
2. Ehrenfest, P. (1927). "Bemerkung über die angenäherte Gültigkeit der klassischen Mechanik innerhalb der Quantenmechanik". Zeitschrift für Physik 45 (7–8): 455–457. Bibcode:1927ZPhy...45..455E. doi:10.1007/BF01329203.
3. Smith, Henrik (1991). Introduction to Quantum Mechanics. World Scientific Pub Co Inc. pp. 108–109. ISBN 978-9810204754.
4. In bra–ket notation, $\frac{\partial}{\partial t}\langle \phi |x\rangle =\frac{-1}{i\hbar}\langle \phi |\hat{H}|x\rangle =\frac{-1}{i\hbar}\langle \phi |x \rangle H=\frac{-1}{i\hbar}\Phi^*H,$ where $\hat{H}$ is the Hamiltonian operator, and H is the Hamiltonian represented in coordinate space (as is the case in the derivation above). In other words, we applied the adjoint operation to the entire Schrödinger equation, which flipped the order of operations for H and Φ.
5. Although the expectation value of the momentum ⟨p⟩, which is a real-number-valued function of time, will have time dependence, the momentum operator itself, p does not, in this picture: Rather, the momentum operator is a constant linear operator on the Hilbert space of the system. The time dependence of the expectation value, in this picture, is due to the time evolution of the wavefunction for which the expectation value is calculated. An Ad hoc example of an operator which does have time dependence is ⟨xt²⟩, where x is the ordinary position operator and t is just the (non-operator) time, a parameter.
6. Bondar, D.; Cabrera, R.; Lompay, R.; Ivanov, M.; Rabitz, H. (2012). "Operational Dynamic Modeling Transcending Quantum and Classical Mechanics". Physical Review Letters 109 (19). arXiv:1105.4014. Bibcode:2012PhRvL.109s0403B. doi:10.1103/PhysRevLett.109.190403.
7. Transtrum, M. K.; Van Huele, J. F. O. S. (2005). "Commutation relations for functions of operators". Journal of Mathematical Physics 46 (6): 063510. Bibcode:2005JMP....46f3510T. doi:10.1063/1.1924703.