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

The competing $^{22}$Ne($\alpha,\gamma$)$^{26}$Mg and $^{22}$Ne($\alpha,n$)$^{25}$Mg reactions control the production of neutrons for the weak $s$-process in massive and AGB stars. In both systems, the ratio between the corresponding reaction rates strongly impacts the total neutron budget and strongly influences the final nucleosynthesis. The $^{22}$Ne($\alpha,\gamma$)$^{26}$Mg and $^{22}$Ne($\alpha,n$)$^{25}$Mg reaction rates was re-evaluated by using newly available information on $^{26}$Mg given by various recent experimental studies. Evaluations of The evaluated $^{22}$Ne($\alpha,\gamma$)$^{26}$Mg reaction rate remains substantially similar to that of Longland {\it et al.} but, including recent results from Texas A\&M, the $^{22}$Ne($\alpha,n$)$^{25}$Mg reaction rate is lower at a range of astrophysically important temperatures. Stellar models computed with NEWTON and MESA predict decreased production of the weak branch $s$-process due to the decreased efficiency of $^{22}$Ne as a neutron source. Using the new reaction rates in the MESA model results in $^{96}$Zr/$^{94}$Zr and $^{135}$Ba/$^{136}$Ba ratios in much better agreement with the measured ratios from presolar SiC grains.

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.

A 70Msun BH was discovered in Milky Way disk in a long period and almost circular detached binary system (LB-1) with a high metallicity 8Msun B star companion. Current consensus on the formation of BHs from high metallicity stars limits the black hole mass to be below 20Msun. Using simple evolutionary model, we show that the formation of a 70Msun BH in high metallicity environment is possible if stellar wind mass loss rates are reduced by factor of 5. As observations indicate, a fraction of massive stars have surface magnetic fields which may quench the wind mass-loss, independently of stellar mass and metallicity. We also computed detailed stellar evolution models and we confirm such a scenario. A non-rotating 85Msun model at Z=0.014 with decreased winds ends up as a 71Msun star prior core-collapse with a 32Msun helium core and a 28Msun CO core. Such star avoids pair-instability pulsation supernova mass loss and may form a 70Msun BH in the direct collapse. Stars that can form such BHs expand to significant size with radius of R>600Rsun, exceeding the size of LB-1 orbit. Therefore, we can explain the formation of BHs upto 70Msun at high metallicity and this result is independent from LB-1. However, if LB-1 hosts a massive BH we are unable to explain how such a binary star system could have formed without invoking some exotic scenarios.

The first massive stars triggered the onset of chemical evolution by releasing the first metals (elements heavier than helium) in the Universe. The nature of these stars and how the early chemical enrichment took place is still largely unknown. Rotational-induced mixing in the stellar interior can impact the nucleosynthesis during the stellar life of massive stars and lead to stellar ejecta having various chemical compositions. We present low and zero-metallicity 20, 25 and 40 $M_{\odot}$ stellar models with various initial rotation rates and assumptions for the nuclear reactions rates. With increasing initial rotation, the yields of light (from $\sim$ C to Al) and trans-iron elements are boosted. The trans-iron elements (especially elements heavier than Ba) are significantly affected by the nuclear reaction uncertainties. The chemical composition of the observed CEMP (carbon-enhanced metal-poor) stars CS29528-028 and HE0336+0113 are consistent with the chemical composition of the material ejected by a fast rotating 40~$M_{\odot}$ model.

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

The $\nu p$ process appears in proton-rich, hot matter which is expanding in a neutrino wind and may be realised in explosive environments such as core-collapse supernovae or in outflows from accretion disks. The impact of uncertainties in nuclear reaction cross sections on the finally produced abundances has been studied by applying Monte Carlo variation of all astrophysical reaction rates in a large reaction network. As the detailed astrophysical conditions of the $\nu p$ process still are unknown, a parameter study was performed, with 23 trajectories covering a large range of entropies and $Y_\mathrm{e}$. The resulting abundance uncertainties are given for each trajectory. The $\nu p$ process has been speculated to contribute to the light $p$ nuclides but it was not possible so far to reproduce the solar isotope ratios. It is found that it is possible to reproduce the solar $^{92}$Mo/$^{94}$Mo abundance ratio within nuclear uncertainties, even within a single trajectory. The solar values of the abundances in the Kr-Sr region relative to the Mo region, however, cannot be achieved within a single trajectory. They may still be obtained from a weighted superposition of different trajectories, though, depending on the actual conditions in the production site. For a stronger constraint of the required conditions, it would be necessary to reduce the uncertainties in the 3$\alpha$ and $^{56}$Ni(n,p)$^{56}$Co rates at temperatures $T>3$ GK.

Most neutron capture elements have a double production by r- and s-processes, but the question of production sites is complex and still open. Recent studies show that including stellar rotation can have a deep impact on nucleosynthesis. We studied the evolution of Sr and Ba in the Milky Way. A chemical evolution model was employed to reproduce the Galactic enrichment. We tested two different nucleosynthesis prescriptions for s-process in massive stars, adopted from the Geneva group and the Rome group. Rotation was taken into account, studying the effects of stars without rotation or rotating with different velocities. We also tested different production sites for the r-process: magneto rotational driven supernovae and neutron star mergers. The evolution of the abundances of Sr and Ba is well reproduced. The comparison with the most recent observations shows that stellar rotation is a good assumption, but excessive velocities result in overproduction of these elements. In particular, the predicted evolution of the [Sr/Ba] ratio at low metallicity does not explain the data at best if rotation is not included. Adopting different rotational velocities for different stellar mass and metallicity better explains the observed trends. Despite the differences between the two sets of adopted stellar models, both show a better agreement with the data assuming an increase of rotational velocity toward low metallicity. Assuming different r-process sources does not alter this conclusion.

Aims. We investigate the s-process during the AGB phase of stellar models whose cores are enforced to rotate at rates consistent with asteroseismology observations of their progenitors and successors. Methods. We calculated new 2M$_{\odot}$, Z=0.01 models, rotating at 0, 125, and 250 km/s at the start of main sequence. An artificial, additional viscosity was added to enhance the transport of angular momentum in order to reduce the core rotation rates to be in agreement with asteroseismology observations. We compared rotation rates of our models with observed rotation rates during the MS up to the end of core He burning, and the white dwarf phase. Results. We present nucleosynthesis calculations for these rotating AGB models that were enforced to match the asteroseismic constraints on rotation rates of MS, RGB, He-burning, and WD stars. In particular, we calculated one model that matches the upper limit of observed rotation rates of core He-burning stars and we also included a model that rotates one order of magnitude faster than the upper limit of the observations. The s-process production in both of these models is comparable to that of non-rotating models. Conclusions. Slowing down the core rotation rate in stars to match the above mentioned asteroseismic constraints reduces the rotationally induced mixing processes to the point that they have no effect on the s-process nucleosynthesis. This result is independent of the initial rotation rate of the stellar evolution model. However, there are uncertainties remaining in the treatment of rotation in stellar evolution, which need to be reduced in order to confirm our conclusions, including the physical nature of our approach to reduce the core rotation rates of our models, and magnetic processes.

The propagation of uncertainties in reaction cross sections and rates of neutron-, proton-, and alpha-induced reactions into the final isotopic abundances obtained in nucleosynthesis models is an important issue in studies of nucleosynthesis and Galactic Chemical Evolution. We developed a Monte Carlo method to allow large-scale postprocessing studies of the impact of nuclear uncertainties on nucleosynthesis. Temperature-dependent rate uncertainties combining realistic experimental and theoretical uncertainties are used. From detailed statistical analyses uncertainties in the final abundances are derived as probability density distributions. Furthermore, based on rate and abundance correlations an automated procedure identifies the most important reactions in complex flow patterns from superposition of many zones or tracers. The method so far was already applied to a number of nucleosynthesis processes. Here we focus on the production of p-nuclei in white dwarfs exploding as thermonuclear (type Ia) supernovae. We find generally small uncertainties in the final abundances despite of the dominance of theoretical nuclear uncertainties. A separate analysis of low- and high-density regions indicates that the total uncertainties are dominated by the high-density regions.

It has been suggested that a $\nu$p process can occur when hot, dense, and proton-rich matter is expanding within a strong flux of anti-neutrinos. In such an environment, proton-rich nuclides can be produced in sequences of proton captures and (n,p) reactions, where the free neutrons are created in situ by $\overline{\nu}_\mathrm{e}+\mathrm{p} \rightarrow \mathrm{n}+\mathrm{e}^+$ reactions. The detailed hydrodynamic evolution determines where the nucleosynthesis path turns off from N = Z line and how far up the nuclear chart it runs. In this work, the uncertainties on the final isotopic abundances stemming from uncertainties in the nuclear reaction rates were investigated in a large-scale Monte Carlo approach, simultaneously varying ten thousand reactions. A large range of model conditions was investigated because a definitive astrophysical site for the $\nu$p process has not yet been identified. The present parameter study provides, for each model, identification of the key nuclear reactions dominating the uncertainty for a given nuclide abundance. As all rates appearing in the $\nu$p process involve unstable nuclei, and thus only theoretical rates are available, the final abundance uncertainties are larger than those for nucleosynthesis processes closer to stability. Nevertheless, most uncertainties remain below a factor of three in trajectories with robust nucleosynthesis. More extreme conditions allow production of heavier nuclides but show larger uncertainties because of the accumulation of the uncertainties in many rates and because the termination of nucleosynthesis is not at equilibrium conditions. It is also found that the solar ratio of the abundances of ${}^{92}$Mo and ${}^{94}$Mo could be reproduced within uncertainties.

The production of the neutron-capture isotopes beyond iron that we observe today in the solar system is the result of the combined contribution of the r-process, the s- process and possibly the i-process. Low-mass AGB (2 < M/Msun < 3) and massive (M >10 Msun ) stars have been identified as the sites of the s-process. In this work we consider the evolution and nucleosynthesis of low-mass AGB stars. We provide an update of the NuGrid Set models, adopting the same general physics assumptions but using an updated convective-boundary mixing model accounting for the contribution from internal gravity waves. The combined data set includes the initial masses Mzams/Msun = 2, 3 for Z = 0.03, 0.02, 0.01. These models are computed with the MESA stellar code and the evolution is followed up to the end of the AGB phase. The nucleosynthesis was calculated for all isotopes in post-processing with the NuGrid mppnp code. The convective boundary mixing model leads to the formation of a 13C-pocket three times wider compared to the one obtained in the previous set of models, bringing the simulation results now in closer agreement with observations. We also discuss the potential impact of other processes inducing mixing, like rotation, adopting parametric models compatible with theory and observations. Complete yield data tables, derived data products and online analytic data access are provided.

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.

Transport of angular momentum has been a challenging topic within the stellar evolution community, even more since the recent asteroseismic surveys. All published studies on rotation using asteroseismic observations show a discrepancy between the observed and calculated rotation rates, indicating there is an undetermined process of angular momentum transport active in these stars. We aim to constrain the efficiency of this process by investigating rotation rates of 2.5 M$_{\odot}$ stars. First, we investigated whether the Tayler-Spruit dynamo could be responsible for the extra transport of angular momentum for stars with an initial mass of 2.5 M$_{\odot}$. Then, by computing rotating models including a constant additional artificial viscosity, we determined the efficiency of the missing process of angular momentum transport by comparing the models to the asteroseismic observations of core helium burning stars. Parameter studies were performed to investigate the effect of the stellar evolution code used, initial mass, and evolutionary stage. We evolved our models into the white dwarf phase, and provide a comparison to white dwarf rotation rates. The Tayler-Spruit dynamo is unable to provide enough transport of angular momentum to reach the observed values of the core helium burning stars investigated in this paper. We find that a value for the additional artificial viscosity $\nu_{\rm{add}}$ around 10$^7$ cm$^2$ s$^{-1}$ provides enough transport of angular momentum. However, the rotational period of these models is too high in the white dwarf phase to match the white dwarf observations. From this comparison we infer that the efficiency of the missing process must decrease during the core helium burning phase. When excluding the $\nu_{\rm{add}}$ during core helium burning phase, we can match the rotational periods of both the core helium burning stars and white dwarfs.

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.

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

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

We investigated the impact of uncertainties in neutron-capture and weak reactions (on heavy elements) on the s-process nucleosynthesis in low-mass stars using a Monte-Carlo based approach. We performed extensive nuclear reaction network calculations that include newly evaluated temperature-dependent upper and lower limits for the individual reaction rates. Our sophisticated approach is able to evaluate the reactions that impact more significantly the final abundances. We found that beta-decay rate uncertainties affect typically nuclides near s-process branchings, whereas most of the uncertainty in the final abundances is caused by uncertainties in neutron capture rates, either directly producing or destroying the nuclide of interest. Combined total nuclear uncertainties due to reactions on heavy elements are approximately 50%.

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

We provide here a significant extension of the NuGrid Set 1 models in mass coverage and toward lower metallicity, adopting the same physics assumptions. The combined data set now includes the initial masses M/Msun = 1, 1.65, 2, 3, 4, 5, 6, 7, 12, 15, 20, 25 for Z = 0.02, 0.01, 0.006, 0.001, 0.0001 with alpha-enhanced composition for the lowest three metallicities. These models are computed with the MESA stellar evolution code and are evolved up to the AGB, the white dwarf stage, or until core collapse. The nucleosynthesis was calculated for all isotopes in post-processing with the NuGrid mppnp code. Explosive nucleosynthesis is based on semi-analytic 1D shock models. Metallicity-dependent mass loss, convective boundary mixing in low- and intermediate mass models and H and He core burning massive star models is included. Convective O-C shell mergers in some stellar models lead to the strong production of odd-Z elements P, Cl, K and Sc. In AGB models with hot dredge-up the convective boundary mixing efficiency is reduced to accommodate for its energetic feedback. In both low-mass and massive star models at the lowest metallicity H-ingestion events are observed and lead to i-process nucleosynthesis and substantial N-15 production. Complete yield data tables, derived data products and online analytic data access are provided.

Several thousand tracers from a 2D model of a thermonuclear supernova were used in a Monte Carlo post-processing approach to determine p-nuclide abundance uncertainties originating from nuclear physics uncertainties in the reaction rates.

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.

The main s-process taking place in low mass stars produces about half of the elements heavier than iron. It is therefore very important to determine the importance and impact of nuclear physics uncertainties on this process. We have performed extensive nuclear reaction network calculations using individual and temperature-dependent uncertainties for reactions involving elements heavier than iron, within a Monte Carlo framework. Using this technique, we determined the uncertainty in the main s-process abundance predictions due to nuclear uncertainties link to weak interactions and neutron captures on elements heavier than iron. We also identified the key nuclear reactions dominating these uncertainties. We found that $\beta$-decay rate uncertainties affect only a few nuclides near s-process branchings, whereas most of the uncertainty in the final abundances is caused by uncertainties in neutron capture rates, either directly producing or destroying the nuclide of interest. Combined total nuclear uncertainties due to reactions on heavy elements are in general small (less than 50%). Three key reactions, nevertheless, stand out because they significantly affect the uncertainties of a large number of nuclides. These are $^{56}$Fe(n,$\gamma$), $^{64}$Ni(n,$\gamma$), and $^{138}$Ba(n,$\gamma$). We discuss the prospect of reducing uncertainties in the key reactions identified in this study with future experiments.

The recently discovered bright type II supernova OGLE14-073 evolved very slowly. The light curve rose to maximum for 90 days from discovery and then declined at a rate compatible with the radioactive decay of 56Co. In this study, we show that a pair-instability supernova is a plausible mechanism for this event. We calculate explosion models and light curves with the radiation hydrodynamics code STELLA starting from two M(ZAMS) = 150 solar masses, Z=0.001 progenitors. We obtain satisfactory fits to OGLE14-073 broadband light curves by including additional 56Ni in the centre of the models and mixing hydrogen down into the inner layers of the ejecta to a radial mass coordinate of 10 solar masses. The extra 56Ni required points to a slightly more massive progenitor star. The mixing of hydrogen could be due to large scale mixing during the explosion. We also present synthetic spectra for our models simulated with the Monte Carlo radiative transfer code ARTIS. The synthetic spectra reproduce the main features of the observed spectra of OGLE14-073. We conclude that OGLE14-073 is one of the most promising candidates for a pair-instability explosion.

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

The s-process in massive stars, producing nuclei up to $A\approx 90$, has a different behaviour at low metallicity if stellar rotation is significant. This enhanced s-process is distinct from the s-process in massive stars around solar metallicity, and details of the nucleosynthesis are poorly known. We investigated nuclear physics uncertainties in the enhanced s-process in metal-poor stars within a Monte-Carlo framework. We applied temperature-dependent uncertainties of reaction rates, distinguishing contributions from the ground state and from excited states. We found that the final abundance of several isotopes shows uncertainties larger than a factor of 2, mostly due to the neutron capture uncertainties. A few nuclei around branching points are affected by uncertainties in the $\beta$-decay.

The s-process, a production mechanism based on slow-neutron capture during stellar evolution, is the origin of about half the elements heavier than iron. Abundance predictions for s-process nucleosynthesis depend strongly on the relevant neutron-capture and $\beta$-decay rates, as well as on the details of the stellar model being considered. Here, we have used a Monte-Carlo approach to evaluate the nuclear uncertainty in s-process nucleosynthesis. We considered the helium burning of massive stars for the weak s-process and low-mass asymptotic-giant-branch stars for the main s-process. Our calculations include a realistic and general prescription for the temperature dependent uncertainty for the reaction cross sections. We find that the adopted uncertainty for (${\rm n},\gamma$) rates, tens of per cent on average, effects the production of s-process nuclei along the line of $\beta$-stability, and that the uncertainties in $\beta$-decay from excited state contributions, has the strongest impact on branching points.

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

Thermonuclear supernovae originating from the explosion of a white dwarf accreting mass from a companion star have been suggested as a site for the production of $p$ nuclides. Such nuclei are produced during the explosion, in layers enriched with seed nuclei coming from prior strong $s$ processing. These seeds are transformed to proton-richer isotopes mainly by photodisintegration reactions. Several thousand trajectories from a 2D explosion model were used in a Monte Carlo approach. Temperature-dependent uncertainties were assigned individually to thousands of rates varied simultaneously in post-processing in an extended nuclear reaction network. The uncertainties in the final nuclear abundances originating from uncertainties in the astrophysical reaction rates were determined. In addition to the 35 classical $p$ nuclides, abundance uncertainties were also determined for the radioactive nuclides $^{92}$Nb, $^{97,98}$Tc, $^{146}$Sm, and for the abundance ratios $Y$(${}^{92}$Mo)/$Y$(${}^{94}$Mo), $Y$(${}^{92}$Nb)/$Y$(${}^{92}$Mo), $Y$(${}^{97}$Tc)/$Y$(${}^{98}$Ru), $Y$(${}^{98}$Tc)/$Y$(${}^{98}$Ru), and $Y$(${}^{146}$Sm)/$Y$(${}^{144}$Sm), important for Galactic Chemical Evolution studies. Uncertainties found were generally lower than a factor of two, although most nucleosynthesis flows mainly involve predicted rates with larger uncertainties. The main contribution to the total uncertainties comes from a group of trajectories with high peak density originating from the interior of the exploding white dwarf. The distinction between low-density and high-density trajectories allows more general conclusions to be drawn, also applicable to other simulations of white dwarf explosions.

Recently the search for the oldest stars have started to focus on the Bulge region. The Galactic bulge hosts extremely old stars, with ages compatible with the ages of the oldest halo stars. The data coming from these recent observations present new chemical signatures and therefore provide complementary constraints to those already found in the halo. So, the study of the oldest bulge stars can improve dramatically the constraints on the nature of first stars and how they polluted the pristine ISM of our Galaxy. We present our first results regarding the light elements (CNO) and the neutron capture elements. Our findings in the oldest bulge stars support the scenario where the first stellar generations have been fast rotators.

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.

Uncertainties in nucleosynthesis models originating from uncertainties in astrophysical reaction rates were estimated in a Monte Carlo variation procedure. Thousands of rates were simultaneously varied within individual, temperature-dependent errors to calculate their combined effect on final abundances. After a presentation of the method, results from application to three different nucleosynthesis processes are shown: the $\gamma$-process and the s-process in massive stars, and the main s-process in AGB stars (preliminary results). Thermal excitation of nuclei in the stellar plasma and the combined action of several reactions increase the final uncertainties above the level of the experimental errors. The total uncertainty, on the other hand, remains within a factor of two even in processes involving a large number of unmeasured rates, with some notable exceptions for nuclides whose production is spread over several stellar layers and for s-process branchings.

In recent years, the viability of the pair-instability supernova (PISN) scenario for explaining superluminous supernovae has all but disappeared except for a few slowly-evolving examples. However, PISN are not predicted to be superluminous throughout the bulk of their mass range. In fact, it is more likely that the first PISN we see (if we have not seen one already) will not be superluminous. Here, we present hydrodynamic simulations of PISNe for four stellar models with unique envelope properties spanning the PISN mass range. In addition, we compute synthetic light curves for comparison with current and future observations. We also investigate, in the context of our most massive model, the prospect of mixing in the supernova ejecta alleviating discrepancies between current PISN models and the remaining superluminous candidate events. To this end, we present the first published 3D hydrodynamic simulations of PISNe. After achieving convergence between 1D, 2D, and 3D simulations we examine mixing in the supernova ejecta and its affect on the bolometric light curve. We observe slight deviations from spherical symmetry which increase with the number of dimensions. We find no significant effects on the bolometric light curve, however we conclude that mixing between the silicon and oxygen rich layers caused by the Rayleigh-Taylor instability may affect spectra.

We investigated the impact of uncertainties in neutron-capture and weak reactions (on heavy elements) on the s-process nucleosynthesis in low-mass stars and massive stars using a Monte-Carlo based approach. We performed extensive nuclear reaction network calculations that include newly evaluated temperature-dependent upper and lower limits for the individual reaction rates. We found beta-decay rate uncertainties affect only a few nuclei near s-process branchings, whereas most of the uncertainty in the final abundances is caused by uncertainties in the neutron capture rates. We suggest a list of uncertain rates as candidates for improved measurement by future experiments.

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

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

The $s$-process in massive stars produces the weak component of the $s$-process (nuclei up to $A \sim 90$), in amounts that match solar abundances. For heavier isotopes, such as barium, production through neutron capture is significantly enhanced in very metal-poor stars with fast rotation. However, detailed theoretical predictions for the resulting final $s$-process abundances have important uncertainties caused both by the underlying uncertainties in the nuclear physics (principally neutron capture reaction and $\beta$-decay rates) as well as by the stellar evolution modeling. In this work, we investigated the impact of nuclear-physics uncertainties relevant to the $s$-process in massive stars. Using a Monte-Carlo based approach, we performed extensive nuclear reaction network calculations that include newly evaluated upper and lower limits for the individual temperature dependent reaction rates. We found that most of the uncertainty in the final abundances is caused by uncertainties in the neutron capture rates, while $\beta$-decay rate uncertainties affect only a few nuclei near $s$-process branchings. The $s$-process in rotating metal-poor stars shows quantitatively different uncertainties and key reactions, although the qualitative characteristics are similar. We confirmed that our results do not significantly change at different metallicities for fast rotating massive stars in the very low metallicity regime. We highlight which of the identified key reactions are realistic candidates for improved measurement by future experiments.

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

We present a chemical abundance analysis of a metal-poor star, ROA 276, in the stellar system omega Centauri. We confirm that this star has an unusually high [Sr/Ba] abundance ratio. Additionally, ROA 276 exhibits remarkably high abundance ratios, [X/Fe], for all elements from Cu to Mo along with normal abundance ratios for the elements from Ba to Pb. The chemical abundance pattern of ROA 276, relative to a primordial omega Cen star ROA 46, is best fit by a fast-rotating low-metallicity massive stellar model of 20 Msun, [Fe/H] = -1.8, and an initial rotation 0.4 times the critical value; no other nucleosynthetic source can match the neutron-capture element distribution. ROA 276 arguably offers the most definitive proof to date that fast-rotating massive stars contributed to the production of heavy elements in the early Universe.

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

We evaluated the uncertainty relevant to s-process nucleosynthesis using a Monte-Carlo centred approach. We are based on a realistic and general prescription of temperature dependent uncertainty for the reactions. We considered massive stars for the weak s-process and AGB stars for the main s-process. We found that the adopted uncertainty for (n,$\gamma$) rates, tens of per cent on average, affect the production of s-process nuclei along the $\beta$-stability line, while for $\beta$-decay, for which contributions from excited states enhances the uncertainty, has the strongest impact on branching points.

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

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

Proton-rich nuclei, the so-called p-nuclei, are made in photodisintegration processes in outer shells of massive stars in the course of the final supernova explosion. Nuclear uncertainties in the production of these nuclei have been quantified in a Monte Carlo procedure. Bespoke temperature-dependent uncertainties were assigned to different types of reactions involving nuclei from Fe to Bi and all rates were varied randomly within the uncertainties. The resulting total production uncertainties of the p-nuclei are below a factor of two, with few exceptions. Key reactions dominating the final uncertainties have been identified in an automated procedure using correlations between rate and abundance uncertainties. Our results are compared to those of a previous study manually varying reaction rates.

With an increasing number of superluminous supernovae (SLSNe) discovered the question of their origin remains open and causes heated debates in the supernova community. Currently, there are three proposed mechanisms for SLSNe: (1) pair-instability supernovae (PISN), (2) magnetar-driven supernovae, and (3) models in which the supernova ejecta interacts with a circumstellar material ejected before the explosion. Based on current observations of SLSNe, the PISN origin has been disfavoured for a number of reasons. Many PISN models provide overly broad light curves and too reddened spectra, because of massive ejecta and a high amount of nickel. In the current study we re-examine PISN properties using progenitor models computed with the GENEC code. We calculate supernova explosions with FLASH and light curve evolution with the radiation hydrodynamics code STELLA. We find that high-mass models (200 and 250 solar masses) at relatively high metallicity (Z=0.001) do not retain hydrogen in the outer layers and produce relatively fast evolving PISNe Type I and might be suitable to explain some SLSNe. We also investigate uncertainties in light curve modelling due to codes, opacities, the nickel-bubble effect and progenitor structure and composition.

Since the emergence of the new class of extremely bright transients, super-luminous supernovae (SLSNe), three main mechanisms to power their light curves (LCs) have been discussed. They are the spin-down of a magnetar, interaction with circumstellar material, and the decay of large amounts of radioactive nickel in pair-instability supernovae (PISNe). Given the high degree of diversity seen within the class, it is possible that all three mechanisms are at work. PISN models can be self- consistently simulated from the main sequence phase of very massive stars (VMS) through to their explosion. These models robustly predict large amounts of radioactive nickel and thus very luminous SN events. However, PISN model LCs evolve more slowly than even the slowest evolving SLSNe. Multidimensional effects on the ejecta structure, specifically the mixing of radioactive nickel out to large radii, could alleviate this discrepancy with observation. Here we explore the multidimensional effects on the LC evolution by simulating the explosion phase in 1D, 2D, and 3D. We find that the ejecta from the multidimensional simulations have slightly shallower abundance gradients due to mixing at shell boundaries. We compute synthetic LCs whose shapes show no discernible differences due to the multidimensional effects.

Nuclear uncertainties in the production of $p$ nuclei in massive stars have been quantified in a Monte Carlo procedure. Bespoke temperature-dependent uncertainties were assigned to different types of reactions involving nuclei from Fe to Bi. Their simultaneous impact was studied in postprocessing explosive trajectories for three different stellar models. It was found that the grid of mass zones in the model of a 25 $M_\odot$ star, which is widely used for investigations of $p$ nucleosynthesis, is too crude to properly resolve the detailed temperature changes required for describing the production of $p$ nuclei. Using models with finer grids for 15 $M_\odot$ and 25 $M_\odot$ stars with initial solar metallicity, it was found that most of the production uncertainties introduced by nuclear reaction uncertainties are smaller than a factor of two. Since a large number of rates were varied at the same time in the Monte Carlo procedure, possible cancellation effects of several uncertainties could be taken into account. Key rates were identified for each $p$ nucleus, which provide the dominant contribution to the production uncertainty. These key rates were found by examining correlations between rate variations and resulting abundance changes. This method is superior to studying flow patterns, especially when the flows are complex, and to individual, sequential variation of a few rates.

In the ONeMg cores of $8.8-9.5~{\rm M}_\odot$ stars, neon and oxygen burning is ignited off-center. Whether the neon-oxygen flame propagates to the center is critical to determine whether these stars undergo Fe core collapse or electron capture induced ONeMg core collapse. We present more details of stars that ignite neon and oxygen burning off-center. The neon flame is established in a similar manner to the carbon flame of super-AGB stars, albeit with a narrower flame width. The criteria for establishing a flame are able to be met if the strict Schwarzschild criterion for convective instability is adopted. Mixing across the interface of the convective shell disrupts the conditions for the propagation of the burning front and instead the shell burns as a series of inward-moving flashes. While this may not directly affect whether the burning will reach the center (as in super-AGB stars), the core is allowed to contract between each shell flash. Reduction of the electron fraction in the shell reduces the Chandrasekhar mass and the center reaches the threshold density for the URCA process to activate and steer the remaining evolution of the core. This highlights the importance of a more accurate treatment of mixing in the stellar interior for yet another important question in stellar astrophysics - determining the properties of stellar evolution and supernova progenitors at the boundary between electron capture supernova and iron core-collapse supernova.

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.

The $s$-process nucleosynthesis in Asymptotic Giant Branch (AGB) stars depends on the modeling of convective boundaries. We present models and s-process simulations that adopt a treatment of convective boundaries based on the results of hydrodynamic simulations and on the theory of mixing due to gravity waves in the vicinity of convective boundaries. Hydrodynamics simulations suggest the presence of convective boundary mixing (CBM) at the bottom of the thermal pulse-driven convective zone. Similarly, convection-induced mixing processes are proposed for the mixing below the convective envelope during third dredge-up where the 13C pocket for the s process in AGB stars forms. In this work we apply a CBM model motivated by simulations and theory to models with initial mass $M = 2$ and $M = 3M_\odot$, and with initial metal content Z = 0.01 and Z = 0.02. As reported previously, the He-intershell abundance of 12C and 16O are increased by CBM at the bottom of pulse-driven convection zone. This mixing is affecting the $^{22}Ne(\alpha,n)^{25}Mg$ activation and the s-process effciency in the 13C-pocket. In our model CBM at the bottom of the convective envelope during the third dredgeup represents gravity wave mixing. We take further into account that hydrodynamic simulations indicate a declining mixing efficiency already about a pressure scale height from the convective boundaries, compared to mixing-length theory. We obtain the formation of the 13C-pocket with a mass of $\approx 10^{-4}M_\odot$. The final $s$-process abundances are characterized by 0.36 < [s=Fe] < 0.78 and the heavy-to-light s-process ratio is 0.23 < [hs=ls] < 0.45. Finally, we compare our results with stellar observations, pre-solar grain measurements and previous work.

We provide a set of stellar evolution and nucleosynthesis calculations that applies established physics assumptions simultaneously to low- and intermediate-mass and massive star models. Our goal is to provide an internally consistent and comprehensive nuclear production and yield data base for applications in areas such as pre-solar grain studies. Our non-rotating models assume convective boundary mixing where it has been adopted before. We include 8 (12) initial masses for $Z = 0.01$ ($0.02$). Models are followed either until the end of the asymptotic giant branch phase or the end of Si burning, complemented by a simple analytic core-collapse supernova models with two options for fallback and shock velocities. The explosions show which pre-supernova yields will most strongly be effected by the explosive nucleosynthesis. We discuss how these two explosion parameters impacts the light elements and the $s$ and $p$ process. For low- and intermediate-mass models our stellar yields from H to Bi include the effect of convective boundary mixing at the He-intershell boundaries and the stellar evolution feedback of the mixing process that produces the $^{13}$C pocket. All post-processing nucleosynthesis calculations use the same nuclear reaction rate network and nuclear physics input. We provide a discussion of the nuclear production across the entire mass range organized by element group. All our stellar nucleosynthesis profile and time evolution output is available electronically, and tools to explore the data on the NuGrid VOspace hosted by the Canadian Astronomical Data Centre are introduced.

The discovery of the most luminous supernova ASASSN-15lh triggered a shock-wave in the supernova community. The three possible mechanisms proposed for the majority of other superluminous supernovae do not produce a realistic physical model for this particular supernova. In the present study we show the limiting luminosity available from a nickel-powered pair-instability supernova. We computed a few exotic nickel-powered explosions with a total mass of nickel up to 1500 solar masses. We used the hydrostatic configurations prepared with the GENEVA and MESA codes, and the STELLA radiative-transfer code for following the explosion of these models. We show that 1500 solar masses of radioactive nickel is needed to power a luminosity of 2x10^45 erg/s. The resulting light curve is very broad and incompatible with the shorter ASASSN-15lh time-scale. This rules out a nickel-powered origin of ASASSN-15lh. In addition, we derive a simple peak luminosity - nickel mass relation from our data, which may serve to estimate of nickel mass from observed peak luminosities.

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

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

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.

Nucleosynthesis of heavy elements requires the use of different experimental and theoretical methods to determine astrophysical reaction rates than light element nucleosynthesis. Additionally, there are also larger uncertainties involved in the astrophysical models, both because the sites are not well known and because of differing numerical treatments in different models. As an example for the latter, the production of p-nuclei is compared in two different stellar models, demonstrating that a model widely used for postproduction calculations may have a zone grid too coarse to follow the synthesis of p-nuclei in detail.

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.

Carbon-rich grains with isotopic anomalies compared to the Sun are found in primitive meteorites. They were made by stars, and carry the original stellar nucleosynthesis signature. Silicon carbide grains of Type X and C, and low-density graphites condensed in the ejecta of core-collapse supernovae. We present a new set of models for the explosive He shell and compare them with the grains showing 12C/13C and 14N/15N ratios lower than solar. In the stellar progenitor H was ingested into the He shell and not fully destroyed before the explosion. Different explosion energies and H concentrations are considered. If the SN shock hits the He-shell region with some H still present, the models can reproduce the C and N isotopic signatures in C-rich grains. Hot-CNO cycle isotopic signatures are obtained, including a large production of 13C and 15N. The short-lived radionuclides 22Na and 26Al are increased by orders of magnitude. The production of radiogenic 22Ne from the decay of 22Na in the He shell might solve the puzzle of the Ne-E(L) component in low-density graphite grains. This scenario is attractive for the SiC grains of type AB with 14N/15N ratios lower than solar, and provides an alternative solution for SiC grains originally classified as nova grains. Finally, this process may contribute to the production of 14N and 15N in the Galaxy, helping to produce the 14N/15N ratio in the solar system.

Aims. We explore the problem of the site of production of Eu. We use also the information present in the observed spread in the Eu abundances in the early Galaxy, not only its average trend. Moreover, we extend to other heavy elements (Ba, Sr, Rb, Zr) our investigations to provide additional constraints to our results. Methods. We adopt a stochastic chemical evolution model taking into account inhomogeneous mixing. The adopted yields of Eu from neutron star mergers (NSM) and from core-collapse supernovae (SNII) are those that are able to explain the average [Eu/Fe]-[Fe/H] trend observed for solar neighborhood stars, in the framework of a well-tested homogeneous model for the chemical evolution of the MilkyWay. Rb, Sr, Zr, and Ba are produced by both the s- and r-process. The s-process contribution by spinstars is the same as in our previous papers. Results. NSM that merge in less than 10 Myr or NSM combined with a source of r-process generated by massive stars can explain the spread of [Eu/Fe] in the Galactic halo. The combination of r-process production by NSM and s-process production by spinstars is able to reproduce the available observational data for Sr, Zr and Ba. We also show the first predictions for Rb in the Galactic halo. Conclusions. We confirm previous results that either NSM with very short time scale or both NSM and at least a fraction of SNII should have contributed to the synthesis of Eu in the Galaxy. The r-process production by NSM - complemented by an s-process production by spinstars - provide results compatible with our previous findings based on other r-process sites. We critically discuss the weak and strong points of both NSM and SNII scenarios for producing Eu and eventually suggest that the best solution is probably a mixed one in which both sources produce Eu. In fact, this scenario better reproduces the scatter observed in all the studied elements. [abridged]

The knowledge of the production of extinct radioactivities like 92Nb and 146Sm by photodisintegration processes in ccSN and SNIa models is essential for interpreting abundances in meteoritic material and for Galactic Chemical Evolution (GCE). The 92Mo/92Nb and 146Sm/144Sm ratios provide constraints for GCE and production sites. We present results for SNIa with emphasis on nuclear uncertainties.

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

The discovery of 150 - 300 M$_{\odot}$ stars in the Local Group and pair-instability supernova candidates at low redshifts has excited interest in this exotic explosion mechanism. Realistic light curves for pair-instability supernovae at near-solar metallicities are key to identifying and properly interpreting these events as more are found. We have modeled pair-instability supernovae of 150 - 500 M$_{\odot}$ Z $\sim$ 0.1 - 0.4 Z$_{\odot}$ stars. These stars lose up to 80% of their mass to strong line-driven winds and explode as bare He cores. We find that their light curves and spectra are quite different from those of Population III pair-instability explosions, which therefore cannot be used as templates for low-redshift events. Although non-zero metallicity pair-instability supernovae are generally dimmer than their Population III counterparts, in some cases they will be bright enough to be detected at the earliest epochs at which they can occur, the formation of the first galaxies at $z \sim$ 10 - 15. Others can masquerade as dim, short duration supernovae that are only visible in the local universe and that under the right conditions could be hidden in a wide variety of supernova classes. We also report for the first time that some pair-instability explosions can create black holes with masses of $\sim$ 100 M$_{\odot}$.

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

In this chapter, after a brief introduction and overview of stellar evolution, we discuss the evolution and nucleosynthesis of very massive stars (VMS: M>100 solar masses) in the context of recent stellar evolution model calculations. This chapter covers the following aspects: general properties, evolution of surface properties, late central evolution, and nucleosynthesis including their dependence on metallicity, mass loss and rotation. 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 VMS at metallicities close to solar end their life as WC(-WO) type Wolf-Rayet stars. Due to very important mass loss through stellar winds, these stars may have luminosities during the advanced phases of their evolution similar to stars with initial masses between 60 and 120 solar masses. A distinctive feature which may be used to disentangle Wolf-Rayet stars originating from VMS from those originating from lower initial masses is the enhanced abundances of neon and magnesium at the surface of WC stars. At solar metallicity, mass loss is so strong that even if a star is born with several hundred solar masses, it will end its life with less than 50 solar masses (using current mass loss prescriptions). At the metallicity of the LMC and lower, on the other hand, mass loss is weaker and might enable star to undergo pair-instability supernovae.

We investigate the enrichment in elements produced by the slow neutron-capture process ($s$-process) in the globular clusters M4 (NGC 6121) and M22 (NGC 6656). Stars in M4 have homogeneous abundances of Fe and neutron-capture elements, but the entire cluster is enhanced in $s$-process elements (Sr, Y, Ba, Pb) relative to other clusters with a similar metallicity. In M22, two stellar groups exhibit different abundances of Fe and $s$-process elements. By subtracting the mean abundances of $s$-poor from $s$-rich stars, we derive $s$-process residuals or empirical $s$-process distributions for M4 and M22. We find that the $s$-process distribution in M22 is more weighted toward the heavy $s$-peak (Ba, La, Ce) and Pb than M4, which has been enriched mostly with light $s$-peak elements (Sr, Y, Zr). We construct simple chemical evolution models using yields from massive star models that include rotation, which dramatically increases $s$-process production at low metallicity. We show that our massive star models with rotation rates of up to 50\% of the critical (break-up) velocity and changes to the preferred $^{17}$O($\alpha$,$\gamma$)$^{21}$Ne rate produce insufficient heavy $s$-elements and Pb to match the empirical distributions. For models that incorporate AGB yields, we find that intermediate-mass yields (with a $^{22}$Ne neutron source) alone do not reproduce the light-to-heavy $s$-element ratios for M4 and M22, and that a small contribution from models with a $^{13}$C pocket is required. With our assumption that $^{13}$C pockets form for initial masses below a transition range between 3.0 and 3.5 M$_\odot$, we match the light-to-heavy s-element ratio in the s-process residual of M22 and predict a minimum enrichment timescale of between 240 and 360 Myr. Our predicted value is consistent with the 300 Myr upper limit age difference between the two groups derived from isochrone fitting.

The [Sr/Ba] and [Y/Ba] scatter observed in some galactic halo stars that are very metal-poor stars and in a few individual stars of the oldest known Milky Way globular cluster NGC 6522,have been interpreted as evidence of early enrichment by massive fast-rotating stars (spinstars). Because NGC 6522 is a bulge globular cluster, the suggestion was that not only the very-metal poor halo stars, but also bulge stars at [Fe/H]~-1 could be used as probes of the stellar nucleosynthesis signatures from the earlier generations of massive stars, but at much higher metallicity. For the bulge the suggestions were based on early spectra available for stars in NGC 6522, with a medium resolution of R~22,000 and a moderate signal-to-noise ratio. The main purpose of this study is to re-analyse the NGC 6522 stars previously reported using new high-resolution (R~45,000) and high signal-to-noise spectra (S/N>100). We aim at re-deriving their stellar parameters and elemental ratios, in particular the abundances of the neutron-capture s-process-dominated elements such as Sr, Y, Zr, La, and Ba, and of the r-element Eu. High-resolution spectra of four giants belonging to the bulge globular cluster NGC 6522 were obtained at the 8m VLT UT2-Kueyen telescope with the UVES spectrograph in FLAMES-UVESconfiguration. The spectroscopic parameters were derived based on the excitation and ionization equilibrium of \ion{Fe}{I} and \ion{Fe}{II}. Our analysis confirms a metallicity [Fe/H] = -0.95+-0.15 for NGC 6522, and the overabundance of the studied stars in Eu (with +~0.2 < [Eu/Fe] < +~0.4) and alpha-elements O and Mg. The neutron-capture s-element-dominated Sr, Y, Zr, Ba, La now show less pronounced variations from star to star. Enhancements are in the range 0.0 < [Sr/Fe] < +0.4, +0.23 < [Y/Fe] < +0.43, 0.0 < [Zr/Fe] < +0.4, 0.0 < [La/Fe] < +0.35,and 0.05 < [Ba/Fe] < +0.55.

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

Carbon-rich grains are observed to condense in the ejecta of recent core-collapse supernovae, within a year after the explosion. Silicon carbide grains of type X are C-rich grains with isotpic signatures of explosive supernova nucleosynthesis have been found in primitive meteorites. Much rarer silicon carbide grains of type C are a special sub-group of SiC grains from supernovae. They show peculiar abundance signatures for Si and S, isotopically heavy Si and isotopically light S, which appear to to be in disagreement with model predictions. We propose that C grains are formed mostly from C-rich stellar material exposed to lower SN shock temperatures than the more common type X grains. In this scenario, extreme $^{32}$S enrichments observed in C grains may be explained by the presence of short-lived $^{32}$Si ($\tau$$_{1/2}$ = 153 years) in the ejecta, produced by neutron capture processes starting from the stable Si isotopes. No mixing from deeper Si-rich material and/or fractionation of Si from S due to molecular chemistry is needed to explain the $^{32}$S enrichments. The abundance of $^{32}$Si in the grains can provide constraints on the neutron density reached during the supernova explosion in the C-rich He shell material. The impact of the large uncertainty of the neutron capture cross sections in the $^{32}$Si region is discussed.

The stellar mass range 8<M/Mo<12 corresponds to the most massive AGB stars and the most numerous massive stars. It is host to a variety of supernova progenitors and is therefore very important for galactic chemical evolution and stellar population studies. In this paper, we study the transition from super-AGB star to massive star and find that a propagating neon-oxygen burning shell is common to both the most massive electron capture supernova (EC-SN) progenitors and the lowest mass iron-core collapse supernova (FeCCSN) progenitors. Of the models that ignite neon burning off-center, the 9.5Mo model would evolve to an FeCCSN after the neon-burning shell propagates to the center, as in previous studies. The neon-burning shell in the 8.8Mo model, however, fails to reach the center as the URCA process and an extended (0.6 Mo) region of low Ye (0.48) in the outer part of the core begin to dominate the late evolution; the model evolves to an EC-SN. This is the first study to follow the most massive EC-SN progenitors to collapse, representing an evolutionary path to EC-SN in addition to that from SAGB stars undergoing thermal pulses. We also present models of an 8.75Mo super-AGB star through its entire thermal pulse phase until electron captures on 20Ne begin at its center and of a 12Mo star up to the iron core collapse. We discuss key uncertainties and how the different pathways to collapse affect the pre-supernova structure. Finally, we compare our results to the observed neutron star mass distribution.

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.

About a year after core collapse supernova, dust starts to condense in the ejecta. In meteorites, a fraction of C-rich presolar grains (e.g., silicon carbide (SiC) grains of Type-X and low density graphites) are identified as relics of these events, according to the anomalous isotopic abundances. Several features of these abundances remain unexplained and challenge the understanding of core-collapse supernovae explosions and nucleosynthesis. We show, for the first time, that most of the measured C-rich grain abundances can be accounted for in the C-rich material from explosive He burning in core-collapse supernovae with high shock velocities and consequent high temperatures. The inefficiency of the $^{12}$C($\alpha$,$\gamma$)$^{16}$O reaction relative to the rest of the $\alpha$-capture chain at $T > 3.5\times10^8 \mathrm{K}$ causes the deepest He-shell material to be carbon rich and silicon rich, and depleted in oxygen. The isotopic ratio predictions in part of this material, defined here as the C/Si zone, are in agreement with the grain data. The high-temperature explosive conditions that our models reach at the bottom of the He shell, can also be representative of the nucleosynthesis in hypernovae or in the high-temperature tail of a distribution of conditions in asymmetric supernovae. Finally, our predictions are consistent with the observation of large $^{44}$Ca/$^{40}$Ca observed in the grains. This is due to the production of $^{44}$Ti together with $^{40}$Ca in the C/Si zone, and/or to the strong depletion of $^{40}$Ca by neutron captures.

Despite much effort in the past decades, the C-burning reaction rate is uncertain by several orders of magnitude, and the relative strength between the different channels 12C(12C,alpha)20Ne, 12C(12C,p)23Na and 12C(12C,n)23Mg is poorly determined. Additionally, in C-burning conditions a high 12C+12C rate may lead to lower central C-burning temperatures and to 13C(alpha,n)16O emerging as a more dominant neutron source than 22Ne(alpha,n)25Mg, increasing significantly the s-process production. This is due to the rapid decrease of the 13N(gamma,p)12C with decreasing temperature, causing the 13C production via 13N(beta+)13C. Presented here is the impact of the 12C+12C reaction uncertainties on the s-process and on explosive p-process nucleosynthesis in massive stars, including also fast rotating massive stars at low metallicity. Using various 12C+12C rates, in particular an upper and lower rate limit of ~ 50000 higher and ~ 20 lower than the standard rate at 5*10^8 K, five 25 Msun stellar models are calculated. The enhanced s-process signature due to 13C(alpha,n)16O activation is considered, taking into account the impact of the uncertainty of all three C-burning reaction branches. Consequently, we show that the p-process abundances have an average production factor increased up to about a factor of 8 compared to the standard case, efficiently producing the elusive Mo and Ru proton-rich isotopes. We also show that an s-process being driven by 13C(alpha,n)16O is a secondary process, even though the abundance of 13C does not depend on the initial metal content. Finally, implications for the Sr-peak elements inventory in the Solar System and at low metallicity are discussed.

We use contemporary evolutionary models for Very Massive Stars (VMS) to assess whether the Eddington limit constrains the upper stellar mass limit. We also consider the interplay between mass and age for the wind properties and spectral morphology of VMS, with reference to the recently modified classification scheme for O2-3.5If*/WN stars. Finally, the death of VMS in the local universe is considered in the context of pair instability supernovae.

(shortened) The first couple of stellar generations may have been massive, of order 100 Msun, and to have played a dominant role in galaxy formation and the chemical enrichment of the early Universe. Some fraction of these objects may have died as pair-instability supernovae or gamma-ray bursts. The winds if these stars may have played an important role in determining these outcomes. As the winds are driven by radiation pressure on spectral lines, their strengths are expected to vary with metallicity. Until now, most mass-loss predictions for metal-poor O-type stars have assumed a scaled-down solar-abundance pattern. However, Population III evolutionary tracks show significant surface enrichment through rotational mixing of CNO-processed material, because even metal-poor stars switch to CNO-burning early on. We address the question of whether the CNO surface enhanced self-enrichment in the first few generations of stars could impact their mass-loss properties. For this, we employ Monte Carlo simulations to establish the local line-force and solve for the momentum equation of the stellar outflow, testing whether an outflow can actually be established by assessing the net acceleration at the sonic point of the flow. Stellar evolution models of rotating metal-poor stars are used to specify the surface chemical composition, focussing on the phases of early enrichment. We find that the mass-loss rates of CNO enhanced metal-poor stars are higher than those of non-enriched stars, but they are much lower than those rates where the CNO abundance is included in the total abundance Z. We present a heuristic formula that provides mass-loss estimates for CNO-dominated winds in relation to scaled-down solar abundances.

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.

[Shortened] The 12C + 12C fusion reaction has been the subject of considerable experimental efforts to constrain uncertainties at temperatures relevant for stellar nucleosynthesis. In order to investigate the effect of an enhanced carbon burning rate on massive star structure and nucleosynthesis, new stellar evolution models and their yields are presented exploring the impact of three different 12C + 12C reaction rates. Non-rotating stellar models were generated using the Geneva Stellar Evolution Code and were later post-processed with the NuGrid Multi-zone Post-Processing Network tool. The enhanced rate causes core carbon burning to be ignited more promptly and at lower temperature. This reduces the neutrino losses, which increases the core carbon burning lifetime. An increased carbon burning rate also increases the upper initial mass limit for which a star exhibits a convective carbon core. Carbon shell burning is also affected, with fewer convective-shell episodes and convection zones that tend to be larger in mass. Consequently, the chance of an overlap between the ashes of carbon core burning and the following carbon shell convection zones is increased, which can cause a portion of the ashes of carbon core burning to be included in the carbon shell. Therefore, during the supernova explosion, the ejecta will be enriched by s-process nuclides synthesized from the carbon core s process. The yields were used to estimate the weak s-process component in order to compare with the solar system abundance distribution. The enhanced rate models were found to produce a significant proportion of Kr, Sr, Y, Zr, Mo, Ru, Pd and Cd in the weak component, which is primarily the signature of the carbon-core s process. Consequently, it is shown that the production of isotopes in the Kr-Sr region can be used to constrain the 12C + 12C rate using the current branching ratio for a- and p-exit channels.

Context. Rotation is known to affect the nucleosynthesis of light elements in massive stars, mainly by rotation-induced mixing. In particular, rotation boosts the primary nitrogen production. Models of rotating stars are able to reproduce the nitrogen observed in low-Z halo stars. Aims. Here we present the first grid of stellar models for rotating massive stars at low Z, where a full s-process network is used to study the impact of rotation-induced mixing on the nucleosynthesis of heavy elements. Methods. We used the Geneva stellar evolution code that includes an enlarged reaction network with nuclear species up to bismuth to calculate 25 M$_\odot$ models at three different Z and with different initial rotation rates. Results. First, we confirm that rotation-induced mixing leads to a production of primary $^{22}$Ne, which is the main neutron source for the s process in massive stars. Therefore rotation boosts the s process in massive stars at all Z. Second, the neutron-to-seed ratio increases with decreasing Z in models including rotation, which leads to the complete consumption of all iron seeds at Z < 1e-3 by the end of core He-burning. Thus at low Z, the iron seeds are the main limitation for this boosted s process. Third, as Z decreases, the production of elements up to the Ba peak increases at the expense of the elements of the Sr peak. We studied the impact of the initial rotation rate and of the uncertain $^{17}$O$(\alpha,\gamma)$ rate (which strongly affects the neutron poison strength of $^{16}$O) on our results. This study shows that rotating models can produce significant amounts of elements up to Ba over a wide range of Z. Fourth, compared to the He-core, the primary $^{22}$Ne production in the He-shell is even higher (> 1% in mass fraction at all Z), which could open the door for an explosive neutron capture nucleosynthesis in the He-shell, with a primary neutron source.

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

We investigate the quantitative constraint on the triple-alpha reaction rate based on stellar evolution theory, motivated by the recent significant revision of the rate proposed by nuclear physics calculations. Targeted stellar models were computed in order to investigate the impact of that rate in the mass range of 0.8 < M / Msun < 25 and in the metallicity range between Z = 0 and Z = 0.02. The revised rate has a significant impact on the evolution of low- and intermediate-mass stars, while its influence on the evolution of massive stars (M >~ 10 Msun) is minimal. We find that employing the revised rate suppresses helium shell flashes on AGB phase for stars in the initial mass range 0.8 < M / Msun < 6, which is contradictory to what is observed. The absence of helium shell flashes is due to the weak temperature dependence of the revised triple-alpha reaction cross section at the temperature involved. In our models, it is suggested that the temperature dependence of the cross section should have at least nu > 10 at T = 1 - 1.2 x 10^8 K where the cross section is proportional to T^{nu}. We also derive the helium ignition curve to estimate the maximum cross section to retain the low-mass first red giants. The semi-analytically derived ignition curves suggest that the reaction rate should be less than ~ 10^{-29} cm^6 s^{-1} mole^{-2} at ~ 10^{7.8} K, which corresponds to about three orders of magnitude larger than that of the NACRE compilation. In an effort to compromise with the revised rates, we calculate and analyze models with enhanced CNO cycle reaction rates to increase the maximum luminosity of the first giant branch. However, it is impossible to reach the typical RGB tip luminosity even if all the reaction rates related to CNO cycles are enhanced by more than ten orders of magnitude.

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)

Depending on mass and metallicity as well as evolutionary phase, stars occasionally experience convective-reactive nucleosynthesis episodes. We specifically investigate the situation when nucleosynthetically unprocessed, H-rich material is convectively mixed with a He-burning zone, for example in convectively unstable shell on top of electron-degenerate cores in AGB stars, young white dwarfs or X-ray bursting neutron stars. Such episodes are frequently encountered in stellar evolution models of stars of extremely low or zero metal content [...] We focus on the convective-reactive episode in the very-late thermal pulse star Sakurai's object (V4334 Sagittarii). Asplund etal. (1999) determined the abundances of 28 elements, many of which are highly non-solar, ranging from H, He and Li all the way to Ba and La, plus the C isotopic ratio. Our simulations show that the mixing evolution according to standard, one-dimensional stellar evolution models implies neutron densities in the He that are too low to obtain a significant neutron capture nucleosynthesis on the heavy elements. We have carried out 3D hydrodynamic He-shell flash convection [...] we assume that the ingestion process of H into the He-shell convection zone leads only after some delay time to a sufficient entropy barrier that splits the convection zone [...] we obtain significantly higher neutron densities (~few 10^15 1/cm^3) and reproduce the key observed abundance trends found in Sakurai's object. These include an overproduction of Rb, Sr and Y by about 2 orders of magnitude higher than the overproduction of Ba and La. Such a peculiar nucleosynthesis signature is impossible to obtain with the mixing predictions in our one-dimensional stellar evolution models. [...] We determine how our results depend on uncertainties of nuclear reaction rates, for example for the C13(\alpha, n)O16 reaction.

We recently determined the mass of the most massive star known to the date, R136a1 with a mass at birth 320 times the mass of our sun, as well as the mass of several other stars that are more massive than 150 M. Such massive stars (~150-300 M) may end their life as pair-instability supernovae (PISN) if they retain enough mass until they die. We have calculated a grid of stellar evolution models in order to investigate the impact of mass loss and rotation on the evolution and fate of these very massive stars. As mass loss is very strong at solar metallicity, our models predict that most of the very massive stars will die as type Ic SNe. Only slowly and non-rotating stars at metallicities below that of the LMC might retain enough mass to produce a PISN. This would mean that the first stellar generations might have produced PISN although their chemical signature is not observed in extremely metal poor stars in the halo of our galaxy.

The contribution by massive stars (M > 9 solar masses) to the weak s-process component of the solar system abundances is primarily due to the 22Ne neutron source, which is activated near the end of helium-core burning. The residual 22Ne left over from helium-core burning is then reignited during carbon burning, initiating further s-processing that modifies the isotopic distribution. This modification is sensitive to the stellar structure and the carbon burning reaction rate. Recent work on the 12C + 12C reaction suggests that resonances located within the Gamow peak may exist, causing a strong increase in the astrophysical S-factor and consequently the reaction rate. To investigate the effect of an increased rate, 25 solar mass stellar models with three different carbon burning rates, at solar metallicity, were generated using the Geneva Stellar Evolution Code (GENEC) with nucleosynthesis post-processing calculated using the NuGrid Multi-zone Post-Processing Network code (MPPNP). The strongest rate caused carbon burning to occur in a large convective core rather than a radiative one. The presence of this large convective core leads to an overlap with the subsequent convective carbon-shell, significantly altering the initial composition of the carbon-shell. In addition, an enhanced rate causes carbon-shell burning episodes to ignite earlier in the evolution of the star, igniting the 22Ne source at lower temperatures and reducing the neutron density.

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.

Massive stars and their supernovae are prominent sources of radioactive isotopes, the observations of which thus can help to improve our astrophysical models of those. Our understanding of stellar evolution and the final explosive endpoints such as supernovae or hypernovae or gamma-ray bursts relies on the combination of magneto-hydrodynamics, energy generation due to nuclear reactions accompanying composition changes, radiation transport, and thermodynamic properties (such as the equation of state of stellar matter). Nuclear energy production includes all nuclear reactions triggered during stellar evolution and explosive end stages, also among unstable isotopes produced on the way. Radiation transport covers atomic physics (e.g. opacities) for photon transport, but also nuclear physics and neutrino nucleon/nucleus interactions in late phases and core collapse. Here we want to focus on the astrophysical aspects, i.e. a description of the evolution of massive stars and their endpoints, with a special emphasis on the composition of their ejecta (in form of stellar winds during the evolution or of explosive ejecta). Low and intermediate mass stars end their evolution as a white dwarf with an unburned C and O composition. Massive stars evolve beyond this point and experience all stellar burning stages from H over He, C, Ne, O and Si-burning up to core collapse and explosive endstages. In this chapter we discuss the nucleosynthesis processes involved and the production of radioactive nuclei in more detail.

Spectroscopic analyses of H-rich WN5-6 stars within the young star clusters NGC 3603 and R136 are presented, using archival HST & VLT spectroscopy, & high spatial resolution near-IR photometry. We derive high T* for the WN stars in NGC 3603 (T*~42+/-2 kK) & R136 (T*~53+/-3 kK) plus clumping-corrected dM/dt ~ 2-5x10^-5 Msun/yr which closely agree with theoretical predictions. These stars make a disproportionate contribution to the global budget of their host clusters. R136a1 alone supplies ~7% of N(LyC) of the entire 30 Dor region. Comparisons with stellar models calculated for the main-sequence evolution of 85-500 Msun suggest ages of ~1.5 Myr & M_init in the range 105 - 170 Msun for 3 systems in NGC 3603, plus 165-320 Msun for 4 stars in R136. Our high stellar masses are supported by dynamical mass determinations for the components of NGC 3603 A1. We consider the predicted L_X of the R136 stars if they were close, colliding wind binaries. R136c is consistent with a colliding wind binary system. However, short period, colliding wind systems are excluded for R136a WN stars if mass ratios are of order unity. Widely separated systems would have been expected to harden owing to early dynamical encounters with other massive stars in such a dense environment. From simulated star clusters, whose constituents are randomly sampled from the Kroupa IMF, both clusters are consistent with a tentative upper mass limit of ~300 Msun. The Arches cluster is either too old, exhibits a deficiency of very massive stars, or more likely stellar masses have been underestimated - M_init for the most luminous stars in the Arches cluster approach 200 Msun according to contemporary stellar & photometric results. The potential for stars greatly exceeding 150 Msun within metal-poor galaxies suggests that such pair-instability SNe could occur within the local universe, as has been claimed for SN 2007bi (abridged).

The aim of this paper is to examine the effects of the horizontal turbulence in differentially rotating stars on the GSF instability and apply our results to pre-supernova models. For this purpose we derive the expression for the GSF instability with account of the thermal transport and smoothing of the mu-gradient by the horizontal turbulence. We apply the new expressions in numerical models of a 20 solar mass star. We show that if N^2_{Omega} < 0 the Rayleigh-Taylor instability cannot be killed by the stabilizing thermal and mu-gradients, so that the GSF instability is always there and we derive the corresponding diffusion coefficient. The GSF instability grows towards the very latest stages of stellar evolution. Close to the deep convective zones in pre-supernova stages, the transport coefficient of elements and angular momentum by the GSF instability can very locally be larger than the shear instability and even as large as the thermal diffusivity. However the zones over which the GSF instability is acting are extremely narrow and there is not enough time left before the supernova explosion for a significant mixing to occur. Thus, even when the inhibiting effects of the mu-gradient are reduced by the horizontal turbulence, the GSF instability remains insignificant for the evolution. We conclude that the GSF instability in pre-supernova stages cannot be held responsible for the relatively low rotation rate of pulsars compared to the predictions of rotating star models.

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

The slow neutron capture process in massive stars (the weak s-process) produces most of the s-only isotopes in the mass region 60 < A < 90. The nuclear reaction rates used in simulations of this process have a profound effect on the final s-process yields. We generated 1D stellar models of a 25 solar mass star varying the 12C + 12C rate by a factor of 10 and calculated full nucleosynthesis using the post-processing code PPN. Increasing or decreasing the rate by a factor of 10 affects the convective history and nucleosynthesis, and consequently the final yields.

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

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

Many nucleosynthesis and mixing processes of low-mass stars as they evolve from the Main Sequence to the thermal-pulse Asymptotic Giant Branch phase (TP-AGB) are well understood (although of course important physics components, e.g. rotation, magnetic fields, gravity wave mixing, remain poorly known). Nevertheless, in the last years presolar grain measurements with high resolution have presented new puzzling problems and strong constraints on nucleosynthesis processes in stars. The goal of the NuGrid collaboration is to present uniform yields for a large range of masses and metallicities, including low$-$mass stars and massive stars and their explosions. Here we present the first calculations of stellar evolution and high-resolution, post-processing simulations of an AGB star with an initial mass of 2 M_sun and solar-like metallicity (Z=0.01), based on the post-processing code PPN. In particular, we analyze the formation and evolution of the radiative 13C-pocket between the 17th TP and the 18th TP. The s-process nucleosynthesis profile of a sample of heavy isotopes is also discussed, before the next convective TP occurrence.

We review some of the uncertainties in calculating nucleosynthetic yields, focusing on the explosion mechanism. Current yield calculations tend to either use a piston, energy injection, or enhancement of neutrino opacities to drive an explosion. We show that the energy injection, or more accurately, an entropy injection mechanism is best-suited to mimic our current understanding of the convection-enhanced supernova engine. The enhanced neutrino-opacity technique is in qualitative disagreement with simulations of core-collapse supernovae and will likely produce errors in the yields. But piston-driven explosions are the most discrepant. Piston-driven explosion severely underestimate the amount of fallback, leading to order-of-magnitude errors in the yields of heavy elements. To obtain yields accurate to the factor of a few level, we must use entropy or energy injection and this has become the NuGrid collaboration approach.

The "collapsar" engine for gamma-ray bursts invokes as its energy source the failure of a normal supernova and the formation of a black hole. Here we present the results of the first three-dimensional simulation of the collapse of a massive star down to a black hole, including the subsequent accretion and explosion. The explosion differs significantly from the axisymmetric scenario obtained in two-dimensional simulations; this has important consequences for the nucleosynthetic yields. We compare the nucleosynthetic yields to those of hypernovae. Calculating yields from three-dimensional explosions requires new strategies in post-process nucleosynthesis; we discuss NuGrid's plan for three-dimensional yields.

We investigate the physical conditions where 44Ti and 56Ni are created in core-collapse supernovae. In this preliminary work we use a series of post-processing network calculations with parametrized expansion profiles that are representative of the wide range of temperatures, densities and electron-to-baryon ratios found in 3D supernova simulations. Critical flows that affect the final yields of 44Ti and 56Ni are assessed.

Simulations of nucleosynthesis in astrophysical environments are at the intersection of nuclear physics reaction rate research and astrophysical applications, for example in the area of galactic chemical evolution or near-field cosmology. Unfortunately, at present the available yields for such applications are based on heterogeneous assumptions between the various contributing nuclear production sites, both in terms of modeling the thermodynamic environment itself as well as the choice of specifc nuclear reaction rates and compilations. On the other side, new nuclear reaction rate determinations are often taking a long time to be included in astrophysical applications. The NuGrid project addresses these issues by providing a set of codes and a framework in which these codes interact. In this contribution we describe the motivation, goals and first results of the NuGrid project. At the core is a new and evolving post-processing nuclesoynthesis code (PPN) that can follow quiescent and explosive nucleosynthesis following multi-zone 1D-stellar evolution as well as multi-zone hydrodynamic input, including explosions. First results are available in the areas of AGB and massive stars.

The s-process production in massive stars at very low metallicities is expected to be negligible due to the low abundance of the neutron source 22Ne, to primary neutron poisons and decreasing iron seed abundances. However, recent models of massive stars including the effects of rotation show that a strong production of 22Ne is possible in the helium core, as a consequence of the primary nitrogen production (observed in halo metal poor stars). Using the PPN post-processing code, we studied the impact of this primary 22Ne on the s process. We find a large production of s elements between strontium and barium, starting with the amount of primary 22Ne predicted by stellar models. There are several key reaction rate uncertainties influencing the s-process efficiency. Among them, 17O(alpha,gamma) may play a crucial role strongly influencing the s process efficiency, or it may play a negligible role, according to the rate used in the calculations. We also report on the development of a new parallel (MPI) post-processing code (MPPNP) designed to follow the complete nucleosynthesis in stars on highly resolved grids. We present here the first post-processing run from the ZAMS up to the end of helium burning for a 15 solar mass model.

We examine observed heavy element abundances in the Cassiopeia A supernova remnant as a constraint on the nature of the Cas A supernova. We compare bulk abundances from 1D and 3D explosion models and spatial distribution of elements in 3D models with those derived from X-ray observations. We also examine the cospatial production of 26Al with other species. We find that the most reliable indicator of the presence of 26Al in unmixed ejecta is a very low S/Si ratio (~0.05). Production of N in O/S/Si-rich regions is also indicative. The biologically important element P is produced at its highest abundance in the same regions. Proxies should be detectable in supernova ejecta with high spatial resolution multiwavelength observations.

We present preliminary results from recent high-resolution double-degenerate merger simulations with the Smooth Particle Hydrodynamics (SPH) technique. We put particular emphasis on verification and validation in our effort and show the importance of details in the initial condition setup for the final outcome of the simulation. We also stress the dynamical importance of including shocks in the simulations. These results represent a first step toward a suite of simulations that will shed light on the question whether double-degenerate mergers are a viable path toward type 1a supernovae. In future simulations, we will make use of the capabilities of the NuGrid collaboration in post-processing SPH particle trajectories with a complete nuclear network to follow the detailed nuclear reactions during the dynamic merger phase.

The nucleosynthetic yield from a supernova explosion depends upon a variety of effects: progenitor evolution, explosion process, details of the nuclear network, and nuclear rates. Especially in studies of integrated stellar yields, simplifications reduce these uncertainties. But nature is much more complex, and to actually study nuclear rates, we will have to understand the full, complex set of processes involved in nucleosynthesis. Here we discuss a few of these complexities and detail how the NuGrid collaboration will address them.

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.

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.

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.

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

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

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.

Mass loss is a very important aspect of the life of massive stars. After briefly reviewing its importance, we discuss the impact of the recently proposed downward revision of mass loss rates due to clumping (difficulty to form Wolf-Rayet stars and production of critically rotating stars). Although a small reduction might be allowed, large reduction factors around ten are disfavoured. We then discuss the possibility of significant mass loss at very low metallicity due to stars reaching break-up velocities and especially due to the metal enrichment of the surface of the star via rotational and convective mixing. This significant mass loss may help the first very massive stars avoid the fate of pair-creation supernova, the chemical signature of which is not observed in extremely metal poor stars. The chemical composition of the very low metallicity winds is very similar to that of the most metal poor star known to date, HE1327-2326 and offer an interesting explanation for the origin of the metals in this star. We also discuss the importance of mass loss in the context of long and soft gamma-ray bursts and pair-creation supernovae. Finally, we would like to stress that mass loss in cooler parts of the HR-diagram (luminous blue variable and yellow and red supergiant stages) are much more uncertain than in the hot part. More work needs to be done in these areas to better constrain the evolution of the most massive stars.

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

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.

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.

A series of fast rotating models at very low metallicity (Z=1e-8) was computed in order to explain the surface abundances observed at the surface of CEMP stars, in particular for nitrogen. The main results are the following: - Strong mixing occurs during He--burning and leads to important primary nitrogen production. - Important mass loss takes place in the RSG stage for the most massive models. The 85 solar mass model loses about three quarter of its initial mass, becomes a WO star and could produce a GRB. - The CNO elements of HE1327-2326 could have been produced in massive rotating stars and ejected by their stellar winds.

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 .

Two series of models were computed. The first series consists of 20 solar mass models with varying initial metallicity (Z=0.02 down to Z=10^{-8}) and rotation (V_{ini}=0-600 km/s). The second one consists of models with an initial metallicity of Z=10^{-8}, masses between 9 and 85 solar masses and fast initial rotation velocities (V_{ini}=600-800 km/s). The most interesting models are the models with Z=10^{-8} ([Fe/H]~-6.6). In the course of helium burning, carbon and oxygen are mixed into the hydrogen burning shell. This boosts the importance of the shell and causes a reduction of the CO core mass. Later in the evolution, the hydrogen shell deepens and produces large amount of primary nitrogen. For the most massive models (M>~60 solar masses), significant mass loss occurs during the red supergiant stage. This mass loss is due to the surface enrichment in CNO elements via rotational and convective mixing. The 85 solar mass model ends up as a WO type Wolf-Rayet star. Therefore the models predict SNe of type Ic and possibly long and soft GRBs at very low metallicities. The rotating 20 solar mass models can best reproduce the observed CNO abundances at the surface of extremely metal poor (EMP) stars and the metallicity trends when their angular momentum content is the same as at solar metallicity (and therefore have an increasing surface velocity with decreasing metallicity). The wind of the massive star models can also reproduce the CNO abundances of the most metal-poor carbon-rich star known to date, HE1327-2326.

The aim of this paper is to look at the evolution of massive stars in order to determine whether or not the progenitor of V838 Mon may be a massive star. In the first part of this paper, the evolution of massive stars around solar metallicity is described, especially the evolution in the Hertzsprung-Russell (HR) diagram. Then, using the observational constraints, the probable progenitors (and their evolution) are described. Using models of single stars, no progenitor can be found amongst massive stars that can satisfy all the observational constraints. Wolf-Rayet stars (stars with initial masses above about 30 solar masses, which have lost their hydrogen rich envelopes) could explain 10 to 100 solar masses of circumstellar material but they are very luminous (L > 100,000 solar luminosities). Main sequence stars crossing the HR diagram and becoming red supergiants (RSG) can have very low effective temperatures but take thousands of years to cross over. Be stars (fast rotating stars with a mass around 10 solar masses), which form disk or B stars accreting matter from a binary companion of a similar mass would need to be compared in detail with the observational constraints. In the future, there will hopefully be further observational constraints on the models coming from the mass and nature (interstellar or circumstellar) of the material producing the light echo and from a frequency estimate of spectacular objects such as V838 Mon.

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

Nine 20 Mo models were computed with metallicities ranging from solar, through $Z=10^{-5}$ ([Fe/H]=~-3.1) down to $Z=10^{-8}$ ([Fe/H]=~-6.1) and with initial rotational velocities between 0 and 600 km/s to study the impact of initial metallicity and rotational velocity. The very large amounts of N14 observed (~0.03 Mo) are only produced at $Z=10^{-8}$ (PopII 1/2). The strong dependence of the N14 yields on rotation and other parameters like the initial mass and metallicity may explain the large scatter in the observations of N14 abundance. The metallicity trends are best reproduced by the models with Omega_ini/Omega_c=~0.75, which is slightly above the mean observed value for OB solar metallicity stars. Indeed, in the model with Vini=600 km/s at $Z=10^{-8}$, the O16 yield is reduced due to strong mixing. This allows in particular to reproduce the upturn for C/O and a slightly decreasing [C/Fe], which are observed below [Fe/H]=~-3.

Two series of models and their yields are presented in this paper. The first series consists of 20 Mo models with varying initial metallicity (solar down to $Z=10^{-8}$) and rotation (Vini=0-600 km/s). The second one consists of models with an initial metallicity of $Z=10^{-8}$, masses between 20 and 85 Mo and average rotation velocities at these metallicities (Vini=600-800 km/s). The most interesting models are the models with $Z=10^{-8}$ ([Fe/H]=~-6.6). In the course of helium burning, carbon and oxygen are mixed into the hydrogen burning shell. This boosts the importance of the shell and causes a reduction of the size of the CO core. Later in the evolution, the hydrogen shell deepens and produces large amount of primary nitrogen. For the most massive models (M>~60 Mo), significant mass loss occurs during the red supergiant stage. This mass loss is due to the surface enrichment in CNO elements via rotational and convective mixing. The yields of the fast rotating 20 Mo models can best reproduce (within our study) the observed abundances at the surface of extremely metal poor (EMP) stars. The wind of the massive models can reproduce the CNO abundances of the carbon--rich UMPs, in particular for the most metal poor star known to date, HE1327-2326.

Recently, measurements of abundances in extremely metal poor (EMP) stars have brought new constraints on stellar evolution models. In an attempt to explain the origin of the abundances observed, we computed pre--supernova evolution models, explosion models and the related nucleosynthesis. In this paper, we start by presenting the pre-SN models of rotating single stars with metallicities ranging from solar metallicity down to almost metal free. We then review key processes in core-collapse and bounce, before we integrate them in a simplistic parameterization for 3D MHD models, which are well underway and allow one to follow the evolution of the magnetic fields during collapse and bounce. Finally, we present explosive nucleosynthesis results including neutrino interactions with matter, which are calculated using the outputs of the explosion models. The main results of the pre-SN models are the following. First, primary nitrogen is produced in large amount in models with an initial metallicity $Z=10^{-8}$. Second, at the same metallicity of $Z=10^{-8}$ and for models with an initial mass larger than about 60 Mo, rotating models may experience heavy mass loss (up to more than half of the initial mass of the star). The chemical composition of these winds can qualitatively reproduce the abundance patterns observed at the surface of carbon-rich EMP stars. Explosive nucleosynthesis including neutrino-matter interactions produce improved abundances for iron group elements, in particular for scandium and zinc. It also opens the way to a new neutrino and proton rich process ($\nu$p-process) able to contribute to the nucleosynthesis of elements with A > 64. (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.

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

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.

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.

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.

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.

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.