The apparent lack of massive O-type stars near the zero-age main sequence (at ages < 2 Myr) is a topic widely discussed. Different explanations for this elusive detection have been proposed, but no firm conclusions have been reached yet. We reassess this empirical result benefiting from the high-quality spectroscopic observations of >400 Galactic O-type stars gathered by the IACOB and OWN surveys. We used temperatures and gravities from a iacob-gbat/fastwind spectroscopic analysis to locate our sample in the Kiel and spectroscopic HR diagrams. We evaluated the completeness of our sample of stars, observational biases using information from the Galactic O star catalog (GOSC), systematics of our methodology, and compare with other recent studies using smaller samples of Galactic O-type stars. We base our discussion on the spectroscopic HR diagram to avoid the use of uncertain distances. We performed a detailed study of the young cluster Trumpler-14 as an example of how Gaia cluster distances can help to construct the associated classical HR diagram. The apparent lack of massive O-type stars near the ZAMS with masses between 30 and 70 Msol persist even when spectroscopic results from a large, non-biased sample of stars are used. We do not find correlation between the dearth of stars and observational biases, limitations of our methodology, or the use of spectroscopic HR diagram instead of the classical one. Investigating the efficiency of mass accretion during the formation process we conclude that an adjustment of the accretion rate towards lower values could reconcile the hotter boundary of detected O-type stars and the theoretical birthline. Last, we discuss that the presence of a small sample of O2-O3.5 stars found closer to the ZAMS might be explained taking into account non-standard star evolution (e.g. binary interaction, mergers, or homogeneous evolution).

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.

Context: Asteroseismic measurements reveal that an unknown efficient angular momentum (AM) transport mechanism is needed for subgiant and red giant stars. A revised prescription for AM transport by the magnetic Tayler instability has been recently proposed as a possible candidate for such a missing mechanism. Results: The revised prescription for the transport by the Tayler instability leads to low core rotation rates after the main sequence that are in better global agreement with asteroseismic measurements than those predicted by models with purely hydrodynamic processes or with the original Tayler-Spruit dynamo. A detailed comparison with asteroseismic data shows that the rotational properties of at most two of the six subgiants can be correctly reproduced by models accounting for this revised magnetic transport process. This result is obtained independently of the value adopted for the calibration parameter in this prescription. We also find that this transport by the Tayler instability faces difficulties in simultaneously reproducing asteroseismic measurements available for subgiant and red giant stars. The low values of the calibration parameter needed to correctly reproduce the rotational properties of two of the six subgiants lead to core rotation rates during the red giant phase that are too high. Inversely, the higher values of this parameter needed to reproduce the core rotation rates of red giants lead to a very low degree of radial differential rotation before the red giant phase, which is in contradiction with the internal rotation of subgiant stars. Conclusions: In its present form, the revised prescription for the transport by the Tayler instability does not provide a complete solution to the missing AM transport revealed by asteroseismology of evolved stars.

We show that it is not possible to determine the final mass $M_{\rm fin}$ of a red supergiant (RSG) at the pre-supernova (SN) stage from its luminosity $L$ and effective temperature $T_{\rm eff}$ alone. Using a grid of stellar models, we demonstrate that for a given value of $L$ and $T_{\rm eff}$, a RSG can have a range of $M_{\rm fin}$ as wide as 3 to $45~\mathrm{M}_{\odot}$. While the probability distribution within these limits is not flat, any individual determination of $M_{\rm fin}$ for a RSG will be degenerate. This makes it difficult to determine its evolutionary history and to map $M_{\rm fin}$ to an initial mass. Single stars produce a narrower range that is difficult to accurately determine without making strong assumptions about mass loss, convection, and rotation. Binaries would produce a wider range of RSG $M_{\rm fin}$. However, the final Helium core mass M$_{\rm He-core}$ is well determined by the final luminosity and we find $\log (\mathrm{M}_{\rm He-core}/M_{\odot}) = 0.659 \log (L/\mathrm{L}_{\odot}) -2.630$ Using this relationship, we derive M$_{\rm He-core}$ for directly imaged SN progenitors and one failed SN candidate. The value of $M_{\rm fin}$ for stripped star progenitors of SNe IIb is better constrained by $L$ and $T_{\rm eff}$ due to the dependence of $T_{\rm eff}$ on the envelope mass $M_{\rm env}$ for $M_{\rm env} \lesssim 1~$M$_{\odot}$. Given the initial mass function, our results apply to the majority of progenitors of core collapse SNe, failed SNe and direct collapse black holes.

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 investigate the effects of mass loss during the main-sequence (MS) and post-MS phases on the final black hole (BH) masses of massive stars. We compute solar metallicity Geneva stellar evolution models of an 85 $M_\odot$ star with mass-loss rate ($M_\odot$) prescriptions for MS and post-MS phases. Such models could lead to massive BHs such as the recently detected 70 $M_\odot$ BH in the LB-1 system. Based on the observational constraints for $M_\odot$ of luminous stars, we discuss two possible scenarios that could produce such a massive BH at high metallicity. First, if the progenitor of LB-1 evolved from the observed population of WNh stars, we show that its average mass loss rate during its post-MS evolution was less than $1\,\times10^{-5}\,M_\odot/$yr. However, this is lower than the typical observed mass-loss rates of LBVs. A second possibility is that the progenitor evolved from a yet undetected population of 80-85~\msun\ stars with strong surface magnetic fields, which could quench mass loss during the MS evolution. In this case, the average mass-loss rate during the post-MS luminous blue variable (LBV) phase has to be less than $5\,\times10^{-5}\,M_\odot/$yr. This value is still low, considering that LBVs such as AG Carinae typically have average mass-loss rates from quiescent stellar winds close or above this level. LBVs directly collapsing to massive BHs are apparently at odds with the evidence that LBVs are the direct progenitors of some supernovae (SNe). To reconcile this, we suggest that LBVs from single stars (or mergers that fully mixed) have large cores and form BHs, while binary LBVs may have smaller cores and larger envelopes (either mergers or mass gainers that do not fully mix) and could produce a SN.

Major mergers of gas-rich galaxies provide promising conditions for the formation of supermassive black holes (SMBHs; $\gtrsim10^5$ M$_\odot$) by direct collapse because they can trigger mass inflows as high as $10^4-10^5$ M$_\odot$ yr$^{-1}$ on sub-parsec scales. However, the channel of SMBH formation in this case, either dark collapse (direct collapse without prior stellar phase) or supermassive star (SMS; $\gtrsim10^4$ M$_\odot$), remains unknown. Here, we investigate the limit in accretion rate up to which stars can maintain hydrostatic equilibrium. We compute hydrostatic models of SMSs accreting at $1-1000$ M$_\odot$ yr$^{-1}$, and estimate the departures from equilibrium a posteriori by taking into account the finite speed of sound. We find that stars accreting above the atomic cooling limit ($\gtrsim10$ M$_\odot$ yr$^{-1}$) can only maintain hydrostatic equilibrium once they are supermassive. In this case, they evolve adiabatically with a hylotropic structure, that is, entropy is locally conserved and scales with the square root of the mass coordinate. Our results imply that stars can only become supermassive by accretion at the rates of atomically cooled haloes ($\sim0.1-10$ M$_\odot$ yr$^{-1}$). Once they are supermassive, larger rates are possible.

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.

In the last decade, the Kepler and CoRoT space-photometry missions have demonstrated the potential of asteroseismology as a novel, versatile and powerful tool to perform exquisite tests of stellar physics, and to enable precise and accurate characterisations of stellar properties, with impact on both exoplanetary and Galactic astrophysics. Based on our improved understanding of the strengths and limitations of such a tool, we argue for a new small/medium space mission dedicated to gathering high-precision, high-cadence, long photometric series in dense stellar fields. Such a mission will lead to breakthroughs in stellar astrophysics, especially in the metal poor regime, will elucidate the evolution and formation of open and globular clusters, and aid our understanding of the assembly history and chemodynamics of the Milky Way's bulge and few nearby dwarf galaxies.

Current burning issues in stellar physics, for both hot and cool stars, concern their magnetism. In hot stars, stable magnetic fields of fossil origin impact their stellar structure and circumstellar environment, with a likely major role in stellar evolution. However, this role is complex and thus poorly understood as of today. It needs to be quantified with high-resolution UV spectropolarimetric measurements. In cool stars, UV spectropolarimetry would provide access to the structure and magnetic field of the very dynamic upper stellar atmosphere, providing key data for new progress to be made on the role of magnetic fields in heating the upper atmospheres, launching stellar winds, and more generally in the interaction of cool stars with their environment (circumstellar disk, planets) along their whole evolution. UV spectropolarimetry is proposed on missions of various sizes and scopes, from POLLUX on the 15-m telescope LUVOIR to the Arago M-size mission dedicated to UV spectropolarimetry.

Over forty years of research suggests that the common envelope phase, in which an evolved star engulfs its companion upon expansion, is the critical evolutionary stage forming short-period, compact-object binary systems, such as coalescing double compact objects, X-ray binaries, and cataclysmic variables. In this work, we adapt the one-dimensional hydrodynamic stellar evolution code, MESA, to model the inspiral of a 1.4M$_{\odot}$ neutron star (NS) inside the envelope of a 12$M_{\odot}$ red supergiant star. We self-consistently calculate the drag force experienced by the NS as well as the back-reaction onto the expanding envelope as the NS spirals in. Nearly all of the hydrogen envelope escapes, expanding to large radii ($\sim$10$^2$ AU) where it forms an optically thick envelope with temperatures low enough that dust formation occurs. We simulate the NS orbit until only 0.8M$_{\odot}$ of the hydrogen envelope remains around the giant star's core. Our results suggest that the inspiral will continue until another $\approx$0.3M$_{\odot}$ are removed, at which point the remaining envelope will retract. Upon separation, a phase of dynamically stable mass transfer onto the NS accretor is likely to ensue, which may be observable as an ultraluminous X-ray source. The resulting binary, comprised of a detached 2.6M$_{\odot}$ helium-star and a NS with a separation of 3.3-5.7R$_{\odot}$, is expected to evolve into a merging double neutron-star, analogous to those recently detected by LIGO/Virgo. For our chosen combination of binary parameters, our estimated final separation (including the phase of stable mass transfer) suggests a very high $\alpha_{\rm CE}$-equivalent efficiency of $\approx$5.

Context. The CoRoT and Kepler missions have paved the way for synergies between exoplanetology and asteroseismology. The use of seismic data helps providing stringent constraints on the stellar properties which directly impact the results of planetary studies. Amongst the most interesting planetary systems discovered by Kepler, Kepler-444 is unique by the quality of its seismic and classical stellar constraints. Its magnitude, age and the presence of 5 small-sized planets orbiting this target makes it an exceptional testbed for exoplanetology. Aims. We aim at providing a detailed characterization of Kepler-444, focusing on the dependency of the results on variations of key ingredients of the theoretical stellar models. This thorough study will serve as a basis for future investigations of the planetary evolution of the system orbiting Kepler-444. Methods. We use local and global minimization techniques to study the internal structure of the exoplanet-host star Kepler-444. We combine seismic observations from the Kepler mission, Gaia DR2 data and revised spectroscopic parameters to precisely constrain its internal structure and evolution. Results. We provide updated robust and precise determinations of the fundamental parameters of Kepler-444 and demonstrate that this low-mass star bore a convective core during a significant portion of its life on the main-sequence. Using seismic data, we are able to estimate the lifetime of the convective core to approximately 8 Gyr out of the 11 Gyr of the evolution of Kepler-444. The revised stellar parameters found by our thorough study are M = 0.754 +- 0.03M_Sun , R = 0.753 +- 0.01R_Sun , Age = 11 +- 1 Gy.

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.

Supermassive stars (SMSs) are candidates for being progenitors of supermassive quasars at high redshifts. However, their formation process requires strong mechanisms that would be able to extract the angular momentum of the gas that the SMSs accrete. We investigate under which conditions the magnetic coupling between an accreting SMS and its winds can remove enough angular momentum for accretion to proceed from a Keplerian disc. We numerically computed the rotational properties of accreting SMSs that rotate at the $\Omega\Gamma$-limit and estimated the magnetic field that is required to maintain the rotation velocity at this limit using prescriptions from magnetohydrodynamical simulations of stellar winds. We find that a magnetic field of 10 kG at the stellar surface is required to satisfy the constraints on stellar rotation from the $\Omega\Gamma$-limit. Magnetic coupling between the envelope of SMSs and their winds could allow for SMS formation by accretion from a Keplerian disc, provided the magnetic field is at the upper end of present-day observed stellar fields. Such fields are consistent with primordial origins.

The Sun is the most studied of stars and a laboratory of fundamental physics. However, the understanding of our star is stained by the solar modelling problem which can stem from various causes. We combine inversions of sound speed, an entropy proxy and the Ledoux discriminant with the position of the base of the convective zone and the photospheric helium abundance to test combinations of ingredients such as equation of state, abundance and opacity tables. We study the potential of the inversions to constrain ad-hoc opacity modifications and additional mixing in the Sun. We show that they provide constraints on these modifications to the ingredients and that the solar problem likely occurs from various sources and using phase shifts with our approach is the next step to take.

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.

Gamma-ray astronomy explores the most energetic photons in nature to address some of the most pressing puzzles in contemporary astrophysics. It encompasses a wide range of objects and phenomena: stars, supernovae, novae, neutron stars, stellar-mass black holes, nucleosynthesis, the interstellar medium, cosmic rays and relativistic-particle acceleration, and the evolution of galaxies. MeV gamma-rays provide a unique probe of nuclear processes in astronomy, directly measuring radioactive decay, nuclear de-excitation, and positron annihilation. The substantial information carried by gamma-ray photons allows us to see deeper into these objects, the bulk of the power is often emitted at gamma-ray energies, and radioactivity provides a natural physical clock that adds unique information. New science will be driven by time-domain population studies at gamma-ray energies. This science is enabled by next-generation gamma-ray instruments with one to two orders of magnitude better sensitivity, larger sky coverage, and faster cadence than all previous gamma-ray instruments. This transformative capability permits: (a) the accurate identification of the gamma-ray emitting objects and correlations with observations taken at other wavelengths and with other messengers; (b) construction of new gamma-ray maps of the Milky Way and other nearby galaxies where extended regions are distinguished from point sources; and (c) considerable serendipitous science of scarce events -- nearby neutron star mergers, for example. Advances in technology push the performance of new gamma-ray instruments to address a wide set of astrophysical questions.

Black hole (BH) spins in low-mass X-ray binaries (LMXBs) cover a range of values that can be explained by accretion after BH birth. In contrast, the three BH spin measurements in high-mass X-ray binaries (HMXBs) show only values near the maximum and likely have a different origin connected to the BH stellar progenitor. We explore here two possible scenarios to explain the high spins of BHs in HMXBs: formation in binaries that undergo mass transfer (MT) during the main sequence (MS; Case-A MT), and very close binaries undergoing chemically homogeneous evolution (CHE). We find that both scenarios are able to produce high-spin BHs if internal angular momentum (AM) transport in the progenitor star after its MS evolution is not too strong (i.e., weak coupling between the stellar core and its envelope). If instead efficient AM transport is assumed, we find that the resulting BH spins are always too low with respect to observations. The Case-A MT model provides a good fit for the BH spins, the masses of the two components, and the final orbital periods for two of the three BHs in HMXBs with measured spins. For one of them, the mass predicted for the BH companion is significantly lower than observed, but this depends strongly on the assumed efficiency of mass transfer. The CHE models predict orbital periods that are too large for all three cases considered here. We expect the Case-A MT to be much more frequent at the metallicities involved, so we conclude that the Case-A MT scenario is preferred. Finally, we predict that the stellar companions of HMXBs formed through the Case-A MT have enhanced nitrogen surface abundances, which can be tested by future observations.

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

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

The Sun is the most studied of all stars. It is a reference for all other observed stars and a laboratory of fundamental physics helping us understand processes occuring in conditions irreproducible on Earth. However, our understanding of the Sun is currently stained by the solar modelling problem which can stem from various causes, such as the opacities, the equation of state and the mixing of chemical elements. We combine inversions of sound speed, an entropy proxy and the Ledoux discriminant with constraints such as the position of the base of the convective zone and the photospheric helium abundance. We test various combinations of standard ingredients for solar modelling such as abundance tables, equation of state, formalism for convection and diffusion and opacity tables and study the diagnostic potential of the inversions to constrain ad-hoc modifications of the opacity profile and additional mixing below the convective envelope. Combining inversions provides stringent constraints on the modifications on the models, far beyond what is achievable only from sound speed inversions. We constrain the form and amplitude of the opacity increase required and show that a 15% increase at log T=6.35 provides a significant improvement but that a more global increase of the opacity, within the uncertainties of current tables, coupled with an additional mixing at the bottom of the convective zone gives the best agreement for low metallicity models. We show that high metallicity models do not satisfy all the inversion results. We conclude that the solar problem likely occurs from various small sources, as ingredients such as the equation of state or the formalism of convection can induce small but significant changes and that using phase shift analyses combined with our approach is the next step for a better understanding of the inaccuracies of solar models just below the convective envelope.

We explore the initial conditions that lead to detached binary evolution of massive pre-MS binaries and ask how large a fraction of the observed binary systems may have been initially formed as low-mass protobinaries and later undergone a significant accretion phase while remaining detached. We develop a family of analytic models to describe the orbital separation, $a$, and mass ratio, $q$, evolution. For a given mass accretion rate onto the binary system, we define a recipe for distributing this mass between the two components. For this we introduce a parameter $\eta$, such that $\frac{\dot{M}_2}{\dot{M}_1}=q^\eta$ at any time, to determine the binary mass ratio evolution. Depending on the choice of $\eta$, any type of mass ratio evolution is possible. Furthermore, we use MESA, a detailed stellar structure code, to calculate an extensive grid of binary sequences where a protobinary undergoes accretion, and we identify the initial conditions that separates detached from non-detached pre-MS binary evolution. A value of $\eta$ around 2 allows accretion growth in detached systems to form close massive binaries on the Zero Age Main-Sequence with minimum orbital periods down to about 1.2 days for $M_{\rm 1, ZAMS}=20-30M_{\odot}$ twin-binaries. $\eta=2$ can also reproduce the observed population of binary systems with primary stars above $6M_{\odot}$. The whole observed range of massive close binaries can form via accretion growth in detached systems, making the binary formation channel of accretion growth a strong contender to explain the formation of massive close binaries, including progenitors of coalescing binary black holes.

Recent studies show that rotation significantly affects the s-process in massive stars. We provide tables of yields for non-rotating and rotating massive stars between 10 and 150 $M_{\odot}$ at $Z=10^{-3}$ ([Fe/H] $=-1.8$). Tables for different mass cuts are provided. The complete s-process is followed during the whole evolution with a network of 737 isotopes, from Hydrogen to Polonium. A grid of stellar models with initial masses of 10, 15, 20, 25, 40, 60, 85, 120 and 150 $M_{\odot}$ and with an initial rotation rate of both 0 or 40$~\%$ of the critical velocity was computed. Three extra models were computed in order to investigate the effect of faster rotation (70$~\%$ of the critical velocity) and of a lower $^{17}$O($\alpha,\gamma$) reaction rate. At the considered metallicity, rotation has a strong impact on the production of s-elements for initial masses between 20 and 60 $M_{\odot}$. In this range, the first s-process peak is boosted by $2-3$ dex if rotation is included. Above 60 $M_{\odot}$, s-element yields of rotating and non-rotating models are similar. Increasing the initial rotation from 40$~\%$ to 70$~\%$ of the critical velocity enhances the production of $40 \lesssim Z \lesssim 60$ elements by $\sim 0.5-1$ dex. Adopting a reasonably lower $^{17}$O($\alpha,\gamma$) rate in the fast rotating model (70$~\%$ of the critical velocity) boosts again the yields of s-elements with $55 \lesssim Z \lesssim 82$ by about 1 dex. In particular, a modest amount of Pb is produced. Together with s-elements, some light elements (particularly fluorine) are strongly overproduced in rotating models.

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

Various binary black hole formation channels have been proposed since the first gravitational event GW150914 was discovered by the Advanced Laser Interferometer Gravitational-Wave Observatory (AdLIGO). For all evolutionary channels based on the evolution of isolated binaries, the immediate progenitor of the binary black hole is a close binary system composed of a black hole and a helium star. We perform detailed binary evolution and systematically explore the parameter space of initial binary properties, including initial black hole and helium star masses, initial rotation of the helium star as well as metallicity. We argue that the spin of the first-born black hole at its birth is negligible ($\lesssim 0.1$), hence the second-born black hole's spin dominates the measured effective spin, $\chi_{\rm eff}$, from gravitational wave events of double black hole mergers. We find that tides can be important only when orbital periods are shorter than 2 days. Upon core collapse, the helium star produces a black hole (the second-born black hole in the system) with a spin that can span the entire range from zero to maximally spinning. We show that the bimodal distribution of the spin of the second-born black hole obtained in recent papers is mainly due to oversimplifying assumptions. We find an anti-correlation between the merging timescale of the two black holes, T$_{\rm merger}$, and the effective spin $\chi_{\rm eff}$. Finally, we provide new prescriptions for the tidal coefficient E$_2$ for both H-rich and the helium-rich stars. We predict that, with future improvements to AdLIGO's sensitivity, the sample of merging binary black hole systems will show an overdensity of sources with positive but small $\chi_{\rm eff}$. originating from lower mass black hole mergers born at low redshift.

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.

The present paper discusses the main physical effects produced by stellar rotation on presupernovae, as well as observations which confirm these effects and their consequences for presupernova models. Rotation critically influences the mass of the exploding cores, the mass and chemical composition of the envelopes and the types of supernovae, as well as the properties of the remnants and the chemical yields. In the formation of gamma-ray bursts, rotation and the properties of rotating stars appear as the key factor. In binaries, the interaction between axial rotation and tidal effects often leads to interesting and unexpected results. Rotation plays a key role in shaping the evolution and nucleosynthesis in massive stars with very low metallicities (metallicity below about the Small Magellanic Cloud metallicity down to Population III stars). At solar and higher metallicities, the effects of rotation compete with those of stellar winds. In close binaries, the synchronisation process can lock the star at a high rotation rate despite strong mass loss and thus both effects, rotation and stellar winds, have a strong impact. In conclusion, rotation is a key physical ingredient of the stellar models and of presupernova stages, and the evolution both of single stars and close binaries. Moreover, important effects are expected along the whole cosmic history.

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.

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

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

Following the renewed interest for axions as a dark matter component, we revisit the effects of energy loss by axion emission on the evolution of the first generation of stars. These stars with zero metallicity are supposed to be massive, more compact, and hotter than subsequent generations. It is hence important to extend previous studies restricted to solar metallicity stars. Our analysis first compares the evolution of solar metallicity 8, 10 and 12 M$_\odot$ stars to previous work. We then calculate the evolution of 8 zero metallicity stars with and without axion losses and with masses ranging from 20 to 150 M$_\odot$. For the solar metallicity models, we confirm the disappearance of the blue loop phase for a value of the axion-photon coupling, $g_{a\gamma}=10^{-10} {\rm GeV}^{-1}$. Regarding the Pop. III stars, we show that for $g_{a\gamma}=10^{-10} {\rm GeV}^{-1}$, their evolution is not much affected by axion losses except within the range of masses 80$-$130 M$_\odot$. Such stars show significant differences both in their tracks within the $T_c$-$\rho_c$ diagram and in their central composition (in particular $^{20}$Ne and $^{24}$Mg). We discuss the origin of these modifications from the stellar physics point of view, and their potential observational signatures.

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.

Context. Globular clusters host stars with chemical peculiarities. The associated helium enrichment is expected to affect the evolution of stars, in general, and of low-mass stars, and in particular the progenitors of white dwarfs (WDs). Aims. We investigate the effects of different initial helium contents on the properties of white dwarfs such as their masses, compositions, and the time since their formation. Results. In a globular cluster at a given age, the He enrichment implies lower initial masses for stars at a given phase. Thus it leads to a decrease of the masses of white dwarfs reaching the cooling sequence. In addition the He enrichment increases the total mass and number of white dwarfs and eventually allows the presence of He white dwarf from single progenitors. Conclusions. The low He enrichment determined in most globular clusters with different methods results in negligible effects on the white dwarf properties. However, in the few globular clusters that display a high He enrichment, this may significantly affect the characteristics of the white dwarfs. In NGC 2808 and {\omega} Centauri the high He enrichment even leads to the formation of He white dwarfs from single He-rich progenitors. Therefore investigating the white dwarf mass domain in globular clusters with a high He enrichment would provide an additional indirect way to measure and constrain the He enrichment degree.

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.

Rotation plays a key role in the star-formation process, from pre-stellar cores to pre-main-sequence (PMS) objects. Understanding the formation of massive stars requires taking into account the accretion of angular momentum during their PMS phase. We study the PMS evolution of objects destined to become massive stars by accretion, focusing on the links between the physical conditions of the environment and the rotational properties of young stars. In particular, we look at the physical conditions that allow the production of massive stars by accretion. We present PMS models computed with a new version of the Geneva Stellar Evolution code self-consistently including accretion and rotation according to various accretion scenarios for mass and angular momentum. We describe the internal distribution of angular momentum in PMS stars accreting at high rates and we show how the various physical conditions impact their internal structures, evolutionary tracks, and rotation velocities during the PMS and the early main sequence. We find that the smooth angular momentum accretion considered in previous studies leads to an angular momentum barrier and does not allow the formation of massive stars by accretion. A braking mechanism is needed in order to circumvent this angular momentum barrier. This mechanism has to be efficient enough to remove more than 2/3 of the angular momentum from the inner accretion disc. Due to the weak efficiency of angular momentum transport by shear instability and meridional circulation during the accretion phase, the internal rotation profiles of accreting stars reflect essentially the angular momentum accretion history. As a consequence, careful choice of the angular momentum accretion history allows circumvention of any limitation in mass and velocity, and production of stars of any mass and velocity compatible with structure equations.

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.

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

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

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

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

Context. Fast rotating red giants in the upper part of the red giant branch have surface velocities that cannot be explained by single star evolution. Aims. We check whether tides between a star and a planet followed by planet engulfment can indeed accelerate the surface rotation of red giants for a sufficient long time in order to produce these fast rotating red giants. Methods. Using rotating stellar models, accounting for the redistribution of the angular momentum inside the star by different transport mechanisms, for the exchanges of angular momentum between the planet orbit and the star before the engulfment and for the deposition of angular momentum inside the star at the engulfment, we study how the surface rotation velocity at the stellar surface evolves. Results. We show that the surface velocities reached at the end of the orbital decay due to tidal forces and planet engulfment can be similar to values observed for fast rotating red giants. This surface velocity then decreases when the star evolves along the red giant branch but at a sufficiently slow pace for allowing stars to be detected with such a high velocity. More quantitatively, star-planet interaction can produce a rapid acceleration of the surface of the star, above values equal to 8 km s$^{-1}$, for periods lasting up to more than 30% the red giant branch phase. The changes of the surface carbon isotopic ratios produced by the dilution of the planetary material into the convective envelope is quite modest. Much more important might be the increase of the lithium abundance due to this effect. However lithium may be affected by many different, still uncertain, processes. Conclusions. The acceleration of the stellar surface to rotation velocities above limits that depend on the surface gravity does appear at the moment as the clearest signature of a star-planet interaction.[abridged]

Context. Globular clusters host multiple populations of long-lived low-mass stars whose origin remains an open question. Several scenarios have been proposed to explain the associated photometric and spectroscopic peculiarities. They differ, for instance, in the maximum helium enrichment they predict for stars of the second population, which these stars can inherit at birth as the result of the internal pollution of the cluster by different types of stars of the first population. Aims. We present the distribution of helium-rich stars in present-day globular clusters as it is expected in the original framework of the fast-rotating massive stars scenario (FRMS) as first-population polluters. We focus on NGC 6752. (to be continued)

The CEMP-no stars are long-lived small mass stars presenting a very low iron content and overabundances of carbon with no sign or only very weak signs for the presence of s- or r-elements. Although the origin of that abundance pattern is still a matter of debate, it was very likely inherited from a previous massive star, that we shall call here the source star. We rely on a recent classification of CEMP-no stars arguing that some of them are made of a material processed by hydrogen burning that was enriched in products of helium burning during the nuclear life of the source star. We examine the possibility of forming CEMP-no stars with this material. We study the nucleosynthesis of the CNO cycle and the Ne-Na Mg-Al chains in a hydrogen burning single zone while injecting the helium burning products $^{12}$C, $^{16}$O, $^{22}$Ne and $^{26}$Mg. We investigate the impact of changing the density and temperature, as well as the injection rate. The nuclear reaction rates involving the creation and destruction of $^{27}$Al are also examined. $^{14}$N, $^{23}$Na, $^{24}$Mg and $^{27}$Al are formed when injecting $^{12}$C, $^{16}$O, $^{22}$Ne and $^{26}$Mg in the hydrogen burning zone. The $^{12}$C/$^{13}$C ratio is constant under various conditions in the hydrogen burning zone. The predicted [Al/Fe] ratio varies up to $\sim 2$ dex depending on the prescription used for the reaction rates involving $^{27}$Al. The experiments we carried out support the view that some CEMP-no stars are made of a material processed by hydrogen burning, coming from a massive star experiencing mild-to-strong rotational mixing. During its burning, this material was likely enriched in helium burning products. No material coming from the carbon-oxygen rich core of the source star should be added to form the daughter star, otherwise the $^{12}$C/$^{13}$C ratio would be largely above the observed range of values.

Context. Efforts to look for signatures of the first stars have concentrated on metal-poor halo objects. However, the low end of the bulge metallicity distribution has been shown to host some of the oldest objects in the Milky Way and hence this Galactic component potentially offers interesting targets to look at imprints of the first stellar generations. As a pilot project, we selected bulge field stars already identified in the ARGOS survey as having [Fe/H] ~ -1 and oversolar [alpha/Fe] ratios, and we used FLAMES-UVES to obtain detailed abundances of key elements that are believed to reveal imprints of the first stellar generations. Aims. The main purpose of this study is to analyse selected ARGOS stars using new high-resolution (R~45,000) and high-signal-to-noise (S/N >100) spectra. We aim to derive their stellar parameters and elemental ratios, in particular the abundances of C, N, the alpha-elements O, Mg, Si, Ca, and Ti, the odd-Z elements Na and Al, the neutron-capture s-process dominated elements Y, Zr, La, and Ba, and the r-element Eu. Methods. High-resolution spectra of five field giant stars were obtained at the 8m VLT UT2-Kueyen telescope with the UVES spectrograph in FLAMES-UVES configuration. Spectroscopic parameters were derived based on the excitation and ionization equilibrium of Fe I and Fe II. The abundance analysis was performed with a MARCS LTE spherical model atmosphere grid and the Turbospectrum spectrum synthesis code.

We have endeavoured to understand the formation and evolution of the black hole (BH) X-ray binary LMC X-3. We estimate the properties of the system at 4 evolutionary stages: 1) at the Zero Age Main Sequence (ZAMS), 2) just prior to the supernova (SN) explosion of the primary, 3) just after the SN, and 4) at the moment of RLO onset.We use a hybrid approach, combining detailed stellar structure and binary evolution calculations with approximate population synthesis models. This allows us to estimate potential natal kicks and the evolution of the BH spin. In the whole analysis we incorporate as model constraints the most up-to-date observational information, encompassing the binary's orbital properties, the companion star mass, effective temperature, surface gravity and radius, as well as the black hole's mass and spin. We find that LMC X-3 began as a ZAMS system with the mass of the primary star in the range $M_{\rm{1,ZAMS}}$ = 22-31 $\rm{M_{\odot}}$ and a secondary star of $M_{\rm{2,ZAMS}} = 5.0-8.3M_{\odot}$, in a wide ($P_{ZAMS} \gtrsim 2.000\, \rm days$) and eccentric ($e_{\rm{ZAMS}} \gtrsim 0.23$) orbit. Just prior to the SN, the primary has a mass of $M_{\mathrm{1,preSN}} = 11.1-18.0\,\rm M_{\odot}$, with the secondary star largely unaffected. The orbital period decreases to $0.6-1.7\, \rm days$, and is still eccentric $0 \leq e_{\rm{preSN}} \leq 0.44$. We find that a symmetric SN explosion with no or small natal kicks (a few tens of $\rm km\, s^{-1}$) imparted on the BH cannot be formally excluded, however, large natal kicks in excess of $\gtrsim 120 \,\rm km\, s^{-1}$ increase the estimated formation rate by an order of magnitude. Following the SN, the system has a BH $M_{\mathrm{BH,postSN}} = 6.4-8.2\,\rm M_{\odot}$ and is put into an eccentric orbit. At the RLO onset the orbit is circularised and it has an orbital period of $P_{\rm{RLO}} = 0.8-1.4\,\rm days$.

Context. As a star evolves, the planet orbits change with time due to tidal interactions, stellar mass losses, friction and gravitational drag forces, mass accretion and evaporation on/by the planet. Stellar rotation modifies the structure of the star and therefore the way these different processes occur. Changes of the orbits, at their turn, have an impact on the rotation of the star. Aims. Models accounting in a consistent way for these interactions between the orbital evolution of the planet and the evolution of the rotation of the star are still missing. The present work is a first attempt to fill this gap. Methods. We compute the evolution of stellar models including a comprehensive treatment of rotational effects together with the evolution of planetary orbits, so that the exchanges of angular momentum between the star and the planetary orbit are treated in a self-consistent way. The evolution of the rotation of the star accounts for the angular momentum exchange with the planet and also follows the effects of the internal transport of angular momentum and chemicals. Results. We show that rotating stellar models without tidal interactions can well reproduce the surface rotations of the bulk of the red giants. However, models without any interactions cannot account for fast rotating red giants in the upper part of the red giant branch, where, such models, whatever the initial rotation considered on the ZAMS, always predict very low velocities. For those stars some interaction with a companion is highly probable and the present rotating stellar models with planets confirm that tidal interaction can reproduce their high surface velocities. We show also that the minimum distance between the planet and the star on the ZAMS that will allow the planet to avoid engulfment and survive is decreased around faster rotating stars. [abridged]

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.

Two main hypotheses for the origin of Galactic cosmic rays are the "supernova" and "superbubble" origin hypotheses. We analyse the evidence for the superbubble hypothesis provided by the measurements of the relativive abundances of isotopes of cobalt and nickel in the cosmic ray flux. We compare the measured upper limit on the abundance of 59Ni in the cosmic ray flux with the 59Ni abundance predictions of the up-to-date stellar evolution models. Non-detection of 59Ni in the cosmic ray flux has previously been attributed to a large time delay of the order of 1e5 yr between the moment of supernova explosion and the onset of particle acceleration process. This large time delay was considered as an argument in favour of the "superbubble" scenario. We show that the recent calculation of the 59Ni yield of massive stars, which takes into account the initial mass range up to 120 solar masses and includes stellar rotation, results in prediction of low 59Ni abundance relative to its decay product 59Co. The predicted abundance is consistent with the upper bound on 59Ni abundance in the cosmic ray flux for the supernova parameters assumed. This result removes the necessity of decay of 59Ni in the time interval between the supernova explosion and the onset of acceleration process and restores the consistency of measurements of 59Ni / 59Co abundances with the "supernova" hypothesis of the CR origin.

The nucleosynthesis in the first massive stars may be constrained by observing the surface composition of long-lived very iron-poor stars born around 10 billion years ago from material enriched by their ejecta. Many interesting clues on physical processes having occurred in the first stars can be obtained just based on nuclear aspects. Two facts are particularly clear, 1) in these first massive stars, mixing must have occurred between the H-burning and the He-burning zone during their nuclear lifetimes; 2) only the outer layers of these massive stars have enriched the material from which the very iron-poor stars, observed today in the halo, have formed. These two basic requirements can be obtained by rotating stellar models at very low metallicity. In the present paper, we discuss the arguments supporting this view and illustrates the sensitivity of the results concerning the [Mg/Al] ratio on the rate of the reaction Na(23)(p,\gamma)Mg(24).

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.

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.

Massive stars likely form by accretion and the evolutionary track of an accreting forming star corresponds to what is called the birthline in the HR diagram. The shape of this birthline is quite sensitive to the evolution of the entropy in the accreting star. We first study the reasons why some birthlines published in past years present different behaviours for a given accretion rate. We then revisit the question of the accretion rate, which allows us to understand the distribution of the observed pre-main-sequence (pre-MS) stars in the Hertzsprung-Russell (HR) diagram. Finally, we identify the conditions needed to obtain a large inflation of the star along its pre-MS evolution that may push the birthline towards the Hayashi line in the upper part of the HR diagram. We present new pre-MS models including accretion at various rates and for different initial structures of the accreting core. From the observed upper envelope of pre-MS stars in the HR diagram, we deduce the accretion law that best matches the accretion history of most of the intermediate-mass stars. In the case of cold disc accretion, the existence of a significant swelling during the accretion phase, which leads to radii $\gtrsim100\,R_\odot$ and brings the star back to the red part of the HR diagram, depends sensitively on the initial conditions. For an accretion rate of $10^{-3}\,M_\odot\rm\,yr^{-1}$, only models starting from a core with a significant radiative region evolve back to the red part of the HR diagram. We also obtain that, in order to reproduce the observed upper envelope of pre-MS stars in the HR diagram with an accretion law deduced from the observed mass outflows in ultra-compact HII regions, the mass effectively accreted onto the star with respect to the total in falling matter decreases when the mass of the star increases.

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.

The material used to form the CEMP-no stars presents signatures of material processed by the CNO cycle and by He-burning from a previous stellar generation called the source stars. In order to reproduce the relative abundance ratios like for instance C/N or $^{12}$C/$^{13}$C, some mixing between the two burning regions must have occured in the source stars and only the outer layers of the stars, with modest amount coming from the CO core, must have been expelled either through stellar winds or at the time of the (faint) supernova event. With new models at low metallicity including rotational mixing, we shall discuss how the variety of abundances observed for CEMP-no stars can be reproduced.

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

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 propose and apply a new classification for the CEMP-no stars, which are "carbon-enhanced metal-poor" stars with no overabundance of s-elements and with [Fe/H] generally inferior or equal to -2.5. This classification is based on the changes in abundances for the elements and isotopes involved in the CNO, Ne-Na, and Mg-Al nuclear cycles. These abundances change very much owing to successive back and forth mixing motions between the He- and H-burning regions in massive stars (the "source stars" responsible for the chemical enrichment of the CEMP-no stars). The wide variety of the ratios [C/Fe], 12C/13C, [N/Fe], [O/Fe], [Na/Fe], [Mg/Fe], [Al/Fe], [Sr/Fe], and [Ba/Fe], which are the main characteristics making the CEMP-no and low s stars so peculiar, is described well in terms of the proposed nucleosynthetic classification. We note that the [(C+N+O)/Fe] ratios significantly increase for lower values of [Fe/H]. The classification of CEMP-no stars and the behavior of [(C+N+O)/Fe] support the presence, in the first stellar generations of the Galaxy, of fast-rotating massive stars experiencing strong mixing and mass loss (spinstars). This result has an impact on the early chemical and spectral evolution of the Galaxy.

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

(Abridged) The CEMP--no stars are "carbon-enhanced-metal-poor" stars that in principle show no evidence of s-- and r--elements from neutron captures. We try to understand the origin and nucleosynthetic site of their peculiar CNO, Ne--Na, and Mg--Al abundances. We compare the observed abundances to the nucleosynthetic predictions of AGB models and of models of rotating massive stars with internal mixing and mass loss. We also analyze the different behaviors of $\alpha$-- and CNO--elements, as well the abundances of elements involved in the Ne--Na and Mg--Al cycles. We show that CEMP-no stars exhibit products of He--burning that have gone through partial mixing and processing by the CNO cycle, producing low $^{12}$C/$^{13}$C and a broad variety of [C/N] and [O/N] ratios. From a $^{12}$C/$^{13}$C vs. [C/N] diagram, we conclude that neither the yields of AGB stars (in binaries or not) nor the yields of classic supernovae can fully account for the observed CNO abundances in CEMP-no stars. Better agreement is obtained once the chemical contribution by stellar winds of fast-rotating massive stars is taken into account, where partial mixing takes place, leading to various amounts of CNO being ejected. CEMP--no stars present a wide variety in the [C/Fe], [N/Fe], [O/Fe], [Na/Fe], [Mg/Fe], [Al/Fe], and [Sr/Fe] ratios. We show that back-and-forth, partial mixing between the He-- and H--regions may account for this variety. Some s--elements, mainly of the first peak, may even be produced by these processes in a small fraction of the CEMP--no stars. We propose a classification scheme for the CEMP--no and low--s stars, based on the changes in composition produced by these successive back-and-forth mixing motions.

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.

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.

Globular clusters are among the oldest structures in the Universe and they host today low-mass stars and no gas. However, there has been a time when they formed as gaseous objects hosting a large number of short-lived, massive stars. Many details on this early epoch have been depicted recently through unprecedented dissection of low-mass globular cluster stars via spectroscopy and photometry. In particular, multiple populations have been identified, which bear the nucleosynthetic fingerprints of the massive hot stars long disappeared. Here we discuss how massive star archeology can been done through the lens of these multiple populations.

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.

We examine some debated points in current discussions about rotating stars: the shape, the gravity darkening, the critical velocities, the mass loss rates, the hydrodynamical instabilities, the internal mixing and N--enrichments. The study of rotational mixing requires high quality data and careful analysis. From recent studies where such conditions are fulfilled, rotational mixing is well confirmed. Magnetic coupling with stellar winds may produce an apparent contradiction, i.e. stars with a low rotation and a high N--enrichment. We point out that it rather confirms the large role of shears in differentially rotating stars for the transport processes. New models of interacting binaries also show how shears and mixing may be enhanced in close binaries which are either spun up or down by tidal interactions.

Magnetic configurations, stable on the long term, appear to exist in various evolutionary phases, from Main-Sequence stars to white dwarfs and neutron stars. The large scale ordered nature of these fields, often approximately dipolar, and their scaling according to the flux conservation scenario favor the model of a fossil field (Duez et al. 2010). We make some first estimates of the magnetic coupling between the stellar cores and the outer layers in red giants and supergiants. Analytical expressions of the truncation radius of the field coupling are established for a convective envelope and for a rotating radiative zone with horizontal turbulence. The timescales of the internal exchanges of angular momentum are considered. Numerical estimates are made on the basis of recent model grids. The direct magnetic coupling of the core to the extended convective envelope of red giants and supergiants appears unlikely. However, we find that the intermediate radiative zone is fully coupled to the core during the He-burning and later phases. This coupling is able to produce a strong spin down of the core of red giants and supergiants, also leading to relatively slowly rotating stellar remnants, like white dwarfs and pulsars. Some angular momentum is also transferred to the outer convective envelope of red giants and supergiants during the He-burning phase and later.

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

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

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

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

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.

(Abridged) Globular clusters host multiple stellar populations showing different sodium enrichments. These various populations can be observed along the main sequence, red giant and horizontal branch phases. Recently it was shown, however, that at least in the globular cluster NGC 6752, no sodium-rich stars are observed along the early asymptotic giant branch, posing an apparent problem for stellar evolution. We present an explanation for this lack of sodium-rich stars in this region of the colour-magnitude diagram.

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]

<Context> Ideally, one would like to know which type of core-collapse SNe is produced by different progenitors and the channels of stellar evolution leading to these progenitors. These links have to be very well known to use the observed frequency of different types of SN events for probing the star formation rate and massive star evolution in different types of galaxies. <Aims> We investigate the link between LBV as SN progenitors and the appearance of episodic radio light curve modulations of the SN event. <Methods> We use the 20Msun and 25Msun models with rotation at solar metallicity, part of an extended grid of stellar models computed by the Geneva team. At their pre-SN stage, these two models have recently been shown to have spectra similar to those of LBV stars and possibly explode as Type IIb SNe. Based on the wind properties before the explosion, we derive the density structure of their circumstellar medium. This structure is used as input for computing the SN radio light curve. <Results> We find that the 20Msun model shows radio light curves with episodic luminosity modulations, similar to those observed in some Type IIb SNe. This occurs because the evolution of the 20Msun model terminates in a region of the HR diagram where radiative stellar winds present strong density variations, caused by the bistability limit. The 25Msun model, ending its evolution in a zone of the HR diagram where no change of the mass-loss rates is expected, presents no such modulations in its radio SN light curve. <Conclusions> Our results reinforce the link between SN progenitors and LBV stars. We also confirm the existence of a physical mechanism for a single star to have episodic radio light curve modulations. In the case of the 25Msun progenitors, we do not obtain modulations in the radio light curve, but our models may miss some outbursting behavior in the late stages of massive stars.

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

The way angular momentum is built up in stars during their formation process may have an impact on their further evolution. In the frame of the cold disc accretion scenario, we study for the first time how angular momentum builds up inside the star during its formation and what are the consequences for its evolution on the main sequence (MS). Computation begins from a hydrostatic core on the Hayashi line of 0.7 Msol at solar metallicity (Z=0.014) rotating as a solid body. Accretion rates depending on the luminosity of the accreting object are considered varying between 1.5e-5 and 1.7e-3 Msol/yr. The accreted matter is assumed to have an angular velocity equal to that of the outer layer of the accreting star. Models are computed for a mass-range on the zero-age main sequence (ZAMS) between 2 and 22 Msol. We study how the internal and surface velocities vary as a function of time during the accretion phase and the evolution towards the ZAMS. Stellar models, whose evolution has been followed along the pre-MS phase, are found to exhibit a shallow gradient of angular velocity on the ZAMS. Interestingly, for masses on the ZAMS larger than 8 Msol, there exists a maximum surface velocity that can be reached through the present scenario of formation. Typically, for 14 Msol models, only stars with surface velocities on the ZAMS lower than about 45% of the critical velocity can be formed. To reach higher velocities would require to start from cores rotating above the critical limit. We find that this upper velocity limit is smaller for higher masses. In contrast, below 8 Msol, there is no restriction and the whole domain of velocities, up to the critical one, can be reached.

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.

Very old halo stars were previously found to show at least four different abundance 'anomalies', which models of fast rotating massive stars (spinstars) can successfully account for: rise of N/O and C/O, low 12C/13C and a primary-like evolution of Be and B. Here we show the impact of these same stars in the enrichment of Sr and Ba in the early Universe. We study if the s-process production of fast rotating massive stars can offer an explanation for the observed spread in [Sr/Ba] ratio in halo stars with metallicity [Fe/H]< -2.5. By means of a chemical inhomogeneous model we compute the enrichment of Sr and Ba by massive stars in the Galactic halo. Our model takes into account, for the first time, the contribution of spinstars. Our model (combining an r-process contribution with a s-process from fast rotating massive stars) is able to reproduce for the first time the observed scatter in the [Sr/Ba] ratio at [Fe/H]< -2.5. Toward higher metallicities, the stochasticity of the star formation fades away due to the increasing number of exploding and enriching stars, and as a consequence the predicted scatter decreases. Our scenario is again based on the existence of spinstars in the early Universe. Very old halo stars were previously found to show at least four other abundance 'anomalies', which rotating models of massive stars can successfully account for. Our results provide a 5th independent signature of the existence of fast rotating massive stars: an early enrichment of the Universe in s-process elements.

Aims. We examine the interactions of various instabilities in rotating stars, which usually are considered as independent. Methods. An analytical study of the problem is performed, account is given to radiative losses, mu-gradients and horizontal turbulence. Results. The diffusion coefficient for an ensemble of instabilities is not given by the sum of the specific coefficients for each instability, but by the solution of a general equation. We find that thermohaline mixing is possible in low-mass red giants only if the horizontal turbulence is very weak. In rotating stars the Rayleigh-Taylor and the shear instabilities need simultaneous treating. This has for consequence that rotation laws of the form 1/r^(alpha) are predicted to be unstable for alpha > 1.6568, while the usual Rayleigh criterion predicts instability only for alpha > 2. Also, the shear instabilities are somehow reduced in Main Sequence stars by the effect of the Rayleigh-Taylor criterion. Various instability criteria should be expressed differently in rotating stars than in simplified geometries.

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

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

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

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

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.

The self-enrichment scenario for globular clusters (GC) requires large amounts of residual gas after the initial formation of the first stellar generation. Recently, Krause et al. (2012) found that supernovae may not be able to expel that gas, as required to explain their present day gas-free state, and suggested that a sudden accretion on to the dark remnants, at a stage when type II supernovae have ceased, may plausibly lead to fast gas expulsion. Here, we explore the consequences of these results for the self-enrichment scenario via fast rotating massive stars (FRMS). We analyse the interaction of FRMS with the intra-cluster medium (ICM), in particular where, when and how the second generation of stars may form. From the results, we develop a timeline of the first approximately 40 Myr of GC evolution. The results of Paper I imply three phases during which the ICM is in a fundamentally different state, namely the wind bubble phase (lasting 3.5 to 8.8 Myr), the supernova phase (lasting 26.2 to 31.5 Myr), and the dark remnant accretion phase (lasting 0.1 to 4 Myr): (i) Quickly after the first generation massive stars have formed, stellar wind bubbles compress the ICM into thin filaments. No stars may form in the normal way during this phase, due to the high Lyman-Werner flux density. If the first generation massive stars have however equatorial ejections, as we proposed in the FRMS scenario, accretion may resume in the shadow of the equatorial ejecta. The second generation stars may then form due to gravitational instability in these discs that are fed by both the FRMS ejecta and pristine gas. [...]

A major paradigm shift has recently revolutionized our picture of globular clusters (GC) that were long thought to be simple systems of coeval stars born out of homogeneous material. Indeed, detailed abundance studies of GC long-lived low-mass stars performed with 8-10m class telescopes, together with high-precision photometry of Galactic GCs obtained with HST,have brought compelling clues on the presence of multiple stellar populations in individual GCs. These stellar subgroups can be recognized thanks to their different chemical properties (more precisely by abundance differences in light elements from carbon to aluminium; see Bragaglia, this volume) and by the appearance of multimodal sequences in the colour-magnitude diagrams (see Piotto, this volume). This has a severe impact on our understanding of the early evolution of GCs, and in particular of the possible role that massive stars played in shaping the intra-cluster medium (ICM) and in inducing secondary star formation. Here we summarize the detailed timeline we have recently proposed for the first 40 Myrs in the lifetime of a typical GC following the general ideas of our so-called "Fast Rotating Massive stars scenario" (FRMS, Decressin et al. 2007b) and taking into account the dynamics of interstellar bubbles produced by stellar winds and supernovae. More details can be found in Krause et al. (2012, 2013).

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

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

We summarize the results presented in Krause et al. (2012) on the impact of supernova-driven shells and dark-remnant accretion on gas expulsion in globular cluster infancy.

The picture of the early evolution of globular clusters has been significantly revised in recent years. Current scenarios require at least two generations of stars of which the first generation (1G), and therefore also the protocluster cloud, has been much more massive than the currently predominating second generation (2G). Fast gas expulsion is thought to unbind the majority of the 1G stars. Gas expulsion is also mandatory to remove metal-enriched supernova ejecta, which are not found in the 2G stars. It has long been thought that the supernovae themselves are the agent of the gas expulsion, based on crude energetics arguments. Here, we assume that gas expulsion happens via the formation of a superbubble, and describe the kinematics by a thin-shell model. We find that supernova- driven shells are destroyed by the Rayleigh-Taylor instability before they reach escape speed for all but perhaps the least massive and most extended clusters. More power is required to expel the gas, which might plausibly be provided by a coherent onset of accretion onto the stellar remnants. The resulting kpc-sized bubbles might be observable in Faraday rotation maps with the planned Square Kilometre Array radio telescope against polarised background radio lobes if a globular cluster would happen to form in front of such a radio lobe.

Little is known about the stellar environment and the genealogy of our solar system. Short-lived radionuclides (SLRs, mean lifetime shorter than 100 Myr) that were present in the solar protoplanetary disk 4.56 Gyr ago could potentially provide insight into that key aspect of our history, were their origin understood. Previous models failed to provide a reasonable explanation of the abundance of two key SLRs, 26Al (mean lifetime = 1.1 Myr) and 60Fe (mean lifetime = 3.7 Myr), at the birth of the solar system by requiring unlikely astrophysical conditions. Our aim is to propose a coherent and generic solution based on the most recent understanding of star-forming mechanisms. Iron-60 in the nascent solar system is shown to have been produced by a diversity of supernovae belonging to a first generation of stars in a giant molecular cloud. Aluminum-26 is delivered into a dense collected shell by a single massive star wind belonging to a second star generation. The Sun formed in the collected shell as part of a third stellar generation. Aluminum-26 yields used in our calculation are based on new rotating stellar models in which 26Al is present in stellar winds during the star main sequence rather than during the Wolf-Rayet phase alone. Our scenario eventually constrains the time sequence of the formation of the two stellar generations that just preceded the solar system formation, along with the number of stars born in these two generations. (Continued).

Aims: We study the influence of rotation and disc lifetime on lithium depletion of pre-main sequence (PMS) solar-type stars. Methods: The impact of rotational mixing and of the hydrostatic effects of rotation on lithium abundances are investigated by computing non-rotating and rotating PMS models that include a comprehensive treatment of shellular rotation. The influence of the disc lifetime is then studied by comparing the lithium content of PMS rotating models experiencing different durations of the disc-locking phase between 3 and 9 Myr. Results: The surface lithium abundance at the end of the PMS is decreased when rotational effects are included. During the beginning of the lithium depletion phase, only hydrostatic effects of rotation are at work. This results in a decrease in the lithium depletion rate for rotating models compared to non-rotating ones. When the convective envelope recedes from the stellar centre, rotational mixing begins to play an important role due to differential rotation near the bottom of the convective envelope. This mixing results in a decrease in the surface lithium abundance with a limited contribution from hydrostatic effects of rotation, which favours lithium depletion during the second part of the PMS evolution. The impact of rotation on PMS lithium depletion is also found to be sensitive to the duration of the disc-locking phase. When the disc lifetime increases, the PMS lithium abundance of a solar-type star decreases owing to the higher efficiency of rotational mixing in the radiative zone. A relationship between the surface rotation and lithium abundance at the end of the PMS is then obtained: slow rotators on the zero-age main sequence are predicted to be more lithium-depleted than fast rotators due to the increase in the disc lifetime.

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.

Supernovae Ib/c are likely to be associated to long GRBs, therefore it is important to compare the SN rate in galaxies with the GRB rate. To do that we computed Type Ib/c SN rates in galaxies of different morphological type by assuming different histories of star formation and different supernova Ib/c progenitors. We included some recent suggestions about the dependence of the minimum mass of single Wolf-Rayet (WR) stars upon the stellar metallicity and therefore upon galactic chemical evolution. We adopted several cosmic star formation rates as functions of cosmic time, either observationally or theoretically derived, including the one computed with our galaxy models. Then we computed the cosmic Type Ib/c SN rates. We derived the following conclusions: i) the ratio cosmic GRB - Type Ib/c rate varies in the range 10^{-2}-10^{-4} in the whole redshift range, thus suggesting that only a small fraction of all the Type Ib/c SNe gives rise to GRBs. ii) The metallicity dependence of Type Ib/c SN progenitors produces lower cosmic SN Ib/c rates at early times, for any chosen cosmic star formation rate. iii) Different theoretical cosmic star formation rates, computed under different scenarios of galaxy formation, produce SN Ib/c cosmic rates which differ mainly at very high redshift. However, it is difficult to draw firm conclusions on the high redshift trend because of the large uncertainties in the data. iv) GRBs can be important tracers of star formation at high redshift if their luminosity function does not vary with redshift and they can help in discriminating among different galaxy formation models.

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

Yellow and red supergiants are evolved massive stars whose numbers and locations on the HR diagram can provide a stringent test for models of massive star evolution. Previous studies have found large discrepancies between the relative number of yellow supergiants observed as a function of mass and those predicted by evolutionary models, while a disagreement between the predicted and observed locations of red supergiants on the HR diagram was only recently resolved. Here we extend these studies by examining the yellow and red supergiant populations of M33. Unfortunately, identifying these stars is difficult as this portion of the color-magnitude diagram is heavily contaminated by foreground dwarfs. We identify the red supergiants through a combination of radial velocities and a two-color surface gravity discriminant and, after re-characterizing the rotation curve of M33 with our newly selected red supergiants, we identify the yellow supergiants through a combination of radial velocities and the strength of the OI $\lambda$7774 triplet. We examine ~1300 spectra in total and identify 121 yellow supergiants (a sample which is unbiased in luminosity above log(L/L\odot) ~ 4.8) and 189 red supergiants. After placing these objects on the HR diagram, we find that the latest generation of Geneva evolutionary tracks show excellent agreement with the observed locations of our red and yellow supergiants, the observed relative number of yellow supergiants with mass and the observed red supergiant upper mass limit. These models therefore represent a drastic improvement over previous generations.

Due to their transitionary nature, yellow supergiants provide a critical challenge for evolutionary modeling. Previous studies within M31 and the SMC show that the Geneva evolutionary models do a poor job at predicting the lifetimes of these short-lived stars. Here we extend this study to the LMC while also investigating the galaxy's red supergiant content. This task is complicated by contamination by Galactic foreground stars that color and magnitude criteria alone cannot weed out. Therefore, we use proper motions and the LMC's large systemic radial velocity (\sim278 km/s) to separate out these foreground dwarfs. After observing nearly 2,000 stars, we identified 317 probable yellow supergiants, 6 possible yellow supergiants and 505 probable red supergiants. Foreground contamination of our yellow supergiant sample was \sim80%, while that of the the red supergiant sample was only 3%. By placing the yellow supergiants on the H-R diagram and comparing them against the evolutionary tracks, we find that new Geneva evolutionary models do an exemplary job at predicting both the locations and the lifetimes of these transitory objects.

We study the effects of feebly or non-annihilating weakly interacting Dark Matter (DM) particles on stars that live in DM environments denser than that of our Sun. We find that the energy transport mechanism induced by DM particles can produce unusual conditions in the core of Main Sequence stars, with effects which can potentially be used to probe DM properties. We find that solar mass stars placed in DM densities of rhochi>= e2 GeV/cm3 are sensitive to Spin-Dependent scattering cross-section sigmsd >= e-37 cm2 and a DM particle mass as low as mchi=5 GeV, accessing a parameter range weakly constrained by current direct detection experiments.

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.

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

This article first reviews the basic physics of rotating stars and their evolution. We examine in particular the changes of the mechanical and thermal equilibrium of rotating stars. An important (predicted and observed) effect is that rotating stars are hotter at the poles and cooler at the equator. We briefly discuss the mass loss by stellar winds, which are influenced by the anisotropic temperature distribution. These anisotropies in the interior are also driving circulation currents, which transports the chemical elements and the angular momentum in stars. Internal differential rotation, if present, creates instabilities and mixing, in particular the shear mixing, the horizontal turbulence and their interactions. A major check of the model predictions concerns the changes of the surface abundances, which are modified by mass loss in the very massive stars and by rotational mixing in O- and B-type stars. We show that the observations confirm the existence of rotational mixing, with much larger effects at lower metallicities. We discuss the predictions of stellar models concerning the evolution of the surface velocities, the evolutionary tracks in the HR diagram and lifetimes, the populations of blue, red supergiants and Wolf-Rayet stars, and the progenitors of type Ibc supernovae. We show, that in many aspects, rotating models provide a much better fit than non-rotating ones. Using the same physical ingredients as those which fit the best the observations of stars at near solar metallicities, we explore the consequences of rotating models for the status of Be stars, the progenitors of Gamma Ray Bursts, the evolution of Pop III stars and of very metal poor stars, the early chemical evolution of galaxies, the origin of the C-enhanced Metal Poor stars (CEMP) and of the chemical anomalies in globular clusters.

Massive stars are crucial building blocks of galaxies and the universe, as production sites of heavy elements and as stirring agents and energy providers through stellar winds and supernovae. The field of magnetic massive stars has seen tremendous progress in recent years. Different perspectives -- ranging from direct field measurements over dynamo theory and stellar evolution to colliding winds and the stellar environment -- fruitfully combine into a most interesting and still evolving overall picture, which we attempt to review here. Zeeman signatures leave no doubt that at least some O- and early B-type stars have a surface magnetic field. Indirect evidence, especially non-thermal radio emission from colliding winds, suggests many more. The emerging picture for massive stars shows similarities with results from intermediate mass stars, for which much more data are available. Observations are often compatible with a dipole or low order multi-pole field of about 1 kG (O-stars) or 300 G to 30 kG (Ap / Bp stars). Weak and unordered fields have been detected in the O-star zeta-Ori A and in Vega, the first normal A-type star with a magnetic field. Theory offers essentially two explanations for the origin of the observed surface fields: fossil fields, particularly for strong and ordered fields, or different dynamo mechanisms, preferentially for less ordered fields. Numerical simulations yield the first concrete stable (fossil) field configuration, but give contradictory results as to whether dynamo action in the radiative envelope of massive main sequence stars is possible. Internal magnetic fields, which may not even show up at the stellar surface, affect stellar evolution as they lead to a more uniform rotation, with more slowly rotating cores and faster surface rotation. Surface metallicities may become enhanced, thus affecting the mass-loss rates.

(abridged) A detailed study of the blue supergiant UIT005 (B2-2.5Ia+) in M33 is presented. The results of our quantitative spectral analysis indicate that the star is a very luminous, log(L/Lsun)~5.9 dex, and massive, M~50 Msun, object, showing a very high nitrogen-to-oxygen ratio in its surface (N/O~8, by mass). Based on the derived Mg and Si abundances, we argue that this high N/O ratio cannot be the result of an initial low O content due to its location on the disk of M33, known to present a steep metallicity gradient. In combination with the He abundance, the most plausible interpretation is that UIT005 is in an advanced stage of evolution, showing in its surface N enrichment and O depletion resulting from mixing with CNO processed material from the stellar interior. A comparison with the predictions of current stellar evolutionary models indicates that there are significant discrepancies, in particular with regard to the degree of chemical processing, with the models predicting a much lower degree of O depletion than observed. At the same time, the mass-loss rate derived in our analysis is an order of magnitude lower than the values considered in the evolutionary calculations. Based on a study of the surrounding stellar population and the nearby cluster NGC588, using WFPC2 photometry, we suggest that UIT005 could be in fact a runaway star from this cluster.

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)

We examine the nature and role of mass loss via an equatorial decretion disk in massive stars with near-critical rotation induced by evolution of the stellar interior. In contrast to the usual stellar wind mass loss set by exterior driving from the stellar luminosity, such decretion-disk mass loss stems from the angular momentum loss needed to keep the star near and below critical rotation, given the interior evolution and decline in the star's moment of inertia. Because the specific angular momentum in a Keplerian disk increases with the square root of the radius, the decretion mass loss associated with a required level of angular momentum loss depends crucially on the outer radius for viscous coupling of the disk, and can be significantly less than the spherical mass loss the spherical, wind-like mass loss commonly assumed in evolutionary calculations. We discuss the physical processes that affect the outer disk radius, including thermal disk outflow, and ablation of the disk material via a line-driven wind induced by the star's radiation. We present parameterized scaling laws for taking account of decretion-disk mass loss in stellar evolution codes, including how these are affected by metallicity, or by presence within a close binary and/or a dense cluster. Effects similar to those discussed here should also be present in accretion disks during star formation, and may play an important role in shaping the distribution of rotation speeds on the ZAMS.

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.

Magnetic fields at the surface of a few early-type stars have been directly detected. These fields have magnitudes between a few hundred G up to a few kG. In one case, evidence of magnetic braking has been found. We investigate the effects of magnetic braking on the evolution of rotating ($\upsilon_{\rm ini}$=200 km s$^{-1}$) 10 M$_\odot$ stellar models at solar metallicity during the main-sequence (MS) phase. The magnetic braking process is included in our stellar models according to the formalism deduced from 2D MHD simulations of magnetic wind confinement by ud-Doula and co-workers. Various assumptions are made regarding both the magnitude of the magnetic field and of the efficiency of the angular momentum transport mechanisms in the stellar interior. When magnetic braking occurs in models with differential rotation, a strong and rapid mixing is obtained at the surface accompanied by a rapid decrease in the surface velocity. Such a process might account for some MS stars showing strong mixing and low surface velocities. When solid-body rotation is imposed in the interior, the star is slowed down so rapidly that surface enrichments are smaller than in similar models with no magnetic braking. In both kinds of models (differentially or uniformly rotating), magnetic braking due to a field of a few 100 G significantly reduces the angular momentum of the core during the MS phase. This reduction is much greater in solid-body rotating models.

Aims: The effects of rotation and magnetic fields on the surface abundances of solar-type stars are studied in order to investigate whether the reported difference in lithium content of exoplanet-host stars can be related to their rotational history. Moreover, the asteroseismic properties predicted for stars with and without exoplanets are compared to determine how such a scenario, which relates the lithium abundances and the rotational history of the star, can be further challenged by observations of solar-like oscillations. Methods: Based on observations of rotational periods of solar-type stars, slow rotators on the zero age main sequence (ZAMS) are modelled with a comprehensive treatment of only the shellular rotation, while fast rotators are modelled including both shellular rotation and magnetic fields. Assuming a possible link between low rotation rates on the ZAMS and the presence of planets as a result of a longer disc-locking phase during the pre-main sequence (PMS), we compare the surface abundances and asteroseismic properties of slow and fast rotating models, which correspond to exoplanet-host stars and stars without detected planets, respectively. Results: We confirm previous suggestions that the difference in the lithium content of stars with and without detected planets can be related to their different rotational history. The larger efficiency of rotational mixing predicted in exoplanet-host stars explains their lithium depletion and also leads to changes in the structure and chemical composition of the central stellar layers. Asteroseismic observations can reveal these changes and can help us distinguish between different possible explanations for the lower lithium content of exoplanet-host 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.

The influence of rotational mixing on the evolution and asteroseismic properties of solar-type stars is studied. Rotational mixing changes the global properties of a solar-type star with a significant increase of the effective temperature resulting in a shift of the evolutionary track to the blue part of the HR diagram. These differences are related to changes of the chemical composition, because rotational mixing counteracts the effects of atomic diffusion leading to larger helium surface abundances for rotating models than for non-rotating ones. Higher values of the large frequency separation are then found for rotating models than for non-rotating ones at the same evolutionary stage, because the increase of the effective temperature leads to a smaller radius and hence to an increase of the stellar mean density. Rotational mixing also has a considerable impact on the structure and chemical composition of the central stellar layers by bringing fresh hydrogen fuel to the core, thereby enhancing the main-sequence lifetime. The increase of the central hydrogen abundance together with the change of the chemical profiles in the central layers result in a significant increase of the values of the small frequency separations and of the ratio of the small to large separations for models including shellular rotation. This increase is clearly seen for models with the same age sharing the same initial parameters except for the inclusion of rotation as well as for models with the same global stellar parameters and in particular the same location in the HR diagram. By computing rotating models of solar-type stars including the effects of a dynamo that possibly occurs in the radiative zone, we find that the efficiency of rotational mixing is strongly reduced when the effects of magnetic fields are taken into account, in contrast to what happens in massive stars.

We study the effect of dark matter (DM) particles in the Sun, focusing in particular on the possible reduction of the solar neutrinos flux due to the energy carried away by DM particles from the innermost regions of the Sun, and to the consequent reduction of the temperature of the solar core. We find that in the very low-mass range between 4 and 10 GeV, recently advocated to explain the findings of the DAMA and CoGent experiments, the effects on neutrino fluxes are detectable only for DM models with very small, or vanishing, self-annihilation cross section, such as the so-called asymmetric DM models, and we study the combination of DM masses and Spin Dependent cross sections which can be excluded with current solar neutrino data. Finally, we revisit the recent claim that DM models with large self-interacting cross sections can lead to a modification of the position of the convective zone, alleviating or solving the solar composition problem. We show that when the `geometric' upper limit on the capture rate is correctly taken into account, the effects of DM are reduced by orders of magnitude, and the position of the convective zone remains unchanged.

The yellow supergiant content of nearby galaxies provides a critical test of massive star evolutionary theory. While these stars are the brightest in a galaxy, they are difficult to identify because a large number of foreground Milky Way stars have similar colors and magnitudes. We previously conducted a census of yellow supergiants within M31 and found that the evolutionary tracks predict a yellow supergiant duration an order of magnitude longer than we observed. Here we turn our attention to the SMC, where the metallicity is 10x lower than that of M31, which is important as metallicity strongly affects massive star evolution. The SMC's large radial velocity (~160 km/s) allows us to separate members from foreground stars. Observations of ~500 candidates yielded 176 near-certain SMC supergiants, 16 possible SMC supergiants, along with 306 foreground stars and provide good relative numbers of yellow supergiants down to 12Mo. Of the 176 near-certain SMC supergiants, the kinematics predicted by the Besancon model of the Milky Way suggest a foreground contamination of >4%. After placing the SMC supergiants on the H-R diagram and comparing our results to the Geneva evolutionary tracks, we find results similar to those of the M31 study: while the locations of the stars on the H-R diagram match the locations of evolutionary tracks well, the models over-predict the yellow supergiant lifetime by a factor of ten. Uncertainties about the mass-loss rates on the main-sequence thus cannot be the primary problem with the models.

Aims: We test predictions of evolution models on mixing of CNO-cycled products in massive stars from a fundamental perspective. Relative changes within the theoretical C:N:O abundance ratios and the buildup of helium are compared with observational results. Methods: A sample of well-studied Galactic massive stars is presented. High-quality optical spectra are carefully analysed using improved NLTE line-formation and comprehensive analysis strategies. The results are put in the context of the existing literature data. Results: A tight trend in the observed N/C vs. N/O ratios and the buildup of helium is found from the self-consistent analysis of main-sequence to supergiant stars for the first time. The catalytic nature of the CNO-cycles is confirmed quantitatively, though further investigations are required to derive a fully consistent picture. Our observational results support the case of strong mixing, as predicted e.g. by evolution models that consider magnetic fields or by models that have gone through the first dredge-up in the case of many supergiants.

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

The Cygnus region harbours a huge complex of massive stars at a distance of 1.0-2.0kpc from us. About 170 O stars are distributed over several OB associations, among which the Cyg OB2 cluster is by far the most important with about 100-120 O stars. These massive stars inject large quantities of radioactive nuclei into the interstellar medium, such as 26Al and 60Fe, and their gamma-ray line decay signals can provide insight into the physics of massive stars and core-collapse supernovae. Past studies of the nucleosynthesis activity of Cygnus have concluded that the level of 26Al decay emission as deduced from CGRO/COMPTEL observations was a factor 2-3 above the predictions based on the theoretical yields available at that time and on the observed stellar content of the Cygnus region. We reevaluate the situation from new measurements of the gamma-ray decay fluxes with INTEGRAL/SPI and new predictions based on recently improved stellar models including some of the effects of stellar rotation for the higher mass stars and a coherent estimate of the contribution from SNIb/c. We developed a population synthesis code to predict the nucleosynthesis activity and corresponding decay fluxes of a given stellar population of massive stars. The observed decay fluxes from the Cygnus complex are found to be consistent with the values predicted by population synthesis at solar metallicity. The observed extent of the 1809keV emission from Cygnus is found to be consistent with the result of a numerical simulation of the diffusion of 26Al inside the superbubble blown by Cyg OB2. Our work indicates that the past dilemma regarding the gamma-ray line emission from Cygnus resulted from an overestimate of the 1809keV flux of the Cygnus complex, combined with an underestimate of the nucleosynthesis yields.

The Cygnus region harbours a huge complex of massive stars at a distance of 1.0-2.0kpc from us. About 170 O stars are distributed over several OB associations, among which the Cyg OB2 cluster is by far the largest with about 100-120 O stars. As a consequence of their successive nuclear-burning episodes, these massive stars inject large quantities of radioactive nuclei into the interstellar medium such as 26Al and 60Fe. The gamma-ray line signal from the latter is a solid tracer of ongoing nucleosynthesis. We want to compare the decay emission from the Cygnus region with predictions of recently improved stellar models. As a first step, we establish observational constraints upon the gamma-ray line emission from 26Al and 60Fe, with particular emphasis placed on separating the emission due to the Cygnus complex from the foreground and background mean Galactic contributions. We used the high-resolution gamma-ray spectrometer INTEGRAL/SPI to analyse the 26Al and 60Fe decay signal from the Cygnus region. The weak gamma-ray line emissions at 1809 and 1173/1332keV have been characterised in terms of photometry, spectrometry, and source morphology. The 1809keV emission from Cygnus is centred on the position of the Cyg OB2 cluster and has a typical size of 9 degrees or more. The total 1809keV flux from the Cygnus region is (6.0 +/- 1.0) x10e-5 ph/cm2/s, but the flux really attributable to the Cygnus complex reduces to (3.9 +/- 1.1) x10e-5 ph/cm2/s. The 1809keV line centroid is agrees with expectations considering the direction and distance of the Cygnus complex, and the line width is consistent with typical motions of the interstellar medium. No decay emission from 60Fe has been observed from the Cygnus region and an upper limit of 1.6 x10e-5 ph/cm2/s was derived.

The influence of rotation on the properties of red giants is studied in the context of the asteroseismic modelling of these stars. While red giants exhibit low surface rotational velocities, we find that the rotational history of the star has a large impact on its properties during the red giant phase. In particular, for stars massive enough to ignite He burning in non-degenerate conditions, rotational mixing induces a significant increase of the stellar luminosity and shifts the location of the core helium burning phase to a higher luminosity in the HR diagram. This of course results in a change of the seismic properties of red giants at the same evolutionary state. As a consequence the inclusion of rotation significantly changes the fundamental parameters of a red giant star as determined by performing an asteroseismic calibration. In particular rotation decreases the derived stellar mass and increases the age. Depending on the rotation law assumed in the convective envelope and on the initial velocity of the star, non-negligible values of rotational splitting can be reached, which may complicate the observation and identification of non-radial oscillation modes for red giants exhibiting moderate surface rotational velocities. By comparing the effects of rotation and overshooting, we find that the main-sequence widening and the increase of the H-burning lifetime induced by rotation (Vini=150 km/s) are well reproduced by non-rotating models with an overshooting parameter of 0.1, while the increase of luminosity during the post-main sequence evolution is better reproduced by non-rotating models with overshooting parameters twice as large. This is due to the fact that rotation not only increases the size of the convective core but also changes the chemical composition of the radiative zone.

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

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

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

The yellow supergiant content of nearby galaxies can provide a critical test of stellar evolution theory, bridging the gap between the hot, massive stars and the cool red supergiants. But, this region of the color-magnitude diagram is dominated by foreground contamination, requiring membership to somehow be determined. Fortunately, the large negative systemic velocity of M31, coupled to its high rotation rate, provides the means for separating the contaminating foreground dwarfs from the bona fide yellow supergiants within M31. Using the MMT, we obtained spectra of about 2900 stars, selected using the color and magnitude range to be yellow supergiants. Comparing the velocities to that of M31's rotation curve, we identified 54 certain, and 66 probable yellow supergiants from among the sea of foreground dwarfs. We find excellent agreement between the location of yellow supergiants in the H-R diagram and that predicted by the latest Geneva evolutionary tracks which include rotation. However, the relative number of yellow supergiants seen as a function of mass varies from that predicted by the models by a factor of more than 10, in the sense that more high mass yellow supergiants are predicted than are actually observed. Comparing the total number of yellow supergiants with masses above 20Mo with the estimated number of unevolved O stars indicates that the duration of the yellow supergiant phase is about 3000 years. This is consistent with what the 12Mo and 15Mo evolutionary tracks predict, but disagrees with the 20,000-80,000 year time scales predicted by the models for higher masses.

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 developed a new population synthesis code for groups of massive stars, where we model the emission of different forms of energy and matter from the stars of the association. In particular, the ejection of the two radioactive isotopes 26Al and 60Fe is followed, as well as the emission of hydrogen ionizing photons, and the kinetic energy of the stellar winds and supernova explosions. We investigate various alternative astrophysical inputs and the resulting output sensitivities, especially effects due to the inclusion of rotation in stellar models. As the aim of the code is the application to relatively small populations of massive stars, special care is taken to address their statistical properties. Our code incorporates both analytical statistical methods applicable to small populations, as well as extensive Monte Carlo simulations. We find that the inclusion of rotation in the stellar models has a large impact on the interactions between OB associations and their surrounding interstellar medium. The emission of 26Al in the stellar winds is strongly enhanced, compared to non-rotating models with the same mass-loss prescription. This compensates the recent reductions in the estimates of mass-loss rates of massive stars due to the effects of clumping. Despite the lower mass-loss rates, the power of the winds is actually enhanced for rotating stellar models. The supernova power (kinetic energy of their ejecta) is decreased due to longer lifetimes of rotating stars, and therefore the wind power dominates over supernova power for the first 6 Myr after a burst of star-formation. For populations typical of nearby star-forming regions, the statistical uncertainties are large and clearly non-Gaussian.

Red supergiants are a short-lived stage in the evolution of moderately massive stars (10-25Mo), and as such their location in the H-R diagram provides an exacting test of stellar evolutionary models. Since massive star evolution is strongly affected by the amount of mass-loss a star suffers, and since the mass-loss rates depend upon metallicity, it is highly desirable to study the physical properties of these stars in galaxies of various metallicities. Here we identify a sample of red supergiants in M31 (the most metal-rich of the Local Group galaxies) and derive their physical properties by fitting MARCS atmosphere models to moderate resolution optical spectroscopy, and from V-K photometry.

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.

We propose to advance investigations of electromagnetic radiation originating in atomic nuclei beyond its current infancy to a true astronomy. This nuclear emission is independent from conditions of gas, thus complements more traditional stronomical methods used to probe the nearby universe. Radioactive gamma-rays arise from isotopes which are made in specific locations inside massive stars, their decay in interstellar space traces an otherwise not directly observable hot and tenuous phase of the ISM, which is crucial for feedback from massive stars. Its intrinsic clocks can measure characteristic times of processes within the ISM. Frontier questions that can be addressed with studies in this field are the complex interiors of massive stars and supernovae which are key agents in galactic dynamics and chemical evolution, the history of star-forming and supernova activity affecting our solar-system environment, and explorations of occulted and inaccessible regions of young stellar nurseries in our Galaxy.

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.

A significant fraction of stars in globular clusters (about 70%-85%) exhibit peculiar chemical patterns with strong abundance variations in light elements along with constant abundances in heavy elements. These abundance anomalies can be created in the H-burning core of a first generation of fast rotating massive stars and the corresponding elements are convoyed to the stellar surface thanks to rotational induced mixing. If the rotation of the stars is fast enough this matter is ejected at low velocity through a mechanical wind at the equator. It then pollutes the ISM from which a second generation of chemically anomalous stars can be formed. The proportion of anomalous to normal star observed today depends on at least two quantities : (1) the number of polluter stars; (2) the dynamical history of the cluster which may lose during its lifetime first and second generation stars in different proportions. Here we estimate these proportions based on dynamical models for globular clusters. When internal dynamical evolution and dissolution due to tidal forces are accounted for, starting from an initial fraction of anomalous stars of 10% produces a present day fraction of about 25%, still too small with respect to the observed 70-85%. In case gas expulsion by supernovae is accounted for, much higher fraction is expected to be produced. In this paper we also address the question of the evolution of the second generation stars that are He-rich, and deduce consequences for the age determination of globular clusters.

Starting from a few topical astrophysical questions which require the knowledge of the age of Pop I stars, we discuss the needed precision on the age in order to make progresses in these areas of research. Then we review the effects of various inputs of the stellar models on the age determination and try to identify those affecting the most the lifetimes of stars.

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.

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.

The goal of this paper is to analyze the impact of a primary neutron source on the s-process nucleosynthesis in massive stars at halo metallicity. Recent stellar models including rotation at very low metallicity predict a strong production of primary N14. Part of the nitrogen produced in the H-burning shell diffuses by rotational mixing into the He core where it is converted to Ne22 providing additional neutrons for the s process. We present nucleosynthesis calculations for a 25 Msun star at [Fe/H] = -3, -4, where in the convective core He-burning about 0.8 % in mass is made of primary Ne22. The usual weak s-process shape is changed by the additional neutron source with a peak between Sr and Ba, where the s-process yields increase by orders of magnitude with respect to the yields obtained without rotation. Iron seeds are fully consumed and the maximum production of Sr, Y and Zr is reached. On the other hand, the s-process efficiency beyond Sr and the ratio Sr/Ba are strongly affected by the amount of Ne22 and by nuclear uncertainties, first of all by the Ne22(alpha,n)Mg25 reaction. Finally, assuming that Ne22 is primary in the considered metallicity range, the s-process efficiency decreases with metallicity due to the effect of the major neutron poisons Mg25 and Ne22. This work represents a first step towards the study of primary neutron source effect in fast rotating massive stars, and its implications are discussed in the light of spectroscopic observations of heavy elements at halo metallicity.

The Gaia satellite will be launched at the end of 2011. It will observe at least 1 billion stars, and among them several million emission line stars and hot stars. Gaia will provide parallaxes for each star and spectra for stars till V magnitude equal to 17. After a general description of Gaia, we present the codes and methods, which are currently developed by our team. They will provide automatically the astrophysical parameters and spectral classification for the hot and emission line stars in the Milky Way and other close Local Group galaxies such as the Magellanic Clouds.

Massive stars played a key role in the early evolution of the Universe. They formed with the first halos and started the re-ionisation. It is therefore very important to understand their evolution. In this paper, we describe the strong impact of rotation induced mixing and mass loss at very low $Z$. The strong mixing leads to a significant production of primary nitrogen 14, carbon 13 and neon 22. Mass loss during the red supergiant stage allows the production of Wolf-Rayet stars, type Ib,c supernovae and possibly gamma-ray bursts (GRBs) down to almost Z=0 for stars more massive than 60 solar masses. Galactic chemical evolution models calculated with models of rotating stars better reproduce the early evolution of N/O, C/O and C12/C13. We calculated the weak s-process production induced by the primary neon 22 and obtain overproduction factors (relative to the initial composition, Z=1.e-6) between 100-1000 in the mass range 60-90.

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.

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

(Abridged) Rotation has been shown to play a determinant role at very low metallicity, bringing heavy mass loss where almost none was expected. Is this still true when the metallicity strictly equals zero? The aim of our study is to get an answer to this question, and to determine how rotation changes the evolution and the chemical signature of the primordial stars. We have calculated 14 differentially-rotating and non-rotating stellar models at zero metallicity, with masses between 9 and 200 Msol. The evolution has been followed up to the pre-supernova stage. We find that Z=0 models rotate with an internal profile Omega(r) close to local angular momentum conservation, because of a very weak core-envelope coupling. Rotational mixing drives a H-shell boost due to a sudden onset of CNO cycle in the shell. This boost leads to a high 14N production. Generally, the rotating models produce much more metals than their non-rotating counterparts. The mass loss is very low, even for the models that reach the critical velocity during the main sequence. Due to the low mass loss and the weak coupling, the core retains a high angular momentum at the end of the evolution. The high rotation rate at death probably leads to a much stronger explosion than previously expected, changing the fate of the models. The inclusion of our yields in a chemical evolution model of the Galactic halo predicts log values of N/O, C/O and 12C/13C ratios of -2.2, -0.95 and 50 respectively at log O/H +12 = 4.2.

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

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

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

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

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

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

In these lecture notes, we present the equations presently used in stellar interior models in order to compute the effects of axial rotation. We discuss the hypotheses made. We suggest that the effects of rotation might play a key role at low metallicity.

Globular clusters as $\omega$ Cen and NGC 2808 appear to have a population of very He-rich stars. From a theoretical point of view, one expects the presence of He-rich stars in all globular clusters showing an oxygen-sodium (O-Na) anticorrelation. In this paper, we briefly recall how fast rotating massive stars could be the main source of the material from which He-rich stars have formed. We speculate that the UV-upturn phenomenon observed in all elliptical galaxies might be due to He-rich stars. If this hypothesis is correct then fast rotating massive stars may have played in the early evolutionary phases of these systems a similar role as the one they played in the nascent globular clusters.

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.

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

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

Aims. We analyse the effects of a first generation of fast rotating massive stars on the dynamical and chemical properties of globular clusters. Methods. We use stellar models of fast rotating massive stars, losing mass through a slow mechanical equatorial winds to produce material rich in H-burning products. We propose that stars with high Na and low O abundances (hereafter anomalous stars) are formed from matter made of slow winds of individual massive stars and of interstellar matter. The proportion of slow wind and of interstellar material is fixed in order to reproduce the observed Li-Na anticorrelation in NGC 6752. Results. In the case that globular clusters, during their lifetime, did not lose any stars, we found that to reproduce the observed ratio of normal to anomalous stars, a flat initial mass function (IMF) is needed, with typically a slope x=0.55 (a Salpeter's IMF has x=1.35). In the case that globular clusters suffer from an evaporation of normal stars, the IMF slope can be steeper: to have x=1.35, about 96% of the normal stars would be lost. We make predictions for the distribution of stars as a function of their [O/Na] and obtain quite reasonable agreement with that one observed for NGC 6752. Predictions for the number fraction of stars with different values of helium, of the 12C/13C and 16O/17O ratios are discussed, as well as the expected relations between values of [O/Na] and those of helium, of [C/N], of 12C/13C and of 16O/17O. Future observations might test these predictions. We also provide predictions for the present day mass of the clusters expressed in units of mass of the gas used to form stars, and for the way the present day mass is distributed between the first and second generation of stars and the stellar remnants.

We recall the basic physical principles governing the evolution of stars with some emphasis on the role played by the nuclear reactions. We argue that in general it is not possible from observations of stars to deduce constraints on the nuclear reaction rates. This is the reason why precise measurements of nuclear reaction rates are a necessity in order to make progresses in stellar physics, nucleosynthesis and chemical evolution of galaxies. There are however some stars which provides useful constraint on nuclear processes. The Wolf-Rayet stars of the WN type present at their surface CNO equilibrium patterns. There is also the particular case of the abundance of $^{22}$Ne at the surface of WC stars. The abundance of this element is a measure of the initial CNO content. Very interestingly, recent determinations of its abundance at the surface of WC stars tend to confirm that massive stars in the solar neighborhood have initial metallicities in agreement with the Asplund et al. (2005) solar abundances.

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

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

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

We review some important observed properties of massive stars. Then we discuss how mass loss and rotation affect their evolution and help in giving better fits to observational constraints. Consequences for nucleosynthesis at different metallicities are discussed. Mass loss appear to be the key feature at high metallicity, while rotation is likely dominant at low and very low metallicities. We discuss various indications supporting the view that very metal poor stars had their evolution strongly affected by rotational mixing. Many features, like the origin of primary nitrogen at low metallicity, that of the C-rich extremely metal poor halo stars, of He-rich stars in massive globular clusters, of the O-Na anticorrelation in globular clusters may be related to the existence of a population of very fast rotating metal poor stars that we tentatively call the {\it spinstars}. A fraction of these {\it spinstars} may also be the progenitors of GRB in very metal poor regions. They may avoid pair instability explosion due to the heavy mass loss undergone during their early evolutionary phases and be, dependent on their frequency, important sources of ionising photons in the early Universe.

We have measured the isotopic abundances of neon and a number of other species in the galactic cosmic rays (GCRs) using the Cosmic Ray Isotope Spectrometer (CRIS) aboard the ACE spacecraft. Our data are compared to recent results from two-component Wolf-Rayet (WR) models. The three largest deviations of galactic cosmic ray isotope ratios from solar-system ratios predicted by these models, 12C/16O, 22Ne/20Ne, and 58Fe/56Fe, are very close to those observed. All of the isotopic ratios that we have measured are consistent with a GCR source consisting of ~20% of WR material mixed with ~80% material with solar-system composition. Since WR stars are evolutionary products of OB stars, and most OB stars exist in OB associations that form superbubbles, the good agreement of our data with WR models suggests that OB associations within superbubbles are the likely source of at least a substantial fraction of GCRs. In previous work it has been shown that the primary 59Ni (which decays only by electron-capture) in GCRs has decayed, indicating a time interval between nucleosynthesis and acceleration of >10^5 yr.In this paper we suggest a scenario that should allow much of the 59Ni to decay in the OB association environment and conclude that the hypothesis of the OB association origin of cosmic rays appears to be viable.

We investigate the effect of new stellar models, which take rotation into account, computed for very low metallicities on the chemical evolution of the earliest phases of the Milky Way. We check the impact of these new stellar yields on a model for the halo of the Milky Way that can reproduce the observed halo metallicity distribution. In this way we try to better constrain the ISM enrichment timescale, which was not done in our previous work. The stellar models adopted in this work were computed under the assumption that the ratio of the initial rotation velocity to the critical velocity of stars is roughly constant with metallicity. This naturally leads to faster rotation at lower metallicity, as metal poor stars are more compact than metal rich ones. We find that the new Z = 10-8 stellar yields computed for large rotational velocities have a tremendous impact on the interstellar medium nitrogen enrichment for log(O/H)+12 < 7 (or [Fe/H]< -3). We show that upon the inclusion of the new stellar calculations in a chemical evolution model for the galactic halo with infall and outflow, both high N/O and C/O ratios are obtained in the very-metal poor metallicity range in agreement with observations. Our results give further support to the idea that stars at very low metallicities could have initial rotational velocities of the order of 600-800kms-1. An important contribution to N from AGB stars is still needed in order to explain the observations at intermediate metallicities. One possibility is that AGB stars at very low metallicities also rotate fast. This could be tested in the future, once stellar evolution models for fast rotating AGB stars will be available.

We present and discuss evolutionary synthesis models for massive stellar populations generated with the Starburst99 code in combination with a new set of stellar evolution models accounting for rotation. The new stellar evolution models were compiled from several data releases of the Geneva group and cover heavy-element abundances ranging from twice solar to one fifth solar. The evolution models were computed for rotation velocities on the zero-age main-sequence of 0 and 300 km/s and with the latest revision of stellar mass-loss rates. Since the mass coverage is incomplete, in particular at non-solar chemical composition, our parameter study is still preliminary and must be viewed as exploratory. Stellar population properties computed with Starburst99 and the new evolution models show some marked differences in comparison with models obtained using earlier tracks. Since individual stars now tend to be more luminous and bluer when on the blue side of the Hertzsprung-Russell diagram, the populations mirror this trend. For instance, increases by factors of two or more are found for the light-to-mass ratios at ultraviolet to near-infrared wavelengths, as well as for the output of hydrogen ionizing photons. If these results are confirmed once the evolution models have matured, recalibrations of certain star-formation and initial mass function indicators will be required.

After a brief review of the observational evidences indicating how the populations of Be stars, red/blue supergiants, Wolf-Rayet stars vary as a function of metallicity, we discuss the implications of these observed trend for our understanding of the massive star evolution. We show how the inclusion of the effects of rotation in stellar models improves significantly the correspondence between theory and observation.

We study the effect of wind anisotropies on the stellar evolution leading to collapsars. Rotating models of a 60 M$_\odot$ star with $\Omega/\Omega_{\rm crit}=0.75$ on the ZAMS, accounting for shellular rotation and a magnetic field, with and without wind anisotropies, are computed at $Z$=0.002 until the end of the core He-burning phase. Only the models accounting for the effects of the wind anisotropies retain enough angular momentum in their core to produce a Gamma Ray Burst (GRB). The chemical composition is such that a type Ic supernova event occurs. Wind anisotropies appear to be a key physical ingredient in the scenario leading to long GRBs.

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

After a review of the many effects of metallicity on the evolution of rotating and non-rotating stars, we discuss the consequences of a high metallicity on massive star populations and on stellar nucleosynthesis. The most striking effect of high metallicity is to enhance the amount of mass lost by stellar winds. Typically at a metallicity $Z=0.001$ only 9% of the total mass returned by non-rotating massive stars is ejected by winds (91% by supernovae explosion), while at solar metallicity this fraction may amount to more than 40%. High metallicity favors the formation of Wolf-Rayet stars and of type Ib supernovae. It however disfavors the occurrence of type Ic supernovae. We estimate empirical yields of carbon based on the observed population of WC stars in the solar neighborhood, and obtain that WC stars eject between 0.2 and 0.4% of the mass initially locked into stars under the form of new synthesized carbon. Models give values well in agreement with these empirical yields. Chemical evolution models indicate that such carbon yields may have important impacts on the abundance of carbon at high metallicity.

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

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

We calculated pre-supernova evolution models of single rotating massive stars. These models reproduce observations during the early stages of the evolution very well, in particular Wolf--Rayet (WR) populations and ratio between type II and type Ib,c supernovae at different metallicities (Z). Using these models we found the following results concerning long and soft gamma--ray burst (GRB) progenitors: - GRBs coming from WO--type (SNIc) WR stars are only produced at low Z (LMC or lower). - The upper metallicity limit for GRBs is reduced to Z ~ 0.004 (SMC) when the effects of magnetic fields are included. - GRBs are predicted from the second (and probably the first) stellar generation onwards.

Some indirect observations, as the high fraction of Be stars at low metallicity, or the necessity for massive stars to be important sources of primary nitrogen, seem to indicate that very metal poor stars were fast rotators. As a consequence of this fast rotation, these stars, contrarily to current wisdom, might lose large amounts of mass during their lifetime. In this paper, we review various mechanisms triggered by rotation which may induce strong mass loss at very low metallicity. The most efficient process comes from surface enrichments in CNO elements which then drive mass loss by stellar winds. Due to this process, a fast rotating 60 M$_\odot$ with metallicities in the range of $Z=10^{-8}$ and $10^{-5}$, can lose between 30 and 55% of its initial mass. This rotationally wind ejected material participates to the chemical evolution of the interstellar medium, enriching it exclusively in H- and He-burning products. In particular, metal poor fast rotating stars may play a key role for explaining the origin of the peculiar abundance pattern observed at the surface of the extremely metal-poor C-rich stars, for explaining the chemical inhomogeneities observed in globular clusters, and the presence of stars in $\omega$ Cen with a very high helium content .

We use the MARCS stellar atmosphere to derive the physical properties of 36 red supergiants (RSGs) in the LMC, and 39 RSGs in the SMC using moderate-resolution optical spectrophotometry (4000-9000A) and broad-band colors (V-R, V-K). The results from the dereddened V-R colors are in good agreement with those derived from the spectrophotometry, but the dereddened V-K colors give temperatures that are 3-4% warmer for the SMC data, with the LMC and Milky Way showing a smaller but similar effect. We conclude that this discrepancy is due to the limitations of 1D models. Our newly derived effective temperatures and bolometric luminosities bring the Magellanic Cloud RSGs into good agreement with stellar evolutionary models that include the effects of rotation. A typical M2~I in the SMC is about 150 K cooler than its Galactic counterpart; one in the LMC is about 50 K cooler. This is in the sense expected due to the lower chemical abundances in the SMC and LMC, although it is not sufficient to explain the shift in average RSG spectral type seen between the SMC, LMC, and Milky Way. Instead, that is due primarily to the change in Hayashi limit with metallicity, as first proposed by Elias et al. (1985). Finally, our study confirms that many RSGs in the Magellanic Clouds are significantly more reddened than OB stars, consistent with our recent findings for Galactic stars that circumstellar dust may contribute several magnitudes of extra visual extinction.

It has been speculated that WR winds may have contaminated the forming solar system, in particular with short-lived radionuclides (half-lives in the approximate 10^5 - 10^8 y range) that are responsible for a class of isotopic anomalies found in some meteoritic materials. We revisit the capability of the WR winds to eject these radionuclides using new models of single non-exploding WR stars with metallicity Z = 0.02. The earlier predictions for non-rotating WR stars are updated, and models for rotating such stars are used for the first time in this context. We find that (1) rotation has no significant influence on the short-lived radionuclide production by neutron capture during the core He-burning phase, and (2) 26Al, 36Cl, 41Ca, and 107Pd can be wind-ejected by a variety of WR stars at relative levels that are compatible with the meteoritic analyses for a period of free decay of around 10^5 y between production and incorporation into the forming solar system solid bodies. We confirm the previously published conclusions that the winds of WR stars have a radionuclide composition that can meet the necessary condition for them to be a possible contaminating agent of the forming solar system. Still, it remains to be demonstrated from detailed models that this is a sufficient condition for these winds to have provided a level of pollution that is compatible with the observations.

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

The blue Main Sequence (bMS) of $\omega$ Cen implies a ratio of helium to metal enrichment $\Delta Y/\Delta Z \approx 70$, which is a major enigma. We show that rotating models of low metallicity stars, which account for the anomalous abundance ratios of extremely metal poor stars, are also useful for understanding the very high $\Delta Y/\Delta Z$ ratio in $\omega$ Cen. Models of massive stars with moderate initial rotation velocities produce stellar winds with large He-- and N--excesses, but without the large C-- (and O--) excesses made by very fast rotation, in agreement with the observed chemical abundance ratios in $\omega$ Cen. It is still uncertain whether the abundance peculiarities of $\omega$ Cen result from the fact that the high velocity contributions of supernovae escaped the globular cluster, usually considered as a tidally stripped core of a dwarf galaxy. Another possibility is a general dominance of wind ejecta at very low $Z$, due to the formation of black holes. Some abundance and isotopic ratios like $Mg/Al$, $Na/Mg$, $Ne/N$, $^{12}C/^{13}C$, $^{16}O/^{18}O$ and $^{17}O/^{18}O$ may allow us to further discriminate between these scenarios and between the AGB and massive star contributions.

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

We have re-evaluated empirical expressions for the abundance determination of N, O, Ne, S, Cl, Ar and Fe taking into account the latest atomic data and constructing an appropriate grid of photoionization models with state-of-the art model atmospheres. Using these expressions we have derived heavy element abundances in the $\sim$ 310 emission-line galaxies from the Data Release 3 of the Sloan Digital Sky Survey (SDSS)with an observed Hbeta flux F(Hbeta)> 1E-14 erg s^{-1} cm^{-2} and for which the [O III] 4363 emission line was detected at least at a 2sigma level, allowing abundance determination by direct methods. The oxygen abundance 12 + log O/H of the SDSS galaxies lies in the range from ~ 7.1 (Zsun/30) to 8.5 (0.7 Zsun). The SDSS sample is merged with a sample of 109 blue compact dwarf (BCD) galaxies with high quality spectra, which contains extremely low-metallicity objects. We use the merged sample to study the abundance patterns of low-metallicity emission-line galaxies. We find that extremely metal-poor galaxies (12 + log O/H < 7.6, i.e. Z < Zsun/12) are rare in the SDSS sample. The alpha element-to-oxygen abundance ratios do not show any significant trends with oxygen abundance, in agreement with previous studies, except for a slight increase of Ne/O with increasing metallicity, which we interpret as due to a moderate depletion of O onto grains in the most metal-rich galaxies. The Fe/O abundance ratio is smaller than the solar value, by up to 1 dex at the high metallicity end. (abridged)

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

Using single star models including the effects of shellular rotation with and without magnetic fields, we show that massive stars at solar metallicity with initial masses lower than about 20-25 M$_\odot$ and with an initial rotation above $\sim 350$ km s$^{-1}$ likely reach the critical velocity during their Main-Sequence phase. This results from the efficient outwards transport of angular momentum by the meridional circulation. This could be a scenario for explaining the Be stars. After the Main-Sequence phase, single star in this mass range can again reach the critical limit when they are on a blue loop after a red supergiant phase \citep{HL98}. This might be a scenario for the formation of B[e] stars, however as discussed by Langer & Heger (1998), this scenario would predict a short B[e] phase (only some 10$^4$ years) with correspondingly small amounts of mass lost.

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

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

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

We examine the generation of a magnetic field in a solar-like star and its effects on the internal distribution of the angular velocity. We suggest that the evolution of a rotating star with magnetic fields leads to an equilibrium value of the differential rotation. This equilibrium is determined by the magnetic coupling, which favours a constant rotation profile, and meridional circulation which tends to build differential rotation. The global equilibrium stage is close to solid body rotation between about 0.7 and 0.2 R_sun, in good agreement with helioseismic measurements.

The abundances of neon isotopes in the galactic cosmic rays (GCRs) are reported using data from the Cosmic Ray Isotope Spectrometer (CRIS) aboard the Advanced Composition Explorer (ACE). We compare our ACE-CRIS data for neon and refractory isotope ratios, and data from other experiments, with recent results from two-component Wolf-Rayet (WR) models. The three largest deviations of GCR isotope ratios from solar-system ratios predicted by these models are indeed present in the GCRs. Since WR stars are evolutionary products of OB stars, and most OB stars exist in OB associations that form superbubbles, the good agreement of these data with WR models suggests that superbubbles are the likely source of at least a substantial fraction of GCRs.

New models of rotating and non-rotating stars are computed for initial masses between 25 and 120 Msun and for metallicities Z = 0.004, 0.008, 0.020 and 0.040 with the aim of reexamining the wind contribution of Wolf-Rayet (WR) stars to the F19 enrichment of the interstellar medium. Models with an initial rotation velocity vini = 300 km/s are found to globally eject less F19 than the non-rotating models. We compare our new predictions with those of Meynet & Arnould (2000), and demonstrate that the F19 yields are very sensitive to the still uncertain F19(alpha,p)Ne22 rate and to the adopted mass loss rates. Using the recommended mass loss rate values that take into account the clumping of the WR wind and the NACRE reaction rates when available, we obtain WR F19 yields that are significantly lower than predicted by Meynet & Arnould (2000), and that would make WR stars non-important contributors to the galactic F19 budget. In view, however, of the large nuclear and mass loss rate uncertainties, we consider that the question of the WR contribution to the galactic F19 remains quite largely open.

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

We examine the effects of the magnetic field created by the Tayler--Spruit dynamo in differentially rotating stars. Magnetic fields of the order of a few $10^4$ G are present through most of the stellar envelope, with the exception of the outer layers. The diffusion coefficient for the transport of angular momentum is very large and it imposes nearly solid body rotation during the MS phase. In turn, solid body rotation drives meridional circulation currents which are much faster than usual and leads to much larger diffusion coefficients than the magnetic diffusivity for the chemical species. The consequence is that the interplay of the thermal and magnetic instabilities favours the chemical transport of elements, while there would be no transport in models with magnetic field only. We also discuss the effects on the stellar interior, lifetimes and HR diagram.

We use moderate-resolution optical spectrophotometry and the new MARCS stellar atmosphere models to determine the effective temperatures of 74 Galactic red supergiants. From these we find a new effective temperature scale that is significantly warmer than those in the literature. We show that this temperature scale, along with the newly derived bolometric corrections, gives much better agreement between our red supergiants and stellar evolutionary tracks. This agreement provides an independent verification of our new temperature scale. The combination of effective temperature and bolometric luminosities allows us to calculate stellar radii; the coolest and most luminous stars have radii of roughly 1500 solar radii (7 AU), in excellent accordance with the largest stellar radii predicted from current evolutionary theory. We find that similar results are obtained for the effective temperatures and bolometric luminosities using only the de-reddened V-K colors, providing a powerful demonstration of the self-consistency of the MARCS models.

We present a new set of stellar yields obtained from rotating stellar models at solar metallicity covering the massive star range (12-60 solar masses). The stellar models were calculated with the latest version of the Geneva stellar evolution code described in Hirschi et al (2004). Evolution and nucleosynthesis are in general followed up to silicon burning. The yields of our non-rotating models are consistent with other calculations and differences can be understood in the light of the treatment of convection and the rate used for C12(a,g)O16. This verifies the accuracy of our calculations and gives a safe basis for studying the effects of rotation on the yields. The contributions from stellar winds and supernova explosions to the stellar yields are presented separately. We then add the two contributions to compute the total stellar yields. Below about 30 solar masses, rotation increases the total metal yields, Z, and in particular the yields of carbon and oxygen by a factor of 1.5-2.5. As a rule of thumb, the yields of a rotating 20 solar masses star are similar to the yields of a non-rotating 30 solar masses star, at least for the light elements considered in this work. For very massive stars (around 60 solar masses), rotation increases the yield of helium but does not significantly affect the yields of heavy elements.

We present a new set of stellar yields obtained from rotating stellar models at solar metallicity covering the massive star range (9-120 solar masses). The stellar models were calculated with the latest version of the Geneva stellar evolution code described in Hirschi et al 2004a. The contributions from stellar winds and from supernova explosions to the stellar yields were calculated separately. The two contributions were then added to compute the total stellar yields (Hirschi et al 2004b). The effects of rotation on pre-supernova models are significant between 15 and 30 solar masses. Above 20 solar masses, rotation may change the radius or colour of the supernova progenitors (blue instead of red supergiant) and the supernova type (IIb or Ib instead of II). Rotation increases the helium and CO core sizes by a factor about 1.5. Thus, rotation increases the total metallic yields and in particular the yields for carbon and oxygen by a factor 1.5-2.5. For Wolf-Rayet stars (M>30 solar masses), rotation increases the yields of helium and other hydrogen burning products but does not significantly affect the yields of heavy elements.

We present new yields of 26Al from Wolf--Rayet stellar winds based on rotating stellar models which account well for numerous observed properties of massive stars. We study the impacts on the yields of a change of initial mass, metallicity and initial rotation velocity. We also consider the effects of a change of mass loss rates during the Wolf--Rayet phase. We show that for surface rotation velocities during the core H--burning phase matching to the observed ones, the quantity of 26Al ejected by a star of a given initial mass and metallicity is roughly doubled when the effects of rotation are taken into account. The metallicity dependence of the yield is, on the other hand, very similar to that obtained from non--rotating models. We estimate that at least about 20% to 50% (e.g. 0.6 --1.4 Msun) of the live 26Al detected in the Milky--Way originates from Wolf--Rayet stellar winds. We show the importance of a good knowledge of the present metallicity gradient and star formation rate in our galaxy for modeling both the variation of the 26Al surface density with the galactocentric distance and the global contribution of the Wolf--Rayet stellar winds to the present galactic mass of 26Al.

The main instabilities induced by rotation in stellar interiors are described. We derive from first principles the general equation describing the transport of the angular momentum. The case of the transport of the chemical species is also discussed. As long as the mass loss rates are not too important, meridional currents, by advecting angular momentum from the inner regions to the outer layers, accelerate the stellar surface during the Main Sequence phase. A 9 M$_\odot$ stellar model at solar metallicity with an equatorial velocity at the beginning of the core H--burning phase equal to 340 km s$^{-1}$ reaches the break--up limit during the MS phase. The model with an initial velocity of 290 km s$^{-1}$ approaches this limit without reaching it. The models with 290 km s$^{-1}$ and 340 km s$^{-1}$ predict enhancements of the N/C ratio at the end of the MS phase equal to 2.8 and 3.2 times the initial value respectively.

Stellar rotation produces an internal mixing of the elements due to shear instability and meridional circulation. This leads to observable $N/C$ enhancements in massive stars above about 7--9 $M_{\odot}$. Rotation also favours mass loss by stellar winds. Mass loss effects dominate for masses above 30 $M_{\odot}$, while mixing dominates below that limit. The effects of mixing are also much larger at lower metallicity $Z$, because the internal $\Omega$--gradients are steeper. This appears to be in agreement with observations in the SMC. At very low $Z$ and Z=0, mixing between the He--burning core and the H--burning shell leads to the production of primary N in intermediate mass stars. Such enrichments increase the metallicity of the rotating star, also massive Z=0 stars with moderate initial velocities currently reach break--up velocity during a fraction of the MS phase. Both effects favour mass loss in Z=0 stars, which have ejecta with abundance anomalies very similar to those of C--rich very metal poor stars.

We show that rotation strongly affects the nature of the supernova progenitor (blue/red supergiant or Wolf--Rayet star), and thus the supernova types. In particular our models well reproduce the variations of the number ratio SNIb/Ic to SNII with metallicity. Rotation also produces envelope enrichments of the N/C ratio, and increases the size of the CO cores. We show the evolution of the specific angular momentum up to the preSN stage and make comparison with neutron stars. We suggest that the rare WO stars, preferentially formed at low metallicity, are the progenitors of GRB.

We discuss the physics of the $\Omega \Gamma$-- Limit, i.e. when the star is unbound as a result of {\emph{both}} rotation and radiation pressure. We suggest that the $\Omega \Gamma$-- Limit is what makes the Humphreys--Davidson Limit. Stellar filiations are discussed, with an emphasis on the final stages, in particular on the final masses, their angular momentum and the chemical yields with account of rotation. A possible relation between WO stars and GRB is emphasized.

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

Grids of models of massive stars ($M \ge$ 20 $M_\odot$) with rotation are computed for metallicities $Z$ ranging from that of the Small Magellanic Cloud (SMC) to that of the Galactic Centre. The hydrostatic effects of rotation, the rotational mixing and the enhancements of the mass loss rates by rotation are included. The evolution of the surface rotational velocities of the most massive O--stars mainly depends on the mass loss rates and thus on the initial $Z$ value. The minimum initial mass for a star for entering the Wolf--Rayet (WR) phase is lowered by rotation. For all metallicities, rotating stars enter the WR phase at an earlier stage of evolution and the WR lifetimes are increased, mainly as a result of the increased duration of the eWNL phase. Models of WR stars predict in general rather low rotation velocities ($ < 50$ km s$^{-1}$) with a few possible exceptions, particularly at metallicities lower than solar where WR star models have in general faster rotation and more chance to reach the break--up limit.The properties of the WR populations as predicted by the rotating models are in general in much better agreement with the observations in nearby galaxies. The observed variation with metallicity of the fractions of type Ib/Ic supernovae with respect to type II supernovae as found by Prantzos & Boissier (\cite{Pr03}) is very well reproduced by the rotating models, while non--rotating models predict much too low ratios.

The 26Al radionuclide can be detected through its decay emission line at 1.809 MeV, as was first observed by Mahoney et al. (1982). Since then, COMPTEL on board of the CGRO satellite, performed a sky survey in this energy range, and provided maps of the 26Al distribution in the Galaxy. These results revealed that the main contributors to the synthesis of 26Al are most likely the massive stars, which contribute through their winds (Wolf-Rayet stars) and through their supernova explosion. Comparison between these observations (in particular observations in localized regions such as the Vela region and the Cygnus region) and the models available at that moment, showed however the need for improvements from both theoretical and observational points of view, in order to improve our understanding of the 26Al galactic distribution as well as that of its synthesis. With the launch of the INTEGRAL satellite in October 2002, the observational part will hopefully be improved, and the construction of better resolution maps at 1.809 MeV is one of the main aims of the mission. From a theoretical point of view, we need the most up-to-date predictions in order to be able to interpret the forthcoming data. In this paper, we address this latter part, and present new results for 26Al production by rotating Wolf-Rayet stars and their contribution to the total amount observed in the Galaxy.

We describe the latest developments of the Geneva stellar evolution code in order to model the pre-supernova evolution of rotating massive stars. Rotating and non-rotating stellar models at solar metallicity with masses equal to 12, 15, 20, 25, 40 and 60 solar masses were computed from the ZAMS until the end of the core silicon burning phase. We took into account meridional circulation, secular shear instabilities, horizontal turbulence and dynamical shear instabilities. Most of the differences between the pre-supernova structures obtained from rotating and non-rotating stellar models have their origin in the effects of rotation during the core hydrogen and helium burning phases. The effects of rotation on pre-supernova models are significant between 15 and 30 solar masses. Indeed, rotation increases the core sizes (and the yields) by a factor ~ 1.5. Above 20 solar masses, rotation may change the colour of the supernova progenitors (blue instead of red supergiant) and the supernova type (Ib instead of II). Rotation affects the lower mass limits for radiative core carbon burning, for iron core collapse and for black hole formation. For Wolf-Rayet stars (M > 30 solar masses), the pre-supernova structures are mostly affected by the intensities of the stellar winds and less by rotational mixing. Finally, the core of our rotating WR stars contain enough angular momentum to produce GRBs.

We further develop the Tayler--Spruit dynamo theory, based on the most efficient instability for generating magnetic fields in radiative layers of differentially rotating stars. We avoid the simplifying assumptions that either the $\mu$-- or the $T$--gradient dominates, but we treat the general case and we also account for the nonadiabatic effects, which favour the growth of the magnetic field. Stars with a magnetic field rotate almost as a solid body. Several of their properties (size of the core, MS lifetimes, tracks, abundances) are closer to those of models without rotation than with rotation only. In particular, the observed N/C or N/H excesses in OB stars are better explained by our previous models with rotation only than by the present models with magnetic fields that predict no nitrogen excesses. We show that there is a complex feedback loop between the magnetic instability and the thermal instability driving meridional circulation. This opens the possibility for further magnetic models, but at this stage we do not know the relative importance of the magnetic fields due to the Tayler instability in stellar interiors.

Using a high resolution spectrum of the secondary star in the black hole binary A0620-00, we have derived the stellar parameters and veiling caused by the accretion disk in a consistent way. We have used a chi^2 minimization procedure to explore a grid of 800 000 LTE synthetic spectra computed for a plausible range of both stellar and veiling parameters. Adopting the best model parameters found, we have determined atmospheric abundances of Fe, Ca, Ti, Ni and Al. The Fe abundance of the star is [Fe/H]=0.14 +- 0.20. Except for Ca, we found the other elements moderately over-abundant as compared with stars in the solar neighborhood of similar iron content. Taking into account the small orbital separation, the mass transfer rate and the mass of the convection zone of the secondary star, a comparison with element yields in supernova explosion models suggests a possible explosive event with a mass cut comparable to the current mass of the compact object. We have also analyzed the Li abundance, which is unusually high for a star of this spectral type and relatively low mass.

Detailed models of alpha Cen A and B based on new seismological data for alpha Cen B by Carrier & Bourban (2003) have been computed using the Geneva evolution code including atomic diffusion. Taking into account the numerous observational constraints now available for the alpha Cen system, we find a stellar model which is in good agreement with the astrometric, photometric, spectroscopic and asteroseismic data. The global parameters of the alpha Cen system are now firmly constrained to an age of t=6.52+-0.30 Gyr, an initial helium mass fraction Y_i=0.275+-0.010 and an initial metallicity (Z/X)_i=0.0434+-0.0020. Thanks to these numerous observational constraints, we confirm that the mixing-length parameter alpha of the B component is larger than the one of the A component, as already suggested by many authors (Noels et al. 1991, Fernandes & Neuforge 1995 and Guenther & Demarque 2000): alpha_B is about 8% larger than alpha_A (alpha_A=1.83+-0.10 and alpha_B=1.97+-0.10). Moreover, we show that asteroseismic measurements enable to determine the radii of both stars with a very high precision (errors smaller than 0.3%). The radii deduced from seismological data are compatible with the new interferometric results of Kervella et al. (2003) even if they are slightly larger than the interferometric radii (differences smaller than 1%).

We describe and discuss the properties of three numerical methods for solving the diffusion equation for the transport of the chemical species and of the angular momentum in stellar interiors. We study through numerical experiments both their accuracy and their ability to provide physical solutions. On the basis of new tests and analyses applied to the stellar astrophysical context, we show that the most robust method to follow the secular evolution is the implicit finite differences method. The importance of correctly estimating the diffusion coefficient between mesh points is emphasized and a procedure for estimating the average diffusion coefficient between a convective and a radiative zone is described.

We present new $^{26}$Al stellar yields from rotating Wolf--Rayet stellar models which, at solar metallicity, well reproduce the observed properties of the Wolf-Rayet populations. These new yields are enhanced with respect to non--rotating models, even with respect to non--rotating models computed with enhanced mass loss rates. We briefly discuss some implications of the use of these new yields for estimating the global contribution of Wolf-Rayet stars to the quantity of $^{26}$Al now present in the Milky Way.

The Geneva evolutionary code has been modified to study the advanced stages (Ne, O, Si burnings) of rotating massive stars. Here we present the results of four 20 solar mass stars at solar metallicity with initial rotational velocities of 0, 100, 200 and 300 km/s in order to show the crucial role of rotation in stellar evolution. As already known, rotation increases mass loss and core masses (Meynet and Maeder 2000). A fast rotating 20 solar mass star has the same central evolution as a non-rotating 26 solar mass star. Rotation also increases strongly net total metal yields. Furthermore, rotation changes the SN type so that more SNIb are predicted (see Meynet and Maeder 2003 and N. Prantzos and S. Boissier 2003). Finally, SN1987A-like supernovae progenitor colour can be explained in a single rotating star scenario.

We compare the current effects of rotation in stellar evolution to those of the magnetic field created by the Tayler instability. In stellar regions, where magnetic field can be generated by the dynamo due to differential rotation (Spruit 2002), we find that the growth rate of the magnetic instability is much faster than for the thermal instability. Thus, meridional circulation is negligible with respect to the magnetic fields, both for the transport of angular momentum and of chemical elements. Also, the horizontal coupling by the magnetic field, which reaches values of a few $10^5$ G, is much more important than the effects of the horizontal turbulence. The field, however, is not sufficient to distort the shape of the equipotentials. We impose the condition that the energy of the magnetic field created by the Tayler--Spruit dynamo cannot be larger than the energy excess present in the differential rotation. This leads to a criterion for the existence of the magnetic field in stellar interiors. Numerical tests are made in a rotating star model of 15 M$_{\odot}$ rotating with an initial velocity of 300 km$\cdot$s$^{-1}$.

We compute the evolution of different abundance ratios in the Milky Way (MW) for two different sets of stellar yields. In one of them stellar rotation is taken into account and we investigate its effects on the chemical evolution model predictions. Moreover, we show that some abundance ratios offer an important tool to investigate the halo-disk discontinuity. For the first time it is shown that the effect of a halt in the star formation between the halo/thick disk and thin disk phases, already suggested from studies based both on Fe/O vs. O/H and Fe/Mg vs. Mg/H, should also be seen in a C/O versus O/H plot if C is produced mainly by low- and intermediate-mass stars (LIMS). The idea that C originates mainly from LIMS is suggested by the flat behavior of the [C/Fe] ratio as a function of metallicity, from [Fe/H]=-2.2 to solar, and by the fact that very recent C/O measurements for stars in the MW halo and disk seem to show a discontinuity around log(O/H)+12=8.4. Finally, a more gentle increase of N abundance with metallicity (or time), relative to models adopting the yields of van den Hoek and Groenewegen (1997), is predicted by using the stellar yields of Meynet and Maeder (2002 - which include stellar rotation but not hot-bottom burning) for intermediate mass stars. This fact has some implications for the timescales of N enrichment and thus for the interpretation of the nature of Damped Lyman Alpha Systems.

We examine the properties of Wolf--Rayet (WR) stars predicted by models of rotating stars taking account of the new mass loss rates for O--type stars and WR stars (Vink et al. \cite{Vink00}, \cite{Vink01}; Nugis & Lamers \cite{NuLa00}) and of the wind anisotropies induced by rotation. We find that the rotation velocities $v$ of WR stars are modest, i.e. about 50 km s$^{-1}$, not very dependant on the initial $v$ and masses. For the most massive stars, the evolution of $v$ is very strongly influenced by the values of the mass loss rates; below $\sim$12 M$_\odot$ the evolution of rotation during the MS phase and later phases is dominated by the internal coupling. Massive stars with extreme rotation may skip the LBV phase. Models having a typical $v$ for the O--type stars have WR lifetimes on the average two times longer than for non--rotating models. The increase of the WR lifetimes is mainly due to that of the H--rich eWNL phase. Rotation allows a transition WN/WC phase to be present for initial masses lower than 60 M$_\odot$. The durations of the other WR subphases are less affected by rotation. The mass threshold for forming WR stars is lowered from 37 to 22 M$_\odot$ for typical rotation. The comparisons of the predicted number ratios WR/O, WN/WC and of the number of transition WN/WC stars show very good agreement with models with rotation, while this is not the case for models with the present--day mass loss rates and no rotation. As to the chemical abundances in WR stars, rotation brings only very small changes for WN stars, since they have equilibrium CNO values. However, WC stars with rotation have on average lower C/He and O/He ratios. The luminosity distribution of WC stars is also influenced by rotation.

The dynamical shear instability is an important mixing process in the advanced stages of the evolution of massive stars. We calculated different models of 15 Mo with an initial rotational velocity, v(rot,ini)= 300 km/s to investigate its efficiency. We found that the dynamical shear instability has a timescale shorter than Oxygen burning timescale and that it slightly enlarges the convective zones and smoothens the omega gradients throughout the evolution. However, its effect is too localized to slow down the core of the star.

We use the rotating stellar models described in the paper by A. Maeder & G. Meynet in this volume to consider the effects of rotation on the evolution of the most massive stars into and during the Wolf--Rayet phase, and on the post-Main Sequence evolution of intermediate mass stars. The two main results of this discussion are the following. First, we show that rotating models are able to account for the observed properties of the Wolf--Rayet stellar populations at solar metallicity. Second, at low metallicities, the inclusion of stellar rotation in the calculation of chemical yields can lead to a longer time delay between the release of oxygen and nitrogen into the interstellar medium following an episode of star formation, since stars of lower masses (compared to non-rotating models) can synthesize primary N. Qualitatively, such an effect may be required to explain the relative abundances of N and O in extragalactic metal--poor environments, particularly at high redshifts.

For $^{12}$C and $^{16}$O, rotating models predict in general enhanced yields. At high metallicity, the carbon and oxygen yields from the very high mass stars, which go through a WR phase, are little affected by rotation. For $^{14}$N, rotation allows the production of important amounts of primary nitrogen in {\it intermediate mass stars} at very low metallicity. The process invoked for this production is different, from the classically accepted scenario {\it i.e.} the Hot Bottom Burning (HBB) in Asymptotic Giant Branch (AGB) stars. Rotating models also predict important productions of primary $^{13}$C, $^{17}$O and $^{22}$Ne at very low metallicity.

Rotation appears as a dominant effect in massive star evolution. It largely affects all the model outputs: inner structure, tracks, lifetimes, isochrones, surface compositions, blue to red supergiant ratios, etc. At lower metallicities, the effects of rotational mixing are larger; also, more stars may reach critical velocity, even if the initial distribution of rotational velocities is the same.

Gamma-ray and microwave observations of the Cygnus region reveal an intense signal of 1.809 MeV line emission, attributed to radioactive decay of 26Al, that is closely correlated with 53 GHz free-free emission, originating from the ionised interstellar medium. We modelled both emissions using a multi-wavelength evolutionary synthesis code for massive star associations that we applied to the known massive star populations in Cygnus. For all OB associations and young open clusters in the field, we determined the population age, distance, and richness as well as the uncertainties in all these quantities from published photometric and spectroscopic data. We propagate the population uncertainties in model uncertainties by means of a Bayesian method. The young globular cluster Cyg OB2 turns out to be the dominant 26Al nucleosynthesis and ionisation source in Cygnus. Our model reproduces the ionising luminosity of the Cygnus region very well, yet it underestimates 26Al production by about a factor of 2. We attribute this underestimation to shortcomings of current nucleosynthesis models, and suggest the inclusion of stellar rotation as possible mechanism to enhance 26Al production. We also modelled 60Fe nucleosynthesis in the Cygnus region, yet the small number of recent supernova events suggests only little 60Fe production. Consequently, a detection of the 1.137 MeV and 1.332 MeV decay lines of 60Fe from Cygnus by the upcoming INTEGRAL observatory is not expected.

We calculate a grid of star models with and without the effects of axial rotation for stars in the mass range between 2 and 60 M$_{\odot}$ for the metallicity $Z = 10^{-5}$. Star models with initial masses superior or equal to 9 M$_\odot$ were computed up to the end of the carbon--burning phase. Star models with masses between 2 and 7 M$_\odot$ were evolved beyond the end of the He--burning phase through a few thermal pulses during the AGB phase. Compared to models at $Z=0.02$, the low $Z$ models show faster rotating cores and stronger internal $\Omega$--gradients, which favour an important mixing of the chemical elements. In very low $Z$ models, primary nitrogen is produced during the He--burning phase by rotational diffusion of $^{12}$C into the H--burning shell. The intermediate mass stars of very low $Z$ are the main producers of primary $^{14}$N, but massive stars also contribute to this production; no significant primary nitrogen is made in models at metallicity $Z$=0.004 or above. We calculate the chemical yields in He, C, N, O and heavy elements and discuss the chemical evolution of the CNO elements at very low Z. Remarkably, the C/O vs O/H diagram is mainly sensitive to the interval of stellar masses, while the N/O vs O/H diagram is mainly sensitive to the average rotation of the stars contributing to the element synthesis. The presently available observations in these diagrams seem to favour contributions either from stars down to about 2 M$_{\odot}$ with normal rotation velocities or from stars above 8 M$_{\odot}$ but with very fast rotation.

We present new determinations of the blue to red supergiant ratio ($B/R$) in young open clusters at various metallicities. For this purpose, we examine the HR diagrams of 45 clusters in the Galaxy and of 4 clusters in the Magellanic Clouds. The identification of supergiants is based on spectroscopic measurements (with photometric counts to check the results). The new counts confirm the increase of the $B/R$ ratio when the metallicity increases with the following normalized relation : $\frac{B/R}{(B/R)_{\odot}} \cong 0.05\cdot e^{3 \frac{Z}{Z_{\odot}}}$, where $Z_{\odot}=0.02$ and $(B/R)_{\odot}$ is the value of $B/R$ at $Z_{\odot}$ which depends on the definition of $B$ and $R$ and on the age interval considered (e.g. for spectroscopic counts including clusters with $\log age$ between 6.8 and 7.5, $(B/R)_{\odot}\cong3$ when $B$ includes O, B and A supergiants).

We investigate the role of stellar axial rotation on the nitrogen nucleosynthesis at low metallicities Z. For this purpose, we have calculated models with initial masses between 2 and 60 M$_\odot$ at Z=0.00001 from the zero age sequence to the phase of thermal pulses for models below or equal to 7 M$_\odot$, and up to the end of central C--burning for the more massive stars. The models include all the main physical effects of rotation. We show that intermediate mass stars with rotation naturally reproduce the occurrence and amount of primary nitrogen in the early star generations in the Universe. We identify two reasons why rotating models at low Z produce primary $^{14}$N: 1) Since the stars lose less angular momentum, they rotate faster. Simultaneously, they are more compact, thus differential rotation and shear mixing are stronger. 2) The H--burning shell has a much higher temperature and is thus closer to the core, which favours mixing between the two.

This report discusses the possibilities for observational tests of current models and ideas regarding the origin of the Galactic Cosmic Rays using facilities which are already available.

We calculate a grid of models with and without the effects of axial rotation for massive stars in the range of 9 to 60 M$_{\odot}$ and metallicity $Z$ = 0.004 appropriate for the SMC. Remarkably, the ratios $\Omega/\Omega_{\mathrm{crit}}$ of the angular velocity to the break-up angular velocity grow strongly during the evolution of high mass stars, contrary to the situation at $Z$ = 0.020. The reason is that at low $Z$, mass loss is smaller and the removal of angular momentum during evolution much weaker, also there is an efficient outward transport of angular momentum by meridional circulation. Thus, a much larger fraction of the stars at lower $Z$ reach break-up velocities and rotation may thus be a dominant effect at low $Z$. The models with rotation well account for the long standing problem of the large numbers of red supergiants observed in low $Z$ galaxies, while current models with mass loss were predicting no red supergiants. We discuss in detail the physical effects of rotation which favour a redwards evolution in the HR diagram. The models also predict large N enrichments during the evolution of high mass stars. The predicted relative N-enrichments are larger at $Z$ lower than solar and this is in very good agreement with the observations for A-type supergiants in the SMC.

The detection of 1.809 MeV gamma-ray line emission from the Cygnus X complex by the COMPTEL telescope is one of the most convincing proves of massive star nucleosynthesis in our Galaxy. The Cygnus X complex is an extremely active nearby region, containing several star forming regions, OB associations and young open star clusters. It houses some of the most massive stars known in our Galaxy and concentrates a large number of Wolf-Rayet stars. Thermal radio continuum emission and intense [C II] line emission reveals widespread ionisation, and at least 60 individual H II regions have been identified. In order to understand the 1.809 MeV gamma-ray line emission from the Cygnus X complex, and to compare the observations to theoretical nucleosynthesis calculations, we modelled the multi-wavelength spectrum of the region by means of an evolutionary synthesis model. Our investigation leads us to the following conclusions: 1. Stellar wind ejection is the dominant mechanism for the observed 26Al enrichment in the Cygnus X region. 2. Cyg OB2 is by far the dominating massive star association in Cygnus X and 1.809 MeV line emission from 26Al produced in this association should be detectable by the spectrometer SPI on INTEGRAL. 3. There is only low supernova activity in the Cygnus X complex and the 60Fe lines should be below INTEGRAL's detection sensitivity. The detectability of an individual massive star cluster (Cyg OB2) by INTEGRAL provides a fantastic opportunity for future nucleosynthesis studies using gamma-ray line spectroscopy. We will explore the scientific potential of such observations, and provide estimates for gamma-ray line intensity distributions based on the massive star census of the Cygnus X region.

We have developed a new diagnostic tool for the study of gamma-ray emission lines from radioactive isotopes (such as 26Al and 60Fe in conjunction with other multi-wavelength observables of Galactic clusters, associations, and alike objects. Our evolutionary synthesis models are based on the code of Cervino & Mas-Hesse (1994), which has been updated to include recent stellar evolution tracks, new stellar atmospheres for OB and WR stars, and nucleosynthetic yields from massive stars during hydrostatic burning phases and explosive SN II and SN Ib events. The temporal evolution of 26Al and 60Fe production, the equivalent yield of 26Al per ionising O7 V star (yal), and other observables are predicted for a coeval population. We also employed our model to estimate the steady state nucleosynthesis of a population of solar metallicity. To allow for a fully quantitative analysis of existing and future multi-wavelength observations, we propose a Bayesian approach that allows the inclusion of IMF richness effects and observational uncertainties in the analysis. In particular, a Monte Carlo technique is adopted to estimate probability distributions for all observables of interest. We outline the procedure of exploiting these distributions by applying our model to a fictive massive star association. Applications to existing observations of the Cygnus and Vela regions will be discussed in companion papers.

Several properties of massive stars with large effects of rotation and radiation are studied. We show that there are 2 roots for the equation giving the rotational velocities at break-up: 1) The usual solution, which is shown to apply when the Eddington ratio $\Gamma$ of the star is smaller than formally 0.639. 2) Above this value of $\Gamma$, there is a second root, inferior to the first one, for the break-up velocity. This second solution tends to zero, when $\Gamma$ tends towards 1. This second root results from the interplay of radiation and rotation, and in particular from the reduction by rotation of the effective mass in the local Eddington factor. The analysis made here should hopefully clarify a recent debate between Langer (\cite{La97,La98}) and Glatzel (\cite{Gla98}). The expression for the global mass loss-rates is a function of both $\Omega$ and $\Gamma$, and this may give raise to extreme mass loss-rates ($\Omega \Gamma $-limit). In particular, for O-type stars, LBV stars, supergiants and Wolf-Rayet stars, even slow rotation may dramatically enhance the mass loss rates. Numerical examples in the range of 9 to 120 M$_{\odot}$ at various $T_\mathrm{eff}$ are given.

Grids of models for rotating stars are constructed in the range of 9 to 120 M$_{\odot}$ at solar metallicity. The following effects of rotation are included: shellular rotation, new structure equations for non-conservative case, surface distorsions, increase of mass loss with rotation, meridional circulation and interaction with horizontal turbulence, shear instability and coupling with thermal effects, advection and diffusion of angular momentum treated in the non-stationary regime, transport and diffusion of the chemical elements. Globally we find that for massive stars the effects of rotation have an importance comparable to those of mass loss. For an average rotation, the tracks in the HR diagram are modified like a moderate overshoot would do. In general, an average rotation may increase the MS lifetime up to about 30 %; for the helium-burning phase the effects are smaller and amount to at most 10 %. Rotation also brings significant surface He- and N-enhancements, they are higher for higher masses and rotation.

First, we review the main physical effects to be considered in the building of evolutionary models of rotating stars on the Upper Main-Sequence (MS). The internal rotation law evolves as a result of contraction and expansion, meridional circulation, diffusion processes and mass loss. In turn, differential rotation and mixing exert a feedback on circulation and diffusion, so that a consistent treatment is necessary. We review recent results on the evolution of internal rotation and the surface rotational velocities for stars on the Upper MS, for red giants, supergiants and W-R stars. A fast rotation is enhancing the mass loss by stellar winds and reciprocally high mass loss is removing a lot of angular momentum. The problem of the ``break-up'' or $\Omega$-limit is critically examined in connection with the origin of Be and LBV stars. The effects of rotation on the tracks in the HR diagram, the lifetimes, the isochrones, the blue to red supergiant ratios, the formation of W-R stars, the chemical abundances in massive stars as well as in red giants and AGB stars, are reviewed in relation to recent observations for stars in the Galaxy and Magellanic Clouds. The effects of rotation on the final stages and on the chemical yields are examined, as well as the constraints placed by the periods of pulsars. On the whole, this review points out that stellar evolution is not only a function of mass M and metallicity Z, but of angular velocity $\Omega$ as well.

Meynet and Arnould (1993) have suggested that Wolf-Rayet (WR) stars could significantly contaminate the Galaxy with 19F. In their scenario, 19F is synthesized at the beginning of the He-burning phase from the 14N left over by the previous CNO-burning core, and is ejected in the interstellar medium when the star enters its WC phase. Recourse to CNO seeds makes the 19F yields metallicity-dependent. These yields are calculated on grounds of detailed stellar evolutionary sequences for an extended range of initial masses (from 25 to 120 Msol) and metallicities (Z = 0.008, 0.02 and 0.04). The adopted mass loss rate prescription enables to account for the observed variations of WR populations in different environments. The 19F abundance in the WR winds of 60 Msol model stars is found to be about 10 to 70 times higher than its initial value, depending on the metallicity. This prediction is used in conjunction with a very simple model for the chemical evolution of the Galaxy to predict that WR stars could be significant (dominant?) contributors to the solar system fluorine content. We also briefly discuss the implications of our model on the possible detection of fluorine at high redshift.

Radioisotopes are natural clocks which can be used to estimate the age of the solar system. They also influence the shape of supernova light curves. In addition, the diffuse emission at 1.8 MeV from the decay of 26Al may provide a measure of the present day nucleosynthetic activity in the Galaxy. Therefore, even if radionuclides represent only a tiny fraction of the cosmic matter, they carry a unique piece of information. A large number of radioisotopes are produced by massive stars at the time of their supernova explosion. A more or less substantial fraction of them are also synthesized during the previous hydrostatic burning phases. These nuclides are then ejected either at the time of the supernova event, or through stellar winds during their hydrostatic burning phases. This paper focusses of the non explosive ejection of radionuclides by non-rotating or rotating Wolf-Rayet stars.

Radionuclides with half-lives ranging from some years to billions of years presumably synthesized outside of the solar system are now recorded in `live' or `fossil' form in various types of materials, like meteorites or the galactic cosmic rays. They bring specific astrophysical messages the deciphering of which is briefly reviewed here, with special emphasis on the contribution of Jerry Wasserburg and his collaborators to this exciting field of research. A brief mention is also made to $\gamma$-ray line astrophysics. It provides an interesting information on live short-lived radionuclides in the present interstellar medium, and thus complements in a very important way the study of extinct radionuclides in meteorites.

We present observations with the Hubble Space Telescope Wide Field Planetary Camera 2 of three fields centered on super-metal-rich globular clusters in the bulge of M31. Our (I,V-I) color-magnitude diagrams reach as faint as I ~ 26.5 mag and clearly reveal the magnitude of the first ascent red giant branch (RGB) tip. We find that the apparent I magnitude of the RGB tip does not become brighter near the center of M31 as concluded by previous investigators. Our observations and artificial star experiments presented in this study strongly support the idea that previous very bright stars were likely the result of spurious detections of blended stars due to crowding in lower resolution images. On the contrary, our observations indicate that, at a mean projected galactocentric distance of 1.1 kpc, the RGB tip is some 1.3 magnitudes fainter than it is at 7 kpc. An analysis of this difference in RGB tip magnitude suggests that the M31 bulge stellar population has a mean metallicity close to that of the Sun.

We present stellar evolutionary models covering the mass range from 0.4 to 1 Msun calculated for metallicities Z=0.020 and 0.001 with the MHD equation of state (Hummer & Mihalas, 1988; Mihalas et al. 1988; D\"appen et al. 1988). A parallel calculation using the OPAL (Rogers et al. 1996) equation of state has been made to demonstrate the adequacy of the MHD equation of state in the range of 1.0 to 0.8 Msun (the lower end of the OPAL tables). Below, down to 0.4 Msun, we have justified the use of the MHD equation of state by theoretical arguments and the findings of Chabrier & Baraffe (1997). We use the radiative opacities by Iglesias & Rogers (1996), completed with the atomic and molecular opacities by Alexander & Fergusson (1994). We follow the evolution from the Hayashi fully convective configuration up to the red giant tip for the most massive stars, and up to an age of 20 Gyr for the less massive ones. We compare our solar-metallicity models with recent models computed by other groups and with observations. The present stellar models complete the set of grids computed with the same up-to-date input physics by the Geneva group [Z=0.020 and 0.001, Schaller et al. (1992), Bernasconi (1996), and Charbonnel et al. (1996); Z=0.008, Schaerer et al. (1992); Z=0.004, Charbonnel et al. (1993); Z=0.040, Schaerer et al. (1993); Z=0.10, Mowlavi et al. (1998); enhanced mass loss rate evolutionary tracks, Meynet et al. (1994)].

The F8Ia supergiant V810 Centauri is part of a long-term high-precision photometric monitoring program on long period variables started twenty years ago. Time series analysis of this unique set of 500 data points, spanning almost fifteen years in the homogeneous Geneva photometric system, is presented. Cluster membership, physical parameters and evolutionary status of the star are reinvestigated. Radial velocity data do not support the cluster membership to Stock 14}. Ultraviolet and optical spectrophotometry is combined with optical and infrared photometry to evaluate the physical parameters of the yellow supergiant (Teff = 5970 K, M_bol = -8.5, R = 420 R_sun) and of its B0III companion. From theoretical stellar evolutionary tracks, an initial mass of 25 M_sun is estimated for V810 Cen, which is actually at the end of its first redward evolution. V810 Cen is a multi-periodic small amplitude variable star, whose amplitudes are variable with time. The period of the main mode, 156 d, is in agreement with the Period--Luminosity--Colour relation for supergiants. This mode is most probably the fundamental radial one. According to the theoretical pulsation periods for the radial modes, calculated from a linear non-adiabatic analysis, the period of the observed second mode, 107 d, is much too long to correspond to the first radial overtone . Thus, this second mode could be a non-radial p-mode. Other transient periods are observed, in particular at 187 d. The l ength of this period suggests a non-radial g-mode. Then, the complex variability of V810 Cen could be due to a mixing of unstable radial and non-radial p- and g-modes.

An analysis of some properties of stellar models as a function of metallicity Z (and helium content Y) is presented, with special attention to those stars with metallicities higher than twice or three times solar. It is shown that the stellar properties as a function of Z are mainly determined by the effects of the opacities at sub-solar metallicities, and by the effects of the mean molecular weight and stellar mass loss at higher metallicities. As a result, very metal-rich stars (Z >~0.05) exhibit properties that deviate from what is expected from the known characteristics at lower metallicities. In particular, they are more luminous and hotter than those at Z <~ 0.05 due to the effect of the mean molecular weight. They have main sequence lifetimes much shorter (60 % shorter at Z=0.1 than at Z=0.02) than those at solar metallicity due to their lower initial hydrogen content. Finally, the high mass loss rates at high metallicities affect significantly the population synthesis of massive stars in very metal-rich regions. An analysis of expected Wolf-Rayet and supernovae populations in such conditions is briefly presented.

We present spectroscopic observations of the central star clusters in NGC 5253 the aim of which is to search for WC stars. Our observations show the presence of Wolf-Rayet (WR) stars not only of WN but also of WC subtype in two star forming regions corresponding to the maximum optical and UV emission. The massive star population we derive is consistent with young bursts of ~3 and 4 Myr. The region of maximum optical emission is found to provide the dominant contribution of the ionizing flux as opposed to the less extinguished region of maximum UV brightness. The presence of WR stars near the N-enriched regions found by Walsh & Roy (1987, 1989) and Kobulnicky et al. (1997) suggests they are a possible source of N. It is presently unclear whether or not our detection of WC stars is compatible with the normal observed He/O and C/O abundance ratios.

We describe the main-sequence evolution of a rotating 9 $M_\odot$ star. Its interior rotation profile is determined by the redistribution of angular momentum through the meridian circulation and through the shear turbulence generated by the differential rotation; the possible effect of internal waves is neglected. We examine the mixing of chemicals produced by the same internal motions. Our modelization is based on the set of equations established by Zahn (1992) and completed in Matias, Talon & Zahn (1996). Our calculations show that the amount of mixing associated with a typical rotation velocity of $\sim 100 km/s$ yields stellar models whose global parameters are very similar to those obtained with the moderate overshooting ($d/H_P \simeq 0.2$) which has been invoked until now to fit the observations. Fast rotation ($\sim 300 km/s$) leads to significant changes of the C/N and O/N surface ratios, but the abundance of He is barely increased. The modifications of the internal composition profile due to such rotational mixing will certainly affect the post--main-sequence evolution.