Related papers: N-body integrators for planets in binary star syst…
Context. Due to our increasing knowledge on the Galactic and stellar neighborhood of the Solar System, modern long-period comet motion studies have to take into account both stellar perturbations and the overall Galactic potential. Aims.…
In the standard model of terrestrial planet formation, planets are formed through giant impacts of planetary embryos after the dispersal of the protoplanetary gas disc. Traditionally, $N$-body simulations have been used to investigate this…
The integration of the equations of motion in gravitational dynamical systems -- either in our Solar System or for extra-solar planetary system -- being non integrable in the global case, is usually performed by means of numerical…
Terrestrial planet formation theory is at a bottleneck, with the growing realization that pairwise collisions are treated far too simply. Here, and in our companion paper (Cambioni et al. 2019) that introduces the training methodology, we…
We discuss the performance of direct summation codes used in the simulation of astrophysical stellar systems on highly distributed architectures. These codes compute the gravitational interaction among stars in an exact way and have an…
We present WHFast, a fast and accurate implementation of a Wisdom-Holman symplectic integrator for long-term orbit integrations of planetary systems. WHFast is significantly faster and conserves energy better than all other Wisdom-Holman…
We investigate the range of applicability of a model for the real-space power spectrum based on N-body dynamics and a (quadratic) Lagrangian bias expansion. This combination uses the highly accurate particle displacements that can be…
We construct explicit integrators of arbitrary even orders of accuracy for massive point vortex dynamics in binary mixture of Bose--Einstein condensates proposed by Richaud et al. The integrators are symplectic and preserve the angular…
We develop a fourth-order numerical integrator to simulate the coupled spin and orbital motions of two rigid bodies having arbitrary mass distributions under the influence of their mutual gravitational potential. We simulate the dynamics of…
We give a short and elementary introduction to Lie group methods. A selection of applications of Lie group integrators are discussed. Finally, a family of symplectic integrators on cotangent bundles of Lie groups is presented and the notion…
Two families of symplectic methods specially designed for second-order time-dependent linear systems are presented. Both are obtained from the Magnus expansion of the corresponding first-order equation, but otherwise they differ in…
A survey of currently known planet-hosting stars indicates that approximately 25% of extrasolar planetary systems are within dual-star environments. Several of these systems contain stellar companions on moderately close orbits, implying…
We use a hybrid N-body program to study the evolution of massive black hole binaries in the centers of galaxies, mainly to understand the factors affecting the binary eccentricity, the response of the galaxy to the binary merger, and the…
This paper explores the stability of an Earth-like planet orbiting a solar mass star in the presence of an outer-lying intermediate mass companion. The overall goal is to estimate the fraction of binary systems that allow Earth-like planets…
The modern study of the dynamics of stellar systems requires the use of high-performance computers. Indeed, an accurate modelization of the structure and evolution of self-gravitating systems like planetary systems, open clusters, globular…
We present a new, simple, fast algorithm to numerically evolve disks of inelastically colliding particles surrounding a central star. Our algorithm adds negligible computational cost to the fastest existing collisionless N-body codes, and…
(abridged) This thesis presents theoretical and statistical techniques broadly related to systems of dynamically-interacting particles composed of several different types of populations. They are applied to observations of dense star…
Probabilistic solvers provide a flexible and efficient framework for simulation, uncertainty quantification, and inference in dynamical systems. However, like standard solvers, they suffer performance penalties for certain stiff systems,…
In this work, the benefits of the phase fitting technique are embedded in high order discrete Lagrangian integrators. The proposed methodology creates integrators with zero phase lag in a test Lagrangian in a similar way used in phase…
Multiple stellar systems are common especially among O and B stars. In order to accurately describe their dynamics, interactions among components must be accounted for. In this work, we describe the new dynamical model in Phoebe, which…