Public code

The programs listed below are posted on GitHub, and you are welcome to use these codes for your own analyses. To download a program, follow the link provided below to access the GitHub page, then press the green ‘Code’ button, then ‘Download ZIP’. Unzip the downloaded files, read the README and you’re ready to go!

SculptingPlanet
Constrain the parameters of a single planet that sculpts the inner edge of a debris disc, using the model of Pearce et al. (2022). Download the Python code here.

GetDebrisStirringLevelInReboundSim
Measure the level of debris-disc stirring in a Rebound simulation using collisional analyses, as described in Costa, Pearce & Krivov (2024). A previous version contained a bug, which caused the collision velocity to be overestimated roughly 50% of the time. This bug was corrected on 26/01/2024. Download the Python code here.

GetSecularAndDiffusionTimes
Calculate the secular and diffusion timescales of a test particle, due to perturbations from a single planet. Download the Python code here.

GetForcedEAndIInTwoPlanetSystem
Calculate the forced eccentricity and forced longitude of pericentre on a test particle in a two-planet system. Download the Python code here.

GetTauMax
Calculate tauMax, the collisional lifetime of the largest bodies in a debris disc, based on the model of Löhne, Krivov & Rodmann (2008). By default the code uses the assumptions in Pearce et al. (2024). Download the Python code here.

GetLaplaceCoefficient
Calculate the Laplace coefficient b_s^(j)(α), which is used in dynamical calculations involving secular interactions or mean-motion resonances. Download the Python code here.

MinSelfStirringDiscMass
Calculate the absolute minimum mass that a debris disc requires to self-stir, using the model of Pearce et al. (2022). Download the Python code here.

PlanetMassToCarveDebrisGap
Calculate the gap carved by a migrating planet embedded in a massive debris disc, using the model of Friebe, Pearce & Löhne (2022). Download the Python code here.