In this work, we investigate the impact of uncertainties due to convective boundary mixing (CBM), commonly called `overshoot', namely the boundary location and the amount of mixing at the convective boundary, on stellar structure and evolution. For this we calculated two grids of stellar evolution models with the MESA code, each with the Ledoux and the Schwarzschild boundary criterion, and vary the amount of CBM. We calculate each grid with the initial masses $15$, $20$ and $25\,\rm{M}_\odot$. We present the stellar structure of the models during the hydrogen and helium burning phases. In the latter, we examine the impact on the nucleosynthesis. We find a broadening of the main-sequence with more CBM, which is more in agreement with observations. Furthermore during the core hydrogen burning phase there is a convergence of the convective boundary location due to CBM. The uncertainties of the intermediate convective zone remove this convergence. The behaviour of this convective zone strongly affects the surface evolution of the model, i.e. how fast it evolves red-wards. The amount of CBM impacts the size of the convective cores and the nucleosynthesis, e.g. the $^{12}$C to $^{16}$O ratio and the weak s-process. Lastly, we determine the uncertainty that the range of parameter values investigated introduce and we find differences of up to $70\%$ for the core masses and the total mass of the star.

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.

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 time evolution of angular momentum and surface rotation of massive stars is strongly influenced by fossil magnetic fields via magnetic braking. We present a new module containing a simple, comprehensive implementation of such a field at the surface of a massive star within the Modules for Experiments in Stellar Astrophysics (MESA) software instrument. We test two limiting scenarios for magnetic braking: distributing the angular momentum loss throughout the star in the first case, and restricting the angular momentum loss to a surface reservoir in the second case. We perform a systematic investigation of the rotational evolution using a grid of OB star models with surface magnetic fields ($M_\star=5-60$ M$_\odot$, $\Omega/\Omega_{\rm crit} =0.2-1.0$, $B_{\rm p} =1-20$ kG). We then employ a representative grid of B-type star models ($M_\star=5, 10, 15$ M$_\odot$, $\Omega/\Omega_{\rm crit} =0.2 , 0.5, 0.8$, $B_{\rm p} = 1, 3 ,10, 30$ kG) to compare to the results of a recent self-consistent analysis of the sample of known magnetic B-type stars. We infer that magnetic massive stars arrive at the zero age main sequence with a range of rotation rates, rather than with one common value. In particular, some stars are required to have close-to-critical rotation at the ZAMS. However, magnetic braking yields surface rotation rates converging to a common low value, making it difficult to infer the initial rotation rates of evolved, slowly-rotating stars.

We present a detailed study of stellar rotation in the massive 1.5 Gyr old cluster NGC 1846 in the Large Magellanic Cloud. Similar to other clusters at this age, NGC 1846 shows an extended main sequence turn-off (eMSTO), and previous photometric studies have suggested it could be bimodal. In this study, we use MUSE integral-field spectroscopy to measure the projected rotational velocities (vsini) of around 1400 stars across the eMSTO and along the upper main sequence of NGC 1846. We measure vsini values up to ~250 km/s and find a clear relation between the vsini of a star and its location across the eMSTO. Closer inspection of the distribution of rotation rates reveals evidence for a bimodal distribution, with the fast rotators centred around vsini = 140 km/s and the slow rotators centred around vsini = 60 km/s. We further observe a lack of fast rotating stars along the photometric binary sequence of NGC 1846, confirming results from the field that suggest that tidal interactions in binary systems can spin down stars. However, we do not detect a significant difference in the binary fractions of the fast and slowly rotating sub-populations. Finally, we report on the serendipitous discovery of a planetary nebula associated with NGC 1846.

We find an asymptotic limit for the integral scale dissipation length in strongly stratified stellar convection; this is not adjustable, and we identify it with the "mixing length parameter", $\alpha \sim 5/3$, a theoretical prediction which agrees to within 10% with the fitted empirical values quoted by eight independent stellar evolutionary groups, as well as several 3D stellar atmosphere groups. For strong stratification, the dissipation length approaches the density scale height; in weak stratification it shrinks to the depth for such thin convective regions. Simulations with different zoning give a monotonic decrease in resolution errors with added zones, down to barely detectable for the most difficult case (lower boundary layer), and an estimate of the dissipation due to the turbulent cascade. No explicit viscosity is required; stellar dissipation occurs from nonlinear fluid effects (Taylor, Onsager). These implicit large eddy simulations of the Euler equations, with simultaneous analysis by Reynolds averaging (RA-ILES), resolve the energy-containing eddies, develop a turbulent cascade down to the grid scale, and produce surfaces of separation at boundary layers. Like experiments (Warhaft), and direct numerical simulations of the Navier-Stokes equations (DNS, Sreenivasan), our simulations develop "anomalous terms" for dissipation in the turbulent cascade (due to intermittency, anisotropy, and interactions between coherent structures). They have more realistic behavior at high order than the K41 theory of Kolmogorov. Mixing length theory (MLT, Boehm-Vitense) overestimates the convective enthalpy flux for strong stratification, which has consequences for the theory of luminous blue variables, for the solar metal abundance, and deep solar convection.

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

We present an extended grid of state-of-the art stellar models for low-mass stars including updated physics (nuclear reaction rates, surface boundary condition, mass-loss rate, angular momentum transport, torque and rotation-induced mixing prescriptions). We aim at evaluating the impact of wind braking, realistic atmospheric treatment, rotation and rotation-induced mixing on the structural and rotational evolution from the pre-main sequence to the turn-off. Using the STAREVOL code, we provide an updated PMS grid. We compute stellar models for 7 different metallicities, from [Fe/H] = -1 dex to [Fe/H] = +0.3 dex with a solar composition corresponding to $Z=0.0134$. The initial stellar mass ranges from 0.2 to 1.5\Ms with extra grid refinement around one solar mass. We also provide rotating models for three different initial rotation rates (slow, median and fast) with prescriptions for the wind braking and disc-coupling timescale calibrated on observed properties of young open clusters. The rotational mixing includes an up-to-date description of the turbulence anisotropy in stably stratified regions. The overall behaviour of our models at solar metallicity -- and its constitutive physics -- is validated through a detailed comparison with a variety of distributed evolutionary tracks. The main differences arise from the choice of surface boundary conditions and initial solar composition. The models including rotation with our prescription for angular momentum extraction and self-consistent formalism for angular momentum transport are able to reproduce the rotation period distribution observed in young open clusters over a broad mass-range. These models are publicly available and may be used to analyse data coming from present and forthcoming asteroseismic and spectroscopic surveys such as Gaia, TESS and PLATO.

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.

Surface magnetic fields have a strong impact on stellar mass loss and rotation and, as a consequence, on the evolution of massive stars. In this work we study the influence of an evolving dipolar surface fossil magnetic field with an initial field strength of 4 kG on the characteristics of 15 M$_{\odot}$ solar metallicity models using the Geneva stellar evolution code. Non-rotating and rotating models considering two different scenarios for internal angular momentum transport are computed, including magnetic field evolution, mass-loss quenching, and magnetic braking. Magnetic field evolution results in weakening the initially strong magnetic field, however, in our models an observable magnetic field is still maintained as the star evolves towards the red supergiant phase. At the given initial mass of the models, mass-loss quenching is modest. Magnetic braking greatly enhances chemical element mixing if radial differential rotation is allowed for, on the other hand, the inclusion of surface magnetic fields yields a lower surface enrichment in the case of near solid-body rotation. Models including surface magnetic fields show notably different trends on the Hunter diagram (plotting nitrogen abundance vs $v \sin i$) compared to those that do not. The magnetic models agree qualitatively with the anomalous `Group 2 stars', showing slow surface rotation and high surface nitrogen enhancement on the main sequence.

Convective boundary mixing is one of the major uncertainties in stellar evolution. In order to study its dependence on boundary properties and turbulence strength in a controlled way, we computed a series of 3D hydrodynamical simulations of stellar convection during carbon burning with a varying boosting factor of the driving luminosity. Our 3D implicit large eddy simulations were computed with the PROMPI code. We performed a mean field analysis of the simulations within the Reynolds-averaged Navier-Stokes framework. Both the vertical RMS velocity within the convective region and the bulk Richardson number of the boundaries are found to scale with the driving luminosity as expected from theory. The positions of the convective boundaries were estimated through the composition profiles across them, and the strength of convective boundary mixing was determined by analysing the boundaries within the framework of the entrainment law. We find that the entrainment is approximately inversely proportional to the bulk Richardson number. Although the entrainment law does not encompass all the processes occurring at boundaries, our results support the use of the entrainment law to describe convective boundary mixing in 1D models, at least for the advanced phases. The next steps and challenges ahead are also discussed.

The semi-empirical initial-final mass relation (IFMR) connects spectroscopically analyzed white dwarfs in star clusters to the initial masses of the stars that formed them. Most current stellar evolution models, however, predict that stars will evolve to white dwarfs $\sim$0.1 M$_\odot$ less massive than that found in the IFMR. We first look at how varying theoretical mass-loss rates, third dredge-up efficiencies, and convective-core overshoot may help explain the differences between models and observations. These parameters play an important role at the lowest masses (M$_{\rm initial}$ $<$ 3 M$_\odot$). At higher masses, only convective-core overshoot meaningfully affects white dwarf mass, but alone it likely cannot explain the observed white dwarf masses nor why the IFMR scatter is larger than observational errors predict. These higher masses, however, are also where rotational mixing in main sequence stars begins to create more massive cores, and hence more massive white dwarfs. This rotational mixing also extends a star's lifetime, making faster rotating progenitors appear like less massive stars in their semi-empirical age analysis. Applying the observed range of young B-dwarf rotations to the MIST or SYCLIST rotational models demonstrates a marked improvement in reproducing both the observed IFMR data and its scatter. The incorporation of both rotation and efficient convective-core overshoot significantly improves the match with observations. This work shows that the IFMR provides a valuable observational constraint on how rotation and convective-core overshoot affect the core evolution of a star.

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.

Renzini (1987) wrote an influential critique of mixing-length theory (MLT) as used in stellar evolution codes, and concluded that three-dimensional (3D) fluid dynamical simulations were needed to clarify several important issues. We have critically explored the limitations of the numerical methods and conclude that they are approaching the required accuracy. Implicit large eddy simulations (ILES) automatically connect large scale turbulence to a Kolmogorov cascade below the grid scale, allowing turbulent boundary layers to remove singularities that appear in the theory. Interactions between coherent structures give multi-modal behavior, driving intermittency and fluctuations. Reynolds averaging (RA) allows us to abstract the essential features of this dynamical behavior of boundaries which are appropriate to stellar evolution, and consider how they relate static boundary conditions (Richardson, Schwarzschild or Ledoux). We clarify several questions concerning when and why MLT works, and does not work, using both analytical theory and 3D high resolution numerical simulations. The composition gradients and boundary layer structure which are produced by our simulations suggest a self-consistent approach to boundary layers, removing the need for ad hoc procedures for 'convective overshooting' and `semi-convection'. In a companion paper we quantify the adequacy of our numerical resolution, determine of the length scale of dissipation (the `mixing length') without astronomical calibration, quantify agreement with the four-fifths law of Kolmogorov for weak stratification, and extend MLT to deal with strong stratification.

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.

Convection plays a key role in the evolution of stars due to energy transport and mixing of composition. Despite its importance, this process is still not well understood. One longstanding conundrum in all 1D stellar evolution codes is the treatment of convective boundaries. In this study we compare two convective uncertainties, the boundary location (Ledoux versus Schwarzschild) and the amount of extra mixing, and their impact on the early evolution of massive stars. With increasing convective boundary mixing (CBM), we find a convergence of the two different boundary locations, a decreasing blue to red super giant ratio and a reduced importance of semiconvection.

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

We recently discovered a yellow supergiant (YSG) in the Small Magellanic Cloud (SMC) with a heliocentric radial velocity of ~300 km/s which is much larger than expected for a star in its location in the SMC. This is the first runaway YSG ever discovered and only the second evolved runaway star discovered in a different galaxy than the Milky Way. We classify the star as G5-8I, and use de-reddened broad-band colors with model atmospheres to determine an effective temperature of 4700+/-250K, consistent with what is expected from its spectral type. The star's luminosity is then L/Lo ~ 4.2+/-0.1, consistent with it being a ~30Myr 9Mo star according to the Geneva evolution models. The star is currently located in the outer portion of the SMC's body, but if the star's transverse peculiar velocity is similar to its peculiar radial velocity, in 10Myr the star would have moved 1.6 degrees across the disk of the SMC, and could easily have been born in one of the SMC's star-forming regions. Based on its large radial velocity, we suggest it originated in a binary system where the primary exploded as a supernovae thus flinging the runaway star out into space. Such stars may provide an important mechanism for the dispersal of heavier elements in galaxies given the large percentage of massive stars that are runaways. In the future we hope to look into additional evolved runaway stars that were discovered as part of our other past surveys.

Massive stars have a strong impact on their surroundings, in particular when they produce a core-collapse supernova at the end of their evolution. In these proceedings, we review the general evolution of massive stars and their properties at collapse as well as the transition between massive and intermediate-mass stars. We also summarise the effects of metallicity and rotation. We then discuss some of the major uncertainties in the modelling of massive stars, with a particular emphasis on the treatment of convection in 1D stellar evolution codes. Finally, we present new 3D hydrodynamic simulations of convection in carbon burning and list key points to take from 3D hydrodynamic studies for the development of new prescriptions for convective boundary mixing in 1D stellar evolution codes.

Be stars are main-sequence massive stars with emission features in their spectrum, which originates in circumstellar gaseous discs. Even though the viscous decretion disc (VDD) model can satisfactorily explain most observations, two important physical ingredients, namely the magnitude of the viscosity ($\alpha$) and the disk mass injection rate, remain poorly constrained. The light curves of Be stars that undergo events of disc formation and dissipation offer an opportunity to constrain these quantities. A pipeline was developed to model these events that uses a grid of synthetic light curves, computed from coupled hydrodynamic and radiative transfer calculations. A sample of 54 Be stars from the OGLE survey of the Small Magellanic Cloud (SMC) was selected for this study. Because of the way our sample was selected (bright stars with clear disc events), it likely represents the densest discs in the SMC. Like their siblings in the Galaxy, the mass of the disc in the SMC increases with the stellar mass. The typical mass and angular momentum loss rates associated with the disk events are of the order of $\sim$$10^{-10}\, M_\odot\,\mathrm{yr^{-1}}$ and $\sim$$5\times 10^{36}\, \mathrm{g\, cm^{2}\, s^{-2}}$, respectively. The values of $\alpha$ found in this work are typically of a few tenths, consistent with recent results in the literature and with the ones found in dwarf novae, but larger than current theory predicts. Considering the sample as a whole, the viscosity parameter is roughly two times larger at build-up ($\left\langle\alpha_\mathrm{bu}\right\rangle = 0.63$) than at dissipation ($\left\langle\alpha_\mathrm{d}\right\rangle = 0.26$). Further work is necessary to verify whether this trend is real or a result of some of the model assumptions.

We present the initial results of the Large Impact of magnetic Fields on the Evolution of hot stars (LIFE) project. The focus of this project is the search for magnetic fields in evolved OBA giants and supergiants with visual magnitudes between 4 and 8, with the aim to investigate how the magnetic fields observed in upper main sequence (MS) stars evolve from the MS until the late post-MS stages. In this paper, we present spectropolarimetric observations of 15 stars observed using the ESPaDOnS instrument of the CFHT. For each star, we have determined the fundamental parameters and have used stellar evolution models to calculate their mass, age and radius. Using the LSD technique, we have produced averaged line profiles for each star. From these profiles, we have measured the longitudinal magnetic field strength and have calculated the detection probability. We report the detection of magnetic fields in two stars of our sample: a weak field of $B_l=1.0\pm0.2\,$G is detected in the post-MS A5 star 19$\,$Aur and a stronger field of $B_l=-230\,\pm\,10\,$G is detected in the MS/post-MS B8/9 star HR$\,$3042.

We present the first detailed three-dimensional hydrodynamic implicit large eddy simulations of turbulent convection for carbon burning. The simulations start with an initial radial profile mapped from a carbon burning shell within a 15 solar mass stellar evolution model. We considered 4 resolutions from 128^3 to 1024^3 zones. These simulations confirm that convective boundary mixing (CBM) occurs via turbulent entrainment as in the case of oxygen burning. The expansion of the boundary into the surrounding stable region and the entrainment rate are smaller at the bottom boundary because it is stiffer than the upper boundary. The results of this and similar studies call for improved CBM prescriptions in 1D stellar evolution models.

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.

Magnetic confinement of the winds of hot, massive stars has far-reaching consequences on timescales ranging from hours to Myr. Understanding the long-term effects of this interplay has already led to the identification of two new evolutionary pathways to form `heavy' stellar mass black holes and pair-instability supernova even at galactic metallicity. We are performing 1D stellar evolution model calculations that, for the first time, account for the surface effects and the time evolution of fossil magnetic fields. These models will be thoroughly confronted with observations and will potentially lead to a significant revision of the derived parameters of observed magnetic massive stars.

The physics of massive stars depends (at least) on convection, mass loss by stellar winds, rotation, magnetic fields and multiplicity. We briefly discuss the impact of the first three processes on the stellar yields trying to identify some guidelines for future works.

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

In the context of the high resolution, high signal-to-noise ratio, high sensitivity, spectropolarimetric survey BritePol, which complements observations by the BRITE constellation of nanosatellites for asteroseismology, we are looking for and measuring the magnetic field of all stars brighter than V=4. In this paper, we present circularly polarised spectra obtained with HarpsPol at ESO in La Silla (Chile) and ESPaDOnS at CFHT (Hawaii) for 3 hot evolved stars: $\iota$ Car, HR 3890, and $\epsilon$ CMa. We detected a magnetic field in all 3 stars. Each star has been observed several times to confirm the magnetic detections and check for variability. The stellar parameters of the 3 objects were determined and their evolutionary status was ascertained employing evolution models computed with the Geneva code. $\epsilon$ CMa was already known and is confirmed to be magnetic, but our modeling indicates that it is located near the end of the main sequence, i.e. it is still in a core hydrogen burning phase. $\iota$ Car and HR 3890 are the first discoveries of magnetic hot supergiants located well after the end of the main sequence on the HR diagram. These stars are probably the descendants of main sequence magnetic massive stars. Their current field strength (a few G) is compatible with magnetic flux conservation during stellar evolution. These results provide observational constraints for the development of future evolutionary models of hot stars including a fossil magnetic field.

Aims: Past observations of fast-rotating massive stars exhibiting normal nitrogen abundances at their surface have raised questions about the rotational mixing paradigm. We revisit this question thanks to a spectroscopic analysis of a sample of bright fast-rotating OB stars, with the goal of quantifying the efficiency of rotational mixing at high rotation rates. Method: Our sample consists of 40 fast rotators on the main sequence, with spectral types comprised between B0.5 and O4. We compare the abundances of some key element indicators of mixing (He, CNO) with the predictions of evolutionary models for single objects and for stars in interacting binary systems. Results: The properties of half of the sample stars can be reproduced by single evolutionary models, even in the case of probable or confirmed binaries that can therefore be true single stars in a pre-interaction configuration. The main problem for the rest of the sample is a mismatch for the [N/O] abundance ratio (we confirm the existence of fast rotators with a lack of nitrogen enrichment) and/or a high helium abundance that cannot be accounted for by models. Modifying the diffusion coefficient implemented in single- star models does not solve the problem as it cannot simultaneously reproduce the helium abundances and [N/O] abundance ratios of our targets. Since part of them actually are binaries, we also compared their chemical properties with predictions for post-mass transfer systems. We found that these models can explain the abundances measured for a majority of our targets, including some of the most helium-enriched, but fail to reproduce them in other cases. Our study thus reveals that some physical ingredients are still missing in current models.

We present the first detailed three-dimensional (3D) hydrodynamic implicit large eddy simulations of turbulent convection of carbon burning in massive stars. Simulations begin with radial profiles mapped from a carbon burning shell within a 15$\,\textrm{M}_\odot$ one-dimensional stellar evolution model. We consider models with $128^3$, $256^3$, $512^3$ and $1024^3$ zones. The turbulent flow properties of these carbon burning simulations are very similar to the oxygen burning case. We performed a mean field analysis of the kinetic energy budgets within the Reynolds-averaged Navier-Stokes framework. For the upper convective boundary region, we find that the numerical dissipation is insensitive to resolution for linear mesh resolutions above 512 grid points. For the stiffer, more stratified lower boundary, our highest resolution model still shows signs of decreasing sub-grid dissipation suggesting it is not yet numerically converged. We find that the widths of the upper and lower boundaries are roughly 30% and 10% of the local pressure scale heights, respectively. The shape of the boundaries is significantly different from those used in stellar evolution models. As in past oxygen-shell burning simulations, we observe entrainment at both boundaries in our carbon-shell burning simulations. In the large P\'eclet number regime found in the advanced phases, the entrainment rate is roughly inversely proportional to the bulk Richardson number, Ri$_{\rm B}$ ($\propto $Ri${\rm_B}^{-\alpha}$, $0.5\lesssim \alpha \lesssim 1.0$). We thus suggest the use of Ri$_{\rm B}$ as a means to take into account the results of 3D hydrodynamics simulations in new 1D prescriptions of convective boundary mixing.

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.

The surface rotations of some red giants are so fast that they must have been spun up by tidal interaction with a close companion, either another star, a brown dwarf, or a planet. We focus here on the case of red giants that are spun up by tidal interaction with a planet. When the distance between the planet and the star decreases, the spin period of the star decreases, the orbital period of the planet decreases, and the reflex motion of the star increases. We study the change rate of these three quantities when the circular orbit of a planet of 15 M$_{J}$ that initially orbits a 2 M$_\odot$ star at 1 au shrinks under the action of tidal forces during the red giant phase. We use stellar evolution models coupled with computations of the orbital evolution of the planet, which allows us to follow the exchanges of angular momentum between the star and the orbit in a consistent way. We obtain that the reflex motion of the red giant star increases by more than 1 m s$^{-1}$ per year in the last $\sim$40 years before the planet engulfment. During this phase, the reflex motion of the star is between 660 and 710 m s$^{-1}$. The spin period of the star increases by more than about 10 minutes per year in the last 3000 y before engulfment. During this period, the spin period of the star is shorter than 0.7 year. During this same period, the variation in orbital period, which is shorter than 0.18 year, is on the same order of magnitude. Changes in reflex-motion and spin velocities are very small and thus most likely out of reach of being observed. The most promising way of detecting this effect is through observations of transiting planets, that is, through{\it } changes of the beginning or end of the transit. A space mission like PLATO might be of great interest for detecting planets that are on the verge of being engulfed by red giants.

The treatment of mixing processes is still one of the major uncertainties in 1D stellar evolution models. This is mostly due to the need to parametrize and approximate aspects of hydrodynamics in hydrostatic codes. In particular, the effect of hydrodynamic instabilities in rotating stars, for example, dynamical shear instability, evades consistent description. We intend to study the accuracy of the diffusion approximation to dynamical shear in hydrostatic stellar evolution models by comparing 1D models to a first-principle hydrodynamics simulation starting from the same initial conditions. We chose an initial model calculated with the stellar evolution code GENEC that is just at the onset of a dynamical shear instability but does not show any other instabilities (e.g., convection). This was mapped to the hydrodynamics code SLH to perform a 2D simulation in the equatorial plane. We compare the resulting profiles in the two codes and compute an effective diffusion coefficient for the hydro simulation. Shear instabilities develop in the 2D simulation in the regions predicted by linear theory to become unstable in the 1D model. Angular velocity and chemical composition is redistributed in the unstable region, thereby creating new unstable regions. Eventually the 2D simulation settles in a symmetric, steady state, which is Richardson stable everywhere, whereas the instability remains for longer in the 1D model due to current limitations in the 1D code. A spatially resolved diffusion coefficient is extracted by comparing the initial and final profiles of mean atomic mass. The presented simulation gives a first insight on hydrodynamics of shear instabilities in a real stellar environment and even allows us to directly extract an effective diffusion coefficient. We see evidence for a critical Richardson number of 0.25 as regions above this value remain stable for the course of the simulation.

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.

Over the last decade, tremendous strides have been achieved in our understanding of magnetism in main sequence hot stars. In particular, the statistical occurrence of their surface magnetism has been established (~10%) and the field origin is now understood to be fossil. However, fundamental questions remain: how do these fossil fields evolve during the post-main sequence phases, and how do they influence the evolution of hot stars from the main sequence to their ultimate demise? Filling the void of known magnetic evolved hot (OBA) stars, studying the evolution of their fossil magnetic fields along stellar evolution, and understanding the impact of these fields on the angular momentum, rotation, mass loss, and evolution of the star itself, is crucial to answering these questions, with far reaching consequences, in particular for the properties of the precursors of supernovae explosions and stellar remnants. In the framework of the BRITE spectropolarimetric survey and LIFE project, we have discovered the first few magnetic hot supergiants. Their longitudinal surface magnetic field is very weak but their configuration resembles those of main sequence hot stars. We present these first observational results and propose to interpret them at first order in the context of magnetic flux conservation as the radius of the star expands with evolution. We then also consider the possible impact of stellar structure changes along evolution.

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.

Recent photometric analysis of the colour-magnitude diagrams (CMDs) of young massive clusters (YMCs) have found evidence for splitting in the main sequence and extended main sequence turn-offs, both of which have been suggested to be caused by stellar rotation. Comparison of the observed main sequence splitting with models has led various authors to suggest a rather extreme stellar rotation distribution, with a minority ($10-30$\%) of stars with low rotational velocities and the remainder ($70-90$\%) of stars rotating near the critical rotation (i.e., near break-up). We test this hypothesis by searching for Be stars within two YMCs in the LMC (NGC 1850 and NGC 1856), which are thought to be critically rotating stars with decretion disks that are (partially) ionised by their host stars. In both clusters we detect large populations of Be stars at the main sequence turn-off ($\sim30-60$\% of stars), which supports previous suggestions of large populations of rapidly rotating stars within massive clusters.

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.

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.

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.

Stellar evolution models of massive stars are important for many areas of astrophysics, for example nucleosynthesis yields, supernova progenitor models and understanding physics under extreme conditions. Turbulence occurs in stars primarily due to nuclear burning at different mass coordinates within the star. The understanding and correct treatment of turbulence and turbulent mixing at convective boundaries in stellar models has been studied for decades but still lacks a definitive solution. This paper presents initial results of a study on convective boundary mixing (CBM) in massive stars. The 'stiffness' of a convective boundary can be quantified using the bulk Richardson number ($\textrm{Ri}_B$), the ratio of the potential energy for restoration of the boundary to the kinetic energy of turbulent eddies. A 'stiff' boundary ($\textrm{Ri}_B \sim 10^4$) will suppress CBM, whereas in the opposite case a 'soft' boundary ($\textrm{Ri}_B \sim 10$) will be more susceptible to CBM. One of the key results obtained so far is that lower convective boundaries (closer to the centre) of nuclear burning shells are 'stiffer' than the corresponding upper boundaries, implying limited CBM at lower shell boundaries. This is in agreement with 3D hydrodynamic simulations carried out by Meakin and Arnett [The Astrophysical Journal 667:448-475, 2007]. This result also has implications for new CBM prescriptions in massive stars as well as for nuclear burning flame front propagation in Super-Asymptotic Giant Branch stars and also the onset of novae.

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.

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.

We present late-time Hubble Space Telescope (HST) ultraviolet (UV) and optical observations of the site of SN 2011dh in the galaxy M51, ~1164 days post-explosion. At the SN location, we observe a point source that is visible at all wavelengths, that is significantly fainter than the spectral energy distribution (SED) of the Yellow Supergiant progenitor observed prior to explosion. The previously reported photometry of the progenitor is, therefore, completely unaffected by any sources that may persist at the SN location after explosion. In comparison with the previously reported late-time photometric evolution of SN 2011dh, we find that the light curve has plateaued at all wavelengths. The SED of the late-time source is clearly inconsistent with a SED of stellar origin. Although the SED is bright at UV wavelengths, there is no strong evidence that the late-time luminosity originates solely from a stellar source corresponding to the binary companion, although a partial contribution to the observed UV flux from a companion star can not be ruled out.

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.

(abridged) The flux-weighted gravity-luminosity relationship (FGLR) of blue supergiant stars (BSG) links their absolute magnitude to the spectroscopically determined flux-weighted gravity log g = Teff^4. BSG are the brightest stars in the universe at visual light and the application of the FGLR has become a powerful tool to determine extragalactic distances. Observationally, the FGLR is a tight relationship with only small scatter. It is, therefore, ideal to be used as a constraint for stellar evolution models. The goal of this work is to investigate whether stellar evolution can reproduce the observed FGLR and to develop an improved foundation of the FGLR as an extragalactic distance indicator. We use different grids of stellar models for initial masses between 9 and 40 Msun, for metallicities between Z = 0.002 and 0.014, with and without rotation, computed with various mass loss rates during the red supergiant phase. For each of these models we discuss the details of post-main sequence evolution and construct theoretical FGLRs by means of population synthesis models which we then compare with the observed FGLR. In general, the stellar evolution model FGLRs agree reasonably well with the observed one. There are, however, differences between the models, in particular with regard to the shape and width (scatter) in the flux-weighted gravity-luminosity plane. The best agreement is obtained with models which include the effects of rotation and assume that the large majority, if not all the observed BSG evolve towards the red supergiant phase and only a few are evolving back from this stage. The effects of metallicity on the shape and scatter of the FGLR are small. (...)

Massive stars burn hydrogen through the CNO cycle during most of their evolution. When mixing is efficient, or when mass transfer in binary systems happens, chemically processed material is observed at the surface of O and B stars. ON stars show stronger lines of nitrogen than morphologically normal counterparts. Whether this corresponds to the presence of material processed through the CNO cycle or not is not known. Our goal is to answer this question. We perform a spectroscopic analysis of a sample of ON stars with atmosphere models. We determine the fundamental parameters as well as the He, C, N, and O surface abundances. We also measure the projected rotational velocities. We compare the properties of the ON stars to those of normal O stars. We show that ON stars are usually helium-rich. Their CNO surface abundances are fully consistent with predictions of nucleosynthesis. ON stars are more chemically evolved and rotate - on average - faster than normal O stars. Evolutionary models including rotation cannot account for the extreme enrichment observed among ON main sequence stars. Some ON stars are members of binary systems, but others are single stars as indicated by stable radial velocities. Hence, mass transfer is not a simple explanation for the observed chemical properties. We conclude that ON stars show extreme chemical enrichment at their surface, consistent with nucleosynthesis through the CNO cycle. Its origin is not clear at present.

Massive stars are key sources of radiative, kinetic, and chemical feedback in the universe. Grids of massive star models computed by different groups each using their own codes, input physics choices and numerical approximations, however, lead to inconsistent results for the same stars. We use three of these 1D codes---GENEC, KEPLER and MESA---to compute non-rotating stellar models of $15~\mathrm{M}_\odot$, $20~\mathrm{M}_\odot$, and $25~\mathrm{M}_\odot$ and compare their nucleosynthesis. We follow the evolution from the main sequence until the end of core helium burning. The GENEC and KEPLER models hold physics assumptions used in large grids of published models. The MESA code was set up to use convective core overshooting such that the CO core masses are consistent with those obtained by GENEC. For all models, full nucleosynthesis is computed using the NuGrid post-processing tool MPPNP. We find that the surface abundances predicted by the models are in reasonable agreement. In the helium core, the standard deviation of the elemental overproduction factors for Fe to Mo is less than $30\,\%$---smaller than the impact of the present nuclear physics uncertainties. For our three initial masses, the three stellar evolution codes yield consistent results. Differences in key properties of the models, e.g., helium and CO core masses and the time spent as a red supergiant, are traced back to the treatment of convection and, to a lesser extent, mass loss. The mixing processes in stars remain the key uncertainty in stellar modelling. Better constrained prescriptions are thus necessary to improve the predictive power of stellar evolution models.

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.

In this contribution we present initial results of a study on convective boundary mixing (CBM) in massive stellar models using the GENEVA stellar evolution code. Before undertaking costly 3D hydrodynamic simulations, it is important to study the general properties of convective boundaries, such as the: composition jump; pressure gradient; and `stiffness'. Models for a 15Mo star were computed. We found that for convective shells above the core, the lower (in radius or mass) boundaries are `stiffer' according to the bulk Richardson number than the relative upper (Schwarzschild) boundaries. Thus, we expect reduced CBM at the lower boundaries in comparison to the upper. This has implications on flame front propagation and the onset of novae.

Recent stellar evolution computations indicate that massive stars in the range ~ 20 - 30 Msun are located in the blue supergiant (BSG) region of the Hertzsprung-Russell diagram at two different stages of their life: immediately after the main sequence (MS, group 1) and during a blueward evolution after the red supergiant phase (group 2). From the observation of the pulsationnal properties of a subgroup of variable BSGs (alpha Cyg variables), one can deduce that these stars belongs to group 2. It is however difficult to simultaneously fit the observed surface abundances and gravity for these stars, and this allows to constrain the physical processes of chemical species transport in massive stars. We will show here that the surface abundances are extremely sensitive to the physics of convection, particularly the location of the intermediate convective shell that appears at the ignition of the hydrogen shell burning after the MS. Our results show that the use of the Ledoux criterion to determine the convective regions in the stellar models leads to a better fit of the surface abundances for alpha Cyg variables than the Schwarzschild one.

Stability boundaries of radial pulsations in massive stars are compared with positions of variable and non-variable blue-supergiants in the spectroscopic HR (sHR) diagram (Langer & Kudritzki 2014), whose vertical axis is $4\log T_{\rm eff}-\log g \,(= \log L/M)$. Observational data indicate that variables tend to have higher $L/M$ than non-variables in agreement with the theoretical prediction. However, many variable blue-supergiants are found to have values of $L/M$ below the theoretical stability boundary; i.e., surface gravities seem to be too high by around 0.2-0.3 dex.

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.

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]

Recent stellar evolution computations show that the blue supergiant (BSG) stars could come from two distinct populations: a first group arising from massive stars that just left the main sequence (MS) and are crossing the Hertzsprung-Russell diagram (HRD) towards the red supergiant (RSG) branch, and a second group coming from stars that have lost considerable amount of mass during the RSG stage and are crossing the HRD for a second time towards the blue region. Due to very different luminosity-to- mass ratio, only stars from the second group are expected to have excited pulsations observable at the surface. In a previous work, we have shown that our models were able to reproduce the pulsational properties of BSGs. However, these models failed to reproduce the surface chemical composition of stars evolving back from a RSG phase. In this paper, we show how the use of the Ledoux criterion instead of the Schwarzschild one for convection allows to significantly improve the agreement with the observed chemical composition, while keeping the agreement with the pulsation periods. This gives some support to the Ledoux criterion.

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.

Aims. Observations show nebulae around some massive stars but not around others. If observed, their chemical composition is far from homogeneous. Our goal is to put these observational features into the context of the evolution of massive stars and their circumstellar medium (CSM) and, more generally, to quantify the role of massive stars for the chemical and dynamical evolution of the ISM. Methods. Using the A-MAZE code, we perform 2d-axisymmetric hydrodynamical simulations of the evolution of the CSM, shaped by stellar winds, for a whole grid of massive stellar models from 15 to 120 Msun and following the stellar evolution from the zero-age main-sequence to the time of supernova explosion. In addition to the usual quantities, we also follow five chemical species: H, He, C, N, and O. Results. We show how various quantities evolve as a function of time: size of the bubble, position of the wind termination shock, chemical composition of the bubble, etc. The chemical composition of the bubble changes considerably compared to the initial composition, particularly during the red-supergiant (RSG) and Wolf-Rayet (WR) phases. In some extreme cases, the inner region of the bubble can be completely depleted in hydrogen and nitrogen, and is mainly composed of carbon, helium and oxygen. We argue why the bubble typically expands at a lower rate than predicted by self-similarity theory. In particular, the size of the bubble is very sensitive to the density of the ISM, decreasing by a factor of around 2.5 for each additional dex in ISM density. The bubble size also decreases with the metallicity of the central star, as low-metallicity stars have weaker winds. Our models qualitatively fit the observations of WR ejecta nebulae.

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?

(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 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]

If a massive star has lost significant mass during its red-supergiant stage, it would return to blue region in the HR diagram and spend a part of the core-He burning stage as a blue supergiant having a luminosity to mass ratio (L/M) considerably larger than about 10^4 (in solar units); the duration depends on the degree of internal mixing and on the metallicity. Then, various stellar pulsations are excited by enhanced \kappa-mechanism and strange mode instability. Assuming these pulsations to be responsible for (at least some of) the quasi-periodic light and radial-velocity variations in \alpha Cygni variables including luminous blue variables (LBVs; or S Dor variables), we can predict masses and surface compositions for these variables, and compare them with observed ones to constrain the evolutionary models. We discuss radial pulsations excited in evolutionary models of an initial mass of 40 M_\odot with solar metallicity of Z=0.014, and compare them to micro-variations in the two Galactic LBVs, HR Car and HD 160529. We have found that these stars should has lost more than half of the initial mass and their surface CNO abundances should be significantly modified from the original ones showing partial H-burning products.

A massive star can enter the blue supergiant region either evolving directly from the main-sequence, or evolving from a previous red supergiant stage. The fractions of the blue supergiants having different histories depend on the internal mixing and mass-loss during the red supergiant stage. We study the possibility to use diagnostics based on stellar pulsation to discriminate blue supergiants having different evolution histories. For this purpose we have studied the pulsation property of massive star models calculated with the Geneva stellar evolution code for initial masses ranging from 8 to 50 M$_\odot$ with a solar metallicity of $Z=0.014$. We have found that radial pulsations are excited in the blue-supergiant region only in the models that had been red-supergiants before. This would provide us with a useful mean to diagnose the history of evolution of each blue-supergiant. At a given effective temperature, much more nonradial pulsations are excited in the model after the red-supergiant stage than in the model evolving towards the red-supergiant. The properties of radial and nonradial pulsations in blue supergiants are discussed. Predicted periods are compared with period ranges observed in some \alpha-Cygni variables in the Galaxy and NGC 300. We have found that blue supergiant models after the red-supergiant stage roughly agree with observed period ranges in most cases. However, we are left with the puzzle that the predicted surface N/C and N/O ratios seem to be too high compared with those of Deneb and Rigel.

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.

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.

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.

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

Wolf-Rayet stars are evolved massive stars, and the relative number of WC-type and WN-type WRs should vary with the metallicity of the host galaxy, providing a sensitive test of stellar evolutionary theory. However, past studies of the WR content of M31 have been biased towards detecting WC stars, as their emission line signatures are much stronger than those of WNs. Here we present the results of a survey covering all of M31's optical disk (2.2 deg^2), with sufficient sensitivity to detect the weaker-lined WN-types. We identify 107 newly found WR stars, mostly of WN-type. This brings the total number of spectroscopically confirmed WRs in M31 to 154, a number we argue is complete to ~95%, except in regions of unusually high reddening. This number is consistent with what we expect from the integrated Halpha luminosity compared to that of M33. The majority of these WRs formed in OB associations around the Population I ring, although 5% are truly isolated. Both the relative number of WC to WN-type stars as well as the WC subtype distribution suggest that most WRs exist in environments with higher-than-solar metallicities, which is consistent with studies of M31's metallicity. Although the WC-to-WN ratio we find for M31 is much lower than that found by previous studies, it is still higher than what the Geneva evolutionary models predict. This may suggest that Roche-lobe overflow produces the excess of WC stars observed at high metallicity, or that the assumed rotational velocities in the models are too high.

The research of exoplanets has entered an era in which we characterize extrasolar planets. This has become possible with measurements of radii and luminosities. Meanwhile, radial velocity surveys discover also very low-mass planets. Uniting all this observational data into one coherent picture to better understand planet formation is an important, but difficult undertaking. Our approach is to develop a model which can make testable predictions for all these observational methods. We continue to describe how we have extended our formation model into a self-consistently coupled formation and evolution model. We show how we calculate the internal structure of the solid core and radiogenic heating. We also improve the protoplanetary disk model. Finally, we conduct population synthesis calculations. We present how the planetary mass-radius relationship of planets with primordial H/He envelopes forms and evolves in time. The basic shape of the M-R relation can be understood from the core accretion model. Low-mass planets cannot bind massive envelopes, while super-critical cores necessarily trigger runway gas accretion, leading to "forbidden" zones in the M-R plane. For a given mass, there is a considerable diversity of radii. We compare the synthetic M-R relation with the observed one, finding good agreement for a>0.1 AU. The synthetic radius distribution is characterized by a strong increase towards small R, and a second, lower local maximum at ~1 Jovian radius. The increase towards small radii reflects the increase of the mass function towards low M. The second local maximum is due to the fact that radii are nearly independent of mass for giant planets. A comparison of the synthetic radius distribution with Kepler data shows agreement for R>2 Earth radii, but divergence for smaller radii. We predict that in the next few years, Kepler should find the second, local maximum at ~1 Jovian radius.

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.

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.

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 increasing observed number of supernova events allows for finding ever more frequently the progenitor star in archive images. In a few cases, the progenitor star is a yellow supergiant star. The estimated position in the Hertzsprung-Russell diagram of these stars is not compatible with the theoretical tracks of classical single star models. According to several authors, the mass-loss rates during the red supergiant phase could be underestimated. We study the impact of an increase of these mass-loss rates on the position of 12 to 15 M\odot stars at the end of their nuclear life, in order to reconcile the theoretical tracks with the observed yellow supergiant progenitors. We perform calculations of 12 to 15 M\odot rotating stellar models using the Geneva stellar evolution code. To account for the uncertainties in the mass-loss rates during the RSG phase, we increase the mass-loss rate of the star (between 3 and 10 times the standard one) during that phase and compare the evolution of stars undergoing such high mass-loss rates with models computed with the standard mass-loss prescription. We show that the final position of the models in the Hertzsprung-Russell diagram depends on the mass loss they undergo during the red supergiant phase. With an increased mass-loss rate, we find that some models end their nuclear life at positions that are compatible with the observed position of several supernova progenitors. We conclude that an increased mass-loss rate (whom physical mechanism still need to be clarified) allows single star models to reproduce simultaneously the estimated position in the HRD of the YSG SN progenitors, as well as the SN type.

[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 medium around massive stars is strongly shaped by the stellar winds. Those winds depend on various stellar parameters (effective temperature, luminosity, chemical composition, rotation, ...), which are varying as a function of the time. Using the wind properties obtained by classical stellar evolution code allows for the multi- D hydrodynamical simulation of the circum-stellar medium accounting for the time variations of the wind. We present here the preliminary results of the simulation of the medium around a fast rotating star. Comparing the results of such simulations with the observed properties of the circum-stellar medium will allow to better understand the interactions of the stellar winds with the interstellar medium, and could give hints on the past history of the mass loss around massive stars, providing constraints for stellar models.

We discuss, in the context of the single star scenario, the nature of the progenitors of Red Supergiants (RSG), of Luminous Blue Variables (LBV) and of Wolf-Rayet (WR) stars. These three different populations correspond to evolved phases of Main Sequence (MS) OB stars. Axial rotation and mass loss have a great influence on massive star evolution in general and more specifically on the durations of these different phases. Moderate rotation and mass loss, during the MS phase, favor the evolution towards the RSG stage. Fast rotation and strong mass loss during the MS phase, in contrast, prevent the star from becoming a RSG and allow the star to pass directly from the OB star phase into the WR phase. Mass loss during the RSG stage may make the star evolve back in the blue part of the HR diagram. We argue that such an evolution may be more common than presently accounted for in stellar models. This might be the reason for the lack of type IIP SNe with RSG progenitors having initial masses between 18 and 30 M$_\odot$. The LBVs do appear as a possible transition phase between O and WR stars or between WNL and WNE stars. Fast rotation and/or strong mass loss during the Main-Sequence phase prevent the formation of LBV stars. The mechanisms driving the very strong ejections shown by LBV stars are still unknown. (abridged)

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.

Whenever stars are rotating very fast (Omega/Omega_crit > 0.7, with Omega_crit the Keplerian angular velocity of the star accounting for its deformation) radiative stellar winds are enhanced in polar regions. This theoretical prediction is now confirmed by interferometric observations of fast rotating stars.} Polar winds remove less angular momentum than spherical winds and thus allow the star to keep more angular momentum. We quantitatively assess the importance of this effect. First we use a semi-analytical approach to estimate the variation of the angular momentum loss when the rotation parameter increases. Then we compute complete 9 MSun stellar models at very high angular velocities (starting on the ZAMS with Omega/Omega_crit = 0.8 and reaching the critical velocity during the Main Sequence) with and without radiative wind anisotropies. When wind anisotropies are accounted for, the angular momentum loss rate is reduced by less than 4% for Omega/Omega_crit < 0.9 with respect to the case of spherical winds. The reduction amounts to at most 30% when the star is rotating near the critical velocity. These values result from two counteracting effects: on the one hand polar winds reduce the loss of angular momentum, on the other hand, surface deformations imply that the mass which is lost at high co-latitude is lost at a larger distance from the rotational axis and thus removes more angular momentum. In contrast with previous studies, which neglected surface deformations, we show that the radiative wind anisotropies have a relatively modest effect on the evolution of the angular momentum content of fast rotating stars.

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.

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.

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

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.

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.

AGB stars have long been held responsible for the important star-to-star variations in light elements observed in Galactic globular clusters. We analyse the main impacts of a first generation of rotating intermediate-mass stars on the chemical properties of second-generation globular cluster stars. The rotating models were computed without magnetic fields and without the effects of internal gravity waves. They account for the transports by meridional currents and turbulence. We computed the evolution of both standard and rotating stellar models with initial masses between 2.5 and 8 Msun within the metallicity range covered by Galactic globular clusters. During central He-burning, rotational mixing transports fresh CO-rich material from the core towards the hydrogen-burning shell, leading to the production of primary 14N. In stars more massive than M > 4 Msun, the convective envelope reaches this reservoir during the second dredge-up episode, resulting in a large increase in the total C+N+O content at the stellar surface and in the stellar wind. The corresponding pollution depends on the initial metallicity. At low- and intermediate-metallicity, it is at odds with the constancy of C+N+O observed among globular cluster low-mass stars. With the given input physics, our models suggest that massive rotating AGB stars have not shaped the abundance patterns observed in low- and intermediate-metallicity globular clusters. Our non-rotating models, on the other hands, do not predict surface C+N+O enhancements, hence are in a better position as sources of the chemical anomalies in globular clusters showing the constancy of the C+N+O. However at the moment, there is no reason to think that intermediate mass stars were not rotating.

We discuss the properties of the progenitors of core collapse SNe, as they can be deduced from rotating stellar models of single stars. The type of the SN progenitor was determined from the surface abundances at the pre-SN stage. The type of the SN event was obtained from the masses of hydrogen and helium ejected at the time of the core-collapse SN event. We find that the minimum amount of helium ejected by a core-collapse SN (of whatever type) is around 0.3 MSun. There is no difference between the WC and WO stars in the ejected masses of helium, CNO elements, and heavy elements. Also no difference is expected between the chemical composition of a WC star resulting from a normal or a homogeneous evolution. The progenitors of type Ib SNe are WNL, WNE, or less massive WC stars. Those of type Ic are WC and WO stars. WO stars are produced in a limited mass range (around 60 MSun) and only at low metallicity (for <= 0.010) as already found. The WO stars are the progenitors of only a small fraction of type Ic. Present stellar models indicate that, at solar metallicity, there is about 1 type Ib SN for 1 type Ic, and this ratio rises to 3 type Ic for 1 type Ib SN at twice solar metallicity. At this metallicity, type Ic's are more frequent than type Ib's because most massive stars that go through a WNE stage evolve further into a WC/WO phase. Current models can account for the observed number ratios SN Ib/SN II and SN Ic/SN II and for their observed variation with the metallicity. In case no supernova occurs when a black hole is formed, single-star models can still account for more than half of the observed (SN Ib+SN Ic)/SN II ratio for Z >= ZSun. For the GRB rate, our models produce too large a number for such an event, even if we restrict the progenitor to the WO stars.

The hot Neptune orbiting around GJ 436 is a unique example of an intermediate mass planet. Its close-in orbit suggests that the planet has undergone migration and its study is fundamental to understanding planet formation and evolution. As it transits its parent star, it is the only Neptune-mass extrasolar planet of known mass and radius, being slightly larger and more massive than Neptune (M=22.6 M_Earth, R=4.19R_Earth). In this regime, several bulk compositions are possible: from an Earth-like core with a thick hydrogen envelope to a water-rich planet with a thin hydrogen envelope comprising a Neptune-like structure. We combine planet-structure modeling with an advanced planet-formation model to assess the likelihood of the different possible bulk compositions of GJ 436 b. We find that both an envelope-free water planet ("Ocean planet") as well as a diminute version of a gaseous giant planet are excluded. Consisting of a rocky core with a thick hydrogen/helium envelope, a "dry" composition produces not only too small a radius but is also a very unlikely outcome of planet formation around such a low-mass star. We conclude that GJ 436 b is probably of much higher rock content than Neptune (more than 45% in mass), with a small H-He envelope (10 - 20% in mass). This is the expected outcome of the gathering of materials during the migration process in the inner disk, creating a population of which the hot Neptune is representative.

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.

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.

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 is a determinant factor which strongly affects the evolution and the fate of massive stars. At low metallicity, stars are supposed to rotate faster than at the solar one. This favors the existence of stars near the critical velocity. In this rotation regime, the deformation of the stellar surface becomes important, and wind anisotropy develops. Polar winds are expected to be dominant for fast rotating hot stars. These polar winds allow the star to lose large quantities of mass and still retain a high angular momentum, and they modifie the evolution of the surface velocity and the final angular momentum kept in the star's core. We show here how these winds affect the final stages of massive stars, according to our knowledge about Gamma Ray Bursts. Computation of theoretical Gamma Ray Bursts rate indicates that our models have too fast rotating cores, and that we need to include an additional effect to spin them down. Magnetic fields in stars act in this direction, and we show how they modify the evolution of massive star up to the final stages.

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.

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.

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

We study the effects of rotation on the outer convective zones of massive stars. We examine the effects of rotation on the thermal gradient and on the Solberg--Hoiland term by analytical developments and by numerical models. Writing the criterion for convection in rotating envelopes, we show that the effects of rotation on the thermal gradient are much larger and of opposite sign to the effect of the Solberg-Hoiland criterion. On the whole, rotation favors convection in stellar envelopes at the equator and to a smaller extent at the poles. In a rotating 20 Msun star at 94% of the critical angular velocity, there are two convective envelopes, with the bigger one having a thickness of 13.2% of the equatorial radius. In the non-rotating model, the corresponding convective zone has a thickness of only 4.6% of the radius. The occurrence of outer convection in massive stars has many consequences.

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.