Step Potentials


The transmission and reflection coefficients are defined in terms of the probability current, J, also known as current density

qm13-eq-01.gif (1281 bytes)

Consider a one-dimensional a beam of particles that are plane wave states. The beam is moves in x-direction and impinges on a barrier as shown in Fig. 1 below 

qm13-im-01.gif (3309 bytes)

 

The transmission coefficient T and reflection coefficient R are defined in terms of the components of the probability current. Let Ji represent the incoming current, Jt represent the incoming transmitted current, and Let Jr represent the reflected current

qm13-eq-02.gif (1150 bytes)


General Case of Step Potential for E > V

 The magnitude of the momentum of the particles in the incoming, transmitted and reflected beams are, respectively

 qm13-eq-03.gif (1351 bytes)

where, in general, k1 k2. The potential is constant for x < 0 and for x > 0 so the wavefunction in each domain is that of a free particle. Since energy is conserved it follows that Ei = Et = Er . The wavefunction for the incoming, transmitted and reflected beams are, respectively

qm13-eq-04.gif (2555 bytes) 

Using Eq. (4a) in the relation for the current in Eq. (1) gives the x-component of the current density

 qm13-eq-05.gif (1457 bytes)

 Substituting Eq. (4a) into Eq. (5)

 qm13-eq-06.gif (5047 bytes)

 A similarly procedure for the transmitted and reflected current densities

qm13-eq-07.gif (1701 bytes)

We these values for the currents the transmission and reflection coefficients become

qm13-eq-08.gif (1476 bytes)

The time-independent Schrodinger equation is

qm13-eq-09.gif (1988 bytes)

For the step potential shown in Fig. 1 can broken up into two domains. Region I is the domain x < 0 and Region II is the domain x > 0.  The wavefunction will be partitioned into two parts as well corresponding to each domain. In Region I the Hamiltonian consists entirely of the kinetic energy term so Eq. (9b) reduces to 

qm13-eq-10.gif (1878 bytes) 

In Region II the Hamiltonian has a kinetic energy term plus a constant potential 

qm13-eq-11.gif (2227 bytes)

The solution to Eq. (9) and (10) is

qm13-eq-12.gif (1638 bytes) 

Since the time-independent Schrodinger equation is second order there must be two boundary conditions. First note that we can set since it is the coefficient of a wave coming from +infinity and since there is no such wave we can set D = 0. They are determined by the condition that the wavefunction and its first derivative must be continuous. Applying this criteria at x = 0 means

qm13-eq-13.gif (1571 bytes) 

Taking the derivative of Eq. (12) gives

qm13-eq-14.gif (1874 bytes) 

Applying the boundary condition in Eq. (13a) gives

 qm13-eq-15.gif (1250 bytes)

Appling the boundary condition in Eq. (12b) gives

 qm13-eq-16.gif (1553 bytes) 

Eqs. (15) and (16) can be written simply as 

 qm13-eq-17.gif (1651 bytes) 

Eq. (17) can be solved for the ratios B1/A1 and C1/A1 to give

qm13-eq-18.gif (1609 bytes) 

The transmission coefficient can now be written as 

qm13-eq-19.gif (1793 bytes)


General Case of Step Potential for E < V  

In Region I we have the same situation as above and therefore the same solution as in Eq. (12a). The relationship for the reflection coefficient will therefore have the same relation as above, i.e. R will still be given as in Eq. (8b).   However in Region II the quantity E – V is negative. Therefore

qm13-eq-20.gif (2205 bytes)

There are two solutions to Eq. (20), one of which is an exponentially increasing function. This must be discarded since it leads to a diverging probability density. Therefore the solution is

qm13-eq-21.gif (1553 bytes) 

The derivatives are

qm13-eq-22.gif (1863 bytes) 

Just as above, applying the boundary condition defined in Eq. (13a) gives

qm13-eq-23.gif (1251 bytes) 

Applying the boundary condition defined in Eq. (13b) gives

qm13-eq-24.gif (1556 bytes) 

Eqs. (23) and (24) can be written simply as

qm13-eq-25.gif (1558 bytes) 

Eq. (25) can be solved for the ratios B2/A2 and C2/A2 to give 

qm13-eq-26.gif (1641 bytes) 

Since the right side of Eq. (26) is of the form z*/z and therefore | z*/z| = 1. This gives

qm13-eq-27.gif (1216 bytes)

The transmission coefficient is found by using the relation R + T = 1. It follows that T = 0. This means that the transmitted current density is zero. The wavefunction in Region II is multiplied by a decreasing exponential. The wavefunction therefore does not represent a propagating wave. Such a wave is called an evanescent wave.


Rectangular Barrier and Tunneling: Case E > V 

Consider now a rectangular barrier for which E > V as shown in Fig. 2  

qm13-im-03.gif (4243 bytes)

 Just as above, in each region the potential is constant. The solutions are again similar to that in E. (4) above. We have

 qm13-eq-28.gif (2519 bytes)

 Conservation of energy implies

 qm13-eq-29.gif (1637 bytes)

The boundary conditions are

 qm13-eq-30.gif (1410 bytes) 

qm13-eq-31.gif (1741 bytes)

Applying the boundary conditions in Eq. (30) gives

 qm13-eq-32.gif (2046 bytes)

 To apply the boundary condition in Eq. (31) we need to first take the derivative of each term in Eq. (30)

 qm13-eq-33.gif (2323 bytes) 

Applying the boundary conditions in Eq. (31) gives

 qm13-eq-34.gif (2665 bytes)

 Dividing Eq. (32) by A gives

 qm13-eq-35.gif (2481 bytes)

Dividing Eq. (34) by A gives

 qm13-eq-36.gif (2784 bytes)

Summary of results

qm13-eq-37.gif (2747 bytes)

After some tedious algebra these equations can be put into the form

qm13-eq-38.gif (2877 bytes)

From a calculation similar to Eq. (7) it can be shown that the transmission and reflection currents are

qm13-eq-39.gif (1767 bytes)

Therefore the transmission and reflection currents are, respectively

qm13-eq-40.gif (1730 bytes)

Substituting Eq. (38a) into 1/T gives, using the values or the energy in Eq. (28)

qm13-eq-41.gif (1405 bytes)

To find R use the relation  T + R = 1 to give R = 1 – T.


Hosted by www.Geocities.ws

1