We introduce SNAPSHOT, a technique to systematically compute stellar structure models in hydrostatic and thermal equilibrium based on 3 structural properties - core mass $M_{\rm core}$, envelope mass $M_{\rm env}$ and core composition. This approach allows us to connect these properties of stellar interiors to the luminosity and effective temperature $T_{\rm eff}$ in a more systematic way than with stellar evolution models. For MS models, we derive an analytical relationship between $M_{\rm core}$, $M_{\rm total}$ and central H abundance that can be used in rapid stellar evolution algorithms. Core-He burning models with $M_{\rm core}/M_{\rm total}$ from 0.2 to 0.8 have convective envelopes, low $T_{\rm eff}$ and will appear as red supergiants. For a given $M_{\rm core}$, they exhibit a small variation in luminosity (0.02 dex) and $T_{\rm eff}$ ($\sim 400\,\mathrm{K}$) over a wide range of $M_{\rm env}$ ($\sim 2 - 20\,\mathrm{M}_{\odot}$). This means that it is not possible to derive red supergiant masses from luminosities and $T_{\rm eff}$ alone. We derive the following relationship between $M_{\rm core}$ and the total luminosity of a red supergiant during core He burning: $\log M_{\rm core} \simeq 0.44 \log L/L_{\odot} - 1.38$. At $M_{\rm core}$/$M_{\rm total} \approx 0.2$, our models exhibit a bi-stability and jump from a RSG to a BSG structure. Our models with $M_{\rm core}/M_{\rm total} > 0.8$, which correspond to stripped stars produced by mass loss or binary interaction, show that $T_{\rm eff}$ has a strong dependence on $M_{\rm env}$, $M_{\rm core}$ and the core composition. We find the mass of one of these stripped stars in a binary system, HD 45166, to be less than its dynamical mass. When a large observational sample of stripped stars becomes available, our results can be used to constrain their masses and the physics of binary interaction.

Massive stars are the drivers of the chemical evolution of dwarf galaxies. We review here the basics of massive star evolution and the specificities of stellar evolution in low-Z environment. We discuss nucleosynthetic aspects and what observations could constrain our view on the first generations of stars.

The effects of rotation on stellar evolution are particularly important at low metallicity, when mass loss by stellar winds diminishes and the surface enrichment due to rotational mixing becomes relatively more pronounced than at high metallicities. Here we investigate the impact of rotation and metallicity on stellar evolution. Using a similar physics as in our previous large grids of models at Z=0.002 and Z=0.014, we compute stellar evolution models with the Geneva code for rotating and nonrotating stars with initial masses (Mini) between 1.7 and 120 Msun and Z=0.0004 (1/35 solar). This is comparable to the metallicities of the most metal poor galaxies observed so far, such as I Zw 18. Concerning massive stars, both rotating and nonrotating models spend most of their core-helium burning phase with an effective temperature higher than 8000 K. Stars become red supergiants only at the end of their lifetimes, and few RSGs are expected. Our models predict very few to no classical Wolf-Rayet stars as a results of weak stellar winds at low metallicity. The most massive stars end their lifetimes as luminous blue supergiants or luminous blue variables, a feature that is not predicted by models with higher metallicities. Interestingly, due to the behavior of the intermediate convective zone, the mass domain of stars producing pair-instability supernovae is smaller at Z=0.0004 than at Z=0.002. We find that during the main sequence phase, the ratio between nitrogen and carbon abundances (N/C) remains unchanged for nonrotating models. However, N/C increases by factors of 10-20 in rotating models at the end of the MS. Cepheids coming from stars with Mini > 4-6 Msun are beyond the core helium burning phase and spend little time in the instability strip. Since they would evolve towards cooler effective temperatures, these Cepheids should show an increase of the pulsation period as a function of age.

Grids of stellar models are useful tools to derive the properties of stellar clusters, in particular young clusters hosting massive stars, and to provide information on the star formation process in various mass ranges. Because of their short evolutionary timescale, massive stars end their life while their low-mass siblings are still on the pre-main sequence (pre-MS) phase. Thus the study of young clusters requires consistent consideration of all the phases of stellar evolution. But despite the large number of grids that are available in the literature, a grid accounting for the evolution from the pre-MS accretion phase to the post-MS phase in the whole stellar mass range is still lacking. We build a grid of stellar models at solar metallicity with masses from 0.8 $M_\odot$ to 120 $M_\odot$, including pre-MS phase with accretion. We use the {\sc genec} code to run stellar models on this mass range. The accretion law is chosen to match the observations of pre-MS objects on the Hertzsprung-Russell diagram. We describe the evolutionary tracks and isochrones of our models. The grid is connected to previous MS and post-MS grids computed with the same numerical method and physical assumptions, which provides the widest grid in mass and age to date. Numerical tables of our models and corresponding isochrones are available online.

Extended MS turn-offs are features commonly found in the colour-magnitude diagrams of young and intermediate age (less than about 2 Gyr) massive star clusters, where the MS turn-off is broader than can be explained by photometric uncertainties, crowding, or binarity. Rotation is suspected to be the cause of this feature, by accumulating fast rotating stars, strongly affected by gravity darkening and rotation-induced mixing, near the MS turn-off. This scenario successfully reproduces the tight relation between the age and the actual extent in luminosity of the extended MS turn-off of observed clusters. Below a given mass (dependent on the metallicity), stars are efficiently braked early on the MS due to the interaction of stellar winds and the surface magnetic field, making their tracks converge towards those of non-rotating tracks in the HRD. When these stars are located at the turn-off of a cluster, their slow rotation causes the extended MS turn-off feature to disappear. We investigate the maximal mass for which this braking occurs at different metallicities, and determine the age above which no extended MS turn-off is expected in clusters. Our models predict that the extended MS turn-off phenomenon disappears at ages older than about 2 Gyr. There is a trend with the metallicity, the age at which the disappearance occurs becoming older at higher metallicity. These results are robust between the two codes used in this work, despite some differences in the input physics and in particular in the detailed description of rotation-induced internal processes and of angular momentum extraction by stellar winds. Comparing our results with clusters in the LMC and Galaxy shows a very good fit to the observations. This strengthens the rotation scenario to explain the cause of the extended MS turn-off phenomenon.

Context: The observations of solar-like oscillations in evolved stars have brought important constraints on their internal rotation rates. To correctly reproduce these data, an efficient transport mechanism is needed in addition to meridional circulation and shear instability. Aims: We study the efficiency of the transport of angular momentum during the subgiant phase. Results: The precise asteroseismic measurements of both core and surface rotation rates available for the six Kepler targets enable a precise determination of the efficiency of the transport of angular momentum needed for each of these subgiants. These results are found to be insensitive to all the uncertainties related to the modelling of rotational effects before the post-main sequence phase. An interesting exception in this context is the case of young subgiants (typical values of log(g) close to 4), because their rotational properties are sensitive to the degree of radial differential rotation on the main sequence. These young subgiants constitute therefore perfect targets to constrain the transport of angular momentum on the main sequence from asteroseismic observations of evolved stars. As for red giants, we find that the efficiency of the additional transport process increases with the mass of the star during the subgiant phase. However, the efficiency of this undetermined mechanism decreases with evolution during the subgiant phase, contrary to what is found for red giants. Consequently, a transport process with an efficiency that increases with the degree of radial differential rotation cannot account for the core rotation rates of subgiants, while it correctly reproduces the rotation rates of red giant stars. This suggests that the physical nature of the additional mechanism needed for the internal transport of angular momentum may be different in subgiant and red giant stars.

A large fraction of massive stars evolve in interacting binary systems, which dramatically modifies the outcome of stellar evolution. We investigated the properties of blue supergiants in binary systems and whether they are suitable for extragalactic distance determinations using the flux-weighted gravity luminosity relationship (FGLR). This is a relationship between the absolute bolometric magnitude $M_{\rm bol}$ and the spectroscopically determined flux-weighted gravity $g_{\rm F} = g/T^4_{\rm eff}$. We computed a grid of binary stellar evolution models with MESA and use the v2.1 BPASS models to examine whether they are compatible with the relatively small scatter shown by the observed relationship. Our models have initial primary masses of 9 - 30 $M_\odot$, initial orbital periods of 10 - 2511 days, mass ratio $q$ = 0.9, and metallicity $Z$ = 0.02. We find that the majority of primary stars that produce blue supergiant stages are consistent with the observed FGLR, with a small offset towards brighter bolometric magnitudes. In between 1\% - 24\% of cases, binary evolution may produce blue supergiants after a mass transfer episode, that lie below the observed FGLR. A small number of such stars have been found in extragalactic FGLR studies, suggesting that they may have evolved through binary interaction. Some models with shorter periods could resemble blue hypergiants and LBVs. We used CMFGEN radiative transfer models to investigate the effects of unresolved secondaries on diagnostics for $T_{\rm eff}$ and $g$, and the biases on the determination of interstellar reddening and $M_{\rm bol}$. We find that the effects are small and within the observed scatter, but could lead to a small overestimate of the luminosity, of $T_{\rm eff}$ and of $g$ for extreme cases. We conclude that the observed FGLR can, in principle, be well reproduced by close binary evolution models.

Binary stars are of course more than two stars, but they are also at least two stars. In this chapter we will review some aspects of the physics governing the evolution of single massive stars. We will also review the uncertainties of key physical ingredients: mass loss, rotation and convection.

We present an analysis of the relatively low mass ($\sim2400$~M$_{\odot}$), $\sim800$~Myr, Galactic open cluster, NGC~2818, using Gaia DR2 results combined with VLT/FLAMES spectroscopy. Using Gaia DR2 proper motions and parallax measurements we are able to select a clean sample of cluster members. This cluster displays a clear extended main sequence turn-off (eMSTO), a phenomenon previously studied mainly in young and intermediate age massive clusters in the Magellanic clouds. The main sequence of NGC~2818 is extremely narrow, with a width of $\sim0.01$ magnitudes (G$_{\rm BP} - $ G$_{\rm RP}$), suggesting very low levels of differential extinction. Using VLT/FLAMES spectroscopy of 60 cluster members to measure the rotational velocity of the stars (Vsini) we find that stars on the red side of the eMSTO have high Vsini ($>160$~km/s) while stars on the blue side have low Vsini ($<160$~km/s), in agreement with model predictions. The cluster also follows the previously discovered trend between the age of the cluster and the extent of the eMSTO. We conclude that stellar rotation is the likely cause of the eMSTO phenomenon.

When planets are formed from the protoplanetary disk and after the disk has dissipated, the evolution of their orbits is governed by tidal interactions, friction, and gravitational drag, and also by changes in the mass of the star and planet. These interactions may change the initial distribution of the distances between the planets and their host star by expanding the original orbit, by contracting it (which may cause an engulfment of the planet by the star), or by destroying the planet. We study the evolution of the orbit of a planet orbiting its host star under the effects of equilibrium tides, dynamical tides, drag (frictional and gravitational), and stellar mass loss. We used the Geneva stellar evolution code to compute the evolution of stars with initial masses of 1 and 1.5 solar mass with different rotation rates at solar metallicity. The star is evolved from the pre-main-sequence (PMS) up to the tip of the red giant branch. We used these models as input for computing the evolution of the planetary orbits. We explored the effects of changing the planet masses (of 1 Earth mass up to 20 Jupiter masses), the distance between the planet and the star (of 0.015 and more than 3 au), the mass, and the spin of the star. We present results when only the equilibrium tide was accounted for and when both equilibrium and dynamical tides were accounted for. The expression for the dynamical tide is a frequency-averaged dissipation of tidally excited inertial waves, obtained from a piecewise homogeneous two-layer stellar model. Gravity wave damping was neglected. Dynamical tides in convective zones have a significant effect on planetary orbits only during the PMS phase and only for fast-rotating stars. They have no significant effects during the PMS phase for initially slow-rotating stars and during the red giant branch phase, regardless of the initial rotation (abridged).

Through the modelling of near-infrared photometry of star plus disk systems with the codes BEDISK/BERAY, we successfully describe the WISE photometric characteristics of Be stars in five young open clusters, NGC 663, NGC 869, NGC 884, NGC 3766 and NGC 4755, broadly studied in the literature. WISE photometry allows previously known Be stars to be detected and to find new Be candidates which could be confirmed spectroscopically. The location of Be stars in the WISE colour-magnitude diagram, separates them in two groups; active (Be stars hosting a developed circumstellar disk) and quiescent objects (Be stars in a diskless phase), and this way, we can explore how often stars are observed in these different stages. The variability observed in most active variable Be stars is compatible with a disk dissipation phase. We find that 50 percent of Be stars in the studied open clusters are in an active phase. We can interpret this as Be stars having a developed circumstellar disk half of the time. The location of Be stars with a developed disk in the CMD require mass loss rates in agreement with values recently reported in the literature. For these objects, we expect to have a tight relation between the equivalent width of the H{\alpha} line and the mass of the disk, if the inclination is known. Also, near-infrared photometry of Be stars in stellar clusters has the potential of being useful to test whether there is a preferential viewing angle.

CEMP-s stars are long-lived low-mass stars with a very low iron content as well as overabundances of carbon and s-elements. Their peculiar chemical pattern is often explained by pollution from a AGB star companion. Recent observations have shown that most of the CEMP-s stars are in a binary system, providing support to the AGB companion scenario. A few CEMP-s stars, however, appear to be single. We inspect four apparently single CEMP-s stars and discuss the possibility that they formed from the ejecta of a previous-generation massive star, referred to as the "source" star. In order to investigate this scenario, we computed low-metallicity massive star models with and without rotation and including complete s-process nucleosynthesis. We find that non-rotating source stars cannot explain the observed abundance of any of the four CEMP-s stars. Three out of the four CEMP-s stars can be explained by a $25$ $M_{\odot}$ source star with $v_{\rm ini} \sim 500$ km s$^{-1}$ (spinstar). The fourth CEMP-s star has a high Pb abundance that cannot be explained by any of the models we computed. Since spinstars and AGB predict different ranges of [O/Fe] and [ls/hs], these ratios could be an interesting way to further test these two scenarios.

Context. CEMP-no stars are long-lived low-mass stars with a very low iron content, overabundances of carbon and no or minor signs for the presence of s- or r-elements. Although their origin is still a matter of debate, they are often considered as being made of a material ejected by a previous stellar generation (source stars). Aims. We place constraints on the source stars from the observed abundance data of CEMP-no stars. Methods. We computed source star models of 20, 32, and 60 M$_{\odot}$ at Z = 10$^{-5}$ with and without fast rotation. For each model we also computed a case with a late mixing event occurring between the hydrogen and helium-burning shell $\sim$ 200 yr before the end of the evolution. This creates a partially CNO-processed zone in the source star. We use the 12C/13C and C/N ratios observed on CEMP-no stars to put constraints on the possible source stars (mass, late mixing or not). Then, we inspect more closely the abundance data of six CEMP-no stars and select their preferred source star(s). Results. Four out of the six CEMP-no stars studied cannot be explained without the late mixing process in the source star. Two of them show nucleosynthetic signatures of a progressive mixing (due e.g. to rotation) in the source star. We also show that a 20 M$_{\odot}$ source star is preferred compared to one of 60 M$_{\odot}$ and that likely only the outer layers of the source stars were expelled to reproduce the observed 12C/13C. Conclusions. The results suggest that (1) a late mixing process could operate in some source stars, (2) a progressive mixing, possibly achieved by fast rotation, is at work in several source stars, (3) $\sim$ 20 M$_{\odot}$ source stars are preferred compared to $\sim$ 60 M$_{\odot}$ ones, and (4) the source star might have preferentially experienced a low energetic supernova with large fallback.

We give a brief overview of where we stand with respect to some old and new questions bearing on how massive stars evolve and end their lifetime. We focus on the following key points that are further discussed by other contributions during this conference: convection, mass losses, rotation, magnetic field and multiplicity. For purpose of clarity, each of these processes are discussed on its own but we have to keep in mind that they are all interacting between them offering a large variety of outputs, some of them still to be discovered.

Classical Cepheid variable stars are both sensitive astrophysical laboratories and accurate cosmic distance tracers. We have recently investigated how the evolutionary effects of rotation impact the properties of these important stars and here provide an accessible overview of some key elements as well as two important consequences. Firstly, rotation resolves the long-standing Cepheid mass discrepancy problem. Second, rotation increases main sequence lifetimes, i.e, Cepheids are approximately twice as old as previously thought. Finally, we highlight the importance of the short-period ends of Cepheid period distributions as indicators for model adequacy.

We study the pulsational properties of rapidly rotating main-sequence B-type stars using linear non-adiabatic analysis of non-radial low-frequency modes taking into account the effect of rotation. We compare the properties of prograde sectoral $g$ and retrograde $r$ modes excited by the $\kappa$ mechanism at the Fe opacity peak with the newly discovered period-luminosity relation that is obeyed by a group of fast-rotating B-type stars in the young open cluster NGC 3766. The observed relation consists of two sequences in the period versus magnitude diagram, at periods shorter than 0.5 days. We find that this property is consistent with similar period-luminosity relations predicted for excited sectoral prograde $g$-modes of azimuthal orders $m=-1$ and $m=-2$ in fast-rotating stars along an isochrone. We further show that some of the rapidly rotating stars that have photometric variability with periods longer than a day may be caused by $r$-mode pulsation predicted to be excited in these stars. One fast-rotating star, in particular, shows both short and long periods that can be explained by the simultaneous excitation of $g$- and $r$-mode pulsations in models of fast-rotating stars.

Near-solar metallicity (and low-redshift) Pair-Instability Supernova (PISN) candidates challenge stellar evolution models. Indeed, at such a metallicity, even an initially very massive star generally loses so much mass by stellar winds that it will avoid the electron-positron pair-creation instability. We use recent results showing that a magnetic field at the surface of a massive star can significantly reduce its effective mass-loss rate to compute magnetic models of very massive stars (VMSs) at solar metallicity and explore the possibility that such stars end as PISNe. We implement the quenching of the mass loss produced by a surface dipolar magnetic field into the Geneva stellar evolution code and compute new stellar models with an initial mass of $200\,M_\odot$ at solar metallicity, with and without rotation. It considerably reduces the total amount of mass lost by the star during its life. For the non-rotating model, the total (CO-core) mass of the models is $72.8\,M_\odot$ ($70.1\,M_\odot$) at the onset of the electron-positron pair-creation instability. For the rotating model, we obtain $65.6\,M_\odot$ ($62.4\,M\odot$). In both cases, a significant fraction of the internal mass lies in the region where pair instability occurs in the $\log(T)-\log(\rho)$ plane. The interaction of the reduced mass loss with the magnetic field efficiently brakes the surface of the rotating model, producing a strong shear and hence a very efficient mixing that makes the star evolve nearly homogeneously. The core characteristics of our models indicate that solar metallicity models of magnetic VMSs may evolve to PISNe (and pulsation PISNe).

The modelling of massive star evolution is a complex task, and is very sensitive to the way physical processes (such as convection, rotation, mass loss, etc.) are included in stellar evolution code. Moreover, the very high observed fraction of binary systems among massive stars makes the comparison with observations difficult. In this paper, we focus on discussing the uncertainties linked to the modelling of convection and rotation in single massive stars.

Context: Constraints on the internal rotation of red giants are now available thanks to asteroseismic observations. Preliminary comparisons with rotating stellar models indicate that an undetermined additional process for the internal transport of angular momentum is required in addition to purely hydrodynamic processes. Aims: We investigate how asteroseismic measurements of red giants can help us characterize the additional transport mechanism. Methods: We first determine the efficiency of the missing transport mechanism for the low-mass red giant KIC 7341231 by computing rotating models that include an additional viscosity corresponding to this process. We then discuss the change in the efficiency of this transport of angular momentum with the mass, metallicity and evolutionary stage. Results: In the case of the low-mass red giant KIC 7341231, we find that the viscosity corresponding to the additional mechanism is constrained to the range 1 x 10^3 - 1.3 x 10^4 cm^2/s. This constraint on the efficiency of the unknown additional transport mechanism during the post-main sequence is obtained independently of any specific assumption about the modelling of rotational effects during the pre-main sequence and the main sequence (in particular, the braking of the surface by magnetized winds and the efficiency of the internal transport of angular momentum before the post-main-sequence phase). When we assume that the additional transport mechanism is at work during the whole evolution of the star together with a solar-calibrated braking of the surface by magnetized winds, the range of nu_add is reduced to 1 - 4 x 10^3 cm^2/s. In addition to being sensitive to the evolutionary stage of the star, we show that the efficiency of the unknown process for internal transport of angular momentum increases with the stellar mass.

In this paper, we discuss some consequences of rotation and mass loss on the evolved stages of massive star evolution. The physical reasons of the time evolution of the surface velocity are explained, and then we show how the late-time evolution of massive stars are impacted in combination with the effects of mass loss. The most interesting result is that in some cases, a massive star can have a blue-red-blue evolution, opening the possibility that Blue Supergiants are composed by two distinct populations of stars: one just leaving the main sequence and crossing the HRD for the first time, and the other one evolving back to the blue side of the HRD after a Red Supergiant phase. We discuss a few possible observational tests that can allow to distinguish these two populations, and how supergiant B[e] stars fit in this context.

$Context.$ Pulsating stars are windows to the physics of stars enabling us to see glimpses of their interior. Not all stars pulsate, however. On the main sequence, pulsating stars form an almost continuous sequence in brightness, except for a magnitude range between $\delta$ Scuti and slowly pulsating B stars. Against all expectations, 36 periodic variables were discovered in 2013 in this luminosity range in the open cluster NGC 3766, the origins of which was a mystery. $Aims.$ We investigate the properties of those new variability class candidates in relation to their stellar rotation rates and stellar multiplicity. $Methods.$ We took multi-epoch spectra over three consecutive nights using ESO's Very Large Telescope. $Results.$ We find that the majority of the new variability class candidates are fast-rotating pulsators that obey a new period-luminosity relation. We argue that the new relation discovered here has a different physical origin to the period-luminosity relations observed for Cepheids. $Conclusions.$ We anticipate that our discovery will boost the relatively new field of stellar pulsation in fast-rotating stars, will open new doors for asteroseismology, and will potentially offer a new tool to estimate stellar ages or cosmic distances.

Context. Red-giant stars may engulf planets. This may increase the rotation rate of their convective envelope, which could lead to strong dynamo-triggered magnetic fields. Aims. We explore the possibility of generating magnetic fields in red giants that have gone through the process of a planet engulfment. We compare them with similar models that evolve without any planets. We discuss the impact of stellar wind magnetic braking on the evolution of the surface velocity of the parent star. Methods. With rotating stellar models with and without planets and an empirical relation between the Rossby number and the surface magnetic field, we deduce the evolution of the surface magnetic field along the red-giant branch. The effects of wind magnetic braking is explored using a relation deduced from MHD simulations. Results. The stellar evolution model of a 1.7 M$_\odot$ without planet engulfment and that has a time-averaged rotation velocity during the Main-Sequence equal to 100 km s$^{-1}$, shows a surface magnetic field triggered by convection larger than 10 G only at the base of the red giant branch, that means for gravities log $g > 3$. When a planet engulfment occurs, such magnetic field can also appear at much lower gravities, i.e. at much higher luminosities along the red giant branch. Typically the engulfment of a 15 M$_J$ planet produces a dynamo triggered magnetic field larger than 10 G for gravities between 2.5 and 1.9. We show that for reasonable wind magnetic braking laws, the high surface velocity reached after a planet engulfment may be maintained sufficiently long for being observable. Conclusions. High surface magnetic fields for red giants in the upper part of the red giant branch is a strong indication of a planet engulfment or of an interaction with a companion. Our theory can be tested by observing fast rotating red giants and check whether they show magnetic fields.

Classical Cepheid variable stars are high-sensitivity probes of stellar evolution and fundamental tracers of cosmic distances. While rotational mixing significantly affects the evolution of Cepheid progenitors (intermediate-mass stars), the impact of the resulting changes in stellar structure and composition on Cepheids on their pulsational properties is hitherto unknown. Here we present the first detailed pulsational instability analysis of stellar evolution models that include the effects of rotation, for both fundamental mode and first overtone pulsation. We employ Geneva evolution models spanning a three-dimensional grid in mass (1.7 - 15 $M_\odot$), metallicity (Z = 0.014, 0.006, 0.002), and rotation (non-rotating, average & fast rotation). We determine (1) hot and cool instability strip (IS) boundaries taking into account the coupling between convection and pulsation, (2) pulsation periods, and (3) rates of period change. We investigate relations between period and (a) luminosity, (b) age, (c) radius, (d) temperature, (e) rate of period change, (f) mass, (g) the flux-weighted gravity-luminosity relation (FWGLR). We confront all predictions aside from those for age with observations, finding generally excellent agreement. We tabulate period-luminosity relations (PLRs) for several photometric pass-bands and investigate how the finite IS width, different IS crossings, metallicity, and rotation affect PLRs. We show that a Wesenheit index based on H, V, and I photometry is expected to have the smallest intrinsic PLR dispersion. We confirm that rotation resolves the Cepheid mass discrepancy. Period-age relations depend significantly on rotation (rotation increases Cepheid ages), offering a straightforward explanation for evolved stars in binary systems that cannot be matched by conventional isochrones assuming a single age. Finally, we show that Cepheids obey a tight FWGLR.

Key physical ingredients governing the evolution of massive stars are mass losses, convection and mixing in radiative zones. These effects are important both in the frame of single and close binary evolution. The present paper addresses two points: 1) the differences between two families of rotating models, i.e. the family of models computed with and without an efficient transport of angular momentum in radiative zones; 2) The impact of the mass losses in single and in close binary models.

Rotation was shown to have a strong impact on the structure and light element nucleosynthesis in massive stars. In particular, models including rotation can reproduce the primary nitrogen observed in halo extremely metal-poor (EMP) stars. Additional exploratory models showed that rotation may enhance $s$-process production at low metallicity. Here we present a large grid of massive star models including rotation and a full $s$-process network to study the impact of rotation on the weak $s$-process. We explore the possibility of producing significant amounts of elements beyond the strontium peak, which is where the weak $s$-process usually stops. We used the Geneva stellar evolution code coupled to an enlarged reaction network with 737 nuclear species up to bismuth to calculate $15-40\,\text{M}_\odot$ models at four metallicities ($Z = 0.014,10^{-3}$, $10^{-5}$, and $10^{-7}$) from the main sequence up to the end of oxygen burning. We confirm that rotation-induced mixing between the convective H-shell and He-core enables an important production of primary $^{14}$N and $^{22}$Ne and $s$-process at low metallicity. At low metallicity, even though the production is still limited by the initial number of iron seeds, rotation enhances the $s$-process production, even for isotopes heavier than strontium, by increasing the neutron to seed ratio. The increase in this ratio is a direct consequence of the primary production of $^{22}$Ne. Despite nuclear uncertainties affecting the $s$-process production and stellar uncertainties affecting the rotation-induced mixing, our results show a robust production of $s$ process at low metallicity when rotation is taken into account. Considering models with a distribution of initial rotation rates enables to reproduce the observed large range of the [Sr/Ba] ratios in (carbon-enhanced and normal) EMP stars.

During the last few years, the Geneva stellar evolution group has released new grids of stellar models, including the effect of rotation and with updated physical inputs (Ekstr\"om et al. 2012; Georgy et al. 2013a,b). To ease the comparison between the outputs of the stellar evolution computations and the observations, a dedicated tool was developed: the Syclist toolbox (Georgy et al. 2014). It allows to compute interpolated stellar models, isochrones, synthetic clusters, and to simulate the time-evolution of stellar populations.

Mass-loss rates during the red supergiant phase are very poorly constrained from an observational or theoretical point of view. However, they can be very high, and make a massive star lose a lot of mass during this phase, influencing considerably the final evolution of the star: will it end as a red supergiant? Will it evolve bluewards by removing its hydrogen-rich envelope? In this paper, we briefly summarise the effects of this mass loss and of the related uncertainties, particularly on the population of blue supergiant stars.

Wolf-Rayet (WR) stars, as they are advanced stages of the life of massive stars, provide a good test for various physical processes involved in the modelling of massive stars, such as rotation and mass loss. In this paper, we show the outputs of the latest grids of single massive stars computed with the Geneva stellar evolution code, and compare them with some observations. We present a short discussion on the shortcomings of single stars models and we also briefly discuss the impact of binarity on the WR populations.

We use the Geneva Syclist isochrone models that include the effects of stellar rotation to investigate the role that rotation has on the resulting colour-magnitude diagram (CMD) of young and intermediate age clusters. We find that if a distribution of rotation velocities exists within the clusters, rotating stars will remain on the main sequence (MS) for longer, appearing to be younger than non-rotating stars within the same cluster. This results in an extended main sequence turn-off (eMSTO) that appears at young ages ($\sim30$~Myr) and lasts beyond 1~Gyr. If this eMSTO is interpreted as an age spread, the resulting age spread is proportional to the age of the cluster, i.e. young clusters ($<100$~Myr) appear to have small age spreads (10s of Myr) whereas older clusters ($\sim1$~Gyr) appear to have much larger spreads, up to a few hundred Myr. We compare the predicted spreads for a sample of rotation rates to observations of young and intermediate age clusters, and find a strong correlation between the measured 'age spread' and the age of the cluster, in good agreement with models of stellar rotation. This suggests that the 'age spreads' reported in the literature may simply be the result of a distribution of stellar rotation velocities within clusters.

The post main-sequence evolution of massive stars is very sensitive to many parameters of the stellar models. Key parameters are the mixing processes, the metallicity, the mass-loss rate and the effect of a close companion. We study how the red supergiant lifetimes, the tracks in the Hertzsprung-Russel diagram (HRD), the positions in this diagram of the pre-supernova progenitor as well as the structure of the stars at that time change for various mass-loss rates during the red supergiant phase (RSG), and for two different initial rotation velocities. The surface abundances of RSGs are much more sensitive to rotation than to the mass-loss rates during that phase. A change of the RSG mass-loss rate has a strong impact on the RSG lifetimes and therefore on the luminosity function of RSGs. At solar metallicity, the enhanced mass-loss rate models do produce significant changes on the populations of blue, yellow and red supergiants. When extended blue loops or blue ward excursions are produced by enhanced mass-loss, the models predict that a majority of blue (yellow) supergiants are post RSG objects. These post RSG stars are predicted to show much smaller surface rotational velocities than similar blue supergiants on their first crossing of the HR gap. The position in the HRD of the end point of the evolution depends on the mass of the hydrogen envelope. More precisely, whenever, at the pre-supernova stage, the H-rich envelope contains more than about 5\% of the initial mass, the star is a red supergiant, and whenever the H-rich envelope contains less than 1\% of the total mass the star is a blue supergiant. For intermediate situations, intermediate colors/effective temperatures are obtained. Yellow progenitors for core collapse supernovae can be explained by the enhanced mass-loss rate models, while the red progenitors are better fitted by the standard mass-loss rate models.

We recently showed that rotation significantly affects most observable Cepheid quantities, and that rotation, in combination with the evolutionary status of the star, can resolve the long-standing Cepheid mass discrepancy problem. We therefore provide a brief overview of our results regarding the problem of Cepheid masses. We also briefly mention the impact of rotation on the Cepheid period-luminosity(-color) relation, which is crucial for determining extragalactic distances, and thus for calibrating the Hubble constant.

(abridged) We present a new tool for building synthetic colour-magnitude diagrams of coeval stellar populations. We study, from a theoretical point of view, the impact of axial rotation of stars on various observed properties of single-aged stellar populations: magnitude at the turnoff, photometric properties of evolved stars, surface velocities, surface abundances, and the impact of rotation on the age determination of clusters by an isochrone fitting. Stellar models for different initial masses, metallicities, and zero-age main sequence (ZAMS) rotational velocities are used for building interpolated stellar tracks, isochrones, and synthetic clusters for various ages and metallicities. The synthetic populations account for the effects of the initial distribution of the rotational velocities on the ZAMS, the impact of the inclination angle and the effects of gravity and limb darkening, unresolved binaries and photometric errors. Interpolated tracks, isochrones, and synthetic clusters can be computed through a public web interface. For clusters with a metallicity in the range [0.002, 0.014] and an age between 30 Myr and 1 Gyr, the fraction of fast rotators on the main sequence (MS) band is the largest just below the turnoff. This is a natural consequence of the increase in the MS lifetime due to rotation. The fraction of fast rotators one magnitude below the turnoff also increases with the age of the cluster between 30 Myr and 1 Gyr. There is an increase in the fraction of enriched stars when the metallicity decreases. We show that the use of isochrones computed from rotating stellar models with an initial rotation that is representative of the average initial rotation of the stars in clusters provides a reasonable estimate of the age, even though stars in a real cluster did not start their evolution with an identical initial rotation.

The nitrogen to carbon (N/C) and nitrogen to oxygen (N/O) ratios are the most sensitive quantities to mixing in stellar interiors of intermediate and massive stars. We further investigate the theoretical properties of these ratios as well as put in context recent observational results obtained by the VLT-FLAMES Survey of massive stars in the Galaxy and the Magellanic Clouds. We consider analytical relations and numerical models of stellar evolution as well as our own stellar atmosphere models, and we critically re-investigate observed spectra. On the theoretical side, the N/C vs N/O plot shows little dependence on the initial stellar masses, rotation velocities, and nature of the mixing processes up to relative enrichment of N/O by a factor of about four, thus this plot constitutes an ideal quality test for observational results. The comparison between the FLAMES Survey and theoretical values shows overall agreement, despite the observational scatter of the published results. The existence of some mixing of CNO products is clearly confirmed, however the accuracy of the data is not sufficient for allowing a test of the significant differences between different models of rotating stars and the Geneva models. We discuss reasons (for the most part due to observational bias) why part of the observational data points should not be considered for this comparison. When these observational data points are not considered, the scatter is reduced. Finally, the N/C vs N/O plot potentially offers a powerful way for discriminating blue supergiants before the red supergiant stage from those after it. Also, red supergiants of similar low velocities may exhibit different N enrichments, depending on their initial rotation during the main-sequence phase.

[Abridged] We aim to improve the understanding of Cepheids from an evolutionary perspective and establish the role of rotation in the Cepheid paradigm. In particular, we are interested in the contribution of rotation to the problem of Cepheid masses, and explore testable predictions of quantities that can be confronted with observations. Evolutionary models including a homogeneous and self-consistent treatment of rotation are studied in detail during the crossings of the classical instability strip (IS). The dependence of several parameters on initial rotation is studied. These parameters include mass, luminosity, temperature, lifetimes, equatorial velocity, surface abundances, and rates of period change. Several key results are obtained: i) mass-luminosity (M-L) relations depend on rotation, particularly during the blue loop phase; ii) luminosity increases between crossings of the IS. Hence, Cepheid M-L relations at fixed initial rotation rate depend on crossing number (faster rotation yields greater luminosity difference between crossings); iii) the Cepheid mass discrepancy problem vanishes when rotation and crossing number are taken into account, without a need for high core overshooting values or enhanced mass loss; iv) rotation creates dispersion around average parameters predicted at fixed mass and metallicity. This is of particular importance for the period-luminosity-relation, for which rotation is a source of intrinsic dispersion; v) enhanced surface abundances do not unambiguously distinguish Cepheids occupying the Hertzsprung gap from ones on blue loops (after dredge-up), since rotational mixing can lead to significantly enhanced Main Sequence (MS) abundances; vi) rotating models predict greater Cepheid ages than non-rotating models due to longer MS lifetimes. Rotation has a significant evolutionary impact on classical Cepheids and should no longer be neglected in their study.

The force that governs the evolution of stars is gravity. Indeed this force drives star formation, imposes thermal and density gradients into stars at hydrostatic equilibrium and finally plays the key role in the last phases of their evolution. Nuclear power in stars governs their lifetimes and of course the stellar nucleosynthesis. The nuclear reactions are at the heart of the changes of composition of the baryonic matter in the Universe. This change of composition, in its turn, has profound consequences on the evolution of stars and galaxies. The energy extracted from the gravitational, respectively nuclear reservoirs during the lifetimes of stars of different masses are estimated. It is shown that low and intermediate mass stars (M < 8 Msol) extract roughly 90 times more energy from their nuclear reservoir than from their gravitational one, while massive stars (M > 8 Msol), which explode in a supernova explosion, extract more than 5 times more energy from the gravitational reservoir than from the nuclear one. We conclude by discussing a few important nuclear reactions and their link to topical astrophysical questions.

(shortened) We provide a grid of single star models covering a mass range from 0.8 to 120 Msun 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). 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/Lsun) > 5.5, is very sensitive to rotational mixing. Strong mixing significantly reduces the MS width. We confirm, but here for the first time on the whole mass range, 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/Lsun) 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 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.

We study how tides in a binary system induce some specific internal shear mixing, able to substantially modify the evolution of close binaries prior to mass transfer. We construct numerical models accounting for tidal interactions, meridional circulation, transport of angular momentum, shears and horizontal turbulence and consider a variety of orbital periods and initial rotation velocities. Depending on orbital periods and rotation velocities, tidal effects may spin down (spin down Case) or spin up (spin up Case) the axial rotation. In both cases, tides may induce a large internal differential rotation. The resulting tidally induced shear mixing (TISM) is so efficient that the internal distributions of angular velocity and chemical elements are greatly influenced. The evolutionary tracks are modified, and in both cases of spin down and spin up, large amounts of nitrogen can be transported to the stellar surfaces before any binary mass transfer. Meridional circulation, when properly treated as an advection, always tends to counteract the tidal interaction, tending to spin up the surface when it is braked down and vice versa. As a consequence, the times needed for the axial angular velocity to become equal to the orbital angular velocity may be larger than given by typical synchronization timescales. Also, due to meridional circulation some differential rotation remains in tidally locked binary systems.

B-type stars are known to rotate at various velocities, including very fast rotators near the critical velocity as the Be stars. In this paper, we provide stellar models covering the mass range between 1.7 to 15 Msun, which includes the typical mass of known Be stars, at Z = 0.014, 0.006, and 0.002 and for an extended range of initial velocities on the zero-age main sequence. We used the Geneva stellar-evolution code, including the effects of shellular rotation, with a numerical treatment that has been improved so the code can precisely track the variation in the angular momentum content of the star as it changes under the influence of radiative winds and/or mechanical mass loss. We discuss the impact of the initial rotation rate on the tracks in the Hertzsprung-Russell diagram, the main-sequence (MS) lifetimes, the evolution of the surface rotation and abundances, as well as on the ejected masses of various isotopes. Among the new results obtained from the present grid we find that 1) fast-rotating stars with initial masses around 1.7 Msun present at the beginning of the core hydrogen-burning phase quite small convective cores with respect to their slowly rotating counterparts. This fact may be interesting to keep in mind in the framework of the asteroseismic studies of such stars. 2) The contrast between the core and surface angular velocity is higher in slower rotating stars. The values presently obtained are in agreement with the very few values obtained for B-type stars from asteroseismology. 3) At Z = 0.002, the stars in the mass range of 1.7 to 3 Msun with a mean velocity on the MS of the order of 150 km/s show N/H enhancement superior to 0.2 dex at mid-MS, and superior to 0.4 dex at the end of the MS phase. At solar metallicity the corresponding values are below 0.2 dex at any time in the MS.

Even though it is broadly accepted that single Be stars are rapidly rotating stars surrounded by a flat rotating circumstellar disk, there is still a debate about how fast these stars rotate and also about the mechanisms involved in the angular-momentum and mass input in the disk. We study the properties of stars that rotate near their critical-rotation rate and investigate the properties of the disks formed by equatorial mass ejections. We used the most recent Geneva stellar evolutionary tracks for rapidly rotating stars that reach the critical limit and used a simple model for the disk structure. We obtain that for a 9 Msun star at solar metallicity, the minimum average velocity during the Main Sequence phase to reach the critical velocity is around 330 km/s, whereas it would be 390 km/s at the metallicity of the Small Magellanic Cloud (SMC). Red giants or supergiants originating from very rapid rotators rotate six times faster and show N/C ratios three times higher than those originating from slowly rotating stars. This difference becomes stronger at lower metallicity. It might therefore be very interesting to study the red giants in clusters that show a large number of Be stars on the MS band. On the basis of our single-star models, we show that the observed Be-star fraction with cluster age is compatible with the existence of a temperature-dependent lower limit in the velocity rate required for a star to become a Be star. The mass, extension, and diffusion time of the disks produced when the star is losing mass at the critical velocity, obtained from simple parametrized expressions, are not too far from those estimated for disks around Be-type stars. At a given metallicity, the mass and the extension of the disk increase with the initial mass and with age on the MS phase. Denser disks are expected in low-metallicity regions.

Mass-loss rates are one of the most relevant parameters determining the evolution of massive stars. In particular, the rates at which the star loses mass during the red-supergiant (RSG) phase is the least constrained by the observations or theory. In this paper, we show how the mass loss during the RSG phase affects the later evolution of the star, as well as the final type of supernova towards which it leads. We also discuss some possibilities to discriminate between blue stars that went through a RSG phase and those which remained in the blue part of the Hertzsprung-Russell diagram.

We review the significant role played by red supergiants (RSGs) in stellar populations, and some challenges and questions they raise for theoretical stellar evolution. We present how metallicity and rotation modify the way stars go to the red part of the Hertzsprung- Russell diagram or come back from it, and how RSGs might keep a trace of their main-sequence evolution. We compare theoretical popu- lation ratios with observed ones.

The rotation of stars has many interesting and important consequences for the photometric and chemical evolution of galaxies. Many of the predictions of models of stellar rotation are now compared with observations of surface abundances and velocities, with interferometric studies of fast rotating stars, with internal rotation profiles as they can be deduced by asteroseismology, to cite just a few observational constraints. In this paper, we investigate how the outputs of models depend on the prescriptions used for the diffusion coefficients included in the shellular rotating models. After recalling the various prescriptions found in the literature, we discuss their impacts on the evolutionary tracks and lifetimes of the Main-Sequence (MS) phase, the changes of the surface composition and velocities during the MS phase, the distribution of the core helium lifetime in the blue and the red part of the HR diagram, the extensions of the blue loops, the evolution of the angular momentum of the core, and the synthesis of primary nitrogen in fast-rotating metal-poor massive stars. While some of these outputs depend only slightly on the prescriptions used (for instance, the evolution of the surface velocities), most of them show a significant dependence. The models which best fit the changes of the surface abundances are those computed with the vertical shear diffusion coefficient of Maeder (1997) and the horizontal shear diffusion coefficient by Zahn (1992).

Stars more massive than about 8 Msun end their lives as a Supernova (SN), an event of fundamental importance Universe-wide. Theoretically, these stars have been expected to be either at the red supergiant, blue supergiant, or Wolf-Rayet stage before the explosion. We performed coupled stellar evolution and atmospheric modeling of stars with initial masses between 20 Msun and 120 Msun. We found that the 20 Msun and 25 Msun rotating models, before exploding as SN, have spectra that do not resemble any of the aforementioned classes of massive stars. Rather, they have remarkable similarities with rare, unstable massive stars known as Luminous Blue Variables (LBV). While observations show that some SNe seem to have had LBVs as progenitors, no theoretical model had yet predicted that a star could explode at this stage. Our models provide theoretical support for relatively low-luminosity LBVs exploding as SN in the framework of single stellar evolution. This is a significant shift in paradigm, meaning that a fraction of LBVs could be the end stage of massive star evolution, rather than a transitory evolutionary phase. We suggest that type IIb SN could have LBV as progenitors, and a prime example could be SN 2008ax.

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. 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 \sim 19.0 MSun for the non rotating models and \sim 16.8 MSun 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\odot, 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.

The availability of asteroseismic constraints for a large sample of stars from the missions CoRoT and Kepler paves the way for various statistical studies of the seismic properties of stellar populations. In this paper, we evaluate the impact of rotation-induced mixing and thermohaline instability on the global asteroseismic parameters at different stages of the stellar evolution from the Zero Age Main Sequence to the Thermally Pulsating Asymptotic Giant Branch to distinguish stellar populations. We present a grid of stellar evolutionary models for four metallicities (Z = 0.0001, 0.002, 0.004, and 0.014) in the mass range between 0.85 to 6.0 Msun. The models are computed either with standard prescriptions or including both thermohaline convection and rotation-induced mixing. For the whole grid we provide the usual stellar parameters (luminosity, effective temperature, lifetimes, ...), together with the global seismic parameters, i.e. the large frequency separation and asymptotic relations, the frequency corresponding to the maximum oscillation power {\nu}_{max}, the maximal amplitude A_{max}, the asymptotic period spacing of g-modes, and different acoustic radii. We discuss the signature of rotation-induced mixing on the global asteroseismic quantities, that can be detected observationally. Thermohaline mixing whose effects can be identified by spectroscopic studies cannot be caracterized with the global seismic parameters studied here. But it is not excluded that individual mode frequencies or other well chosen asteroseismic quantities might help constraining this mixing.

We examine some properties of stars evolving close to the classical Eddington limit for electron-scattering opacity, when these stars maintain a chemically homogeneous structure as a result of mixing and/or mass loss. We consider analytical relations and models computed with the Geneva code. Homologous, chemically homogeneous stars evolving with a constant Eddington factor obey a relation of the form mu^2 M = const. This applies, for example, to Wolf-Rayet (WR) stars in stages without hydrogen. The value of the constant may depend on the metallicity, initial mass, evolutionary stage, and physical processes included in the considered homologous evolutionary sequence. An average value of the constant between 20 and 40 in solar units is consistent with the masses of Galactic WR stars.

[abridged] 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. 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 Msun. 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 ZAMS with a rotation rate Vini/Vcrit=0.4. The evolution is computed until the end of the central carbon-burning phase, the early 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. [...]

The Be star phenomenon is related to fast rotation, although the cause of this fast rotation is not yet clearly established. The basic effects of fast rotation on the stellar structure are reviewed: oblateness, mixing, anisotropic winds. The processes governing the evolution of the equatorial velocity of a single star (transport mechanisms and mass loss) are presented, as well as their metallicity dependence. The theoretical results are compared to observations of B and Be stars in the Galaxy and the Magellanic Clouds.

Be stars are thought to be fast rotating stars surrounded by an equatorial disc. The formation, structure and evolution of the disc are still not well understood. In the frame of single star models, it is expected that the surface of an initially fast rotating star can reach its keplerian velocity (critical velocity). The Geneva stellar evolution code has been recently improved, in order to obtain some estimates of the total mass loss and of the mechanical mass loss rates in the equatorial disc during the whole critical rotation phase. We present here the first results of the computation of a grid of fast rotating B stars evolving towards the Be phase, and discuss the first estimates we obtained.

Magnetic field is playing an important role at all stages of star evolution from star formation to the endpoints. The main effects are briefly reviewed. We also show that O-type stars have large convective envelopes, where convective dynamo could work. There, fields in magnetostatic balance have intensities of the order of 100 G. A few OB stars with strong polar fields (Henrichs et al. 2003a) show large N-enhancements indicating a strong internal mixing. We suggest that the meridional circulation enhanced by an internal rotation law close to uniform in these magnetic stars is responsible for the observed mixing. Thus, it is not the magnetic field itself which makes the mixing, but the strong thermal instability associated to solid body rotation. A critical question for evolution is whether a dynamo is at work in radiative zones of rotating stars. The Tayler-Spruit (TS) dynamo is the best candidate. We derive some basic relations for dynamos in radiative layers. Evolutionary models with TS dynamo show important effects: internal rotation coupling and enhanced mixing, all model outputs being affected.

Metal-free stars are assumed to evolve at constant mass because of the very low stellar winds. This leads to large CO-core mass at the end of the evolution, so primordial stars with an initial mass between 25 and 85 Msol are expected to end as direct black holes, the explosion energy being too weak to remove the full envelope. We show that when rotation enters into play, some mass is lost because the stars are prone to reach the critical velocity during the main sequence evolution. Contrarily to what happens in the case of very low- but non zero-metallicity stars, the enrichment of the envelope by rotational mixing is very small and the total mass lost remains modest. The compactness of the primordial stars lead to a very inefficient transport of the angular momentum inside the star, so the profile of Omega(r) is close to Omega r^2 = const. As the core contracts, the rotation rate increases, and the star ends its life with a fast spinning core. Such a configuration has been shown to modify substantially the dynamics of the explosion. Where one expected a weak explosion or none at all, rotation might boost the explosion energy and drive a robust supernova. This will have important consequences in the way primordial stars enriched the early Universe.

(Abridged) 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 equals 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. We have calculated 14 differentially-rotating and non-rotating stellar models at zero metallicity, with masses between 9 and 200 Msol. The evolution has been followed up to the pre-supernova stage. We find that Z=0 models rotate with an internal profile Omega(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 14N production. 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 sequence. 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 12C/13C ratios of -2.2, -0.95 and 50 respectively at log O/H +12 = 4.2.

Very massive primordial stars (140 Msol < M < 260 Msol) are supposed to end their lives as PISN. Such an event can be traced by a typical chemical signature in low metallicity stars, but at the present time, this signature is lacking in the extremely metal-poor stars we are able to observe. Does it mean that those very massive objects were not formed, contrarily to the primordial star formation scenarios ? Could it be possible that they avoided this tragical fate ? We explore the effects of rotation, anisotropical mass loss and magnetic field on the core size of very massive Population III models. We find that magnetic fields provide the strong coupling that is lacking in standard evolution metal-free models and our 150 Msol Population III model avoids indeed the pair-instability explosion.

(abridged) Fast stellar rotation is currently the most promising mechanism for producing primary nitrogen in metal-poor massive stars. Chemical evolution models computed with the inclusion of the yields of fast rotating models at a metallicity $Z=10^{-8}$ can account for the high N/O abundances observed in normal metal-poor halo stars. If, as believed, intermediate mass stars did not have enough time to contribute to the interstellar medium enrichment at such low metallicities, the above result constitutes a strong case for the existence of fast rotators in the primordial Universe. An important result of stellar models of fast rotators is that large quantities of primary 13C are produced. Hence, our goal is to investigate the consequence of fast rotation on the evolution of the 12C/13C ratio in the interstellar medium at low metallicity. We predict that, if fast rotating massive stars were common phenomena in the early Universe, the primordial interstellar medium of galaxies with a star formation history similar to the one inferred for our galactic halo should have 12C/13C ratios between 30-300. Without fast rotators, the predicted 12C/13C ratios would be $\sim$ 4500 at [Fe/H]=-3.5, increasing to $\sim$ 31000 at around [Fe/H]=-5.0. Current data on very metal-poor giant normal stars in the galactic halo agree better with chemical evolution models including fast rotators. To test our predictions, challenging measurements of the 12C/13C in more extremely metal-poor giants and turnoff stars are required.

We determine which are the mechanisms accelerating the surface of single stars during the Main Sequence evolution. We have computed 112 stellar models of four different initial masses between 3 and 60 M_sun, at four different metallicities between 0 and 0.020, and with seven different values of the ratio Omega/Omega_crit between 0.1 and 0.99. For all the models, computations were performed until either the end of the Main Sequence evolution or the reaching of the critical limit. The evolution of surface velocities during the Main Sequence lifetime results from an interplay between meridional circulation (bringing angular momentum to the surface) and mass loss by stellar winds (removing it). The dependence on metallicity of these two mechanisms plays a key role in determining for each metallicity, a limiting range of initial masses (spectral types) for stars able to reach or at least approach the critical limit. Present models predict a higher frequency of fast rotating stars in clusters with ages between 10 and 25 Myr. This is the range of ages where most of Be stars are observed.

Very massive primordial stars ($140 M_{\odot} < M < 260 M_{\odot}$) are supposed to end their lives as pair-instability supernovae. Such an event can be traced by a typical chemical signature in low metallicity stars, but at the present time, this signature is lacking in the extremely metal-poor stars we are able to observe. Does it mean that those very massive objects did not form, contrarily to the primordial star formation scenarios? Could they avoid this tragical fate? We explore the effects of rotation, anisotropic mass loss and magnetic fields on the core size of a very massive Population III model, in order to check if its mass is sufficiently modified to prevent the pair instability. We obtain that a Population III model of $150 M_{\odot}$ with $\upsilon/\upsilon_{\rm crit}=0.56$ computed with the inclusion of wind anisotropy and Tayler-Spruit dynamo avoids the pair instability explosion.

We propose the Wind of Fast Rotating Massive Stars scenario to explain the origin of the abundance anomalies observed in globular clusters. We compute and present models of fast rotating stars with initial masses between 20 and 120 Msun for an initial metallicity Z=0.0005 ([Fe/H]=-1.5). We discuss the nucleosynthesis in the H-burning core of these objects and present the chemical composition of their ejecta. We consider the impact of uncertainties in the relevant nuclear reaction rates. Fast rotating stars reach the critical velocity at the beginning of their evolution and remain near the critical limit during the rest of the main sequence and part of the He-burning phase. As a consequence they lose large amounts of material through a mechanical wind which probably leads to the formation of a slow outflowing disk. The material in this slow wind is enriched in H-burning products and presents abundance patterns similar to the chemical anomalies observed in globular cluster stars. In particular, the C, N, O, Na and Li variations are well reproduced by our model. However the rate of the 24Mg(p,gamma) has to be increased by a factor 1000 around 50 MK in order to reproduce the whole amplitude of the observed Mg-Al anticorrelation. We discuss how the long-lived low-mass stars currently observed in globular clusters could have formed out of the slow wind material ejected by massive stars.

Although the theoretical study of very low metallicity (Z) and metal-free stars is not new, their importance has recently greatly increased since two related fields have been developing rapidly. The first is cosmological simulations of the formation of the first stars and of the reionisation period. The second is the observations of extremely metal poor stars. In this paper, we present pre-supernova evolution models of massive rotating stars at very low Z (Z=1e-8) and at Z=0. Rotation has a strong impact on mass loss and nucleosynthesis. Models reaching break-up velocities lose up to ten percents of their initial mass. In very low Z models, rotational and convective mixing enhances significantly the surface content in carbon, nitrogen and oxygen (CNO) when the star becomes a red supergiant. This induces a strong mass loss for stars more massive than about 60 solar masses. Our models predict type Ib,c supernovae and gamma-ray bursts at very low Z. Rotational mixing also induces a large production of CNO elements, in particular of primary nitrogen. The stellar wind chemical composition is compatible with the most metal-poor star know to date, HE 1327-2326, for CNO elements. Our models reproduce the early evolution of nitrogen in the Milky Way.

Big Bang nucleosynthesis produces only light elements and the very first generation stars are thus formed from metal-free clouds. They start the production of heavy elements during their life, and enrich the interstellar medium through their explosive death. Stellar evolution models show that the treatment of rotation has important effects on the evolution of those metal-free stars: for example, rotating models produce up to five orders of magnitude more primary nitrogen than non rotating models, due to internal mixing. This will have an impact in the composition of the second generation stars, some of which may now be observed in the Galactic halo. In the case Population III stars were very massive and would end up as direct black holes, rotation again have an interesting effect of enhancing mass loss through centrifugal force and surface enrichment. CNO composition patterns observed in ultra metal-poor halo stars may be explained by a 'wind only' contribution.

We investigate the effect of new stellar models, which take rotation into account, computed for a metallicity Z = 10^{-8} on the chemical evolution of the earliest phases of the Milky Way. These models are 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 interstellar medium nitrogen enrichment for log(O/H)+12 < 7 (or [Fe/H]< -3).We show that upon the inclusion of the 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}.

The aim of this paper is to determine the initial rotational velocities required on the ZAMS for single stars to reach the critical velocity, sometimes called the break-up velocity, during the Main-Sequence (MS) phase. Some useful relations between $\Omega/\Omega_{\rm crit}$, $\upsilon/\upsilon_{\rm crit}$ ($\upsilon$ is the velocity at the equator), the moments of inertia, the angular momenta, the kinetic energy in the rotation and various other basic physical quantities are obtained.

The effects of rotation on low-metallicity stellar models are twofold: first, the models reach break-up during main sequence and may lose mass by mechanical process; second, strong internal mixing brings freshly synthesized elements towards the surface and raises the effective metallicity to higher values, so that the initially very low radiative winds are enhanced. Those two effects are also found in Z=0 models, but to a lesser degree because of structural differences. This weak mass loss becomes interesting though in the case of the black hole-doomed stars (25-140 M_\sun and > 260 M_\sun) because it allows these stars to contribute to the enrichment of the medium by stellar winds.

Rotation deeply affects the evolution of very metal poor massive stars. Indeed, even moderately rotating stars reach the break--up limit during the Main--Sequence (MS) phase, they evolve rapidly to the red after the core H--burning phase and important surface enrichment in CNO elements occurs at the supergiant stage. High rotation and the self enrichment in heavy elements $Z$ tend to enhance the quantity of mass lost by stellar winds. As a numerical example, we obtain that a 60 M$_\odot$ at $Z=0.00001$, with an initial rotation corresponding to 0.7 of the critical velocity, loses half of its initial mass during its lifetime. The ejected material, enriched in helium and CNO elements, has interesting consequences for the early chemical evolution of galaxies.