This is the repository associated to the paper "Stellar Evolution in Real Time II: R Hydrae and an Open-Source Grid of >3000 Seismic TP-AGB Models Computed with MESA"
The most recent version of this paper will be available on arxiv on Friday, June 21st 2024. A pdf copy is included in this repository.
Due to the volume of data, the tracks are spread across four different Zenodo listings and several different .tar.gz files. Each Zenodo listing has its own DOI.
To use the smallest number of separate Zenodo listings possible, not all data are grouped in the same way.
The tar.gz files labeled with a string of the form "AGB_to_end_M1.XX_y-fixed.tar.gz" include all models having a mass beginning with 1 (1.00, 1.10, 1.20, etc) and adopting the fixed helium assumption.
Files labeled in the form "AGB_to_end_M1.XX_y-varied.tar.gz" include the same, but using the helium-varied assumption.
Files labeled with a string of the form "AGB_to_end_M3.10_all.tar.gz" include both the fixed-helium and varied-helium tracks for all masses that start with 3 (3.00, 3.10, 3.20, etc).
These data are not required to use the grid visualizer!! However, the grid visualizer uses lists of filenames generated based on these data, uploaded separately. If you wish to modify the underlying dataset, the list files must be regenerated.
visualization tools for grids of AGB models accompanying RHya paper
which loads the visualization tool in the browser
- three folders of realizations: yi-fixed_realizations, yi-varied_realizations, non-linear_realizations
all_best_fits*.dat
data files contain the best fit statistics as a function of mass, metallicity, and mode ID assumption (FM vs O1), computed according to different assumptions regarding agreement with period, luminosity, effective temperature, and radius. For example, the file all_best_fits_FM_Feb20_nonlin_yi_hardness75.dat
located in the non-linear_realizations/
folder will show a realization of the fundamental mode (FM) computed using the non-linear (nonlin) scaling relations on MESA models that have a fixed helium assumption and fit the pulse morphology according to a hardness parameter of 75.
See paper for more detail on the statistics.
a list of all models run to at least the onset of the thermally pulsing AGB phase, needed to generate certain aspects of the visualizer. The file names correspond verbatim to the MESA output data files that were used to make the pulse spectra figures in the visualizer and are available at https://meridithjoyce.com/pulse_data/
a list of png files needed in the visualizer. The filenames correspond verbatim to the files used in the visualizer.
This file can be reproduced by running the associated_plots.py
script in a directory containing the desired history data files
a crude average age in Myr taken by summing the model star's age at the onset of the TP-AGB phase and its age at the termination of the run and dividing by two.
script used to generate the pngs
file containing the composite observational measurements of period and period uncertainty as a function of (relative) time