Related papers: Long-term planetary integration with individual ti…
Symplectic N-body integrators are widely used to study problems in celestial mechanics. The most popular algorithms are of 2nd and 4th order, requiring 2 and 6 substeps per timestep, respectively. The number of substeps increases rapidly…
Symplectic integration algorithms have become popular in recent years in long-term orbital integrations because these algorithms enforce certain conservation laws that are intrinsic to Hamiltonian systems. For problems with large variations…
Reliable studies of the long-term dynamics of planetary systems require numerical integrators that are accurate and fast. The challenge is often formidable because the chaotic nature of many systems requires relative numerical error bounds…
Modern N-body techniques for planetary dynamics are generally based on symplectic algorithms specially adapted to the Kepler problem. These methods have proven very useful in studying planet formation, but typically require the timestep for…
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…
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 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…
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…
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…
Computational efficiency demands discretised, hierarchically organised, and individually adaptive time-step sizes (known as the block-step scheme) for the time integration of N-body models. However, most existing N-body codes adapt…
Symplectic integrators are the tool of choice for many researchers studying dynamical systems because of their good long-term energy conservation properties. For systems with a dominant central mass, symplectic integrators are also highly…
$N$-body integrations are used to model a wide range of astrophysical dynamics, but they suffer from errors which make their orbits diverge exponentially in time from the correct orbits. Over long time-scales, their reliability needs to be…
Due to the chaotic nature of the Solar System, the question of its dynamic long-term stability can only be answered in a statistical sense, e.g. based on numerical ensemble integrations of nearby orbits. Destabilization, including…
Due to the chaotic nature of the Solar System, the question of its long-term stability can only be answered in a statistical sense, for instance, based on numerical ensemble integrations of nearby orbits. Destabilization of the inner…
I report results from accurate numerical integrations of Solar System orbits over the past 100Myr with the integrator package HNBody. The simulations used different integrator algorithms, step sizes, initial conditions, and included effects…
Orbital parameters of planets are fitted directly to an appropriate set of observations. It is shown how to use the rigorous Deming method combined with a numerical integration of gravitation equations. In all, 65 parameters of the nine…
The symplectic Wisdom-Holman map revolutionized long-term integrations of planetary systems. There is freedom in such methods of how to split the Hamiltonian and which coordinate system to employ, and several options have been proposed in…
We present an analytical proof assisted by computer calculations for the dynamical stability of the eight main planets and Pluto for the next 100,000 years. It means that the semi-major axes of the planets will not change significantly…
Symplectic integration algorithms are well-suited for long-term integrations of Hamiltonian systems because they preserve the geometric structure of the Hamiltonian flow. However, this desirable property is generally lost when adaptive…
We show that short-term perturbations among massive planets in multiple planet systems can result in radial velocity variations of the central star which differ substantially from velocity variations derived assuming the planets are…