Source author record

Marco Pignatari

Marco Pignatari 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

36works
8topics
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

36 published item(s)

preprint2026arXiv

Impact of newly measured 26Al(n, p)26Mg and 26Al(n, α)23Na reaction rates on the nucleosynthesis of 26Al in stars

The cosmic production of the short-lived radioactive nuclide 26Al is crucial for our understanding of the evolution of stars and galaxies. However, simulations of the stellar sites producing 26Al are still weakened by significant nuclear uncertainties. We re-evaluate the 26Al(n, p)26Mg, and 26Al(n, α)23Na ground state reactivities from 0.01 GK to 10 GK, based on the recent n TOF measurement combined with theoretical predictions and a previous measurement at higher energies, and test their impact on stellar nucleosynthesis. We computed the nucleosynthesis of low- and high-mass stars using the Monash nucleosynthesis code, the NuGrid mppnp code, and the FUNS stellar evolutionary code. Our low-mass stellar models cover the 2-3 Msun mass range with metallicities between Z = 0.01 and 0.02, their predicted 26Al/27Al ratios are compared to 62 meteoritic SiC grains. For high-mass stars, we test our reactivities on two 15 Msun models with Z = 0.006 and 0.02. The new reactivities allow low-mass AGB stars to reproduce the full range of 26Al/27Al ratios measured in SiC grains. The final 26Al abundance in high-mass stars, at the point of highest production, varies by a factor of 2.4 when adopting the upper, or lower limit of our rates. However, stellar uncertainties still play an important role in both mass regimes. The new reactivities visibly impact both low- and high-mass stars nucleosynthesis and allow a general improvement in the comparison between stardust SiC grains and low-mass star models. Concerning explosive nucleosynthesis, an improvement of the current uncertainties between T9~0.3 and 2.5 is needed for future studies.

preprint2026arXiv

The NuGrid AGB Evolution and Nucleosynthesis Data Set

Asymptotic Giant Branch (AGB) stars play a key role in the chemical evolution of galaxies. These stars are the fundamental stellar site for the production of light elements such as C, N and F, and half of the elements heavier than Fe via the slow neutron capture process (s-process). Hence, detailed computational models of AGB stars' evolution and nucleosynthesis are essential for galactic chemical evolution. In this work, we discuss the progress in updating the NuGrid data set of AGB stellar models and abundance yields. All stellar models have been computed using the MESA stellar evolution code, coupled with the post-processing mppnp code to calculate the full nucleosynthesis. The final data set will include the initial masses Mini/Msun = 1, 1.65, 2, 3, 4, 5, 6 and 7 for initial metallicities Z = 0.0001, 0.001, 0.006, 0.01, 0.02 and 0.03. Observed s-process abundances on the surfaces of evolved stars as well as the typical light elements in the composition of H-deficient post-AGB stars are reproduced. A key short-term goal is to complete and expand the AGB stars data set for the full metallicity range. Chemical yield tables are provided for the available models.

preprint2024arXiv

Spectacular nucleosynthesis from early massive stars

Stars formed with initial mass over 50 Msun are very rare today, but they are thought to be more common in the early universe. The fates of those early, metal-poor, massive stars are highly uncertain. Most are expected to directly collapse to black holes, while some may explode as a result of rotationally powered engines or the pair-creation instability. We present the chemical abundances of J0931+0038, a nearby low-mass star identified in early followup of SDSS-V Milky Way Mapper, which preserves the signature of unusual nucleosynthesis from a massive star in the early universe. J0931+0038 has relatively high metallicity ([Fe/H] = -1.76 +/- 0.13) but an extreme odd-even abundance pattern, with some of the lowest known abundance ratios of [N/Fe], [Na/Fe], [K/Fe], [Sc/Fe], and [Ba/Fe]. The implication is that a majority of its metals originated in a single extremely metal-poor nucleosynthetic source. An extensive search through nucleosynthesis predictions finds a clear preference for progenitors with initial mass > 50 Msun, making J0931+0038 one of the first observational constraints on nucleosynthesis in this mass range. However the full abundance pattern is not matched by any models in the literature. J0931+0038 thus presents a challenge for the next generation of nucleosynthesis models and motivates study of high-mass progenitor stars impacted by convection, rotation, jets, and/or binary companions. Though rare, more examples of unusual early nucleosynthesis in metal-poor stars should be found in upcoming large spectroscopic surveys.

preprint2022arXiv

Comparison between core-collapse supernova nucleosynthesis and meteoric stardust grains: investigating magnesium, aluminium, and chromium

Isotope variations of nucleosynthetic origin among Solar System's solid samples are well documented, yet the origin of these variations is still uncertain. The observed variability of \iso{54}Cr among materials formed in different regions of the proto-planetary disk has been attributed to variable amounts of presolar chromium-rich oxide (chromite) grains, which exist within the meteoritic stardust inventory and most likely originated from some type of supernova explosions. To investigate if core-collapse supernovae (CCSNe) could be the site of origin of these grains, we analyse yields of CCSN models of stars with initial mass 15, 20 and 25 M$_{\odot}$, and solar metallicity. We present an extensive abundance data set of the Cr, Mg, and Al isotopes as a function of enclosed mass. We find cases in which the explosive C-ashes produce a composition in good agreement with the observed \iso{54}Cr/\iso{52}Cr and \iso{53}Cr/\iso{52}Cr ratios as well as the \iso{50}Cr/\iso{52}Cr ratios. Taking into account that the signal at atomic mass 50 could also originate from \iso{50}Ti, the ashes of explosive He-burning also match the observed ratios. Addition of material from the He ashes (enriched in Al and Cr relative to Mg to simulate the make-up of chromite grains) to the Solar System composition may reproduce the observed correlation between Mg and Cr anomalies, while material from the C-ashes does not present significant Mg anomalies together with Cr isotopic variations. In all cases, non-radiogenic, stable Mg isotope variations dominate over the variations expected from \iso{26}Al.

preprint2022arXiv

Origin of Plutonium-244 in the Early Solar System

We investigate the origin in the early Solar System of the short-lived radionuclide 244Pu (with a half life of 80 Myr) produced by the rapid (r) neutron-capture process. We consider two large sets of r-process nucleosynthesis models and analyse if the origin of 244Pu in the ESS is consistent with that of the other r and slow (s) neutron-capture process radioactive nuclei. Uncertainties on the r-process models come from both the nuclear physics input and the astrophysical site. The former strongly affects the ratios of isotopes of close mass (129I/127I, 244Pu/238U, and 247Pu/235U). The 129I/247Cm ratio, instead, which involves isotopes of a very different mass, is much more variable than those listed above and is more affected by the physics of the astrophysical site. We consider possible scenarios for the evolution of the abundances of these radioactive nuclei in the galactic interstellar medium and verify under which scenarios and conditions solutions can be found for the origin of 244Pu that are consistent with the origin of the other isotopes. Solutions are generally found for all the possible different regimes controlled by the interval ($δ$) between additions from the source to the parcel of interstellar medium gas that ended up in the Solar System, relative to decay timescales. If r-process ejecta in interstellar medium are mixed within a relatively small area (leading to a long $δ$), we derive that the last event that explains the 129I and 247Cm abundances in the early Solar System can also account for the abundance of 244Pu. Due to its longer half life, however, 244Pu may have originated from a few events instead of one only. If r-process ejecta in interstellar medium are mixed within a relatively large area (leading to a short $δ$), we derive that the time elapsed from the formation of the molecular cloud to the formation of the Sun was 9-16 Myr.

preprint2022arXiv

Radioactive nuclei in the early Solar System: analysis of the 15 isotopes produced by core-collapse supernovae

Short-lived radioactive isotopes (SLRs) with half-lives between 0.1 to 100 Myr can be used to probe the origin of the Solar System. In this work, we examine the core-collapse supernovae production of the 15 SLRs produced: $^{26}$Al, $^{36}$Cl, $^{41}$Ca, $^{53}$Mn, $^{60}$Fe, $^{92}$Nb, $^{97}$Tc, $^{98}$Tc, $^{107}$Pd, $^{126}$Sn, $^{129}$I, $^{135}$Cs, $^{146}$Sm, $^{182}$Hf, and $^{205}$Pb. We probe the impact of the uncertainties of the core-collapse explosion mechanism by examining a collection of 62 core-collapse models with initial masses of 15, 20, and 25M$_{\odot}$, explosion energies between 3.4$\times$10$^{50}$ and 1.8$\times$10$^{52}$ ergs and compact remnant masses between 1.5M$_{\odot}$and 4.89M$_{\odot}$. We identify the impact of both explosion energy and remnant mass on the final yields of the SLRs. Isotopes produced within the innermost regions of the star, such as $^{92}$Nb and $^{97}$Tc, are the most affected by the remnant mass, $^{92}$Nb varying by five orders of magnitude. Isotopes synthesised primarily in explosive C-burning and explosive He-burning, such as $^{60}$Fe, are most affected by explosion energies. $^{60}$Fe increases by two orders of magnitude from the lowest to the highest explosion energy in the 15M$_{\odot}$model. The final yield of each examined SLR is used to compare to literature models.

preprint2022arXiv

The RADIOSTAR Project

Radioactive nuclei are the key to understanding the circumstances of the birth of our Sun because meteoritic analysis has proven that many of them were present at that time. Their origin, however, has been so far elusive. The ERC-CoG-2016 RADIOSTAR project is dedicated to investigating the production of radioactive nuclei by nuclear reactions inside stars, their evolution in the Milky Way Galaxy, and their presence in molecular clouds. So far, we have discovered that: (i) radioactive nuclei produced by slow ($^{107}$Pd and $^{182}$Hf) and rapid ($^{129}$I and $^{247}$Cm) neutron captures originated from stellar sources - asymptotic giant branch (AGB) stars and compact binary mergers, respectively - within the galactic environment that predated the formation of the molecular cloud where the Sun was born; (ii) the time that elapsed from the birth of the cloud to the birth of the Sun was of the order of 10$^7$ years, and (iii) the abundances of the very short-lived nuclei $^{26}$Al, $^{36}$Cl, and $^{41}$Ca can be explained by massive star winds in single or binary systems, if these winds directly polluted the early Solar System. Our current and future work, as required to finalise the picture of the origin of radioactive nuclei in the Solar System, involves studying the possible origin of radioactive nuclei in the early Solar System from core-collapse supernovae, investigating the production of $^{107}$Pd in massive star winds, modelling the transport and mixing of radioactive nuclei in the galactic and molecular cloud medium, and calculating the galactic chemical evolution of $^{53}$Mn and $^{60}$Fe and of the p-process isotopes $^{92}$Nb and $^{146}$Sm.

preprint2021arXiv

129I and 247Cm in Meteorites Constrain the Last Astrophysical Source of Solar r-process Elements

The composition of the early Solar System can be inferred from meteorites. Many elements heavier than iron were formed by the rapid neutron-capture process (r process), but the astrophysical sources where this occurred remain poorly understood. We demonstrate that the near-identical half-lives ($\simeq$ 15.6 Myr) of the radioactive r-process nuclei 129I and 247Cm preserve their ratio, irrespective of the time between production and incorporation into the Solar System. We constrain the last r-process source by comparing the measured meteoritic 129I / 247Cm = 438 $\pm$ 184 to nucleosynthesis calculations based on neutron star merger and magneto-rotational supernova simulations. Moderately neutron-rich conditions, often found in merger disk ejecta simulations, are most consistent with the meteoritic value. Uncertain nuclear physics data limit our confidence in this conclusion.

preprint2021arXiv

A Large-scale Approach to Modelling Molecular Biosignatures: The Diatomics

This work presents the first steps to modelling synthetic rovibrational spectra for all molecules of astrophysical interest using a new approach implemented in the Prometheus code. The goal is to create a new comprehensive source of first-principles molecular spectra, thus bridging the gap for missing data to help drive future high-resolution studies. Our primary application domain is for molecules identified as signatures of life in planetary atmospheres (biosignatures), but our approach is general and can be applied to other systems. In this work we evaluate the accuracy of our method by studying four diatomic molecules H$_2$, O$_2$, N$_2$ and CO, all of which have well-known spectra. Prometheus uses the Transition-Optimised Shifted Hermite (TOSH) theory to account for anharmonicity for the fundamental $ν=0 \rightarrow ν=1$ band, along with thermal profile modeling for the rotational transitions. To this end, we expand TOSH theory to enable the modeling of rotational constants. We show that our simple model achieves results that are a better approximation of the real spectra than those produced through a harmonic approach. We compare our results with high-resolution HITRAN and ExoMol spectral data. We find that modelling accuracy tends to diminish for rovibrational transition away from the band origin, thus highlighting the need for the theory to be further adapted.

preprint2020arXiv

Heavy Elements Nucleosynthesis On Accreting White Dwarfs: building seeds for the p-process

The origin of the proton-rich trans-iron isotopes in the solar system is still uncertain. Single-degenerate thermonuclear supernovae (SNIa) with n-capture nucleosynthesis seeds assembled in the external layers of the progenitor's rapidly accreting white dwarf phase may produce these isotopes. We calculate the stellar structure of the accretion phase of five white dwarf models with initial masses >~ 0.85Msun using the stellar code MESA. The near-surface layers of the 1, 1.26, 1.32 and 1.38Msun models are most representative of the regions in which the bulk of the p nuclei are produced during SNIa explosions, and for these models we also calculate the neutron-capture nucleosynthesis in the external layers. Contrary to previous rapidly-accreting white dwarf models at lower mass, we find that the H-shell ashes are the main site of n-capture nucleosynthesis. We find high neutron densities up to several 10^15 cm^-3 in the most massive WDs. Through the recurrence of the H-shell ashes these intermediate neutron densities can be sustained effectively for a long time leading to high neutron exposures with a strong production up to Pb. Both the neutron density and the neutron exposure increase with increasing the mass of the accreting WD. Finally, the SNIa nucleosynthesis is calculated using the obtained abundances as seeds. We obtain solar to super-solar abundances for p-nuclei with A>96. Our models show that SNIa are a viable p-process production site.

preprint2020arXiv

Origin of large meteoritic SiC stardust grains in metal-rich AGB stars

Stardust grains that originated in ancient stars and supernovae are recovered from meteorites and carry the detailed composition of their astronomical sites of origin. We present evidence that the majority of large ($μ$m-sized) meteoritic silicon carbide (SiC) grains formed in C-rich asymptotic giant branch (AGB) stars that were more metal-rich than the Sun. In the framework of the slow neutron-captures (the s process) that occurs in AGB stars the lower-than-solar 88Sr/86Sr isotopic ratios measured in the large SiC grains can only be accompanied by Ce/Y elemental ratios that are also lower than solar, and predominately observed in metal-rich barium stars - the binary companions of AGB stars. Such an origin suggests that these large grains represent the material from high-metallicity AGB stars needed to explain the s-process nucleosynthesis variations observed in bulk meteorites (Ek et al. 2020). In the outflows of metal-rich, C-rich AGB stars SiC grains are predicted to be small ($\simeq$ 0.2 $μ$m-sized); large ($\simeq$ $μ$m-sized) SiC grains can grow if the number of dust seeds is two to three orders of magnitude lower than the standard value of $10^{-13}$ times the number of H atoms. We therefore predict that with increasing metallicity the number of dust seeds might decrease, resulting in the production of larger SiC grains.

preprint2019arXiv

Chromium Nucleosynthesis and Silicon-Carbon Shell Mergers in Massive Stars

We analyze the production of the element Cr in galactic chemical evolution (GCE) models using the NuGrid nucleosynthesis yields set. We show that the unusually large [Cr/Fe] abundance at [Fe/H] $\approx 0$ reported by previous studies using those yields and predicted by our Milky Way model originates from the merging of convective Si-burning and C-burning shells in a 20 $M_\odot$ model at metallicity $Z=0.01$, about an hour before the star explodes. This merger mixes the incomplete burning material in the Si shell, including $^{51}$V and $^{52}$Cr, out to the edge of the carbon/oxygen (CO) core. The adopted supernova model ejects the outer 2 $M_\odot$ of the CO core, which includes a significant fraction of the Cr-rich material. When including this 20 $M_\odot$ model at $Z=0.01$ in the yields interpolation scheme of our GCE model for stars in between 15 and 25 $M_\odot$, we overestimate [Cr/Fe] by an order of magnitude at [Fe/H] $\approx$ 0 relative to observations in the Galactic disk. This raises a number of questions regarding the occurrence of Si-C shell mergers in nature, the accuracy of different simulation approaches, and the impact of such mergers on the pre-supernova structure and explosion dynamics. According to the conditions in this 1D stellar model, the substantial penetration of C-shell material into the Si-shell could launch a convective-reactive global oscillation, if a merger does take place. In any case, GCE provides stringent constraints on the outcome of this stellar evolution phase.

preprint2019arXiv

The Nucleosynthetic Yields of Core-Collapse Supernovae, prospects for the Next Generation of Gamma-Ray Astronomy

Though the neutrino-driven convection model for the core-collapse explosion mechanism has received strong support in recent years, there are still many uncertainties in the explosion parameters -- such as explosion energy, remnant mass, and end-of-life stellar abundances as initial conditions. Using a broad set of spherically symmetric core-collapse simulations we examine the effects of these key parameters on explosive nucleosynthesis and final explosion yields. Post-bounce temperature and density evolution of ZAMS 15, 20, and 25 solar mass progenitors are post-processed through the Nucleosynthesis Grid (NuGrid) nuclear network to obtain detailed explosive yields. In particular, this study focuses on radio-isotopes that are of particular interest to the next generation of gamma-ray astronomical observations; 43K, 47Ca, 44Sc, 47Sc, 48V, 48Cr, 51Cr, 52Mn, 59Fe, 56Co, 57Co, 57Ni. These nuclides may be key in advancing our understanding of the inner workings of core-collapse supernovae by probing the parameters of the explosion engine. We find that the isotopes that are strong indicators of explosion energy are 43K, 47Ca, 44Sc, 47Sc, and 59Fe, those that are dependent on the progenitor structure are 48V, 51Cr, and 57Co, and those that probe neither are 48Cr, 52Mn, 57Ni, and 56Co. We discuss prospects of observing these radionuclides in supernova remnants.

preprint2018arXiv

H-He Shell Interactions and Nucleosynthesis in Massive Population III Stars

We report on our ongoing investigation into the nucleosynthetic and hydrodynamic nature of mixing at the interface between the H- and He-convection zones in massive Pop III stars. Studying a grid of 26 1D stellar evolution simulations with different mixing assumptions, we find that H-He interactions occur in 23/26 cases. We demonstrate the nucleosynthesis expected in a H-He interaction in an 80M$_\odot$. Finally, we describe our progress in simulating a Pop III double convection zone in the PPMStar hydrodynamics code.

preprint2016arXiv

Origin of the p-process radionuclides 92Nb and 146Sm in the early Solar System and inferences on the birth of the Sun

The abundances of 92Nb and 146Sm in the early Solar System are determined from meteoritic analysis and their stellar production is attributed to the p process. We investigate if their origin from thermonuclear supernovae deriving from the explosion of white dwarfs with mass above the Chandrasekhar limit is in agreement with the abundance of 53Mn, another radionuclide present in the early Solar System and produced in the same events. A consistent solution for 92Nb and 53Mn cannot be found within the current uncertainties and requires that the 92Nb/92Mo ratio in the early Solar System is at least 50% lower than the current nominal value, which is outside its present error bars. A different solution is to invoke another production site for 92Nb, which we find in the alpha-rich freezeout during core-collapse supernovae from massive stars. Whichever scenario we consider, we find that a relatively long time interval of at least ~10 Myr must have elapsed from when the star-forming region where the Sun was born was isolated from the interstellar medium and the birth of the Sun. This is in agreement with results obtained from radionuclides heavier than iron produced by neutron captures and lends further support to the idea that the Sun was born in a massive star-forming region together with many thousands of stellar siblings.

preprint2016arXiv

Stellar Origins of Extremely $^{\text{13}}C$- and $^{15}N$-enriched Presolar SiC Grains: Novae or Supernovae?

Extreme excesses of $^{13}C$ ($^{12}C$/$^{13}C$<10) and $^{15}N$ ($^{14}N$/$^{15}N$<20) in rare presolar SiC grains have been considered diagnostic of an origin in classical novae, though an origin in core collapse supernovae (CCSNe) has also been proposed. We report C, N, and Si isotope data for 14 submicron- to micron-sized $^{13}C$- and $^{15}N$-enriched presolar SiC grains ($^{12}C$/$^{13}C$<16 and $^{14}N$/$^{15}N$<~100) from Murchison, and their correlated Mg-Al, S, and Ca-Ti isotope data when available. These grains are enriched in $^{13}C$ and $^{15}N$, but with quite diverse Si isotopic signatures. Four grains with $^{29,30}Si$ excesses similar to those of type C SiC grains likely came from CCSNe, which experienced explosive H burning occurred during explosions. The independent coexistence of proton- and neutron-capture isotopic signatures in these grains strongly supports heterogeneous H ingestion into the He shell in pre-supernovae. Two of the seven putative nova grains with $^{30}Si$ excesses and $^{29}Si$ depletions show lower-than-solar $^{34}S$/$^{32}S$ ratios that cannot be explained by classical nova nucleosynthetic models. We discuss these signatures within the CCSN scenario. For the remaining five putative nova grains, both nova and supernova origins are viable because explosive H burning in the two stellar sites could result in quite similar proton-capture isotopic signatures. Three of the grains are sub-type AB grains that are also $^{13}C$ enriched, but have a range of higher $^{14}N$/$^{15}N$. We found that $^{15}N$-enriched AB grains (~50<$^{14}N$/$^{15}N$<~100) have distinctive isotopic signatures compared to putative nova grains, such as higher $^{14}N$/$^{15}N$, lower $^{26}Al$/$^{27}Al$, and lack of $^{30}Si$ excess, indicating weaker proton-capture nucleosynthetic environments.

preprint2016arXiv

The Diverse Origins of Neutron-Capture Elements in the Metal-Poor Star HD 94028: Possible Detection of Products of i-process Nucleosynthesis

We present a detailed analysis of the composition and nucleosynthetic origins of the heavy elements in the metal-poor ([Fe/H]=-1.62+/-0.09) star HD94028. Previous studies revealed that this star is mildly enhanced in elements produced by the slow neutron-capture process (s-process; e.g., [Pb/Fe]=+0.79+/-0.32) and rapid neutron-capture process (r-process; e.g., [Eu/Fe]=+0.22+/-0.12), including unusually large molybdenum ([Mo/Fe]=+0.97+/-0.16) and ruthenium ([Ru/Fe]=+0.69+/-0.17) enhancements. However, this star is not enhanced in carbon ([C/Fe]=-0.06+/-0.19). We analyze an archival near-ultraviolet spectrum of HD94028, collected using the Space Telescope Imaging Spectrograph on board the Hubble Space Telescope, and other archival optical spectra collected from ground-based telescopes. We report abundances or upper limits derived from 64 species of 56 elements. We compare these observations with s-process yields from low-metallicity AGB evolution and nucleosynthesis models. No combination of s- and r-process patterns can adequately reproduce the observed abundances, including the super-solar [As/Ge] ratio (+0.99+/-0.23) and the enhanced [Mo/Fe] and [Ru/Fe] ratios. We can fit these features when including an additional contribution from the intermediate neutron-capture process (i process), which perhaps operated by the ingestion of H in He-burning convective regions in massive stars, super-AGB stars, or low-mass AGB stars. Currently, only the i process appears capable of consistently producing the super-solar [As/Ge] ratios and ratios among neighboring heavy elements found in HD94028. Other metal-poor stars also show enhanced [As/Ge] ratios, hinting that operation of the i process may have been common in the early Galaxy.

preprint2016arXiv

The production of proton-rich isotopes beyond iron: The $γ$ process in stars

Beyond iron, a small fraction of the total abundances in the Solar System is made of proton-rich isotopes, the p nuclei. The clear understanding of their production is a fundamental challenge for nuclear astrophysics. The p nuclei constrain the nucleosynthesis in core-collapse and thermonuclear supernovae. The $γ$ process is the most established scenario for the production of the p nuclei, which are produced via different photodisintegration paths starting on heavier nuclei. A large effort from nuclear physics is needed to access the relevant nuclear reaction rates far from the valley of stability. This review describes the production of the heavy proton-rich isotopes by the $γ$ process in stars, and explores the state of the art of experimental nuclear physics to provide nuclear data for stellar nucleosynthesis.

preprint2016arXiv

Uncertainties in Galactic Chemical Evolution Models

We use a simple one-zone galactic chemical evolution model to quantify the uncertainties generated by the input parameters in numerical predictions, for a galaxy with properties similar to those of the Milky Way. We compiled several studies from the literature to gather the current constraints for our simulations regarding the typical value and uncertainty of seven basic parameters, which are: the lower and upper mass limits of the stellar initial mass function (IMF), the slope of the high-mass end of the stellar IMF, the slope of the delay-time distribution function of Type Ia supernovae (SNe Ia), the number of SNeIa per Msun formed, the total stellar mass formed, and the final mass of gas. We derived a probability distribution function to express the range of likely values for every parameter, which were then included in a Monte Carlo code to run several hundred simulations with randomly selected input parameters. This approach enables us to analyze the predicted chemical evolution of 16 elements in a statistical way by identifying the most probable solutions along with their 68% and 95% confidence levels. Our results show that the overall uncertainties are shaped by several input parameters that individually contribute at different metallicities, and thus at different galactic ages. The level of uncertainty then depends on the metallicity and is different from one element to another. Among the seven input parameters considered in this work, the slope of the IMF and the number of SNe Ia are currently the two main sources of uncertainty.

preprint2015arXiv

$s$-process production in rotating massive stars at solar and low metallicities

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.

preprint2015arXiv

Iron-Group Abundances in the Metal-Poor Main Sequence Turnoff Star HD~84937

We have derived new very accurate abundances of the Fe-group elements Sc through Zn (Z = 21-30) in the bright main-sequence turnoff star HD 84937, based on high-resolution spectra covering the visible and ultraviolet spectral regions. New or recent laboratory transition data for 14 species of seven elements have been used. Abundances from more than 600 lines of non-Fe species have been combined with about 550 Fe lines in HD 84937 to yield abundance ratios of high precision. The abundances have been determined from both neutral and ionized transitions, which generally are in agreement with each other. We find no substantial departures from standard LTE Saha ionization balance in this [Fe/H] = -2.32 star. Noteworthy among the abundances are: [Co/Fe] = 0.14 and [Cu/Fe] = -0.83, in agreement with past studies abundance trends in this and other low metallicity stars; and <[Sc,Ti,V/Fe]> = 0.31, which has not been noted previously. A detailed examination of scandium, titanium, and vanadium abundances in large-sample spectroscopic surveys reveals that they are positively correlated in stars with [Fe/H] < -2; HD 84937 lies at the high end of this correlation. These trends constrain the synthesis mechanisms of Fe-group elements. We also examine the GCE abundance trends of the Fe-group elements, including a new nucleosynthesis model with jet-like explosion effects.

preprint2015arXiv

Nucleosynthesis simulations for the production of the p-nuclei $^{\text{92}}$Mo and $^{\text{94}}$Mo in a Supernova type II model

We present a nucleosynthesis sensitivity study for the $γ$-process in a Supernova type II model within the NuGrid research platform. The simulations aimed at identifying the relevant local production and destruction rates for the p-nuclei of molybdenum and at determining the sensitivity of the final abundances to these rates. We show that local destruction rates strongly determine the abundance of $^{92}$Mo and $^{94}$Mo, and quantify the impact.

preprint2014arXiv

Barium Isotopic Composition of Mainstream Silicon Carbides from Murchison: Constraints for s-Process Nucleosynthesis in AGB Stars

We present barium, carbon, and silicon isotopic compositions of 38 acid-cleaned presolar SiC grains from Murchison. Comparison with previous data shows that acid washing is highly effective in removing barium contamination. Strong depletions in $δ$($^{138}$Ba/$^{136}$Ba) values are found, down to $-$400 permil, which can only be modeled with a flatter $^{13}$C profile within the $^{13}$C pocket than is normally used. The dependence of $δ$($^{138}$Ba/$^{136}$Ba) predictions on the distribution of $^{13}$C within the pocket in AGB models allows us to probe the $^{13}$C profile within the $^{13}$C pocket and the pocket mass in asymptotic giant branch (AGB) stars. In addition, we provide constraints on the $^{22}$Ne$(α,n)^{25}$Mg rate in the stellar temperature regime relevant to AGB stars, based on $δ$($^{134}$Ba/$^{136}$Ba) values of mainstream grains. We found two nominally mainstream grains with strongly negative $δ$($^{134}$Ba/$^{136}$Ba) values that cannot be explained by any of the current AGB model calculations. Instead, such negative values are consistent with the intermediate neutron capture process ($i$-process), which is activated by the Very Late Thermal Pulse (VLTP) during the post-AGB phase and characterized by a neutron density much higher than the $s$-process. These two grains may have condensed around post-AGB stars. Finally, we report abundances of two $p$-process isotopes, $^{130}$Ba and $^{132}$Ba, in single SiC grains. These isotopes are destroyed in the $s$-process in AGB stars. By comparing their abundances with respect to that of $^{135}$Ba, we conclude that there is no measurable decay of $^{135}$Cs ($t_{1/2}$= 2.3 Ma) to $^{135}$Ba in individual SiC grains, indicating condensation of barium, but not cesium into SiC grains before $^{135}$Cs decayed.

preprint2014arXiv

Chemical abundances in bright giants of the globular cluster M62 (NGC 6266)

With the exception of Terzan 5, all the Galactic globular clusters that possess significant metallicity spreads, such as omega Cen and M22, are preferentially the more luminous clusters with extended horizontal branches. Here we present radial velocities and chemical abundances for seven bright giants in the globular cluster M62, a previously little-studied cluster. With M_V = -9.18, M62 is the ninth most luminous Galactic globular cluster and has an extended horizontal branch. Within our sample, we find (i) no evidence for a dispersion in metallicity, [Fe/H], beyond the measurement uncertainties, (ii) star-to-star abundance variations for C, O, Na and Al with the usual correlations between these elements as seen in other globular clusters, and (iii) a global enrichment for the elements Zr, Ba and La at the level [X/Fe] = +0.4 dex. For elements heavier than La, the abundance ratios are consistent with the scaled-solar $r$-process distribution. Below La, the abundances are anomalous when compared to the scaled-solar s-process or r-process distributions. For these elements, the abundance signature in M62 is in agreement with predictions of the s-process from fast-rotating massive stars, although the high [Rb/Y] ratio we measure may be a challenge to this scenario.

preprint2014arXiv

Code dependencies of pre-supernova evolution and nucleosynthesis in massive stars: Evolution to the end of core helium burning

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.

preprint2014arXiv

MESA and NuGrid simulations of classical novae: CO and ONe nova nucleosynthesis

Classical novae are the result of thermonuclear flashes of hydrogen accreted by CO or ONe white dwarfs, leading eventually to the dynamic ejection of the surface layers. These are observationally known to be enriched in heavy elements, such as C, O and Ne that must originate in layers below the H-flash convection zone. Building on our previous work, we now present stellar evolution simulations of ONe novae and provide a comprehensive comparison of our models with published ones. Some of our models include exponential convective boundary mixing to account for the observed enrichment of the nova ejecta even when accreted material has a solar abundance distribution. Our models produce maximum temperature evolution profiles and nucleosynthesis yields in good agreement with models that generate enriched ejecta by assuming that the accreted material was pre-mixed. We confirm for ONe novae the result we reported previously, i.e.\ we found that $^3$He could be produced {\it in situ} in solar-composition envelopes accreted with slow rates ($\dot{M} < 10^{-10}\,M_\odot/\mbox{yr}$) by cold ($T_{\rm WD} < 10^7$ K) CO WDs, and that convection was triggered by $^3$He burning before the nova outburst in that case. In addition, we now find that the interplay between the $^3$He production and destruction in the solar-composition envelope accreted with an intermediate rate, e.g.\ $\dot{M} = 10^{-10}\,M_\odot/\mbox{yr}$, by the $1.15\,M_\odot$ ONe WD with a relatively high initial central temperature, e.g.\ $T_{\rm WD} = 15\times 10^6$ K, leads to the formation of a thick radiative buffer zone that separates the bottom of the convective envelope from the WD surface. (Abridged)

preprint2013arXiv

Evidence for Radiogenic Sulfur-32 in Type AB Presolar Silicon Carbide Grains?

We report C, Si, and S isotope measurements on 34 presolar silicon carbide grains of Type AB, characterized by 12C/13C < 10. Nitrogen, Mg-Al-, and Ca-Ti-isotopic compositions were measured on a subset of these grains. Three grains show large 32S excesses, a signature that has been previously observed for grains from supernovae (SNe). Enrichments in 32S may be due to contributions from the Si/S zone and the result of S molecule chemistry in still unmixed SN ejecta or due to incorporation of radioactive 32Si from C-rich explosive He shell ejecta. However, a SN origin remains unlikely for the three AB grains considered here, because of missing evidence for 44Ti, relatively low 26Al/27Al ratios (a few times 10-3), and radiogenic 32S along with low 12C/13C ratios. Instead, we show that born-again asymptotic giant branch (AGB) stars that have undergone a very-late thermal pulse (VLTP), known to have low 12C/13C ratios and enhanced abundances of the light s-process elements, can produce 32Si, which makes such stars attractive sources for AB grains with 32S excesses. This lends support to the proposal that at least some AB grains originate from born-again AGB stars, although uncertainties in the born-again AGB star models and possible variations of initial S-isotopic compositions in the parent stars of AB grains make it difficult to draw a definitive conclusion.

preprint2013arXiv

On the internal pollution mechanisms in the globular cluster NGC 6121 (M4): heavy-element abundances and AGB models

Globular clusters display significant variations in their light-element content, pointing to the existence of a second stellar generation formed from the ejecta of an earlier generation. The nature of these internal polluters is still a matter of debate: the two most popular scenarios indicate intermediate-mass AGB stars and fast rotating massive stars. Abundances determination for some key elements can help distinguish between these competitor candidates. We present in this paper Y abundances for a sample of 103 red giant branch stars in NGC 6121. Within measurement errors, we find that the [Y/Fe] is constant in this cluster contrary to a recent suggestion. For a subsample of six stars we also find [Rb/Fe] to be constant, consistent with previous studies showing no variation in other s-process elements. We also present a new set of stellar yields for intermediate-mass AGB stellar models of 5 and 6 solar masses, including heavy element s-process abundances. The uncertainties on the mass-loss rate, the mixing-length parameter, and the nuclear reaction rates have a major impact on the stellar abundances. Within the IM-AGB pollution scenario, the constant abundance of heavy elements inside the cluster requires a marginal s-process efficiency in IM-AGB stars. Such a constrain could still be satisfied by the present models assuming a stronger mass-loss rate. The uncertainties mentioned above are limiting the predictive power of intermediate-mass AGB models. For these reasons, at the moment we are not able to clearly rule out their role as main polluters of the second population stars in globular clusters.

preprint2013arXiv

Relics of ancient post-AGB stars in a primitive meteorite

Graphite is one of the many presolar circumstellar condensate species found in primitive meteorites. While the isotopic compositions of low-density graphite grains indicate an origin in core-collapse supernovae, some high-density grains have extreme isotopic anomalies in C, Ca and Ti, which cannot be explained by envelope predictions of asymptotic giant branch (AGB) stars or theoretical supernova models. The Ca and Ti isotopic anomalies, however, match the predictions of He-shell abundances in AGB stars. In this study, we show that the C, Ca, and Ti isotopic anomalies are consistent with nucleosynthesis predictions of the H-ingestion phase during a very late thermal pulse (VLTP) event in post-AGB stars. The low $^{12}$C/$^{13}$C isotopic ratios in these grains are a result of abundant $^{12}$C efficiently capturing the protons that are being ingested during the VLTP. Very high neutron densities of $\sim 10^{15}$ cm$^{-3}$, typical of the $i$-process, are achieved during this phase in post-AGB stars. The large $^{42,43,44}$Ca excesses in some graphite grains are indicative of neutron capture nucleosynthesis during VLTP. The comparison of VLTP nucleosynthesis calculations to the graphite data also indicate that apparent anomalies in the Ti isotopic ratios are due to large contributions from $^{46,48}$Ca, which cannot be resolved from the isobars $^{46,48}$Ti during the measurements. We conclude that presolar graphite grains with moderate to extreme Ca and Ti isotopic anomalies originate in post-AGB stars that suffer a very late thermal pulse.

preprint2013arXiv

Reproducing the observed abundances in RCB and HdC stars with post-double degenerate merger models - constraints on merger and post-merger simulations and physics processes

The R Coronae Borealis (RCB) stars are hydrogen-deficient, variable stars that are most likely the result of He-CO WD mergers. They display extremely low oxygen isotopic ratios, 16O/18O ~ 1 - 10, 12C/13C>=100, and enhancements up to 2.6dex in F and in s-process elements from Zn to La, compared to solar. These abundances provide stringent constraints on the physical processes during and after the double-degenerate merger. As shown before O-isotopic ratios observed in RCB stars cannot result from the dynamic double-degenerate merger phase, and we investigate now the role of the long-term 1D spherical post-merger evolution and nucleosynthesis based on realistic hydrodynamic merger progenitor models. We adopt a model for extra envelope mixing to represent processes driven by rotation originating in the dynamical merger. Comprehensive nucleosynthesis post-processing simulations for these stellar evolution models reproduce, for the first time, the full range of the observed abundances for almost all the elements measured in RCB stars: 16O/18O ratios between 9 and 15, C-isotopic ratios above 100, and ~1.4 - 2.35dex F enhancements, along with enrichments in s-process elements. The nucleosynthesis processes in our models constrain the length and temperature in the dynamic merger shell-of-fire feature as well as the envelope mixing in the post-merger phase. s-process elements originate either in the shell-of-fire merger feature or during the post-merger evolution, but the contribution from the AGB progenitors is negligible. The post-merger envelope mixing must eventually cease ~ 10^6yr after the dynamic merger phase, before the star enters the RCB phase.

preprint2012arXiv

Do R Coronae Borealis Stars Form from Double White Dwarf Mergers?

A leading formation scenario for R Coronae Borealis (RCB) stars invokes the merger of degenerate He and CO white dwarfs (WD) in a binary. The observed ratio of 16O/18O for RCB stars is in the range of 0.3-20 much smaller than the solar value of ~500. In this paper, we investigate whether such a low ratio can be obtained in simulations of the merger of a CO and a He white dwarf. We present the results of five 3-dimensional hydrodynamic simulations of the merger of a double white dwarf system where the total mass is 0.9 Mdot and the initial mass ratio (q) varies between 0.5 and 0.99. We identify in simulations with $q\lesssim0.7$ a feature around the merged stars where the temperatures and densities are suitable for forming 18O. However, more 16O is being dredged-up from the C- and O-rich accretor during the merger than the amount of 18O that is produced. Therefore, on a dynamical time scale over which our hydrodynamics simulation runs, a 16O/18O ratio of ~2000 in the "best" case is found. If the conditions found in the hydrodynamic simulations persist for 10^6 seconds the oxygen ratio drops to 16 in one case studied, while in a hundred years it drops to ~4 in another case studied, consistent with the observed values in RCB stars. Therefore, the merger of two white dwarfs remains a strong candidate for the formation of these enigmatic stars.

preprint2012arXiv

Nova Framework: A New Tool For Modeling of Nova Outbursts and Nucleosynthesis

Classical novae are the results of surface thermonuclear explosions of H-rich material accreted by white dwarfs (WDs) from their low-mass main-sequence or red-giant binary companions. Chemical composition analysis of their ejecta shows that nova outbursts occur on both carbon-oxygen (CO) and oxygen-neon (ONe) WDs, and that there is cross-boundary mixing between the accreted envelope and underlying WD. We have combined the stellar evolution code MESA and post-processing nucleosynthesis tools of NuGrid into a framework that allows to produce up-to-date models of nova outbursts and compute detailed nucleosynthesis in novae occurring on CO and ONe WDs. The convective boundary mixing (CBM) in our 1D numerical simulations is implemented using a diffusion coefficient that is exponentially decreasing with a distance below the bottom of the convective envelope. This MESA CBM prescription is based on the findings in 3D hydrodynamic simulations that the velocity field, and along with it the mixing expressed in terms of a diffusion coefficient, decays exponentially in the stable layer adjacent to a convective boundary. The framework can also use the commonly adopted 1D nova model in which the CBM is mimicked by assuming that the accreted envelope has been pre-enriched with WD's material. In this preliminary report, we present the most interesting new results related to CO and ONe nova outbursts that have been obtained with the Nova Framework.

preprint2011arXiv

Convective-reactive proton-C12 combustion in Sakurai's object (V4334 Sagittarii) and implications for the evolution and yields from the first generations of stars

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(α, n)O16 reaction.

preprint2010arXiv

Fluorine and Sodium in C-rich low-metallicity stars

We present the N, O, F and Na abundance and 12C/13C isotopic ratio measurements or upper limits for a sample of 10 C-rich, metal-poor giant stars, eight enhanced in s-process (CEMP-s) elements and two poor in n-capture elements (CEMP-no). The abundances are derived from IR, K-band, high-resolution CRIRES@VLT spectra obtained. The metallicity of our sample ranges from [Fe/H]=-3.4 to -1.3. F abundance could be measured only in two CEMP-s stars. With [F/Fe]=0.64, one is mildly F-overabundant, while the other is F-rich, at [F/Fe]=1.44. For the remaining eight objects, including both CEMP-no in our sample, only upper limits on F abundance could be placed. Our measurements and upper limits show that there is a spread in [F/C+N] ratio in CEMP-s stars as predicted by theory. Predictions from nucleosynthetic models for low-mass, low-metallicity Asymptotic Giant Branch stars, account for the derived F abundances, while the upper limits on F content derived for most of the stars are lower than the predicted values. The measured Na content is accounted for by AGB models in the 1.25 to 1.7 Msun range, confirming that the stars responsible for the peculiar abundance pattern observed in CEMP-s stars are low-mass, low-metallicity AGB stars, in agreement with the most accepted astrophysical scenario. We conclude that the mechanism of F production in current state-of-the-art low-metallicity low-mass AGB models needs further scrutiny and that F measurements in a larger number of metal-poor stars are needed to better constraint the models.

preprint2010arXiv

The effect of 12C + 12C rate uncertainties on the weak s-process component

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.

preprint2008arXiv

Evolution and chemical and dynamical effects of high-mass stars

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.