Pseudospectral methods on a semi-infinite interval with application to the Hydrogen atom: a comparison of the mapped Fourier-sine method with Laguerre series and rational Chebyshev expansions

  • Authors:

  • John P. Boyd;C. Rangan;P. H. Bucksbaum

  • Affiliations:

  • Department of Atmospheric, Oceanic and Space Science, University of Michigan, 2455 Hayward Avenue, Ann Arbor, MI;Department of Physics, FOCUS Center, University of Michigan, Ann Arbor, MI;Department of Physics, FOCUS Center, University of Michigan, Ann Arbor, MI

  • Venue:
  • Journal of Computational Physics
  • Year:
  • 2003

Quantified Score

Hi-index 31.47

Visualization

Abstract

The Fourier-sine-with-mapping pseudospectral algorithm of Fattal et al. [Phys. Rev. E 53 (1996) 1217] has been applied in several quantum physics problems. Here, we compare it with pseudospectral methods using Laguerre functions and rational Chebyshev functions. We show that Laguerre and Chebyshev expansions are better suited for solving problems in the interval r ∈ [0, ∞] (for example, the Coulomb-Schrödinger equation), than the Fourier-sinemapping scheme. All three methods give similar accuracy for the hydrogen atom when the scaling parameter L is optimum, but the Laguerre and Chebyshev methods are less sensitive to variations in L. We introduce a new variant of rational Chebyshev functions which has a more uniform spacing of grid points for large r, and gives somewhat better results than the rational Chebyshev functions of Boyd [J. Comp. Phys. 70 (1987) 63].