Code dependencies of pre-supernova evolution and nucleosynthesis in massive stars: Evolution to the end of core helium burning

TL;DR

This study compares massive star evolution models from three different codes up to core helium burning, highlighting how differences in convection treatment and mass loss affect predictions, and emphasizing the need for better modeling of stellar mixing processes.

Contribution

It provides a direct comparison of stellar evolution and nucleosynthesis predictions from GENEC, KEPLER, and MESA codes, focusing on the impact of physical assumptions on model outcomes.

Findings

Surface abundances agree reasonably across models.

Elemental overproduction factors vary less than 30%.

Differences mainly due to convection and mass loss treatments.

Abstract

Massive stars are key sources of radiative, kinetic, and chemical feedback in the universe. Grids of massive star models computed by different groups each using their own codes, input physics choices and numerical approximations, however, lead to inconsistent results for the same stars. We use three of these 1D codes---GENEC, KEPLER and MESA---to compute non-rotating stellar models of $15~\mathrm{M}_\odot$, $20~\mathrm{M}_\odot$, and $25~\mathrm{M}_\odot$ and compare their nucleosynthesis. We follow the evolution from the main sequence until the end of core helium burning. The GENEC and KEPLER models hold physics assumptions used in large grids of published models. The MESA code was set up to use convective core overshooting such that the CO core masses are consistent with those obtained by GENEC. For all models, full nucleosynthesis is computed using the NuGrid post-processing tool MPPNP. We find that the surface abundances predicted by the models are in reasonable agreement. In the helium core, the standard deviation of the elemental overproduction factors for Fe to Mo is less than $30\,\%$---smaller than the impact of the present nuclear physics uncertainties. For our three initial masses, the three stellar evolution codes yield consistent results. Differences in key properties of the models, e.g., helium and CO core masses and the time spent as a red supergiant, are traced back to the treatment of convection and, to a lesser extent, mass loss. The mixing processes in stars remain the key uncertainty in stellar modelling. Better constrained prescriptions are thus necessary to improve the predictive power of stellar evolution models.