AccretR is a planetary body accretion Monte Carlo code that calculates mass, radius and bulk composition along a specified growth track, for orderly/hierarchical, runaway, and random particle accretion models. Elements in the model include concentrations of: H, C, N, O, Na, Mg, Al, Si, S, Cl, K, Ca, and Fe. Maximal water is also computed, assuming all H goes into forming water. Accretional heat is also calculated.
Currently, the code is optimized to build Jupiter's moon Europa, and Saturn's moons Titan and Enceladus, from CI, CM, CR, CK, CO and CV carbonaceous chondrite meteorites, cometary material (using comet 67P/Churyumov-Gerasimenko), and pure water ice. Literature sources for each of these compositions are found in the code.
Some care must be taken when selecting to accrete either Europa, Titan or Enceladus - simply comment out the unnecessary lines, and run in R. RStudio may be used. The code can run on a HPC, and has been run successfully on the Texas Advanced Computing Center (TACC) Stampede2 HPC (https://www.tacc.utexas.edu/systems/stampede2).
To do:
- Timing of accretion needs to be tweaked to allow a change in radiative heat dissipation rates (i.e., shedding accretional heat).
- Evaluation and quantification of melting and vaporization of maximal water reservoir after each particle impact (computationally intensive).
- Various quality of life improvements, such as making the selection of bulk compositions and growth tracks more intuitive.
- Rewrite as a R notebook.
- Change the multi-plotting functions to tidyverse plots.
- Adding new constraints on accretion material hydration and temperature derived from literature (e.g. hydrodynamic escape constraints, Bierson & Nimmo (2020). ApJ, 897(2), L43. https://doi.org/10.3847/2041-8213/aba11a).
Desiderata:
- Changing isotope reservoirs of stable isotopes.
- A pretty GUI.