Source author record

R. Hirschi

R. Hirschi appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

28works
10topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

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

Published work

28 published item(s)

preprint2022arXiv

Dynamics in a stellar convective layer and at its boundary: Comparison of five 3D hydrodynamics codes

Our ability to predict the structure and evolution of stars is in part limited by complex, 3D hydrodynamic processes such as convective boundary mixing. Hydrodynamic simulations help us understand the dynamics of stellar convection and convective boundaries. However, the codes used to compute such simulations are usually tested on extremely simple problems and the reliability and reproducibility of their predictions for turbulent flows is unclear. We define a test problem involving turbulent convection in a plane-parallel box, which leads to mass entrainment from, and internal-wave generation in, a stably stratified layer. We compare the outputs from the codes FLASH, MUSIC, PPMSTAR, PROMPI, and SLH, which have been widely employed to study hydrodynamic problems in stellar interiors. The convection is dominated by the largest scales that fit into the simulation box. All time-averaged profiles of velocity components, fluctuation amplitudes, and fluxes of enthalpy and kinetic energy are within $\lesssim 3σ$ of the mean of all simulations on a given grid ($128^3$ and $256^3$ grid cells), where $σ$ describes the statistical variation due to the flow's time dependence. They also agree well with a $512^3$ reference run. The $128^3$ and $256^3$ simulations agree within $9\%$ and $4\%$, respectively, on the total mass entrained into the convective layer. The entrainment rate appears to be set by the amount of energy that can be converted to work in our setup and details of the small-scale flows in the boundary layer seem to be largely irrelevant. Our results lend credence to hydrodynamic simulations of flows in stellar interiors. We provide in electronic form all outputs of our simulations as well as all information needed to reproduce or extend our study.

preprint2022arXiv

Realistic 3D hydrodynamics simulations find significant turbulent entrainment in massive stars

Our understanding of stellar structure and evolution coming from one-dimensional (1D) stellar models is limited by uncertainties related to multi-dimensional processes taking place in stellar interiors. 1D models, however, can now be tested and improved with the help of detailed three-dimensional (3D) hydrodynamics models, which can reproduce complex multi-dimensional processes over short timescales, thanks to the recent advances in computing resources. Among these processes, turbulent entrainment leading to mixing across convective boundaries is one of the least understood and most impactful. Here we present the results from a set of hydrodynamics simulations of the neon-burning shell in a massive star, and interpret them in the framework of the turbulent entrainment law from geophysics. Our simulations differ from previous studies in their unprecedented degree of realism in reproducing the stellar environment. Importantly, the strong entrainment found in the simulations highlights the major flaws of the current implementation of convective boundary mixing in 1D stellar models. This study therefore calls for major revisions of how convective boundaries are modelled in 1D, and in particular the implementation of entrainment in these models. This will have important implications for supernova theory, nucleosynthesis, neutron stars and black holes physics.

preprint2022arXiv

The impact of $^{17}$O$+α$ reaction rate uncertainties on the s-process in rotating massive stars

Massive stars are crucial to galactic chemical evolution for elements heavier than iron. Their contribution at early times in the evolution of the Universe, however, is unclear due to poorly constrained nuclear reaction rates. The competing $^{17}$O($α,γ$)$^{21}$Ne and $^{17}$O($α,n$)$^{20}$Ne reactions strongly impact weak s-process yields from rotating massive stars at low metallicities. Abundant $^{16}$O absorbs neutrons, removing flux from the s-process, and producing $^{17}$O. The $^{17}$O($α,n$)$^{20}$Ne reaction releases neutrons, allowing continued s-process nucleosynthesis, if the $^{17}$O($α,γ$)$^{21}$Ne reaction is sufficiently weak. While published rates are available, they are based on limited indirect experimental data for the relevant temperatures and, more importantly, no uncertainties are provided. The available nuclear physics has been evaluated, and combined with data from a new study of astrophysically relevant $^{21}$Ne states using the $^{20}$Ne($d,p$)$^{21}$Ne reaction. Constraints are placed on the ratio of the ($α,n$)/($α,γ$) reaction rates with uncertainties on the rates provided for the first time. The new rates favour the ($α,n$) reaction and suggest that the weak s-process in rotating low-metallicity stars is likely to continue up to barium and, within the computed uncertainties, even to lead.

preprint2022arXiv

The p-process in exploding rotating massive stars

The p-process nucleosynthesis can explain proton-rich isotopes that are heavier than iron, which are observed in the Solar System, but discrepancies still persist and important questions concerning the astrophysical site(s) of the p-process remain unanswered. We investigate how the p-process operates in exploding rotating massive stars that have experienced an enhanced s-process nucleosynthesis during their life through rotational mixing. We computed 25 $M_{\odot}$ stellar models at a metallicity of $Z=10^{-3}$ with different initial rotation velocities and rates for the uncertain $^{17}$O($α$,$γ$)$^{21}$Ne reaction. The nucleosynthesis calculation, followed with a network of 737 isotopes, was coupled to stellar evolution, and the p-process nucleosynthesis was calculated in post-processing during both the final evolutionary stages and spherical explosions of various energies. In our models, the p-nuclides are mainly synthesized during the explosion, but not much during the ultimate hydrostatic burning stages. The p-process yields mostly depend on the initial number of trans-iron seeds, which in turn depend on the initial rotation. We found that the impact of rotation on the p-process is comparable to the impact of rotation on the s-process. From no to fast rotation, the s-process yields of nuclides with mass number $A<140$ increase by $3-4$ dex, and so do the p-process yields. Fast rotation with a lower $^{17}$O($α,γ$) rate significantly produces s- and p-nuclides with $A\geq140$. Our results suggest that the contribution of core-collapse supernovae from massive stars to the solar (and Galactic) p-nuclei has been underestimated in the past, and more specifically, that the contribution from massive stars with sub-solar metallicities may even dominate. A more detailed study including stellar models with a wide range of masses and metallicities remains to be performed.

preprint2022arXiv

Very massive star winds as sources of the short-lived radioactive isotope $^{26}$Al

The $^{26}$Al short-lived radioactive nuclide is the source of the observed galactic diffuse $γ$-ray emission at 1.8 MeV. While different sources of $^{26}$Al have been explored, such as AGB stars, massive stars winds, and supernovae, the contribution of very massive stars has never been studied. We study the stellar wind contribution of very massive stars, i.e stars with initial masses between 150 and 300 M$_\odot$, to the enrichment in $^{26}$Al of the galactic interstellar medium. We discuss the production of $^{26}$Al by studying rotating and non-rotating very massive stellar models with initial masses between 150 and 300 M$_\odot$ for metallicities Z=0.006, 0.014, and 0.020. We confront this result to a simple Milky Way model taking into account both the metallicity and the star formation rate gradients. We obtain that very massive stars in the Z=0.006-0.020 metallicity range might be very significant contributors to the $^{26}$Al enrichment of the interstellar medium. Typically, the contribution of the winds of massive stars to the total quantity of $^{26}$Al in the Galaxy increases by 150\% when very massive stars are considered. Very massive stars, despite their rarity, might be important contributors to $^{26}$Al and overall very important actors for nucleosynthesis in the Galaxy.

preprint2021arXiv

Constraints on stellar rotation from the evolution of Sr and Ba in the Galactic halo

Recent studies show that the chemical evolution of Sr and Ba in the Galaxy can be explained if different production sites, hosting r- and s-processes, are taken into account. However, the question of unambiguously identifying these sites is still unsolved. Massive stars are shown to play an important role in the production of s-material if rotation is considered. In this work, we study in detail the contribution of rotating massive stars to the production of Sr and Ba, in order to explain their chemical evolution, but also to constrain the rotational behaviour of massive stars. A stochastic chemical evolution model was employed to reproduce the enrichment of the Galactic halo. We developed new methods for model-data comparison which help to objectively compare the stochastic results to the observations. We employed these methods to estimate the value of free parameters which describe the rotation of massive stars, assumed to be dependent on the stellar metallicity. We constrain the parameters using the observations for Sr and Ba. Employing these parameters for rotating massive stars in our stochastic model, we are able to correctly reproduce the chemical evolution of Sr and Ba, but also Y, Zr and La. The data supports a decrease of both the mean rotational velocities and their dispersion with increasing metallicity. Our results show that a metallicity-dependent rotation is a necessary assumption to explain the s-process in massive stars. Our novel methods of model-data comparison represent a promising tool for future galactic chemical evolution studies.

preprint2020arXiv

The evolutionary roads leading to low effective spins, high black hole masses, and O1/O2 rates of LIGO/Virgo binary black holes

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.

preprint2020arXiv

The Formation of a 70 Msun Black Hole at High Metallicity

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.

preprint2019arXiv

Impact of Uncertainties in Astrophysical Reaction Rates on Nucleosynthesis in the $νp$ Process

The $ν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 $ν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 $ν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$α$ and $^{56}$Ni(n,p)$^{56}$Co rates at temperatures $T>3$ GK.

preprint2018arXiv

The s-process nucleosynthesis in low mass stars: impact of the uncertainties in the nuclear physics determined by Monte Carlo variations

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

preprint2016arXiv

Looking for imprints of the first stellar generations in metal-poor bulge field stars

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.

preprint2016arXiv

NuGrid stellar data set I. Stellar yields from H to Bi for stars with metallicities Z = 0.02 and Z = 0.01

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.

preprint2016arXiv

Uncertainties in the production of $p$ nuclei in massive stars obtained from Monte Carlo variations

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.

preprint2015arXiv

Carbon-rich presolar grains from massive stars. Subsolar 12C/13C and 14N/15N ratios and the mystery of 15N

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.

preprint2015arXiv

The role of neutron star mergers in the chemical evolution of the Galactic halo

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]

preprint2015arXiv

Theory Considerations For Nucleosynthesis Beyond Fe With Special Emphasis On p-Nuclei In Massive 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.

preprint2014arXiv

High-resolution abundance analysis of red giants in the globular cluster NGC 6522

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.

preprint2014arXiv

Production of 92Nb, 92Mo, and 146Sm in the gamma-process in SNIa

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.

preprint2013arXiv

Production of carbon-rich presolar grains from massive stars

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($α$,$γ$)$^{16}$O reaction relative to the rest of the $α$-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.

preprint2013arXiv

Silicon carbide grains of type C provide evidence for the production of the unstable isotope $^{32}$Si in supernovae

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 ($τ$$_{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.

preprint2012arXiv

Mass-loss predictions for evolved very metal-poor massive stars

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

preprint2012arXiv

The 12C + 12C reaction and the impact on nucleosynthesis in massive stars

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.

preprint2012arXiv

The effect of 12C + 12C rate uncertainties on the evolution and nucleosynthesis of massive stars

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

preprint2012arXiv

Very Massive Stars and the Eddington Limit

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.

preprint2011arXiv

Non-standard s process in low metallicity massive rotating stars

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$(α,γ)$ 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.

preprint2006arXiv

GRB progenitors at low metallicities

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.

preprint2006arXiv

Stellar evolution of massive stars at very low metallicities

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 ($ν$p-process) able to contribute to the nucleosynthesis of elements with A > 64. (Abridged)