Key points are not available for this paper at this time.
Helium-burning stars, in particular Cepheids, are especially difficult to model, as the choice of free parameters can greatly impact the shape of the blue loops - the part of the evolutionary track at which instability strip is crossed. Contemporary one-dimensional stellar evolution codes, like MESA (Modules for Experiments in Stellar Astrophysics), come with a large number of free parameters that allow to model the physical processes in stellar interiors under many assumptions. The uncertainties that arise from this freedom are rarely discussed in the literature despite their impact on the evolution of the model. We calculate a grid of evolutionary models with MESA, varying several controls, like solar mixture of heavy elements, mixing length theory prescription, nuclear reaction rates, scheme to determine convective boundaries, atmosphere model, temporal and spatial resolution, and quantify their impact on age and location of the evolutionary track on the HR diagram from the main sequence till the end of core-helium burning. Our investigation was conducted for a full range of masses and metallicities expected for classical Cepheids. The uncertainties are significant, especially during core-helium burning, reaching or exceeding the observational uncertainties of log Teff and log L for detached eclipsing binary systems. For models below 9 solar masses, thin convective shells develop and evolve erratically, not allowing the models to converge. A careful inspection of Kippenhahn diagrams and convergence study is advised for a given mass and metallicity, to assess how severe this problem is and to what extent it may affect the evolution.
Ziółkowska et al. (Tue,) studied this question.