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.