Researcher profile

Stephane Goriely

Stephane Goriely contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate

Trust snapshot

Quick read

Trust 21 - Emerging
19works
0followers
7topics
4close collaborators

Actions

Decide how to stay connected

Follow researcher0

Research graph

See the researcher in context

Open full explorer

Inspect adjacent work, topics, institutions and collaborators without jumping out to a separate graph page.

Building this graph slice

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

19 published item(s)

preprint2022arXiv

Dynamical ejecta of neutron star mergers with nucleonic weak processes II: Kilonova emission

The majority of existing results for the kilonova (or macronova) emission from material ejected during a neutron-star (NS) merger is based on (quasi-)one-zone models or manually constructed toy-model ejecta configurations. In this study we present a kilonova analysis of the material ejected during the first ~10ms of a NS merger, called dynamical ejecta, using directly the outflow trajectories from general relativistic smoothed-particle hydrodynamics simulations including a sophisticated neutrino treatment and the corresponding nucleosynthesis results, which have been presented in Part I of this study. We employ a multi-dimensional two-moment radiation transport scheme with approximate M1 closure to evolve the photon field and use a heuristic prescription for the opacities found by calibration with atomic-physics based reference results. We find that the photosphere is generically ellipsoidal but augmented with small-scale structure and produces emission that is about 1.5-3 times stronger towards the pole than the equator. The kilonova typically peaks after 0.7-1.5days in the near-infrared frequency regime with luminosities between 3-7x10^40erg/s and at photospheric temperatures of 2.2-2.8x10^3K. A softer equation of state or higher binary-mass asymmetry leads to a longer and brighter signal. Significant variations of the light curve are also obtained for models with artificially modified electron fractions, emphasizing the importance of a reliable neutrino-transport modeling. None of the models investigated here, which only consider dynamical ejecta, produces a transient as bright as AT2017gfo. The near-infrared peak of our models is incompatible with the early blue component of AT2017gfo.

preprint2021arXiv

Skyrme-Hartree-Fock-Bogoliubov mass models on a 3D Mesh: Effect of triaxial shape

The modeling of nuclear reactions and radioactive decays in astrophysical or earth-based conditions requires detailed knowledge of the masses of essentially all nuclei. Microscopic mass models based on nuclear energy density functionals (EDFs) can be descriptive and used to provide this information. The concept of intrinsic symmetry breaking is central to the predictive power of EDF approaches, yet is generally not exploited to the utmost by mass models because of the computational demands of adjusting up to about two dozen parameters to thousands of nuclear masses. We report on a first step to bridge the gap between what is presently feasible for studies of individual nuclei and large-scale models: we present a new Skyrme-EDF-based model that was adjusted using a three-dimensional coordinate-space representation, for the first time allowing for both axial and triaxial deformations during the adjustment process. To compensate for the substantial increase in computational cost brought by the latter, we have employed a committee of multilayer neural networks to model the objective function in parameter space and guide us towards the overall best fit. The resulting mass model BSkG1 is computed with the EDF model independently of the neural network. It yields a root mean square (rms) deviation on the 2457 known masses of 741 keV and an rms deviation on the 884 measured charge radii of 0.024 fm.

preprint2016arXiv

The temperature and chronology of heavy-element synthesis in low-mass stars

Roughly half of the heavy elements (atomic mass greater than that of iron) are believed to be synthesized in the late evolutionary stages of stars with masses between 0.8 and 8 solar masses. Deep inside the star, nuclei (mainly iron) capture neutrons and progressively build up (through the slow-neutron-capture process, or s-process) heavier elements that are subsequently brought to the stellar surface by convection. Two neutron sources, activated at distinct temperatures, have been proposed: 13C and 22Ne, each releasing one neutron per alpha-particle (4He) captured. To explain the measured stellar abundances, stellar evolution models invoking the 13C neutron source (which operates at temperatures of about one hundred million kelvin) are favoured. Isotopic ratios in primitive meteorites, however, reflecting nucleosynthesis in the previous generations of stars that contributed material to the Solar System, point to higher temperatures (more than three hundred million kelvin), requiring at least a late activation of 22Ne. Here we report a determination of the s-process temperature directly in evolved low-mass giant stars, using zirconium and niobium abundances, independently of stellar evolution models. The derived temperature supports 13C as the s-process neutron source. The radioactive pair 93Zr-93Nb used to estimate the s-process temperature also provides, together with the pair 99Tc-99Ru, chronometric information on the time elapsed since the start of the s-process, which we determine to be one million to three million years.

preprint2015arXiv

Comprehensive nucleosynthesis analysis for ejecta of compact binary mergers

We present the first comprehensive study of r-process element nucleosynthesis in the ejecta of compact binary mergers (CBMs) and their relic black-hole (BH)-torus systems. The evolution of the BH-accretion tori is simulated for seconds with a Newtonian hydrodynamics code including viscosity effects, pseudo-Newtonian gravity for rotating BHs, and an energy-dependent two-moment closure scheme for the transport of electron neutrinos and antineutrinos. The investigated cases are guided by relativistic double neutron star (NS-NS) and NS-BH merger models, producing ~3-6 Msun BHs with rotation parameters of A~0.8 and tori of 0.03-0.3 Msun. Our nucleosynthesis analysis includes the dynamical (prompt) ejecta expelled during the CBM phase and the neutrino and viscously driven outflows of the relic BH-torus systems. While typically ~20-25% of the initial accretion-torus mass are lost by viscously driven outflows, neutrino-powered winds contribute at most another ~1%, but neutrino heating enhances the viscous ejecta significantly. Since BH-torus ejecta possess a wide distribution of electron fractions (0.1-0.6) and entropies, they produce heavy elements from A~80 up to the actinides, with relative contributions of A>130 nuclei being subdominant and sensitively dependent on BH and torus masses and the exact treatment of shear viscosity. The combined ejecta of CBM and BH-torus phases can reproduce the solar abundances amazingly well for A>90. Varying contributions of the torus ejecta might account for observed variations of lighter elements with 40<Z<56 relative to heavier ones, and a considerable reduction of the prompt ejecta compared to the torus ejecta, e.g. in highly asymmetric NS-BH mergers, might explain the composition of heavy-element deficient stars.

preprint2015arXiv

Detailed homogeneous abundance studies of 14 Galactic s-process enriched post-AGB stars: In search of lead (Pb)

This paper is part of a larger project in which we systematically study the chemical abundances of Galactic and extragalactic post-asymptotic giant branch (post-AGB) stars. Lead (Pb) is the final product of the s-process nucleosynthesis and is predicted to have large overabundances with respect to other s-process elements in AGB stars of low metallicities. However, Pb abundance studies of enriched post-AGB stars in the Magellanic Clouds show the Pb overabundance is not observed. We used high-resolution UVES and HERMES spectra for detailed spectral abundance studies of our sample of 14 Galactic post-AGB stars. We do not find any clear evidence of Pb overabundances in our sample. Stars with T(eff) > 7500 K do not provide strong constraints on the Pb abundance. We conclude that the discrepancy between theory and observation increases towards lower metallicities. All objects, except IRAS 17279-1119, confirm the relation between neutron exposure [hs/ls] and third dredge-up efficiency [s/Fe], whereas no relation between metallicity and neutron exposure is detected within the metallicity range of our total sample (-1.4 < [Fe/H] < -0.2). The mild enrichment of IRAS 17279-1119 can probably be attributed to a cut-off of the AGB evolution due to binary interactions. To our knowledge, IRAS 17279-1119 is the first s-process enhanced Galactic post-AGB star known in a binary system and is a possible precursor of the extrinsic Ba dwarf stars. Lead-rich stars have yet to be found in post-AGB stars. extrinsic Ba dwarf stars. Lead-rich stars are yet to be found in post-AGB stars.

preprint2015arXiv

Impact of the phonon coupling on the photon strength function

The pygmy dipole resonances and photon strength functions in stable and unstable Ni and Sn isotopes are calculated within the microscopic self-consistent version of the extended theory of finite fermi systems which includes phonon coupling effects.The Skyrme forces SLy4 is used. A pygmy dipole resonance in $^{72}$Ni is predicted at the mean energy of 12.4 MeV exhausting 25.7\% of the total energy-weighted sum rule. The microscopically obtained photon E1 strength functions are compared with available experimental data and used to calculate nuclear reaction properties. Average radiative widths and radiative neutron capture cross sections have been calculated taking the phonon coupling into account as well as the uncertainties caused by various microscopic level density models. In all three quantities considered, the contribution of phonon coupling turned out to be significant and is found necessary to explain available experimental data.

preprint2015arXiv

Impact of weak interactions of free nucleons on the r-process in dynamical ejecta from neutron-star mergers

We investigate beta-interactions of free nucleons and their impact on the electron fraction (Y_e) and r-process nucleosynthesis in ejecta characteristic of binary neutron star mergers (BNSMs). For that we employ trajectories from a relativistic BNSM model to represent the density-temperature evolutions in our parametric study. In the high-density environment, positron captures decrease the neutron richness at the high temperatures predicted by the hydrodynamic simulation. Circumventing the complexities of modelling three-dimensional neutrino transport, (anti)neutrino captures are parameterized in terms of prescribed neutrino luminosities and mean energies, guided by published results and assumed as constant in time. Depending sensitively on the adopted neutrino-antineutrino luminosity ratio, neutrino processes increase Y_e to values between 0.25 and 0.40, still allowing for a successful r-process compatible with the observed solar abundance distribution and a significant fraction of the ejecta consisting of r-process nuclei. If the electron neutrino luminosities and mean energies are relatively large compared to the antineutrino properties, the mean Y_e might reach values >0.40 so that neutrino captures seriously compromise the success of the r-process. In this case, the r-abundances remain compatible with the solar distribution, but the total amount of ejected r-material is reduced to a few percent, because the production of iron-peak elements is favored. Proper neutrino physics, in particular also neutrino absorption, have to be included in BNSM simulations before final conclusions can be drawn concerning r-processing in this environment and concerning observational consequences like kilonovae, whose peak brightness and color temperature are sensitive to the composition-dependent opacity of the ejecta.

preprint2015arXiv

Nucleosynthesis in dynamical and torus ejecta of compact binary mergers

We present a comprehensive study of r-process element nucleosynthesis in the ejecta of compact binary mergers (CBMs) and their relic black-hole (BH)-torus systems. The evolution of the BH-accretion tori is simulated for seconds with a Newtonian hydrodynamics code including viscosity effects, pseudo-Newtonian gravity for rotating BHs, and an energy-dependent two-moment closure scheme for the transport of electron neutrinos and antineutrinos. The investigated cases are guided by relativistic double neutron star (NS-NS) and NS-BH merger models, producing ~3-6 Msun BHs with rotation parameters of A~0.8 and tori of 0.03-0.3 Msun. Our nucleosynthesis analysis includes the dynamical (prompt) ejecta expelled during the CBM phase and the neutrino and viscously driven outflows of the relic BH-torus systems. While typically ~20-25% of the initial accretion-torus mass are lost by viscously driven outflows, neutrino-powered winds contribute at most another ~1%, but neutrino heating enhances the viscous ejecta significantly. Since BH-torus ejecta possess a wide distribution of electron fractions (0.1-0.6) and entropies, they produce heavy elements from A~80 up to the actinides, with relative contributions of A>130 nuclei being subdominant and sensitively dependent on BH and torus masses and the exact treatment of shear viscosity. The combined ejecta of CBM and BH-torus phases can reproduce the solar abundances amazingly well for A>90. Varying contributions of the torus ejecta might account for observed variations of lighter elements with 40<Z<56 relative to heavier ones, and a considerable reduction of the prompt ejecta compared to the torus ejecta, e.g. in highly asymmetric NS-BH mergers, might explain the composition of heavy-element deficient stars.

preprint2015arXiv

On microscopic theory of radiative nuclear reaction characteristics

A survey of some results in the modern microscopic theory of properties of nuclear reactions with gamma-rays is given. First of all, we discuss the impact of phonon coupling (PC) on the photon strength function (PSF) because it represents the most natural physical source of additional strength found for Sn isotopes in recent experiments that could not be explained within the stan- dard HFB+QRPA approach. The self-consistent version of the Extended Theory of Finite Fermi Systems in the Quasiparticle Time Blocking Approximation, or simply QTBA, is applied. It uses the HFB mean field and includes both the QRPA and PC effects on the basis of the SLy4 Skyrme force. With our microscopic E1 PSFs, the following properties have been calculated for many stable and unstable even-even semi-magic Sn and Ni isotopes as well as for double-magic 132Sn and 208Pb using the reaction codes EMPIRE and TALYS with several nuclear level density (NLD) models: 1) the neutron capture cross sections, 2) the corresponding neutron capture gamma spectra, 3) the av- erage radiative widths of neutron resonances. In all the properties considered, the PC contribution turned out to be significant, as compared with the standard QRPA one, and necessary to explain the available experimental data. The results with the phenomenological so-called generalized super- fluid NLD model turned out to be worse, on the whole, than those obtained with the microscopic HFB+combinatorial NLD model. Finally, we also discuss the modern microscopic NLD models based on the self-consistent HFB method and show their relevance to explain experimental data as compared with the phenomeno- logical models. The use of these self-consistent microscopic approaches is of particular relevance for nuclear astrophysics, but also for the study of double-magic nuclei.

preprint2014arXiv

Neutron-powered precursors of kilonovae

The merger of binary neutron stars (NSs) ejects a small quantity of neutron rich matter, the radioactive decay of which powers a day to week long thermal transient known as a kilonova. Most of the ejecta remains sufficiently dense during its expansion that all neutrons are captured into nuclei during the r-process. However, recent general relativistic merger simulations by Bauswein and collaborators show that a small fraction of the ejected mass (a few per cent, or ~1e-4 Msun) expands sufficiently rapidly for most neutrons to avoid capture. This matter originates from the shocked-heated interface between the merging NSs. Here we show that the beta-decay of these free neutrons in the outermost ejecta powers a `precursor&#39; to the main kilonova emission, which peaks on a timescale of a few hours following merger at U-band magnitude ~22 (for an assumed distance of 200 Mpc). The high luminosity and blue colors of the neutron precursor render it a potentially important counterpart to the gravitational wave source, that may encode valuable information on the properties of the merging binary (e.g. NS-NS versus NS-black hole) and the NS equation of state. Future work is necessary to assess the robustness of the fast moving ejecta and the survival of free neutrons in the face of neutrino absorptions, although the precursor properties are robust to a moderate amount of leptonization. Our results provide additional motivation for short latency gravitational wave triggers and rapid follow-up searches with sensitive ground based telescopes.

preprint2014arXiv

Nucleosynthesis constraints on the neutron star-black hole merger rate

We derive constraints on the time-averaged event rate of neutron star-black hole (NS-BH) mergers by using estimates of the population-integrated production of heavy rapid neutron-capture (r-process) elements with nuclear mass numbers A > 140 by such events in comparison to the Galactic repository of these chemical species. Our estimates are based on relativistic hydrodynamical simulations convolved with theoretical predictions of the binary population. This allows us to determine a strict upper limit of the average NS-BH merger rate of ~6*10^-5 per year. We quantify the uncertainties of this estimate to be within factors of a few mostly because of the unknown BH spin distribution of such systems, the uncertain equation of state of NS matter, and possible errors in the Galactic content of r-process material. Our approach implies a correlation between the merger rates of NS-BH binaries and of double NS systems. Predictions of the detection rate of gravitational-wave signals from such compact-object binaries by Advanced LIGO and Advanced Virgo on the optimistic side are incompatible with the constraints set by our analysis.

preprint2014arXiv

Stellar origin of the 182Hf cosmochronometer and the presolar history of solar system matter

Among the short-lived radioactive nuclei inferred to be present in the early solar system via meteoritic analyses, there are several heavier than iron whose stellar origin has been poorly understood. In particular, the abundances inferred for 182Hf (half-life = 8.9 million years) and 129I (half-life = 15.7 million years) are in disagreement with each other if both nuclei are produced by the rapid neutron-capture process. Here, we demonstrate that contrary to previous assumption, the slow neutron-capture process in asymptotic giant branch stars produces 182Hf. This has allowed us to date the last rapid and slow neutron-capture events that contaminated the solar system material at roughly 100 million years and 30 million years, respectively, before the formation of the Sun.

preprint2013arXiv

NACRE II: an update of the NACRE compilation of charged-particle-induced thermonuclear reaction rates for nuclei with mass number $A < 16$

An update of the NACRE compilation [Angulo et al., Nucl. Phys. A 656 (1999) 3] is presented. This new compilation, referred to as NACRE II, reports thermonuclear reaction rates for 34 charged-particle induced, two-body exoergic reactions on nuclides with mass number $A<16$, of which fifteen are particle-transfer reactions and the rest radiative capture reactions. When compared with NACRE, NACRE II features in particular (1) the addition to the experimental data collected in NACRE of those reported later, preferentially in the major journals of the field by early 2013, and (2) the adoption of potential models as the primary tool for extrapolation to very low energies of astrophysical $S$-factors, with a systematic evaluation of uncertainties. As in NACRE, the rates are presented in tabular form for temperatures in the $10^{6}$ $\simeq\leq$ T $\leq$ $10^{10}$ K range. Along with the &#39;adopted&#39; rates, their low and high limits are provided. The new rates are available in electronic form as part of the Brussels Library (BRUSLIB) of nuclear data. The NACRE II rates also supersede the previous NACRE rates in the Nuclear Network Generator (NETGEN) for astrophysics. [http://www.astro.ulb.ac.be/databases.html.]

preprint2012arXiv

Databases and tools for nuclear astrophysics applications BRUSsels Nuclear LIBrary (BRUSLIB), Nuclear Astrophysics Compilation of REactions II (NACRE II) and Nuclear NETwork GENerator (NETGEN)

An update of a previous description of the BRUSLIB+NACRE package of nuclear data for astrophysics and of the web-based nuclear network generator NETGEN is presented. The new version of BRUSLIB contains the latest predictions of a wide variety of nuclear data based on the most recent version of the Brussels-Montreal Skyrme-HFB model. The nuclear masses, radii, spin/parities, deformations, single-particle schemes, matter densities, nuclear level densities, E1 strength functions, fission properties, and partition functions are provided for all nuclei lying between the proton and neutron drip lines over the 8<=Z<=110 range, whose evaluation is based on a unique microscopic model that ensures a good compromise between accuracy, reliability, and feasibility. In addition, these various ingredients are used to calculate about 100000 Hauser-Feshbach n-, p-, a-, and gamma-induced reaction rates based on the reaction code TALYS. NACRE is superseded by the NACRE II compilation for 15 charged-particle transfer reactions and 19 charged-particle radiative captures on stable targets with mass numbers A < 16. NACRE II features the inclusion of experimental data made available after the publication of NACRE in 1999 and up to 2011. In addition, the extrapolation of the available data to the very low energies of astrophysical relevance is improved through the systematic use of potential models. Uncertainties in the rates are also evaluated on this basis. The latest release of the web-based tool NETGEN is presented. It contains in a fully documented form the new BRUSLIB and NACRE II data, as well as new experiment-based radiative neutron capture cross sections. The full new versions of BRUSLIB, NACRE II, and NETGEN are available electronically at http://www.astro.ulb.ac.be/NuclearData. The nuclear material is presented in an extended tabular form complemented with a variety of graphical interfaces.

preprint2012arXiv

Post-AGB stars in the SMC as tracers of stellar evolution: the extreme s-process enrichment of the 21 micrometer star J004441.04-732136.4

By virtue of their spectral types, favourable bolometric corrections as well as their constrained distances, post-AGB stars in external galaxies offer unprecedented tests to AGB nucleosynthesis and dredge-up predictions. We focus here on one object J004441.04-732136.4, which is the only known 21 micrometer source of the SMC. Spectral abundance results reveal J004441.04-732136.4 to be one of the most s-process enriched objects found up to date, while the photospheric C/O ratio of 1.9 +- 0.7, shows the star is only modestly C-rich. J004441.04-732136.4 also displays a low [Fe/H] = -1.34 +- 0.32, which is significantly lower than the mean metallicity of the SMC. From the SED, a luminosity of 7600 +- 200 solar luminosities is found, together with E(B-V) = 0.64 +- 0.02. According to evolutionary post-AGB tracks, the initial mass should be approximately 1.3 solar masses. The photometric variability shows a clear period of 97.6 +- 0.3 days. The detected C/O as well as the high s-process overabundances (e.g. [Y/Fe] = 2.15, [La/Fe] = 2.84) are hard to reconcile with the predictions. The chemical models also predict a high Pb abundance, which is not compatible with the detected spectrum, and a very high 12C/13C, which is not yet constrained by observations. The predictions are only marginally dependent on the evolution codes used. We show that our theoretical predictions match the s-process distribution, but fail in reproducing the detected high overabundances and predict a high Pb abundance which is not detected. Additionally, there remain serious problems in explaining the observed pulsational properties of this source.

preprint2011arXiv

R-Process Nucleosynthesis in Dynamically Ejected Matter of Neutron Star Mergers

Although the rapid neutron-capture process, or r-process, is fundamentally important for explaining the origin of approximately half of the stable nuclei with A > 60, the astrophysical site of this process has not been identified yet. Here we study r-process nucleosynthesis in material that is dynamically ejected by tidal and pressure forces during the merging of binary neutron stars (NSs) and within milliseconds afterwards. For the first time we make use of relativistic hydrodynamical simulations of such events, defining consistently the conditions that determine the nucleosynthesis, i.e., neutron enrichment, entropy, early density evolution and thus expansion timescale, and ejecta mass. We find that 10^{-3}-10^{-2} solar masses are ejected, which is enough for mergers to be the main source of heavy (A > 140) galactic r-nuclei for merger rates of some 10^{-5} per year. While asymmetric mergers eject 2-3 times more mass than symmetric ones, the exact amount depends weakly on whether the NSs have radii of ~15 km for a &#34;stiff&#34; nuclear equation of state (EOS) or ~12 km for a &#34;soft&#34; EOS. R-process nucleosynthesis during the decompression becomes largely insensitive to the detailed conditions because of efficient fission recycling, producing a composition that closely follows the solar r-abundance distribution for nuclei with mass numbers A > 140. Estimating the light curve powered by the radioactive decay heating of r-process nuclei with an approximative model, we expect high emission in the B-V-R bands for 1-2 days with potentially observable longer duration in the case of asymmetric mergers because of the larger ejecta mass.

preprint2011arXiv

Standard Big-Bang Nucleosynthesis up to CNO with an improved extended nuclear network

Primordial or Big Bang nucleosynthesis (BBN) is one of the three strong evidences for the Big- Bang model together with the expansion of the Universe and the Cosmic Microwave Background radiation. In this study, we improve the standard BBN calculations taking into account new nuclear physics analyses and we enlarge the nuclear network until Sodium. This is, in particular, important to evaluate the primitive value of CNO mass fraction that could affect Population III stellar evolution. For the first time we list the complete network of more than 400 reactions with references to the origin of the rates, including \approx 270 reaction rates calculated using the TALYS code. Together with the cosmological light elements, we calculate the primordial Beryllium, Boron, Carbon, Nitrogen and Oxygen nuclei. We performed a sensitivity study to identify the important reactions for CNO, 9Be and Boron nucleosynthesis. We reevaluated those important reaction rates using experimental data and/or theoretical evaluations. The results are compared with precedent calculations: a primordial Beryllium abundance increase by a factor of 4 compared to its previous evaluation, but we note a stability for B/H and for the CNO/H abundance ratio that remains close to its previous value of 0.7 \times 10-15. On the other hand, the extension of the nuclear network has not changed the 7Li value, so its abundance is still 3-4 times greater than its observed spectroscopic value.

preprint2004arXiv

Thorium-rich halo star HD221170: further evidence against the universality of the r-process

We report the abundance determination in the atmosphere of the bright halo star HD221170. The spectra were taken with the Terskol Observatory&#39;s 2.0-m telescope with a resolution R=45000 and signal-to-noise ratio up to 250 in the wavelength region 3638-10275 Å. The adopted atmospheric parameters correspond to an effective temperature \Tef=4475 K, a surface gravity \lgg=1.0, a microturbulent velocity \vmi=1.7 \kms, and a macroturbulent velocity \vma=4 \kms. The abundances of 43 chemical elements were determined with the method of spectrum synthesis. The large overabundances (by 1 dex relative to iron) of elements with Z$>38$ are shown to follow the same pattern as the solar r-abundances. The present HD221170 analysis confirms the non-universality of the r-process, or more exactly the observation that the astrophysical sites hosting the r-process do not always lead to a unique relative abundance distribution for the bulk Ba to Hg elements, the Pb-peak elements, and the actinides.

preprint2003arXiv

Pre-suprenova evolution of rotating massive stars

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.