An Implementation of Stochastic Forces for the N-body Code REBOUND

Abstract

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.