-
Reversible time-step adaptation for the integration of few-body systems
Authors:
Tjarda C. N. Boekholt,
Timothee Vaillant,
Alexandre C. M. Correia
Abstract:
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 suffers from discontinuities in the time evolution of the time step. This has a large impact on the functioning of time step symmetrisation algorithms. We provide new…
▽ More
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 suffers from discontinuities in the time evolution of the time step. This has a large impact on the functioning of time step symmetrisation algorithms. We provide new demonstrations of previous findings that a smooth and weighted average over all pairwise time steps in the N-body system, improves the level of energy conservation. Furthermore, we compare the performance of 27 different time step criteria, by considering 3 methods for weighting time steps and 9 symmetrisation methods. We present performance tests for strongly chaotic few-body systems, including unstable triples, giant planets in a resonant chain, and the current Solar System. We find that the harmonic symmetrisation methods (methods A3 and B3 in our notation) are the most robust, in the sense that the symmetrised time step remains close to the time step function. Furthermore, based on our Solar System experiment, we find that our new weighting method based on direct pairwise averaging (method W2 in our notation), is slightly preferred over the other methods.
△ Less
Submitted 19 December, 2022;
originally announced December 2022.
-
A direct N-body integrator for modelling the chaotic, tidal dynamics of multi-body extrasolar systems: TIDYMESS
Authors:
Tjarda C. N. Boekholt,
Alexandre C. M. Correia
Abstract:
Tidal dissipation plays an important role in the dynamical evolution of moons, planets, stars and compact remnants. The interesting complexity originates from the interplay between the internal structure and external tidal forcing. Recent and upcoming observing missions of exoplanets and stars in the Galaxy help to provide constraints on the physics of tidal dissipation. It is timely to develop ne…
▽ More
Tidal dissipation plays an important role in the dynamical evolution of moons, planets, stars and compact remnants. The interesting complexity originates from the interplay between the internal structure and external tidal forcing. Recent and upcoming observing missions of exoplanets and stars in the Galaxy help to provide constraints on the physics of tidal dissipation. It is timely to develop new N-body codes, which allow for experimentation with various tidal models and numerical implementations. We present the open-source N-body code TIDYMESS, which stands for ``TIdal DYnamics of Multi-body ExtraSolar Systems''. This code implements a creep deformation law for the bodies, parametrized by their fluid Love numbers and fluid relaxation times. Due to tidal and centrifugal deformations, we approximate the general shape of a body to be an ellipsoid. We calculate the associated gravitational field to quadruple order, from which we derive the gravitational accelerations and torques. The equations of motion for the orbits, spins and deformations are integrated directly using a fourth-order integration method based on a symplectic composition. We implement a novel integration method for the deformations, which allows for a time step solely dependent on the orbits, and not on the spin periods or fluid relaxation times. This feature greatly speeds up the calculations, while also improving the consistency when comparing different tidal regimes. We demonstrate the capabilities and performance of TIDYMESS, particularly in the niche regime of parameter space where orbits are chaotic and tides become non-linear.
△ Less
Submitted 8 September, 2022;
originally announced September 2022.
-
The paradox of infinitesimal granularity: Chaos and the reversibility of time in Newton's theory of gravity
Authors:
Simon Portegies Zwart,
Tjarda Boekholt
Abstract:
The fundamental laws of physics are time-symmetric, but our macroscopic experience contradicts this. The time reversibility paradox is partly a consequence of the unpredictability of Newton's equations of motion. We measure the dependence of the fraction of irreversible, gravitational N-body systems on numerical precision and find that it scales as a power law. The stochastic wave packet reduction…
▽ More
The fundamental laws of physics are time-symmetric, but our macroscopic experience contradicts this. The time reversibility paradox is partly a consequence of the unpredictability of Newton's equations of motion. We measure the dependence of the fraction of irreversible, gravitational N-body systems on numerical precision and find that it scales as a power law. The stochastic wave packet reduction postulate then introduces fundamental uncertainties in the Cartesian phase space coordinates that propagate through classical three-body dynamics to macroscopic scales within the triple's lifetime. The spontaneous collapse of the wave function then drives the global chaotic behavior of the Universe through the superposition of triple systems (and probably multi-body systems). The paradox of infinitesimal granularity then arises from the superposition principle, which states that any multi-body system is composed of an ensemble of three-body problems.
△ Less
Submitted 6 September, 2022;
originally announced September 2022.
-
On the Jacobi capture origin of binaries with applications to the Earth-Moon system and black holes in galactic nuclei
Authors:
Tjarda C. N. Boekholt,
Connar Rowan,
Bence Kocsis
Abstract:
Close encounters between two bodies in a disc often result in a single orbital deflection. However, within their Jacobi volumes, where the gravitational forces between the two bodies and the central body become competitive, temporary captures with multiple close encounters become possible outcomes: a Jacobi capture. We perform 3-body simulations in order to characterise the dynamics of Jacobi capt…
▽ More
Close encounters between two bodies in a disc often result in a single orbital deflection. However, within their Jacobi volumes, where the gravitational forces between the two bodies and the central body become competitive, temporary captures with multiple close encounters become possible outcomes: a Jacobi capture. We perform 3-body simulations in order to characterise the dynamics of Jacobi captures in the plane. We find that the phase space structure resembles a Cantor-like set with a fractal dimension of about 0.4. The lifetime distribution decreases exponentially, while the distribution of the closest separation follows a power law with index 0.5. In our first application, we consider the Jacobi capture of the Moon. We demonstrate that both tidal captures and giant impacts are possible outcomes. The impact speed is well approximated by a parabolic encounter, while the impact angles follow that of a uniform beam on a circular target. Jacobi captures at larger heliocentric distances are more likely to result in tidal captures. In our second application, we find that Jacobi captures with gravitational wave dissipation can result in the formation of binary black holes in galactic nuclei. The eccentricity distribution is approximately super-thermal and includes both prograde and retrograde orientations. We conclude that dissipative Jacobi captures form an efficient channel for binary formation, which motivates further research into establishing the universality of Jacobi captures across multiple astrophysical scales.
△ Less
Submitted 24 November, 2022; v1 submitted 17 March, 2022;
originally announced March 2022.
-
Chaos in self-gravitating many-body systems: Lyapunov time dependence of $N$ and the influence of general relativity
Authors:
Simon F. Portegies Zwart,
Tjarda C. N. Boekholt,
Emiel Por,
Adrian S. Hamers,
Steve L. W. McMillan
Abstract:
In self-gravitating $N$-body systems, small perturbations introduced at the start, or infinitesimal errors that are produced by the numerical integrator or are due to limited precision in the computer, grow exponentially with time. For Newton's gravity, we confirm earlier results that for relatively homogeneous systems, this rate of growth per crossing time increases with $N$ up to $N \sim 30$, bu…
▽ More
In self-gravitating $N$-body systems, small perturbations introduced at the start, or infinitesimal errors that are produced by the numerical integrator or are due to limited precision in the computer, grow exponentially with time. For Newton's gravity, we confirm earlier results that for relatively homogeneous systems, this rate of growth per crossing time increases with $N$ up to $N \sim 30$, but that for larger systems, the growth rate has a weaker scaling with $N$. For concentrated systems, however, the rate of exponential growth continues to scale with $N$. In relativistic self-gravitating systems, the rate of growth is almost independent of $N$. This effect, however, is only noticeable when the system's mean velocity approaches the speed of light to within three orders of magnitude. The chaotic behavior of systems with more than a dozen bodies for the usually adopted approximation of only solving the pairwise interactions in the Einstein-Infeld-Hoffmann equation of motion is qualitatively different than when the interaction terms (or cross terms) are taken into account. This result provides a strong motivation for follow-up studies on the microscopic effect of general relativity on orbital chaos, and on the influence of higher-order cross-terms in the Taylor-series expansion of the Einstein-Infeld-Hoffmann equations of motion.
△ Less
Submitted 17 November, 2021; v1 submitted 22 September, 2021;
originally announced September 2021.
-
Radio astronomy and Space science in Azores: enhancing the Atlantic VLBI infrastructure cluster
Authors:
Domingos Barbosa,
Bruno Coelho,
Sonia Antón,
Miguel Bergano,
Tjarda Boekholt,
Alexandre C. M. Correia,
Dalmiro Maia,
João Pandeirada,
Valério Ribeiro,
Jason Adams,
João Paulo Barraca,
Diogo Gomes,
Bruno Morgado
Abstract:
Radio astronomy and Space Infrastructures in the Azores have a great scientific and industrial interest because they benefit from a unique geographical location in the middle of the North Atlantic allowing a vast improvement in the sky coverage. This fact obviously has a very high added value for: i) the establishment of space tracking and communications networks for the emergent global small sate…
▽ More
Radio astronomy and Space Infrastructures in the Azores have a great scientific and industrial interest because they benefit from a unique geographical location in the middle of the North Atlantic allowing a vast improvement in the sky coverage. This fact obviously has a very high added value for: i) the establishment of space tracking and communications networks for the emergent global small satellite fleets ii) it is invaluable to connect the radio astronomy infrastructure networks in Africa, Europe and America continents using Very Large Baseline Interferometry (VLBI) techniques, iii) it allows excellent potential for monitoring space debris and Near Earth Objects (NEOs). There is in S. Miguel island a 32-metre SATCOM antenna that could be integrated in advanced VLBI networks and be capable of additional Deep Space Network ground support. This paper explores the space science opportunities offered by the upgrade of the S. Miguel 32-metre SATCOM antenna into a world-class infrastructure for radio astronomy and space exploration: it would enable a Deep Space Network mode and would constitute a key space facility for data production, promoting local digital infrastructure investments and the testing of cutting-edge information technologies. Its Atlantic location also enables improvements in angular resolution, provides many baseline in East-West and North-South directions connecting the emergent VLBI stations in America to Europe and Africa VLBI arrays therefore contributing for greater array imaging capabilities especially for sources or well studied fields close to or below the celestial equator, where ESO facilities, ALMA, SKA and its precursors do or will operate and observe in the coming decades.
△ Less
Submitted 7 July, 2021;
originally announced July 2021.
-
Gargantuan chaotic gravitational three-body systems and their irreversibility to the Planck length
Authors:
Tjarda Boekholt,
Simon Portegies Zwart,
Mauri Valtonen
Abstract:
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 consequence is that integrations of the N-body problem unavoidably magnify truncation and rounding errors to macroscopic scales. Hitherto, a quantitative relation betwee…
▽ More
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 consequence is that integrations of the N-body problem unavoidably magnify truncation and rounding errors to macroscopic scales. Hitherto, a quantitative relation between chaos in stellar dynamical systems and the level of irreversibility remained undetermined. In this work we study chaotic three-body systems in free fall initially using the accurate and precise N-body code Brutus, which goes beyond standard double-precision arithmetic. We demonstrate that the fraction of irreversible solutions decreases as a power law with numerical accuracy. This can be derived from the distribution of amplification factors of small initial perturbations. Applying this result to systems consisting of three massive black holes with zero total angular momentum, we conclude that up to five percent of such triples would require an accuracy of smaller than the Planck length in order to produce a time-reversible solution, thus rendering them fundamentally unpredictable.
△ Less
Submitted 10 February, 2020;
originally announced February 2020.
-
Newton vs the machine: solving the chaotic three-body problem using deep neural networks
Authors:
Philip G. Breen,
Christopher N. Foley,
Tjarda Boekholt,
Simon Portegies Zwart
Abstract:
Since its formulation by Sir Isaac Newton, the problem of solving the equations of motion for three bodies under their own gravitational force has remained practically unsolved. Currently, the solution for a given initialization can only be found by performing laborious iterative calculations that have unpredictable and potentially infinite computational cost, due to the system's chaotic nature. W…
▽ More
Since its formulation by Sir Isaac Newton, the problem of solving the equations of motion for three bodies under their own gravitational force has remained practically unsolved. Currently, the solution for a given initialization can only be found by performing laborious iterative calculations that have unpredictable and potentially infinite computational cost, due to the system's chaotic nature. We show that an ensemble of solutions obtained using an arbitrarily precise numerical integrator can be used to train a deep artificial neural network (ANN) that, over a bounded time interval, provides accurate solutions at fixed computational cost and up to 100 million times faster than a state-of-the-art solver. Our results provide evidence that, for computationally challenging regions of phase-space, a trained ANN can replace existing numerical solvers, enabling fast and scalable simulations of many-body systems to shed light on outstanding phenomena such as the formation of black-hole binary systems or the origin of the core collapse in dense star clusters.
△ Less
Submitted 16 October, 2019;
originally announced October 2019.
-
A Keplerian-based Hamiltonian Splitting for Gravitational $N$-body Simulations
Authors:
G. Gonçalves Ferrari,
T. Boekholt,
S. F. Portegies Zwart
Abstract:
We developed a Keplerian-based Hamiltonian splitting for solving the gravitational $N$-body problem. This splitting allows us to approximate the solution of a general $N$-body problem by a composition of multiple, independently evolved $2$-body problems. While the Hamiltonian splitting is exact, we show that the composition of independent $2$-body problems results in a non-symplectic non-time-symm…
▽ More
We developed a Keplerian-based Hamiltonian splitting for solving the gravitational $N$-body problem. This splitting allows us to approximate the solution of a general $N$-body problem by a composition of multiple, independently evolved $2$-body problems. While the Hamiltonian splitting is exact, we show that the composition of independent $2$-body problems results in a non-symplectic non-time-symmetric first-order map. A time-symmetric second-order map is then constructed by composing this basic first-order map with its self-adjoint. The resulting method is precise for each individual $2$-body solution and produces quick and accurate results for near-Keplerian $N$-body systems, like planetary systems or a cluster of stars that orbit a supermassive black hole. The method is also suitable for integration of $N$-body systems with intrinsic hierarchies, like a star cluster with primordial binaries. The superposition of Kepler solutions for each pair of particles makes the method excellently suited for parallel computing; we achieve $\gtrsim 64\%$ efficiency for only $8$ particles per core, but close to perfect scaling for $16384$ particles on a $128$ core distributed-memory computer. We present several implementations in \texttt{Sakura}, one of which is publicly available via the AMUSE framework.
△ Less
Submitted 13 February, 2014;
originally announced February 2014.