Very Massive Star Models: I. Impact of Rotation and Metallicity and Comparisons with Observations

In addition to being spectacular objects, Very Massive Stars (VMS) are suspected to have a tremendous impact on their environment and on the whole cosmic evolution. The nucleosynthesis both during their advanced stages and their final explosion may contribute greatly to the overall enrichment of the Universe. Their resulting supernovae are candidates for the most superluminous events and their extreme conditions also lead to very important radiative and mechanical feedback effects, from local to cosmic scale. We explore the impact of rotation and metallicity on the evolution of very massive stars across cosmic times. With the recent implementation of an equation of state in the GENEC stellar evolution code, appropriate for describing the conditions in the central regions of very massive stars in the advanced phases, we present new results on VMS evolution from Population III to solar metallicity. Low metallicity VMS models are highly sensitive to rotation, while the evolution of higher metallicity models is dominated by mass loss effects. The mass loss affects strongly their surface velocity evolution, breaking quickly at high metallicity while reaching the critical velocity for low metallicity models. The comparison to observed VMS in the LMC shows that the mass loss prescriptions used for these models are compatible with observed mass loss rates. In our framework for modelling rotation, our models of VMS need a high initial velocity to reproduce the observed surface velocities. The surface enrichment of these VMS is difficult to explain with only one initial composition, and could suggest multiple populations in the R136 cluster. At a metallicity typical of R136, only our non- or slowly rotating VMS models may produce Pair Instability supernovae. The most massive black holes that can be formed are less massive than about 60 M$_\odot$.