Search

Presolar grains are stardust particles that condensed in the ejecta or in the outflows of dying stars and can today be extracted from meteorites. They recorded the nucleosynthetic fingerprint of their parent stars and thus serve as valuable probes of these astrophysical sites. The most common types of presolar silicon carbide grains (called mainstream SiC grains) condensed in the outflows of asymptotic giant branch stars. Their measured silicon isotopic abundances are not significantly influenced by nucleosynthesis within the parent star, but rather represents the pristine stellar composition. Silicon isotopes can thus be used as a proxy for galactic chemical evolution. However, the measured correlation of $^{29}$Si/$^{28}$Si versus $^{30}$Si/$^{28}$Si does not agree with any current chemical evolution model. Here, we use a Monte Carlo model to vary nuclear reaction rates within their theoretical or experimental uncertainties and process them through stellar nucleosynthesis and galactic chemical evolution models to study the variation of silicon isotope abundances based on these nuclear reaction rate uncertainties. We find that these uncertainties can indeed be responsible for the discrepancy between measurements and models and that the slope of the silicon isotope correlation line measured in mainstream SiC grains agrees with chemical evolution models within the nuclear reaction rate uncertainties. Our result highlights the importance of future precision reaction rate measurements for resolving the apparent data-model discrepancy.

We investigate how stellar feedback from the first stars (Population III) distributes metals through the interstellar and intergalactic medium using the star-by-star cosmological hydrodynamics simulation, Aeos. We find that energy injected from the supernovae of the first stars is enough to expel a majority of gas and injected metals beyond the virial radius of halos with mass $M_* \lesssim 10^7$ M$_\odot$, regardless of the number of supernovae. This prevents self-enrichment and results in a non-monotonic increase in metallicity at early times. Most minihalos ($M \gtrsim 10^5 \, \rm M_\odot$) do not retain significant fractions of the yields produced within their virial radii until they have grown to halo masses of $M \gtrsim 10^7 \, \rm M_\odot$. The loss of metals to regions well beyond the virial radius delays the onset of enriched star formation and extends the period that Population III star formation can persist. We also explore the contributions of different nucleosynthetic channels to 10 individual elements. On the timescale of the simulation (lowest redshift $z=14.3$), enrichment is dominated by core-collapse supernovae for all elements, but with a significant contribution from asymptotic giant branch winds to the s-process elements, which are normally thought to only be important at late times. In this work, we establish important mechanisms for early chemical enrichment which allows us to apply Aeos in later epochs to trace the evolution of enrichment during the complete transition from Population III to Population II stars., Comment: 16 pages, 12 figures, 1 animated figure (see source files), submitted to ApJ

We study the sensitivity of presupernova evolution and supernova nucleosynthesis yields of massive stars to variations of the initial composition. We use the solar abundances from Lodders (2009), and compute two different initial stellar compositions: i) scaled solar abundances, and ii) the isotopic galactic chemical history model (GCH) developed by West and Heger (2013b). We run a grid of models using the KEPLER stellar evolution code, with 7 initial stellar masses, 12 initial metallicities, and two for each scaling method to explore the effects on nucleosynthesis over a metallicity range of $-4.0\leq[Z]\leq+0.3$. We find that the compositions from the GCH model better reproduce the weak \emph{s}-process peak than the scaled solar models. The model yields are then used in the OMEGA Galactic Chemical Evolution (GCE) code to assess this result further. We find that initial abundances used in computing stellar structure have more of an impact on GCE results than initial abundances used in the burn network, with the GCH model again being favored when compared to observations. Lastly, a machine learning algorithm was used to verify the free parameter values of the GCH model, which were previously found by West and Heger (2013b) using a stochastic fitting process. The updated model is provided as an accessible tool for further nucleosynthesis studies.

The AEOS project introduces a series of high-resolution cosmological simulations that model star-by-star chemical enrichment and galaxy formation in the early Universe, achieving 1 pc resolution. These simulations capture the complexities of galaxy evolution within the first ~300 Myr by modeling individual stars and their feedback processes. By incorporating chemical yields from individual stars, AEOS generates galaxies with diverse stellar chemical abundances, linking them to hierarchical galaxy formation and early nucleosynthetic events. These simulations underscore the importance of chemical abundance patterns in ancient stars as vital probes of early nucleosynthesis, star formation histories, and galaxy formation. We examine the metallicity floors of various elements resulting from Pop III enrichment, providing best-fit values for eight different metals (e.g., [O/H] = -4.0) to guide simulations without Pop III models. Additionally, we identify galaxies that begin star formation with Pop II after external enrichment and investigate the frequency of CEMP stars at varying metallicities. The AEOS simulations offer detailed insights into the relationship between star formation, feedback, and chemical enrichment. Future work will extend these simulations to later epochs to interpret the diverse stellar populations of the Milky Way and its satellites., Comment: 19 pages, 12 figures, accepted for publication in ApJ

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.

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$\sigma$ 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., Comment: 14 pages, 11 figures, accepted by MNRAS

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., Comment: 36 pages, 26 figures, 1 Table, 1 Appendix, Accepted for publication in MNRAS

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., Comment: 22 pages, 11 figures

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., Comment: 14 pages, 10 figures, Accepted for publication in MNRAS

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 ($\delta$) 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 $\delta$), 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 $\delta$), we derive that the time elapsed from the formation of the molecular cloud to the formation of the Sun was 9-16 Myr., Comment: 17 pages, 3 figures, published on Universe as part of the Special Issue "Nuclear Astrophysics in the Era of High Precision Astronomy"

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., Comment: Accepted for publication in ApJL 28/7/2022

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., Comment: accepted for publication on Universe Special Issue AGB Stars - In Honor of Professor Maurizio Busso on the Occasion of His 70th Birthday

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 $\alpha$-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., Comment: accepted for publication by A&A

$\textit{Thailand-UK Python+Astronomy Summer School}$ (ThaiPASS), a collaborative project comprising UK and Thai institutions and assess its impact and possible application to schools in the United Kingdom. Since its inception in 2018, the annual ThaiPASS have trained around 60 Thai high-school students in basic data handling skills using Python in the context of various astronomy topics, using current research from the teaching team. Our impact assessment of the 5 day summer schools show an overwhelmingly positive response from students in both years, with over $80\%$ of students scoring the activities above average in all activities but one. We use this data to suggest possible future improvements. We also discuss how ThaiPASS may inspire further outreach and engagement activities within the UK and beyond., Comment: Accepted for publication into Physics Education; 21 pages, 5 figures and 3 tables

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., Comment: 28 pages, 17 figures, accepted for publication in A&A

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., Comment: 36 pages, 7 figures, 7 tables

Short-lived radionuclides (SLRs) with mean-lives $\tau$ 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 $\gamma$. 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 $\tau$, $\gamma$, 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 $\tau/\gamma\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 $\tau/\gamma=1$. For $\tau/\gamma\lesssim1$, we provide the probability for the SLR abundance to carry the signature of only one enrichment event, which is greater than 50% when $\tau/\gamma\lesssim0.3$. For $0.3\lesssim\tau/\gamma\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 $\tau/\gamma\sim3$ for these isotopes., Comment: 25 pages, 11 figures, 4 tables, submitted to ApJ

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., Comment: 8 pages, 2 figures, NPA-IX proceedings

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., Comment: White paper submitted to the Canadian Long Range Plan 2020. Minor formatting changes relative to submitted version

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 \AA~ for 209 stars and Ru I lines at 4080, 4584, 4757 \AA~ 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., Comment: 23 pages, 12 figures, 7 tables, accepted in MNRAS

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., Comment: 8 pages, 5 figures, submitted to ApJL

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., Comment: 23 pages, 9 figures, 3 tables, accepted in ApJ

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., Comment: 14 pages including 3 figures

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., Comment: 21 pages, 14 figures

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., Comment: Accepted by ApJ

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., Comment: 29 pages, 10 figures, 1 table, re-submitted to ApJ (revised version)

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. $\alpha$ 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., Comment: 18 pages, 7 figures (plus 2 page, 2 figure appendix). Accepted to ApJ

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., Comment: Proceedings for the FRIB Theory Alliance workshop "FRIB and the GW170817 kilonova", held 16-27 July 2018 at the Facility for Rare Isotope Beams, Michigan State University, East Lansing, MI USA; 44 pages, 14 figures

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., Comment: 19 pages, 9 figures, ApJ Supplement Special Issue on Data, accepted

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., Comment: 19 pages, 13 figures, 3 tables. Submitted to ApJ (version 2)

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., Comment: 18 pages, 10 figures, 3 tables, published in ApJS

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., Comment: 26 pages, 13 figures, 2 tables, submitted to ApJ (version 2)

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., Comment: 13 pages, 3 figures, 3 tables, submitted to ApJ (version 2)

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., Comment: 7 pages, 3 figures, MNRAS Accepted

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., Comment: Accepted in ApJ, 22 pages, 14 figures, 1 table

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., Comment: 3 pages, 2 figures, submitted to JPS Conference Proceedings, Nuclei in the Cosmos XIV

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., Comment: 16 pages, 9 figures, 1 table, accepted in ApJ

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., Comment: 14 pages, 2 tables, 14 figures, submitted to MNRAS

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., Comment: 20 pages, 12 figures, 7 tables, accepted in ApJ (v.3 - 11/05/2016)

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., Comment: 17 pages, 14 color figures. Accepted for publication in MNRAS