Aims. Many topical astrophysical research areas, such as the properties of planet host stars, the nature of the progenitors of different types of supernovae and gamma ray bursts, and the evolution of galaxies, require complete and homogeneous sets of stellar models at different metallicities in order to be studied during the whole of cosmic history. We present here a first set of models for solar metallicity, where the effects of rotation are accounted for in a homogeneous way. Methods. We computed a grid of 48 different stellar evolutionary tracks, both rotating and non-rotating, at Z = 0.014, spanning a wide mass range from 0.8 to 120 M . For each of the stellar masses considered, electronic tables provide data for 400 stages along the evolutionary track and at each stage, a set of 43 physical data are given. These grids thus provide an extensive and detailed data basis for comparisons with the observations. The rotating models start on the zero-age main sequence (ZAMS) with a rotation rate υ ini /υ crit = 0.4. The evolution is computed until the end of the central carbon-burning phase, the early asymptotic giant branch (AGB) phase, or the core helium-flash for, respectively, the massive, intermediate, and both low and very low mass stars. The initial abundances are those deduced by Asplund and collaborators, which best fit the observed abundances of massive stars in the solar neighbourhood. We update both the opacities and nuclear reaction rates, and introduce new prescriptions for the mass-loss rates as stars approach the Eddington and/or the critical velocity. We account for both atomic diffusion and magnetic braking in our low-mass star models.Results. The present rotating models provide a good description of the average evolution of non-interacting stars. In particular, they reproduce the observed main-sequence width, the positions of the red giant and supergiant stars in the Hertzsprung-Russell (HR) diagram, the observed surface compositions and rotational velocities. Very interestingly, the enhancement of the mass loss during the red-supergiant stage, when the luminosity becomes supra-Eddington in some outer layers, help models above 15−20 M to lose a significant part of their hydrogen envelope and evolve back into the blue part of the HR diagram. This result has interesting consequences for the blue to red supergiant ratio, the minimum mass for stars to become Wolf-Rayet stars, and the maximum initial mass of stars that explode as type II−P supernovae.

Spectroscopic analyses of hydrogen-rich WN 5-6 stars within the young star clusters NGC 3603 and R136 are presented, using archival Hubble Space Telescope and Very Large Telescope spectroscopy, and high spatial resolution near-IR photometry, including Multi-Conjugate Adaptive Optics Demonstrator (MAD) imaging of R136. We derive high stellar temperatures for the WN stars in NGC 3603 (T * ∼ 42 ± 2 kK) and R136 (T * ∼ 53 ± 3 kK) plus clumping-corrected mass-loss rates of 2-5 × 10 −5 M yr −1 which closely agree with theoretical predictions from Vink et al. These stars make a disproportionate contribution to the global ionizing and mechanical wind power budget of their host clusters. Indeed, R136a1 alone supplies ∼7 per cent of the ionizing flux of the entire 30 Doradus region. Comparisons with stellar models calculated for the main-sequence evolution of 85-500 M accounting for rotation suggest ages of ∼1.5 Myr and initial masses in the range 105-170 M for three systems in NGC 3603, plus 165-320 M for four stars in R136. Our high stellar masses are supported by consistent spectroscopic and dynamical mass determinations for the components of NGC 3603A1. We consider the predicted X-ray luminosity of the R136 stars if they were close, colliding wind binaries. R136c is consistent with a colliding wind binary system. However, short period, colliding wind systems are excluded for R136a WN stars if mass ratios are of order unity. Widely separated systems would have been expected to harden owing to early dynamical encounters with other massive stars within such a high-density environment. From simulated star clusters, whose constituents are randomly sampled from the Kroupa initial mass function, both NGC 3603 and R136 are consistent with an tentative upper mass limit of ∼300 M . The Arches cluster is either too old to be used to diagnose the upper mass limit, exhibits a deficiency of very massive stars, or more likely stellar masses have been underestimated -initial masses for the most luminous stars in the Arches cluster approach 200 M according to contemporary stellar and photometric results. The potential for stars greatly exceeding 150 M within metal-poor galaxies suggests that such pair-instability supernovae could occur within the local universe, as has been claimed for SN 2007bi.

Aims. We study the impact of a subsolar metallicity on various properties of non-rotating and rotating stars, such as surface velocities and abundances, lifetimes, evolutionary tracks, and evolutionary scenarios. Methods. We provide a grid of single star models covering a mass range of 0.8 to 120 M with an initial metallicity Z = 0.002 with and without rotation. We discuss the impact of a change in the metallicity by comparing the current tracks with models computed with exactly the same physical ingredients but with a metallicity Z = 0.014 (solar). Results. We show that the width of the main-sequence (MS) band in the upper part of the Hertzsprung-Russell diagram (HRD), for luminosity above log (L/L ) > 5.5, is very sensitive to rotational mixing. Strong mixing significantly reduces the MS width. Here for the first time over the whole mass range, we confirm that surface enrichments are stronger at low metallicity provided that comparisons are made for equivalent initial mass, rotation, and evolutionary stage. We show that the enhancement factor due to a lowering of the metallicity (all other factors kept constant) increases when the initial mass decreases. Present models predict an upper luminosity for the red supergiants (RSG) of log (L/L ) around 5.5 at Z = 0.002 in agreement with the observed upper limit of RSG in the Small Magellanic Cloud. We show that models using shear diffusion coefficient, which is calibrated to reproduce the surface enrichments observed for MS B-type stars at Z = 0.014, can also reproduce the stronger enrichments observed at low metallicity. In the framework of the present models, we discuss the factors governing the timescale of the first crossing of the Hertzsprung gap after the MS phase. We show that any process favouring a deep localisation of the H-burning shell (steep gradient at the border of the H-burning convective core, low CNO content), and/or the low opacity of the H-rich envelope favour a blue position in the HRD for the whole, or at least a significant fraction, of the core He-burning phase.

Context. In recent years, many very interesting observations have appeared concerning the positions of Wolf-Rayet (WR) stars in the Hertzsprung-Russell diagram (HRD), the number ratios of WR stars, the nature of Type Ibc supernova (SN) progenitors, long and soft gamma ray bursts (LGRB), and the frequency of these various types of explosive events. These observations represent key constraints on massive star evolution. Aims. We study, in the framework of the single-star evolutionary scenario, how rotation modifies the evolution of a given initial mass star towards the WR phase and how it impacts the rates of Type Ibc SNe. We also discuss the initial conditions required to obtain collapsars and LGRB. Methods. We used a recent grid of stellar models computed with and without rotation to make predictions concerning the WR populations and the frequency of different types of core-collapse SNe. Current rotating models were checked to provide good fits to the following features: solar luminosity and radius at the solar age, main-sequence width, red-giant and red-supergiant (RSG) positions in the HRD, surface abundances, and rotational velocities. Results. Rotating stellar models predict that about half of the observed WR stars and at least half of the Type Ibc SNe may be produced through the single-star evolution channel. Rotation increases the duration of the WNL and WNC phases, while reducing those of the WNE and WC phases, as was already shown in previous works. Rotation increases the frequency of Type Ic SNe. The upper mass limit for Type II-P SNe is ∼19.0 M for the non rotating models and ∼16.8 M for the rotating ones. Both values agree with observations. Moreover, present rotating models provide a very good fit to the progenitor of SN 2008ax. We discuss future directions of research for further improving the agreement between the models and the observations. We conclude that the mass-loss rates in the WNL and RSG phases are probably underestimated at present. We show that up to an initial mass of 40 M , a surface magnetic field inferior to about 200 G may be sufficient to produce some braking. Much lower values are needed at the red supergiant stage. We suggest that the presence/absence of any magnetic braking effect may play a key role in questions regarding rotation rates of young pulsars and the evolution leading to LGRBs.

All ten LIGO/Virgo binary black hole (BH-BH) coalescences reported following the O1/O2 runs have near-zero effective spins. There are only three potential explanations for this. If the BH spin magnitudes are large, then: (i) either both BH spin vectors must be nearly in the orbital plane or (ii) the spin angular momenta of the BHs must be oppositely directed and similar in magnitude. Then there is also the possibility that (iii) the BH spin magnitudes are small. We consider the third hypothesis within the framework of the classical isolated binary evolution scenario of the BH-BH merger formation. We test three models of angular momentum transport in massive stars: a mildly efficient transport by meridional currents (as employed in the Geneva code), an efficient transport by the Tayler-Spruit magnetic dynamo (as implemented in the MESA code), and a very-efficient transport (as proposed by Fuller et al.) to calculate natal BH spins. We allow for binary evolution to increase the BH spins through accretion and account for the potential spin-up of stars through tidal interactions. Additionally, we update the calculations of the stellar-origin BH masses, including revisions to the history of star formation and to the chemical evolution across cosmic time. We find that we can simultaneously match the observed BH-BH merger rate density and BH masses and BH-BH effective spins. Models with efficient angular momentum transport are favored. The updated stellar-mass weighted gas-phase metallicity evolution now used in our models appears to be key for obtaining an improved reproduction of the LIGO/Virgo merger rate estimate. Mass losses during the pair-instability pulsation supernova phase are likely to be overestimated if the merger GW170729 hosts a BH more massive than 50 M⊙. We also estimate rates of black hole-neutron star (BH-NS) mergers from recent LIGO/Virgo observations. If, in fact. angular momentum transport in massive stars is efficient, then any (electromagnetic or gravitational wave) observation of a rapidly spinning BH would indicate either a very effective tidal spin up of the progenitor star (homogeneous evolution, high-mass X-ray binary formation through case A mass transfer, or a spin- up of a Wolf-Rayet star in a close binary by a close companion), significant mass accretion by the hole, or a BH formation through the merger of two or more BHs (in a dense stellar cluster).

Abstract.We describe the latest developments of the Geneva stellar evolution code in order to model the pre-supernova evolution of rotating massive stars. Rotating and non-rotating stellar models at solar metallicity with masses equal to 12, 15, 20, 25, 40 and 60 M were computed from the ZAMS until the end of the core silicon burning phase. We took into account meridional circulation, secular shear instabilities, horizontal turbulence and dynamical shear instabilities. We find that dynamical shear instabilities mainly smoothen the sharp angular velocity gradients but do not transport angular momentum or chemical species over long distances. Most of the differences between the pre-supernova structures obtained from rotating and non-rotating stellar models have their origin in the effects of rotation during the core hydrogen and helium burning phases. The advanced stellar evolutionary stages appear too short in time to allow the rotational instabilities considered in this work to have a significant impact during the late stages. In particular, the internal angular momentum does not change significantly during the advanced stages of the evolution. We can therefore have a good estimate of the final angular momentum at the end of the core helium burning phase. The effects of rotation on pre-supernova models are significant between 15 and 30 M . Indeed, rotation increases the core sizes (and the yields) by a factor ∼1.5. Above 20 M , rotation may change the radius or colour of the supernova progenitors (blue instead of red supergiant) and the supernova type (IIb or Ib instead of II). Rotation affects the lower mass limits for radiative core carbon burning, for iron core collapse and for black hole formation. For Wolf-Rayet stars (M > ∼ 30 M ), the pre-supernova structures are mostly affected by the intensities of the stellar winds and less by rotational mixing.

Depending on mass and metallicity as well as evolutionary phase, stars occasionally experience convectivereactive nucleosynthesis episodes. We specifically investigate the situation when nucleosynthetically unprocessed, H-rich material is convectively mixed with a He-burning zone, for example in convectively unstable shell on top of electron-degenerate cores in AGB stars, young white dwarfs or X-ray bursting neutron stars. Such episodes are frequently encountered in stellar evolution models of stars of extremely low or zero metal content, such as the first stars. We have carried out detailed nucleosynthesis simulations based on stellar evolution models and informed by hydrodynamic simulations. We focus on the convective-reactive episode in the very-late thermal pulse star Sakurai's object (V4334 Sagittarii). Asplund et al. (1999) determined the abundances of 28 elements, many of which are highly non-solar, ranging from H, He and Li all the way to Ba and La, plus the C isotopic ratio. Our simulations show that the mixing evolution according to standard, one-dimensional stellar evolution models implies neutron densities in the He intershell ( few 10 11 cm −3 ) that are too low to obtain a significant neutron capture nucleosynthesis on the heavy elements. We have carried out 3D hydrodynamic He-shell flash convection simulations in 4π geometry to study the entrainment of H-rich material. Guided by these simulations we assume that the ingestion process of H into the He-shell convection zone leads only after some delay time to a sufficient entropy barrier that splits the convection zone into the original one driven by He-burning and a new one driven by the rapid burning of ingested H. By making such mixing assumptions that are motivated by our hydrodynamic simulations we obtain significantly higher neutron densities (∼ few 10 15 cm −3 ) and reproduce the key observed abundance trends found in Sakurai's object. These include an overproduction of Rb, Sr and Y by about 2 orders of magnitude higher than the overproduction of Ba and La. Such a peculiar nucleosynthesis signature is impossible to obtain with the mixing predictions in our one-dimensional stellar evolution models. The simulated Li abundance and the isotopic ratio 12 C/ 13 C are as well in agreement with observations. Details of the observed heavy element abundances can be used as a sensitive diagnostic tool for the neutron density, for the neutron exposure and, in general, for the physics of the convective-reactive phases in stellar evolution. For example, the high elemental ratio Sc/Ca and the high Sc production indicate high neutron densities. The diagnostic value of such abundance markers depends on uncertain nuclear physics input. We determine how our results depend on uncertainties of nuclear reaction rates, for example for the 13 C(α, n) 16 O reaction. 10 Although even in this case multi-dimensional effects of convection have to be taken into account eventually as simulations by indicate that the velocity profile at the bottom of the convective shell is ...

The stellar mass range 8 M/M 12 corresponds to the most massive asymptotic giant branch (AGB) stars and the most numerous massive stars. It is host to a variety of supernova (SN) progenitors and is therefore very important for galactic chemical evolution and stellar population studies. In this paper, we study the transition from super-AGB (SAGB) star to massive star and find that a propagating neon-oxygen-burning shell is common to both the most massive electron capture supernova (EC-SN) progenitors and the lowest mass iron-core-collapse supernova (FeCCSN) progenitors. Of the models that ignite neon-burning off-center, the 9.5 M star would evolve to an FeCCSN after the neon-burning shell propagates to the center, as in previous studies. The neon-burning shell in the 8.8 M model, however, fails to reach the center as the URCA process and an extended (0.6 M ) region of low Y e (0.48) in the outer part of the core begin to dominate the late evolution; the model evolves to an EC-SN. This is the first study to follow the most massive EC-SN progenitors to collapse, representing an evolutionary path to EC-SN in addition to that from SAGB stars undergoing thermal pulses (TPs). We also present models of an 8.75 M SAGB star through its entire TP phase until electron captures on 20 Ne begin at its center and of a 12 M star up to the iron core collapse. We discuss key uncertainties and how the different pathways to collapse affect the pre-SN structure. Finally, we compare our results to the observed neutron star mass distribution.