Self-Consistent Nonlinear Classical Cepheid Pulsations During Stellar Evolution with MESA
Abstract
We extend the time-dependent convection treatment in MESA by introducing eddy-viscous damping. This software change brings MESA-TDC into closer alignment with the radial stellar pulsation framework of MESA-RSP. We demonstrate that the inclusion of the eddy viscosity in hydrodynamic stellar models remains stable on evolutionary timescales. We then present the first self-consistent integration of large-amplitude, nonlinear Classical Cepheid pulsations directly within a MESA-star evolutionary run, demonstrating that the time-dependent convection formalism implemented in MESA-star and the MESA radial stellar pulsation (RSP) module are physically identical. Starting from a 6~\ blue-loop stellar evolution model, we demonstrate evolving the entire stellar model through pulsations as well as pausing the evolution, excising the core, and remeshing the envelope to match the grid used by MESA-RSP. We compare the pulsation properties (e.g., period, light and radius curves, and growth rate) with a matched MESA-RSP run, and find reasonable agreement between the two modules. This unified approach eliminates the reliance on separate post-processing workflows and enables fully coupled evolution-pulsation simulations. This approach enables future studies of stellar pulsations with the inclusion of composition gradients, mass loss, or rotation. It also enables future studies of the ε mechanism as well as providing a physical source of viscosity for other science cases explored using MESA's hydrodynamics solver. We have integrated these modifications into the MESA-star module, enabling open-source use by the community.
Turn this paper into a full lesson
ArcXiv compiles a staged curriculum from this paper: 8-12 lessons across beginner → advanced, synthesised section guides, visuals, flashcards, a quiz, exercises, and on-demand deep dives per section. Grounded in the abstract, never invented.