Related papers: Long-term planetary integration with individual ti…
We study the long term orbital evolution of a terrestrial planet under the gravitational perturbations of a giant planet. In particular, we are interested in situations where the two planets are in the same plane and are relatively close.…
The possibility of measuring the post-Newtonian gravitoelectric correction to the orbital period of a test particle freely orbiting a spherically symmetric mass in the Solar System is analyzed. It should be possible, in principle, to detect…
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…
Recent observational missions have uncovered a significant number of compact multi-exoplanet systems. The tight orbital spacing of these systems has led to much effort being applied to the understanding of their stability; however, a key…
We present an explicit multiscale algorithm for solving differential equations for problems with high-frequency modes that can be averaged over by separating and scaling the fast and slow dynamics within a single equation. We introduce a…
The Kepler space mission has detected a large number of exoplanets in multiple transiting planet systems. Previous studies found that these Kepler multiple planet systems exhibit an intra-system uniformity, namely planets in the same system…
The time step criterion plays a crucial role in direct N-body codes. If not chosen carefully, it will cause a secular drift in the energy error. Shared, adaptive time step criteria commonly adopt the minimum pairwise time step, which…
Periodic solutions of the three body problem are very important for understanding its dynamics either in a theoretical framework or in various applications in celestial mechanics. In this paper we discuss the computation and continuation of…
Modern lunar-planetary ephemerides are numerically integrated on the observational timespan of more than 100 years (with the last 20 years having very precise astrometrical data). On such long timespans, not only finite difference…
If an orbit is fitted from combined RV and astrometric data, the orbit should be physically consistent with both data sets. The Keplerian orbit of a planet is a highly nonlinear function of seven parameters. The astrometric orbit problem…
We present an algorithm for the rapid numerical integration of a time-periodic ODE with a small dissipation term that is $C^1$ in the velocity. Such an ODE arises as a model of spin-orbit coupling in a star/planet system, and the motivation…
Planetary systems with more than two bodies will experience orbital crossings at a time related to the initial orbital separations of the planets. After a crossing, the system enters a period of chaotic evolution ending in the reshaping of…
We explore the orbital dynamics of systems consisting of three planets, each as massive as the Earth, on coplanar, initially circular, orbits about a star of one solar mass. The initial semimajor axes of the planets are equally spaced in…
In order to perform numerical studies of long-term stability in nonlinear Hamiltonian systems, one needs a numerical integration algorithm which is symplectic. Further, this algorithm should be fast and accurate. In this paper, we propose…
We demonstrate that the system of Jovian planets (Sun+Jupiter+Saturn+Uranus+Neptune), integrated for 200 million years as an isolated 5-body system using many sets of initial conditions all within the uncertainty bounds of their currently…
Direct N-body simulations and symplectic integrators are effective tools to study the long-term evolution of planetary systems. The Wisdom-Holman (WH) integrator in particular has been used extensively in planetary dynamics as it allows for…
A planetary instability occurring at time $<100$ My after formation of the giant planets in our solar system can be responsible for some characteristics of the inner solar system. However, the actual influence of the instability on the…
State-of-the-art long-term solar system integrations include several second order effects such as the Sun's quadrupole moment J2 and a contribution from asteroids (plus the Moon and general relativity). We recently showed that including 10…
A pure two-body problem has seven integrals including the Kepler energy, the Laplace vector, and the angular momentum vector. However, only five of them are independent. When the five independent integrals are preserved, the two other…
An optimization program is used to re-adjust the initial conditions, in order to reproduce as closely as possible the predictions of a complete ephemeris by using simplified equations of motion in the numerical integration. The adjustment…