Numerov Solutions for Single- and Double-Well Potentials

In quantum mechanics, energy states of a one-particle system are given by solutions of the time-independent Schrödinger equation
where is the potential energy and is an energy eigenfunction. In principle, one can always find a which satisfies this equation, however, not all are valid solutions since the boundary conditions must also be satisfied. In fact, it is imposition of the boundary conditions on which gives rise to a discrete spectrum of energy states (for a bound system). In the examples here, we use the Numerov approach to generate solutions of the time-independent Schrödinger equation for bound systems. In the examples here, as , thus must vanish as becomes large. As shown here, only certain energies will give rise to solutions that satisfy this condition. The Numerov solution will always diverge as becomes large; however, by carefully adjusting the energy until it is very close to the exact eigenvalue, the point of divergence can be pushed to greater and greater values of .
In the first example, we consider a particle in a harmonic potential . With energy units such that , the exact energy eigenvalues are at . Using the slider, one can scan the energy and find that at only for these discrete values will the computed solutions obey the boundary conditions.
In the second example, we consider the case of a particle in a double-well potential which has two states below the tunneling barrier. The slider is initially at the energy corresponding to the ground state. By carefully adjusting the slider, one can determine the energy of the first excited state, which shows a splitting of the energy levels due to tunneling.



  • [Snapshot]
  • [Snapshot]
  • [Snapshot]
  • [Snapshot]


Snapshot 1: ground state of double-well potential
In this case quantum mechanical tunneling produces a delocalization of the wavefunction between the two wells. One can clearly see that a considerable portion of is "underneath the barrier".
Snapshot 2: Numerov solution for the first excited state of a harmonic oscillator at
Notice that stationary solutions occur only when when obeys the boundary conditions that it vanish as .
In the next two snapshots, when is slightly above or below this value, diverges as becomes large.
For an elementary overview of quantum mechanics and the use of the Numerov approach to compute stationary solutions, the interested viewer is referred to I. Levine, Quantum Chemistry, 5th ed., Upper Saddle River, New Jersey: Prentice Hall, 1999.


    • Share:

Embed Interactive Demonstration New!

Just copy and paste this snippet of JavaScript code into your website or blog to put the live Demonstration on your site. More details »

Files require Wolfram CDF Player or Mathematica.

Mathematica »
The #1 tool for creating Demonstrations
and anything technical.
Wolfram|Alpha »
Explore anything with the first
computational knowledge engine.
MathWorld »
The web's most extensive
mathematics resource.
Course Assistant Apps »
An app for every course—
right in the palm of your hand.
Wolfram Blog »
Read our views on math,
science, and technology.
Computable Document Format »
The format that makes Demonstrations
(and any information) easy to share and
interact with.
STEM Initiative »
Programs & resources for
educators, schools & students.
Computerbasedmath.org »
Join the initiative for modernizing
math education.
Step-by-Step Solutions »
Walk through homework problems one step at a time, with hints to help along the way.
Wolfram Problem Generator »
Unlimited random practice problems and answers with built-in step-by-step solutions. Practice online or make a printable study sheet.
Wolfram Language »
Knowledge-based programming for everyone.
Powered by Wolfram Mathematica © 2018 Wolfram Demonstrations Project & Contributors  |  Terms of Use  |  Privacy Policy  |  RSS Give us your feedback
Note: To run this Demonstration you need Mathematica 7+ or the free Mathematica Player 7EX
Download or upgrade to Mathematica Player 7EX
I already have Mathematica Player or Mathematica 7+