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.
hub
REBOUND: An open-source multi-purpose N-body code for collisional dynamics
35 Pith papers cite this work, alongside 1,022 external citations. Polarity classification is still indexing.
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.
hub tools
citation-role summary
citation-polarity summary
years
2026 35roles
background 3representative citing papers
Discovery of β Pic d, a 2-4 M_Jup planet at semi-major axis >30 au, via spectral template matching in JWST NIRSpec and MIRI data, confirmed by radial velocity and orbital stability.
SHARD is a collision-resolution framework that interpolates SPH results for remnants and clusters unresolved debris into N-body fragments while enforcing mass and water conservation for planet formation simulations.
Secular chaos regulated by the inner-to-outer ZLK timescale ratio R enables high-eccentricity tidal migration in 3+1 systems at sub-critical inclinations, producing polar hot Jupiters.
Nereid is proposed as Neptune's sole surviving regular satellite rather than a captured irregular satellite, based on new spectroscopic data and dynamical modeling.
Coordinated observations of a ~45 g earthgrazing fireball localize thermospheric shock generation to volatile-enhanced hydrodynamic shielding that enables continuum-like flow and detectable infrasound.
Realistic ISM fluctuations from TIGRESS-NCR MHD simulations drive stellar radial heating with σ_R ∝ t^{1/2} (early, cold) and t^{1/5} (late, warm) plus substantial migration explained by quasilinear theory with λ_* ~600 pc and τ_* ~70 Myr.
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.
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.
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 ratio distribution.
New Gemini and JWST data show 450P/LONEOS has a 1.8 km nucleus, CO2 gas production of ~7e24 molec/s, water ice absorption, and activity consistent with CO2 release from crystallizing amorphous ice.
Scattering between close-in super-Earths and secularly perturbed cold planets can launch some planets into free-floating orbits.
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 matching IAS15 accuracy at up to 9x speed.
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 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.
Two extremely low-density Jupiter-sized planets on long-period orbits around TOI-791 were confirmed via ground-based photometry and TTV-derived masses.
Simulations indicate EKL interactions excite TTV phases and disrupt resonances in a significant fraction of near-resonant super-Earth systems over 16 Myr.
Numerical integrations favor interpreting the Quaoar occultation feature as a satellite rather than a triangular-point confined arc, while identifying stable regions for additional moons.
New 2025 transit timing for HIP 41378 f confirms large TTVs and is combined with prior data on planets d and e in an N-body model to update ephemerides and predict future transits.
Dust particles inspiraling under PR drag are temporarily trapped in MMRs, escape due to overstability of librations, then undergo stochastic scattering with predicted fractions for planetary collision, sublimation, or ejection via an analytic epicycle model and Monte Carlo method.
Orbital calculations indicate that 1998 KY26 could be Phobos 1 if the probe executed two burns totaling 1.9 km/s after its 1988 loss.
N-body simulations show that 14 Herculis's orbital architecture requires primordial ejection of an additional massive planet.
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 destabilization near certain spacings.
N-body simulations demonstrate that post-capture chaotic planet-planet interactions around pulsars can produce stable low-eccentricity orbits after ejections.
citing papers explorer
-
Interstellar Medium-Driven Orbital Transport -- I. Radial Heating and Migration
Realistic ISM fluctuations from TIGRESS-NCR MHD simulations drive stellar radial heating with σ_R ∝ t^{1/2} (early, cold) and t^{1/5} (late, warm) plus substantial migration explained by quasilinear theory with λ_* ~600 pc and τ_* ~70 Myr.
-
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.
-
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.
-
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.