Recognition: unknown
REBOUND: An open-source multi-purpose N-body code for collisional dynamics
read the original abstract
REBOUND is a new multi-purpose N-body code which is freely available under an open-source license. It was designed for collisional dynamics such as planetary rings but can also solve the classical N-body problem. It is highly modular and can be customized easily to work on a wide variety of different problems in astrophysics and beyond. REBOUND comes with three symplectic integrators: leap-frog, the symplectic epicycle integrator (SEI) and a Wisdom-Holman mapping (WH). It supports open, periodic and shearing-sheet boundary conditions. REBOUND can use a Barnes-Hut tree to calculate both self-gravity and collisions. These modules are fully parallelized with MPI as well as OpenMP. The former makes use of a static domain decomposition and a distributed essential tree. Two new collision detection modules based on a plane-sweep algorithm are also implemented. The performance of the plane-sweep algorithm is superior to a tree code for simulations in which one dimension is much longer than the other two and in simulations which are quasi-two dimensional with less than one million particles. In this work, we discuss the different algorithms implemented in REBOUND, the philosophy behind the code's structure as well as implementation specific details of the different modules. We present results of accuracy and scaling tests which show that the code can run efficiently on both desktop machines and large computing clusters.
This paper has not been read by Pith yet.
Forward citations
Cited by 16 Pith papers
-
Enhancement of the Rate of Tidal Disruption Events in Active Galactic Nuclei due to the Sweeping Secular Resonance Mechanism
Sweeping secular resonance from an intermediate-mass companion and depleting disk enhances AGN TDE rates to 10^{-3}-10^{-2} per galaxy per year on Myr timescales.
-
Trajectory-Agnostic Asteroid Detection in TESS with Deep Learning
A W-Net deep learning model detects asteroids in TESS data independently of trajectory by rotating training image cubes and using adaptive normalization for data scaling.
-
Suppression of Resonant Overstability at Sharp Migration Gradients
Sharp migration-rate gradients in protoplanetary disks quench resonant overstability when the dimensionless steepness parameter β exceeds the ratio of semi-major axis to eccentricity evolution timescales.
-
Impact of Cold Jupiter Scattering on the Mean-Motion Resonance of Inner Small Planets
Scattering cold Jupiters disrupt inner mean-motion resonances via secular perturbations from their orbital history, driving resonance circulation in most 2:1 and 3:2 configurations and explaining the Kepler period rat...
-
A Coordinate System for Dynamical Instabilities in Hierarchical Systems in REBOUND
A new wide-binary coordinate system in REBOUND's TRACE integrator produces qualitatively correct results for planet-planet scattering, stellar flybys, and ZLK oscillations where prior hybrid integrators fail, while ma...
-
Eccentricity as a signature of hierarchical subsolar-mass mergers in collapsar disks
Numerical relativity simulations of hierarchical mergers in collapsar-disk geometry show eccentricity building to 0.6 and partially surviving (up to 0.1) to merger, offering a gravitational-wave signature for this for...
-
Fates of the sub-stellar objects (FOSSO) II. Evidence for Suppression of Metal Pollution in White Dwarfs by Close Substellar Companions
Close substellar companions suppress metal pollution in white dwarfs by a factor of 5.75 with 87% protection efficiency, backed by observations of 17 systems and N-body simulations.
-
Planetesimal-Driven Instabilities in Resonant Chains of Cold Neptunes and Their Dynamical Outcomes
Planetesimal disks with 1-4% of the planetary mass disrupt resonant Neptune chains, triggering instabilities that scatter planets to ~0.1 au orbits and enable hot Neptune formation on 10-100 Myr timescales.
-
A Robust Launching Mechanism for Freely-Floating Planets from Host Stars with Close-in Planets
Scattering between close-in super-Earths and secularly perturbed cold planets can launch some planets into free-floating orbits.
-
Formation of stable exoplanetary systems around pulsars by capture: An exercise in computational classical mechanics
N-body simulations demonstrate that post-capture chaotic planet-planet interactions around pulsars can produce stable low-eccentricity orbits after ejections.
-
Life is But a Stream: The Distribution of Planetary Systems Along Stellar Streams and their Properties
Simulations show stars at the edges of stellar streams are more likely to retain unperturbed planetary systems than stars near the stream center.
-
Viscously Stirring Particle Disks into Lorentzians and Gaussians to Infer Dynamical and Collisional Masses (ARKS XIII)
Viscous stirring via gravitational scattering produces lognormal inclination distributions that yield Lorentzian vertical density profiles, which relax to Gaussians after equipartition, enabling estimates of perturber...
-
Terrestrial planet formation in the era of GPU computing
GPU-accelerated N-body simulations show that the common acceleration factor f distorts planetary chemical compositions and that terrestrial planets can form resonant chains without gas-driven orbital migration.
-
Orbital Stability of Closely-Spaced Four-planet Systems
Four-planet systems exhibit exponentially increasing lifetimes with orbital spacing, intermediate between three- and five-planet systems, with resonances causing shorter lifetimes and third-order MMRs adding destabili...
-
Observation Timelines for the Potential Lunar Impact of Asteroid 2024 YR4
Asteroid 2024 YR4 has a 4.3% chance of striking the Moon in 2032, producing a 1 km crater with a visible optical flash, hours of infrared afterglow, magnitude-5 seismic reverberation, and escaping debris that could re...
-
Comet 1P/Halley Completes 15 Orbits in 1,151 Years: Commensurability with the Solar System Quasi-Period and Evidence for Jupiter-Saturn Dynamical Coupling
Comet Halley exhibits a precise 15:1 commensurability with a 1,151-year solar-system quasi-period, with Jupiter and Saturn providing coherent perturbation cancellation over that baseline.
discussion (0)
Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.