Mollified Impulse Methods for Highly Oscillatory Differential Equations

We introduce a family of impulselike methods for the integration of highly oscillatory second-order differential equations whose forces can be split into a fast part and a slow part. Methods of this family are specified by two weight functions $\phi$, $\psi$; one is used to average positions and the other to mollify the force. When the fast forces are conservative and $\phi=\psi$, the methods here coincide with the mollified impulse methods introduced by Garcia-Archilla, Sanz-Serna, and Skeel. On the other hand, the methods here extend to nonlinear situations a well-known class of exponential integrators introduced by Hairer and Lubich for cases of linear fast forces. A convergence analysis is presented that provides insight into the role played by the processes of mollification and averaging in avoiding order reduction. A simple condition on the weight functions is shown to be both necessary and sufficient to avoid order reduction.

[1]  A. Iserles On the numerical quadrature of highly‐oscillating integrals I: Fourier transforms , 2004 .

[2]  Antonio G. García A Brief Walk Through Sampling Theory , 2002 .

[3]  E. Hairer,et al.  Numerical Energy Conservation for Multi-Frequency Oscillatory Differential Equations , 2005 .

[4]  Ernst Hairer,et al.  Simulating Hamiltonian dynamics , 2006, Math. Comput..

[5]  Jack Wisdom,et al.  The Origin of the Kirkwood Gaps: A Mapping for Asteroidal Motion Near the 3/1 Commensurability , 1982 .

[6]  W. Arendt Vector-valued laplace transforms and cauchy problems , 2002 .

[7]  Ernst Hairer,et al.  Solving Ordinary Differential Equations I: Nonstiff Problems , 2009 .

[8]  J. Verwer,et al.  Stability of Runge-Kutta Methods for Stiff Nonlinear Differential Equations , 1984 .

[9]  Marlis Hochbruck,et al.  Error analysis of exponential integrators for oscillatory second-order differential equations , 2006 .

[10]  E. Hairer,et al.  Geometric Numerical Integration , 2022, Oberwolfach Reports.

[11]  Mark E. Tuckerman,et al.  Reversible multiple time scale molecular dynamics , 1992 .

[12]  E. Hairer,et al.  Solving Ordinary Differential Equations II: Stiff and Differential-Algebraic Problems , 2010 .

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

[14]  F. Krogh,et al.  Solving Ordinary Differential Equations , 2019, Programming for Computations - Python.

[15]  E. Hairer,et al.  Solving Ordinary Differential Equations II , 2010 .

[16]  Klaus Schulten,et al.  Generalized Verlet Algorithm for Efficient Molecular Dynamics Simulations with Long-range Interactions , 1991 .

[17]  Ernst Hairer,et al.  Long-Time Energy Conservation of Numerical Methods for Oscillatory Differential Equations , 2000, SIAM J. Numer. Anal..

[18]  Robert D. Skeel,et al.  Long-Time-Step Methods for Oscillatory Differential Equations , 1998, SIAM J. Sci. Comput..