Stellar Evolution and Tidal Dissipation in REBOUNDx

Abstract

To study the post-main sequence evolution of the Solar system and exoplanetary systems more accurately and efficiently, we introduce two new features to REBOUNDx, an extended library for the N-body integrator REBOUND. The first is a convenient parameter interpolator for coupling different physics and integrators using numerical splitting schemes. The second implements a constant time lag model for tides without evolving spins. We demonstrate various uses of these features using stellar evolution data from MESA (Modules for Experiments in Stellar Astrophysics) as an example. The results of our tests agree with several studies in the literature on post-main sequence orbital evolution, and our convergence and performance studies respectively demonstrate our implementations' accuracy and efficiency. These additional effects are publicly available as of REBOUNDx's latest release.

arXiv commentAccepted by MNRAS, 9 pages, 6 figures, source code available at https://github.com/dtamayo/reboundx, Jupyter notebooks to reproduce figures available at https://github.com/sabaronett/REBOUNDxPaper

References in corpus (33)

• Modules for Experiments in Stellar Astrophysics (MESA)
• Modules for Experiments in Stellar Astrophysics (MESA): Giant Planets, Oscillations, Rotation, and Massive Stars
• Modules for Experiments in Stellar Astrophysics (MESA): Binaries, Pulsations, and Explosions
• MESA Isochrones and Stellar Tracks (MIST). I: Solar-Scaled Models
• Modules for Experiments in Stellar Astrophysics (MESA): Convective Boundaries, Element Diffusion, and Massive Star Explosions
• Modules for Experiments in Stellar Astrophysics (MESA): Pulsating Variable Stars, Rotation, Convective Boundaries, and Energy Conservation
• Low Temperature Opacities
• REBOUND: An open-source multi-purpose N-body code for collisional dynamics
• IAS15: A fast, adaptive, high-order integrator for gravitational dynamics, accurate to machine precision over a billion orbits
• Shell-model calculations of stellar weak interaction rates: II. Weak rates for nuclei in the mass range A=45-65 in supernovae environment
• Approximate input physics for stellar modelling
• Updated Electron-Conduction Opacities: The Impact on Low-Mass Stellar Models
• WHFast: A fast and unbiased implementation of a symplectic Wisdom-Holman integrator for long term gravitational simulations
• Tidal Decay of Close Planetary Orbits
• A New Version of Reimers' law of Mass Loss Based on a Physical Approach
• Thermodynamic functions of dense plasmas: analytic approximations for astrophysical applications
• REBOUNDx: A Library for Adding Conservative and Dissipative Forces to Otherwise Symplectic N-body Integrations
• Hot Jupiters and Cool Stars
• Distant future of the Sun and Earth revisited
• Simulations of two-planet systems through all phases of stellar evolution: implications for the instability boundary and white dwarf pollution
• Coulomb tunneling for fusion reactions in dense matter: Path integral Monte Carlo versus mean field
• Hybrid Symplectic Integrators for Planetary Dynamics
• Mercury-T: A new code to study tidally evolving multi-planet systems. Applications to Kepler-62
• The Solar System's Post-Main Sequence Escape Boundary
• Star-planet interactions. V. Dynamical and equilibrium tides in convective zones
• High order symplectic integrators for planetary dynamics and their implementation in REBOUND
• Second-order variational equations for N-body simulations
• An estimate of the $k_2$ Love number of WASP-18Ab from its radial velocity measurements
• The fates of Solar system analogues with one additional distant planet
• Dynamical Measurements of the Interior Structure of Exoplanets
• Non-intrusive hierarchical coupling strategies for multi-scale simulations in gravitational dynamics
• Computational astrophysics for the future: An open, modular approach with agreed standards would facilitate astrophysical discovery
• JANUS: A bit-wise reversible integrator for N-body dynamics