Related papers: Integration of Few Body Celestial Systems Implemen…
Small bodies are time-capsules of different eras of solar system history from the most primitive materials within the solar system to evolved pieces of larger bodies. A small body sample return program is an essential component of small…
We describe the astrophysical and numerical basis of N-body simulations, both of collisional stellar systems (dense star clusters and galactic centres) and collisionless stellar dynamics (galaxies and large-scale structure). We explain and…
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…
Much of standard galaxy dynamics rests on the implicit assumption that the corresponding N-body problem is (near) integrable. This notion although leading to great simplification is by no means a fact. It is therefore important to develop…
Various solutions are displayed and analyzed (both analytically and numerically) of arecently-introduced many-body problem in the plane which includes both integrable and nonintegrable cases (depending on the values of the coupling…
We introduce here our new approach to modeling particle cloud evolution off surface of small bodies (asteroids and comets), following the evolution of ejected particles requires dealing with various time and spatial scales, in an efficient,…
It is widely known that numerically integrated orbits are more precise than analytical theories for celestial bodies. However, calculation of the positions of celestial bodies via numerical integration at time $t$ requires the amount of…
(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 two-body problem in General Relativity has been the subject of many analytical investigations. After reviewing some of the methods used to tackle this problem (and, more generally, the N-body problem), we focus on a new, recently…
The first integrals of the Kepler problem are used to compute preliminary orbits starting from two short observed arcs of a celestial body, which may be obtained either by optical or radar observations. We write polynomial equations for…
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…
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…
Modeling interacting galaxies to reproduce observed systems is still a challenge due to the extended parameter space (among other problems). Orbit and basic galaxy parameters can be tackled by fast simulation techniques like the restricted…
Chaos is present in most stellar dynamical systems and manifests itself through the exponential growth of small perturbations. Exponential divergence drives time irreversibility and increases the entropy in the system. A numerical…
We present a new time-stepping criterion for N-body simulations that is based on the true dynamical time of a particle. This allows us to follow the orbits of particles correctly in all environments since it has better adaptivity than…
We present a fully covariant transport framework for Molecular Dynamics that enables a consistent description of the evolution of relativistic N-body systems. For the first time, we derive relativistic equations of motion incorporating both…
We develop a formalism for General Relativistic N-body simulations in the weak field regime, suitable for cosmological applications. The problem is kept tractable by retaining the metric perturbations to first order, the first derivatives…
Given the high-precision modern space mission, a precise relativistic modeling of observations is required. By solving the eikonal equation with the post-Newtonian approximation, the light propagation is determined by the iterative method…
Conservative symmetric second-order one-step schemes are derived for dynamical systems describing various many-body systems using the Discrete Multiplier Method. This includes conservative schemes for the $n$-species Lotka-Volterra system,…
In this work, we present the hitherto most efficient and accurate method for the numerical integration of post-Newtonian equations of motion. We first transform the Poisson system as given by the post-Newtonian approximation to canonically…