Revisiting gravitational instability in protostellar discs with improved radiative cooling models

This dataset contains the setup files and a subset of the simulation output data for the paper "Revisiting gravitational instability in protostellar discs with improved radiative cooling models" (Young et al. 2026). Files for the simulations using the modified Lombardi method are provided. The tar file name refers to the simulation parameters e.g. "M0.5_R100_Mu0.6.tar.gz" is for a stellar mass = 0.5 Msun, disc outer radius = 100au and mu=mdisc/mstar = 0.6. Data from the other simulations are available from the corresponding author. The simulations were produced with a modified version of Phantom (Price+ 2018) which is hosted on GitHub: https://github.com/alisonkyoung1/phantom Binary files are written in little endian format and can be read with Splash (see here: https://splash-viz.readthedocs.io/en/latest/. The remaining files are plain text files.