Raphael Hirschi

Grids of stellar models with rotation III. Models from 0.8 to 120 M at a metallicity Z = 0.002 ?

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 R136 star cluster hosts several stars whose individual masses greatly exceed the accepted 150 M⊙ stellar mass limit

Spectroscopic analyses of hydrogen-rich WN5‐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 MultiConjugate 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% of the ionizing flux of the entire 30 Doradus region. Compar isons 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 3603 A1. 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 mo re 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-instab ility supernovae could occur within the local universe, as has been claimed for SN 2007bi.

Stellar evolution with rotation - XII. Pre-supernova models

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 Mwere 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.

Grids of stellar models with rotation - II. WR populations and supernovae/GRB progenitors at Z = 0.014

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 –

Effects of rotation on the evolution of primordial stars

Context. Though still beyond our observational abilities, Populati on III stars are interesting objects in many perspectives. T hey are responsible for the re-ionisation of the inter-galactic me dium. They also left their chemical imprint in the early Universe, imprint which can be deciphered in the most metal-poor stars in the halo of our Galaxy. Aims. Rotation has been shown to play a determinant role at very low metallicity, bringing heavy mass loss where almost none was expected. Is this still true when the metallicity strictly e quals zero? The aim of our study is to get an answer to this question, and to determine how rotation changes the evolution and the chemical signature of the primordial stars. Methods. We have calculated seven differentially-rotating stellar models at zero metallicity, w ith masses between 9 and 200 M⊙. For each mass, we have also calculated a corresponding model without rotation. The evolution has been followed up to the pre-supernova stage. Results. We find that Z = 0 models rotate with an internal profile (r) close to local angular momentum conservation, because of a very weak core-envelope coupling. Rotational mixing drives a H-shell boost due to a sudden onset of CNO cycle in the shell. This boost leads to a high 14 N production, which can be as much as 10 6 times higher than the production of the non-rotating models. Generally, the rotating models produce much more metals than their non-rotating counterparts. The mass loss is very low, even for the models that reach the critical velocity during the main s equence. It may however have an impact on the chemical enrichment of the Universe, because some of the stars are supposed to collapse directly into black holes. They would contribute to the enrichment only through their winds. While in that case non-rotating st ars would not contribute at all, rotating stars may leave an i mprint in their surrounding. Due to the low mass loss and the weak coupling, the core retains a high angular momentum at the end of the evolution. The high rotation rate at death probably leads to a much stronger explosion than previously expected, changing the fate of the models. The inclusion of our yields in a chemical evolution model of the Galactic halo predicts log values of N/O, C/O and 12 C/ 13 C ratios of -2.2, -0.95 and 50 respectively at log O/H +12 = 4.2.

A strong case for fast stellar rotation at very low metallicities

We investigate the effect of new stellar models taking rotation into account and computing for a metallicity Z = 10 -8 on the chemical evolution of the earliest phases of the Milky Way. These models were computed under the assumption that the ratio of the initial rotation velocity to the critical velocity of stars is roughly constant with metallicity. This naturally leads to faster rotation at lower metallicity, as metal-poor stars are more compact than metal-rich ones. We find that the new Z = 10 -8 stellar yields have a tremendous impact on the nitrogen enrichment of the interstellar medium for log(O/H) + 12 < 7 (or [Fe/H] < -3). We show that including Z = 10 -8 stellar yields in chemical evolution models, both high N/O and C/O ratios are obtained in the very-metal poor metallicity range, in agreement with observations. Our results give further support to the idea that stars at very low metallicities could have rotational velocities of the order of 600-800 km s -1 .

Stellar evolution with rotation. XIII. Predicted GRB rates at various Z

We present the evolution of rotation in models of massive single stars covering a wide range of masses and metallicities. These models reproduce observations during the early stages of the evolution very well, in particular Wolf-Rayet (WR) populations and ratio between type II and type Ib,c supernovae at different metallicities. Our models predict the production of fast-rotating black holes. Models with large initial masses or high metallicity end their lives with less angular momentum in their central remnant with respect to the break-up limit for the remnant. Many WR star models satisfy the three main criteria (black hole formation, loss of hydrogen-rich envelope, and enough angular momentum to form an accretion disk around the black hole) for gamma-ray bursts (GRB) production via Woosleys collapsar model. If we consider all types of WR stars as GRB progenitors, there would be too many GRBs compared to observations but if we consider only WO stars (type Ic supernovae as is the case for SN2003dh/GRB030329) as GRB progenitors, the GRB production rates are in much better agreement with observations. WO stars are produced only at low metallicities in the present series of models. This prediction can be tested by future observations.

Evolution and fate of very massive stars

There is observational evidence that supports the existence of very massive stars (VMS) in the local universe. First, VMS (Mini ≲ 320 M⊙) have been observed in the Large Magellanic Clouds (LMC). Secondly, there are observed supernovae (SNe) that bear the characteristics of pair creation supernovae (PCSNe, also referred to as pair instability SN) which have VMS as progenitors. The most promising candidate to date is SN 2007bi. In order to investigate the evolution and fate of nearby VMS, we calculated a new grid of models for such objects, for solar, LMC and Small Magellanic Clouds (SMC) metallicities, which covers the initial mass range from 120 to 500 M⊙. Both rotating and non-rotating models were calculated using the GENEVA stellar evolution code and evolved until at least the end of helium burning and for most models until oxygen burning. Since VMS have very large convective cores during the main-sequence phase, their evolution is not so much affected by rotational mixing, but more by mass loss through stellar winds. Their evolution is never far from a homogeneous evolution even without rotational mixing. All the VMS, at all the metallicities studied here, end their life as WC(WO)-type Wolf-Rayet stars. Because of very important mass losses through stellar winds, these stars may have luminosities during the advanced phases of their evolution similar to stars with initial masses between 60 and 120 M⊙. A distinctive feature which may be used to disentangle Wolf-Rayet stars originating from VMS from those originating from lower initial masses would be the enhanced abundances of Ne and Mg at the surface of WC stars. This feature is however not always apparent depending on the history of mass loss. At solar metallicity, none of our models is expected to explode as a PCSN. At the metallicity of the LMC, only stars more massive than 300 M⊙ are expected to explode as PCSNe. At the SMC metallicity, the mass range for the PCSN progenitors is much larger and comprises stars with initial masses between about 100 and 290 M⊙. All VMS in the metallicity range studied here produce either a Type Ib SN or a Type Ic SN but not a Type II SN. We estimate that the progenitor of SN 2007bi, assuming a SMC metallicity, had an initial mass between 160 and 175 M⊙. None of models presented in this grid produces gamma-ray bursts or magnetars. They lose too much angular momentum by mass loss or avoid the formation of a black hole by producing a completely disruptive PCSN.

CONVECTIVE-REACTIVE PROTON-12C COMBUSTION IN SAKURAI'S OBJECT (V4334 SAGITTARII) AND IMPLICATIONS FOR THE EVOLUTION AND YIELDS FROM THE FIRST GENERATIONS OF STARS

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. Subject headings: stars: AGB and post-AGB — stars: abundances — stars: evolution — stars: interior — stars: individual (V4334 Sagittarii) — physical data and processes: hydrodynamics — physical data and processes: nuclear reactions, nucleosynthesis, abundances

Non-standard s-process in low metallicity massive rotating stars

Context. Rotation is known to affect the nucleosynthesis of light elements in massive stars, mainly by rotation-induced mixing. In particular, rotation boosts the primary nitrogen production. Models of rotating stars are able to reproduce the nitrogen observed in low-Z halo stars. Aims. Here we present the first grid of stellar models for rotating massive stars at low Z, where a full s-process network is used to study the impact of rotation-induced mixing on the nucleosynthesis of heavy elements. Methods. We used the Geneva stellar evolution code that includes an enlarged reaction network with nuclear species up to bismuth to calculate 25 M