Solution of Schrödinger equation for a step potential

In quantum mechanics and scattering theory, the one dimensional step potential is an idealized system used to model incident, reflected and transmitted matter waves. The problem consists of solving the time-independent Schrödinger equation for a particle with a step-like potential in one dimension. Typically, the potential is modelled as a Heaviside step function.

Calculation

Schrödinger equation and potential function

Scattering at a finite potential step of height V₀, shown in green. The amplitudes and direction of left and right moving waves are indicated. Yellow is the incident wave, blue are reflected and transmitted, red does not occur. E > V₀ for this figure.

The time-independent Schrödinger equation for the wave function $\psi(x)$ is

$H\psi(x)=\left[-\frac{\hbar^2}{2m}\frac{d^2}{dx^2}+V(x)\right]\psi(x)=E\psi(x),$

where H is the Hamiltonian, ħ is the reduced Planck constant, m is the mass, E the energy of the particle. The step potential is simply the Heaviside step function:

$V(x) = \Theta(x)=\left\langle \begin{align} 0 & \quad x<0 \\ 1 & \quad x>0 \end{align}\right.$

where V₀ is the height of the barrier.

The barrier is positioned at x = 0, though any position x₀ may be chosen without changing the results, simply by shifting position of the step by −x₀.

The first term in the Hamiltonian, $-\frac{\hbar^2}{2m}\frac{d^2}{dx^2}\psi$ is the kinetic energy of the particle.

Solution

The step divides space in two parts: x < 0 and x > 0. In any of these parts the potential is constant, meaning the particle is quasi-free, and the solution of the Schrödinger equation can be written as a superposition of left and right moving waves (see free particle)

$\psi_1(x)= \frac{1}{\sqrt{k_1}} \left(A_\rightarrow e^{i k_1 x} + A_\leftarrow e^{-ik_1x}\right)\quad x<0$,

$\psi_2(x)= \frac{1}{\sqrt{k_2}} \left(B_\rightarrow e^{i k_2 x} + B_\leftarrow e^{-ik_2x}\right)\quad x>0$

where subscripts 1 and 2 denote the regions x < 0 and x > 0 respcetivley, the subscripts (→) and (←) on the amplitudes A and B denote the direction of the particle's velocity vector: right and left respectivley.

The coefficients 1√k₁ and 1√k₂ are normalization constants. The wave vectors related to the energy by

$k_1=\sqrt{2m E/\hbar^2}$,

$k_2=\sqrt{2m (E-V_0)/\hbar^2}$

both of which have the same form as the De Broglie relation (in one dimension)

$p=\hbar k$.

Boundary conditions

The coefficients A, B have to be found from the boundary conditions of the wave function at x = 0. The wave function and its derivative have to be continuous everywhere, so:

$\psi_1(0)=\psi_2(0)$,

$\frac{d}{dx}\psi_1(0) = \frac{d}{dx}\psi_2(0)$.

Inserting the wave functions, the boundary conditions give the following restrictions on the coefficients

$\sqrt{k_2}(A_\rightarrow+A_\leftarrow)=\sqrt{k_1}(B_\rightarrow+B_\leftarrow)$

$\sqrt{k_1}(A_\rightarrow-A_\leftarrow)=\sqrt{k_2}(B_\rightarrow-B_\leftarrow)$

Transmission and reflection

It is useful to compare the situation to the classical case. In both cases, the particle behaves as a free particle outside of the barrier region. A classical particle with energy E larger than the barrier height V₀ will be slowed down but never reflected by the barrier, while a classical particle with E < V₀ incident on the barrier from the left would always be reflected.

To study the quantum case, let us consider the following situation: a particle incident on the barrier from the left side A_→. It may be reflected (A_←) or transmitted B_→. Here and in the following assume E > V₀.

To find the amplitudes for reflection and transmission for incidence from the left, we set in the above equations A_→ = 1 (incoming particle), A_← = R (reflection), B_← = 0 (no incoming particle from the right) and B_→ = T (transmission). We then solve the set of two linear equations for T and R.

The result is:

$T=\frac{2\sqrt{k_1k_2}}{k_1+k_2}$

$R=\frac{k_1-k_2}{k_1+k_2}.$

The model is symmetric with respect to a parity transformation and at the same time interchange k₁ and k₂. For incidence from the right we have therefore the amplitudes for transmission and reflection

$T'=T=\frac{2\sqrt{k_1k_2}}{k_1+k_2}$

$R'=-R=\frac{k_2-k_1}{k_1+k_2}.$

Analysis of the expressions

Reflection and transmission probability at a Heaviside-step potential. Dashed: classical result. Solid lines: quantum mechanics. For E < V₀ the classical and quantum problem give the same result.

Energy less than step height (E < V₀)

For energies E < V₀, the wave function to the right of the step is exponentially decaying over a distance $1/(ik_2)$.

Energy greater than step height (E > V₀)

In this energy range the transmission and reflection coefficient differ from the classical case. They are the same for incidence from the left and right:

$T=|t|^2=|T'|^2=\frac{4k_1k_2}{(k_1+k_2)^2}$

$R=|r|^2=|R'|^2=1-T=\frac{(k_1-k_2)^2}{(k_1+k_2)^2}$

In the limit of large energies E ≫ V₀, we have k₁ ≈ k₂ and the classical result T = 1, R = 0 is recovered.

Thus there is a finite probability for a particle with an energy larger than the step height to be reflected.

Applications

The Heaviside step potential mainly serves as an exercise in introductory quantum mechanics, as the solution requires understanding of a variety of quantum mechanical concepts: wavefunction normalization, continuiuty, incident/reflection/transmission amplitudes, and probabilities.

A similar problem to the one considered appears in the physics of normal-metal superconductor interfaces. Quasiparticles are scattered at the pair potential which in the simplest model may be assumed to have a step-like shape. The solution of the Bogoliubov-de Gennes equation resembles that of the discussed Heaviside-step potential. In the superconductor normal-metal case this gives rise to Andreev reflection.

See also

• Rectangular potential barrier
• Finite potential well
• Infinite potential well
• Delta potential barrier (QM)
• Finite potential barrier (QM)

Sources

• Quantum Mechanics Demystified, D. McMahon, Mc Graw Hill (USA), 2006, ISBN(10-) 0-07-145546 9
• Quantum Physics of Atoms, Molecules, Solids, Nuclei, and Particles (2nd Edition), R. Eisberg, R. Resnick, John Wiley & Sons, 1985, ISBN 978-0-471-873730
• Quantum Mechanics, E. Abers, Pearson Ed., Addision Wesley, Prentice Hall Inc, 2004, ISBN 9780131461000
• Elementary Quantum Mechanics, N.F. Mott, Wykeham Science, Wykeham Press (Taylor & Francis Group), 1972, ISBN 0-85109-270-5
• Stationary States, A. Holden, College Physics Monographs (USA), Oxford University Press, 1971, ISBN 0-19-851121-3
• Quantum mechanics, E. Zaarur, Y. Peleg, R. Pnini, Schaum’s Oulines, Mc Graw Hill (USA), 1998, ISBN (10-) 007-0540187

Further reading

• The New Quantum Universe, T.Hey, P.Walters, Cambridge University Press, 2009, ISBN 978-0-521-56457-2.
• Quantum Field Theory, D. McMahon, Mc Graw Hill (USA), 2008, ISBN 978-0-07-154382-8
• Quantum mechanics, E. Zaarur, Y. Peleg, R. Pnini, Schaum’s Easy Oulines Crash Course, Mc Graw Hill (USA), 2006, ISBN (10-)007-145533-7 ISBN (13-)978-007-145533-6