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 product of V₀, the height of the barrier, and the Heaviside step function:

$V(x) = \begin{cases} 0, & x < 0 \\ V_0, & x \ge 0 \end{cases}$

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) with the wave vectors in the respective regions being

$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

$(A_\rightarrow+A_\leftarrow)=(B_\rightarrow+B_\leftarrow)$

$k_1(A_\rightarrow-A_\leftarrow)=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. Once we have found the quantum-mechanical result we will return to the question of how to recover the classical limit.

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 k₁/k₂) (transmission). We then solve for T and R.

The result is:

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

$\sqrt{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

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

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

Analysis of the expressions

File:Steppotdiag.svg

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/(k_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'|=\frac{4k_1k_2}{(k_1+k_2)^2}$

$R=|R'|=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.

Classical limit

The result obtained for R depends only on the ratio E/V₀. This seems superficially to violate the correspondence principle, since we obtain a finite probability of reflection regardless of the value of Planck's constant or the mass of the particle. For example, we seem to predict that when a marble rolls to the edge of a table, there can be a large probability that it is reflected back rather than falling off. Consistency with classical mechanics is restored by eliminating the unphysical assumption that the step potential is discontinuous. When the step function is replaced with a ramp that spans some finite distance w, the probability of reflection approaches zero in the limit $wk\rightarrow0$, where k is the wavenumber of the particle.

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, continuity, 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-87373-0
• Quantum Mechanics, E. Abers, Pearson Ed., Addison Wesley, Prentice Hall Inc, 2004, ISBN 978-0-13-146100-0
• 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-1.
• 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 Outlines Crash Course, Mc Graw Hill (USA), 2006, ISBN (10-)007-145533-7 ISBN (13-)978-007-145533-6