Evolving massive stars to core collapse with GENEC: Extension of equation of state, opacities and effective nuclear network

TL;DR

This paper enhances the GENEC stellar evolution code to accurately model massive stars up to core collapse by upgrading physical ingredients, enabling improved pre-supernova models for supernova and remnant studies.

Contribution

The authors extend GENEC with improved equations of state, nuclear networks, and opacities, producing realistic pre-supernova models comparable to other codes like MESA.

Findings

Generated pre-supernova models for 15, 20, 25 solar mass stars.

Models are consistent with MESA and literature, with differences in shell structure.

Progenitors suitable for supernova and collapse simulations.

Abstract

Stars with initial mass above roughly 8 solar masses will evolve to form a core made of iron group elements at which point no further exothermic nuclear reactions between charged nuclei may prevent the core collapse. Electron captures, neutrino losses, and the photo-disintegration of heavy nuclei trigger the collapse of these stars. Models at the brink of core collapse are produced using stellar evolution codes and these pre-collapse models may be used in the study of the subsequent dynamical evolution (including their explosion as supernovae and the formation of compact remnants such as neutron stars or black holes). We upgrade the physical ingredients employed by the GENeva stellar Evolution Code, GENEC, so that it may cover the regime of high temperatures and high densities required to produce progenitors of core-collapse. We have improved GENEC in three directions, equation of state, the nuclear reaction network and the radiative and conductive opacities adapted for the computation of the advanced phases of evolution. We produce a small grid of pre-supernova models of stars with zero-age main sequence masses of 15, 20 and 25 solar masses at solar and less than half solar metallicities. The results are compared with analogous models produced with the MESA code. The global properties of our new models, particularly of their inner cores, are comparable to models computed with MESA and pre-existing progenitors in the literature. Between codes the exact shell structure varies impacting explosion predictions. Using GENEC with state-of-the-art physics, we have produced massive stellar progenitors prior to collapse. These progenitors are suitable for follow-up studies, including the dynamical collapse and supernova phases. Larger grids of supernova progenitors are now feasible, with potential for further dynamical evolution.