Related papers: Significant improvement in planetary system simula…
Symplectic integrators separate a problem into parts that can be solved in isolation, alternately advancing these sub-problems to approximate the evolution of the complete system. Problems with a single, dominant mass can use mixed-variable…
This paper is a review of the dynamics of a system of planets. It includes the study of averaged equations in both non-resonant and resonant systems and shows the great deal of situations in which the angle between the two semi-major axes…
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…
A symplectic integrator algorithm suitable for hierarchical triple systems is formulated and tested. The positions of the stars are followed in hierarchical Jacobi coordinates, whilst the planets are referenced purely to their primary. The…
Symplectic integrators are widely used for long-term integration of conservative astrophysical problems due to their ability to preserve the constants of motion; however, they cannot in general be applied in the presence of nonconservative…
Symplectic integrators are the preferred method of solving conservative $N$-body problems in cosmological, stellar cluster, and planetary system simulations because of their superior error properties and ability to compute orbital…
Symplectic numerical methods have become a widely-used choice for the accurate simulation of Hamiltonian systems in various fields, including celestial mechanics, molecular dynamics and robotics. Even though their characteristics are…
We derive, and discuss the properties of, a symplectic map for the dynamics of bodies on nearly parabolic orbits. The orbits are perturbed by a planet on a circular, coplanar orbit interior to the pericenter of the parabolic orbit. The map…
We present new splitting methods designed for the numerical integration of near-integrable Hamiltonian systems, and in particular for planetary N-body problems, when one is interested in very accurate results over a large time span. We…
(abbreviated) We use a semi-numerical approach to study the secular behavior of a system composed of a central star and two massive planets in eccentric co-planar orbits. We show that the secular dynamics of this system can be described…
The interval approach to computation of dynamics of celestial bodies in the planetary problem has been considered. It is based on the refusal from idealization of infinitely high resolving capacity of measuring tools, and forms an…
Most physical systems are modelled by an ordinary or a partial differential equation, like the n-body problem in celestial mechanics. In some cases, for example when studying the long term behaviour of the solar system or for complex…
We use numerical $N$-body experiments to explore the statistics of multiple systems formed in small-$N$ subclusters, i.e. the distributions of orbital semi-major axis, $a$, orbital eccentricity, $e$, mass ratio, $q$, mutual orbital…
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…
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…
Explicit symplectic integrators have been important tools for accurate and efficient approximations of mechanical systems with separable Hamiltonians. For the first time, the article proposes for arbitrary Hamiltonians similar integrators,…
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…
The long-term variations in the orbit of the Earth govern the insolation on its surface and hence its climate. The use of the astronomical signal, whose imprint has been recovered in the geological records, has revolutionized the…
Hamiltonian systems are differential equations which describe systems in classical mechanics, plasma physics, and sampling problems. They exhibit many structural properties, such as a lack of attractors and the presence of conservation…
This article reviews some integrators particularly suitable for the numerical resolution of differential equations on a large time interval. Symplectic integrators are presented. Their stability on exponentially large time is shown through…