A P-stable exponentially fitted method is developed in this paper for the numerical integration of the Schrödinger equation. An application to the bound-states problem (we solve the radial Schrödinger equation in order to find eigenvalues for which the wavefunction and its derivative are continuous
Exponential fitted methods for the numerical integration of the Schrödinger equation
✍ Scribed by T.E. Simos
- Publisher
- Elsevier Science
- Year
- 1992
- Tongue
- English
- Weight
- 500 KB
- Volume
- 71
- Category
- Article
- ISSN
- 0010-4655
No coin nor oath required. For personal study only.
✦ Synopsis
A new sixth-order Runge-Kutta type method is developed for the numerical integration of the one-dimensional Schrodinger equation. The formula developed contains certain free parameters which allows it to be fitted automatically to exponential functions. We give a comparative error analysis with other sixth-order exponentially fitted methods. The theoretical and numerical results indicate that the new method is more accurate than the other exponentially fitted methods.
📜 SIMILAR VOLUMES
This algorithm uses a high-order, variable step Runge-Kutta like method in the region where the potential term dominates, and an exponential or Bessel fitted method in the asymptotic region. This approach can be used to compute scattering phase shifts in an efficient and reliable manner. A Fortran p
A new integral equation method for the numerical solution of the radial Schrödinger equation in one dimension, developed by the authors (1997, J. Comput. Phys. 134, 134), is extended to systems of coupled Schrödinger equations with both positive and negative channel energies. The method, carried out
By using the methods of perturbation theory it is possible to construct simple formulae for the numerical integration of the Schrodinger equation, and also to calculate expectation values solely by means of simple eigenvalue calculations. 1. The basic theory may be a multiplicative operator instead