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

TL;DR

This paper introduces a new module for the REBOUND N-body simulation code that models stochastic forces affecting particle motion, useful for simulating turbulent astrophysical disks with correlated noise.

Contribution

It provides a flexible, memory-efficient implementation of stochastic force modeling with customizable amplitude and auto-correlation time within REBOUND.

Findings

Enables simulation of turbulent disk environments with stochastic forces.

Offers correlated noise generation with user-defined parameters.

Available as a freely accessible module in REBOUNDx.

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.