### Ehrenfest's Theorem and Derivation

Ehrenfest's Theorem Statement:

The theorem states that

Proof of theorem:

The proof of the theorem for one-dimensional motion of a particle by showing that

1) $\frac{d \left < x \right >}{dt} = \frac{\left < p_{x} \right > }{m}$

2) $\frac{d \left < p_{x} \right >}{dt} = \left < F_{x} \right >$

1.) To Show that: $\frac{d \left < x \right > }{dt} = \frac{\left < p_{x} \right > }{m}$

Let $x$ is the position coordinate of a particle of mass $m$, at time $t$

The expectation value of position $x$ of a particle is given by

$\left < x \right > = \int_{- \infty}^{+ \infty} \psi^{*} (x,t) . x \: \psi (x,t) dx \qquad (1)$

Differentiating the above equation $(1)$ with respect to $t$

$\frac{d \left < x \right > }{dt} = \int_{- \infty}^{+ \infty} x \frac{\partial (\psi \psi^{*})}{\partial t} dx \qquad(2)$

We know the probablity current density

$\frac{\partial (\psi \psi^{*})}{\partial t} = \frac{i \hbar}{2m} \frac{\partial}{\partial x} \left[ \psi^{*} \frac{\partial \psi}{\partial x} - \psi \frac{\partial \psi^{*} }{\partial x} \right] \qquad(3)$

Now substitute the above eqaution$(3)$ in eqaution $(2)$

$\frac{d \left < x \right > }{dt} = \frac{i \hbar}{2m} \int_{- \infty}^{+ \infty} x \frac{\partial}{\partial x} \left[ \psi^{*} \frac{\partial \psi}{\partial x} - \psi \frac{\partial \psi^{*} }{\partial x} \right] dx $

Integrating the right-hand side by parts of the above equation, we get

$\frac{d \left < x \right > }{dt} = \frac{i \hbar}{2m} \left[ x \left( \psi^{*} \frac{\partial \psi}{\partial x} - \psi \frac{\partial \psi^{*} }{\partial x} \right) \right]^{+\infty}_{-\infty} \\
\qquad
- \frac{i \hbar}{2m} \int_{- \infty}^{+ \infty} \left( \psi^{*} \frac{\partial \psi}{\partial x} - \psi \frac{\partial \psi^{*} }{\partial x} \right) dx $

As $x$ approaches either $+ \infty$ or $-\infty$, $\psi$ and $\frac{\partial \psi}{\partial x}$ approach zero, and therefore the first term becomes zero.

Hence we get

$\frac{d \left < x \right > }{dt} = - \frac{i \hbar}{2m} \int_{- \infty}^{+ \infty} \left( \psi^{*} \frac{\partial \psi}{\partial x} - \psi \frac{\partial \psi^{*} }{\partial x} \right) dx
\\
\qquad\qquad\qquad\qquad\qquad\qquad ---(4)$

The expectation value of $p_{x}$ is given by

$ \left < p_{x} \right > = \int_{- \infty}^{+ \infty} \psi^{*} \frac{\hbar}{i} \frac{\partial \psi}{\partial x} $

$ \int_{- \infty}^{+ \infty} \psi^{*} \frac{\partial \psi}{\partial x} dx = \frac{i}{\hbar}\left < p_{x} \right > \qquad(5)$

Similarly

$ \int_{- \infty}^{+ \infty} \psi \frac{\partial \psi^{*}}{\partial x} dx = - \frac{i}{\hbar}\left < p_{x} \right > \qquad(6)$

Substituting the values of these integrals in equation $(4)$

$\frac{d \left < x \right > }{dt} = - \frac{i \hbar}{2m} \left[ \frac{i}{\hbar}\left < p_{x} \right > + \frac{i}{\hbar}\left < p_{x} \right >\right] $

$\frac{d \left < x \right > }{dt} = - \frac{\left < p_{x} \right >}{m} \qquad(7)$

This is the first result of Ehrenfest's Theorem.

2) To show that: $\frac{d \left < p_{x} \right >}{dt} = \left < F_{x} \right >$

We know that the expectation value of the momentum $p_{x}$ is given by

$ \left < p_{x} \right > = \int_{- \infty}^{+ \infty} \psi^{*} \frac{\hbar}{i} \frac{\partial \psi}{\partial x} $

$ \left < p_{x} \right > =\frac{\hbar}{i} \int_{- \infty}^{+ \infty} \psi^{*} \frac{\partial \psi}{\partial x} \qquad(8)$

Differentiating the equation $(8)$ with respect to $t$, we get

$\frac{d \left < p_{x} \right >}{dt} = \frac{\hbar}{i} \int_{- \infty}^{+ \infty} \left[ \frac{\partial \psi^{*}}{\partial t} \frac{\partial \psi}{\partial x} + \psi^{*} \frac{\partial^{2} \psi}{\partial x \partial t} \right] $

$\frac{d \left < p_{x} \right >}{dt} = \int_{- \infty}^{+ \infty} \left[-i \hbar \frac{\partial \psi^{*}}{\partial t} \frac{\partial \psi}{\partial x} - i\hbar \psi^{*} \frac{\partial^{2} \psi}{\partial x \partial t} \right]
\\
\qquad\qquad\qquad\qquad\qquad\qquad ---(9)$

Now the time-dependent Schrodinger equations for $\psi$ and $\psi^{*}$ are

$i \hbar \frac{\partial \psi}{\partial t} =- \frac{\hbar^{2}}{2m} \frac{\partial^{2} \psi}{\partial x^{2}} + V \psi \qquad(10)$

The complex conjugate of Schrodinger function

$-i \hbar \frac{\partial \psi^{*}}{\partial t} = -\frac{\hbar^{2}}{2m} \frac{\partial^{2} \psi^{*}}{\partial x^{2}} + V \psi^{*} \qquad(11)$

Differentiating the equation $(10)$ with respect to $x$

$i \hbar \frac{\partial^{2} \psi}{\partial x \partial t} = - \frac{\hbar^{2}}{2m} \frac{\partial^{3} \psi}{\partial x^{3}} + \frac{\partial (V \psi)}{\partial x} \qquad(12)$

Now substitute the value of $-i \hbar \frac{\partial \psi^{*}}{\partial t}$ and $i \hbar \frac{\partial^{2} \psi}{\partial x \partial t}$ in the equation $(9)$, we get

Now put $\frac{\partial (V \psi)}{\partial x}= \left\{ \psi \frac{\partial V }{\partial x}+ V\frac{\partial \psi}{\partial x} \right\}$ in above equation:

$\frac{d \left < p_{x} \right >}{dt} = -\frac{\hbar^{2}}{2m} \int_{- \infty}^{+ \infty} \frac{\partial}{\partial x} \left( \frac{\partial \psi^{*}}{\partial x}\frac{\partial \psi}{\partial x} - \psi^{*} \frac{\partial^{2} \psi}{\partial x^{2}} \right) dx + \int_{- \infty}^{+ \infty} \left( V \psi^{*} \frac{\partial \psi}{\partial x} - \psi^{*} \left\{ \psi \frac{\partial V }{\partial x}+ V\frac{\partial \psi}{\partial x} \right\} \right) dx$

As $x$ approaches either $+ \infty $ or $-\infty$ and $\frac{\partial \psi}{\partial x}$ is zero. Therefore the first term of the above equation on the right-hand side will be zero.

$\frac{d \left < p_{x} \right >}{dt} = \int_{- \infty}^{+ \infty} - \psi \frac{\partial V }{\partial x} \psi^{*} dx$

$\frac{d \left < p_{x} \right >}{dt} = - \int_{- \infty}^{+ \infty} \psi \frac{\partial V }{\partial x} \psi^{*} dx$

$\frac{d \left < p_{x} \right >}{dt} = -\left < \frac{\partial V }{\partial x} \right > $

Here the $\left < \frac{\partial V }{\partial x} \right >$ is the average value or expectation value of potential gradient and the negative value of the potential gradient is equal to the average value or expectation value of force $\left < F_{x} \right >$ along the $x$ direction.

$\frac{d \left < p_{x} \right >}{dt} = \left < F_{x} \right > $

This is the second result of Ehrenfest theorem and it represents Newton's second law of motion. Thus if the expectation values of dynamical quantities for a particle are, considered, quantum mechanics given the equations of classical mechanics.

Quantum mechanics gives the same results as classical mechanics for a particle for which average or expectation values of dynamical quantities are involved.

$\frac{d \left < p_{x} \right >}{dt} = \int_{- \infty}^{+ \infty} \left[ \left( -\frac{\hbar^{2}}{2m} \frac{\partial^{2} \psi^{*}}{\partial x^{2}} + V \psi^{*} \right) \frac{\partial \psi}{\partial x} - \psi^{*} \left( - \frac{\hbar^{2}}{2m} \frac{\partial^{3} \psi}{\partial x^{3}} + \frac{\partial (V \psi)}{\partial x} \right) \right]$

$\frac{d \left < p_{x} \right >}{dt} = \int_{- \infty}^{+ \infty} \left[-\frac{\hbar^{2}}{2m} \left( \frac{\partial^{2} \psi^{*}}{\partial x^{2}}\frac{\partial \psi}{\partial x} - \psi^{*} \frac{\partial^{3} \psi}{\partial x^{3}} \right) - \left( V \psi^{*} \frac{\partial \psi}{\partial x} - \psi^{*} \frac{\partial (V \psi)}{\partial x} \right) \right] dx$

$\frac{d \left < p_{x} \right >}{dt} = -\frac{\hbar^{2}}{2m} \int_{- \infty}^{+ \infty} \left( \frac{\partial^{2} \psi^{*}}{\partial x^{2}}\frac{\partial \psi}{\partial x} - \psi^{*} \frac{\partial^{3} \psi}{\partial x^{3}} \right) dx + \int_{- \infty}^{+ \infty} \left( V \psi^{*} \frac{\partial \psi}{\partial x} - \psi^{*} \frac{\partial (V \psi)}{\partial x} \right) dx$

$\frac{d \left < p_{x} \right >}{dt} = -\frac{\hbar^{2}}{2m} \int_{- \infty}^{+ \infty} \frac{\partial}{\partial x} \left( \frac{\partial \psi^{*}}{\partial x}\frac{\partial \psi}{\partial x} - \psi^{*} \frac{\partial^{2} \psi}{\partial x^{2}} \right) dx + \int_{- \infty}^{+ \infty} \left( V \psi^{*} \frac{\partial \psi}{\partial x} - \psi^{*} \frac{\partial (V \psi)}{\partial x} \right) dx$

$\frac{d \left < p_{x} \right >}{dt} = -\frac{\hbar^{2}}{2m} \left[ \frac{\partial \psi^{*}}{\partial x}\frac{\partial \psi}{\partial x} - \psi^{*} \frac{\partial^{2} \psi}{\partial x^{2}} \right]_{- \infty}^{+ \infty} + \int_{- \infty}^{+ \infty} \left( V \psi^{*} \frac{\partial \psi}{\partial x} - \psi^{*} \left\{ \psi \frac{\partial V }{\partial x}+ V\frac{\partial \psi}{\partial x} \right\} \right) dx$

$\frac{d \left < p_{x} \right >}{dt} = \int_{- \infty}^{+ \infty} \left( V \psi^{*} \frac{\partial \psi}{\partial x} - \psi^{*} \left\{ \psi \frac{\partial V }{\partial x}+ V\frac{\partial \psi}{\partial x} \right\} \right) dx$

$\frac{d \left < p_{x} \right >}{dt} = \int_{- \infty}^{+ \infty} \left( V \psi^{*} \frac{\partial \psi}{\partial x} - \psi \frac{\partial V }{\partial x} \psi^{*} - V \psi^{*} \frac{\partial \psi}{\partial x} \right) dx$