-
Nova contributions to the chemical evolution of the Milky Way
Authors:
Alex J. Kemp,
Amanda I. Karakas,
Andrew R. Casey,
Benoit Cote,
Robert G. Izzard,
Zara Osborn
Abstract:
Context. The explosive burning that drives nova eruptions results in unique nucleosynthesis that heavily over-produces certain isotopes relative to the solar abundance. However, novae are often ignored when considering the chemical evolution of our Galaxy due to their low ejecta masses. Aims. In this work, we use previously computed synthetic nova populations and the galactic chemical evolution co…
▽ More
Context. The explosive burning that drives nova eruptions results in unique nucleosynthesis that heavily over-produces certain isotopes relative to the solar abundance. However, novae are often ignored when considering the chemical evolution of our Galaxy due to their low ejecta masses. Aims. In this work, we use previously computed synthetic nova populations and the galactic chemical evolution code OMEGA+ to assess the impact that novae have on the evolution of stable elemental and isotopic abundances. Methods. We combine populations of novae computed using the binary population synthesis code binary_c with the galactic chemical evolution code OMEGA+ and detailed, white dwarf mass-dependent nova yields to model the nucleosynthetic contributions of novae to the evolution of the Milky Way. We consider three different nova yield profiles, each corresponding to a different set of nova yield calculations. Results. Despite novae from low-mass white dwarfs (WDs) dominating nova ejecta contributions, we find that novae occurring on massive WDs are still able to contribute significantly to many isotopes, particularly those with high mass numbers. We find that novae can produce up to 35% of the Galactic 13C and 15N mass by the time the model Galaxy reaches [Fe/H] = 0, and earlier in the evolution of the Galaxy (between [Fe/H] = -2 and -1) novae may have been the dominant source of 15N. Predictions for [13C/Fe], [15N/Fe], 12C/13C, and 14N/15N abundances ratios vary by up to 0.2 dex at [Fe/H] = 0 and by up to 0.7 dex in [15N/Fe] and 14N/15N between [Fe/H] = -2 and -1 (corresponding approximately to Galactic ages of 170 Myr and 1 Gyr in our model). The Galactic evolution of other stable isotopes (excluding Li) is not noticeably affected by including novae.
△ Less
Submitted 26 July, 2024;
originally announced July 2024.
-
Observational constraints on the origin of the elements. VI. Origin and evolution of neutron-capture elements as probed by the Gaia-ESO survey
Authors:
Jianhui Lian,
Nicholas Storm,
Guillaume Guiglion,
Aldo Serenelli,
Benoit Cote,
Amanda I. Karakas,
Nick Boardman,
Maria Bergemann
Abstract:
Most heavy elements beyond the iron peak are synthesized via neutron capture processes. The nature of the astrophysical sites of neutron capture processes is still very unclear. In this work we explore the observational constraints of the chemical abundances of s-process and r-process elements on the sites of neutron-capture processes by applying Galactic chemical evolution (GCE) models to the dat…
▽ More
Most heavy elements beyond the iron peak are synthesized via neutron capture processes. The nature of the astrophysical sites of neutron capture processes is still very unclear. In this work we explore the observational constraints of the chemical abundances of s-process and r-process elements on the sites of neutron-capture processes by applying Galactic chemical evolution (GCE) models to the data from Gaia-ESO large spectroscopic stellar survey. For the r-process, the [Eu/Fe]-[Fe/H] distribution suggests a short delay time of the site that produces Eu. Other independent observations (e.g., NS-NS binaries), however, suggest a significant fraction of long delayed ($>1$Gyr) neutron star mergers (NSM). When assuming NSM as the only r-process sites, these two observational constraints are inconsistent at above 1$σ$ level. Including short delayed r-process sites like magneto-rotational supernova can resolve this inconsistency. For the s-process, we find a weak metallicity dependence of the [Ba/Y] ratio, which traces the s-process efficiency. Our GCE model with up-to-date yields of AGB stars qualitatively reproduces this metallicity dependence, but the model predicts a much higher [Ba/Y] ratio compared to the data. This mismatch suggests that the s-process efficiency of low mass AGB stars in the current AGB nucleosynthesis models could be overestimated.
△ Less
Submitted 21 September, 2023; v1 submitted 2 August, 2023;
originally announced August 2023.
-
The chemical evolution of the solar neighbourhood for planet-hosting stars
Authors:
Marco Pignatari,
Thomas C. L. Trueman,
Kate A. Womack,
Brad K. Gibson,
Benoit Côté,
Diego Turrini,
Christopher Sneden,
Stephen J. Mojzsis,
Richard J. Stancliffe,
Paul Fong,
Thomas V. Lawson,
James D. Keegans,
Kate Pilkington,
Jean-Claude Passy,
Timothy C. Beers,
Maria Lugaro
Abstract:
Theoretical physical-chemical models for the formation of planetary systems depend on data quality for the Sun's composition, that of stars in the solar neighbourhood, and of the estimated "pristine" compositions for stellar systems. The effective scatter and the observational uncertainties of elements within a few hundred parsecs from the Sun, even for the most abundant metals like carbon, oxygen…
▽ More
Theoretical physical-chemical models for the formation of planetary systems depend on data quality for the Sun's composition, that of stars in the solar neighbourhood, and of the estimated "pristine" compositions for stellar systems. The effective scatter and the observational uncertainties of elements within a few hundred parsecs from the Sun, even for the most abundant metals like carbon, oxygen and silicon, are still controversial. Here we analyse the stellar production and the chemical evolution of key elements that underpin the formation of rocky (C, O, Mg, Si) and gas/ice giant planets (C, N, O, S). We calculate 198 galactic chemical evolution (GCE) models of the solar neighbourhood to analyse the impact of different sets of stellar yields, of the upper mass limit for massive stars contributing to GCE ($M_{\rm up}$) and of supernovae from massive-star progenitors which do not eject the bulk of the iron-peak elements (faint supernovae). Even considering the GCE variation produced via different sets of stellar yields, the observed dispersion of elements reported for stars in the Milky Way disk is not reproduced. Among others, the observed range of super-solar [Mg/Si] ratios, sub-solar [S/N], and the dispersion of up to 0.5 dex for [S/Si] challenge our models. The impact of varying $M_{\rm up}$ depends on the adopted supernova yields. Thus, observations do not provide a constraint on the M$_{\rm up}$ parametrization. When including the impact of faint supernova models in GCE calculations, elemental ratios vary by up to 0.1-0.2 dex in the Milky Way disk; this modification better reproduces observations.
△ Less
Submitted 23 July, 2023;
originally announced July 2023.
-
Inhomogeneous enrichment of radioactive nuclei in the Galaxy: Deposition of live Mn-53, Fe-60, Hf-182, and Pu-244 into deep-sea archives. Surfing the wave?
Authors:
Benjamin Wehmeyer,
Andrés Yagüe López,
Benoit Côté,
Maria K. Pető,
Chiaki Kobayashi,
Maria Lugaro
Abstract:
While modelling the galactic chemical evolution (GCE) of stable elements provides insights to the formation history of the Galaxy and the relative contributions of nucleosynthesis sites, modelling the evolution of short-lived radioisotopes (SLRs) can provide supplementary timing information on recent nucleosynthesis. To study the evolution of SLRs, we need to understand their spatial distribution.…
▽ More
While modelling the galactic chemical evolution (GCE) of stable elements provides insights to the formation history of the Galaxy and the relative contributions of nucleosynthesis sites, modelling the evolution of short-lived radioisotopes (SLRs) can provide supplementary timing information on recent nucleosynthesis. To study the evolution of SLRs, we need to understand their spatial distribution. Using a 3-dimensional GCE model, we investigated the evolution of four SLRs: Mn-53, Fe-60, Hf-182, and Pu-244 with the aim of explaining detections of recent (within the last $\approx$1-20 Myr) deposition of live Mn-53, Fe-60, and Pu-244 of extrasolar origin into deep-sea reservoirs. We find that core-collapse supernovae (CCSNe) are the dominant propagation mechanism of SLRs in the Galaxy. This results in the simultaneously arrival of these four SLRs on Earth, although they could have been produced in different astrophysical sites, which can explain why live extrasolar Mn-53, Fe-60, and Pu-244 are found within the same, or similar, layers of deep-sea sediments. We predict that Hf-182 should also be found in such sediments at similar depths.
△ Less
Submitted 21 February, 2023; v1 submitted 11 January, 2023;
originally announced January 2023.
-
Chemical Evolution of Fluorine in the Milky Way
Authors:
Kate A. Womack,
Fiorenzo Vincenzo,
Brad K. Gibson,
Benoit Côté,
Marco Pignatari,
Hannah E. Brinkman,
Paolo Ventura,
Amanda Karakas
Abstract:
Fluorine has many different potential sites and channels of production, making narrowing down a dominant site of fluorine production particularly challenging. In this work, we investigate which sources are the dominant contributors to the galactic fluorine by comparing chemical evolution models to observations of fluorine abundances in Milky Way stars covering a metallicity range -2$<$[Fe/H]$<$0.4…
▽ More
Fluorine has many different potential sites and channels of production, making narrowing down a dominant site of fluorine production particularly challenging. In this work, we investigate which sources are the dominant contributors to the galactic fluorine by comparing chemical evolution models to observations of fluorine abundances in Milky Way stars covering a metallicity range -2$<$[Fe/H]$<$0.4 and upper limits in the range -3.4$<$[Fe/H]$<$-2.3. In our models, we use a variety of stellar yield sets in order to explore the impact of varying both AGB and massive star yields on the chemical evolution of fluorine. In particular, we investigate different prescriptions for initial rotational velocity in massive stars as well as a metallicity dependent mix of rotational velocities. We find that the observed [F/O] and [F/Fe] abundance ratios at low metallicity and the increasing trend of [F/Ba] at [Fe/H]$\gtrsim$-1 can only be reproduced by chemical evolution models assuming, at all metallicities, a contribution from rapidly rotating massive stars with initial rotational velocities as high as 300km s$^{-1}$. A mix of rotational velocities may provide a more physical solution than the sole use of massive stars with $v_{\text{rot}}$=300$\text{km s}^{-1}$, which are predicted to overestimate the fluorine and average s-process elemental abundances at [Fe/H]$\gtrsim$-1. The contribution from AGB stars is predicted to start at [Fe/H]$\approx$-1 and becomes increasingly important at high metallicity, being strictly coupled to the evolution of the nitrogen abundance. Finally, by using modern yield sets, we investigate the fluorine abundances of Wolf-Rayet winds, ruling them out as dominant contributors to the galactic fluorine.
△ Less
Submitted 1 November, 2022;
originally announced November 2022.
-
Enrichment of the Galactic disc with neutron-capture elements: Gd, Dy, and Th
Authors:
T. Mishenina,
M. Pignatari,
T. Gorbaneva,
B. Cote,
A. Yague Lopez,
F. -K. Thielemann,
C. Soubiran
Abstract:
The study of the origin of heavy elements is one of the main goals of nuclear astrophysics. In this paper, we present new observational data for the heavy $r$-process elements gadolinium (Gd, Z=64), dysprosium (Dy, Z=66) and thorium (Th, Z=90) in a sample of 276 Galactic disc stars ( --1.0$<$[Fe/H]$<$+0.3). The stellar spectra have a high resolution of 42,000 and 75,000, and the signal-to-noise ra…
▽ More
The study of the origin of heavy elements is one of the main goals of nuclear astrophysics. In this paper, we present new observational data for the heavy $r$-process elements gadolinium (Gd, Z=64), dysprosium (Dy, Z=66) and thorium (Th, Z=90) in a sample of 276 Galactic disc stars ( --1.0$<$[Fe/H]$<$+0.3). The stellar spectra have a high resolution of 42,000 and 75,000, and the signal-to-noise ratio higher than 100. The LTE abundances of Gd, Dy and Th have been determined by comparing the observed and synthetic spectra for three Gd lines (149 stars), four Dy lines (152 stars) and the Th line at 4019.13 A (170 stars). For about 70% of the stars in our sample Gd and Dy are measured for the first time, and Th for 95% of the stars. Typical errors vary from 0.07 to 0.16 dex. This paper provides the first extended set of Th observations in the Milky Way disc. Together with europium (Eu, Z = 63) data from our previous studies, we have compared these new observations with nucleosynthesis predictions and Galactic Chemical Evolution simulations. We confirm that [Gd/Fe] and [Dy/Fe] show the same behavior of Eu. We study with GCE simulations the evolution of [Th/Fe] in comparison with [Eu/Fe], showing that unlike Eu either the Th production is metallicity dependent in case of a unique source of the r-process in the Galaxy, or the frequency of the Th-rich r-process source is decreasing with the increasing of [Fe/H].
△ Less
Submitted 24 August, 2022;
originally announced August 2022.
-
Origin of Plutonium-244 in the Early Solar System
Authors:
Maria Lugaro,
Andrés Yagüe López,
Benjámin Soós,
Benoit Côté,
Mária Pető,
Nicole Vassh,
Benjamin Wehmeyer,
Marco Pignatari
Abstract:
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-pr…
▽ More
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.
△ Less
Submitted 3 August, 2022;
originally announced August 2022.
-
Viability of novae as sources of Galactic lithium
Authors:
Alex J. Kemp,
Amanda I. Karakas,
Andrew R. Casey,
Benoit Cote,
Robert G. Izzard,
Zara Osborn
Abstract:
Of all the light elements, the evolution of lithium (Li) in the Milky Way is perhaps the most difficult to explain. Li is difficult to synthesize and is easily destroyed, making most stellar sites unsuitable for producing Li in sufficient quantities to account for the proto-solar abundance. For decades, novae have been proposed as a potential explanation to this 'Galactic Li problem', and the rece…
▽ More
Of all the light elements, the evolution of lithium (Li) in the Milky Way is perhaps the most difficult to explain. Li is difficult to synthesize and is easily destroyed, making most stellar sites unsuitable for producing Li in sufficient quantities to account for the proto-solar abundance. For decades, novae have been proposed as a potential explanation to this 'Galactic Li problem', and the recent detection of 7Be in the ejecta of multiple nova eruptions has breathed new life into this theory. In this work, we assess the viability of novae as dominant producers of Li in the Milky Way. We present the most comprehensive treatment of novae in a galactic chemical evolution code to date, testing theoretical- and observationally-derived nova Li yields by integrating metallicity-dependent nova ejecta profiles computed using the binary population synthesis code binary c with the galactic chemical evolution code OMEGA+. We find that our galactic chemical evolution models which use observationally-derived Li yields account for the proto-solar Li abundance very well, while models relying on theoretical nova yields cannot reproduce the proto-solar observation. A brief exploration of physical uncertainties including single-stellar yields, the metallicity resolution of our nova treatment, common-envelope physics, and nova accretion efficiencies indicates that this result is robust to physical assumptions. Scatter within the observationally-derived Li yields in novae is identified as the primary source of uncertainty, motivating further observations of 7Be in nova ejecta.
△ Less
Submitted 27 June, 2022;
originally announced June 2022.
-
Observational constraints on the origin of the elements. V. Non-LTE abundance ratios of [Ni/Fe] in Galactic stars and enrichment by sub-Chandrasekhar mass SNe
Authors:
P. Eitner,
M. Bergemann,
A. J. Ruiter,
O. Avril,
I. R. Seitenzahl,
M. R. Gent,
B. Côté
Abstract:
We constrain the role of different SN Ia channels in the chemical enrichment of the Galaxy by studying the abundances of nickel in Galactic stars. We investigate four different SN Ia sub-classes, including the classical single-degenerate near-Chandrasekhar mass SN Ia, the fainter SN Iax systems associated with He accretion from the companion, as well as two sub-Ch mass SN Ia channels. The latter i…
▽ More
We constrain the role of different SN Ia channels in the chemical enrichment of the Galaxy by studying the abundances of nickel in Galactic stars. We investigate four different SN Ia sub-classes, including the classical single-degenerate near-Chandrasekhar mass SN Ia, the fainter SN Iax systems associated with He accretion from the companion, as well as two sub-Ch mass SN Ia channels. The latter include the double-detonation of a white dwarf accreting helium-rich matter and violent white dwarf mergers. NLTE models of Fe and Ni are used in the abundance analysis. In the GCE models, we include new delay time distributions arising from the different SN Ia channels, as well as recent yields for core-collapse supernovae and AGB stars. The data-model comparison is performed using a Markov chain Monte Carlo framework that allows us to explore the entire parameter space allowed by the diversity of explosion mechanisms and the Galactic SN Ia rate, taking into account the uncertainties of the observed data. We show that NLTE effects have a non-negligible impact on the observed [Ni/Fe] ratios in the Galactic stars. The NLTE corrections to Ni abundances are not large, but strictly positive, lifting the [Ni/Fe] ratios by ~+0.15 dex at [Fe/H] =-2. We find that that the distributions of [Ni/Fe] in LTE and in NLTE are very tight, with a scatter of < 0.1 dex at all metallicities, supporting earlier work. In LTE, most stars have scaled-solar Ni abundances, [Ni/Fe] = 0, with a slight tendency for sub-solar [Ni/Fe] ratios at lower [Fe/H]. In NLTE, however, we find a mild anti-correlation between [Ni/Fe] and metallicity, and a slightly elevated [Ni/Fe] ratios at [Fe/H] < -1.0. The NLTE data can be explained by the GCE models calculated with a substantial, ~ 75%, fraction of sub-Ch SN Ia.
△ Less
Submitted 14 July, 2023; v1 submitted 21 June, 2022;
originally announced June 2022.
-
Horizons: Nuclear Astrophysics in the 2020s and Beyond
Authors:
H. Schatz,
A. D. Becerril Reyes,
A. Best,
E. F. Brown,
K. Chatziioannou,
K. A. Chipps,
C. M. Deibel,
R. Ezzeddine,
D. K. Galloway,
C. J. Hansen,
F. Herwig,
A. P. Ji,
M. Lugaro,
Z. Meisel,
D. Norman,
J. S. Read,
L. F. Roberts,
A. Spyrou,
I. Tews,
F. X. Timmes,
C. Travaglio,
N. Vassh,
C. Abia,
P. Adsley,
S. Agarwal
, et al. (140 additional authors not shown)
Abstract:
Nuclear Astrophysics is a field at the intersection of nuclear physics and astrophysics, which seeks to understand the nuclear engines of astronomical objects and the origin of the chemical elements. This white paper summarizes progress and status of the field, the new open questions that have emerged, and the tremendous scientific opportunities that have opened up with major advances in capabilit…
▽ More
Nuclear Astrophysics is a field at the intersection of nuclear physics and astrophysics, which seeks to understand the nuclear engines of astronomical objects and the origin of the chemical elements. This white paper summarizes progress and status of the field, the new open questions that have emerged, and the tremendous scientific opportunities that have opened up with major advances in capabilities across an ever growing number of disciplines and subfields that need to be integrated. We take a holistic view of the field discussing the unique challenges and opportunities in nuclear astrophysics in regards to science, diversity, education, and the interdisciplinarity and breadth of the field. Clearly nuclear astrophysics is a dynamic field with a bright future that is entering a new era of discovery opportunities.
△ Less
Submitted 16 May, 2022;
originally announced May 2022.
-
The RADIOSTAR Project
Authors:
Maria Lugaro,
Benoit Côté,
Marco Pignatari,
Andrés Yagüe López,
Hannah Brinkman,
Borbála Cseh,
Jacqueline Den Hartogh,
Carolyn Louise Doherty,
Amanda Irene Karakas,
Chiaki Kobayashi,
Thomas Lawson,
Mária Pető,
Benjámin Soós,
Thomas Trueman,
Blanka Világos
Abstract:
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,…
▽ More
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.
△ Less
Submitted 16 February, 2022;
originally announced February 2022.
-
Galactic Chemical Evolution of Radioactive Isotopes with an s-process Contribution
Authors:
T. C. L. Trueman,
B. Côté,
A. Yagüe López,
J. den Hartogh,
M. Pignatari,
B. Soós,
A. I. Karakas,
M. Lugaro
Abstract:
Analysis of inclusions in primitive meteorites reveals that several short-lived radionuclides (SLRs) with half-lives $0.1-100$ Myr existed in the early Solar System (ESS). We investigate the ESS origin of $^{107}$Pd, $^{135}$Cs, and $^{182}$Hf, which are produced by $slow$ neutron captures (the $s$-process) in asymptotic giant branch (AGB) stars. We modelled the galactic abundances of these SLRs u…
▽ More
Analysis of inclusions in primitive meteorites reveals that several short-lived radionuclides (SLRs) with half-lives $0.1-100$ Myr existed in the early Solar System (ESS). We investigate the ESS origin of $^{107}$Pd, $^{135}$Cs, and $^{182}$Hf, which are produced by $slow$ neutron captures (the $s$-process) in asymptotic giant branch (AGB) stars. We modelled the galactic abundances of these SLRs using the OMEGA+ galactic chemical evolution (GCE) code and two sets of mass- and metallicity-dependent AGB nucleosynthesis yields (Monash and FRUITY). Depending on the ratio of the mean life $τ$ of the SLR to the average length of time between the formation of AGB progenitor $γ$, we calculate timescales relevant for the birth of the Sun. If $τ/γ\gtrsim2$, we predict self-consistent isolation times between 9 and 26 Myr by decaying the GCE predicted $^{107}$Pd/$^{108}$Pd, $^{135}$Cs/$^{133}$Cs, and $^{182}$Hf/$^{180}$Hf ratios to their respective ESS ratios. The predicted $^{107}$Pd/$^{182}$Hf ratio indicates that our GCE models are missing $9-73\%$ of $^{107}$Pd and $^{108}$Pd in the ESS. This missing component may have come from AGB stars of higher metallicity than those that contributed to the ESS in our GCE code. If $τ/γ\lesssim0.3$, we calculate instead the time ($T_{\rm LE}$) from the last nucleosynthesis event that added the SLRs into the presolar matter to the formation of the oldest solids in the ESS. For the 2 M$_\odot$, $Z=0.01$ Monash model we find a self-consistent solution of $T_{\rm LE}=25.5$ Myr.
△ Less
Submitted 21 October, 2021;
originally announced October 2021.
-
Monte Carlo Investigation of the Ratios of Short-Lived Radioactive Isotopes in the Interstellar Medium
Authors:
A. Yagüe López,
B. Côté,
M. Lugaro
Abstract:
Short-lived radioactive nuclei (SLR) with mean lives below 100 Myr provide us with unique insights into current galactic nucleosynthetic events, as well as events that contributed to the material of our Solar System more that 4.6 Gyr ago. Here we present a statistical analysis of the ratios of these radioactive nuclei at the time of early Solar System (ESS) using both analytical derivations and Mo…
▽ More
Short-lived radioactive nuclei (SLR) with mean lives below 100 Myr provide us with unique insights into current galactic nucleosynthetic events, as well as events that contributed to the material of our Solar System more that 4.6 Gyr ago. Here we present a statistical analysis of the ratios of these radioactive nuclei at the time of early Solar System (ESS) using both analytical derivations and Monte Carlo methods. We aim to understand the interplay between the production frequency and the mean lives of these isotopes, and its impact on their theoretically predicted ratios in the interstellar medium (ISM). We find that when the ratio of two SRLs, instead of the ratios of each single SLR relative to its stable or long-lived isotope, is considered, not only the uncertainties related to the galactic chemical evolution of the stable isotope are completely eliminated, but also the statistical uncertainties are much lower. We identify four ratios, 247Cm/129I, 107Pd/182Hf, 97Tc/98Tc, and 53Mn/97Tc, that have the potential to provide us with new insights into the r-, s-, and p-process nucleosynthesis at the time of the formation of the Sun, and need to be studied using variable stellar yields. Additionally, the latter two ratios need to be better determined in the ESS to allow us to fully exploit them to investigate the galactic sites of the p process.
△ Less
Submitted 19 May, 2021;
originally announced May 2021.
-
Metallicity-dependent nucleosynthetic yields of Type Ia supernovae originating from double detonations of sub-M$_{\text{Ch}}$ white dwarfs
Authors:
Sabrina Gronow,
Benoit Cote,
Florian Lach,
Ivo R. Seitenzahl,
Christine E. Collins,
Stuart A. Sim,
Friedrich K. Roepke
Abstract:
Double detonations in sub-Chandrasekhar mass carbon-oxygen white dwarfs with helium shell are a potential explosion mechanism for a Type Ia supernova (SNe Ia). It comprises a shell detonation and subsequent core detonation. The focus of our study is on the effect of the progenitor metallicity on the nucleosynthetic yields. For this, we compute and analyse a set of eleven different models with vary…
▽ More
Double detonations in sub-Chandrasekhar mass carbon-oxygen white dwarfs with helium shell are a potential explosion mechanism for a Type Ia supernova (SNe Ia). It comprises a shell detonation and subsequent core detonation. The focus of our study is on the effect of the progenitor metallicity on the nucleosynthetic yields. For this, we compute and analyse a set of eleven different models with varying core and shell masses at four different metallicities each. This results in a total of 44 models at metallicities between 0.01$Z_\odot$ and 3$Z_\odot$. Our models show a strong impact of the metallicity in the high density regime. The presence of $^{22}$Ne causes a neutron-excess which shifts the production from $^{56}$Ni to stable isotopes such as $^{54}$Fe and $^{58}$Ni in the $α$-rich freeze-out regime. The isotopes of the metallicity implementation further serve as seed nuclei for additional reactions in the shell detonation. Most significantly, the production of $^{55}$Mn increases with metallicity confirming the results of previous work. A comparison of elemental ratios relative to iron shows a relatively good match to solar values for some models. Super-solar values are reached for Mn at 3$Z_\odot$ and solar values in some models at $Z_\odot$. This indicates that the required contribution of SNe Ia originating from Chandrasekhar mass WDs can be lower than estimated in orevious work to reach solar values of [Mn/Fe] at [Fe/H]$=0$. Our galactic chemical evolution models suggest that SNe Ia from sub-Chandrasekhar mass white dwarfs, along with core-collapse supernovae, could account for more than 80% of the solar Mn abundance. Using metallicity-dependent SN Ia yields helps to reproduce the upward trend of [Mn/Fe] as a function of metallicity for the solar neighborhood. These chemical evolution predictions, however, depend on the massive star yields adopted in the calculations.
△ Less
Submitted 1 September, 2021; v1 submitted 25 March, 2021;
originally announced March 2021.
-
Nucleosynthesis imprints from different Type Ia Supernova explosion scenarios and implications for galactic chemical evolution
Authors:
F. Lach,
F. K. Roepke,
I. R. Seitenzahl,
B. Coté,
S. Gronow,
A. J. Ruiter
Abstract:
We analyze the nucleosynthesis yields of various Type Ia supernova explosion simulations including pure detonations in sub- Chandrasekhar mass white dwarfs, double detonations and pure helium detonations of sub-Chandrasekhar mass white dwarfs with an accreted helium envelope, a violent merger model of two white dwarfs and deflagrations as well as delayed detonations in Chandrasekhar mass white dwa…
▽ More
We analyze the nucleosynthesis yields of various Type Ia supernova explosion simulations including pure detonations in sub- Chandrasekhar mass white dwarfs, double detonations and pure helium detonations of sub-Chandrasekhar mass white dwarfs with an accreted helium envelope, a violent merger model of two white dwarfs and deflagrations as well as delayed detonations in Chandrasekhar mass white dwarfs. We focus on the iron peak elements Mn, Zn and Cu. To this end, we also briefly review the different burning regimes and production sites of these elements as well as the results of abundance measurements and several galactic chemical evolution studies. We find that super-solar values of [Mn/Fe] are not restricted to Chandrasekhar mass explosion models. Scenarios including a helium detonation can significantly contribute to the production of Mn, in particular the models proposed for calcium-rich transients. Although Type Ia supernovae are often not accounted for as production sites of Zn and Cu, our models involving helium shell detonations can produce these elements in super-solar ratios relative to Fe. Our results suggest a re-consideration of Type Ia supernova yields in galactic chemical evolution models. A detailed comparison with observations can provide new insight into the progenitor and explosion channels of these events.
△ Less
Submitted 27 October, 2020;
originally announced October 2020.
-
The gas, metal and dust evolution in low-metallicity local and high-redshift galaxies
Authors:
Ambra Nanni,
Denis Burgarella,
Patrice Theulé,
Benoit Côté,
Hiroyuki Hirashita
Abstract:
The chemical enrichment in the interstellar medium (ISM) of galaxies is regulated by several physical processes: stellar evolution, grain formation and destruction, galactic inflows and outflows. Understanding such processes is essential to follow the chemical enrichment of galaxies through the cosmic epochs, and to interpret the observations. Despite the importance of such topics, the efficiency…
▽ More
The chemical enrichment in the interstellar medium (ISM) of galaxies is regulated by several physical processes: stellar evolution, grain formation and destruction, galactic inflows and outflows. Understanding such processes is essential to follow the chemical enrichment of galaxies through the cosmic epochs, and to interpret the observations. Despite the importance of such topics, the efficiency of the different processes driving the evolution of baryons in galaxies, remain controversial. We revise the current description of metal and dust evolution in local low-metallicity dwarf galaxies and we develop a description for Lyman Break Galaxies. Our main goal is to reproduce i) the peak in the mass of dust over the mass of stars (sMdust) observed within few hundred Myrs; ii) the decrease of the sMdust at later time. The spectral energy distribution of the galaxies is fitted with the "Code Investigating GALaxies Emission" (CIGALE), through which the stellar and dust masses, and the star formation rate are estimated. For some of the dwarf galaxies, the metal and gas content are also available. We run different calculations of chemical evolution in galaxies, and we fit the observed properties through the model predictions. We show that i) a top-heavy initial mass function that favours massive stars and a dust condensation fraction for Type II Supernovae (SNe II) of 50% or more help to reproduce the peak of sMdust observed after 100 Myrs since the beginning of the cycle; ii) galactic outflows play a crucial role in reproducing the decline in sMdust with age, and they are more efficient than grain destruction from SNe II; iii) a star formation efficiency (mass of gas converted into stars) of few per cent is required to explain the metallicity of local dwarf galaxies; iv) dust growth in the ISM is not necessary to reproduce the sMdust and, if present, its effect is erased by galactic outflows.
△ Less
Submitted 26 June, 2020;
originally announced June 2020.
-
129I and 247Cm in Meteorites Constrain the Last Astrophysical Source of Solar r-process Elements
Authors:
Benoit Côté,
Marius Eichler,
Andrés Yagüe,
Nicole Vassh,
Matthew R. Mumpower,
Blanka Világos,
Benjámin Soós,
Almudena Arcones,
Trevor M. Sprouse,
Rebecca Surman,
Marco Pignatari,
Maria K. Pető,
Benjamin Wehmeyer,
Thomas Rauscher,
Maria Lugaro
Abstract:
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…
▽ More
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.
△ Less
Submitted 2 March, 2021; v1 submitted 8 June, 2020;
originally announced June 2020.
-
Stochastic Chemical Evolution of Radioactive Isotopes with a Monte Carlo Approach
Authors:
Benoit Côté,
Andrés Yagüe,
Blanka Világos,
Maria Lugaro
Abstract:
Short-lived radionuclides (SLRs) with mean-lives $τ$ of a few to hundreds Myr provide unique opportunities to probe recent nucleosynthesis events in the interstellar medium, and the physical conditions in which the Sun formed. Here we quantify the uncertainty in the predicted evolution of SLRs within a parcel of interstellar gas given the stochastic nature of stellar enrichment events. We assume t…
▽ More
Short-lived radionuclides (SLRs) with mean-lives $τ$ of a few to hundreds Myr provide unique opportunities to probe recent nucleosynthesis events in the interstellar medium, and the physical conditions in which the Sun formed. Here we quantify the uncertainty in the predicted evolution of SLRs within a parcel of interstellar gas given the stochastic nature of stellar enrichment events. We assume that an enrichment progenitor is formed at every time interval $γ$. For each progenitor, we randomly sample the delay time between its formation and its enrichment event, based on several delay-time distribution (DTD) functions that cover a wide range of astrophysical sites. For each set of $τ$, $γ$, and DTD function, we follow the abundances of SLRs for 15 Gyr, and repeat this process thousands of times to derive their probability distributions. For $τ/γ\gtrsim2$, the distributions depend on the DTD function and we provide tabulated values and analytical expressions to quantify the spread. The relative abundance uncertainty reaches a maximum of $\sim$ 60% for $τ/γ=1$. For $τ/γ\lesssim1$, we provide the probability for the SLR abundance to carry the signature of only one enrichment event, which is greater than 50% when $τ/γ\lesssim0.3$. For $0.3\lesssimτ/γ\lesssim 2$, a small number of events contributed to the SLR abundance. This case needs to be investigated with a separate statistical method. We find that an isolation time for the birth of the Sun of roughly $9-13$ Myr is consistent with the observed abundances of $^{60}$Fe, $^{107}$Pd, and $^{182}$Hf in the early Solar System, when assuming $τ/γ\sim3$ for these isotopes.
△ Less
Submitted 13 November, 2019; v1 submitted 4 November, 2019;
originally announced November 2019.
-
The Impact of Nuclear Physics Uncertainties on Galactic Chemical Evolution Predictions
Authors:
Benoit Côté,
Pavel Denissenkov,
Falk Herwig,
Chris L. Fryer,
Krzysztof Belczynski,
Nicole Vassh,
Matthew R. Mumpower,
Jonas Lippuner,
Marco Pignatari,
Ashley J. Ruiter
Abstract:
Modeling the evolution of the elements in the Milky Way is a multidisciplinary and challenging task. In addition to simulating the 13 billion years evolution of our Galaxy, chemical evolution simulations must keep track of the elements synthesized and ejected from every astrophysical site of interest (e.g., supernova, compact binary merger). The elemental abundances of such ejecta, which are a fun…
▽ More
Modeling the evolution of the elements in the Milky Way is a multidisciplinary and challenging task. In addition to simulating the 13 billion years evolution of our Galaxy, chemical evolution simulations must keep track of the elements synthesized and ejected from every astrophysical site of interest (e.g., supernova, compact binary merger). The elemental abundances of such ejecta, which are a fundamental input for chemical evolution codes, are usually taken from theoretical nucleosynthesis calculations performed by the nuclear astrophysics community. Therefore, almost all chemical evolution predictions rely on the nuclear physics behind those calculations. In this proceedings, we highlight the impact of nuclear physics uncertainties on galactic chemical evolution predictions. We demonstrate that nuclear physics and galactic evolution uncertainties both have a significant impact on interpreting the origin of neutron-capture elements in our Solar System. Those results serve as a motivation to create and maintain collaborations between the fields of nuclear astrophysics and galaxy evolution.
△ Less
Submitted 31 October, 2019;
originally announced November 2019.
-
LRP2020: The cosmic origin and evolution of the elements
Authors:
Rodrigo Fernández,
Falk Herwig,
Samar Safi-Harb,
Iris Dillmann,
Kim A. Venn,
Benoit Côté,
Craig O. Heinke,
Erik Rosolowsky,
Tyrone E. Woods,
Daryl Haggard,
Luis Lehner,
John J. Ruan,
Daniel M. Siegel,
Jo Bovy,
Alan A. Chen,
Andrew Cumming,
Barry Davids,
Maria R. Drout,
Reiner Krüecken
Abstract:
The origin of many elements of the periodic table remains an unsolved problem. While many nucleosynthetic channels are broadly understood, significant uncertainties remain regarding certain groups of elements such as the intermediate and rapid neutron-capture processes, the p-process, or the origin of odd-Z elements in the most metal-poor stars. Canada has a long tradition of leadership in nuclear…
▽ More
The origin of many elements of the periodic table remains an unsolved problem. While many nucleosynthetic channels are broadly understood, significant uncertainties remain regarding certain groups of elements such as the intermediate and rapid neutron-capture processes, the p-process, or the origin of odd-Z elements in the most metal-poor stars. Canada has a long tradition of leadership in nuclear astrophysics, dating back to the work of Alastair Cameron in the 1950s. Recent faculty hires have further boosted activity in the field, including transient observation and theory, survey science on galactic nucleosynthesis, and nuclear experiments. This white paper contains a brief overview of recent activity in the community, highlighting strengths in each sub-field, and provides recommendations to improve interdisciplinary collaboration. Sustaining Canadian leadership in the next decade will require, on the observational side, access to transient and non-transient surveys like LSST, SKA, or MSE, support for target-of-opportunity observing in current and future Canadian telescopes, and participation in next-generation X-ray telescopes such as ATHENA. State-of-the-art theoretical predictions will require an ambitious succession plan for the Niagara supercomputer to support large parallel jobs. We propose a funding instrument for postdoctoral training that reflects the interdisciplinary nature of nuclear astrophysics research, and the creation of a national collaborative funding program that allows for joint projects and workshop organization.
△ Less
Submitted 21 October, 2019;
originally announced October 2019.
-
Aluminium-26 from massive binary stars I: non-rotating models
Authors:
H. E. Brinkman,
C. L. Doherty,
O. R. Pols,
E. T. Li,
B. Côté,
M. Lugaro
Abstract:
Aluminium-26 is a short-lived radionuclide with a half-life of 0.72Myr, which is observed today in the Galaxy via gamma-ray spectroscopy and is inferred to have been present in the early Solar System via analysis of meteorites. Massive stars are considered the main contributors of Al26. Although most massive stars are found in binary systems, the effect, however, of binary interactions on the Al26…
▽ More
Aluminium-26 is a short-lived radionuclide with a half-life of 0.72Myr, which is observed today in the Galaxy via gamma-ray spectroscopy and is inferred to have been present in the early Solar System via analysis of meteorites. Massive stars are considered the main contributors of Al26. Although most massive stars are found in binary systems, the effect, however, of binary interactions on the Al26 yields have not been investigated since Braun & Langer (1995). Here we aim to fill this gap. We have used the MESA stellar evolution code to compute massive (10Msun<=M<=80Msun), non-rotating, single and binary stars of solar metallicity (Z=0.014). We computed the wind yields for the single stars and for the binary systems where mass transfer plays a major role. Depending on the initial mass of the primary star and orbital period, the Al26 yield can either increase or decrease in a binary system. For binary systems with primary masses up to ~35-40Msun, the yield can increase significantly, especially at the lower mass-end, while above ~45Msun the yield becomes similar to the single star yield or even decreases. Our preliminary results show that compared to supernova explosions, the contribution of mass-loss in binary systems to the total Al26 abundance produced by a stellar population is minor. On the other hand, if massive star mass-loss is the origin of Al26 in the early Solar System, our results will have significant implications for the identification of the potential stellar, or stellar population, source.
△ Less
Submitted 10 September, 2019;
originally announced September 2019.
-
Using failed supernovae to constrain the Galactic r-process element production
Authors:
B. Wehmeyer,
C. Frohlich,
B. Côté,
M. Pignatari,
F. -K. Thielemann
Abstract:
Rapid neutron capture process (r-process) elements have been detected in a large fraction of metal-poor halo stars, with abundances relative to iron (Fe) that vary by over two orders of magnitude. This scatter is reduced to less than a factor of 3 in younger Galactic disc stars. The large scatter of r-process elements in the early Galaxy suggests that the r-process is made by rare events, like com…
▽ More
Rapid neutron capture process (r-process) elements have been detected in a large fraction of metal-poor halo stars, with abundances relative to iron (Fe) that vary by over two orders of magnitude. This scatter is reduced to less than a factor of 3 in younger Galactic disc stars. The large scatter of r-process elements in the early Galaxy suggests that the r-process is made by rare events, like compact binary mergers and rare sub-classes of supernovae. Although being rare, neutron star mergers alone have difficulties to explain the observed enhancement of r-process elements in the lowest metallicity stars compared to Fe. The supernovae producing the two neutron stars already provide a substantial Fe abundance where the r-process ejecta from the merger would be injected. In this work we investigate another complementary scenario, where the r-process occurs in neutron star-black hole mergers in addition to neutron star mergers. Neutron star-black hole mergers would eject similar amounts of r-process matter as neutron star mergers, but only the neutron star progenitor would have produced Fe. Furthermore, a reduced efficiency of Fe production from single stars significantly alters the age-metallicity relation, which shifts the onset of r-process production to lower metallicities. We use the high-resolution [(20 pc)3/cell] inhomogeneous chemical evolution tool `ICE' to study the outcomes of these effects. In our simulations, an adequate combination of neutron star mergers and neutron star-black hole mergers qualitatively reproduces the observed r-process abundances in the Galaxy.
△ Less
Submitted 15 August, 2019;
originally announced August 2019.
-
Enrichment of the Galactic disc with neutron-capture elements: Mo and Ru
Authors:
Tamara Mishenina,
Marco Pignatari,
Tatiana Gorbaneva,
Claudia Travaglio,
Benoit Côté,
Friedrich-Karl Thielemann,
Caroline Soubiran
Abstract:
We present new observational data for the heavy elements molybdenum (Mo, Z = 42) and ruthenium (Ru, Z = 44) in F-, G-, and K-stars belonging to different substructures of the Milky Way. The range of metallicity covered is --1.0 $<$ [Fe/H] $<$ +0.3. The spectra of Galactic disc stars have a high resolution of 42,000 and 75,000 and signal-to-noise ratio better than 100. Mo and Ru abundances were der…
▽ More
We present new observational data for the heavy elements molybdenum (Mo, Z = 42) and ruthenium (Ru, Z = 44) in F-, G-, and K-stars belonging to different substructures of the Milky Way. The range of metallicity covered is --1.0 $<$ [Fe/H] $<$ +0.3. The spectra of Galactic disc stars have a high resolution of 42,000 and 75,000 and signal-to-noise ratio better than 100. Mo and Ru abundances were derived by comparing the observed and synthetic spectra in the region of Mo I lines at 5506, 5533 Å~ for 209 stars and Ru I lines at 4080, 4584, 4757 Å~ for 162 stars using the LTE approach. For all the stars, the Mo and Ru abundance determinations are obtained for the first time with an average error of 0.14 dex. This is the first extended sample of stellar observations for Mo and Ru in the Milky Way disc, and together with earlier observations in halo stars it is pivotal in providing a complete picture of the evolution of Mo and Ru across cosmic timescales.
The Mo and Ru abundances were compared with those of the neutron-capture elements (Sr, Y, Zr, Ba, Sm, Eu). The complex nucleosynthesis history of Mo and Ru is compared with different Galactic Chemical Evolution (GCE) simulations. In general, present theoretical GCE simulations show underproduction of Mo and Ru at all metallicities compared to observations. This highlights a significant contribution of nucleosynthesis processes not yet considered in our simulations. A number of possible scenarios are discussed.
△ Less
Submitted 8 August, 2019;
originally announced August 2019.
-
A New Model for Electron-Capture Supernovae in Galactic Chemical Evolution
Authors:
S. Jones,
B. Côté,
F. K. Roepke,
S. Wanajo
Abstract:
We examine the contribution of electron-capture supernovae (ECSNe), low-mass SNe from collapsing Fe cores (FeCCSNe), and rotating massive stars to the chemical composition of the Galaxy. Our model includes contributions to chemical evolution from both thermonuclear ECSNe (tECSNe) and gravitational collapse ECSNe (cECSNe). We show that if ECSNe are predominantly gravitational collapse SNe but about…
▽ More
We examine the contribution of electron-capture supernovae (ECSNe), low-mass SNe from collapsing Fe cores (FeCCSNe), and rotating massive stars to the chemical composition of the Galaxy. Our model includes contributions to chemical evolution from both thermonuclear ECSNe (tECSNe) and gravitational collapse ECSNe (cECSNe). We show that if ECSNe are predominantly gravitational collapse SNe but about 15% are partial thermonuclear explosions, the model is able to reproduce the solar abundances of several important and problematic isotopes including $^{48}$Ca, $^{50}$Ti and $^{54}$Cr together with $^{58}$Fe, $^{64}$Ni, $^{82}$Se and $^{86}$Kr and several of the Zn--Zr isotopes. A model in which no cECSNe occur, only tECSNe with low-mass FeCCSNe or rotating massive stars, proves also very successful at reproducing the solar abundances for these isotopes. Despite the small mass range for the progenitors of ECSNe and low-mass FeCCSNe, the large production factors suffice for the solar inventory of the above isotopes. Our model is compelling because it introduces no new tensions with the solar abundance distribution for a Milky Way model -- only tending to improve the model predictions for several isotopes. The proposed astrophysical production model thus provides a natural and elegant way to explain one of the last uncharted territories on the periodic table of astrophysical element production.
△ Less
Submitted 6 August, 2019;
originally announced August 2019.
-
Chromium Nucleosynthesis and Silicon-Carbon Shell Mergers in Massive Stars
Authors:
Benoit Côté,
Samuel Jones,
Falk Herwig,
Marco Pignatari
Abstract:
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 metallic…
▽ More
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.
△ Less
Submitted 17 June, 2019;
originally announced June 2019.
-
Galactic Chemical Evolution of Radioactive Isotopes
Authors:
Benoit Côté,
Maria Lugaro,
Rene Reifarth,
Marco Pignatari,
Blanka Világos,
Andrés Yagüe,
Brad K. Gibson
Abstract:
The presence of short-lived ($\sim$\,Myr) radioactive isotopes in meteoritic inclusions at the time of their formation represents a unique opportunity to study the circumstances that led to the formation of the Solar System. To interpret these observations we need to calculate the evolution of radioactive-to-stable isotopic ratios in the Galaxy. We present an extension of the open-source galactic…
▽ More
The presence of short-lived ($\sim$\,Myr) radioactive isotopes in meteoritic inclusions at the time of their formation represents a unique opportunity to study the circumstances that led to the formation of the Solar System. To interpret these observations we need to calculate the evolution of radioactive-to-stable isotopic ratios in the Galaxy. We present an extension of the open-source galactic chemical evolution codes NuPyCEE and JINAPyCEE that enables to track the decay of radioactive isotopes in the interstellar medium. We show how the evolution of isotopic ratio depends on the star formation history and efficiency, star-to-gas mass ratio, and galactic outflows. Given the uncertainties in the observations used to calibrate our model, our predictions for isotopic ratios at the time of formation of the Sun are uncertain by a factor of 3.6. At that time, to recover the actual radioactive-to-stable isotopic ratios predicted by our model, one can multiply the steady-state solution (see Equation~1) by $2.3^{+3.4}_{-0.7}$. However, in the cases where the radioactive isotope has a half-life longer than $\sim$\,200\,Myr, or the target radioactive or stable isotopes have mass- and/or metallicity-depended production rates, or they originate from different sources with different delay-time distributions, or the reference isotope is radioactive, our codes should be used for more accurate solutions. Our preliminary calculations confirm the dichotomy between radioactive nuclei in the early Solar System with $r$- and $s$-process origin, and that $^{55}$Mn and $^{60}$Fe can be explained by galactic chemical evolution, while $^{26}$Al cannot.
△ Less
Submitted 19 May, 2019;
originally announced May 2019.
-
Catching Element Formation In The Act
Authors:
Chris L. Fryer,
Frank Timmes,
Aimee L. Hungerford,
Aaron Couture,
Fred Adams,
Wako Aoki,
Almudena Arcones,
David Arnett,
Katie Auchettl,
Melina Avila,
Carles Badenes,
Eddie Baron,
Andreas Bauswein,
John Beacom,
Jeff Blackmon,
Stephane Blondin,
Peter Bloser,
Steve Boggs,
Alan Boss,
Terri Brandt,
Eduardo Bravo,
Ed Brown,
Peter Brown,
Steve Bruenn. Carl Budtz-Jorgensen,
Eric Burns
, et al. (194 additional authors not shown)
Abstract:
Gamma-ray astronomy explores the most energetic photons in nature to address some of the most pressing puzzles in contemporary astrophysics. It encompasses a wide range of objects and phenomena: stars, supernovae, novae, neutron stars, stellar-mass black holes, nucleosynthesis, the interstellar medium, cosmic rays and relativistic-particle acceleration, and the evolution of galaxies. MeV gamma-ray…
▽ More
Gamma-ray astronomy explores the most energetic photons in nature to address some of the most pressing puzzles in contemporary astrophysics. It encompasses a wide range of objects and phenomena: stars, supernovae, novae, neutron stars, stellar-mass black holes, nucleosynthesis, the interstellar medium, cosmic rays and relativistic-particle acceleration, and the evolution of galaxies. MeV gamma-rays provide a unique probe of nuclear processes in astronomy, directly measuring radioactive decay, nuclear de-excitation, and positron annihilation. The substantial information carried by gamma-ray photons allows us to see deeper into these objects, the bulk of the power is often emitted at gamma-ray energies, and radioactivity provides a natural physical clock that adds unique information. New science will be driven by time-domain population studies at gamma-ray energies. This science is enabled by next-generation gamma-ray instruments with one to two orders of magnitude better sensitivity, larger sky coverage, and faster cadence than all previous gamma-ray instruments. This transformative capability permits: (a) the accurate identification of the gamma-ray emitting objects and correlations with observations taken at other wavelengths and with other messengers; (b) construction of new gamma-ray maps of the Milky Way and other nearby galaxies where extended regions are distinguished from point sources; and (c) considerable serendipitous science of scarce events -- nearby neutron star mergers, for example. Advances in technology push the performance of new gamma-ray instruments to address a wide set of astrophysical questions.
△ Less
Submitted 7 February, 2019;
originally announced February 2019.
-
Nucleosynthetic Yields from Neutron Stars Accreting in Binary Common Envelopes
Authors:
James D. Keegans,
Chris L. Fryer,
Samuel W. Jones,
Benoit Cote,
Krzysztof Belczynski,
Falk Herwig,
Marco Pignatari,
Alison M. Laird,
Christian Aa. Diget
Abstract:
Massive-star binaries can undergo a phase where one of the two stars expands during its advanced evolutionary stage as a giant and envelops its companion, ejecting the hydrogen envelope and tightening its orbit. Such a common envelope phase is required to tighten the binary orbit in the formation of many of the observed X-ray binaries and merging compact binary systems. In the formation scenario f…
▽ More
Massive-star binaries can undergo a phase where one of the two stars expands during its advanced evolutionary stage as a giant and envelops its companion, ejecting the hydrogen envelope and tightening its orbit. Such a common envelope phase is required to tighten the binary orbit in the formation of many of the observed X-ray binaries and merging compact binary systems. In the formation scenario for neutron star binaries, the system might pass through a phase where a neutron star spirals into the envelope of its giant star companion. These phases lead to mass accretion onto the neutron star. Accretion onto these common-envelope-phase neutron stars can eject matter that has undergone burning near to the neutron star surface. This paper presents nucleosynthetic yields of this ejected matter, using population synthesis models to study the importance of these nucleosynthetic yields in a galactic chemical evolution context. Depending on the extreme conditions in temperature and density found in the accreted material, both proton-rich and neutron-rich nucleosynthesis can be obtained, with efficient production of neutron rich isotopes of low Z material at the most extreme conditions, and proton rich isotopes, again at low Z, in lower density models. Final yields are found to be extremely sensitive to the physical modeling of the accretion phase. We show that neutron stars accreting in binary common envelopes might be a new relevant site for galactic chemical evolution, and therefore more comprehensive studies are needed to better constrain nucleosynthesis in these objects.
△ Less
Submitted 5 February, 2019;
originally announced February 2019.
-
Following the Cosmic Evolution of Pristine Gas III: The Observational Consequences of the Unknown Properties of Population III Stars
Authors:
Richard Sarmento,
Evan Scannapieco,
Benoit Côté
Abstract:
We study the observational consequences of several unknown properties of Population III (Pop III) stars using large-scale cosmological simulations that include a subgrid model to track the unresolved mixing of pollutants. Varying the value of the critical metallicity that marks the boundary between Pop III and Population II (Pop II) star formation across 2 dex has a negligible effect on the fracti…
▽ More
We study the observational consequences of several unknown properties of Population III (Pop III) stars using large-scale cosmological simulations that include a subgrid model to track the unresolved mixing of pollutants. Varying the value of the critical metallicity that marks the boundary between Pop III and Population II (Pop II) star formation across 2 dex has a negligible effect on the fraction of Pop III stars formed and the subsequent fraction of Pop III flux from high-redshift galaxies. However, adopting a log normal initial mass function (IMF) for Pop III stars, in place of a baseline Salpeter IMF, results in a Pop III star formation rate density (SFRD) that is 1/4 of the baseline rate. The flux from high-redshift galaxies modeled with this IMF is highly bimodal, resulting in a tiny fraction of $z \leq 8$ galaxies with more than 75\% of their flux coming from Pop III stars. However, at $z=9$, right before reionization in our simulations, $\approx$ 20\% of galaxies are Pop III-bright with $m_{\rm UV} \le 31.4$ mag and at least 75\% of their flux generated by Pop III stars . Additionally, the log normal Pop III IMF results in a population of carbon enhanced, metal poor stars in reasonable agreement with MW halo observations. Our analysis supports the conclusion that the Pop III IMF was dominated by stars in the 20-120$M_{\odot}$ range that generate SN with carbon-enhanced ejecta.
△ Less
Submitted 11 January, 2019;
originally announced January 2019.
-
Neutron Star Mergers Might not be the Only Source of r-Process Elements in the Milky Way
Authors:
Benoit Côté,
Marius Eichler,
Almudena Arcones,
Camilla J. Hansen,
Paolo Simonetti,
Anna Frebel,
Chris L. Fryer,
Marco Pignatari,
Moritz Reichert,
Krzysztof Belczynski,
Francesca Matteucci
Abstract:
Probing the origin of r-process elements in the universe represents a multi-disciplinary challenge. We review the observational evidence that probe the properties of r-process sites, and address them using galactic chemical evolution simulations, binary population synthesis models, and nucleosynthesis calculations. Our motivation is to define which astrophysical sites have significantly contribute…
▽ More
Probing the origin of r-process elements in the universe represents a multi-disciplinary challenge. We review the observational evidence that probe the properties of r-process sites, and address them using galactic chemical evolution simulations, binary population synthesis models, and nucleosynthesis calculations. Our motivation is to define which astrophysical sites have significantly contributed to the total mass of r-process elements present in our Galaxy. We found discrepancies with the neutron star (NS-NS) merger scenario. Assuming they are the only site, the decreasing trend of [Eu/Fe] at [Fe/H]\,$>-1$ in the disk of the Milky Way cannot be reproduced while accounting for the delay-time distribution (DTD) of coalescence times ($\propto~t^{-1}$) derived from short gamma-ray bursts and population synthesis models. Steeper DTD functions ($\propto~t^{-1.5}$) or power laws combined with a strong burst of mergers before the onset of Type~Ia supernovae can reproduce the [Eu/Fe] trend, but this scenario is inconsistent with the similar fraction of short gamma-ray bursts and Type~Ia supernovae occurring in early-type galaxies, and reduces the probability of detecting GW170817 in an early-type galaxy. One solution is to assume an extra production site of Eu that would be active in the early universe, but would fade away with increasing metallicity. If this is correct, this extra site could be responsible for roughly 50% of the Eu production in the early universe, before the onset of Type~Ia supernovae. Rare classes of supernovae could be this additional r-process source, but hydrodynamic simulations still need to ensure the conditions for a robust r-process pattern.
△ Less
Submitted 4 March, 2019; v1 submitted 10 September, 2018;
originally announced September 2018.
-
Metal Mixing and Ejection in Dwarf Galaxies is Dependent on Nucleosynthetic Source
Authors:
Andrew Emerick,
Greg L. Bryan,
Mordecai-Mark Mac Low,
Benoit Côté,
Kathryn V. Johnston,
Brian W. O'Shea
Abstract:
Using a high resolution simulation of an isolated dwarf galaxy, accounting for multi-channel stellar feedback and chemical evolution on a star-by-star basis, we investigate how each of 15 metal species are distributed within our multi-phase interstellar medium (ISM) and ejected from our galaxy by galactic winds. For the first time, we demonstrate that the mass fraction probability distribution fun…
▽ More
Using a high resolution simulation of an isolated dwarf galaxy, accounting for multi-channel stellar feedback and chemical evolution on a star-by-star basis, we investigate how each of 15 metal species are distributed within our multi-phase interstellar medium (ISM) and ejected from our galaxy by galactic winds. For the first time, we demonstrate that the mass fraction probability distribution functions (PDFs) of individual metal species in the ISM are well described by a piecewise log-normal and power-law distribution. The PDF properties vary within each ISM phase. Hot gas is dominated by recent enrichment, with a significant power-law tail to high metal fractions, while cold gas is predominately log-normal. In addition, elements dominated by asymptotic giant branch (AGB) wind enrichment (e.g. N and Ba) mix less efficiently than elements dominated by supernova enrichment (e.g. $α$ elements and Fe). This result is driven by the differences in source energetics and source locations, particularly the higher chance compared to massive stars for AGB stars to eject material into cold gas. Nearly all of the produced metals are ejected from the galaxy (only 4% are retained), but over 20% of metals dominated by AGB enrichment are retained. In dwarf galaxies, therefore, elements synthesized predominately through AGB winds should be both overabundant and have a larger spread compared to elements synthesized in either core collapse or Type Ia supernovae. We discuss the observational implications of these results, their potential use in developing improved models of galactic chemical evolution, and their generalization to more massive galaxies.
△ Less
Submitted 7 November, 2018; v1 submitted 4 September, 2018;
originally announced September 2018.
-
FRIB and the GW170817 Kilonova
Authors:
A. Aprahamian,
R. Surman,
A. Frebel,
G. C. McLaughlin,
A. Arcones,
A. B. Balantekin,
J. Barnes,
Timothy C. Beers,
Erika M. Holmbeck,
Jinmi Yoon,
Maxime Brodeur,
T. M. Sprouse,
Nicole Vassh,
Jolie A. Cizewski,
Jason A. Clark,
Benoit Cote,
Sean M. Couch,
M. Eichler,
Jonathan Engel,
Rana Ezzeddine,
George M. Fuller,
Samuel A. Giuliani,
Robert Grzywacz,
Sophia Han,
C. J. Horowitz
, et al. (23 additional authors not shown)
Abstract:
In July 2018 an FRIB Theory Alliance program was held on the implications of GW170817 and its associated kilonova for r-process nucleosynthesis. Topics of discussion included the astrophysical and nuclear physics uncertainties in the interpretation of the GW170817 kilonova, what we can learn about the astrophysical site or sites of the r process from this event, and the advances in nuclear experim…
▽ More
In July 2018 an FRIB Theory Alliance program was held on the implications of GW170817 and its associated kilonova for r-process nucleosynthesis. Topics of discussion included the astrophysical and nuclear physics uncertainties in the interpretation of the GW170817 kilonova, what we can learn about the astrophysical site or sites of the r process from this event, and the advances in nuclear experiment and theory most crucial to pursue in light of the new data. Here we compile a selection of scientific contributions to the workshop, broadly representative of progress in r-process studies since the GW170817 event.
△ Less
Submitted 3 September, 2018;
originally announced September 2018.
-
r-Process Nucleosynthesis: Connecting Rare-Isotope Beam Facilities with the Cosmos
Authors:
C. J. Horowitz,
A. Arcones,
B. Côté,
I. Dillmann,
W. Nazarewicz,
I. U. Roederer,
H. Schatz,
A. Aprahamian,
D. Atanasov,
A. Bauswein,
J. Bliss,
M. Brodeur,
J. A. Clark,
A. Frebel,
F. Foucart,
C. J. Hansen,
O. Just,
A. Kankainen,
G. C. McLaughlin,
J. M. Kelly,
S. N. Liddick,
D. M. Lee,
J. Lippuner,
D. Martin,
J. Mendoza-Temis
, et al. (13 additional authors not shown)
Abstract:
This is an exciting time for the study of r-process nucleosynthesis. Recently, a neutron star merger GW170817 was observed in extraordinary detail with gravitational waves and electromagnetic radiation from radio to gamma rays. The very red color of the associated kilonova suggests that neutron star mergers are an important r-process site. Astrophysical simulations of neutron star mergers and core…
▽ More
This is an exciting time for the study of r-process nucleosynthesis. Recently, a neutron star merger GW170817 was observed in extraordinary detail with gravitational waves and electromagnetic radiation from radio to gamma rays. The very red color of the associated kilonova suggests that neutron star mergers are an important r-process site. Astrophysical simulations of neutron star mergers and core collapse supernovae are making rapid progress. Detection of both, electron neutrinos and antineutrinos from the next galactic supernova will constrain the composition of neutrino-driven winds and provide unique nucleosynthesis information. Finally FRIB and other rare-isotope beam facilities will soon have dramatic new capabilities to synthesize many neutron-rich nuclei that are involved in the r-process. The new capabilities can significantly improve our understanding of the r-process and likely resolve one of the main outstanding problems in classical nuclear astrophysics. However, to make best use of the new experimental capabilities and to fully interpret the results, a great deal of infrastructure is needed in many related areas of astrophysics, astronomy, and nuclear theory. We will place these experiments in context by discussing astrophysical simulations and observations of r-process sites, observations of stellar abundances, galactic chemical evolution, and nuclear theory for the structure and reactions of very neutron-rich nuclei. This review paper was initiated at a three-week International Collaborations in Nuclear Theory program in June 2016 where we explored promising r-process experiments and discussed their likely impact, and their astrophysical, astronomical, and nuclear theory context.
△ Less
Submitted 11 May, 2018;
originally announced May 2018.
-
Cyberhubs: Virtual Research Environments for Astronomy
Authors:
Falk Herwig,
Robert Andrassy,
Nic Annau,
Ondrea Clarkson,
Benoit Cote,
Aaron D'Sa,
Sam Jones,
Belaid Moa,
Jericho O'Connell,
David Porter,
Christian Ritter,
Paul Woodward
Abstract:
Collaborations in astronomy and astrophysics are faced with numerous cyber infrastructure challenges, such as large data sets, the need to combine heterogeneous data sets, and the challenge to effectively collaborate on those large, heterogeneous data sets with significant processing requirements and complex science software tools. The cyberhubs system is an easy-to-deploy package for small to med…
▽ More
Collaborations in astronomy and astrophysics are faced with numerous cyber infrastructure challenges, such as large data sets, the need to combine heterogeneous data sets, and the challenge to effectively collaborate on those large, heterogeneous data sets with significant processing requirements and complex science software tools. The cyberhubs system is an easy-to-deploy package for small to medium-sized collaborations based on the Jupyter and Docker technology, that allows web-browser enabled, remote, interactive analytic access to shared data. It offers an initial step to address these challenges. The features and deployment steps of the system are described, as well as the requirements collection through an account of the different approaches to data structuring, handling and available analytic tools for the NuGrid and PPMstar collaborations. NuGrid is an international collaboration that creates stellar evolution and explosion physics and nucleosynthesis simulation data. The PPMstar collaboration performs large-scale 3D stellar hydrodynamics simulation of interior convection in the late phases of stellar evolution. Examples of science that is presently performed on cyberhubs, in the areas 3D stellar hydrodynamic simulations, stellar evolution and nucleosynthesis and Galactic chemical evolution, are presented.
△ Less
Submitted 6 February, 2018;
originally announced February 2018.
-
i-Process Contribution of Rapidly Accreting White Dwarfs to the Solar Composition of First-Peak Neutron-Capture Elements
Authors:
Benoit Côté,
Pavel Denissenkov,
Falk Herwig,
Ashley J. Ruiter,
Christian Ritter,
Marco Pignatari,
Krzysztof Belczynski
Abstract:
Rapidly accreting white dwarfs (RAWDs) have been proposed as contributors to the chemical evolution of heavy elements in the Galaxy. Here, we test this scenario for the first time and determine the contribution of RAWDs to the solar composition of first-peak neutron-capture elements. We add the metallicity-dependent contribution of RAWDs to the one-zone galactic chemical evolution code OMEGA accor…
▽ More
Rapidly accreting white dwarfs (RAWDs) have been proposed as contributors to the chemical evolution of heavy elements in the Galaxy. Here, we test this scenario for the first time and determine the contribution of RAWDs to the solar composition of first-peak neutron-capture elements. We add the metallicity-dependent contribution of RAWDs to the one-zone galactic chemical evolution code OMEGA according to RAWD rates from binary stellar population models combined with metallicity-dependent $i$-process stellar yields calculated following the models of Denissenkov et al. (2017). With this approach we find that the contribution of RAWDs to the evolution of heavy elements in the Galaxy could be responsible for a significant fraction of the solar composition of Kr, Rb, Sr, Y, Zr, Nb, and Mo ranging from $2$ to $45\%$ depending on the element, the enrichment history of the Galactic gas, and the total mass ejected per RAWD. This contribution could explain the missing solar Lighter Element Primary Process for some elements (e.g., Sr, Y, and Zr). We do not overproduce any isotope relative to the solar composition, but $^{96}$Zr is produced in a similar amount. The $i$ process produces efficiently the Mo stable isotopes $^{95}$Mo and $^{97}$Mo. When nuclear reaction rate uncertainties are combined with our GCE uncertainties, the upper limits for the predicted RAWD contribution increase by a factor of $1.5-2$ for Rb, Sr, Y, and Zr, and by 3.8 and 2.4 for Nb and Mo, respectively. We discuss the implication of the RAWD stellar evolution properties on the single degenerate Type Ia supernova scenario.
△ Less
Submitted 21 January, 2018; v1 submitted 20 December, 2017;
originally announced December 2017.
-
SYGMA: Stellar Yields for Galactic Modeling Applications
Authors:
Christian Ritter,
Benoit Côté,
Falk Herwig,
Julio F. Navarro,
Chris L. Fryer
Abstract:
The stellar yields for galactic modeling applications (SYGMA) code is an open-source module that models the chemical ejecta and feedback of simple stellar populations (SSPs). It is intended for use in hydrodynamical simulations and semi-analytic models of galactic chemical evolution. The module includes the enrichment from asymptotic giant branch (AGB) stars, massive stars, SNIa and neutron-star m…
▽ More
The stellar yields for galactic modeling applications (SYGMA) code is an open-source module that models the chemical ejecta and feedback of simple stellar populations (SSPs). It is intended for use in hydrodynamical simulations and semi-analytic models of galactic chemical evolution. The module includes the enrichment from asymptotic giant branch (AGB) stars, massive stars, SNIa and neutron-star mergers. An extensive and extendable stellar yields library includes the NuGrid yields with all elements and many isotopes up to Bi. Stellar feedback from mechanic and frequency-dependent radiative luminosities are computed based on NuGrid stellar models and their synthetic spectra. The module further allows for customizable initial-mass functions and supernova Ia (SNIa) delay-time distributions to calculate time-dependent ejecta based on stellar yield input. A variety of r-process sites can be included. A comparison of SSP ejecta based on NuGrid yields with those from Portinari et al. (1998) and Marigo (2001) reveals up to a factor of 3.5 and 4.8 less C and N enrichment from AGB stars at low metallicity, a result we attribute to NuGrid's modeling of hot-bottom burning. Different core-collapse supernova explosion and fallback prescriptions may lead to substantial variations for the accumulated ejecta of C, O and Si in the first $10^7\, \mathrm{yr}$ at $Z=0.001$. An online interface of the open-source SYGMA module enables interactive simulations, analysis and data extraction of the evolution of all species formed by the evolution of simple stellar populations.
△ Less
Submitted 17 September, 2018; v1 submitted 24 November, 2017;
originally announced November 2017.
-
Validating Semi-Analytic Models of High-Redshift Galaxy Formation using Radiation Hydrodynamical Simulations
Authors:
Benoit Côté,
Devin Silvia,
Brian W. O'Shea,
Britton Smith,
John H. Wise
Abstract:
We use a cosmological hydrodynamic simulation calculated with Enzo and the semi-analytic galaxy formation model (SAM) GAMMA to address the chemical evolution of dwarf galaxies in the early universe. The long-term goal of the project is to better understand the origin of metal-poor stars and the formation of dwarf galaxies and the Milky Way halo by cross-validating these theoretical approaches. We…
▽ More
We use a cosmological hydrodynamic simulation calculated with Enzo and the semi-analytic galaxy formation model (SAM) GAMMA to address the chemical evolution of dwarf galaxies in the early universe. The long-term goal of the project is to better understand the origin of metal-poor stars and the formation of dwarf galaxies and the Milky Way halo by cross-validating these theoretical approaches. We combine GAMMA with the merger tree of the most massive galaxy found in the hydrodynamic simulation and compare the star formation rate, the metallicity distribution function (MDF), and the age-metallicity relationship predicted by the two approaches. We found that the SAM can reproduce the global trends of the hydrodynamic simulation. However, there are degeneracies between the model parameters and more constraints (e.g., star formation efficiency, gas flows) need to be extracted from the simulation to isolate the correct semi-analytic solution. Stochastic processes such as bursty star formation histories and star formation triggered by supernova explosions cannot be reproduced by the current version of GAMMA. Non-uniform mixing in the galaxy's interstellar medium, coming primarily from self-enrichment by local supernovae, causes a broadening in the MDF that can be emulated in the SAM by convolving its predicted MDF with a Gaussian function having a standard deviation of ~0.2 dex. We found that the most massive galaxy in the simulation retains nearby 100% of its baryonic mass within its virial radius, which is in agreement with what is needed in GAMMA to reproduce the global trends of the simulation.
△ Less
Submitted 28 March, 2018; v1 submitted 17 October, 2017;
originally announced October 2017.
-
The Origin of r-Process Elements in the Milky Way
Authors:
Benoit Côté,
Chris L. Fryer,
Krzysztof Belczynski,
Oleg Korobkin,
Martyna Chruślińska,
Nicole Vassh,
Matthew R. Mumpower,
Jonas Lippuner,
Trevor M. Sprouse,
Rebecca Surman,
Ryan Wollaeger
Abstract:
Some of the heavy elements, such as gold and europium (Eu), are almost exclusively formed by the rapid neutron capture process (r-process). However, it is still unclear which astrophysical site between core-collapse supernovae and neutron star - neutron star (NS-NS) mergers produced most of the r-process elements in the universe. Galactic chemical evolution (GCE) models can test these scenarios by…
▽ More
Some of the heavy elements, such as gold and europium (Eu), are almost exclusively formed by the rapid neutron capture process (r-process). However, it is still unclear which astrophysical site between core-collapse supernovae and neutron star - neutron star (NS-NS) mergers produced most of the r-process elements in the universe. Galactic chemical evolution (GCE) models can test these scenarios by quantifying the frequency and yields required to reproduce the amount of europium (Eu) observed in galaxies. Although NS-NS mergers have become popular candidates, their required frequency (or rate) needs to be consistent with that obtained from gravitational wave measurements. Here we address the first NS-NS merger detected by LIGO/Virgo (GW170817) and its associated Gamma-ray burst and analyze their implication on the origin of r-process elements. The range of NS-NS merger rate densities of 320-4740 Gpc$^{-3}$ yr$^{-1}$ provided by LIGO/Virgo is remarkably consistent with the range required by GCE to explain the Eu abundances in the Milky Way with NS-NS mergers, assuming the solar r-process abundance pattern for the ejecta. Under the same assumption, this event has produced about 1-5 Earth masses of Eu, and 3-13 Earth masses of gold. When using theoretical calculations to derive Eu yields, constraining the role of NS-NS mergers becomes more challenging because of nuclear astrophysics uncertainties. This is the first study that directly combines nuclear physics uncertainties with GCE calculations. If GW170817 is a representative event, NS-NS mergers can produce Eu in sufficient amounts and are likely to be the main r-process site.
△ Less
Submitted 15 January, 2018; v1 submitted 16 October, 2017;
originally announced October 2017.
-
On the Impact of Neutron Star Binaries Natal-Kick Distribution on the Galactic r-process Enrichment
Authors:
Mohammadtaher Safarzadeh,
Benoit Côté
Abstract:
We study the impact of the neutron star binaries' (NSBs) natal kick distribution on the Galactic r-process enrichment. We model the growth of a Milky Way type halo based on N-body simulation results and its star formation history based on multi epoch abundance matching techniques. We consider the NSBs that merge well beyond the galaxy's effective radius ($>2\times R_\mathrm{eff}$) do not contribut…
▽ More
We study the impact of the neutron star binaries' (NSBs) natal kick distribution on the Galactic r-process enrichment. We model the growth of a Milky Way type halo based on N-body simulation results and its star formation history based on multi epoch abundance matching techniques. We consider the NSBs that merge well beyond the galaxy's effective radius ($>2\times R_\mathrm{eff}$) do not contribute to Galactic r-process enrichment. Assuming a power-law delay-time distribution (DTD) function ($\propto t^{-1}$) with $t_\mathrm{min}=30$ Myr for binaries' coalescence timescales, and an exponential profile for their natal kick distribution with an average value of 180 km s$^{-1}$, we show that up to $\sim$ 40% of all formed NSBs do not contribute to r-process enrichment by $z=0$, either because they merge far from the galaxy at a given redshift (up to $\sim$ 25%) or have not yet merged by today ($\sim$15%). Our result is largely insensitive to the details of the DTD function. Assuming a constant coalescence timescale of 100 Myr well approximates the adopted DTD with 30% of the NSBs not contributing to r-process enrichment. Our results, although rather dependent on the adopted natal kick distribution, represent a first step towards estimating the impact of natal kicks and DTD functions on r-process enrichment of galaxies that would need to be incorporated in the hydrodynamical simulations.
△ Less
Submitted 24 July, 2017;
originally announced July 2017.
-
Convective-reactive nucleosynthesis of K, Sc, Cl and p-process isotopes in O-C shell mergers
Authors:
C. Ritter,
R. Andrassy,
B. Côté,
F. Herwig,
P. R. Woodward,
M. Pignatari,
S. Jones
Abstract:
We address the deficiency of odd-Z elements P, Cl, K and Sc in galactic chemical evolution models through an investigation of the nucleosynthesis of interacting convective O- and C shells in massive stars. 3D hydrodynamic simulations of O-shell convection with moderate C-ingestion rates show no dramatic deviation from spherical symmetry. We derive a spherically averaged diffusion coefficient for 1…
▽ More
We address the deficiency of odd-Z elements P, Cl, K and Sc in galactic chemical evolution models through an investigation of the nucleosynthesis of interacting convective O- and C shells in massive stars. 3D hydrodynamic simulations of O-shell convection with moderate C-ingestion rates show no dramatic deviation from spherical symmetry. We derive a spherically averaged diffusion coefficient for 1D nucleosynthesis simulations which show that such convective-reactive ingestion events can be a production site for P, Cl, K and Sc. An entrainment rate of $10^{-3}{{\rm M}_\odot}$/s features overproduction factors $OP_\mathrm{s} \approx 7$. Full O-C shell mergers in our 1D stellar evolution massive star models have overproduction factors $OP_\mathrm{m}>1 \mathrm{dex}$ but for such cases 3D hydrodynamic simulations suggest deviations from spherical symmetry. p-process species can be produced with overproduction factors of $OP_\mathrm{m}>1 \mathrm{dex}$, e.g. for $^{130,132}Ba$. Using the uncertain prediction of the $15{{\rm M}_\odot}$, $Z=0.02$ massive star model ($OP_\mathrm{m} \approx 15$) as representative for merger or entrainment convective-reactive events involving O- and C-burning shells, and assume that such events occur in more than 50% of all stars, our chemical evolution models reproduce the observed Galactic trends of the odd-Z elements.
△ Less
Submitted 25 September, 2017; v1 submitted 19 April, 2017;
originally announced April 2017.
-
Advanced LIGO Constraints on Neutron Star Mergers and R-Process Sites
Authors:
Benoit Côté,
Krzysztof Belczynski,
Chris L. Fryer,
Christian Ritter,
Adam Paul,
Benjamin Wehmeyer,
Brian W. O'Shea
Abstract:
The role of compact binary mergers as the main production site of r-process elements is investigated by combining stellar abundances of Eu observed in the Milky Way, galactic chemical evolution (GCE) simulations, binary population synthesis models, and Advanced LIGO gravitational wave measurements. We compiled and reviewed seven recent GCE studies to extract the frequency of neutron star - neutron…
▽ More
The role of compact binary mergers as the main production site of r-process elements is investigated by combining stellar abundances of Eu observed in the Milky Way, galactic chemical evolution (GCE) simulations, binary population synthesis models, and Advanced LIGO gravitational wave measurements. We compiled and reviewed seven recent GCE studies to extract the frequency of neutron star - neutron star (NS-NS) mergers that is needed in order to reproduce the observed [Eu/Fe] vs [Fe/H] relationship. We used our simple chemical evolution code to explore the impact of different analytical delay-time distribution (DTD) functions for NS-NS mergers. We then combined our metallicity-dependent population synthesis models with our chemical evolution code to bring their predictions, for both NS-NS mergers and black hole - neutron star mergers, into a GCE context. Finally, we convolved our results with the cosmic star formation history to provide a direct comparison with current and upcoming Advanced LIGO measurements. When assuming that NS-NS mergers are the exclusive r-process sites, and that the ejected r-process mass per merger event is 0.01 Msun, the number of NS-NS mergers needed in GCE studies is about 10 times larger than what is predicted by standard population synthesis models. These two distinct fields can only be consistent with each other when assuming optimistic rates, massive NS-NS merger ejecta, and low Fe yields for massive stars. Population synthesis models and GCE simulations are in agreement with the current upper limit (O1) established by Advanced LIGO during their first run of observations. Upcoming measurements will provide important constraints on the actual local NS-NS merger rate, will provide insights on the plausibility of the GCE requirement, and will help to define whether or not compact binary mergers can be the dominant source of r-process elements in the Universe.
△ Less
Submitted 12 February, 2017; v1 submitted 7 October, 2016;
originally announced October 2016.
-
JINA-NuGrid Galactic Chemical Evolution Pipeline
Authors:
Benoit Côté,
Christian Ritter,
Falk Herwig,
Brian W. O'Shea,
Marco Pignatari,
Devin Silvia,
Samuel Jones,
Chris L. Fryer
Abstract:
Galactic chemical evolution is a topic that involves nuclear physics, stellar evolution, galaxy evolution, observation, and cosmology. Continuous communication and feedback between these fields is a key component in improving our understanding of how galaxies form and how elements are created and recycled in galaxies and intergalactic space. In this proceedings, we present the current state of the…
▽ More
Galactic chemical evolution is a topic that involves nuclear physics, stellar evolution, galaxy evolution, observation, and cosmology. Continuous communication and feedback between these fields is a key component in improving our understanding of how galaxies form and how elements are created and recycled in galaxies and intergalactic space. In this proceedings, we present the current state of the JINA-NuGrid chemical evolution pipeline. It is designed to probe the impact of nuclear astrophysics uncertainties on galactic chemical evolution, to improve our knowledges regarding the origin of the elements in a cosmological context, and to create the required interdisciplinary connections.
△ Less
Submitted 29 September, 2016;
originally announced September 2016.
-
The Impact of Modeling Assumptions in Galactic Chemical Evolution Models
Authors:
Benoit Côté,
Brian W. O'Shea,
Christian Ritter,
Falk Herwig,
Kim A. Venn
Abstract:
We use the OMEGA galactic chemical evolution code to investigate how the assumptions used for the treatment of galactic inflows and outflows impact numerical predictions. The goal is to determine how our capacity to reproduce the chemical evolution trends of a galaxy is affected by the choice of implementation used to include those physical processes. In pursuit of this goal, we experiment with th…
▽ More
We use the OMEGA galactic chemical evolution code to investigate how the assumptions used for the treatment of galactic inflows and outflows impact numerical predictions. The goal is to determine how our capacity to reproduce the chemical evolution trends of a galaxy is affected by the choice of implementation used to include those physical processes. In pursuit of this goal, we experiment with three different prescriptions for galactic inflows and outflows and use OMEGA within a Markov Chain Monte Carlo code to recover the set of input parameters that best reproduces the chemical evolution of nine elements in the dwarf spheroidal galaxy Sculptor. This provides a consistent framework for comparing the best-fit solutions generated by our different models. Despite their different degrees of intended physical realism, we found that all three prescriptions can reproduce in an almost identical way the stellar abundance trends observed in Sculptor. While the three models have the same capacity to fit the data, the best values recovered for the parameters controlling the number of Type Ia supernovae and the strength of galactic outflows, are substantially different and in fact mutually exclusive from one model to another. For the purpose of understanding how a galaxy evolves, we conclude that only reproducing the evolution of a limited number of elements is insufficient and can lead to misleading conclusions. More elements or additional constraints such as the galaxy's star formation efficiency and the gas fraction are needed in order to break the degeneracy between the different modeling assumptions. Our results show that the successes and failures of chemical evolution models are predominantly driven by the input stellar yields, rather than by the complexity of the galaxy model itself. Simple models such as OMEGA are therefore sufficient to test and validate stellar yields.
△ Less
Submitted 7 December, 2016; v1 submitted 26 April, 2016;
originally announced April 2016.
-
Mass and Metallicity Requirement in Stellar Models for Galactic Chemical Evolution Applications
Authors:
Benoit Côté,
Christopher West,
Alexander Heger,
Christian Ritter,
Brian W. O'Shea,
Falk Herwig,
Claudia Travaglio,
Sara Bisterzo
Abstract:
We used a one-zone chemical evolution model to address the question of how many masses and metallicities are required in grids of massive stellar models in order to ensure reliable galactic chemical evolution predictions. We used a set of yields that includes seven masses between 13 and 30 Msun, 15 metallicities between 0 and 0.03 in mass fraction, and two different remnant mass prescriptions. We…
▽ More
We used a one-zone chemical evolution model to address the question of how many masses and metallicities are required in grids of massive stellar models in order to ensure reliable galactic chemical evolution predictions. We used a set of yields that includes seven masses between 13 and 30 Msun, 15 metallicities between 0 and 0.03 in mass fraction, and two different remnant mass prescriptions. We ran several simulations where we sampled subsets of stellar models to explore the impact of different grid resolutions. Stellar yields from low- and intermediate-mass stars and from Type Ia supernovae have been included in our simulations, but with a fixed grid resolution. We compared our results with the stellar abundances observed in the Milky Way for O, Na, Mg, Si, Ca, Ti, and Mn. Our results suggest that the range of metallicity considered is more important than the number of metallicities within that range, which only affects our numerical predictions by about 0.1 dex. We found that our predictions at [Fe/H] < -2 are very sensitive to the metallicity range and the mass sampling used for the lowest metallicity included in the set of yields. Variations between results can be as high as 0.8 dex, for any remnant mass prescription. At higher [Fe/H], we found that the required number of masses depends on the element of interest and on the remnant mass prescription. With a monotonic remnant mass prescription where every model explodes as a core-collapse supernova, the mass resolution induces variations of 0.2 dex on average. But with a remnant mass prescription that includes islands of non-explodability, the mass resolution can cause variations of about 0.2 to 0.7 dex depending on the choice of metallicity range. With such a prescription, explosive or non-explosive models can be missed if not enough masses are selected, resulting in over- or under-estimations of the mass ejected by massive stars.
△ Less
Submitted 19 August, 2016; v1 submitted 15 February, 2016;
originally announced February 2016.
-
Uncertainties in Galactic Chemical Evolution Models
Authors:
Benoit Côté,
Christian Ritter,
Brian W. O'Shea,
Falk Herwig,
Marco Pignatari,
Samuel Jones,
Chris Fryer
Abstract:
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…
▽ More
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.
△ Less
Submitted 11 May, 2016; v1 submitted 21 September, 2015;
originally announced September 2015.
-
Galactic Outflows and Evolution of the Interstellar Medium
Authors:
Benoit Côté,
Hugo Martel,
Laurent Drissen,
Carmelle Robert
Abstract:
We present a model to self-consistently describe the joint evolution of starburst galaxies and the galactic wind resulting from this evolution. We combine the population synthesis code Starburst99 with a semi-analytical model of galactic outflows and a model for the distribution and abundances of chemical elements inside the outflows. Starting with a galaxy mass, formation redshift, and adopting a…
▽ More
We present a model to self-consistently describe the joint evolution of starburst galaxies and the galactic wind resulting from this evolution. We combine the population synthesis code Starburst99 with a semi-analytical model of galactic outflows and a model for the distribution and abundances of chemical elements inside the outflows. Starting with a galaxy mass, formation redshift, and adopting a particular form for the star formation rate, we describe the evolution of the stellar populations in the galaxy, the evolution of the metallicity and chemical composition of the interstellar medium (ISM), the propagation of the galactic wind, and the metal-enrichment of the intergalactic medium (IGM). In this paper, we study the properties of the model, by varying the mass of the galaxy, the star formation rate, and the efficiency of star formation. Our main results are the following: (1) For a given star formation efficiency f*, a more extended period of active star formation tends to produce a galactic wind that reaches a larger extent. If f* is sufficiently large, the energy deposited by the stars completely expels the ISM. Eventually, the ISM is being replenished by mass loss from supernovae and stellar winds. (2) For galaxies with masses above 10^11 Msun, the material ejected in the IGM always falls back onto the galaxy. Hence lower-mass galaxies are the ones responsible for enriching the IGM. (3) Stellar winds play a minor role in the dynamical evolution of the galactic wind, because their energy input is small compared to supernovae. However, they contribute significantly to the chemical composition of the galactic wind. We conclude that the history of the ISM enrichment plays a determinant role in the chemical composition and extent of the galactic wind, and therefore its ability to enrich the IGM.
△ Less
Submitted 9 December, 2011;
originally announced December 2011.