Radial pulsation runaway in massive red supergiants in late evolutionary stage and implications to hydrogen-rich supernovae

This repository contains the simulation setup files associated with the paper:"Radial pulsation runaway in massive red supergiants in late evolutionary stage and implications for hydrogen-rich supernovae", published in Monthly Notices of the Royal Astronomical Society.The archive includes source code, inlist files, and related configuration files required to reproduce the simulations presented in the paper.In the original computations, the source codes are compiled by GNU Fortran (GCC) 12.1.0 on linux PC.All computations were performed using: - MESA version 24.08.1 - MESA SDK version 22.6.1 =DIRECTORY STRUCTUREThe archive consists of 11 directories, named M0130 through M0180, corresponding to initial stellar masses of 13–18 solar masses.Each directory contains the following subdirectories: 1.baseThis directory contains the source code and inlist files for the base models presented in the paper. To run a simulation: ./rn This produces the standard MESA outputs: history.data profile*.data Additional contents: - photos/ : original photo files produced in the simulations - nohup.out : saved standard output from the original run 2.model*These directories contain source and inlist files for short time-step simulations.Each directory includes: - A restart photo file in photos/ (e.g., M0130/model/photos/x00010640) This file is copied from the corresponding base model. These simulations restart from the photo with a reduced timestep limit specified in inlist_common.To run a simulation (e.g, in case of M0130/model1/): ./re x00010640 Outputs: - Standard MESA outputs (history.data, profile*.data) Notes: - Except for M0140 and M0170, each M0??? directory contains 8 models with different initial epochs. - Higher mass models likely to stop before reaching 100 years after the start of the simulation, owing to the developement of large-amplitude pulsations. Check the original standard output (nohup.out) in each directory. =RESOLUTION STUDY MODELS (M0140 and M0170 only)These directories include additional models to test numerical resolution: - dt2_model* Time-step increased by a factor of 2 relative to the fiducial setup - dt05_model* Time-step decreased by a factor of 2 - dm25_model* Higher spatial resolution with: delta_mesh_coefficient = 0.25 (fiducial value: 0.5) =PYTHON SCRIPTThe respository also contains the following plotting script written in Python. - evolution.py: reads and plots the time evolution of the properties of a base model - plot_surface_properties.py : plots the time evolutions of the properties of a short time-step model. - pulse_profile.py : performs period analysis and plots results =NOTES - Restart runs rely on provided photo files and modified timestep controls. - The included nohup.out files document the original execution logs for reproducibility.本仓库包含与发表于《皇家天文学会月报》（Monthly Notices of the Royal Astronomical Society）的论文《晚期演化阶段大质量红超巨星的径向脉动失控及其对富氢超新星的启示》相关的模拟配置文件。本归档包包含复现论文中所述模拟所需的源代码、inlist配置文件（inlist）及相关辅助配置文件。本次原始计算中，源代码于Linux个人计算机上通过GNU Fortran (GCC) 12.1.0编译完成。所有计算均基于以下工具完成： - MESA (MESA) 24.08.1版本 - MESA SDK (MESA SDK) 22.6.1版本 =目录结构 本归档包包含11个目录，命名为M0130至M0180，分别对应13至18倍太阳质量的初始恒星质量。每个目录均包含如下子目录： 1. 基础模型（base） 本目录包含论文中所述基础模型的源代码与inlist配置文件。 运行模拟的命令为： ./rn 运行后将生成标准MESA输出文件： history.data、profile*.data 额外内容包括： - photos/：模拟过程中生成的原始照片文件 - nohup.out：原始运行时保存的标准输出日志 2. 模型组（model*） 此类目录包含短时步模拟所需的源代码与inlist配置文件。每个目录均包含： - photos/目录下的重启照片文件（例如M0130/model/photos/x00010640），该文件源自对应基础模型。 此类模拟将从该照片文件处重启，并采用inlist_common中指定的缩减时步限制。以M0130/model1/为例，运行模拟的命令为： ./re x00010640 输出文件包括： - 标准MESA输出文件（history.data、profile*.data） 注意事项： - 除M0140与M0170外，每个M0???目录均包含8个不同初始时刻的模型。 - 质量更大的模型可能因大幅脉动的发展，在模拟启动后100年之前终止运行，请查阅各目录下的原始标准输出文件（nohup.out）以确认情况。 =分辨率测试模型（仅适用于M0140与M0170） 此类目录包含用于测试数值分辨率的额外模型： - dt2_model*：时步为基准配置的2倍 - dt05_model*：时步为基准配置的1/2 - dm25_model*：空间分辨率更高，设置为delta_mesh_coefficient = 0.25（基准值为0.5） =Python脚本 本仓库还包含以下Python编写的绘图脚本： - evolution.py：读取并绘制基础模型的属性随时间演化曲线 - plot_surface_properties.py：绘制短时步模型的属性随时间演化曲线 - pulse_profile.py：执行周期分析并绘制结果 =补充说明 - 重启运行依赖于提供的照片文件与修改后的时步控制参数。 - 附带的nohup.out文件记录了原始运行日志，以确保结果可复现。