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.

To more thoroughly study the effects of radiative forces on the orbits of small, astronomical bodies, we introduce the Yarkovsky effect into reboundx, an extensional library for the N-body integrator rebound. Two different versions of the Yarkovsky effect (the “Full Version” and the “Simple Version”) are available for use, depending on the needs of the user. We provide demonstrations for both versions of the effect and compare their computational efficiency with another previously implemented radiative force. In addition, we show how this effect can be used in tandem with other features in reboundx by simulating the orbits of asteroids during the asymptotic giant branch phase of a 2 M ☉ star. This effect is made freely available for use with the latest release of reboundx.