We recast the well-known Numerov method for solving Schrödinger’s equation into a representation of the kinetic energy operator on a discrete lattice. With just a few lines of code in a high-level programming environment such as mathematica, it is simple to calculate and plot accurate eigenvalues and eigenvectors for a variety of potential problems. We illustrate the method by calculating high-accuracy solutions for the potential.
REFERENCES
1.
J. M.
Blatt
, “Practical points concerning the solution of the Schrödinger equation
,” J. Comput. Phys.
1
, 382
–396
(1967
).2.
P. C.
Chow
, “Computer solutions to the Schrödinger equation
,” Am. J. Phys.
40
, 730
–734
(1972
).3.
R. H.
Landau
, M. J.
Paez
, and C. C.
Bordeianu
, A Survey of Computational Physics
(Princeton U.P.
, Princeton, NJ
, 2008
).4.
M.
Belloni
and W.
Christian
, “Time development in quantum mechanics using a reduced Hilbert space approach
,” Am. J. Phys.
76
, 385
–392
(2008
).5.
J. J.
Sakurai
and J. J.
Napolitano
, Modern Quantum Mechanics
, 2nd ed. (Addison Wesley
, San Francisco
, 2010
).6.
A.
Gangopadhyaya
, J. V.
Mallow
, and C.
Rasinariu
, Supersymmetric Quantum Mechanics: An Introduction
(World Scientific
, Hackensack, NJ
, 2011
).7.
See supplementary material at http://dx.doi.org/10.1119/1.4748813 or <http://www.physics.wisc.edu/~tgwalker/NumerovExamples/> for the code of Fig. 1 and some worked examples.
8.
A.
Derevianko
, E.
Luc-Koenig
, and F.
Masnou-Seeuws
, “Application of B-splines in determining the eigenspectrum of diatomic molecules: Robust numerical description of halo-state and Feshbach molecules
,” Can. J. Phys.
87
, 67
–74
(2009
).9.
C. C.
Marston
and G. G.
Balint-Kurti
, “The Fourier grid Hamiltonian method for bound state eigenvalues and eigenfunctions
,” J. Chem. Phys.
91
, 3571
–3576
(1989
).© 2012 American Association of Physics Teachers.
2012
American Association of Physics Teachers
AAPT members receive access to the American Journal of Physics and The Physics Teacher as a member benefit. To learn more about this member benefit and becoming an AAPT member, visit the Joining AAPT page.
