Physlets run in a Java-enabled browser on the latest Windows & Mac operating systems.

If Physlets do not run, click here for help on updating Java and setting Java security.

# Section 9.4: Plane Wave Scattering: Potential Energy Steps

Show Real and Imaginary Parts | Show Phase as Color

Please wait for the animation to completely load.

We begin our study of quantum-mechanical scattering by looking at a general case where we have two regions (Regions I and II) and a change in the potential energy at the interface between the two regions. Typically we choose *V* = 0 in Region I and *V* = *V*_{0} in Region II, and the change occurs at *x* = 0. This situation is shown in Step Up Animation 1 and Step Up Animation 2. The wave functions when *E* > *V*_{0} are

Ψ_{inc }= *A* exp{*i*(*k*_{1}*x* − ω*t*)} Ψ_{refl } = *B* exp{−*i*(*k*_{1}*x *+ ω*t*)} Ψ_{trans } = *C* exp{*i*(*k*_{2}x − ω*t*)}, (9.9)

where *k*_{1}^{2 }= 2*mE*/*ħ*^{2} and *k*_{2}^{2 }= 2*m*(*E *−* V*_{0})/*ħ*^{2}, while the momenta are just: *p*_{inc }= *ħk*_{1}, *p*_{refl } = −*ħk*_{1}, and *p*_{trans }= *ħk*_{2}. In addition, we note that the ω's are the same. Why? The ω's are related to the *total energy* of the plane wave. The *total energy * must be the same in each region (the entire wave function represents one state of the Schrödinger equation) even when the *k*'s are different in each region. Therefore, since *E* = *ħ*ω, the ω's must be the same. In addition, if the ω's were different, the solution to the boundary-value problem we will shortly consider would not necessarily be valid for all times.

To construct transmission and reflection coefficients, we must first calculate *J*_{x} for the incident, reflected, and transmitted waves from the equation^{4}

*J*_{x }= (*ħ*/2*mi*)**{**ψ* [∂ Ψ/∂*x*] − [∂ Ψ*/∂*x*]Ψ**} **. (9.12)

We find that

*J*_{inc }= (*ħk*_{1}/*m*) |*A*|^{2} *J*_{refl }= −(*ħk*_{1}/*m*) |*B*|^{2} *J*_{trans }= (*ħk*_{2}/*m*) |*C*|^{2} . (9.13)

For *E* > *V* we may use the above *parts of the wave function*, but we must match the wave function parts at the change in potential energy (the boundary at *x* = 0) since the parts in Eq. (9.7) represent just one wave function.^{5 }We require that ψ_{I}(0) = ψ_{II}(0) and that ψ_{I}'(0) = ψ_{I}'(0). We can use time-independent wave functions (since the time dependence for all wave function parts is identical, it cancels):

Ψ_{inc }= *A* exp(*ik*_{1}*x*) Ψ_{refl } = *B* exp(−*ik*_{1}*x*) Ψ_{trans} = *C* exp(*ik*_{2}*x*),

where Ψ_{I }= Ψ_{inc }+ Ψ_{refl}. When we match Ψ_{I}(0) = Ψ_{II}(0), this gives *A* + *B* = *C*. When we match Ψ_{I}'(0) = Ψ_{I}'(0), we get *k*_{1}*A* −* k*_{1}*B* = *k*_{2}*C*. From these conditions we solve for the fractions |*B*/*A*|^{2} and |*C*/*A*|^{2}

|*B*/*A*|^{2 }= (1 − *k*_{2}/*k*_{1})^{2}/(1 + *k*_{2}/*k*_{1})^{2} and |*C*/*A*|^{2 }= 4/ (1 + *k*_{2}/*k*_{1})^{2},

since they occur in the the current densities, Eq.(9.12), which are necessary in order to calculate the transmission and reflection coefficients from Eq.(9.10). We therefore find that

*R* = (1 − *k*_{2}/*k*_{1})^{2}/(1 + *k*_{2}/*k*_{1})^{2} (9.14)

and

*T* = (*k*_{2}/*k*_{1}) [4/(1 + *k*_{2}/*k*_{1})^{2}]. (9.15)

Recall *k*_{1}^{2 }= 2*mE*/*ħ*^{2} and *k*_{2}^{2 }= 2*m*(*E* − *V*_{0})/*ħ*^{2}. Therefore, *k*_{2}^{2}/*k*_{1}^{2} = (*E* − *V*_{0})/*E* = 1 − *V*_{0}/*E*. As *E* increases, *k*_{2 }→ *k*_{1}, and therefore T → 1, as expected. With a bit of algebra, you can also convince yourself that *T* + *R* = 1 by adding Eqs. (9.14) and (9.15) together.

If we had the opposite circumstance, as shown in Step Down Animation 1 and Step Down Animation 2, with an incident plane wave subject to a potential energy *V*_{0} and then at *x *> 0 experiencing no potential energy, this would amount to *k*_{2 }→ *k*_{1} and *k*_{1 }→ *k*_{2}. Multiplying this result by the square of the ratio of the *k*'s gets us back to the original *T* and *R* in Eqs. (9.14) and (9.15). Therefore the transmission and reflection coefficients for these two different problems are actually the same.

Look at both Step Up Animation 1 and Step Up Animation 2. Step Up Animation 1 shows you the quantum-mechanical plane wave in both regions and the top panel shows the energy diagram. Change the energy of the incident plane wave by moving the slider and note the effect on the plane wave and also the transmission and reflection coefficients. Step Up Animation 2 shows the same situation but with a graph of the transmission and reflection coefficients as a function of energy. Click-drag the plane wave and drag up to increase the energy or drag down to decrease the energy. If you do this slowly enough, you can create a very nice graph of the transmission and reflection coefficients versus *E* for a fixed *V*_{0}.

For *E* < *V*_{0}, we can use the time-independent wave functions

Ψ_{inc }= *A* exp(*ik*_{1}*x*) Ψ_{refl }= *B* exp(_{ }−*ik*_{1}*x*) Ψ_{trans} = *C* exp(_{ }−κ_{2}*x*),

where *k*_{1}^{2 }= 2*mE*/*ħ*^{2} and κ_{2}^{2 }= 2*m*(|*E* − *V*_{0}|)/*ħ*^{2}. We find for this case

*J*_{inc}= (*ħk*_{1}/*m*) |*A*|^{2} *J*_{refl }= _{ }−(*ħk*_{1}/*m*) |*B*|^{2} *J*_{trans }= 0.

Now we must again match the parts of the wave function. We require that Ψ_{I}(0) = Ψ_{II}(0) and that Ψ_{I}'(0) = Ψ_{I}'(0). When we match Ψ_{I}(0) = Ψ_{II}(0), this gives *A* + *B* = *C*. When we match Ψ_{I}'(0) = Ψ_{I}'(0), we get *ik*_{1}*A*_{ } − *ik*_{1}*B* =_{ }− κ_{2}*C*. We eventually arrive at the relationship

|*B*/*A*|^{2 }= [(*k*_{1 }+ *i*κ_{2})/(*k*_{1 } − *i*κ_{2})] / [(*k*_{1 }− *i*κ_{2})(*k*_{1 }+ *i*κ_{2})] = 1.

Note in this case, *R* = |*B*/*A*|^{2 }= 1. From *R* = 1 we also expect *T* = 0, which is also clear from the fact that *J*_{trans}= 0. In addition, if *E* = *V*_{0}, we find that *T* = 0. These results agree with the classical scattering case *except* that there is a finite probability of finding the particle in the *classically-forbidden* region.

^{4}You may be wondering why we care about the probability current density. The conservation of probability in its integral form in one dimension is

*dP*/*dt* + ∫ (*dJ*_{x}/*dx*) *dx* = 0 [integral from −∞ to +∞] ,

and for a scattering process, the total probability of finding the particle over all space does not change with time, or that *dP*/*dt* = 0. Therefore: ∫ (*dJ*_{x}/*dx*) *dx* = 0 or that *J*_{x}(∞) = *J*_{x}(−∞). We can identify *J*_{x}(∞) ≡ *J*_{trans }and *J*_{x}(−∞) ≡ *J*_{inc}+ *J*_{refl}. We find that *J*_{trans} = *J*_{inc}+ *J*_{refl} and therefore *J*_{inc}= *J*_{trans} − *J*_{refl} , which leads directly to

1 = *J*_{trans}/*J*_{inc} − *J*_{refl} /*J*_{inc}.

The sign in front of the *J*_{refl} term is due to the fact that this expression is always *negative*, therefore −*J*_{refl} will be *positive*. We incorporate this sign in the definition of the reflection and transmission coefficients as

*T* = *J*_{trans}/*J*_{inc} ≡ |*J*_{trans}/*J*_{inc}| (9.10)

and

*R* = −*J*_{refl} /*J*_{inc} ≡ |*J*_{refl} /*J*_{inc}|. (9.11)

^{5}Also look for the *E* = *V*_{0} case in the animations. The solution in Region II should be a straight line as discussed.

« previous

next »