Dismiss Notice
Join Physics Forums Today!
The friendliest, high quality science and math community on the planet! Everyone who loves science is here!

Transmission and Reflection Coefficients

  1. Nov 14, 2006 #1
    Consider the step potential defined by
    V(x) = 0 if x < 0
    Vb > 0 if x=> 0

    a) Compute te relfection and tranmission coefficients for a particle of energy E > Vb approaching from th left

    For x < 0
    Schrodinger equaion since V = 0
    [tex] -\frac{\hbar^2}{2m} \frac{\partial^2 \psi}{\partial x^2} = E \psi(x) [/tex]

    [tex] \frac{\partial^2 \psi}{\partial x^2} = -\frac{2mE}{\hbar^2} \psi(x) [/tex]

    define [tex] k_{1}^2 = \frac{2mE}{\hbar^2} [/tex]

    [tex] \frac{\partial^2 \psi}{\partial x^2} = -k^2 \psi(x) [/tex]

    [tex] \psi(x) = A \cos(k_{1}x) + B \sin(k_{1} x) [/tex]

    for x=> 0
    [tex] -\frac{\hbar^2}{2m} \frac{\partial^2 \psi}{\partial x^2} = (E-V_{B}) \psi(x) [/tex]
    [tex] k_{2}^2 = \frac{2m}{\hbar^2} (E - V_{B}) [/tex]
    [tex] \frac{\partial^2 \psi}{\partial x^2} = k_{2}^2 \psi(x) [/tex]
    [tex] \psi_(x) = C \exp(ik_{2} x) + D\exp(-ik_{2}x) [/tex]
    but we only care about the one with the negative exponential
    hence C = 0
    [tex] \psi_(x) = D\exp(-ik_{2}x) [/tex]

    are these solutions correct? Im just afraid i may have made a stupid mistake by switching the negative sign somewhere :(

    ok to proceed i apply the boundary conditions s.t.

    [tex] \psi_{I}(0) = \psi_{T}(0) [/tex]
    [tex] \psi_{I}'(0) = \psi_{T}'(0) [/tex]

    we calculate probability current density j

    [tex] j_{trans} = \frac{\hbar k_{2}}{m} |D|^2 [/tex]
    [tex] j _{inc} = \frac{\hbar k_{1}}{m} |A|^2 [/tex]
    [tex] j_{refl} = \frac{\hbar k_{1}}{m} |B|^2 [/tex]

    and transmission coefficient is calculated like this

    [tex] T = \left| \frac{j_{trans}}{j_{inc}} \right| [/tex]
    [tex] R = \left| \frac{j_{refl}}{j_{inc}} \right| [/tex]

    thank you for your input
  2. jcsd
  3. Nov 15, 2006 #2


    User Avatar
    Science Advisor
    Homework Helper
    Gold Member

    You want a wave traveling to the right so you must keep the C wavefunction and set D =0.

  4. Jul 4, 2010 #3
    Sorry to resurrect a long dead thread, but if i wanted to impose the condition D = 0, ie incoming waves only from one direction on standing wave computed numerically, how would i go about doing that? I know Fourier analysis can give me the momentum components, but thats about it. Would i have to do it in the time dependent case?
  5. Jul 5, 2010 #4


    User Avatar
    Science Advisor
    Homework Helper
    Gold Member

    Say you have a wave incident from the left. The wavefunction would be

    ψI = A*Exp(ik1x)+B*Exp(-ik1x) in region I (x ≤ 0)
    ψII = C*Exp(ik2x) in region II (x ≥ 0)

    You have an incident and reflected wave in region I (constants A and B respectively) and a transmitted wave in region II. Of course, you have to impose the proper boundary conditions at x=0.

    For a lovely simulation of this and potential barriers in general, you may wish to go to
  6. Jul 5, 2010 #5
    I understand that perfectly from an analytic point of view, but what i've got is the equation (actually want to do the KG equation but good to start with Schrod) with a set of boundary conditions, mass of the particle and potential and thats it. I then solve numerically and get a standing wave. It is in this process that i am unsure as to how to impose the correct condition to kill the other incoming wave. In fact, i'm really lost because essentially mathematica will spit out what are essentially shifted cos waves in either regime, which when you take a fourier transform will have momentum components in both directions, both with equal magnitude (as you expect from a cosine/sine wave). I thought perhaps a fourier decomposotion would pull out the momentum components and i could kill them off accordingly but it didnt work really :(

Share this great discussion with others via Reddit, Google+, Twitter, or Facebook