MESA input/output for Evolutionary and Observational Consequences of Dyson Sphere Feedback

MESA inlists, extras, and output for Evolutionary and Observational Consequences of Dyson Sphere Feedback, published in the Astrophysical Journal, 2022. Paper DOI: 10.3847/1538-4357/ac3421 Files included for Irradiated star experiments (recreating the work of Tout et al. (1989) in MESA) in irradiatedstar.zip: src/ run_star_extras.f: MESA extension allowing for a constant temperature bath to irradiate a star. inlists/ inlist_zams: Create a zero age main sequence star. inlist_irrad_tams: Evolve a star from ZAMS to TAMS in a constant temperature irradiation bath. inlist_irrad_hefl: Attempt to evolve a star from TAMS to the helium flash in a constant temperature irradiation bath. Note: This model crashes before the helium flash for many parameter choices. plots/ tout_recreation_plot.ipynb: Python Jupyter notebook for generating the Figure 1 in the paper. lifetimes.pdf: My MESA recreation of Figure 5 from Tout et al. (1989). LOGS/ *_history: MESA output for the evolution of a given stellar mass to a given age in a given temperature bath. Files included for Dyson Sphere Feedback experiments in dsfeedback.zip: src/ run_star_extras.f: MESA extension allowing for a fraction of a star's luminosity to be dumped back onto the outer layer as extra heat. inlists/ inlist_sphere_tams: Evolve a selected model to the end of the main sequence with a given DS feedback fraction. inlist_zams: Create a zero age main sequence star. inlist_zams_lowmass: Create a zero age main sequence star with a timestep limit that helps with lower mass stars. inlist_zamstoage: Evolve a star normally from ZAMS to a selected age (0.6Gyr used here for 2Msun). inlist_zamstosolarage: Evolve a star normally from ZAMS to 4.6 Gyr (used before DS application for 1Msun and lower stars). plots/ ds_feedback_plots.ipynb: Python Jupyter notebook for generating the figures 2-9 in the paper. newplots/*: Plots generated from the notebook LOGS/ *_history: MESA output for the evolution of a given stellar mass, to a given age, and DS feedback level *_profile: MESA output for the structure of a given stellar mass, at one given age, and DS feedback level Email me at [e-mail] if you have any comments or questions!