An upper limit on the spins of merging binary black holes formed through binary evolution

Input files to reproduce simulations of the paper "An upper limit on the spins of merging binary black holes formed through binary evolution". Mesa version used was r23.05.01 with the MESA SDK x86_64-linux-22.6.1 for stripped star models and MESA version 11701 with MESA SDK x86_64-linux-20190503 for binary models (to reproduce old results). Files for MESA simulations are contained in the folder 'mesa_models', containing the following:

bin_models: Binary evolution models for a chemically homogeneous binary evolved with different models for angular momentum transport.

original: Original simulation from du Buisson et al. 2020 using the TS dynamo
TSF_Fuller_Lu2021: Same simulation but using the model for the TS dynamo of Fuller & Lu 2021
solid_body: Simulation done with solid body rotation
Z10Msun: stripped star models with metallicity Zsun/10. Each folder inside contains one model labeled by log10(M/Msun)
Z50Msun: stripped star models with metallicity Zsun/50. Each folder inside contains one model labeled by log10(M/Msun)
Z250Msun: stripped star models with metallicity Zsun/250. Each folder inside contains one model labeled by log10(M/Msun)
In addition, scripts to reproduce figures are included. These were done with julia version 1.9.2 and the Makie plotting package. A Manifest.toml and Project.toml are included to reproduce the environment used to produce the figures and analyze results. The included files are:

waveform.jl: Recalculates detection probabilities of du Buisson et al. 2020 including the effect of spin.
simplified_model.jl: Produces figures related to the semi-analytical model presented in the paper. Used for figures 1-4 of the paper
merger_rates.jl: Produces merger rate predictions. Used for figures 5-10 of the paper.
Extra folders included are:

ComputerModernFont: Fonts used for plotting
duBuisson_data: Original data from du Buisson et al. 2020. This includes both the outcomes of their mesa grids and their population synthesis calculations (see https://zenodo.org/records/3667546 for details).
sensitivity_curves: Sensitivity curves for GW detectors used to recompute detection probabilities