Secular dynamics of hierarchical multiple systems composed of nested binaries, with an arbitrary number of bodies and arbitrary hierarchical structure. III. Suborbital effects: hybrid integration techniques and orbit-averaging corrections

Abstract

The SecularMultiple code, presented in two previous papers of this series, integrates the long-term dynamical evolution of multiple systems with any number of bodies and hierarchical structure, provided that the system is composed of nested binaries. In the formalism underlying SecularMultiple, we previously averaged over all orbits in the system. This approximation significantly speeds up numerical integration of the equations of motion, making large population synthesis studies possible. However, the orbit averaging approximation can break down when the secular evolution timescale of the system is comparable to or shorter than any of the orbital periods in the system. Here, we present an update to SecularMultiple in which we incorporate hybrid integration techniques, and orbit-averaging corrections. With this update, the user can specify which orbits should be integrated directly (without averaging), or assuming averaged orbits. For orbits that are integrated directly, we implemented two integration techniques, one which is based on the regularised Kustaanheimo-Stiefel equations of motion in element form. We also implemented analytical orbit-averaging corrections for pairwise interactions to quadrupole order. The updates presented here provide more flexibility for integrating the long-term dynamical evolution of hierarchical multiple systems. By effectively combining direct integration and orbit averaging the long-term evolution can be accurately computed, but with significantly lower computational cost compared to existing direct N-body codes. We give a number of examples in which the new features are beneficial. Our updated code, which is written in C++ supplemented with a user-friendly interface in Python, is freely available.

Accepted for publication in MNRAS. 15 pages, 4 figures. Code can be found at https://github.com/hamers/secularmultiple

References in corpus (18)

• Lidov-Kozai Cycles with Gravitational Radiation: Merging Black Holes in Isolated Triple Systems
• Mergers and Obliquities in Stellar Triples
• Implementing Few-Body Algorithmic Regularization with Post-Newtonian Terms
• Satellite dynamics on the Laplace surface
• Steady-state planet migration by the Kozai-Lidov mechanism in stellar binaries
• The Kozai-Lidov Mechanism in Hydrodynamical Disks
• Black hole mergers from quadruples
• The triple evolution dynamical instability: Stellar collisions in the field and the formation of exotic binaries
• Algorithmic regularization with velocity-dependent forces
• Black hole-neutron star mergers from triples
• Spin-Orbit Misalignment of Merging Black-Hole Binaries with Tertiary Companions
• Induced Ellipticity for Inspiraling Binary Systems
• Secular chaotic dynamics in hierarchical quadruple systems, with applications to hot Jupiters in stellar binaries and triples
• A slow-down time-transformed symplectic integrator for solving the few-body problem
• MSTAR -- a fast parallelised algorithmically regularised integrator with minimum spanning tree coordinates
• Chaotic quadruple secular evolution and the production of misaligned exomoons and Warm Jupiters in stellar multiples
• Misaligned accretion disc formation via Kozai-Lidov oscillations
• Massive binary star mergers in galactic nuclei: implications for blue stragglers, binary S-stars and gravitational waves

Cited by in corpus (6)

• Multiple Stellar Evolution: a population synthesis algorithm to model the stellar, binary, and dynamical evolution of multiple-star systems
• WD 1856 b: a close giant planet around a white dwarf that could have survived a common-envelope phase
• Was GW190412 born from a hierarchical 3+1 quadruple configuration?
• Properties of von Zeipel-Lidov-Kozai oscillations in triple systems at the quadrupole order: relaxing the test particle approximation
• Forming Hot Subluminous Stars from Hierarchical Triples -- I. The Role of an Outer Tertiary on Formation Channels
• Laplace surface dynamics, revisited: satellites, exo-planets and debris with distant, eccentric companions