Search

The use of computer simulations to model physical systems has gained significant traction in recent years. A key factor in ensuring the accuracy of these models is the proper calibration of model parameters based on real-world observations or experimental data. Inevitably, uncertainties arise, and Bayesian methods provide a robust framework for quantifying and propagating these uncertainties to model predictions. However, predictions can become inaccurate if model errors are neglected. A promising approach to address this issue involves embedding a bias term in the inference parameters, allowing the quantified bias to influence non-observed Quantities of Interest (QoIs). This paper introduces a more interpretable framework for bias embedding compared to existing methods. Current likelihood formulations that incorporate embedded bias often fail when measurement noise is present. To overcome these limitations, we adapt the existing likelihood models to properly account for noise and propose two new formulations designed to address the shortcomings of the previous approaches. Moreover, we evaluate the performance of this bias-embedding approach in the presence of discrepancies between measurements and model predictions, including noise and outliers. Particular attention is given to how the uncertainty associated with the bias term propagates to the QoIs, enabling a more comprehensive statistical analysis of prediction reliability. Finally, the proposed embedded bias model is applied to estimate the uncertainty in the predicted heat flux from a transient thermal simulation, using temperature observations to illustrate its effectiveness., Comment: 37 pages, 24 figures, 5 tables

We provide yields from 189 neutrino-driven core-collapse supernova (CCSN) simulations covering zero-age main sequence masses between 11 and 75 solar masses and three different metallicities. Our CCSN simulations have two main advantages compared to previous methods used for applications in Galactic chemical evolution (GCE). Firstly, the mass cut between remnant and ejecta evolves naturally. Secondly, the neutrino luminosities and thus the electron fraction are not modified. Both is key to obtain an accurate nucleosynthesis. We follow the composition with an in-situ nuclear reaction network including the 16 most abundant isotopes and use the yields as input in a GCE model of the Milky Way. We adopt a GCE which takes into account infall of gas as well as nucleosynthesis from a large variety of stellar sources. The GCE model is calibrated to reproduce the main features of the solar vicinity. For the CCSN models, we use different calibrations and propagate the uncertainty. We find a big impact of the CCSN yields on our GCE predictions. We compare the abundance ratios of C, O, Ne, Mg, Si, S, Ar, Ca, Ti, and Cr with respect to Fe to an observational data set as homogeneous as possible. From this, we conclude that at least half of the massive stars have to explode to match the observed abundance ratios. If the explosions are too energetic, the high amount of iron will suppress the abundance ratios. With this, we demonstrate how GCE models can be used to constrain the evolution and deaths of massive stars., Comment: 13 pages, 10 figures, submitted to MNRAS

Material expelled from binary neutron star (BNS) mergers can harbor r-process nucleosynthesis and power a Kilonova (KN), both intimately related to the astrophysical conditions of the ejection. In turn such conditions indirectly depend on the equation of state (EOS) describing matter inside the neutron star. Therefore, in principle the above observables can hold valuable information on nuclear matter, as the merger gravitational wave signal already does. In this work, we consider the outcome of a set of BNS merger simulations employing different finite-temperature nuclear EOSs. The latter are obtained from a Skyrme-type interaction model where nuclear properties, such as the incompressibility and the nucleon effective mass at saturation density, are systematically varied. We post-process the ejecta using a reaction network coupled with a semi-analytic KN model, to asses the sensitivity on the input EOS of the final yields and the KN light curves. Both of them are found to be non-trivially influenced by the EOS, with the overall outcome being dominated by the heterogeneous outflows from the remnant disk, hosting a variable degree of neutron-rich material. The dynamical ejecta can be more directly related to the EOS parameters considered, however, we find its role in the yields production and the KN emission too entangled with the other ejecta components, in order to infer solid correlations. This result highlights the strong degeneracy that intervenes between the merger outcome and the behaviour of the intrinsic nuclear matter, and places itself as a limit to the employment of EOS-constraining approaches of such kind.

We present here an extended nuclear network, with 90 species, designed for being coupled with hydrodynamic simulations, which includes neutrons, protons, electrons, positrons, and the corresponding neutrino and anti-neutrino emission. This network is also coupled with temperature, making it extremely robust and, together with its size, unique of its kind. The inclusion of electron captures on free protons makes the network very appropriate for multidimensional studies of Type Ia supernova explosions, especially when the exploding object is a massive white dwarf. The results obtained with the proposed medium-sized network compare fairly well, to a few percent, with those computed with the extended network WinNet (> 2000 isotopes) in scenarios reproducing the gross physical conditions of current Type Ia supernova explosion models. In those cases where the carbon and oxygen fuel ignites at high density, the high-temperature plateau typical of the nuclear statistical equilibrium regime is well defined and stable, allowing large integration time steps. We show that the inclusion of electron captures on free protons substantially improves the estimation of the electron fraction of the mixture. Therefore, the pressure is better determined than in networks where electron captures are excluded, which will ultimately lead to more reliable hydrodynamic models. Explosive combustion of helium at low density, occurring near the surface layer of a white dwarf, is also better described with the proposed network, which gives nuclear energy generation rates much closer to WinNet than typical reduced alpha networks., Comment: 19 pages, 21 Figures and 5 Tables. Accepted for publication in A&A

The production of heavy elements is one of the main by-products of the explosive end of massive stars. A long sought goal is finding differentiated patterns in the nucleosynthesis yields, which could permit identifying a number of properties of the explosive core. Among them, the traces of the magnetic field topology are particularly important for \emph{extreme} supernova explosions, most likely hosted by magnetorotational effects. We investigate the nucleosynthesis of five state-of-the-art magnetohydrodynamic models with fast rotation that have been previously calculated in full 3D and that involve an accurate neutrino transport (M1). One of the models does not contain any magnetic field and synthesizes elements around the iron group, in agreement with other CC-SNe models in literature. All other models host a strong magnetic field of the same intensity, but with different topology. For the first time, we investigate the nucleosynthesis of MR-SNe models with a quadrupolar magnetic field and a 90 degree tilted dipole. We obtain a large variety of ejecta compositions reaching from iron nuclei to nuclei up to the third r-process peak. We assess the robustness of our results by considering the impact of different nuclear physics uncertainties such as different nuclear masses, $\beta^{-}$-decays and $\beta^{-}$-delayed neutron emission probabilities, neutrino reactions, fission, and a feedback of nuclear energy on the temperature. We find that the qualitative results do not change with different nuclear physics input. The properties of the explosion dynamics and the magnetic field configuration are the dominant factors determining the ejecta composition.

The elemental abundances between strontium and silver ($Z = 38-47$) observed in the atmospheres of very metal-poor stars (VMP) in the Galaxy may contain the fingerprint of the weak $r$-process and $\nu p$-process occurring in early core-collapse supernovae explosions. In this work, we combine various astrophysical conditions based on a steady-state model to cover the richness of the supernova ejecta in terms of entropy, expansion timescale, and electron fraction. The calculated abundances based on different combinations of conditions are compared with stellar observations with the aim of constraining supernova ejecta conditions. We find that some conditions of the neutrino-driven outflows consistently reproduce the observed abundances of our sample. In addition, from the successful combinations, the neutron-rich trajectories better reproduce the observed abundances of Sr-Zr ($Z= 38-40$), while the proton-rich ones, Mo-Pd ($Z= 42-47$)., Comment: 18 pages, 11 figures, submitted to ApJ

Simulation-based digital twins must provide accurate, robust and reliable digital representations of their physical counterparts. Quantifying the uncertainty in their predictions plays, therefore, a key role in making better-informed decisions that impact the actual system. The update of the simulation model based on data must be then carefully implemented. When applied to complex standing structures such as bridges, discrepancies between the computational model and the real system appear as model bias, which hinders the trustworthiness of the digital twin and increases its uncertainty. Classical Bayesian updating approaches aiming to infer the model parameters often fail at compensating for such model bias, leading to overconfident and unreliable predictions. In this paper, two alternative model bias identification approaches are evaluated in the context of their applicability to digital twins of bridges. A modularized version of Kennedy and O'Hagan's approach and another one based on Orthogonal Gaussian Processes are compared with the classical Bayesian inference framework in a set of representative benchmarks. Additionally, two novel extensions are proposed for such models: the inclusion of noise-aware kernels and the introduction of additional variables not present in the computational model through the bias term. The integration of such approaches in the digital twin corrects the predictions, quantifies their uncertainty, estimates noise from unknown physical sources of error and provides further insight into the system by including additional pre-existing information without modifying the computational model., Comment: 32 pages, 21 figures, 9 tables. Submitted for consideration to Applied Stochastic Models in Business and Industry

The dynamics in mergers of binary neutron star (BNS) systems depend sensitively on the equation of state (EOS) of dense matter. This has profound implications on the emission of gravitational waves (GWs) and the ejection of matter in the merger and post-merger phases and is thus of high interest for multi-messenger astronomy. Today, a variety of nuclear EOSs are available with various underlying microphysical models. This calls for a study to focus on EOS effects from different physical nuclear matter properties and their influence on BNS mergers. We perform simulations of equal-mass BNS mergers with a set of 9 different EOSs based on Skyrme density functionals. In the models, we systematically vary the effective nucleon mass, incompressibility, and symmetry energy at saturation density. This allows us to investigate the influence of specific nuclear matter properties on the dynamics of BNS mergers. We analyze the impact of these properties on the merger dynamics, the fate of the remnant, disk formation, ejection of matter, and gravitational wave emission. Our results indicate that some aspects of the merger are sensitive to the EOS around saturation density while others are sensitive to the behavior towards higher densities, e.g., characterized by the slope of the pressure as a function of density. The detailed density dependence of the EOS thus needs to be taken into account to describe its influence on BNS mergers., Comment: 17 pages, 19 figures, published version

We present the state-of-the-art single-zone nuclear reaction network WinNet that is capable of calculating the nucleosynthetic yields of a large variety of astrophysical environments and conditions. This ranges from the calculation of the primordial nucleosynthesis, where only a few nuclei are considered, to the ejecta of neutron star mergers with several thousands of involved nuclei. Here we describe the underlying physics and implementation details of the reaction network. We additionally present the numerical implementation of two different integration methods, the implicit Euler method and Gears method along with their advantages and disadvantages. We furthermore describe basic example cases of thermodynamic conditions that we provide together with the network and demonstrate the reliability of the code by using simple test cases. With this publication, WinNet is publicly available and open source at GitHub and Zenodo.

We present core-collapse supernova simulations including nuclear reaction networks that impact explosion dynamics and nucleosynthesis. The different composition treatment can lead to changes in the neutrino heating in the vicinity of the shock by modifying the number of nucleons and thus the neutrino-opacity of the region. This reduces the ram pressure outside the shock and allows an easier expansion. The energy released by the nuclear reactions during collapse also slows down the accretion and aids the shock expansion. In addition, nuclear energy generation in the postshocked matter produces up to $20\%$ more energetic explosions. Nucleosynthesis is affected due to the different dynamic evolution of the explosion. Our results indicate that the energy generation from nuclear reactions helps to sustain late outflows from the vicinity of the proto-neutron star, synthesizing more neutron-rich species. Furthermore, we show that there are systematic discrepancies between the ejecta calculated with in-situ and ex-situ reaction networks. These differences stem from the intrinsic characteristics of evolving the composition in hydrodynamic simulations or calculating it with Lagrangian tracer particles. The mass fractions of some Ca, Ti, Cr, and Fe isotopes are consistently underproduced in postprocessing calculations, leading to different nucleosynthesis paths. Our results suggest that large in-situ nuclear reaction networks are important for a realistic feedback of the energy generation, the neutrino heating, and a more accurate ejecta composition., Comment: Submitted to ApJ. Received 2022 October 20; revised 2023 May 13; accepted 2023 May 15

Aims. The Chemical Evolution of R-process Elements in Stars (CERES) project aims to provide a homogeneous analysis of a sample of metal-poor stars ([Fe/H]<-1.5). We present the stellar parameters and the chemical abundances of elements up to Zr for a sample of 52 giant stars.Methods. We relied on a sample of high signal-to-noise UVES spectra. We determined stellar parameters from Gaia photometry and parallaxes. Chemical abundances were derived using spectrum synthesis and model atmospheres.Results. We determined chemical abundances of 26 species of 18 elements: Na, Mg, Al, Si, Ca, Sc, Ti, V, Cr, Mn, Fe, Co, Ni, Cu, Zn, Sr, Y, and Zr. For several stars, we were able to measure both neutral and ionised species, including Si, Sc, Mn, and Zr. We have roughly doubled the number of measurements of Cu for stars at [Fe/H] <= -2.5. The homogeneity of the sample made it possible to highlight the presence of two Zn-rich stars ([Zn/Fe]~+0.7), one r-rich and the other r-poor. We report the existence of two branches in the [Zn/Fe] versus [Ni/Fe] plane and suggest that the high [Zn/Fe] branch is the result of hypernova nucleosynthesis. We discovered two stars with peculiar light neutron-capture abundance patterns: CES1237+1922 (also known as BS 16085-0050), which is ~1 dex underabundant in Sr, Y, and Zr with respect to the other stars in the sample, and CES2250-4057 (also known as HE 2247-4113), which shows a ~1 dex overabundance of Sr with respect to Y and Zr.Conclusions. The high quality of our dataset allowed us to measure hardly detectable ions. This can provide guidance in the development of line formation computations that take deviations from local thermodynamic equilibrium and hydrodynamical effects into account.

Magnetorotational supernovae are a rare type of core-collapse supernovae where the magnetic field and rotation play a central role in the dynamics of the explosion. We present the post-processed nucleosynthesis of state-of-the-art neutrino-MHD supernova models that follow the post explosion evolution for few seconds. We find three different dynamical mechanisms to produce heavy r-process elements: i) a prompt ejection of matter right after core bounce, ii) neutron-rich matter that is ejected at late times due to a reconfiguration of the protoneutronstar shape, iii) small amount of mass ejected with high entropies in the center of the jet. We investigate total ejecta yields, including the ones of unstable nuclei such as $^{26}$Al, $^{44}$Ti, $^{56}$Ni, and $^{60}$Fe. The obtained $^{56}$Ni masses vary between $0.01 - 1\,\mathrm{M_\odot}$. The latter maximum is compatible with hypernova observations. Furthermore, all of our models synthesize Zn masses in agreement with observations of old metal-poor stars. We calculate simplified light curves to investigate whether our models can be candidates for superluminous supernovae. The peak luminosities obtained from taking into account only nuclear heating reach up to a few $\sim 10^{43} \,\mathrm{erg\,s^{-1}}$. Under certain conditions, we find a significant impact of the $^{66}$Ni decay chain that can raise the peak luminosity up to $\sim 38\%$ compared to models including only the $^{56}$Ni decay chain. This work reinforces the theoretical evidence on the critical role of magnetorotational supernovae to understand the occurrence of hypernovae, superluminous supernovae, and the synthesis of heavy elements.

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., Comment: 96 pages. Submitted to Journal of Physics G

A promising astrophysical site to produce the lighter heavy elements of the first $r$-process peak ($Z = 38-47$) is the moderately neutron rich ($0.4 < Y_e < 0.5$) neutrino-driven ejecta of explosive environments, such as core-collapse supernovae and neutron star mergers, where the weak $r$-process operates. This nucleosynthesis exhibits uncertainties from the absence of experimental data from $(\alpha,xn)$ reactions on neutron-rich nuclei, which are currently based on statistical model estimates. In this work, we report on a new study of the nuclear reaction impact using a Monte Carlo approach and improved $(\alpha,xn)$ rates based on the Atomki-V2 $\alpha$ Optical Model Potential ($\alpha$OMP). We compare our results with observations from an up-to-date list of metal-poor stars with [Fe/H] $<$ -1.5 to find conditions of the neutrino-driven wind where the lighter heavy elements can be synthesized. We identified a list of $(\alpha,xn)$ reaction rates that affect key elemental ratios in different astrophysical conditions. Our study aims on motivating more nuclear physics experiments on $(\alpha, xn)$ reactions using current and the new generation of radioactive beam facilities and also more observational studies of metal-poor stars., Comment: 16 pages, 8 figures, submitted to ApJ

We investigate the post-explosion phase in core-collapse supernovae with 2D hydrodynamical simulations and a simple neutrino treatment. The latter allows us to perform 46 simulations and follow the evolution of the 32 successful explosions during several seconds. We present a broad study based on three progenitors (11.2 $M_\odot$, 15 $M_\odot$, and 27 $M_\odot$), different neutrino-heating efficiencies, and various rotation rates. We show that the first seconds after shock revival determine the final explosion energy, remnant mass, and properties of ejected matter. Our results suggest that a continued mass accretion increases the explosion energy even at late times. We link the late-time mass accretion to initial conditions such as rotation strength and shock deformation at explosion time. Only some of our simulations develop a neutrino-driven wind that survives for several seconds. This indicates that neutrino-driven winds are not a standard feature expected after every successful explosion. Even if our neutrino treatment is simple, we estimate the nucleosynthesis of the exploding models for the 15 $M_\odot$ progenitor after correcting the neutrino energies and luminosities to get a more realistic electron fraction., Comment: 16 pages, 15 figures, submitted to ApJ

We study the evolution of Eu and Ba abundances in local group dwarf spheroidal and ultra faint dwarf galaxies by means of detailed chemical evolution models and compare our results with new sets of homogeneous abundances. The adopted models include gas infall and outflow and have been previously tested. We investigate several production scenarios for r-process elements: merging neutron stars and magneto-rotational driven supernovae. Production of Ba through the main s-process acting in low- and intermediate- mass stars is considered as well. We also test different sets of nucleosynthesis yields. For merging neutron stars we adopt either a constant and short delay time for merging or a delay time distribution function. Our simulations show that: i) if r-process elements are produced only by a quick source, it is possible to reproduce the [Eu/Fe] vs [Fe/H], but those models fail in reproducing the [Ba/Fe] vs [Fe/H]. ii) If r-process elements are produced only with longer delays the opposite happens. iii) If both a quick source and a delayed one are adopted, such as magneto-rotational driven supernovae and merging neutron stars with a delay time distribution, the [Eu/Fe] abundance pattern is successfully reproduced, but models still fail in reproducing the [Ba/Fe]. iv) On the other hand, the characteristic abundances of Reticulum II can be reproduced only if both the Eu and the r-process fraction of Ba are produced on short and constant time delays during a single merging event. We discuss also other possible interpretations, including an inhomogeneous mixing of gas which might characterize this galaxy.

We present and discuss three extremely r-process enhanced stars located in the massive dwarf spheriodal galaxy Fornax. These stars are very unique with an extreme Eu enrichment ($1.25\le \mathrm{[Eu/Fe]} \le 1.45$) at high metallicities ($-1.3 \le \mathrm{[Fe/H]} \le -0.8$). They have the largest Eu abundances ever observed in a dwarf galaxy opening new opportunities to further understand the origin of heavy elements formed by the r-process. We derive stellar abundances of Co, Zr, La, Ce, Pr, Nd, Er, and Lu using 1-dimensional, local thermodynamic equilibrium (LTE) codes and model atmospheres in conjunction with state-of-the art yield predictions. We derive Zr in the largest sample of stars (105) known to date in a dwarf galaxy. Accurate stellar abundances combined with a careful assessment of the yield predictions have revealed three metal-rich stars in Fornax showing a pure r-process pattern. We define a new class of stars, namely Eu-stars, as r-II stars (i.e., [Eu/Fe]$>1$) at high metallicities (i.e., $\mathrm{[Fe/H]}\gtrsim -1.5$). The stellar abundance pattern contains Lu, observed for the first time in a dwarf galaxy, and reveals that a late burst of star formation has facilitated extreme r-process enhancement late in the galaxy's history ($<4$\,Gyr ago). Due to the large uncertainties associated with the nuclear physics input in the yield predictions, we cannot yet determine the r-process site leading to the three Eu-stars in Fornax. Our results demonstrate that extremely r-rich stars are not only associated with ultra faint low-mass dwarf galaxies, but can be born also in massive dwarf galaxies.

We present the nucleosynthesis of magneto-rotational supernovae (MR-SNe) including neutrino-driven and magneto-rotational-driven ejecta based, for the first time, on 2D simulations with accurate neutrino transport. The models analysed here have different rotation and magnetic fields, allowing us to explore the impact of these two key ingredients. The accurate neutrino transport of the simulations is critical to analyse the slightly neutron-rich and proton-rich ejecta that are similar to the, also neutrino-driven, ejecta in standard supernovae. In the model with strong magnetic field, the r-process produces heavy elements up to the third r-process peak ($A\sim 195$), in agreement with previous works. This model presents a jet-like explosion with proton-rich jets surrounded by neutron-rich material where the r-process occurs. We have estimated a lower limit for $^{56}$Ni of $2.5\times10^{-2} M_\odot$, which is still well below the expected hypernova value. Longer simulations including the accretion disc evolution are required to get a final prediction. In addition, we have found that the late evolution is critical in a model with weak magnetic field in which late-ejected neutron-rich matter produces elements up to the second r-process peak. Even if we cannot yet provide conclusions for hypernova nucleosynthesis, our results agree with observations of old stars and radioactive isotopes in supernova remnants. This makes MR-SNe a good additional scenario to neutron star mergers for the synthesis of heavy elements and brings us closer to understand their origin and the role of MR-SNe in the early Galaxy nucleosynthesis.

Lighter heavy elements beyond iron and up to around silver can form in neutrino-driven ejecta in core-collapse supernovae and neutron star mergers. Slightly neutron-rich conditions favour a weak r-process that follows a path close to stability. Therefore, the beta decays are slow compared to the expansion time scales, and ($\alpha$,n) reactions become critical to move matter towards heavier nuclei. The rates of these reactions are calculated with the statistical model and their main uncertainty, at energies relevant for the weak r-process, is the $\alpha$+nucleus optical potential. There are several sets of parameters to calculate the $\alpha$+nucleus optical potential leading to large deviations for the reaction rates, exceeding even one order of magnitude. Recently the $^{96}$Zr($\alpha$,n)$^{99}$Mo reaction has been identified as a key reaction that impacts the production of elements from Ru to Cd. Here, we present the first cross section measurement of this reaction at energies (6.22 MeV $\leq$ E$_\mathrm{c.m.}$ $\leq$ 12.47 MeV) relevant for the weak r-process. The new data provide a stringent test of various model predictions which is necessary to improve the precision of the weak r-process network calculations. The strongly reduced reaction rate uncertainty leads to very well-constrained nucleosynthesis yields for $Z = 44 - 48$ isotopes under different neutrino-driven wind conditions.

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

We present a large homogeneous set of stellar parameters and abundances across a broad range of metallicities, involving $13$ classical dwarf spheroidal (dSph) and ultra-faint dSph (UFD) galaxies. In total this study includes $380$ stars in Fornax, Sagittarius, Sculptor, Sextans, Carina, Ursa Minor, Draco, Reticulum II, Bootes I, Ursa Major II, Leo I, Segue I, and Triangulum II. This sample represents the largest, homogeneous, high-resolution study of dSph galaxies to date. With our homogeneously derived catalog, we are able to search for similar and deviating trends across different galaxies. We investigate the mass dependence of the individual systems on the production of $\alpha$-elements, but also try to shed light on the long-standing puzzle of the dominant production site of r-process elements. We use data from the Keck observatory archive and the ESO reduced archive to reanalyze stars from these $13$ dSph galaxies. We automatize the step of obtaining stellar parameters, but run a full spectrum synthesis to derive all abundances except for iron. The homogenized set of abundances yielded the unique possibility to derive a relation between the onset of type Ia supernovae and the stellar mass of the galaxy. Furthermore, we derived a formula to estimate the evolution of $\alpha$-elements. Placing all abundances consistently on the same scale is crucial to answer questions about the chemical history of galaxies. By homogeneously analysing Ba and Eu in the 13 systems, we have traced the onset of the s-process and found it to increase with metallicity as a function of the galaxy's stellar mass. Moreover, the r-process material correlates with the $\alpha$-elements indicating some co-production of these, which in turn would point towards rare core-collapse supernovae rather than binary neutron star mergers as host for the r-process at low [Fe/H] in the investigated dSph systems.

We present results of three-dimensional (3D), radiation-magnetohydrodynamics (MHD) simulations of core-collapse supernovae in full general relativity (GR) with spectral neutrino transport. In order to study the effects of progenitor's rotation and magnetic fields, we compute three models, where the precollapse rotation rate and magnetic fields are included parametrically to a 20 M$_{\odot}$ star. While we find no shock revival in our two non-magnetized models during our simulation times ($\sim500$ ms after bounce), the magnetorotationally (MR) driven shock expansion immediately initiates after bounce in our rapidly rotating and strongly magnetized model. We show that the expansion of the MR-driven flows toward the polar directions is predominantly driven by the magnetic pressure, whereas the shock expansion toward the equatorial direction is supported by neutrino heating. Our detailed analysis indicates that the growth of the so-called kink instability may hinder the collimation of jets, resulting in the formation of broader outflows. Furthermore we find a dipole emission of lepton number, only in the MR explosion model, whose asymmetry is consistent with the explosion morphology. Although it is similar to the lepton-number emission self-sustained asymmetry (LESA), our analysis shows that the dipole emission occurs not from the protoneutron star convection zone but from above the neutrino sphere indicating that it is not associated with the LESA. We also report several unique neutrino signatures, which are significantly dependent on both the time and the viewing angle, if observed, possibly providing a rich information regarding the onset of the MR-driven explosion., Comment: 19 pages, 14 figures, accepted for publication in the ApJ

Neutrino-driven ejecta in core collapse supernovae (CCSNe) offer an interesting astrophysical scenario where lighter heavy elements between Sr and Ag can be synthesized. Previous studies emphasized the important role that ($\alpha,n$) reactions play in the production of these elements, particularly in neutron-rich and alpha-rich environments. In this paper, we have investigated the sensitivity of elemental abundances to specific ($\alpha,n$) reaction-rate uncertainties under different astrophysical conditions. Following a Monte Carlo nucleosynthesis study with over 36 representative astrophysical wind conditions, we have identified the most important reactions based on their impact on the final elemental abundances. Experimental studies of these reactions will reduce the nucleosynthesis uncertainties and make it possible to use observations to understand the origin of lighter heavy elements and the astrophysical conditions where they are formed.

Neutrino-driven winds following core collapse supernovae have been proposed as a suitable site where the so-called light heavy elements (between Sr to Ag) can be synthetized. For moderately neutron-rich winds, ($\alpha,n$) reactions play a critical role in the weak r process, becoming the main mechanism to drive nuclear matter towards heavier elements. In this paper we summarize the sensitivity of network-calculated abundances to the astrophysical conditions, and to uncertainties in the ($\alpha,n$) reaction rates. A list of few ($\alpha,n$) reactions were identified to dominate the uncertainty in the calculated elemental abundances. Measurements of these reactions will allow to identify the astrophysical conditions of the weak r process by comparing calculated/observed abundances in r-limited stars.

Half of all the elements in the universe heavier than iron were created by rapid neutron capture. The theory for this astrophysical `$r$-process' was worked out six decades ago and requires an enormous neutron flux to make the bulk of these elements. Where this happens is still debated. A key piece of missing evidence is the identification of freshly-synthesised $r$-process elements in an astrophysical site. Current models and circumstantial evidence point to neutron star mergers as a probable $r$-process site, with the optical/infrared `kilonova' emerging in the days after the merger a likely place to detect the spectral signatures of newly-created neutron-capture elements. The kilonova, AT2017gfo, emerging from the gravitational-wave--discovered neutron star merger, GW170817, was the first kilonova where detailed spectra were recorded. When these spectra were first reported it was argued that they were broadly consonant with an outflow of radioactive heavy elements, however, there was no robust identification of any element. Here we report the identification of the neutron-capture element strontium in a re-analysis of these spectra. The detection of a neutron-capture element associated with the collision of two extreme-density stars establishes the origin of $r$-process elements in neutron star mergers, and demonstrates that neutron stars contain neutron-rich matter.

Several extremely metal-poor stars are known to have an enhanced thorium abundance. These actinide-boost stars have likely inherited material from an r-process that operated under different conditions than the most frequent r-process in the galaxy. In this article, we explore the sensitivity of actinide production in r-process calculations on the hydrodynamical conditions as well as on the nuclear physics. We find that the initial electron fraction Ye is the most important factor determining the actinide yields and that the abundance ratios between long-lived actinides and lanthanides like europium can vary for different conditions in our calculations. In our setup, conditions with high entropies systematically lead to lower actinide abundances relative to other r-process elements. Furthermore, actinide-enhanced ejecta can be distinguished from the "regular" composition also in other ways, most notably in the second r-process peak abundances., Comment: 20 pages, 14 figures, submitted to 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

Metal-poor stars in the Galactic halo often show strong enhancements in carbon and/or neutron-capture elements. However, the Galactic bulge is notable for its paucity of carbon-enhanced metal-poor (CEMP) and/or CH-stars, with only two such objects known to date. This begs the question whether the processes that produced their abundance distribution were governed by a comparable nucleosynthesis in similar stellar sites as for their more numerous counterparts in the halo. Recently, two contenders of such stars were discovered in the bulge, at [Fe/H] = $-1.5$ and $-$2.5 dex, both of which show enhancements in [C/Fe] of 0.4 and 1.4 dex, [Ba/Fe] in excess of 1.3 dex, and also elevated nitrogen. The more metal-poor of the stars is matched by standard $s$-process nucleosynthesis in low-mass Asymptotic Giant Branch (AGB) polluters. The other star shows an abnormally high [Rb/Fe] ratio. Here, we investigate the origin of the abundance peculiarities in the Rb-rich star by new, detailed measurements of heavy element abundances and by comparing the chemical element ratios of 36 species to models of neutron-capture nucleosynthesis. The $i$-process with intermediate neutron densities between those of the $s$- and $r$-processes has been previously found to provide good matches of CEMP stars with enhancements in both $r$- and $s$-process elements, rather than invoking a superposition of yields from the respective individual processes. However, the peculiar bulge star is incompatible with a pure $i$-process from a single ingestion event. Instead, it can, statistically, be better reproduced by models accounting for two proton ingestion events, or by an $i$-process component in combination with $s$-process nucleosynthesis in low-to-intermediate mass AGB stars, indicating multiple polluters. [abridged], Comment: 11 pages, 6 figures, accepted for publication in Astronomy & Astrophysics

We investigate the impact of different properties of the nuclear equation of state in core-collapse supernovae, with a focus on the proto-neutron-star contraction and its impact on the shock evolution. To this end, we introduce a range of equations of state that vary the nucleon effective mass, incompressibility, symmetry energy, and nuclear saturation point. This allows us to point to the different effects in changing these properties from the Lattimer and Swesty to the Shen et al. equations of state, the two most commonly used equations of state in simulations. In particular, we trace the contraction behavior to the effective mass, which determines the thermal nucleonic contributions to the equation of state. Larger effective masses lead to lower pressures at nuclear densities and a lower thermal index. This results in a more rapid contraction of the proto-neutron star and consequently higher neutrino energies, which aids the shock evolution to a faster explosion., Comment: 6 pages, 4 figures, minor changes, published version

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)

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

Neutrinos play a crucial role in the core-collapse supernova (CCSN) explosion mechanism. The requirement of accurately calculating the transport of neutrinos makes simulations of the CCSN mechanism extremely challenging and computationally expensive. Historically, this stiff challenge has been met by making approximations to the full transport equation. In this work, we compare CCSN simulations in one- and two-dimensions with three approximate neutrino transport schemes, each implemented in the {\tt FLASH} simulation framework. We compare a two moment M1 scheme with an analytic closure (M1), the Isotropic Diffusion Source Approximation (IDSA), and the Advanced Spectral Leakage (ASL) method. We identify and discuss the advantages and disadvantages of each scheme. For each approximate transport scheme, we use identical grid setups, hydrodynamics, and gravity solvers to investigate the transport effects on supernova shock dynamics and neutrino quantities. We find that the transport scheme has a small effect on the evolution of protoneutron star (PNS) radius, PNS mass, and the mass accretion rate. The neutrino luminosities, mean energies, and shock radii have a $\sim$ 10-20\% quantitative difference but the overall qualitative trends are fairly consistent between all three approximations. We find larger differences in the gain region properties, including the gain region mass and the net heating rate in the gain region, as well as the strength of PNS convection in the core. We investigate the progenitor, nuclear equation of state, and stochastic perturbation dependence of our simulations and find similar magnitudes of impact on key quantities. We also compare the computational expense of the various approximations., Comment: 28 pages, 11 figures, accepted by J. Phys. G focus issue on core-collapse supernovae

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., Comment: 132 pages, 18 figures, submitted to Journal of Physics G

The origin of the so-called $p$-isotopes $^{92,94}\mathrm{Mo}$ and $^{96,98}\mathrm{Ru}$ in the solar system remains a mystery as several astrophysical scenarios fail to account for them. In addition, data on presolar silicon carbide grains of type X (SiC X) exhibit peculiar Mo patterns, especially for $^{95,97}\mathrm{Mo}$. We examine production of Mo and Ru isotopes in neutrino-driven winds associated with core-collapse supernovae (CCSNe) over a wide range of conditions. We find that proton-rich winds can make dominant contributions to the solar abundance of $^{98}\mathrm{Ru}$, significant contributions to those of $^{96}$Ru ($\lesssim 40\%$) and $^{92}$Mo ($\lesssim 27\%$), and relatively minor contributions to that of $^{94}$Mo ($\lesssim 14\%$). In contrast, neutron-rich winds make negligible contributions to the solar abundances of $^{92,94}$Mo and cannot produce $^{96,98}$Ru. However, we show that some neutron-rich winds can account for the peculiar Mo patterns in SiC X grains. Our results can be generalized if conditions similar to those studied here are also obtained for other types of ejecta in either CCSNe or neutron star mergers., Comment: 12 pages, 9 figures

Core-collapse supernovae produce elements between Fe and Ag depending on the properties of the ejected matter. Despite the fast progress in supernova simulations in the last decades, there are still uncertainties in the astrophysical conditions. In this paper we investigate the impact of astrophysical uncertainties on the nucleosynthesis. Since a systematic study based on trajectories from hydrodynamic simulations is computationally very expensive, we rely on a steady-state model. By varying the mass and radius of the proto-neutron star as well as electron fraction in the steady-state model, we cover a wide range of astrophysical conditions. In our study, we find four abundance patterns which can be formed in neutron-rich neutrino-driven ejecta. This provides a unique template of trajectories that can be used to investigate the impact of nuclear physics input on the nucleosynthesis for representative astrophysical conditions. Furthermore, we link these four patterns to the neutron-to-seed and alpha-to-seed ratios at $T=3$~GK. Therefore, our results give a good overview of the potential nucleosynthesis evolution which can occur in a supernova simulation., Comment: 9 pages, 10 figures, 1 table

Weak reactions are critical for the neutron richness of the matter dynamically ejected after the merger of two neutron stars. The neutron richness, defined by the electron fraction (Ye), determines which heavy elements are produced by the r-process and thus directly impacts the kilonova light curve. In this work, we have performed a systematic and detailed post-processing study of the impact of weak reactions on the distribution of the electron fraction and of the entropy on the dynamic ejecta obtained from an equal mass neutron star binary merger simulated in full general relativity and with microscopic equation of state. Previous investigations indicated that shocks increase Ye, however our results show that shocks can also decrease Ye, depending on their thermodynamical conditions. Moreover, we have found that neutrino absorption are key and need to be considered in future simulations. We also demonstrated that the angular dependence of the neutrino luminosity and the spatial distribution of the ejecta can lead to significant difference in the electron fraction distribution. In addition to the detailed study of the Ye evolution and its dependences, we have performed nucleosynthesis calculations. They clearly point to the necessity of improving the neutrino treatment in current simulations to be able to predict the contribution of neutron star mergers to the chemical history of the universe and to reliable calculate their kilonova light curves., Comment: 33 pages, 14 figures. Submitted to Classical and Quantum Gravity, focus issue: "Rattle and shine"

When binary systems of neutron stars merge, a very small fraction of their rest mass is ejected, either dynamically or secularly. This material is neutron-rich and its nucleosynthesis could provide the astrophysical site for the production of heavy elements in the universe, together with a kilonova signal confirming neutron-star mergers as the origin of short gamma-ray bursts. We perform full general-relativistic simulations of binary neutron-star mergers employing three different nuclear-physics EOSs, considering both equal- and unequal-mass configurations, and adopting a leakage scheme to account for neutrino radiative losses. Using a combination of techniques, we carry out an extensive and systematic study of the hydrodynamical, thermodynamical, and geometrical properties of the matter ejected dynamically, employing the WinNet nuclear-reaction network to recover the relative abundances of heavy elements produced by each configurations. Among the results obtained, three are particularly important. First, we find that both the properties of the dynamical ejecta and the nucleosynthesis yields are robust against variations of the EOS and masses, and match very well the observed chemical abundances. Second, using a conservative but robust criterion for unbound matter, we find that the amount of ejected mass is $\lesssim 10^{-3}\,M_{\odot}$, hence at least one order of magnitude smaller than what normally assumed in modelling kilonova signals. Finally, using a simplified and gray-opacity model we assess the observability of the infrared kilonova emission finding, that for all binaries the luminosity peaks around $\sim1/2$ day in the $H$-band, reaching a maximum magnitude of $-13$, and decreasing rapidly after one day. These rather low luminosities make the prospects for detecting kilonovae less promising than what assumed so far., Comment: 24 pages, 18 figures, submitted to PRD