Related papers: Fast and reliable symplectic integration for plane…
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…
Given that secular perturbations in a binary system not only excite high orbital eccentricities but also alter the planetary orbital inclination, the classical Keplerian orbital model is no longer applicable for orbital retrieval. The…
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…
We show that symplectic Runge-Kutta methods provide effective symplectic integrators for Hamiltonian systems with index one constraints. These include the Hamiltonian description of variational problems subject to position and velocity…
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…
In recent publications, the construction of explicit symplectic integrators for Schwarzschild and Kerr type spacetimes is based on splitting and composition methods for numerical integrations of Hamiltonians or time-transformed Hamiltonians…
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…
Optimization of low-thrust trajectories that involve a larger number of orbit revolutions is considered a challenging problem. This paper describes a high-precision symplectic method and optimization techniques to solve the minimum-energy…
Many exoplanetary systems are multiplanet configurations whose long-term dynamics are governed by N-body gravitational interactions. Consequently, their detection signatures cannot be adequately described by Keplerian orbits. Accurately…
Using a Newtonian model of the Solar System with all 8 planets, we perform extensive tests on various symplectic integrators of high orders, searching for the best splitting scheme for long term studies in the Solar System. These…
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…
We report our new code (named SACRA) for numerical relativity simulations in which an adaptive mesh refinement algorithm is implemented. In this code, the Einstein equations are solved in the BSSN formalism with a fourth-order finite…
Aims: We investigate the orbital and phase space properties of loss cone stars that interact strongly with a hard, high-redshift binary supermassive black hole (SMBH) system formed in a cosmological scenario. Methods: We use a novel hybrid…
Background: Ab initio many-body methods have been developed over the past ten years to address mid-mass nuclei... As progress in the design of inter-nucleon interactions is made, further efforts must be made to tailor many-body methods.…
Transit timing variations - deviations from strict periodicity between successive passages of a transiting planet - can be used to probe the structure and dynamics of multiple-planet systems. In this paper, we examine prospects for…
We present a family of modified Hermite integrators of arbitrary order possessing superior behaviour for the integration of Keplerian and near-Keplerian orbits. After recounting the derivation of Hermite N-body integrators of arbitrary…
The quantum many-body bound-state problem in its computationally successful coupled cluster method (CCM) representation is reconsidered. In conventional practice one factorizes the ground-state wave functions $|\Psi\rangle= e^S…
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…
A wide variety of outstanding problems in astrophysics involve the motion of a large number of particles ($N\gtrsim 10^{6}$) under the force of gravity. These include the global evolution of globular clusters, tidal disruptions of stars by…
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…