Delta potential barrier (QM)

From Wikipedia, the free encyclopedia

The Delta potential barrier is a textbook problem of quantum mechanics. The problem consists of solving the time-independent Schrödinger equation for a particle in a delta function potential in one dimension.

Scattering at a delta function potential of strength λ. The amplitudes and direction of left and right moving waves are indicated. In red, those waves used for the derivation of the reflection and transmission amplitude.
Scattering at a delta function potential of strength λ. The amplitudes and direction of left and right moving waves are indicated. In red, those waves used for the derivation of the reflection and transmission amplitude.

[edit] Calculation

The time-independent Schrödinger equation for the wave function ψ(x) reads

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

where H is the Hamiltonian, \hbar is the (reduced) Planck constant, m is the mass, E the energy of the particle and

V(x)=\frac{\lambda}{m}\delta(x)

is the delta function potential (barrier) with strength λ > 0. Here we have chosen the potential to be at the origin. Without changing the results, any other shifted position was possible. The first term in the Hamiltonian, -\frac{\hbar^2}{2m}\frac{d^2}{dx^2}\psi is the kinetic energy.

The barrier divides the space in two parts (x < 0,x > 0). In any of these parts the particle is 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_L(x)= A_r e^{i k x} + A_l e^{-ikx}\quad x<0, and
\psi_R(x)= B_r e^{i k x} + B_l e^{-ikx}\quad x>0

where the wave vector is related to the energy via

k=\sqrt{2m E}/\hbar . The index r/l on the coefficients A and B denotes the direction of the velocity vector. Those coefficients have to be found from the boundary conditions of the wave function at x = 0:
ψL = ψR,

\frac{d}{dx}\psi_L = \frac{d}{dx}\psi_R - \frac{2\lambda}{\hbar^2} \psi_R.

The second of these equations follows from integrating the Schrödinger equation with respect to x. The boundary conditions thus give the following restrictions on the coefficients

Ar + Al = Br + Bl
ik(A_r-A_l-B_r+B_l)=-\frac{2\lambda}{\hbar^2}(B_r+B_l).

[edit] Transmission and reflection

At this point, it is instructive to compare the situation to the classical case. In both cases, the particle behaves as a free particle outside of the barrier region. However, a classical particle having a finite energy cannot pass the infinitely high potential barrier and will be reflected from it. To study the quantum case, let us consider the following situation: a particle incident on the barrier from the left side (Ar). It may be reflected (Al) or transmitted (Br).

To find the amplitudes for reflection and transmission for incidence from the left, we put in the above equations Ar = 1 (incoming particle), Al = r (reflection), Bl=0 (no incoming particle from the right) and Br = t (transmission) and solve for r,t. The result is:

t=\frac{1}{\frac{i\lambda}{\hbar^2k}+1}
r=\frac{1}{\frac{i\hbar^2 k}{\lambda}-1}.

Due to the mirror symmetry of the model, the amplitudes for incidence from the right are the same as those from the left. The surprising result, from the classical point of view, is that there is a non-zero probability (given by the transmission coefficient)

T=|t|^2=\frac{1}{1+\frac{\lambda^2}{\hbar^4k^2}}= \frac{1}{1+\frac{\lambda^2}{2m\hbar^2 E}}

for the particle to be transmitted through the barrier. This effect which differs from the classical case is called quantum tunneling.

For completeness, the probability for reflection (given by the reflection coefficient) is:

R=|r|^2=1-T=\frac{1}{1+\frac{\hbar^4k^2}{\lambda^2}}= \frac{1}{1+\frac{2m\hbar^2 E}{\lambda^2}}.

[edit] Remarks, Application

The calculation presented above may at first seem unrealistic and hardly useful. However it has proved to be a suitable model for a variety of real-life systems. One such example regards the interfaces between two conducting materials. In the bulk of the materials, the motion of the electrons is quasi free and can be described by the kinetic term in the above Hamiltonian with an effective mass m. Often the surfaces of such materials are covered with oxide layers or are not ideal for other reasons. This thin, non-conducting layer may then be modeled by a local delta-function potential as above. Electrons may then tunnel from one material to the other giving rise to a current.

The operation of a scanning tunneling microscope (STM) relies on this tunneling effect. In that case, the barrier is due to the air between the tip of the STM and the underlying object. The strength of the barrier is related to the separation being stronger the further apart the two are. For a more general model of this situation, see Finite potential barrier (QM). The delta function potential barrier is the limiting case of the model considered there for very high and narrow barriers.

The above model is one-dimensional while the space around us is three-dimensional. So in fact one should solve the Schrödinger equation in three dimensions. On the other hand, many systems only change along one coordinate direction and are translationally invariant along the others. The Schrödinger equation may then be reduced to the case considered here by an Ansatz for the wave function of the type: Ψ(x,y,z) = ψ(x)φ(y,z).