The Adaptive Verlet Method

We discuss the integration of autonomous Hamiltonian systems via dynamical rescaling of the vector field (reparameterization of time). Appropriate rescalings (e.g., based on normalization of the vector field or on minimum particle separation in an N-body problem) do not alter the time-reversal symmetry of the flow, and it is desirable to maintain this symmetry under discretization. For standard form mechanical systems without rescaling, this can be achieved by using the explicit leapfrog--Verlet method; we show that explicit time-reversible integration of the reparameterized equations is also possible if the parameterization depends on positions or velocities only. For general rescalings, a scalar nonlinear equation must be solved at each step, but only one force evaluation is needed. The new method also conserves the angular momentum for an N-body problem. The use of reversible schemes, together with a step control based on normalization of the vector field (arclength reparameterization), is demonstrated in several numerical experiments, including a double pendulum, the Kepler problem, and a three-body problem.

[1]  Robert D. Skeel,et al.  Symplectic Integration with Variable Stepsize , 1994 .

[2]  Tanner,et al.  Classical mechanics of two-electron atoms. , 1993, Physical review. A, Atomic, molecular, and optical physics.

[3]  Jack Wisdom.  Some Aspects of Chaotic Behavior in the Solar System , 1988 .

[4]  R. Ruth.  A Can0nical Integrati0n Technique , 1983, IEEE Transactions on Nuclear Science.

[5]  G. J. Cooper.  Stability of Runge-Kutta Methods for Trajectory Problems , 1987 .

[6]  E. Hairer,et al.  Accurate long-term integration of dynamical systems , 1995 .

[7]  Archie E. Roy,et al.  Predictability, stability, and chaos in N-body dynamical systems , 1991 .

[8]  Robert D. Skeel,et al.  Explicit canonical methods for Hamiltonian systems , 1992 .

[9]  Mari Paz Calvo,et al.  The Development of Variable-Step Symplectic Integrators, with Application to the Two-Body Problem , 1993, SIAM J. Sci. Comput..

[10]  J. C. Simo,et al.  On the stability of symplectic and energy-momentum algorithms for non-linear Hamiltonian systems with symmetry , 1996 .

[11]  Robert D. Skeel,et al.  Does variable step size ruin a symplectic integrator , 1992 .

[12]  J. M. Sanz-Serna,et al.  Numerical Hamiltonian Problems , 1994 .

[13]  E. Hairer.  Backward analysis of numerical integrators and symplectic methods , 1994 .

[14]  Linda R. Petzold,et al.  Numerical solution of initial-value problems in differential-algebraic equations , 1996, Classics in applied mathematics.

[15]  S. Mikkola,et al.  A Numerical Experimenter’s View of the Few-Body Problem , 1988 .

[16]  Yosi Shibberu.  Time-discretization of Hamiltonian dynamical systems☆ , 1994 .

[17]  S. Aarseth.  Integration Methods for Small N-Body Systems , 1988 .

[18]  Gene H. Golub,et al.  Matrix computations , 1983 .

[19]  P. Hut,et al.  Building a better leapfrog , 1995 .