DOI: 10.5281/zenodo.17945923
title: MESA files, scripts, and data for: Self-Consistent Nonlinear Cepheid Pulsations During Stellar Evolution with MESA

description: MESA model files, data, and plotting scripts for plots seen in 'Self-Consistent Nonlinear Cepheid Pulsations During Stellar Evolution with MESA'. All models were run in the development version of MESA, and have been tested to function in the r25.12.1 release of MESA, with the latest mesasdk-aarch64-macos-25.12.2 +,  sdk available at http://user.astro.wisc.edu/~townsend/resource/download/mesasdk.

All plotting files and data are made available and can be running using python 3. Plots were generated with python 3.9.16.

(Note the colors module used for synthetic photometry in this work has since been replaced with a new colors module in MESA.)

A file tree is provided below:

├── data

│   ├── 6M_evol_histories

│   │   ├── alfam0.data

│   │   ├── alfam10.data

│   │   ├── alfam25.data

│   │   ├── alfam35.data

│   │   ├── alfam45.data

│   │   └── alfam55.data

│   ├── last_few_cycles_for_table

│   │   ├── RSP_6M_env_model.data

│   │   ├── TDC_6M_load_from_evol_env_model.data

│   │   ├── TDC_6M_load_from_evol_full_model.data

│   │   └── TDC_6M_load_from_rsp_env_model.data

│   ├── log_scrubber2.py

│   ├── RSP_Cepheid_6M

│   │   ├── history.data

│   │   ├── interval1_history.data

│   │   └── LOGS_pulsation

│   │       ├── history.data

│   │       ├── pgstar.dat

│           ├── profile*.data ...

│   ├── TDC_Cepheid_6M

│   │   └── history.data

│   ├── TDC_Cepheid_6M_evol

│   │   └── history.data

│   ├── TDC_Cepheid_6M_evol_no_uq

│   │   └── history.data

│   ├── TDC_Cepheid_6M_fm

│   │   ├── full_model.data

│   │   ├── history.data

│   │   ├── last_few_cycles_full_model.data

│   │   └── pgstar.dat

│   └── TDC_Cepheid_6M_load_from_rsp

│       ├── history.data

│       └── LOGS_pulsation

│           ├── history.data

│           ├── pgstar.dat

│           ├── profile*.data ...

├── mesa.mplstyle

├── models

│   ├── dev_TDC_Cepheid_6M

│   │   ├── ck

│   │   ├── clean

│   │   ├── evol.mp4

│   │   ├── gyre.data

│   │   ├── gyre.in

│   │   ├── he_dep.mod

│   │   ├── history_columns.list

│   │   ├── inlist

│   │   ├── inlist_common

│   │   ├── inlist_extra

│   │   ├── inlist_pgstar

│   │   ├── inlist_pulses

│   │   ├── inlist_pulses_header

│   │   ├── inlist_to_he_dep

│   │   ├── inlist_to_he_dep_header

│   │   ├── make

│   │   │   └── makefile

│   │   ├── mk

│   │   ├── profile_columns.list

│   │   ├── pulse.mp4

│   │   ├── re

│   │   ├── re_all

│   │   ├── re_mid_pulse

│   │   ├── re_nomodfiles

│   │   ├── README.rst

│   │   ├── rn

│   │   ├── rn_nomodfiles

│   │   ├── rn_pulse

│   │   ├── rn1

│   │   ├── run_all

│   │   ├── src

│   │   │   ├── run_star_extras_TDC_pulsation_defs.inc

│   │   │   ├── run_star_extras_TDC_pulsation.inc

│   │   │   ├── run_star_extras.f90

│   │   │   └── run.f90

│   │   └── standard_he_dep.mod

│   ├── dev_TDC_Cepheid_6M_alfam0p25_full_model

│   │   ├── cepheid_full_model_pulsation.mp4

│   │   ├── ck

│   │   ├── clean

│   │   ├── gyre.data

│   │   ├── gyre.in

│   │   ├── he_dep.mod

│   │   ├── history_columns.list

│   │   ├── inlist

│   │   ├── inlist_common

│   │   ├── inlist_extra

│   │   ├── inlist_pgstar

│   │   ├── inlist_pulses

│   │   ├── inlist_pulses_header

│   │   ├── inlist_to_he_dep

│   │   ├── inlist_to_he_dep_header

│   │   ├── make

│   │   │   └── makefile

│   │   ├── mk

│   │   ├── profile_columns.list

│   │   ├── re

│   │   ├── re_all

│   │   ├── re_mid_pulse

│   │   ├── re_nomodfiles

│   │   ├── README.rst

│   │   ├── rn

│   │   ├── rn_nomodfiles

│   │   ├── rn_pulse

│   │   ├── rn1

│   │   ├── run_all

│   │   ├── src

│   │   │   ├── run_star_extras_TDC_pulsation_defs.inc

│   │   │   ├── run_star_extras_TDC_pulsation.inc

│   │   │   ├── run_star_extras.f90

│   │   │   └── run.f90

│   │   └── standard_he_dep.mod

│   ├── dev_TDC_Cepheid_6M_load_from_RSP

│   │   ├── ck

│   │   ├── clean

│   │   ├── evol.mp4

│   │   ├── gyre.data

│   │   ├── gyre.in

│   │   ├── history_columns.list

│   │   ├── inlist

│   │   ├── inlist_common

│   │   ├── inlist_extra

│   │   ├── inlist_pgstar

│   │   ├── inlist_pulses

│   │   ├── inlist_pulses_header

│   │   ├── inlist_to_he_dep

│   │   ├── inlist_to_he_dep_header

│   │   ├── load_from_rsp_pulsation.mp4

│   │   ├── make

│   │   │   └── makefile

│   │   ├── mk

│   │   ├── profile_columns.list

│   │   ├── re

│   │   ├── re_all

│   │   ├── re_mid_pulse

│   │   ├── re_nomodfiles

│   │   ├── README.rst

│   │   ├── rn

│   │   ├── rn_nomodfiles

│   │   ├── rn_pulse

│   │   ├── rn1

│   │   ├── rsp_cepheid_6M_cycle0.mod

│   │   ├── run_all

│   │   └── src

│   │       ├── run_star_extras_TDC_pulsation_defs.inc

│   │       ├── run_star_extras_TDC_pulsation.inc

│   │       ├── run_star_extras.f90

│   │       └── run.f90

│   ├── rsp_Cepheid_6M_compare

│   │   ├── ck

│   │   ├── clean

│   │   ├── history_columns.list

│   │   ├── inlist

│   │   ├── inlist_rsp_Cepheid

│   │   ├── inlist_rsp_Cepheid_header

│   │   ├── inlist_rsp_pgstar_default

│   │   ├── make

│   │   │   └── makefile

│   │   ├── mk

│   │   ├── profile_columns.list

│   │   ├── re

│   │   ├── restart_photo

│   │   ├── rn

│   │   ├── rn1

│   │   ├── RSP_ceph_6M.mp4

│   │   └── src

│   │       ├── run_star_extras.f90

│   │       └── run.f90

│   └── rsp_Cepheid_6M_compare_make_mod

│       ├── ck

│       ├── clean

│       ├── history_columns.list

│       ├── inlist

│       ├── inlist_rsp_Cepheid

│       ├── inlist_rsp_Cepheid_header

│       ├── inlist_rsp_pgstar_default

│       ├── LOGS

│       │   ├── history.data

│       │   ├── LINA_eigen1.data

│       │   ├── LINA_eigen2.data

│       │   ├── LINA_eigen3.data

│       │   ├── LINA_period_growth.data

│       │   ├── LINA_work1.data

│       │   ├── LINA_work2.data

│       │   ├── LINA_work3.data

│       │   ├── pgstar.dat

│       │   ├── profile1.data

│       │   └── profiles.index

│       ├── make

│       │   ├── makefile

│       │   ├── run_star_extras.mod

│       │   ├── run_star_extras.o

│       │   ├── run_star_extras.smod

│       │   ├── run_star.mod

│       │   ├── run_star.o

│       │   └── run.o

│       ├── mk

│       ├── photos

│       │   ├── 1000

│       │   └── x192

│       ├── png

│       ├── profile_columns.list

│       ├── re

│       ├── rn

│       ├── rn1

│       ├── rsp_cepheid_6M_cycle0.mod

│       ├── src

│       │   ├── run_star_extras.f90

│       │   └── run.f90

│       ├── star

│       └── testhub.yml

├── plot_evol.py

├── plot_profiles.py

├── plot_pulsation_alfam.py

├── plot_pulsation.py

├── plots

│   ├── alfam_L_lightcurves.pdf

│   ├── alfam_R_lightcurves.pdf

│   ├── Combined_pulsation_panel.pdf

│   ├── compared_hr_evol_to_rsp.pdf

│   ├── compared_hr.pdf

│   ├── growth_rate_panel.pdf

│   ├── hr_evol_pulsation.pdf

│   ├── hr_evol.pdf

│   ├── hr_pulsation.pdf

│   ├── L_evol_lightcurves.pdf

│   ├── L_lightcurves.pdf

│   ├── logrho_phase_vs_radius_last_cycle.png

│   ├── logtau_phase_vs_radius_last_cycle.png

│   ├── opacity_phase_vs_logq_last_cycle.png

│   ├── opacity_phase_vs_radius_last_cycle.png

│   ├── phase_vs_logq_last_cycle.png

│   ├── phase_vs_radius_last_cycle.png

│   ├── profiles.pdf

│   ├── R_lightcurves.pdf

│   ├── RSP_pulsation_panel.pdf

│   ├── TDC_evol_pulsation_panel.pdf

│   ├── TDC_pulsation_panel.pdf

│   ├── TDC_vs_RSP_pulsation_panel.pdf

│   ├── y_face_phase_vs_logq_last_cycle.png

│   └── y_face_phase_vs_radius_last_cycle.png

└── profile_plots

    ├── 6M_load_from_rsp.data

    ├── 6M_rsp_pre_load_1192.data

    ├── 6M_rsp_pre_load_1193.data

    ├── evol_after_cut.data

    ├── evol_before_cut.data

    ├── mesa.mplstyle

    ├── plot_envelope_structure.py

    ├── profiles_logq.pdf

    └── profiles_rho_T.pdf

 

33 directories, 6503 files