Gradient Symplectic Algorithms for Solving the Radial Schrodinger Equation
arXiv:physics/0509229 · doi:10.1063/1.2150831
Abstract
The radial Schrodinger equation for a spherically symmetric potential can be regarded as a one dimensional classical harmonic oscillator with a time-dependent spring constant. For solving classical dynamics problems, symplectic integrators are well known for their excellent conservation properties. The class of {\it gradient} symplectic algorithms is particularly suited for solving harmonic oscillator dynamics. By use of Suzuki's rule for decomposing time-ordered operators, these algorithms can be easily applied to the Schrodinger equation. We demonstrate the power of this class of gradient algorithms by solving the spectrum of highly singular radial potentials using Killingbeck's method of backward Newton-Ralphson iterations.
19 pages, 10 figures
References in corpus (7)
- Gradient Symplectic Algorithms for Solving the Schroedinger Equation with Time-Dependent Potentials
- On the construction of high-order force gradient algorithms for integration of motion in classical and quantum systems
- Fourth Order Gradient Symplectic Integrator Methods for Solving the Time-Dependent Schrödinger Equation
- Fourth-Order Algorithms for Solving the Imaginary Time Gross-Pitaevskii Equation in a Rotating Anisotropic Trap
- Calculation of the spiked harmonic oscillators through a generalized pseudospectral method
- Higher Order Force Gradient Symplectic Algorithms
- Forward Symplectic Integrators and the Long Time Phase Error in Periodic Motions