Source author record

Tjarda C. N. Boekholt

Tjarda C. N. Boekholt appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

astro-ph.EP physics.comp-ph astro-ph.CO astro-ph.HE astro-ph.IM gr-qc nlin.CD

Catalog footprint

What is connected

3works

7topics

4close collaborators

Actions

Connect this record

Open graph Browse works

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

preprint2022arXiv

Reversible time-step adaptation for the integration of few-body systems

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.

preprint2021arXiv

Chaos in self-gravitating many-body systems: Lyapunov time dependence of $N$ and the influence of general relativity

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.

preprint2016arXiv

The Origin of Chaos in the Orbit of Comet 1P/Halley

According to Munoz-Gutierrez et al. (2015) the orbit of comet 1P/Halley is chaotic with a surprisingly small Lyapunov time scale of order its orbital period. In this work we analyse the origin of chaos in Halley's orbit and the growth of perturbations, in order to get a better understanding of this unusually short time scale. We perform N-body simulations to model Halley's orbit in the Solar System and measure the separation between neighbouring trajectories. To be able to interpret the numerical results, we use a semi-analytical map to demonstrate different growth modes, i.e. linear, oscillatory or exponential, and transitions between these modes. We find the Lyapunov time scale of Halley's orbit to be of order 300 years, which is significantly longer than previous estimates in the literature. This discrepancy could be due to the different methods used to measure the Lyapunov time scale. A surprising result is that next to Jupiter, also encounters with Venus contribute to the exponential growth in the next 3000 years. Finally, we note an interesting application of the sub-linear, oscillatory growth mode to an ensemble of bodies moving through the Solar System. Whereas in the absence of encounters with a third body the ensemble spreads out linearly in time, the accumulation of weak encounters can increase the lifetime of such systems due to the oscillatory behaviour.