Related papers: N-body integrators for planets in binary star syst…
The gravitational $N$-body problem, which is fundamentally important in astrophysics to predict the motion of $N$ celestial bodies under the mutual gravity of each other, is usually solved numerically because there is no known general…
Many applications in computational physics that use numerical integrators based on splitting and composition can benefit from the development of optimized algorithms and from choosing the best ordering of terms. The cost in programming and…
Calculating the long term solution of ordinary differential equations, such as those of the $N$-body problem, is central to understanding a wide range of dynamics in astrophysics, from galaxy formation to planetary chaos. Because generally…
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…
The majority of star formation results in binaries or higher multiple systems, and planets in such systems are constrained to a limited range of orbital parameters in order to remain stable against perturbations from stellar companions.…
We present a new symplectic integrator designed for collisional gravitational $N$-body problems which makes use of Kepler solvers. The integrator is also reversible and conserves 9 integrals of motion of the $N$-body problem to machine…
We present a detailed comparison of several integration schemes applied to the dynamic system consisting of a charged particle on the Kerr background endowed with the axisymmetric electromagnetic test field. In particular, we compare the…
Many exoplanets are discovered in binary star systems in internal or in circumbinary orbits. Whether the planet can be habitable or not depends on the possibility to maintain liquid water on its surface, and therefore on the luminosity of…
In this work we propose a new numerical approach to distinguish between regular and chaotic orbits in Hamiltonian systems, based on the simultaneous integration of both the orbit and the deviation vectors using a symplectic scheme, hereby…
Symplectic methods, in particular the Wisdom-Holman map, have revolutionized our ability to model the long-term, conservative dynamics of planetary systems. However, many astrophysically important effects are dissipative. The consequences…
Leapfrog integration has been the method of choice in N-body simulations owing to its low computational cost for a symplectic integrator with second order accuracy. We introduce a new leapfrog integrator that allows for variable timesteps…
The dynamic equation of mass point in rotating coordinates is governed by Coriolis and centrifugal force, besides a corotating potential relative to frame. Such a system is no longer a canonical Hamiltonian system so that the construction…
Symplectic integrators that preserve the geometric structure of Hamiltonian flows and do not exhibit secular growth in energy errors are suitable for the long-term integration of N-body Hamiltonian systems in the solar system. However, the…
Dependable numerical results from long-time simulations require stable numerical integration schemes. For Hamiltonian systems, this is achieved with symplectic integrators, which conserve the symplectic condition and exactly solve for the…
We present a simple algorithm to switch between $N$-body time integrators in a reversible way. We apply it to planetary systems undergoing arbitrarily close encounters and highly eccentric orbits, but the potential applications are broader.…
Recently a new class of numerical integration methods -- ``mixed variable symplectic integrators'' -- has been introduced for studying long-term evolution in the conservative gravitational few-body problem. These integrators are an order of…
Hill's equations are an approximation that is useful in a number of areas of astrophysics including planetary rings and planetesimal disks. We derive a symplectic method for integrating Hill's equations based on a generalized leapfrog. This…
We present a simple choice of integration variables that can be used to exploit the near-integrable character of problems in celestial mechanics. The approach is based on the well-known principle of variation of parameters: instead of…
The formation and evolution of protoplanetary systems, the breeding grounds of planet formation, is a complex dynamical problem that involves many orders of magnitudes. To serve this purpose, we present a new hybrid algorithm that combines…
This paper describes a fourth-order integration algorithm for the gravitational N-body problem based on discrete Lagrangian mechanics. When used with shared timesteps, the algorithm is momentum conserving and symplectic. We generalize the…