## Probability current

Required math: calculus

Required physics: Schrödinger equation

Reference: Griffiths, David J. (2005), Introduction to Quantum Mechanics, 2nd Edition; Pearson Education – Problem 1.14.

Using the Schrödinger equation we can derive an interesting quantity called the probability current. Using the probabilistic interpretation of the wave function, the probability of a particle being between ${x=a}$ and ${x=b}$ is

$\displaystyle P_{ab}=\int_{a}^{b}\left|\Psi^{*}\Psi\right|dx$

The rate of change of this probability can then be expressed in terms of spatial derivatives using the Schrödinger equation:

 $\displaystyle \frac{dP_{ab}}{dt}$ $\displaystyle =$ $\displaystyle \int_{a}^{b}\left[\frac{\partial\Psi^{*}}{\partial t}\Psi+\frac{\partial\Psi}{\partial t}\Psi^{*}\right]dx$ $\displaystyle$ $\displaystyle =$ $\displaystyle \int_{a}^{b}\left\{ -\frac{i\hbar}{2m}\frac{\partial^{2}\Psi^{*}}{\partial x^{2}}-\frac{1}{i\hbar}V\Psi^{*}\right\} \Psi dx$ $\displaystyle$ $\displaystyle$ $\displaystyle +\int_{a}^{b}\left\{ \frac{i\hbar}{2m}\frac{\partial^{2}\Psi}{\partial x^{2}}+\frac{1}{i\hbar}V\Psi\right\} \Psi^{*}dx$ $\displaystyle$ $\displaystyle =$ $\displaystyle \frac{i\hbar}{2m}\int_{a}^{b}\left[\frac{\partial^{2}\Psi}{\partial x^{2}}\Psi^{*}-\frac{\partial^{2}\Psi^{*}}{\partial x^{2}}\Psi\right]dx$

We can now apply integration by parts to each term.

 $\displaystyle \int_{a}^{b}\frac{\partial^{2}\Psi}{\partial x^{2}}\Psi^{*}dx$ $\displaystyle =$ $\displaystyle \left.\frac{\partial\Psi}{\partial x}\Psi^{*}\right|_{a}^{b}-\int_{a}^{b}\frac{\partial\Psi}{\partial x}\frac{\partial\Psi^{*}}{\partial x}dx$ $\displaystyle -\int_{a}^{b}\frac{\partial^{2}\Psi^{*}}{\partial x^{2}}\Psi dx$ $\displaystyle =$ $\displaystyle \left.-\frac{\partial\Psi^{*}}{\partial x}\Psi\right|_{a}^{b}+\int_{a}^{b}\frac{\partial\Psi}{\partial x}\frac{\partial\Psi^{*}}{\partial x}dx$

Adding these terms together, we get

$\displaystyle \frac{dP_{ab}}{dt}=\frac{i\hbar}{2m}\left[\left.\frac{\partial\Psi}{\partial x}\Psi^{*}\right|_{a}^{b}\left.-\frac{\partial\Psi^{*}}{\partial x}\Psi\right|_{a}^{b}\right]$

If we define the probability current as

$\displaystyle J(x,t)\equiv\frac{i\hbar}{2m}\left(\frac{\partial\Psi^{*}}{\partial x}\Psi-\frac{\partial\Psi}{\partial x}\Psi^{*}\right)$

we can write the rate of change of probability as

$\displaystyle \frac{dP_{ab}}{dt}=J(a,t)-J(b,t)$

As an example, if the wave function is given by

$\displaystyle \Psi(x,t)=\left(\frac{2am}{\pi\hbar}\right)^{1/4}e^{-a\left[\left(mx^{2}/\hbar\right)+it\right]}$

(we’ve taken the constant in front so that it normalizes the wave function), then

$\displaystyle \frac{\partial\Psi}{\partial x}=-\frac{1}{\pi^{1/4}}\left(\frac{2am}{\hbar}\right)^{5/4}xe^{-a\left[\left(mx^{2}/\hbar\right)+it\right]}$

So for the probability current we get

 $\displaystyle J(x,t)$ $\displaystyle =$ $\displaystyle -\frac{i\hbar}{2m}\frac{1}{\pi^{1/4}}\left(\frac{2am}{\hbar}\right)^{5/4}xe^{-a\left[\left(mx^{2}/\hbar\right)-it\right]}\left(\frac{2am}{\pi\hbar}\right)^{1/4}e^{-a\left[\left(mx^{2}/\hbar\right)+it\right]}$ $\displaystyle$ $\displaystyle$ $\displaystyle +\frac{i\hbar}{2m}\frac{1}{\pi^{1/4}}\left(\frac{2am}{\hbar}\right)^{5/4}xe^{-a\left[\left(mx^{2}/\hbar\right)+it\right]}\left(\frac{2am}{\pi\hbar}\right)^{1/4}e^{-a\left[\left(mx^{2}/\hbar\right)-it\right]}$ $\displaystyle$ $\displaystyle =$ $\displaystyle 0$

A bit of an anti-climax after all that mathematics.