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.

About these ads

Trackbacks

By The free particle: probability current « Physics tutorials on Tuesday, 24 July 2012 at 16:06

[...] of change of probability of a particle in a given range of can be written as the difference in probability current at the two ends. The current is defined [...]
By Probability current in 3-d « Physics tutorials on Friday, 18 January 2013 at 16:57

[...] seen the concept of the probability current in one dimension. A generalization to 3-d gives the [...]

Physics tutorials