MESA files, scripts, and data for: Self-Consistent Nonlinear Classical Cepheid Pulsations During Stellar Evolution with MESA

MESA model files, data, and plotting scripts for plots seen in 'Self-Consistent Nonlinear Classical 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