Classical Cepheids pulsate — their radii, temperatures, and luminosities oscillate on periods of days to weeks. Stellar evolution codes model how stars change over millions of years. These two timescales differ by a factor of a billion. Until now, no code could run both simultaneously: the pulsation was calculated separately, as a post-processing step on a frozen evolutionary snapshot.
The MESA stellar evolution code now integrates large-amplitude nonlinear Cepheid pulsations directly within an evolutionary simulation. Starting from a 6-solar-mass blue-loop model, the code evolves the entire stellar structure — core nuclear burning, convective envelope, atmospheric oscillation — in a single coupled calculation. The pulsation periods, light curves, and radius variations match standalone pulsation calculations.
The enabling modification is surprisingly specific: adding eddy-viscous damping to MESA's time-dependent convection treatment. Without this term, the convective instabilities on short timescales corrupt the long-term evolutionary calculation. With it, the hydrodynamic models remain stable on evolutionary timescales while correctly capturing pulsation dynamics.
The unification matters because the features separated by the traditional workflow — composition gradients from nuclear burning, mass loss from stellar winds, rotational mixing — all operate on evolutionary timescales but affect pulsation properties. Running them in a single simulation means these couplings are captured naturally rather than approximated. The separated workflow assumed the pulsation didn't affect the evolution and vice versa. The coupled simulation tests that assumption — and enables investigation of cases where it fails.