An Implementation of Stochastic Forces for the N-body code REBOUND
arXiv:2205.06757 · doi:10.3847/2515-5172/ac6e41
Abstract
We describe the implementation of a new module which can be used to simulate physical systems in which the motion of particles is affected by stochastic forces. Such forces are expected to be present in turbulent circumstellar disks or remnant planetesimal disks. Our implementation offers a convenient way to generate correlated noise with a user-specified amplitude and auto-correlation time for each particle. The module has minimal memory requirements and is freely available within the REBOUNDx additional effects library.
arXiv comment3 pages, 1 figure, accepted by RNAAS, source code available at https://github.com/dtamayo/reboundx
References in corpus (10)
- 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
- REBOUNDx: A Library for Adding Conservative and Dissipative Forces to Otherwise Symplectic N-body Integrations
- On the Orbital Evolution of Low Mass Protoplanets in Turbulent, Magnetised Disks
- Turbulent Torques on Protoplanets in a Dead Zone
- On the evolution of mean motion resonances through stochastic forcing: Fast and slow libration modes and the origin of HD128311
- Brownian Motion in Planetary Migration
- General Analysis of Type I Planetary Migration with Stochastic Perturbations
- A new paradigm for reproducing and analyzing N-body simulations of planetary systems