An N-body model built from AMUSE used to calculate the short-lived radioisotopic enrichment from massive stars via winds and supernovae, as well as enrichment through interloping AGB stars. Much of this code is written to be as fast as possible, and has a few kludges here and there, I have done my best to be as explanitory as possible in the code-base. Some issues remain with the code such that it may need a bit of help setting up, feel free to message me by email.
The code is used in the following papers, all papers are open-access.
- Short-Lived Radioisotopic enrichment from AGB interlopers in low-mass star-forming regions - Eatson & Parker (2025) MNRAS, 544, 1 (343–354).
- Toward a unified injection model of short-lived radioisotopes in N-body simulations of star-forming regions - Eatson, Parker & Lichtenberg (2024) ApJ, 977, 1 (13).
Some important requirements, see requirements.txt for a full list, or install through the command pip3 -r requirements.txt.
amuse: N-body simulation (BHTree) and stellar evolution (SeBa).pandas: Results tabulation and processing, saving results data.py_ubjson: Saving simulation states as compressed binary files (much faster thanpickle).numpy: Ubiquitous mathemtatics library.numba: High-speed processing of data at C-like speeds for SLR enrichment in-line with N-body simulation with limited overhead.matplotlib: Data visualisation inal26_plotzstandard: Compression of data files.