Abstract
We present a method whereby the 3D dynamics of the electronic wave packet in a hydrogen atom can be calculated efficiently. The method is constructed so as to satisfy the following two requirements: the wave function is zero at the Coulomb singular point so that the numerical difficulties concerning the singularity are avoided; the coordinate system is chosen so that the differential operators in the Hamiltonian can be well evaluated by the 3-point finite difference formula even near the singular point. The generalized cylindrical coordinate system (ξλ, z, φ) is introduced to satisfy the above conditions, and the value of λ is determined to be 3/2. The Schrödinger equation is discretized in time and space and solved by the Peaceman-Rachford method. The λ = 3/2 coodinate system helps the reduction in the number of grid points. To examine the numerical stability and accuracy of our method, we first apply it to cases where no laser field is turned on. The errors for the ordinary cylindrical coordinate system (λ = 1) are more than ten times as large as those for λ = 3/2. We then apply the method to the case where the atom interacts with a linearly polarized strong laser pulse. Our method is highly reliable and is a powerful tool for analyzing ionization processes of the hydrogen atom.
Original language | English |
---|---|
Pages (from-to) | 148-159 |
Number of pages | 12 |
Journal | Journal of Computational Physics |
Volume | 130 |
Issue number | 1 |
DOIs | |
Publication status | Published - 1997 Jan 1 |
ASJC Scopus subject areas
- Numerical Analysis
- Modelling and Simulation
- Physics and Astronomy (miscellaneous)
- Physics and Astronomy(all)
- Computer Science Applications
- Computational Mathematics
- Applied Mathematics