Reversible adaptive regularization: perturbed Kepler motion and classical atomic trajectories

Open Access
Publisher Website
Google Scholar
Abstract
Reversible and adaptive integration methods based on Kustaanheimo–Stiefel regularization and modified Sundman transformations are applied to simulate general perturbed Kepler motion and to compute classical trajectories of atomic systems (e.g. Rydberg atoms). The new family of reversible adaptive regularization methods also conserve angular momentum and exhibit superior energy conservation and numerical stability in long–time integrations. The schemes are appropriate for scattering, for astronomical calculations of escape time and long–term stability, and for classical and semiclassical studies of atomic dynamics. The components of an algorithm for trajectory calculations are described. Numerical experiments illustrate the effectiveness of the reversible approach.

This publication has 38 references indexed in Scilit: