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

Numerical solution of one dimensional Schrodinger equation

  1. Apr 18, 2014 #1
    I want to solve one dimensional Schrodinger equation for a scattering problem. The potential function is 1/ ( 1+exp(-x) ). So at -∞ it goes to 0 and at ∞ it's 1. The energy level is more than 1. I used Numerov's method and integrated it from +∞ (far enough) backwards with an initial value =1 . But I believe it's wrong b/c squared wave function is oscillating on whole interval and it's supposed to be constant after the jump in potential. I know that I'm doing somewhere wrong in my solution So I would appreciate you if you help me by this or introduce me some sources.
  2. jcsd
  3. Apr 19, 2014 #2
    Sounds like your solution is a linear combination of a particle incoming from the left and a particle incoming from the right. Interference of ##e^{ipx}## and ##e^{-ipx}## terms would produce this oscillating behavior. Maybe you can find two solutions and take a linear combination to eliminate the left-moving component on the right side of the potential jump?
  4. Apr 19, 2014 #3
    Actually we have two wave function A exp(ipx) and B exp(−ipx) at far left and one C exp(iqx) at far right. I need to write my program in a way which the coefficients A, B and C been determined self-consistently.
  5. Apr 19, 2014 #4
    Right, but the most general solution to the Schrodinger equation also has a D exp(-iqx) component at the far right. Are you doing anything to prevent this component from appearing in your numerical solution?
Share this great discussion with others via Reddit, Google+, Twitter, or Facebook