All ten LIGO/Virgo binary black hole (BH-BH) coalescences reported from the O1/O2 runs have near zero effective spins. There are only three potential explanations of this fact. 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. Or, (iii) the BH spin magnitudes are small. We test 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, include revisions to the history of star formation and to the chemical evolution across cosmic time. We find that we can match simultaneously the observed BH-BH merger rate density, BH masses, and 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 a key in better reproducing the LIGO/Virgo merger rate estimate. Mass losses during the pair-instability pulsation supernova phase are likely overestimated if the merger GW170729 hosts a BH more massive than 50 Msun. We also estimate rate of BH-NS mergers from recent LIGO/Virgo observations. Our updated models of BH-BH, BH-NS and NS-NS mergers are now publicly available at www.syntheticuniverse.org.

The mass-loss rates of red supergiant stars (RSGs) are poorly constrained by direct measurements, and yet the subsequent evolution of these stars depends critically on how much mass is lost during the RSG phase. In 2012 the Geneva evolutionary group updated their mass-loss prescription for RSGs with the result that a 20 solar mass star now loses 10x more mass during the RSG phase than in the older models. Thus, higher mass RSGs evolve back through a second yellow supergiant phase rather than exploding as Type II-P supernovae, in accord with recent observations (the so-called "RSG Problem"). Still, even much larger mass-loss rates during the RSG phase cannot be ruled out by direct measurements of their current dust-production rates, as such mass-loss is episodic. Here we test the models by deriving a luminosity function for RSGs in the nearby spiral galaxy M31 which is sensitive to the total mass loss during the RSG phase. We carefully separate RSGs from asymptotic giant branch stars in the color-magnitude diagram following the recent method exploited by Yang and collaborators in their Small Magellanic Cloud studies. Comparing our resulting luminosity function to that predicted by the evolutionary models shows that the new prescription for RSG mass-loss does an excellent job of matching the observations, and we can readily rule out significantly larger values.

Massive stars with solar metallicity lose important amounts of rotational angular momentum through their winds. When a magnetic field is present at the surface of a star, efficient angular momentum losses can still be achieved even when the mass-loss rate is very modest, at lower metallicities, or for lower-initial-mass stars. In a close binary system, the effect of wind magnetic braking also interacts with the influence of tides, resulting in a complex evolution of rotation. We study the interactions between the process of wind magnetic braking and tides in close binary systems. We discuss the evolution of a 10 M$_\odot$ star in a close binary system with a 7 M$_\odot$ companion using the Geneva stellar evolution code. The initial orbital period is 1.2 days. The 10 M$_\odot$ star has a surface magnetic field of 1 kG. Various initial rotations are considered. We use two different approaches for the internal angular momentum transport. In one of them, angular momentum is transported by shear and meridional currents. In the other, a strong internal magnetic field imposes nearly perfect solid-body rotation. The evolution of the primary is computed until the first mass-transfer episode occurs. The cases of different values for the magnetic fields and for various orbital periods and mass ratios are briefly discussed. We show that, independently of the initial rotation rate of the primary and the efficiency of the internal angular momentum transport, the surface rotation of the primary will converge, in a time that is short with respect to the main-sequence lifetime, towards a slowly evolving velocity that is different from the synchronization velocity. (abridged).

We present Hubble Space Telescope photometry of NGC1850, a ~100 Myr, ~10^5 Msun cluster in the Large Magellanic Cloud. The colour magnitude diagram clearly shows the presence of an extended main sequence turnoff (eMSTO). The use of non-rotating stellar isochrones leads to an age spread of ~40 Myr. This is in good agreement with the age range expected when the effects of rotation in MSTO stars are wrongly interpreted in terms of age spread. We also do not find evidence for multiple, isolated episodes of star-formation bursts within the cluster, in contradiction to scenarios that invoke actual age spreads to explain the eMSTO phenomenon. NGC 1850 therefore continues the trend of eMSTO clusters where the inferred age spread is proportional to the age of the cluster. While our results confirm a key prediction of the scenario where stellar rotation causes the eMSTO feature, direct measurements of the rotational rate of MSTO stars is required to definitively confirm or refute whether stellar rotation is the origin of the eMSTO phenomenon or if it is due to an as yet undiscovered effect.

After a brief recall of the main impacts of stellar rotation on the structure and the evolution of stars, four topics are addressed: 1) the links between magnetic fields and rotation; 2) the impact of rotation on the age determination of clusters; 3) the exchanges of angular momentum between the orbit of a planet and the star due to tides; 4) the impact of rotation on the early chemical evolution of the Milky Way and the origin of the Carbon-Enhanced-Metal-Poor stars.

We investigate the impact of tidal interactions, before any mass transfer, on various properties of the stellar models. We study the conditions for obtaining homogeneous evolution triggered by tidal interactions, and for avoiding any Roche lobe overflow during the Main-Sequence phase. We consider the case of rotating stars computed with a strong coupling mediated by an interior magnetic field. In models without any tidal interaction (single stars and wide binaries), homogeneous evolution in solid body rotating models is obtained when two conditions are realized: the initial rotation must be high enough, the loss of angular momentum by stellar winds should be modest. This last point favors metal-poor fast rotating stars. In models with tidal interactions, homogeneous evolution is obtained when rotation imposed by synchronization is high enough (typically a time-averaged surface velocities during the Main-Sequence phase above 250 km s$^{-1}$), whatever the mass losses. In close binaries, mixing is stronger at higher than at lower metallicities. Homogeneous evolution is thus favored at higher metallicities. Roche lobe overflow avoidance is favored at lower metallicities due to the fact that stars with less metals remain more compact. We study also the impact of different processes for the angular momentum transport on the surface abundances and velocities in single and close binaries. In models where strong internal coupling is assumed, strong surface enrichments are always associated to high surface velocities in binary or single star models. In contrast, models computed with mild coupling may produce strong surface enrichments associated to low surface velocities. Close binary models may be of interest for explaining homogeneous massive stars, fast rotating Wolf-Rayet stars, and progenitors of long soft gamma ray bursts, even at high metallicities.

We present a new set of synthesis models for stellar populations obtained with Starburst99, which are based on new stellar evolutionary tracks with rotation. We discuss models with zero rotation velocity and with velocities of 40% of the break-up velocity on the zero-age main-sequence. These values are expected to bracket realistic rotation velocity distributions in stellar populations. The new rotating models for massive stars are more luminous and hotter due to a larger convective core and enhanced surface abundances. This results in pronounced changes in the integrated spectral energy distribution of a population containing massive stars. The changes are most significant at the shortest wavelengths where an increase of the ionizing luminosity by up to a factor of 5 is predicted. We also show that high equivalent widths of recombination lines may not necessarily indicate a very young age but can be achieved at ages as late as 10 Myr. Comparison of these two boundary cases (0 and 40% of the break-up velocity) will allow users to evaluate the effects of rotation and provide guidance for calibrating the stellar evolution models. We also introduce a new theoretical ultraviolet spectral library built from the Potsdam Wolf-Rayet (PoWR) atmospheres. Its purpose is to help identify signatures of Wolf-Rayet stars in the ultraviolet whose strength is sensitive to the particulars of the evolution models. The new models are available for solar and 1/7th solar metallicities. A complete suite of models can be generated on the Starburst99 website (www.stsci.edu/science/starburst99/). The updated Starburst99 package can be retrieved from this website as well.

For the first time, the interior and spectroscopic evolution of a massive star is analyzed from the zero-age main sequence (ZAMS) to the pre-supernova (SN) stage. For this purpose, we combined stellar evolution models using the Geneva code and atmospheric models using CMFGEN. With our approach, we were able to produce observables, such as a synthetic high-resolution spectrum and photometry, aiding the comparison between evolution models and observed data. Here we analyze the evolution of a non-rotating 60 Msun star and its spectrum throughout its lifetime. Interestingly, the star has a supergiant appearance (luminosity class I) even at the ZAMS. We find the following evolutionary sequence of spectral types: O3 I (at the ZAMS), O4 I (middle of the H-core burning phase), B supergiant (BSG), B hypergiant (BHG), hot luminous blue variable (LBV; end of H-core burning), cool LBV (H-shell burning through the beginning of the He-core burning phase), rapid evolution through late WN and early WN, early WC (middle of He-core burning), and WO (end of He-core burning until core collapse). We find the following spectroscopic phase lifetimes: 3.22e6 yr for the O-type, 0.34e5 yr (BSG), 0.79e5 yr (BHG), 2.35e5 yr (LBV), 1.05e5 yr (WN), 2.57e4 yr (WC), and 3.80e4 yr (WO). Compared to previous studies, we find a much longer (shorter) duration for the early WN (late WN) phase, as well as a long-lived LBV phase. We show that LBVs arise naturally in single-star evolution models at the end of the MS when the mass-loss rate increases as a consequence of crossing the bistability limit. We discuss the evolution of the spectra, magnitudes, colors, and ionizing flux across the star's lifetime, and the way they are related to the evolution of the interior. [abridged]

Core-collapse supernova (SN) explosions mark the end of the tumultuous life of massive stars. Determining the nature of their progenitors is a crucial step towards understanding the properties of SNe. Until recently, no progenitor has been directly detected for SN of type Ibc, which are believed to come from massive stars that lose their hydrogen envelope through stellar winds and from binary systems where the companion has stripped the H envelope from the primary. Here we analyze recently reported observations of iPTF13bvn, which could possibly be the first detection of a SN Ib progenitor based on pre-explosion images. Very interestingly, the recently published Geneva models of single stars can reproduce the observed photometry of the progenitor candidate and its mass-loss rate, confirming a recently proposed scenario. We find that a single WR star with initial mass in the range 31-35 Msun fits the observed photometry of the progenitor of iPTF13bvn. The progenitor likely has a luminosity of log (L/Lsun)~5.55, surface temperature ~45000 K, and mass of ~10.9 Msun at the time of explosion. Our non-rotating 32 Msun model overestimates the derived radius of the progenitor, although this could likely be reconciled with a fine-tuned model of a more massive (between 40 and 50 Msun), hotter, and luminous progenitor. Our models indicate a very uncertain ejecta mass of ~8 Msun, which is higher than the average of the SN Ib ejecta mass that is derived from the lightcurve (2-4 Msun). This possibly high ejecta mass could produce detectable effects in the iPTF13bvn lightcurve and spectrum. If the candidate is indeed confirmed to be the progenitor, our results suggest that stars with relatively high initial masses (>30 Msun) can produce visible SN explosions at their deaths and do not collapse directly to a black hole.

We discuss four questions dealing with massive star evolution. The first one is about the origin of slowly rotating, non-evolved, nitrogen rich stars. We propose that these stars may originate from initially fast rotating stars whose surface has been braked down. The second question is about the evolutionary status of alpha-Cygni variables. According to their pulsation properties, these stars should be post red supergiant stars. However, some stars at least present surface abundances indicating that they should be pre red supergiant stars. How to reconcile these two contradictory requirements? The third one concerns the various supernova types which are the end point of the evolution of stars with initial masses between 18 and 30 Msun, i.e. the most massive stars which go through a red supergiant phase during their lifetime. Do they produce types IIP, IIL, IIn, IIb or Ib supernovae or do they end without producing any SN event? Finally, we shall discuss reasons why so few progenitors of type Ibc supernovae have yet been detected?

We investigate the fundamental properties of core-collapse Supernova (SN) progenitors from single stars at solar metallicity. We combine Geneva stellar evolutionary models with initial masses of Mini=20-120 Msun with atmospheric/wind models using CMFGEN. We provide synthetic photometry and high-resolution spectra of hot stars at the pre-SN stage. For Mini=9-20 Msun, we supplement our analysis using publicly available MARCS model atmospheres of RSGs. We employ observational criteria of spectroscopic classification and find that massive stars, depending on Mini and rotation, end their lives as red supergiants (RSG), yellow hypergiants (YHG), luminous blue variables (LBV), and Wolf-Rayet (WR) stars of the WN and WO spectral types. For rotating models, we obtain the following types of SN progenitors: WO1-3 (Mini <= 32 Msun), WN10-11 (25 < Mini < 32 Msun), LBV (20 <= Mini < 25 Msun), G1 Ia+ (18 < Mini < 20 Msun), and RSGs (9 <= Mini <= 18 Msun). For non-rotating models, we find spectral types WO1-3 (Mini>40 Msun), WN7-8 (25 < Mini <= 40 Msun), WN11h/LBV (20 < Mini <= 25 Msun), and RSGs (9 <= Mini <= 20 Msun). Our rotating models indicate that SN IIP progenitors are all RSG, SN IIL/b progenitors are 56% LBVs and 44% YHGs, SN Ib progenitors are 96% WN10-11 and 4% WOs, and SN Ic progenitors are all WO stars. We find that not necessarily the most massive and luminous SN progenitors are the brighter ones in a given filter. We show that SN IIP progenitors (RSGs) are bright in the RIJHK_S filters and faint in the UB filters. SN IIL/b progenitors (LBVs and YHGs), and SN Ib progenitors (WNs) are relatively bright in optical/IR filters, while SN Ic progenitors (WOs) are faint in all optical filters. We argue that SN Ib and Ic progenitors from single stars should be undetectable in the available pre-explosion images with the current magnitude limits, in agreement with observational results.[abridged]

There is observational evidence that supports the existence of Very Massive Stars in the local universe. First, very massive stars (Mini<=320 M) have been observed in the Large Magellanic Cloud . Second, there are observed SNe that bear the characteristics of Pair Creation Supernovae which have very massive stars as progenitors. The most promising candidate to date is SN2007bi. In order to investigate the evolution and fate of nearby very massive stars, we calculated a new grid of models for such objects, for solar, LMC and SMC metallicities, which covers the initial mass range from 120 to 500M. 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 very massive stars 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 stars. 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 stars in the metallicity range studied here produce either a type Ib or a type Ic SN. We estimate that the progenitor of SN2007bi, assuming a SMC metallicity, had an initial mass between 160 and 175 M . None of models presented in this grid produce GRBs or magnetars. They lose too much angular momentum by mass loss or avoid the formation of a BH by producing a completely disruptive PCSN.

We discuss the physics and the evolution of a typical massive star passing through an evolutionary stage similar to that of Betelgeuse. After a brief introduction recalling various observed parameters of Betelgeuse, we discuss the Pre-Main-Sequence phase (PMS), the Main-Sequence (MS) phase, the physics governing the duration of the first crossing of the HR diagram, the red supergiant stage (RSG), the post-red supergiant phases and the final fate of solar metallicity stars with masses between 9 and 25 M$_\odot$. We examine the impact of different initial rotation and of various prescriptions for the mass loss rates during the red supergiant phase. We show that, whatever the initial rotation rate (chosen between 0 and 0.7$\times\upsilon_{\rm crit}$, $\upsilon_{\rm crit}$ being the surface equatorial velocity producing a centrifugal acceleration balancing exactly the gravity) and the mass loss rates during the RSG stage (varied between a standard value and 25 times that value), a 15 M$_\odot$ star always ends its lifetime as a RSG and explodes as a type II-P or II-L supernova.

Optical and infrared interferometers definitively established that the photometric standard Vega (alpha Lyrae) is a rapidly rotating star viewed nearly pole-on. Recent independent spectroscopic analyses could not reconcile the inferred inclination angle with the observed line profiles, preferring a larger inclination. In order to resolve this controversy, we observed Vega using the six-beam Michigan Infrared Combiner on the Center for High Angular Resolution Astronomy Array. With our greater angular resolution and dense (u,v)-coverage, we find Vega is rotating less rapidly and with a smaller gravity darkening coefficient than previous interferometric results. Our models are compatible with low photospheric macroturbulence and also consistent with the possible rotational period of ~0.71 days recently reported based on magnetic field observations. Our updated evolutionary analysis explicitly incorporates rapid rotation, finding Vega to have a mass of 2.15+0.10_-0.15 Msun and an age 700-75+150 Myrs, substantially older than previous estimates with errors dominated by lingering metallicity uncertainties (Z=0.006+0.003-0.002).

We present a sample of synthetic massive stellar populations created using the Starburst99 evolutionary synthesis code and new sets of stellar evolutionary tracks, including one set that adopts a detailed treatment of rotation. Using the outputs of the Starburst99 code, we compare the populations' integrated properties, including ionizing radiation fields, bolometric luminosities, and colors. With these comparisons we are able to probe the specific effects of rotation on the properties of a stellar population. We find that a population of rotating stars produces a much harder ionizing radiation field and a higher bolometric luminosity, changes that are primarily attributable to the effects of rotational mixing on the lifetimes, luminosities, effective temperatures, and mass loss rates of massive stars. We consider the implications of the profound effects that rotation can have on a stellar population, and discuss the importance of refining stellar evolutionary models for future work in the study of extragalactic, and particularly high-redshift, stellar populations.

We present dense grids of stellar models suitable for comparison with observable quantities measured with great precision, such as those derived from binary systems or planet-hosting stars. We computed new Geneva models without rotation at metallicities Z=0.006, 0.01, 0.014, 0.02, 0.03 and 0.04 (i.e. [Fe/H] from -0.33 to +0.54) and with mass in small steps from 0.5 to 3.5 Msun. Great care was taken in the procedure for interpolating between tracks in order to compute isochrones. Several properties of our grids are presented as a function of stellar mass and metallicity. Those include surface properties in the Hertzsprung-Russell diagram, internal properties including mean stellar density, sizes of the convective cores, and global asteroseismic properties. We checked our interpolation procedure and compared interpolated tracks with computed tracks. The deviations are less than 1% in radius and effective temperatures for most of the cases considered. We also checked that the present isochrones provide nice fits to four couples of observed detached binaries and to the observed sequences of the open clusters NGC 3532 and M67. Including atomic diffusion in our models with M<1.1 Msun leads to variations in the surface abundances that should be taken into account when comparing with observational data of stars with measured metallicities. For that purpose, iso-Zsurf lines are computed. These can be requested for download from a dedicated web page together with tracks at masses and metallicities within the limits covered by the grids. The validity of the relations linking Z and FeH is also re-assessed in light of the surface abundance variations in low-mass stars.

The current state-of-the-art of population synthesis is reviewed. The field is currently undergoing major revisions with the recognition of several key processes as new critical ingredients. Stochastic effects can artificially enhance or suppress certain evolutionary phases and/or stellar mass regimes and introduce systematic biases in, e.g., the determination of the stellar initial mass function. Post-main-sequence evolution is often associated with irregular variations of stellar properties on ultra-short time-scales. Examples are asymptotic giant branch stars and luminous blue variables, both of which are poorly treated in the models. Stars rarely form in isolation, and the fraction of truly single stars may be very small. Therefore, stellar multiplicity must be accounted for since many systems will develop tidal interaction over the course of their evolution. Last but not least, stellar rotation can drastically increase stellar temperatures and luminosities, which in turn leads to revised mass-to-light ratios in population synthesis models.

According to recent theoretical studies, the progenitors of Long Gamma Ray Bursts should be very fast rotating stars, massive enough but not so for collapsing into a black hole. In addition, recent observations seem to show that stars of about 20 solar masses could be at the origin of LGRBs. At low metallicity B-type stars rotate faster than at higher metallicity. We found with the ESO-WFI an occurrence of Be/Oe stars, that are quasi critical rotators, 3 to 5 times larger in the SMC than in the Galaxy. According to our results, and using observational clues on the SMC WR stars, as well as the theoretical predictions of the characteristics must have the LGRB progenitors, we have identified the low metallicity massive Be/Oe stars as potential LGRB progenitors. To support this identification, the expected rates and the numbers of LGRB were then calculated and compared to the observed ones: 3 to 6 LGRBs were found in the local universe in 11 years while 8 were actually observed.

Since the external regions of the envelopes of rapidly rotating early-type stars are unstable to convection, a coupling may exist between the convection and the internal rotation. We explore what can be learned from spectroscopic and interferometric observations about the properties of the rotation law in the external layers of these objects. Using simple relations between the entropy and specific rotational quantities, some of which are found to be efficient at accounting for the solar differential rotation in the convective region, we derived analytical solutions that represent possible differential rotations in the envelope of early-type stars. A surface latitudinal differential rotation may not only be an external imprint of the inner rotation, but induces changes in the stellar geometry, the gravitational darkening, the aspect of spectral line profiles, and the emitted spectral energy distribution. By studying the equation of the surface of stars with non-conservative rotation laws, we conclude that objects undergo geometrical deformations that are a function of the latitudinal differential rotation able to be scrutinized both spectroscopically and by interferometry. The combination of Fourier analysis of spectral lines with model atmospheres provides independent estimates of the surface latitudinal differential rotation and the inclination angle. Models of stars at different evolutionary stages rotating with internal conservative rotation laws were calculated to show that the Roche approximation can be safely used to account for the gravitational potential. The surface temperature gradient in rapid rotators induce an acceleration to the surface angular velocity. A non-zero differential rotation parameter may indicate that the rotation is neither rigid nor shellular underneath the stellar surface.

Context. Light elements and nitrogen surface abundances together can constrain the mixing efficiencies in massive stars on the main sequence, because moderate mixing at the surface leads to a depletion of light elements but only later to an enrichment in nitrogen. Aims. We want to test the rotational mixing prescriptions included in the Geneva stellar evolution code (GENEC) by following the evolution of surface abundances of light isotopes in massive stars. Methods. The GENEC is a 1D code containing sophisticated prescriptions for rotational mixing. We implemented an extended reaction network into this code including the light elements Li, Be and B, which allowed us to perform calculations testing the rotation induced mixing. Results. We followed 9, 12 and 15 solar mass models with rotation from the zero age main sequence up to the end of He burning. The calculations show the expected behaviour with faster depletion of light isotopes for faster rotating stars and more massive stars. Conclusions. We find that the mixing prescriptions used in the present rotating models for massive single stars can account for most of the observations. However the uncertainties are quite large making it hard to draw a firm conclusion on the mixing scenario.

At low metallicity the B-type stars rotate faster than at higher metallicity, typically in the SMC. As a consequence, it was expected a larger number of fast rotators in the SMC than in the Galaxy, in particular more Be/Oe stars. With the ESO-WFI in its slitless mode, the SMC open clusters were examined and an occurence of Be stars 3 to 5 times larger than in the Galaxy was found. The evolution of the angular rotational velocity seems to be the main key on the understanding of the specific behaviour and of the stellar evolution of such stars at different metallicities. With the results of this WFI study and using observational clues on the SMC WR stars and massive stars, as well as the theoretical indications of long gamma-ray burst progenitors, we identify the low metallicity massive Be and Oe stars as potential LGRB progenitors. Therefore the expected rates and numbers of LGRB are calculated and compared to the observed ones, leading to a good probability that low metallicity Be/Oe stars are actually LGRB progenitors.

Several studies have shown recently that at low metallicity B-type stars rotate faster than in environments of high metallicity. This is a typical case in the SMC. As a consequence, it is expected that a larger number of fast rotators is found in the SMC than in the Galaxy, in particular a higher fraction of Be/Oe stars. Using the ESO-WFI in its slitless mode, the data from the SMC open clusters were examined and an occurrence of Be stars 3 to 5 times larger than in the Galaxy was found. The evolution of the angular rotational velocity at different metallicities seems to be the main key to understand the specific behavior and evolution of these stars. According to the results from this WFI study, the observational clues obtained from the SMC WR stars and massive stars, and the theoretical predictions of the characteristics must have the long gamma-ray burst progenitors, we have identified the low metallicity massive Be and Oe stars as potential LGRB progenitors. To this end, the ZAMS rotational velocities of the SMC Be/Oe stars were determined and compared to models. The expected rates and the numbers of LGRB were then calculated and compared to the observed ones. Thus, a high probability was found that low metallicity Be/Oe stars can be LGRB progenitors. In this document, we describe the different steps followed in these studies: determination of the number of Be/Oe stars at different metallicities, identification of the clues that lead to suppose the low metallicity Be/Oe stars as LGRB progenitors, comparison of models with observations.

(abridged) The three most iron poor stars presently known ([Fe/H] equal to -5.96, -5.4 and -4.75) are carbon-rich, they are called C-Rich Ultra-Metal Poor Stars (CRUMPS). The origin of their peculiar surface abundances is not understood. We propose a synthetic view of the different models so far proposed to explain the peculiar abundances observed at the surface of the CRUMP stars. We deduce some expected trends based on nucleosynthetic arguments and look for signatures allowing to discriminate among models. We discuss the conditions for having CRUMP stars which are He-rich, i.e. with a mass fraction of helium greater than 0.30 and up to 0.60. We discuss the chemical composition of stars made of interstellar medium mixed with wind material of very metal poor massive stars, with wind plus supernova ejecta and with material extracted from the envelope of early AGB stars. Rotating and non-rotating models are considered. CRUMP stars formed from wind material of massive stars mixed with small amounts of pristine interstellar medium are He-rich (helium mass fraction between 0.30 and 0.60), Li-depleted and present low $^{12}$C/$^{13}$C ratios (inferior to 10 in number). Such He-rich stars, if discovered, would confirm that the most metal poor CRUMPs formed from essentially pure wind/envelope material. They would provide the most direct way to probe the nucleosynthetic outputs of the first generations of stars. We show that rotation is a key ingredient to explain the abundance patterns of CRUMPS stars.

Context: The identification of long-gamma-ray-bursts (LGRBs) is still uncertain, although the collapsar engine of fast-rotating massive stars is gaining a strong consensus. Aims: We propose that low-metallicity Be and Oe stars, which are massive fast rotators, as potential LGRBs progenitors. Methods: We checked this hypothesis by 1) testing the global specific angular momentum of Oe/Be stars in the ZAMS with the SMC metallicity, 2) comparing the ZAMS ($\Omega/\Omega_{\rm c},M/M_{\odot}$) parameters of these stars with the area predicted theoretically for progenitors with metallicity $Z=0.002$, and 3) calculating the expected rate of LGRBs/year/galaxy and comparing them with the observed ones. To this end, we determined the ZAMS linear and angular rotational velocities for SMC Be and Oe stars using the observed vsini parameters, corrected from the underestimation induced by the gravitational darkening effect. Results: The angular velocities of SMC Oe/Be stars are on average $<\Omega/\Omega_{\rm c}>=0.95$ in the ZAMS. These velocities are in the area theoretically predicted for the LGRBs progenitors. We estimated the yearly rate per galaxy of LGRBs and the number of LGRBs produced in the local Universe up to z=0.2. We have considered that the mass range of LGRB progenitors corresponds to stars hotter than spectral types B0-B1 and used individual beaming angles from 5 to 15\degr. We thus obtain $R^{\rm pred}_{\rm LGRB}\sim10^{-7}$ to $\sim10^{-6}$ LGRBs/year/galaxy, which represents on average 2 to 14 LGRB predicted events in the local Universe during the past 11 years. The predicted rates could widely surpass the observed ones [(0.2-3)$\times10^{-7}$ LGRBs/year/galaxy; 8 LGRBs observed in the local Universe during the last 11 years] if the stellar counts were made from the spectral type B1-B2, in accordance with the expected apparent spectral types of the appropriate massive fast rotators. Conclusion: We conclude that the massive Be/Oe stars with SMC metallicity could be LGRBs progenitors. Nevertheless, other SMC O/B stars without emission lines, which have high enough specific angular momentum, can enhance the predicted $R_{\rm LGRB}$ rate.

We show that the peculiar surface abundance patterns of Carbon Enhanced Metal Poor (CEMP) stars has been inherited from material having been processed by H- and He-burning phases in a previous generation of stars (hereafter called the ``Source Stars''). In this previous generation, some mixing must have occurred between the He- and the H-burning regions in order to explain the high observed abundances of nitrogen. In addition, it is necessary to postulate that a very small fraction of the carbon-oxygen core has been expelled (either by winds or by the supernova explosion). Therefore only the most outer layers should have been released by the Source Stars. Some of the CEMP stars may be He-rich if the matter from the Source Star is not too much diluted with the InterStellar Medium (ISM). Those stars formed from nearly pure ejecta would also be Li-poor.

A variation of the fundamental constants is expected to affect the thermonuclear rates important for stellar nucleosynthesis. In particular, because of the very small resonant energies of Be8 and C12, the triple $\alpha$ process is extremely sensitive to any such variations. Using a microscopic model for these nuclei, we derive the sensitivity of the Hoyle state to the nucleon-nucleon potential allowing for a change in the magnitude of the nuclear interaction. We follow the evolution of 15 and 60 solar mass, zero metallicity stellar models, up to the end of core helium burning. These stars are assumed to be representative of the first, Population III stars. We derive limits on the variation of the magnitude of the nuclear interaction and model dependent limits on the variation of the fine structure constant based on the calculated oxygen and carbon abundances resulting from helium burning. The requirement that some C12 and O16 be present are the end of the helium burning phase allows for permille limits on the change of the nuclear interaction and limits of order 10^{-5} on the fine structure constant relevant at a cosmological redshift of z ~ 15-20.

We discuss three effects of axial rotation at low metallicity. The first one is the mixing of the chemical species which is predicted to be more efficient in low metallicity environments. A consequence is the production of important quantities of primary $^{14}$N, $^{13}$C, $^{22}$Ne and a strong impact on the nucleosynthesis of the {\it s}-process elements. The second effect is a consequence of the first. Strong mixing makes possible the apparition at the surface of important quantities of CNO elements. This increases the opacity of the outer layers and may trigger important mass loss by line driven winds. The third effect is the fact that, during the main-sequence phase, stars, at very low metallicity, reach more easily than their more metal rich counterparts, the critical velocity\footnote{The critical velocity is the surface equatorial velocity such that the centrifugal acceleration compensates for the local gravity.}. We discuss the respective importance of these three effects as a function of the metallicity. We show the consequences for the early chemical evolution of the galactic halo and for explaining the CEMP stars. We conclude that rotation is probably a key feature which contributes in an important way to shape the evolution of the first stellar generations in the Universe.

When the surface angular velocity is above about 70% of the critical angular velocity, many interesting features appear which may be tested by interferometric observations, like significant deformation of stars, variation of the effective temperature with the latitude. Also polar winds become important and equatorial disks may appear. Near the critical limit, convection is also favored in the outer layers. In the present paper, we emphasize the need for a proper estimate of the critical velocity since this is the ratio of the actual velocity of the star to that critical velocity which determines the amplitude of the above effects. We recall the existence of two critical velocities. The first one, also called the classical critical velocity is the one to consider when the star has an Eddington factor inferior to 0.639, while the second one is the one to be considered when the Eddington factor is above 0.639. The features of the star at these two critical limits may be very different.

Massive stars are "cosmic engines" (cf the title of the IAU Symposium 250). They drive the photometric and chemical evolution of galaxies, inject energy and momentum through stellar winds and supernova explosions, they modify in this way the physical state of the interstellar gas and have an impact on star formation. The evolution of massive stars depends sensitively on the metallicity which has an impact on the intensity of the line driven stellar winds and on rotational mixing. We can distinguish four metallicity regimes: 1.- the Pop III regime $0 \le Z < \sim 10^{-10}$; 2.- The low metallicity regime $10^{-10} \le Z < 0.001$; 3.- The near solar metallicity regime $0.001 \le Z < 0.020$; 4.- The high metallicity regime $0.020 \le Z$. In each of these metallicity ranges, some specific physical processes occur. In this review we shall discuss these physical processes and their consequences for nucleosynthesis and the massive star populations in galaxies. We shall mainly focus on the effects of axial rotation and mass loss by line driven winds, although of course other processes like binarity, magnetic fields, transport processes by internal waves may also play important roles.

We review general characteristics of massive stars, present the main observable constraints that stellar models should reproduce. We discuss the impact of massive star nucleosynthesis on the early phases of the chemical evolution of the Milky Way (MW). We show that rotating models can account for the important primary nitrogen production needed at low metallicity. Interestingly such rotating models can also better account for other features as the variation with the metallicity of the C/O ratio. Damped Lyman Alpha (DLA) systems present similar characteristics as the halo of the MW for what concern the N/O and C/O ratios. Although in DLAs, the star formation history might be quite different from that of the halo, in these systems also, rotating stars (both massive and intermediate) probably play an important role for explaining these features. The production of primary nitrogen is accompanied by an overproduction of other elements as $^{13}$C, $^{22}$Ne and s-process elements. We show also how the observed variation with the metallicity of the number ratio of type Ibc to type II supernovae may be a consequence of the metallicity dependence of the line-driven stellar winds.

Recently, the concept of rotational mixing has been challenged by some authors (e.g. Hunter et al. 2008). We show that the excess N/H is a multivariate function f(M, age, v, sin i, multiplicity,Z). To find a correlation of a multivariate function with some parameter, it is evidently necessary to limit as much as possible the range of the other involved parameters. When this is done, the concept of rotational mixing is supported by the observations. We also show that the sample data are not free from several biases. A fraction of about 20% of the stars may escape to the relation as a result of binary evolution.

Mass loss and axial rotation are playing key roles in shaping the evolution of massive stars. They affect the tracks in the HR diagram, the lifetimes, the surface abundances, the hardness of the radiation field, the chemical yields, the presupernova status, the nature of the remnant, the mechanical energy released in the interstellar medium, etc... In this paper, after recalling a few characteristics of mass loss and rotation, we review the effects of these two processes at different metallicities. Rotation probably has its most important effects at low metallicities, while mass loss and rotation deeply affect the evolution of massive stars at solar and higher than solar metallicities.

We study the impact of the capture and annihilation of Weakly Interacting Massive Particles (WIMPs) on the evolution of Pop III stars. With a suitable modification of the Geneva stellar evolution code, we study the evolution of 20 and 200 M$_\odot$ stars in Dark Matter Haloes with densities between 10$^{8}$ and $10^{11}$ GeV/cm$^3$ during the core H-burning phase, and, for selected cases, until the end of the core He-burning phase. We find that for WIMP densities higher than 5.3 $10^{10}(\sigma^{SD}_p/10^{-38} {cm}^2)^{-1}$ GeV cm$^{-3}$ the core H-burning lifetime of $20 M_{\odot}$ and $200 M_{\odot}$ stars exceeds the age of the Universe, and stars are sustained only by WIMP annihilations. We determine the observational properties of these `frozen` objects and show that they can be searched for in the local Universe thanks to their anomalous mass-radius relation, which should allow unambiguous discrimination from normal stars.

New constraints on stellar models are provided by large surveys of massive stars, interferometric observations and asteroseismology. After a review of the main results so far obtained, we present new results from rotating models and discuss comparisons with observed features. We conclude that rotation is a key feature of massive star physics.

Massive stars played a key role in the early evolution of the Universe. They formed with the first halos and started the re-ionisation. It is therefore very important to understand their evolution. In this review, we first recall the effect of metallicity (Z) on the evolution of massive stars. We then describe the strong impact of rotation induced mixing and mass loss at very low Z. The strong mixing leads to a significant production of primary nitrogen 14, carbon 13 and neon 22. Mass loss during the red supergiant stage allows the production of Wolf-Rayet stars, type Ib,c supernovae and possibly gamma-ray bursts (GRBs) down to almost Z=0 for stars more massive than 60 solar masses. Galactic chemical evolution models calculated with models of rotating stars better reproduce the early evolution of N/O, C/O and C12/C13. Finally, the impact of magnetic fields is discussed in the context of GRBs.

Some useful developments in the model physics are briefly presented, followed by model results on chemical enrichments and WR stars. We discuss the expected rotation velocities of WR stars. We emphasize that the (C+O)/He ratio is a better chemical indicator of evolution for WC stars than the C/He ratios. With or without rotation, at a given luminosity the (C+O)/He ratios should be higher in regions of lower metallicity Z. Also, for a given (C+O)/He ratio the WC stars in lower Z regions have higher luminosities. The WO stars, which are likely the progenitors of supernovae SNIc and of some GRBs, should preferentially be found in regions of low Z and be the descendants of very high initial masses. Finally, we emphasize the physical reasons why massive rotating low Z stars may also experience heavy mass loss

The main effect of axial rotation on the evolution of massive PopIII stars is to trigger internal mixing processes which allow stars to produce significant amounts of primary nitrogen 14 and carbon 13. Very metal poor massive stars produce much more primary nitrogen than PopIII stars for a given initial mass and rotation velocity. The very metal poor stars undergo strong mass loss induced by rotation. One can distinguish two types of rotationnaly enhanced stellar winds: 1) Rotationally mechanical winds occurs when the surface velocity reaches the critical velocity at the equator, {\it i.e.} the velocity at which the centrifugal acceleration is equal to the gravity; 2) Rotationally radiatively line driven winds are a consequence of strong internal mixing which brings large amounts of CNO elements at the surface. This enhances the opacity and may trigger strong line driven winds. These effects are important for an initial value of $\upsilon/\upsilon_{\rm crit}$ of 0.54 for a 60 M$_\odot$ at $Z=10^{-8}$, {\it i.e.} for initial values of $\upsilon/\upsilon_{\rm crit}$ higher than the one ($\sim$0.4) corresponding to observations at solar $Z$. These two effects, strong internal mixing leading to the synthesis of large amounts of primary nitrogen and important mass losses induced by rotation, occur for $Z$ between about 10$^{-8}$ and 0.001. For metallicities above 0.001 and for reasonable choice of the rotation velocities, internal mixing is no longer efficient enough to trigger these effects.

Mass loss plays a dominant role in the evolution of massive stars at solar metallicity. After discussing different mass loss mechanisms and their metallicity dependence, we present the possibility of strong mass loss at very low metallicity. Our models at Z=1e-8 show that stars more massive than about 60 solar masses may lose a significant fraction of their initial mass in the red supergiant phase. This mass loss is due to the surface enrichment in CNO elements via rotational and convective mixing. Our 85 solar mass model ends its life as a fast rotating WO type Wolf-Rayet star. Therefore the models predict the existence of type Ic SNe and long and soft GRBs at very low metallicities. Such strong mass loss in the red supergiant phase or the Omega-Gamma limit could prevent the most massive stars from ending as pair-creation supernovae. The very low metallicity models calculated are also very interesting from the nucleosynthesis point of view. Indeed, the wind of the massive star models can reproduce the CNO abundances of the most metal-poor carbon-rich star known to date, HE1327-2326. Finally, using chemical evolution models, we are able to reproduce the evolution of CNO elements as observed in the normal extremely metal poor stars.

We investigate the effect of new stellar models, which take rotation into account, computed for very low metallicities on the chemical evolution of the earliest phases of the Milky Way. We check the impact of these new stellar yields on a model for the halo of the Milky Way that can reproduce the observed halo metallicity distribution. In this way we try to better constrain the ISM enrichment timescale, which was not done in our previous work. The stellar models adopted in this work 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 computed for large rotational velocities have a tremendous impact on the interstellar medium nitrogen enrichment for log(O/H)+12 < 7 (or [Fe/H]< -3). We show that upon the inclusion of the new stellar calculations in a chemical evolution model for the galactic halo with infall and outflow, 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 initial rotational velocities of the order of 600-800kms-1. An important contribution to N from AGB stars is still needed in order to explain the observations at intermediate metallicities. One possibility is that AGB stars at very low metallicities also rotate fast. This could be tested in the future, once stellar evolution models for fast rotating AGB stars will be available.

At very low metallicity, the effects of differential rotation have a more important impact on the evolution of stars than at high metallicity. Rotational mixing leads to the production of great quantities of helium and of primary $^{14}$N by massive stars. Rotation induces important mass loss and allows stars to locally strongly enrich the interstellar medium in CNO elements. Stars formed from interstellar clouds enriched by the winds of fast rotating massive stars would present surface abundances similar to those of C-rich extremely metal-poor stars. C-rich stars can also be formed by mass accretion in a binary system where the primary would be a fast rotating intermediate mass star in the early-AGB phase. Fast rotation may also lead to the formation of collapsars even at very low metallicity and make the most massive stars avoid the pair instability.

Rotating massive stars at $Z=10^{-8}$ and $10^{-5}$ lose a great part of their initial mass through stellar winds. The chemical composition of the rotationally enhanced winds of very low $Z$ stars is very peculiar. The winds show large CNO enhancements by factors of $10^3$ to $10^7$, together with large excesses of $^{13}$C and $^{17}$O and moderate amounts of Na and Al. The excesses of primary N are particularly striking. When these ejecta from the rotationally enhanced winds are diluted with the supernova ejecta from the corresponding CO cores, we find [C/Fe], [N/Fe],[O/Fe] abundance ratios very similar to those observed in the C--rich extremely metal poor stars (CEMP). We show that rotating AGB stars and rotating massive stars have about the same effects on the CNO enhancements. Abundances of s-process elements and the $^{12}$C/$^{13}$C ratio could help us to distinguish between contributions from AGB and massive stars. On the whole, we emphasize the dominant effects of rotation for the chemical yields of extremely metal poor stars.

We examine the role of rotation on the evolution and chemical yields of very metal--poor stars. The models include the same physics, which was applied successfully at the solar $Z$ and for the SMC, in particular, shear diffusion, meridional circulation, horizontal turbulence, and rotationally enhanced mass loss. Models of very low $Z$ experience a much stronger internal mixing in all phases than at solar $Z$. Also, rotating models at very low $Z$, contrary to the usual considerations, show a large mass loss, which mainly results from the efficient mixing of the products of the 3$\alpha$ reaction into the H--burning shell. This mixing allows convective dredge--up to enrich the stellar surface in heavy elements during the red supergiant phase, which in turn favours a large loss of mass by stellar winds, especially as rotation also increases the duration of this phase. On the whole, the low $Z$ stars may lose about half of their mass. Massive stars initially rotating at half of their critical velocity are likely to avoid the pair--instability supernova. The chemical composition of the rotationally enhanced winds of very low $Z$ stars show large CNO enhancements by factors of $10^3$ to $10^7$, together with large excesses of $^{13}$C and $^{17}$O and moderate amounts of Na and Al. The excesses of primary N are particularly striking. When these ejecta from the rotationally enhanced winds are diluted with the supernova ejecta from the corresponding CO cores, we find [C/Fe], [N/Fe],[O/Fe] abundance ratios that are very similar to those observed in the C--rich, extremely metal--poor stars (CEMP). We show that rotating AGB stars and rotating massive stars have about the same effects on the CNO enhancements.