The Yarkovsky effect in REBOUNDx

Abstract

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_{\odot }$ star. This effect is made freely available for use with the latest release of REBOUNDx.

arXiv commentAccepted by ApJS, 9 pages, 3 figures, source code available at https://github.com/dtamayo/reboundx, Jupyter notebooks for figures available at https://github.com/Nofe4108/REBOUNDx_Paper

References in corpus (13)

• Comprehensive analytic formulae for stellar evolution as a function of mass and metallicity
• 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
• WHFast: A fast and unbiased implementation of a symplectic Wisdom-Holman integrator for long term gravitational simulations
• Ancient planetary systems are orbiting a large fraction of white dwarf stars
• REBOUNDx: A Library for Adding Conservative and Dissipative Forces to Otherwise Symplectic N-body Integrations
• Eccentric Planets and Stellar Evolution as a Cause of Polluted White Dwarfs
• The orbital evolution of asteroids, pebbles and planets from giant branch stellar radiation and winds
• High order symplectic integrators for planetary dynamics and their implementation in REBOUND
• Speeding past planets? Asteroids radiatively propelled by giant branch Yarkovsky effects
• Detection of the Yarkovsky effect for C-type asteroids in the Veritas family
• JANUS: A bit-wise reversible integrator for N-body dynamics
• Embedded operator splitting methods for perturbed systems

Cited by in corpus (5)

• High-resolution resonant portraits of a single-planet white dwarf system
• Post-main sequence thermal evolution of planetesimals
• Dynamical feasibility of (3) Juno as a parent body of the H chondrites
• Planetesimals drifting through dusty and gaseous white dwarf debris discs: Types I, II and III-like migration
• An ancient L- type family associated to (460) Scania in the Middle Main Belt as revealed by Gaia DR3 spectra