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

<p>MESA inlists, extras, and output for&nbsp;<a href="https://iopscience.iop.org/article/10.3847/1538-4357/ac3421">Evolutionary and Observational Consequences of Dyson Sphere Feedback,</a> published in the Astrophysical Journal, 2022.</p> <p>Paper&nbsp;DOI: 10.3847/1538-4357/ac3421</p> <p>Files included for Irradiated star experiments (recreating&nbsp;the work of Tout et al. (1989) in MESA) in irradiatedstar.zip:</p> <ul> <li>src/ <ul> <li>run_star_extras.f: MESA extension allowing for a constant temperature bath to irradiate a star.</li> </ul> </li> <li>inlists/ <ul> <li>inlist_zams: Create a zero age main sequence star.</li> <li>inlist_irrad_tams: Evolve a star from ZAMS to TAMS in a constant temperature irradiation&nbsp;bath.</li> <li>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.</li> </ul> </li> <li>plots/ <ul> <li>tout_recreation_plot.ipynb: Python Jupyter notebook for generating the Figure&nbsp;1&nbsp;in the paper.</li> <li>lifetimes.pdf: My MESA recreation of Figure 5 from Tout et al. (1989).</li> </ul> </li> <li>LOGS/ <ul> <li>*_history: MESA output for the evolution of a given stellar mass to a given age in a given temperature bath.</li> </ul> </li> </ul> <p>Files included for Dyson Sphere Feedback experiments&nbsp;in dsfeedback.zip:</p> <ul> <li>src/ <ul> <li>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.</li> </ul> </li> <li>inlists/ <ul> <li>inlist_sphere_tams: Evolve a selected model to the end of the main sequence&nbsp;with a given DS feedback fraction.</li> <li>inlist_zams: Create a zero age main sequence star.</li> <li>inlist_zams_lowmass:&nbsp;Create a zero age main sequence star with a timestep limit that helps with lower mass stars.</li> <li>inlist_zamstoage: Evolve a star normally from ZAMS to a selected age (0.6Gyr used here for 2Msun).</li> <li>inlist_zamstosolarage: Evolve a star normally from ZAMS to 4.6 Gyr (used before DS application for 1Msun and lower stars).</li> </ul> </li> <li>plots/ <ul> <li>ds_feedback_plots.ipynb: Python Jupyter notebook for generating the figures 2-9 in the paper.</li> <li>newplots/*: Plots generated from the notebook</li> </ul> </li> <li>LOGS/ <ul> <li>*_history: MESA output for the evolution of a given stellar mass, to a given age, and DS feedback level</li> <li>*_profile: MESA output for the structure of a given stellar mass, at one given age, and DS feedback level</li> </ul> </li> </ul> <p>Email me at mhuston@berkeley.edu if you have any comments or questions!</p> <p>MESA version: 12778</p> <p>MESA SDK version: x86_64-macos-20.4.1</p>