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.

Accepted 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 (10)

• Modules for Experiments in Stellar Astrophysics (MESA)
• Modules for Experiments in Stellar Astrophysics (MESA): Pulsating Variable Stars, Rotation, Convective Boundaries, and Energy Conservation
• Updated Electron-Conduction Opacities: The Impact on Low-Mass Stellar Models
• Hot Jupiters and Cool Stars
• Distant future of the Sun and Earth revisited
• Coulomb tunneling for fusion reactions in dense matter: Path integral Monte Carlo versus mean field
• Hybrid Symplectic Integrators for Planetary Dynamics
• 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
• Non-intrusive hierarchical coupling strategies for multi-scale simulations in gravitational dynamics