Related papers: Are long-term $N$-body simulations reliable?
The inverse problem is studied in multi-body systems with nonlinear dynamics representing, e.g., phase-locked wave systems, standard multimode and random lasers. Using a general model for four-body interacting complex-valued variables we…
The orbital solutions of published multi-planet systems are not necessarily dynamically stable on timescales comparable to the lifetime of the system as a whole. For this reason, dynamical tests of the architectures of proposed exoplanetary…
Spin-orbit coupling of planetary systems plays an important role in the dynamics and habitability of planets. However, symplectic integrators that can accurately simulate not only how orbit affects spin but also how spin affects orbit have…
Simulating the evolution of the gravitational N-body problem becomes extremely computationally expensive as N increases since the problem complexity scales quadratically with the number of bodies. We study the use of Artificial Neural…
In the preliminary design of space missions it can be useful to identify regions of dynamics that drive the system's behaviour or separate qualitatively different dynamics. The Lagrangian Coherent Structure (LCS) has been widely used in the…
This work discusses the main analogies and differences between the deterministic approach underlying most cosmological N-body simulations and the probabilistic interpretation of the problem that is often considered in mathematics and…
We present implementations of a fourth-order symplectic integrator on graphic processing units for three $N$-body models with long-range interactions of general interest: the Hamiltonian Mean Field, Ring and two-dimensional self-gravitating…
A revision of Stodolkiewicz's Monte-Carlo code is used to simulate evolution of star clusters. The new method treats each superstar as a single star and follows the evolution and motion of all individual stellar objects. The first…
In this work we propose a new numerical approach to distinguish between regular and chaotic orbits in Hamiltonian systems, based on the simultaneous integration of both the orbit and the deviation vectors using a symplectic scheme, hereby…
In this fourth article on weighing the Galactic disk using the shape of the phase-space spiral, we have tested our method on a billion particle three-dimensional N-body simulation, comprised of a Milky Way like host galaxy and a merging…
An accurate and efficient method dealing with the few-body dynamics is important for simulating collisional N-body systems like star clusters and to follow the formation and evolution of compact binaries. We describe such a method which…
Binary and multiple star systems are extreme environments for the formation and long-term presence of extrasolar planets. Circumstellar planets are subject to gravitational perturbations from the distant companion star, and this interaction…
We present a qualitative diagnostic based on the continuous wavelet transform, for the detection of irregular behavior in time series of particle simulations. We apply the method to three qualitatively different gravitational 3-body…
The "gravitational million-body problem," to model the dynamical evolution of a self-gravitating, collisional N-body system with ~10^6 particles over many relaxation times, remains a major challenge in computational astrophysics.…
Many exo-solar systems discovered in the last decade consist of planets orbiting in resonant configurations and consequently, their evolution should show long-term stability. However, due to the mutual planetary interactions a multi-planet…
In this paper I focus on three topics related to the dynamical evolution of small galaxy groups, for which the input of N-body simulations has been decisive. These are the merging rates in compact groups, the properties of remnants of…
We demonstrate that the output of a cosmological N-body simulation can, to remarkable accuracy, be scaled to represent the growth of large-scale structure in a cosmology with parameters similar to but different from those originally…
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…
Upcoming Large Scale Structure surveys aim to achieve an unprecedented level of precision in measuring galaxy clustering. However, accurately modeling these statistics may require theoretical templates that go beyond second-order…
(Abridged) The probability of detecting additional companions is essentially unchanged with respect to the single-planet configurations, but after fitting and subtraction of orbits with astrometric signal-to-noise ratio $\alpha/\sigma_d\to…