← Exercise library← All exercises for this theme
Theme 3 — Postulates and Applications

Probability current

Keywords: probability current · continuity equation · plane wave · standing wave · evanescent wave · velocity

Exercise 1 : Probability current

For a particle of mass mm described by a one-dimensional wave function ψ(x)\psi(x), the probability current is defined by

j=ℏm Im⁡(ψ∗ dψdx).j=\frac{\hbar}{m}\,\operatorname{Im}\Bigl(\psi^*\,\frac{d\psi}{dx}\Bigr).

For a stationary state, the continuity equation ∂tρ+∂xj=0\partial_t\rho+\partial_xj=0 requires jj to be independent of xx. Here we calculate the current for several standard wave functions. Throughout the exercise, kk and κ\kappa are strictly positive real numbers, and AA, BB are complex constants.

Question 1
Calculate the current of the plane wave ψ(x)=Aeikx\psi(x)=Ae^{ikx}. Interpret the result using the classical velocity v=p/mv=p/m, where p=ℏkp=\hbar k.

Solution
We have ψ′=ikAeikx\psi'=ikAe^{ikx}, so

ψ∗ψ′=A∗e−ikx⋅ikAeikx=ik∣A∣2,j=ℏmIm⁡(ik∣A∣2)=ℏkm∣A∣2.\psi^*\psi'=A^*e^{-ikx}\cdot ikAe^{ikx}=ik|A|^2, \qquad j=\frac{\hbar}{m}\operatorname{Im}(ik|A|^2)=\frac{\hbar k}{m}|A|^2 .

With p=ℏkp=\hbar k, we recognise j=v ρj=v\,\rho, where ρ=∣A∣2\rho=|A|^2 is the probability density and v=p/mv=p/m the classical velocity. This is the usual relation between flux, density, and velocity in a fluid: the plane wave describes a uniform probability flux moving towards increasing xx at velocity vv.

Question 2
Calculate the current of the superposition ψ(x)=Aeikx+Be−ikx\psi(x)=Ae^{ikx}+Be^{-ikx}. Show that the cross terms do not contribute. Hence find the current of a standing wave ψ(x)=Csin⁡(kx)\psi(x)=C\sin(kx).

Solution
We have ψ′=ik(Aeikx−Be−ikx)\psi'=ik(Ae^{ikx}-Be^{-ikx}), and

ψ∗ψ′=ik(A∗e−ikx+B∗eikx)(Aeikx−Be−ikx)=ik(∣A∣2−∣B∣2+B∗Ae2ikx−A∗Be−2ikx).\psi^*\psi'=ik\bigl(A^*e^{-ikx}+B^*e^{ikx}\bigr)\bigl(Ae^{ikx}-Be^{-ikx}\bigr) =ik\Bigl(|A|^2-|B|^2+B^*Ae^{2ikx}-A^*Be^{-2ikx}\Bigr).

Let z=B∗Ae2ikxz=B^*Ae^{2ikx}. The two cross terms are z−z∗=2iIm⁡zz-z^*=2i\operatorname{Im}z, which is purely imaginary. Multiplication by ikik gives −2kIm⁡z-2k\operatorname{Im}z, which is real and therefore does not contribute to the imaginary part. Thus

j=ℏkm(∣A∣2−∣B∣2).j=\frac{\hbar k}{m}\bigl(|A|^2-|B|^2\bigr).

The total current is the difference between the current of the wave propagating to the right and that of the wave propagating to the left: the two waves transport probability independently of each other. For the standing wave Csin⁡(kx)=C2i(eikx−e−ikx)C\sin(kx)=\frac{C}{2i}(e^{ikx}-e^{-ikx}), we have ∣A∣=∣B∣=∣C∣/2|A|=|B|=|C|/2, and the current vanishes: the two opposing fluxes cancel. This can also be seen directly: ψ∗ψ′=∣C∣2ksin⁡(kx)cos⁡(kx)\psi^*\psi'=|C|^2k\sin(kx)\cos(kx) is real.

Question 3
More generally, show that the current vanishes for any wave function of the form ψ(x)=eiθf(x)\psi(x)=e^{i\theta}f(x), where ff is a real function and θ\theta a constant. Apply this to the evanescent wave ψ(x)=Ae−κx\psi(x)=Ae^{-\kappa x}.

Solution
We have ψ∗ψ′=e−iθf⋅eiθf′=ff′\psi^*\psi'=e^{-i\theta}f\cdot e^{i\theta}f'=ff', which is real: its imaginary part vanishes, and j=0j=0. The constant phase θ\theta plays no role. The evanescent wave Ae−κxAe^{-\kappa x} has this form, with f(x)=∣A∣e−κxf(x)=|A|e^{-\kappa x} and θ=arg⁡A\theta=\arg A. It therefore carries no current: this is why an evanescent wave beyond an infinite potential step corresponds to no transmission.

Question 4
Calculate the current of ψ(x)=Aeκx+Be−κx\psi(x)=Ae^{\kappa x}+Be^{-\kappa x}. Under what condition on AA and BB is it non-zero? Why is this result useful for understanding tunnelling?

Solution
We have ψ′=κ(Aeκx−Be−κx)\psi'=\kappa(Ae^{\kappa x}-Be^{-\kappa x}), and

ψ∗ψ′=κ(A∗eκx+B∗e−κx)(Aeκx−Be−κx)=κ(∣A∣2e2κx−∣B∣2e−2κx+B∗A−A∗B).\psi^*\psi'=\kappa\bigl(A^*e^{\kappa x}+B^*e^{-\kappa x}\bigr)\bigl(Ae^{\kappa x}-Be^{-\kappa x}\bigr) =\kappa\Bigl(|A|^2e^{2\kappa x}-|B|^2e^{-2\kappa x}+B^*A-A^*B\Bigr).

The first two terms are real. The last two can be written B∗A−(B∗A)∗=2iIm⁡(AB∗)B^*A-(B^*A)^*=2i\operatorname{Im}(AB^*). Hence

j=2ℏκmIm⁡(AB∗).j=\frac{2\hbar\kappa}{m}\operatorname{Im}(AB^*).

As required, this current is independent of xx. It is non-zero if and only if AA and BB are both non-zero and do not have the same phase modulo π\pi. A single real exponential carries no current, but a combination of the two exponentials with coefficients of different phases does. This is exactly what occurs inside a potential barrier of finite width: the matching conditions impose such a combination, and the current through the barrier equals the transmitted current.

Question 5
Write an arbitrary wave function in the form ψ(x)=f(x) eiS(x)/ℏ\psi(x)=f(x)\,e^{iS(x)/\hbar}, with ff and SS real and f≥0f\geq0. Show that j=ρ vj=\rho\,v, with ρ=f2\rho=f^2 and v=S′(x)/mv=S'(x)/m. Recover the result of the first question.

Solution
Differentiating the product gives

ψ′=(f′+iℏfS′)eiS/ℏ,ψ∗ψ′=ff′+iℏf2S′.\psi'=\Bigl(f'+\frac{i}{\hbar}fS'\Bigr)e^{iS/\hbar}, \qquad \psi^*\psi'=ff'+\frac{i}{\hbar}f^2S' .

The imaginary part is f2S′/ℏf^2S'/\hbar, and therefore

j=ℏm⋅f2S′ℏ=f2 S′m=ρ v.j=\frac{\hbar}{m}\cdot\frac{f^2S'}{\hbar}=f^2\,\frac{S'}{m}=\rho\,v .

The current depends only on the phase gradient: it is the phase of the wave function that encodes the local velocity of the probability flow. For the plane wave, f=∣A∣f=|A| and S=ℏkx+ℏarg⁡AS=\hbar kx+\hbar\arg A, so S′=ℏkS'=\hbar k and v=ℏk/mv=\hbar k/m, in agreement with the first question.