Related papers: Time stepping N-body simulations
Symplectic integrators are a foundation to the study of dynamical $N$-body phenomena, at scales ranging from from planetary to cosmological. These integrators preserve the Poincar\'e invariants of Hamiltonian dynamics. The $N$-body…
Recently, our group developed explicit symplectic methods for curved spacetimes that are not split into several explicitly integrable parts, but are via appropriate time transformations. Such time-transformed explicit symplectic integrators…
It is well known that symplectic integrators lose their near energy preservation properties when variable step sizes are used. The most common approach to combine adaptive step sizes and symplectic integrators involves the Poincar\'e…
We present new almost time-reversible integrators for solution of planetary systems consisting of "planets" and a dominant mass ("star"). The algorithms can be considered adaptive generalizations of the Wisdom--Holman method, in which all…
A second order explicit one-step numerical method for the initial value problem of the general ordinary differential equation is proposed. It is obtained by natural modifications of the well-known leapfrog method, which is a second order,…
This is the first in a series of papers devoted to fully general-relativistic $N$-body simulations applied to late-time cosmology. The purpose of this paper is to present the combination of a numerical relativity scheme, discretization…
The leapfrog integrator is widely used because of its excellent stability in molecular dynamics simulation. This is recognized as being due to the existence of a discrete variational structure of the equations. We introduce a modified…
A fixed time-step variational integrator cannot preserve momentum, energy, and symplectic form simultaneously for nonintegrable systems. This barrier can be overcome by treating time as a discrete dynamic variable and deriving adaptive…
A basic leapfrog integrator and its energy-preserving and variational / symplectic variants are proposed and studied for the numerical integration of the equations of motion of relativistic charged particles in an electromagnetic field. The…
We derive a fully computable aposteriori error estimator for a Galerkin finite element solution of the wave equation with explicit leapfrog time-stepping. Our discrete formulation accommodates both time evolving meshes and leapfrog based…
We present a multiscale integrator for Hamiltonian systems with slowly varying quadratic stiff potentials that uses coarse timesteps (analogous to what the impulse method uses for constant quadratic stiff potentials). This method is based…
We report an accessible and robust tool for evaluating the effects of Coulomb collisions on a test particle in a plasma that obeys Maxwell-J\"uttner statistics. The implementation is based on the Beliaev-Budker collision integral which…
A new algorithm for numerical integration of the rigid-body equations of motion is proposed. The algorithm uses the leapfrog scheme and the quantities involved are angular velocities and orientational variables which can be expressed in…
The existence of explicit symplectic integrators for general nonseparable Hamiltonian systems is an open and important problem in both numerical analysis and computing in science and engineering, as explicit integrators are usually more…
A drift-kick-drift (DKD) type leapfrog symplectic integrator applied for a time-transformed separable Hamiltonian (or time-transformed symplectic integrator; TSI) has been known to conserve the Kepler orbit exactly. We find that for an…
Most direct N-body integrations of planetary systems use a symplectic integrator with a fixed timestep. A large timestep is desirable in order to speed up the numerical simulations. However, simulations yield unphysical results if the…
Symplectic integrators have made it possible to study the long-term evolution of planetary systems with direct N-body simulations. In this paper we reassess the accuracy of such simulations by running a convergence test on 20Myr…
Astrophysical research in recent decades has made significant progress thanks to the availability of various $N$-body simulation techniques. With the rapid development of high-performance computing technologies, modern simulations have been…
We introduce Hermite-leapfrog methods for first order wave systems. The new Hermite-leapfrog methods pair leapfrog time-stepping with the Hermite methods of Goodrich and co-authors. The new schemes stagger field variables in both time and…
We suggest a numerical integration procedure for solving the equations of motion of certain classical spin systems which preserves the underlying symplectic structure of the phase space. Such symplectic integrators have been successfully…