Bound-State Solutions of the Schrödinger Equation by Numerical Integration
This Demonstration shows the mathematical solution of the time-independent Schrödinger equation for four potentials, the harmonic oscillator , the V-shaped potential , the anharmonic oscillator , and a square-well potential. The wavefunction is arbitrarily fixed at . You can obtain linearly independent solutions by numerical integration for different values of the derivative and the energy level. The vertical dashed lines indicate the locations of the classical turning points.
For any value of ,it is always possible to tune in a value of such that goes to zero either in the limit as or in the limit as . The goal is to find the discrete values of (the eigenvalues of the Hamiltonian) such that an acceptable solution (the eigenfunction) exists that goes to zero for both limits and .
For instance, the eigenvalues of the harmonic oscillator are , , , , …., while the first five eigenvalues corresponding to the V-shaped potential are , , , , and .
Units are such that , , and the proportionality coefficients in the potential functions are as indicated.