Search

We present and discuss three extremely r-process enhanced stars located in the massive dwarf spheroidal galaxy Fornax. These stars are very unique with an extreme Eu enrichment (1.25 ≤ [Eu/Fe]≤1.45) at high metallicities (−1.3 ≤ [Fe/H]≤−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 one-dimensional, local thermodynamic equilibrium 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., [Fe/H] ≳ −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 (

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

We present the nucleosynthesis of magneto-rotational supernovae (MR-SNe) including neutrino-driven and magneto-rotational-driven ejecta based, for the first time, on two-dimensional 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 ∼ 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 56Ni of 2.5 × 10−2M⊙, which is still well below the expected hypernova value. Longer simulations including the accretion disk 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 lately 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.

Context. 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. Aims. 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 α-elements, but also try to shed light on the long-standing puzzle of the dominant production site of r-process elements. Methods. We used data from the Keck observatory archive and the ESO reduced archive to reanalyze stars from these 13 classical dSph and UFD galaxies. We automatized the step of obtaining stellar parameters, but ran a full spectrum synthesis (1D, local thermal equilibrium) to derive all abundances except for iron to which we applied nonlocal thermodynamic equilibrium corrections where possible. Results. The homogenized set of abundances yielded the unique possibility of deriving 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 α-elements. This reveals a universal relation of these systems across a large range in mass. Finally, we show that between stellar masses of 2.1 × 107 M⊙ and 2.9 × 105 M⊙, there is no dependence of the production of heavy r-process elements on the stellar mass of the galaxy. Conclusions. Placing all abundances consistently on the same scale is crucial to answering questions about the chemical history of galaxies. By homogeneously analyzing 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 α-elements indicating some coproduction of these, which in turn would point toward rare core-collapse supernovae rather than binary neutron star mergers as a 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., 19 pages, 14 figures, accepted for publication in the ApJ

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, (α,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 (α,n) reaction rates. A list of few (α,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.

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

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., 28 pages, 11 figures, accepted by J. Phys. G focus issue on core-collapse supernovae

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., 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., 9 pages, 10 figures, 1 table

In this chapter we provide a description of the current state of the art in the field of electromagnetic emission and nucleosynthesis from neutron star binaries. We will discuss binaries composed by two neutron stars or by a neutron star and a black hole. While we took the effort to represent what we believe are some of the most important results in this field, we strongly recommend the interested reader to check also the references cited in the mentioned publications for a more in detail view of each topic. As the reader will see from the following sections, this field has gained considerable attention in the latest years with an increasing number of publications on the topic. For each argument we will also discuss what we believe are the current limitations that will need to be addressed in the future.

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

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"

Binary neutron star mergers are plausible progenitor candidates for short gamma-ray bursts (GRBs); however, a detailed explanation of their central engine is still lacking. The annihilation of neutrino pairs has been proposed as one of the possible powering mechanisms. We present calculations of the energy and momentum deposition operated by neutrino pair annihilation above merger remnants. Starting from the results of a detailed, three-dimensional simulation of the aftermath of a binary neutron star merger, we compute the deposition rates over a time scale comparable to the life time of the disk ($t \approx 0.4$ s), assuming a long-lived massive neutron star (MNS). We model neutrino emission using a spectral leakage scheme and compute the neutrino annihilation rates using a ray-tracing algorithm. We find that the presence of the MNS increases the energy deposition rate by a factor $\sim 2$, mainly due to the annihilation of radiation coming from the MNS with radiation coming from the disk. We compute the impact of relativistic effects and discover that, despite they can significantly change the local rate intensity, the volume-integrated results are only marginally decreased. The cumulative deposited energy, extrapolated to 1 sec, is $\approx 2.2 \times 10^{49} \, {\rm erg}$. A comparison with the inferred short GRB energetics reveals that in most cases this energy is not large enough, even assuming small jet opening angles and a long-lived MNS. Significantly more intense neutrino luminosities (a factor ~5-10 larger) are required to explain most of the observed short GRB. We conclude that it is unlikely that neutrino pair annihilation can explain the central engine of short GRBs alone., 38 pages, 19 figures

We follow the longterm evolution of the dynamic ejecta of neutron star mergers for up to 100 years and over a density range of roughly 40 orders of magnitude. We include the nuclear energy input from the freshly synthesized, radioactively decaying nuclei in our simulations and study its effects on the remnant dynamics. Although the nuclear heating substantially alters the longterm evolution, we find that running nuclear networks over purely hydrodynamic simulations (i.e. without heating) yields actually acceptable nucleosynthesis results. The main dynamic effect of the radioactive heating is to quickly smooth out inhomogeneities in the initial mass distribution, subsequently the evolution proceeds self-similarly and after 100 years the remnant still carries the memory of the initial binary mass ratio. We also explore the nucleosynthetic yields for two mass ejection channels. The dynamic ejecta very robustly produce "strong" r-process elements with $A > 130$ with a pattern that is essentially independent of the details of the merging system. From a simple model we find that neutrino-driven winds yield "weak" r-process contributions with $50 < A < 130$ whose abundance patterns vary substantially between different merger cases. This is because their electron fraction, set by the ratio of neutrino luminosities, varies considerably from case to case. Such winds do not produce any $^{56}{\rm Ni}$, but a range of radioactive isotopes that are long-lived enough to produce a second, radioactively powered electromagnetic transient in addition to the "macronova" from the dynamic ejecta. While our wind model is very simple, it nevertheless demonstrates the potential of such neutrino-driven winds for electromagnetic transients and it motivates further, more detailed neutrino-hydrodynamic studies. The properties of the mentioned transients are discussed in more detail in a companion paper., 13 pages, 11 figures, accepted to MNRAS

After a successful core-collapse supernova, a neutrino-driven wind develops where it is possible to synthesize lighter heavy elements ($30, Comment: 15 pages, 8 figures

The direct impact of white dwarfs has been suggested as a plausible channel for type Ia supernovae. In spite of their (a priori) rareness, in highly populated globular clusters and in galactic centers, where the amount of white dwarfs is considerable, the rate of violent collisions between two of them might be non-negligible. Even more, there are indications that binary white dwarf systems orbited by a third stellar-mass body have an important chance to induce a clean head-on collision. Therefore, this scenario represents a source of contamination for the supernova light-curves sample that it is used as standard candles in cosmology, and it deserves further investigation. Some groups have conducted numerical simulations of this scenario, but their results show several differences. In this paper we address some of the possible sources of these differences, presenting the results of high resolution hydrodynamical simulations jointly with a detailed nuclear post-processing of the nuclear abundances, to check the viability of white dwarf collisions to produce significant amounts of 56Ni. To that purpose, we use a 2D-axial symmetric smoothed particle hydrodynamic code to obtain a resolution considerably higher than in previous studies. In this work, we also study how the initial mass and nuclear composition affect the results. The gravitational wave emission is also calculated, as this is a unique signature of this kind of events. All calculated models produce a significant amount of 56Ni, ranging from 0.1 Msun to 1.1 Msun, compatible not only with normal-Branch type Ia supernova but also with the subluminous and super-Chandrasekhar subset. Nevertheless, the distribution mass-function of white dwarfs favors collisions among 0.6-0.7 Msun objects, leading to subluminous events., Comment: 24 pages, 12 figures, accepted for publication in MNRAS

We have explored the Eu production in the Milky Way by means of a very detailed chemical evolution model. In particular, we have assumed that Eu is formed in merging neutron star (or neutron star black hole) binaries as well as in type II supernovae. We have tested the effects of several important parameters influencing the production of Eu during the merging of two neutron stars, such as: i) the time scale of coalescence, ii) the Eu yields and iii) the range of initial masses for the progenitors of the neutron stars. The yields of Eu from type II supernovae are very uncertain, more than those from coalescing neutron stars, so we have explored several possibilities. We have compared our model results with the observed rate of coalescence of neutron stars, the solar Eu abundance, the [Eu/Fe] versus [Fe/H] relation in the solar vicinity and the [Eu/H] gradient along the Galactic disc. Our main results can be summarized as follows: i) neutron star mergers can be entirely responsible for the production of Eu in the Galaxy if the coalescence time scale is no longer than 1 Myr for the bulk of binary systems, the Eu yield is around $3 \times 10^{-7}$ M$_\odot$, and the mass range of progenitors of neutron stars is 9-50 M$_\odot$; ii) both type II supernovae and merging neutron stars can produce the right amount of Eu if the neutron star mergers produce $2 \times 10^{-7}$ M$_\odot$ per system and type II supernovae, with progenitors in the range 20-50 M$_\odot$, produce yields of Eu of the order of $10^{-8}-10^{-9}$ M$_\odot$; iii) either models with only neutron stars producing Eu or mixed ones can reproduce the observed Eu abundance gradient along the Galactic disc., 9 pages, 8 figures, accepted for publication in MNRAS

An erroneous interpretation of the hydrodynamical results led to an incorrect determination of the fallback masses in Ugliano et al. (2012), which also (on a smaller level) affects the neutron star masses provided in that paper. This problem was already addressed and corrected in the follow-up works by Ertl et al. (2015) and Sukhbold et al. (2015). Therefore, the reader is advised to use the new data of the latter two publications. In the remaining text of this Erratum we present the differences of the old and new fallback results in detail and explain the origin of the mistake in the original analysis by Ugliano et al. (2012)., 3 pages, 2 figures; submitted to The Astrophysical Journal

The comparison between observational abundance features and those obtained from nucleosynthesis predictions of stellar evolution and/or explosion simulations can scrutinize two aspects: (a) the conditions in the astrophysical production site and (b) the quality of the nuclear physics input utilized. Here we test the abundance features of r-process nucleosynthesis calculations using four different fission fragment distribution models. Furthermore, we explore the origin of a shift in the third r-process peak position in comparison with the solar r-process abundances which has been noticed in a number of merger nucleosynthesis predictions. We show that this shift occurs during the r-process freeze-out when neutron captures and β-decays compete and an (n,γ)-(γ,n) equilibrium is not maintained anymore. During this phase neutrons originate mainly from fission of material above A = 240. We also investigate the role of β-decay half-lives from recent theoretical advances, which lead either to a smaller amount of fissioning nuclei during freeze-out or a faster (and thus earlier) release of fission neutrons, which can (partially) prevent this shift and has an impact on the second and rare-earth peak as well.

In this study we explore nucleosynthesis in the dynamic ejecta of compact binary mergers. We are particularly interested in the question how sensitive the resulting abundance patterns are to the parameters of the merging system. Therefore, we systematically investigate combinations of neutron star masses in the range from 1.0 to 2.0 M? and, for completeness, we compare the results with those from two simulations of a neutron star black hole merger. The ejecta masses vary by a factor of 5 for the studied systems, but all amounts are (within the uncertainties of the merger rates) compatible with being a major source of the cosmic r-process. The ejecta undergo robust r-process nucleosynthesis which produces all the elements from the second to the third peak in close-to-solar ratios. Most strikingly, this r-process is extremely robust, and all 23 investigated binary systems yield practically identical abundance patterns. This is mainly the result of the ejecta being extremely neutron rich (Ye similar to 0.04) and the r-process path meandering along the neutron drip line so that the abundances are determined entirely by nuclear rather than astrophysical properties. While further questions related to galactic chemical evolution need to be explored in future studies, we consider this robustness together with the ease with which both the second and third peak are reproduced as strong indications that compact binary mergers are prime candidates for the sources of the observed unique heavy r-process component.

This article addresses three of the four nucleosynthesis processes involved in producing heavy nuclei beyond Fe (with a main focus on the r -process). Opposite to the fourth process (the s -process), which operates in stellar evolution during He- and C-burning, they are all related to explosive burning phases, (presumably) linked to core collapse supernova events of massive stars. The (classical) p -process is identified with explosive Ne/O-burning in outer zones of the progenitor star. It is initiated by the passage of the supernova shock wave and acts via photodisintegration reactions like a spallation process which produces neighboring (proton-rich) isotopes from pre-existing heavy nuclei. The reproduction of some of the so-called lighter p -isotopes with A 100 faces problems in this environment. The only recently discovered ν p -process is related to the innermost ejecta, the neutrino wind expelled from the hot proto-neutron star after core collapse and the supernova explosion. This neutrino wind is proton-rich in its early phase, producing nuclei up to 64Ge. Reactions with neutrinos permit to overcome decay/reaction bottlenecks for the flow beyond 64Ge, thus producing light p -isotopes, which face problems in the classical p -process scenario. The understanding of the r -process, being identified for a long time with rapid neutron captures and passing through nuclei far from stability, is still experiencing major problems. These are on the one hand related to nuclear uncertainties far from stability (masses, half-lives, fission barriers), affecting the process speed and abundance peaks. On the other hand the site is still not definitely located, yet. (i) Later, possibly neutron-rich, high entropy phases of the neutrino wind (if they materialize!) could permit its operation. (ii) Other options include the ejection of very neutron-rich neutron star-like matter, occurring possibly in neutron star mergers or core collapse supernova events with jets, related to prior stellar evolution with high rotation rates and magnetic fields. Two different environments are required for a weak and a main/strong r -process, witnessed by observations of low metallicity stars and meteoritic inclusions, which could possibly be identified with the two options listed above, i.e. the weak r -process could be related to the neutrino wind when changing from p -rich to n -rich conditions.

The most promising astrophysical sources of kHz gravitational waves (GWs) are the inspiral and merger of binary neutron star(NS)/black hole systems. Maximizing the scientific return of a GW detection will require identifying a coincident electro-magnetic (EM) counterpart. One of the most likely sources of isotropic EM emission from compact object mergers is a supernova-like transient powered by the radioactive decay of heavy elements synthesized in ejecta from the merger. We present the first calculations of the optical transients from compact object mergers that self-consistently determine the radioactive heating by means of a nuclear reaction network; using this heating rate, we model the light curve with a one dimensional Monte Carlo radiation transfer calculation. For an ejecta mass ~1e-2 M_sun[1e-3 M_sun] the resulting light curve peaks on a timescale ~ 1 day at a V-band luminosity nu L_nu ~ 3e41[1e41] ergs/s (M_V = -15[-14]); this corresponds to an effective "f" parameter ~3e-6 in the Li-Paczynski toy model. We argue that these results are relatively insensitive to uncertainties in the relevant nuclear physics and to the precise early-time dynamics and ejecta composition. Due to the rapid evolution and low luminosity of NS merger transients, EM counterpart searches triggered by GW detections will require close collaboration between the GW and astronomical communities. NS merger transients may also be detectable following a short-duration Gamma-Ray Burst or "blindly" with present or upcoming optical transient surveys. Because the emission produced by NS merger ejecta is powered by the formation of rare r-process elements, current optical transient surveys can directly constrain the unknown origin of the heaviest elements in the Universe.

We explore the effects of r-process nucleosynthesis on fall-back accretion in neutron star(NS)-NS and black hole-NS mergers, and the resulting implications for short-duration gamma-ray bursts (GRBs). Though dynamically important, the energy released during the r-process is not yet taken into account in merger simulations. We use a nuclear reaction network to calculate the heating (due to beta-decays and nuclear fission) experienced by material on the marginally-bound orbits nominally responsible for late-time fall-back. Since matter with longer orbital periods t_orb experiences lower densities, for longer periods of time, the total r-process heating rises rapidly with t_orb, such that material with t_orb > 1 seconds can become completely unbound. Thus, r-process heating fundamentally changes the canonical prediction of an uninterrupted power-law decline in the fall-back rate dM/dt at late times. When the timescale for r-process to complete is > 1 second, the heating produces a complete cut-off in fall-back accretion after ~ 1 second; if robust, this would imply that fall-back accretion cannot explain the late-time X-ray flaring observed following some short GRBs. However, for a narrow, but physically plausible, range of parameters, fall-back accretion can resume after ~ 10 s, despite having been strongly suppressed for ~ 1-10 s after the merger. This suggests the intriguing possibility that the gap observed between the prompt and extended emission in short GRBs is a manifestation of r-process heating., 7 pages; 4 figures; submitted to MNRAS

We present a nucleosynthesis study on the neutrino-driven wind after neutron star merger. Post-processing tracers from a hydrodynamical simulation, we determine nucleosynthesis yields that depend on both life time of the massive neutron star and polar angle. In comparison with the dynamic ejecta we find a complementary nucleosynthesis, which can explain the r-process pattern beyond the first r-process peak. Additionally, we predict possible forms of light curves for the detection of a kilonova.

We investigate the behavior and consequences of the reverse shock that terminates the supersonic expansion of the baryonic wind which is driven by neutrino heating off the surface of (non-magnetized) new-born neutron stars in supernova cores. To this end we perform long-time hydrodynamic simulations in spherical symmetry. In agreement with previous relativistic wind studies, we find that the neutrino-driven outflow accelerates to supersonic velocities and in case of a compact, about 1.4 solar mass (gravitational mass) neutron star with a radius of about 10 km, the wind reaches entropies of about 100 k_B per nucleon. The wind, however, is strongly influenced by the environment of the supernova core. It is decelerated and shock-heated abruptly by a termination shock that forms when the supersonic outflow collides with the slower preceding supernova ejecta. The radial position of this reverse shock varies with time and depends on the strength of the neutrino wind and the different conditions in progenitor stars with different masses and structure. Its basic properties and behavior can be understood by simple analytic considerations. We demonstrate that the entropy of matter going through the reverse shock can increase to a multiple of the asymptotic wind value. Seconds after the onset of the explosion it therefore can exceed 400 k_B per nucleon. The temperature of the shocked wind has typically dropped to about or less than 10^9 K, and density and temperature in the shock-decelerated matter continue to decrease only very slowly. Such conditions might strongly affect the important phases of supernova nucleosynthesis in a time and progenitor dependent way. (abridged), 23 pages, 13 figures; discussion of results revised and new Sect.4 added; accepted by Astronomy & Astrophysics

Heavy elements like gold, platinum or uranium are produced in the r-process, which needs neutron-rich and explosive environments. Neutron star mergers are a promising candidate for an r-process site. They exhibit three different channels for matter ejection fulfilling these conditions: dynamic ejecta due to tidal torques, neutrino-driven winds and evaporating matter from the accretion disk. We present a first study of the integrated nucleosynthesis for a neutrino-driven wind from a neutron star merger with a hyper-massive neutron star. Trajectories from a recent hydrodynamical simulation are divided into four different angle regions and post-processed with a reaction network. We find that the electron fraction varies around $Y_e \approx 0.1 - 0.4$, but its distribution differs for every angle of ejection. Hence, the wind ejecta do not undergo a robust r-process, but rather possess distinct nucleosynthesis yields depending on the angle range. Compared to the dynamic ejecta, a smaller amount of neutron-rich matter gets unbound, but the production of lighter heavy elements with $A \lesssim 130$ in the neutrino-driven wind can complement the strong r-process of the dynamic ejecta., Comment: 6 pages, 4 figures, accepted for publication on PoS

Context. Strontium has proven itself to be one of the most important neutron-capture elements in the study of metal-poor stars. Th anks to the strong absorption lines of Sr, they can be detected eve n in the most metal-poor stars and also in low-resolution spe ctra. However, we still cannot explain the large star-to-star abundance sc atter we derive for metal-poor stars. Aims. Here we compare Galactic chemical evolution (GCE) predictions with improved abundances for Sri and Srii, including updated atomic data, to evaluate possible explanations for the larg e star-to-star scatter at low metallicities. Methods. We have derived abundances under both local thermodynamic equilibrium (LTE) and non-LTE (NLTE) for stars spanning a large interval of metallicities, as well as a broad range of other stellar parameters. Gravities and metallicities are also determined in NLTE. We employed MARCS stellar atmospheres and MOOG for theLTE spectrum synthesis, while MAFAGS and DETAIL were used to derive the NLTE abundances. We verified the consisten cy of the two methods in LTE. Results. We confirm that the ionisation equilibrium between Sr I and Sr II is satisfied under NLTE but not LTE, where the di fference between neutral and ionised Sr is on average ∼ 0.3 dex. We show that the NLTE corrections are of increasing importance as the metallicity decreases. For the stars with [Fe/H] > −3, the Sri NLTE correction is∼ 0.35/0.55 dex in dwarfs/giants, while the Srii NLTE correction is

Neutrino-driven winds that follow core collapse supernovae were thought to be the site where half of the heavy elements are produced by the r-process. Although recent hydrodynamic simulations show that the conditions in the wind are not enough for the r-process, lighter heavy elements like Sr, Y, and Zr can be produced. However, it is still not clear whether the conditions in the wind are slightly neutron rich or proton rich. Here we investigate the nucleosynthesis in neutrino-driven winds for both these conditions and systematically explore the impact of wind parameters on abundances. Our results show the difficulty of obtaining a robust abundance pattern in neutron-rich winds where also an over production of Sr, Y, and Zr is likely. In proton-rich conditions, the abundances smoothly change when varying wind parameters. Constraints for wind parameters and neutrino energies and luminosities will soon become available by combining nucleosynthesis studies, like the one presented here, with new and future experimental data and observations., Comment: 18 pages, 6figures, to appear in Journal of Physics G Focus issue on "Nucleosynthesis and the role of neutrinos"

Abundances of low-metallicity stars offer a unique opportunity to understand the contribution and conditions of the different processes that synthesize heavy elements. Many old, metal-poor stars show a robust abundance pattern for elements heavier than Ba, and a less robust pattern between Sr and Ag. Here we probe if two nucleosynthesis processes are sufficient to explain the stellar abundances at low metallicity, and we carry out a site independent approach to separate the contribution from these two processes or components to the total observationally derived abundances. Our approach provides a method to determine the contribution of each process to the production of elements such as Sr, Zr, Ba, and Eu. We explore the observed star-to-star abundance scatter as a function of metallicity that each process leads to. Moreover, we use the deduced abundance pattern of one of the nucleosynthesis components to constrain the astrophysical conditions of neutrino-driven winds from core-collapse supernovae., 13 pages, published in ApJ

We present a detailed, 3D hydrodynamics study of the neutrino-driven winds that emerge from the remnant of a NS merger. Our simulations are performed with the Newtonian, Eulerian code FISH, augmented by a detailed, spectral neutrino leakage scheme that accounts for heating due to neutrino absorption in optically thin conditions. Consistent with the 2D study of Dessart et al. (2009), we find that a strong baryonic wind is blown out along the original binary rotation axis within $100$ ms after the merger. We compute a lower limit on the expelled mass of $3.5 \times 10^{-3} M_{\odot}$, large enough to be relevant for heavy element nucleosynthesis. The physical properties vary significantly between different wind regions. For example, due to stronger neutrino irradiation, the polar regions show substantially larger $Y_e$ than those at lower latitudes. This has its bearings on the nucleosynthesis: the polar ejecta produce interesting r-process contributions from $A\sim 80$ to about 130, while the more neutron-rich, lower-latitude parts produce also elements up to the third r-process peak near $A\sim 195$. We also calculate the properties of electromagnetic transients that are powered by the radioactivity in the wind, in addition to the macronova transient that stems from the dynamic ejecta. The high-latitude (polar) regions produce UV/optical transients reaching luminosities up to $10^{41} {\rm erg \, s^{-1}}$, which peak around 1 day in optical and 0.3 days in bolometric luminosity. The lower-latitude regions, due to their contamination with high-opacity heavy elements, produce dimmer and more red signals, peaking after $\sim 2$ days in optical and infrared. Our numerical experiments indicate that it will be difficult to infer the collapse time-scale of the HMNS to a BH based on the wind electromagnetic transient, at least for collapse time-scales larger than the wind production time-scale., Comment: 25 pages, 4 tables, 22 figures. Submitted to MNRAS

Core-collapse supernovae produce a hot protoneutron star that cools emitting huge amounts of neutrinos of all flavors. The interaction of these neutrinos with the outer layers of the protoneutron star produces an outflow of matter whose composition is determined by the luminosities and energies of the emitted neutrinos and antineutrinos. The presence of light nuclei like deuterons and tritons can have a big impact in the average energies of the emitted antineutrinos and consequently in the neutron-richness of the ejected matter. Recent hydrodynamical models show that the ejected matter is in fact proton-rich and constitutes the site of the ν p -process where antineutrino absorption reactions catalyze the nucleosynthesis of nuclei with A > 64 .

The origin of rapid neutron capture (r-process) nuclei remains one of the longest standing mysteries in nuclear astrophysics. Core collapse supernovae (SNe) and neutron star binary mergers are likely r-process sites, but little evidence yet exists for their in situ formation in such environments. Motivated by the advent of sensitive new or planned X-ray telescopes such as the Nuclear Spectroscopic Telescope Array (NuSTAR) and the Large Observatory for X-ray Timing (LOFT), we revisit the prospects for the detection of X-ray decay lines from r-process nuclei in young or nearby supernova remnants. For all remnants planned to be observed by NuSTAR (and several others), we conclude that r-process nuclei are detectable only if the remnant possesses a large overabundance O > 1e3 relative to the average yield per SN. Prospects are better for the next Galactic SN (assumed age of 3 years and distance of 10 kpc), for which an average r-process yield is detectable via the 10.7(9.2) keV line complexes of Os194 by LOFT at 6 sigma (5 sigma) confidence; the 27.3 keV line complex of Sb125 is detectable by NuSTAR at 2 sigma for O > 2. We also consider X-rays lines from the remnants of Galactic magnetars, motivated by the much higher r-process yields of the magneto-rotationally driven SNe predicted to birth magnetars. The ~ 3.6-3.9 keV lines of Sn126 are potentially detectable in the remnants of the magnetars 1E1547.0-5408 and 1E2259+586 by LOFT for an assumed r-process yield predicted by recent simulations. The (non-)detection of these lines can thus probe whether magnetars are indeed born with millisecond periods. Finally, we consider a blind survey of the Galactic plane with LOFT for r-process lines from the most recent binary neutron star merger remnant, concluding that a detection is unlikely without additional information on the merger location., 11 pages, 6 figures, 4 tables, submitted to MNRAS

We perform hydrodynamic supernova simulations in spherical symmetry for over 100 single stars of solar metallicity to explore the progenitor-explosion and progenitor-remnant connections established by the neutrino-driven mechanism. We use an approximative treatment of neutrino transport and replace the high-density interior of the neutron star (NS) by an inner boundary condition based on an analytic proto-NS core-cooling model, whose free parameters are chosen such that explosion energy, nickel production, and energy release by the compact remnant of progenitors around 20 solar masses are compatible with Supernova 1987A. Thus we are able to simulate the accretion phase, initiation of the explosion, subsequent neutrino-driven wind phase for 15-20 s, and the further evolution of the blast wave for hours to days until fallback is completed. Our results challenge long-standing paradigms. We find that remnant mass, launch time, and properties of the explosion depend strongly on the stellar structure and exhibit large variability even in narrow intervals of the progenitors' zero-age-main-sequence mass. While all progenitors with masses below about 15 solar masses yield NSs, black hole (BH) as well as NS formation is possible for more massive stars, where partial loss of the hydrogen envelope leads to weak reverse shocks and weak fallback. Our NS baryonic masses of ~1.2-2.0 solar masses and BH masses >6 solar masses are compatible with a possible lack of low-mass BHs in the empirical distribution. Neutrino heating accounts for SN energies between some 10^{50} erg and about 2*10^{51} erg, but can hardly explain more energetic explosions and nickel masses higher than 0.1-0.2 solar masses. These seem to require an alternative SN mechanism., 10 pages, 6 figures (7 eps files); extended version to account for referee comments and questions; accepted by Astrophys. J

We examine magnetorotationally driven supernovae as sources of $r$-process elements in the early Galaxy. On the basis of thermodynamic histories of tracer particles from a three-dimensional magnetohydrodynamical core-collapse supernova model with approximated neutrino transport, we perform nucleosynthesis calculations with and without considering the effects of neutrino absorption reactions on the electron fraction ($Y_{e}$) during post-processing. We find that the peak distribution of $Y_{e}$ in the ejecta is shifted from $\sim0.15$ to $\sim0.17$ and broadened toward higher $Y_{e}$ due to neutrino absorption. Nevertheless, in both cases the second and third peaks of the solar $r$-process element distribution can be well reproduced. The rare progenitor configuration that was used here, characterized by a high rotation rate and a large magnetic field necessary for the formation of bipolar jets, could naturally provide a site for the strong $r$-process in agreement with observations of the early galactic chemical evolution., Comment: 5 pages, 4 figures

We present nucleosynthesis studies based on hydrodynamical simulations of core-collapse supernovae and their subsequent neutrino-driven winds. Although the conditions found in these simulations are not suitable for the rapid neutron capture (r-process) to produce elements heavier than A$\sim$130, this can be solved by artificially increasing the wind entropy. In this way one can mimic the general behavior of an ejecta where the r-process occurs. We study the impact of the long-time dynamical evolution and of the nuclear physics input on the final abundances and show that different nuclear mass models lead to significant variations in the abundances. These differences can be linked to the behavior of nuclear masses far from stability. In addition, we have analyzed in detail the effect of neutron capture and beta-delayed neutron emission when matter decays back to stability. In all our studied cases, freeze out effects are larger than previously estimated and produce substantial changes in the post freeze out abundances., Comment: 10 pages, 3 figures. Talk at Nuclei in the Cosmos XI (Heidelberg, 2010), to appear in Proceedings of Science

After the initiation of the explosion of core-collapse supernovae, neutrinos emitted from the nascent neutron star drive a supersonic baryonic outflow. This neutrino-driven wind interacts with the more slowly moving, earlier supernova ejecta forming a wind termination shock (or reverse shock), which changes the local wind conditions and their evolution. Important nucleosynthesis processes (alpha-process, charged-particle reactions, r-process, and vp-process) occur or might occur in this environment. The nucleosynthesis depends on the long-time evolution of density, temperature, and expansion velocity. Here we present two-dimensional hydrodynamical simulations with an approximate description of neutrino-transport effects, which for the first time follow the post-bounce accretion, onset of the explosion, wind formation, and the wind expansion through the collision with the preceding supernova ejecta. Our results demonstrate that the anisotropic ejecta distribution has a great impact on the position of the reverse shock, the wind profile, and the long-time evolution. This suggests that hydrodynamic instabilities after core bounce and the consequential asymmetries may have important effects on the nucleosynthesis-relevant conditions in the neutrino-heated baryonic mass flow from proto-neutron stars., 13 pages, 8 figures, accepted for publication in A&A