Meteoritic analysis demonstrates that radioactive nuclei heavier than iron were present in the early Solar System. Among them, $^{129}$I and $^{247}$Cm both have a rapid neutron-capture process ($r$ process) origin and decay on the same timescale ($\simeq$ 15.6 Myr). We show that the $^{129}$I/$^{247}$Cm abundance ratio in the early Solar System (438$\pm$184) is immune to galactic evolution uncertainties and represents the first direct observational constraint for the properties of the last $r$-process event that polluted the pre-solar nebula. We investigate the physical conditions of this event using nucleosynthesis calculations and demonstrate that moderately neutron-rich ejecta can produce the observed ratio. We conclude that a dominant contribution by exceedingly neutron-rich ejecta is highly disfavoured.

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.

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 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.

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.

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.

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.

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]

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.

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.

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.

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.

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.

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.

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.

After a successful core-collapse supernova, a neutrino-driven wind develops where it is possible to synthesize lighter heavy elements ($30<Z<45$). In the early galaxy, the origin of these elements is associated to the r-process and to an additional process. Here we assume the additional process corresponds to the weak r-process or alpha-process taking place in neutrino-driven winds. Based on a trajectory obtained from hydrodynamical simulations we study the astrophysics and nuclear physics uncertainties with our main focus on the $(\alpha,n)$ reactions. These are critical to redistribute the matter and allow to further move it from light to heavy elements after nuclear statistical equilibrium freezes out. In this first sensitivity study, we vary all $(\alpha,n)$ reactions by given constant factors which are justified based on the uncertainties of the statistical model and its nuclear physics input, mainly alpha optical potentials under the conditions found in the wind. Our results show that $(\alpha,n)$ reactions are indeed very important to predict abundances. Therefore, further studies will follow to identify individual critical reactions. Since the nucleosynthesis path is close to stability, these reactions can be measured in the near future. This will uniquely allow to reduce the nuclear physics uncertainties and thus to use observations to constrain and understand the astrophysical conditions in the wind.

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.

Nuclear masses play a fundamental role in understanding how the heaviest elements in the Universe are created in the $r$-process. We predict $r$-process nucleosynthesis yields using neutron capture and photodissociation rates that are based on nuclear density functional theory. Using six Skyrme energy density functionals based on different optimization protocols, we determine for the first time systematic uncertainty bands -- related to mass modeling -- for $r$-process abundances in realistic astrophysical scenarios. We find that features of the underlying microphysics make an imprint on abundances especially in the vicinity of neutron shell closures: abundance peaks and troughs are reflected in trends of neutron separation energy. Further advances in nuclear theory and experiments, when linked to observations, will help in the understanding of astrophysical conditions in extreme $r$-process sites.

This white paper informs the nuclear astrophysics community and funding agencies about the scientific directions and priorities of the field and provides input from this community for the 2015 Nuclear Science Long Range Plan. It summarizes the outcome of the nuclear astrophysics town meeting that was held on August 21-23, 2014 in College Station at the campus of Texas A&M University in preparation of the NSAC Nuclear Science Long Range Plan. It also reflects the outcome of an earlier town meeting of the nuclear astrophysics community organized by the Joint Institute for Nuclear Astrophysics (JINA) on October 9- 10, 2012 Detroit, Michigan, with the purpose of developing a vision for nuclear astrophysics in light of the recent NRC decadal surveys in nuclear physics (NP2010) and astronomy (ASTRO2010). The white paper is furthermore informed by the town meeting of the Association of Research at University Nuclear Accelerators (ARUNA) that took place at the University of Notre Dame on June 12-13, 2014. In summary we find that nuclear astrophysics is a modern and vibrant field addressing fundamental science questions at the intersection of nuclear physics and astrophysics. These questions relate to the origin of the elements, the nuclear engines that drive life and death of stars, and the properties of dense matter. A broad range of nuclear accelerator facilities, astronomical observatories, theory efforts, and computational capabilities are needed. With the developments outlined in this white paper, answers to long standing key questions are well within reach in the coming decade.

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).

Neutrino-driven winds that follow core-collapse supernovae are an exciting astrophysical site for the production of heavy elements. Although hydrodynamical simulations show that the conditions in the wind are not extreme enough for a r-process up to uranium, neutrino-driven winds may be the astrophysical site where lighter heavy elements between Sr an Ag are produced, either by the weak r-process or by the $\nu \!p$-process. However, it is still not clear if the conditions in the wind are slightly neutron-rich or proton-rich. Therefore, we investigate the nucleosynthesis in the wind for neutron- and proton-rich conditions and systematically explore the impact of wind parameters on abundances. Here we focus on molybdenum that has raised attention because several astrophysical scenarios failed to reproduce the solar system (SoS) abundance ratio of $^{92}\mathrm{Mo}$ and $^{94}\mathrm{Mo}$. Moreover, available data of SiC X grains exhibit different isotopic ratios of $^{95}\mathrm{Mo}$ and $^{97}\mathrm{Mo}$ than in the SoS. We have investigated if neutrino-driven winds can reproduce the SoS $Y(^{92}\mathrm{Mo})/Y(^{94}\mathrm{Mo})$ and can explain the origin of the $Y(^{95}\mathrm{Mo})/Y(^{97}\mathrm{Mo})$ found in SiC X.

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.

We present a comprehensive nucleosynthesis study of the neutrino-driven wind in the aftermath of a binary neutron star merger. Our focus is the initial remnant phase when a massive central neutron star is present. Using tracers from a recent hydrodynamical simulation, we determine total masses and integrated abundances to characterize the composition of unbound matter. We find that the nucleosynthetic yields depend sensitively on both the life time of the massive neutron star and the polar angle. Matter in excess of up to $9 \cdot 10^{-3} M_\odot$ becomes unbound until $\sim 200~{\rm ms}$. Due to electron fractions of $Y_{\rm e} \approx 0.2 - 0.4$ mainly nuclei with mass numbers $A < 130$ are synthesized, complementing the yields from the earlier dynamic ejecta. Mixing scenarios with these two types of ejecta can explain the abundance pattern in r-process enriched metal-poor stars. Additionally, we calculate heating rates for the decay of the freshly produced radioactive isotopes. The resulting light curve peaks in the blue band after about $4~{\rm h}$. Furthermore, high opacities due to heavy r-process nuclei in the dynamic ejecta lead to a second peak in the infrared after $3-4~{\rm d}$.

Comparing observational abundance features with nucleosynthesis predictions of stellar evolution 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. We test the abundance features of r-process nucleosynthesis calculations for the dynamical ejecta of neutron star merger simulations based on three different nuclear mass models: The Finite Range Droplet Model (FRDM), the (quenched version of the) Extended Thomas Fermi Model with Strutinsky Integral (ETFSI-Q), and the Hartree-Fock-Bogoliubov (HFB) mass model. We make use of corresponding fission barrier heights and compare the impact of four different fission fragment distribution models on the final r-process abundance distribution. In particular, we explore the abundance distribution in the second r-process peak and the rare-earth sub-peak as a function of mass models and fission fragment distributions, as well as the origin of a shift in the third r-process peak position. The latter has been noticed in a number of merger nucleosynthesis predictions. We show that the shift occurs during the r-process freeze-out when neutron captures and {\beta}-decays compete and an (n,{\gamma})-({\gamma},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 {\beta}-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.

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.

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.

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.

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.

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.

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.

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.

Strontium has proven itself to be one of the most important neutron-capture elements in the study of metal-poor stars. Thanks to the strong absorption lines of Sr, they can be detected even in the most metal-poor stars and also in low-resolution spectra. However, we still cannot explain the large star-to-star abundance scatter we derive for metal-poor stars. Here we contrast Galactic chemical evolution (GCE) with improved abundances for SrI+II including updated atomic data, to evaluate possible explanations for the large star-to-star scatter at low metallicities. We derive abundances under both local thermodynamic equilibrium (LTE) and non-LTE (NLTE) for stars spanning a large interval of stellar parameters. Gravities and metallicities are also determined in NLTE. We confirm that the ionisation equilibrium between SrI and SrII is satisfied under NLTE but not LTE, where the difference between SrI and SrII is on average ~0.3dex. 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.55dex in dwarfs/giants, while the Sr II NLTE correction is +/-0.05dex. On the basis of the large NLTE corrections, SrI should not be applied as a chemical tracer under LTE, while it is a good tracer under NLTE. SrII is a good tracer under both LTE and NLTE (down to [Fe/H]\sim -3), and LTE is a safe assumption for this majority species. However, the Sr abundance from SrII lines is dependent on an accurate surface gravity determination, which can be obtained from NLTE spectroscopy of Fe lines or from parallax measurements. We could not explain the star-to-star scatter (which remains under both LTE and NLTE) by the use of the GCE model, since the Sr yields to date are too uncertain to draw firm conclusions. At least two production sites seem necessary in order to account for this large scatter (abridged).

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.

In this study we explore the 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 \Msun 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 five for the studied systems, but all amounts are (within the uncertainties of the merger rates) compatible with being a major source of cosmic r-process. The ejecta undergo a 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, all 23 investigated binary systems yield practically identical abundance patterns. This is mainly the result of the ejecta being extremely neutron rich (\ye $\approx0.04$) and the r-process path meandering along the neutron drip line so that the abundances are determined entirely by nuclear rather than by astrophysical properties. This robustness together with the ease with which both the second and third peak are reproduced make compact binary mergers the prime candidate for the source of the observed unique heavy r-process component.

Neutrino-driven winds, which follow core-collapse supernova explosions, present a fascinating nuclear astrophysics problem that requires understanding advanced astrophysics simulations, the properties of matter and neutrino interactions under extreme conditions, the structure and reactions of exotic nuclei, and comparisons against forefront astronomical observations. The neutrino-driven wind has attracted vast attention over the last 20 years as it was suggested to be a candidate for the astrophysics site where half of the heavy elements are produced via the r-process. In this review, we summarize our present understanding of neutrino-driven winds from the dynamical and nucleosynthesis perspectives. Rapid progress has been made during recent years in understanding the wind with improved simulations and better micro physics. The current status of the fields is that hydrodynamical simulations do not reach the extreme conditions necessary for the r-process and the proton or neutron richness of the wind remains to be investigated in more detail. However, nucleosynthesis studies and observations point already to neutrino-driven winds to explain the origin of lighter heavy elements, such as Sr, Y, Zr.

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.

We study the impact of the late time dynamical evolution of ejecta from core-collapse supernovae on \nu p-process nucleosynthesis. Our results are based on hydrodynamical simulations of neutrino wind ejecta. Motivated by recent two-dimensional wind simulations, we vary the dynamical evolution during the \nu p-process and show that final abundances strongly depend on the temperature evolution. When the expansion is very fast, there is not enough time for antineutrino absorption on protons to produce enough neutrons to overcome the \beta-decay waiting points and no heavy elements beyond A=64 are produced. The wind termination shock or reverse shock dramatically reduces the expansion speed of the ejecta. This extends the period during which matter remains at relatively high temperatures and is exposed to high neutrino fluxes, thus allowing for further (p,\gamma) and (n,p) reactions to occur and to synthesize elements beyond iron. We find that the \nu p-process starts to efficiently produce heavy elements only when the temperature drops below ~3 GK. At higher temperatures, due to the low alpha separation energy of 60Zn (S_{\alpha} = 2.7 MeV) the reaction 59Cu(p,\alpha)56Ni is faster than the reaction 59Cu(p,\gamma)60Zn. This results in the closed NiCu cycle that we identify and discuss here for the first time. We also investigate the late phase of the \nu p-process when the temperatures become too low to maintain proton captures. Depending on the late neutron density, the evolution to stability is dominated by \beta decays or by (n,\gamma) reactions. In the latter case, the matter flow can even reach the neutron-rich side of stability and the isotopic composition of a given element is then dominated by neutron-rich isotopes.

We show that long-range correlations for nuclear masses have a significant effect on the synthesis of heavy elements by the r-process. As calculated by Delaroche et al. [1], these correlations suppress magic number effects associated with minor shells. This impacts the calculated abundances before the third r-process peak (at mass number A~195), where the abundances are low and form a trough. This trough and the position of the third abundance peak are strongly affected by the masses of nuclei in the transition region between deformed and spherical. Based on different astrophysical environments, our results demonstrate that a microscopic theory of nuclear masses including correlations naturally smoothens the separation energies, thus reducing the trough and improving the agreement with observed solar system abundances.

We present first comparisons between Light Element Primary Process (LEPP) abundances observed in some ultra metal poor (UMP) stars and nucleosynthesis calculations based on long-time hydrodynamical simulations of core-collapse supernovae and their neutrino-driven wind. UMP star observations indicate Z>38 elements include the contributions of at least two nucleosynthesis components: r-process nuclei that are synthesized by rapid neutron capture in a yet unknown site and LEPP elements (mainly Sr, Y, Zr). We show that neutrino-driven wind simulations can explain the observed LEPP pattern. We explore in detail the sensitivity of the calculated abundances to the electron fraction, which is a key nucleosynthesis parameter but poorly known due to uncertainties in neutrino interactions and transport. Our results show that the observed LEPP pattern can be reproduced in proton- and neutron-rich winds.

We use results from long-time core-collapse supernovae simulations to investigate the impact of the late time evolution of the ejecta and of the nuclear physics input on the calculated r-process abundances. Based on the latest hydrodynamical simulations, heavy r-process elements cannot be synthesized in the neutrino-driven winds that follow the supernova explosion. However, by artificially increasing the wind entropy, elements up to A=195 can be made. In this way one can reproduce the typical behavior of high-entropy ejecta where the r-process is expected to occur. We identify which nuclear physics input is more important depending on the dynamical evolution of the ejecta. When the evolution proceeds at high temperatures (hot r-process), an (n,g)-(g,n) equilibrium is reached. While at low temperature (cold r-process) there is a competition between neutron captures and beta decays. In the first phase of the r-process, while enough neutrons are available, the most relevant nuclear physics input are the nuclear masses for the hot r-process and the neutron capture and beta-decay rates for the cold r-process. At the end of this phase, the abundances follow a steady beta flow for the hot r-process and a steady flow of neutron captures and beta decays for the cold r-process. After neutrons are almost exhausted, matter decays to stability and our results show that in both cases neutron captures are key for determining the final abundances, the position of the r-process peaks, and the formation of the rare-earth peak. In all the cases studied, we find that the freeze out occurs in a timescale of several seconds.

The specific mechanism and astrophysical site for the production of half of the elements heavier than iron via rapid neutron capture (r-process) remains to be found. In order to reproduce the abundances of the solar system and of the old halo stars, at least two components are required: the heavy r-process nuclei (A>130) and the weak r-process which correspond to the lighter heavy nuclei (A<130). In this work, we present nucleosynthesis studies based on trajectories of hydrodynamical simulations for core-collapse supernovae and their subsequent neutrino-driven winds. We show that the weak r-process elements can be produced in neutrino-driven winds and we relate their abundances to the neutrino emission from the nascent neutron star. Based on the latest hydrodynamical simulations, heavy r-process elements cannot be synthesized in the neutrino-driven winds. However, by artificially increasing the wind entropy, elements up to A=195 can be made. In this way one can mimic the general behavior of an ejecta where the r-process occurs. We use this to study the impact of the nuclear physics input (nuclear masses, neutron capture cross sections, and beta-delayed neutron emission) and of the long-time dynamical evolution on the final abundances.

The electron-to-nucleon ratio or electron fraction is a key parameter in many astrophysical studies. Its value is determined by weak-interaction rates that are based on theoretical calculations subject to several nuclear physics uncertainties. Consequently, it is important to have a model independent way of constraining the electron fraction value in different astrophysical environments. Here we show that nuclear statistical equilibrium combined with beta equilibrium can provide such a constraint. We test the validity of this approximation in presupernova models and give lower limits for the electron fraction in type Ia supernova and accretion-induced collapse.

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.

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.

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 have performed hydrodynamical simulations of the long-time evolution of proto-neutron stars to study the nucleosynthesis using the resulting wind trajectories. Although the conditions found in the present wind models are not favourable for the production of heavy elements, a small enhancement of the entropy results in the production of r-process elements with A $\approx$ 195. This allows us to explore the sensitivity of their production to the hydrodynamical evolution (wind termination shock) and nuclear physics input used.

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.

The r-process constitutes one of the major challenges in nuclear astrophysics. Its astrophysical site has not yet been identified but there is observational evidence suggesting that at least two possible sites should contribute to the solar system abundance of r-process elements and that the r-process responsible for the production of elements heavier than Z=56 operates quite robustly producing always the same relative abundances. From the nuclear-physics point of view the r-process requires the knowledge of a large number of reaction rates involving exotic nuclei. These include neutron capture rates, beta-decays and fission rates, the latter for the heavier nuclei produced in the r-process. We have developed for the first time a complete database of reaction rates that in addition to neutron-capture rates and beta-decay half-lives includes all possible reactions that can induce fission (neutron-capture, beta-decay and spontaneous fission) and the corresponding fission yields. In addition, we have implemented these reaction rates in a fully implicit reaction network. We have performed r-process calculations for the neutrino-driven wind scenario to explore whether or not fission can contribute to provide a robust r-process pattern.

We study the composition of the outer layers of a protoneutron star and show that light nuclei are present in substantial amounts. The composition is dominated by nucleons, deuterons, tritons and alpha particles; 3He is present in smaller amounts. This composition can be studied in laboratory experiments with new neutron-rich radioactive beams that can reproduce similar densities and temperatures. After including the corresponding neutrino interactions, we demonstrate that light nuclei have a small impact on the average energy of the emitted electron neutrinos, but are significant for the average energy of antineutrinos. During the early post-explosion phase, the average energy of electron antineutrinos is slightly increased, while at later times during the protoneutron star cooling it is reduced by about 1 MeV. The consequences of these changes for nucleosynthesis in neutrino-driven supernova outflows are discussed.

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)

The neutrino-driven wind from a nascent neutron star at the center of a supernova expands into the earlier ejecta of the explosion. Upon collision with this slower matter the wind material is decelerated in a wind termination shock. By means of hydrodynamic simulations in spherical symmetry we demonstrate that this can lead to a large increase of the wind entropy, density, and temperature, and to a strong deceleration of the wind expansion. The consequences of this phenomenon for the possible r-process nucleosynthesis in the late wind still need to be explored in detail. Two-dimensional models show that the wind-ejecta collision is highly anisotropic and could lead to a directional dependence of the nucleosynthesis even if the neutrino-driven wind itself is spherically symmetric.