SpaceHub: A high-performance gravity integration toolkit for few-body problems in astrophysics

Abstract

We present the open source few-body gravity integration toolkit {\tt SpaceHub}. {\tt SpaceHub} offers a variety of algorithmic methods, including the unique algorithms AR-Radau, AR-Sym6, AR-ABITS and AR-chain${}^{+}$ which we show out-perform other methods in the literature and allow for fast, precise and accurate computations to deal with few-body problems ranging from interacting black holes to planetary dynamics. We show that AR-Sym6 and AR-chain${}^{+}$, with algorithmic regularization, chain algorithm, active round-off error compensation and a symplectic kernel implementation, are the fastest and most accurate algorithms to treat black hole dynamics with extreme mass ratios, extreme eccentricities and very close encounters. AR-Radau, the first regularized Radau integrator with round off error control down to 64 bits floating point machine precision, has the ability to handle extremely eccentric orbits and close approaches in long-term integrations. AR-ABITS, a bit efficient arbitrary precision method, achieves any precision with the least CPU cost compared to other open source arbitrary precision few-body codes. With the implementation of deep numerical and code optimization, these new algorithms in {\tt SpaceHub} prove superior to other popular high precision few-body codes in terms of performance, accuracy and speed.

arXiv commentSubmitted to MNRAS. Comments are welcome

References in corpus (13)

• IAS15: A fast, adaptive, high-order integrator for gravitational dynamics, accurate to machine precision over a billion orbits
• Implementing Few-Body Algorithmic Regularization with Post-Newtonian Terms
• A class of symplectic integrators with adaptive timestep for separable Hamiltonian systems
• Black Hole and Neutron Star Binary Mergers in Triple Systems: II. Merger Eccentricity and Spin-Orbit Misalignment
• Binary Mergers near a Supermassive Black Hole: Relativistic Effects in Triples
• Constraining the Black Hole Initial Mass Function with LIGO/VIRGO Observations
• MSTAR -- a fast parallelised algorithmically regularised integrator with minimum spanning tree coordinates
• Planetary Architectures in Interacting Stellar Environments
• Hot Jupiter and ultra-cold Saturn formation in dense star clusters
• The fate of close encounters between binary stars and binary supermassive black holes
• Giant planet swaps during close stellar encounters
• Hypervelocity binaries from close encounters with a SMBH-IMBH binary: orbital properties and diagnostics
• The cosmological distribution of compact object mergers from dynamical interactions with SMBH binaries

Cited by in corpus (9)

• Symmetry Breaking in Dynamical Encounters in the Disks of Active Galactic Nuclei
• Close Encounters of Stars with Stellar-mass Black Hole Binaries
• The two phases of core formation -- orbital evolution in the centres of ellipticals with supermassive black hole binaries
• BIFROST: simulating compact subsystems in star clusters using a hierarchical fourth-order forward symplectic integrator code
• Hot Jupiter formation in dense clusters: secular chaos in multi-planetary systems
• Examining the effects of dark matter spikes on eccentric intermediate mass ratio inspirals using N-body simulations
• Evolution of Massive Black Hole Binaries in Collisionally Relaxed Nuclear Star Clusters -- Impact of Mass Segregation
• Effective two-body scatterings around a massive object
• Short-Range Forces Can Catalyze Extreme Orbital Evolution in Hierarchical Triples