Supernova explosions of massive stars are one of the primary sites for the production of the elements in the universe. Up to now, stars with zero-age main-sequence masses in the range of 35--50~$M_\odot$ had mostly been representing the failed supernova explosion branch. In contrast, it has been demonstrated recently that the appearance of exotic phases of hot and dense matter, associated with a sufficiently strong phase transition from nuclear matter to the quark-gluon plasma at high baryon density, can trigger supernova explosions of such massive supergiant. Here, we present the first results obtained from an extensive nucleosynthesis analysis for material being ejected from the surface of the newly born proto-neutron star of such supernova explosions. These ejecta contain an early neutron-rich component and a late-time high-entropy neutrino-driven wind. The nucleosynthesis robustly overcomes the production of nuclei associated with the second $r$-process peak, at nuclear mass number $A\simeq 130$, and proceeds beyond the formation of the third peak ($A\simeq 195$) to the actinides. These yields may account for metal-poor star observations concerning $r$-process elements such as strontium and europium in the Galaxy at low metalicity, while the actinide yields suggests that this source may be a candidate contributing to the abundances of radioactive $^{244}$Pu measured in deep-sea sediments on Earth.

The production of about half of the heavy elements found in nature is assigned to a specific astrophysical nucleosynthesis process: the rapid neutron capture process (r-process). Although this idea has been postulated more than six decades ago, the full understanding faces two types of uncertainties/open questions: (a) The nucleosynthesis path in the nuclear chart runs close to the neutron-drip line, where presently only limited experimental information is available, and one has to rely strongly on theoretical predictions for nuclear properties. (b) While for many years the occurrence of the r-process has been associated with supernovae, more recent studies have cast substantial doubts on this environment. Alternative scenarios include the mergers of neutron stars, neutron-star black hole mergers, but possibly also rare classes of supernovae as well as hypernovae/collapsars with polar jet ejecta and also accretion disk outflows related to the collapse of fast rotating massive stars with high magnetic fields. Stellar r-process abundance observations, have provided insights into, and constraints on the frequency of and conditions in the responsible stellar production sites. One of them, neutron star mergers, was just identified and related to the Gravitational Wave event GW170817. High resolution observations, increasingly more precise due to improved experimental atomic data, have been particularly important in defining the heavy element abundance patterns of the old halo stars, and thus determining the extent, and nature, of the earliest nucleosynthesis in our Galaxy. Combining new results and important breakthroughs in the related nuclear, atomic and astronomical fields of science, this review attempts to provide an answer to the question "How Were the Elements from Iron to Uranium Made?" (Abridged)

In this fourth paper of the series, we use the parametrized, spherically symmetric explosion method PUSH to perform a systematic study of two sets of non-rotating stellar progenitor models. Our study includes pre-explosion models with metallicities Z=0 and Z=Z$_{\odot}\times 10^{-4}$ and covers a progenitor mass range from 11 up to 75 M$_\odot$. We present and discuss the explosion properties of all models and predict remnant (neutron star or black hole) mass distributions within this approach. We also perform systematic nucleosynthesis studies and predict detailed isotopic yields as function of the progenitor mass and metallicity. We present a comparison of our nucleosynthesis results with observationally derived $^{56}$Ni ejecta from normal core-collapse supernovae and with iron-group abundances for metal-poor star HD~84937. Overall, our results for explosion energies, remnant mass distribution, $^{56}$Ni mass, and iron group yields are consistent with observations of normal CCSNe. We find that stellar progenitors at low and zero metallicity are more prone to BH formation than those at solar metallicity, which allows for the formation of BHs in the mass range observed by LIGO/VIRGO.

In a previously presented proof-of-principle study we established a parametrized spherically symmetric explosion method (PUSH) that can reproduce many features of core-collapse supernovae. The present paper goes beyond a specific application that is able to reproduce observational properties of SN1987A and performs a systematic study of the explosion properties for an extensive set of non-rotating, solar metallicity stellar progenitor models in the mass range from 10.8 to 120 M$_\odot$.This includes the transition from neutron stars to black holes as the final result of the collapse of massive stars, and the relation of the latter to supernovae and faint/failed supernovae. The present paper provides the basis for extended nucleosynthesis predictions in a forthcoming paper to be employed in galactic evolution models.

Rapid neutron capture process (r-process) elements have been detected in a large fraction of metal-poor halo stars, with abundances relative to iron (Fe) that vary by over two orders of magnitude. This scatter is reduced to less than a factor of 3 in younger Galactic disc stars. The large scatter of r-process elements in the early Galaxy suggests that the r-process is made by rare events, like compact binary mergers and rare sub-classes of supernovae. Although being rare, neutron star mergers alone have difficulties to explain the observed enhancement of r-process elements in the lowest metallicity stars compared to Fe. The supernovae producing the two neutron stars already provide a substantial Fe abundance where the r-process ejecta from the merger would be injected. In this work we investigate another complementary scenario, where the r-process occurs in neutron star-black hole mergers in addition to neutron star mergers. Neutron star-black hole mergers would eject similar amounts of r-process matter as neutron star mergers, but only the neutron star progenitor would have produced Fe. Furthermore, a reduced efficiency of Fe production from single stars significantly alters the age-metallicity relation, which shifts the onset of r-process production to lower metallicities. We use the high-resolution [(20 pc)3/cell] inhomogeneous chemical evolution tool `ICE' to study the outcomes of these effects. In our simulations, an adequate combination of neutron star mergers and neutron star-black hole mergers qualitatively reproduces the observed r-process abundances in the Galaxy.

We present new observational data for the heavy elements molybdenum (Mo, Z = 42) and ruthenium (Ru, Z = 44) in F-, G-, and K-stars belonging to different substructures of the Milky Way. The range of metallicity covered is --1.0 $<$ [Fe/H] $<$ +0.3. The spectra of Galactic disc stars have a high resolution of 42,000 and 75,000 and signal-to-noise ratio better than 100. Mo and Ru abundances were derived by comparing the observed and synthetic spectra in the region of Mo I lines at 5506, 5533 \AA~ for 209 stars and Ru I lines at 4080, 4584, 4757 \AA~ for 162 stars using the LTE approach. For all the stars, the Mo and Ru abundance determinations are obtained for the first time with an average error of 0.14 dex. This is the first extended sample of stellar observations for Mo and Ru in the Milky Way disc, and together with earlier observations in halo stars it is pivotal in providing a complete picture of the evolution of Mo and Ru across cosmic timescales. The Mo and Ru abundances were compared with those of the neutron-capture elements (Sr, Y, Zr, Ba, Sm, Eu). The complex nucleosynthesis history of Mo and Ru is compared with different Galactic Chemical Evolution (GCE) simulations. In general, present theoretical GCE simulations show underproduction of Mo and Ru at all metallicities compared to observations. This highlights a significant contribution of nucleosynthesis processes not yet considered in our simulations. A number of possible scenarios are discussed.

The recent detection of 60Fe in the cosmic rays provides conclusive evidence that there is a recently synthesized component (few MY) in the GCRs (Binns et al. 2016). In addition, these nuclei must have been synthesized and accelerated in supernovae near the solar system, probably in the Sco-Cen OB association subgroups, which are about 100 pc distant from the Sun. Recent theoretical work on the production of r-process nuclei appears to indicate that it is difficult for SNe to produce the solar system abundances relative to iron of r-process elements with high atomic number (Z), including the actinides (Th, U, Np, Pu, and Cm). Instead, it is believed by many that the heaviest r-process nuclei, or perhaps even all r-process nuclei, are produced in binary neutron star mergers. Since we now know that there is at least a component of the GCRs that has been recently synthesized and accelerated, models of r-process production by SNe and BNSM can be tested by measuring the relative abundances of these ultra-heavy r-process nuclei, and especially the actinides, since they are radioactive and provide clocks that give the time interval from nucleosynthesis to detection at Earth. Since BNSM are believed to be much less frequent in our galaxy than SNe (roughly 1000 times less frequent, the ratios of the actinides, each with their own half-life, will enable a clear determination of whether the heaviest r-process nuclei are synthesized in SNe or in BNSM. In addition, the r-process nuclei for the charge range from 34 to 82 can be used to constrain models of r-process production in BNSM and SNe. Thus, GCRs become a multi-messenger component in the study of BNSM and SNe.

There is now solid experimental evidence of at least one supernova explosion within 100 pc of Earth within the last few million years, from measurements of the short-lived isotope 60Fe in widespread deep-ocean samples, as well as in the lunar regolith and cosmic rays. This is the first established example of a specific dated astrophysical event outside the Solar System having a measurable impact on the Earth, offering new probes of stellar evolution, nuclear astrophysics, the astrophysics of the solar neighborhood, cosmic-ray sources and acceleration, multi-messenger astronomy, and astrobiology. Interdisciplinary connections reach broadly to include heliophysics, geology, and evolutionary biology. Objectives for the future include pinning down the nature and location of the established near-Earth supernova explosions, seeking evidence for others, and searching for other short-lived isotopes such as 26Al and 244Pu. The unique information provided by geological and lunar detections of radioactive 60Fe to assess nearby supernova explosions make now a compelling time for the astronomy community to advocate for supporting multi-disciplinary, cross-cutting research programs.

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.

The enrichment history of heavy neutron-capture elements in the Milky Way disc provides fundamental information about the chemical evolution of our Galaxy and about the stellar sources that made those elements. In this work we give new observational data for Sr, the element at the first neutron-shell closure beyond iron, N=50, based on the analysis of the high resolution spectra of 276 Galactic disc stars. The Sr abundance was derived by comparing the observed and synthetic spectra in the region of the SrI 4607 A line, making use of the LTE approximation. NLTE corrections lead to an increase of the abundance estimates obtained under LTE, but for these lines they are minor near solar metallicity. The average correction that we find is 0.151 dex. The star that is mostly affected is HD 6582, with a 0.244 dex correction. The behavior of the Sr abundance as a function of metallicity is discussed within a stellar nucleosynthesis context, in comparison with the abundance of the heavy neutron-capture elements Ba (Z=56) and Eu (Z=63). The comparison of the observational data with the current GCE models confirm that the s-process contributions from Asymptotic Giant Branch stars and from massive stars are the main sources of Sr in the Galactic disc and in the Sun, while different nucleosynthesis sources can explain the high [Sr/Ba] and [Sr/Eu] ratios observed in the early Galaxy.

Modeling core-collapse supernovae (CCSNe) with neutrino transport in three dimensions (3D) requires tremendous computing resources and some level of approximation. We present a first comparison study of CCSNe in 3D with different physics approximations and hydrodynamics codes. We aim to assess the impact of the hydrodynamics code, approximations for the neutrino and gravity treatments, and rotation on the simulation of CCSNe in 3D. We use four different hydrodynamics codes in this work (ELEPHANT, FLASH, fGR1, and SPHYNX) in combination with two different neutrino treatments, the isotropic diffusion source approximation (IDSA) and two-moment M1, and three different gravity treatments: Newtonian, 1D General Relativity (GR) correction, and full GR). Additional parameters discussed in this study are the inclusion of neutrino-electron scattering via a parametrized deleptonization (PD) and the influence of rotation. The four codes compared in this work include Eulerian and fully Lagrangian (smoothed particle hydrodynamics) codes for the first time. They show agreement in the overall evolution of the collapse phase and early post-bounce within the range of 10% (20% in some cases). The comparison of the different neutrino treatments highlights the need to further investigate the antineutrino luminosities in IDSA, which tend to be relatively high. We also demonstrate the requirement for a more detailed heavy-lepton neutrino leakage. When comparing with a full GR code, including an M1 transport method, we confirm the influence of neutrino-electron scattering during the collapse phase, which is adequately captured by the PD scheme. Also, the effective GR potential reproduces the overall dynamic evolution correctly in all Newtonian codes. Additionally, we verify that rotation aids the shock expansion and estimate the overall angular momentum losses for each code in rotating scenarios.

In a previously presented proof-of-principle study, we established a parametrized spherically symmetric explosion method (PUSH) that can reproduce many features of core-collapse supernovae for a wide range of pre-explosion models. The method is based on the neutrino-driven mechanism and follows collapse, bounce and explosion. There are two crucial aspects of our model for nucleosynthesis predictions. First, the mass cut and explosion energy emerge simultaneously from the simulation (determining, for each stellar model, the amount of Fe-group ejecta). Second, the interactions between neutrinos and matter are included consistently (setting the electron fraction of the innermost ejecta). In the present paper, we use the successful explosion models from Ebinger et al. (2018) which include two sets of pre-explosion models at solar metallicity, with combined masses between 10.8 and 120 M$_{\odot}$. We perform systematic nucleosynthesis studies and predict detailed isotopic yields. The resulting $^{56}$Ni ejecta are in overall agreement with observationally derived values from normal core-collapse supernovae. The Fe-group yields are also in agreement with derived abundances for metal-poor star HD84937. We also present a comparison of our results with observational trends in alpha element to iron ratios.

In this White Paper we present the potential of the enhanced X-ray Timing and Polarimetry (eXTP) mission for studies related to Observatory Science targets. These include flaring stars, supernova remnants, accreting white dwarfs, low and high mass X-ray binaries, radio quiet and radio loud active galactic nuclei, tidal disruption events, and gamma-ray bursts. eXTP will be excellently suited to study one common aspect of these objects: their often transient nature. Developed by an international Consortium led by the Institute of High Energy Physics of the Chinese Academy of Science, the eXTP mission is expected to be launched in the mid 2020s.

This review touches on historical aspects, going back to the early days of nuclear astrophysics, initiated by B$^2$FH and Cameron, discusses (i) the required nuclear input from reaction rates and decay properties up to the nuclear equation of state, continues (ii) with the tools to perform nucleosynthesis calculations and (iii) early parametrized nucleosynthesis studies, before (iv) reliable stellar models became available for the late stages of stellar evolution. It passes then through (v) explosive environments from core-collapse supernovae to explosive events in binary systems (including type Ia supernovae and compact binary mergers), and finally (vi) discusses the role of all these nucleosynthesis production sites in the evolution of galaxies. The focus is put on the comparison of early ideas and present, very recent, understanding.

Female role models reduce the impact on women of stereotype threat, i.e., of being at risk of conforming to a negative stereotype about one's social, gender, or racial group. This can lead women scientists to underperform or to leave their scientific career because of negative stereotypes such as, not being as talented or as interested in science as men. Sadly, history rarely provides role models for women scientists; instead, it often renders these women invisible. In response to this situation, we present a selection of twelve outstanding women who helped to develop nuclear astrophysics.

We investigate axisymmetric black hole~(BH) formation and its gravitational wave (GW) and neutrino signals with self-consistent core-collapse supernova simulations of a non-rotating $40~M_\odot$ progenitor star using the isotropic diffusion source approximation for the neutrino transport and a modified gravitational potential for general relativistic effects. We consider four different neutron star (NS) equations of state~(EoS): LS220, SFHo, BHB$\Lambda\phi$ and DD2, and study the impact of the EoS on BH formation dynamics and GW emission. We find that the BH formation time is sensitive to the EoS from 460 to $>1300$~ms and is delayed in multiple dimensions for $\sim~100-250$~ms due to the finite entropy effects. Depending on the EoS, our simulations show the possibility that shock revival can occur along with the collapse of the proto-neutron star~(PNS) to a BH. The gravitational waveforms contain four major features that are similar to previous studies but show extreme values: (1)~a low frequency signal ($\sim~300-500$~Hz) from core-bounce and prompt convection, (2)~a strong signal from the PNS g-mode oscillation among other features, (3)~a high frequency signal from the PNS inner-core convection, and (4)~signals from the standing accretion shock instability and convection. The peak frequency at the onset of BH formation reaches to $\sim~2.3$~kHz. The characteristic amplitude of a 10~kpc object at peak frequency is detectable but close to the noise threshold of the Advanced~LIGO and KAGRA, suggesting that the next generation gravitational wave detector will need to improve the sensitivity at the kHz domain to better observe stellar-mass BH formation from core-collapse supernovae or failed supernovae.

Using isochronous mass spectrometry at the experimental storage ring CSRe in Lanzhou, the masses of $^{82}$Zr and $^{84}$Nb were measured for the first time with an uncertainty of $\sim 10$ keV, and the masses of $^{79}$Y, $^{81}$Zr, and $^{83}$Nb were re-determined with a higher precision. %The latter differ significantly from their literature values. The latter are significantly less bound than their literature values. Our new and accurate masses remove the irregularities of the mass surface in this region of the nuclear chart. Our results do not support the predicted island of pronounced low $\alpha$ separation energies for neutron-deficient Mo and Tc isotopes, making the formation of Zr-Nb cycle in the $rp$-process unlikely. The new proton separation energy of $^{83}$Nb was determined to be 490(400)~keV smaller than that in the Atomic Mass Evaluation 2012. This partly removes the overproduction of the $p$-nucleus $^{84}$Sr relative to the neutron-deficient molybdenum isotopes in the previous $\nu p$-process simulations.

We study the final fate of a very massive star by performing full general relativistic (GR), three-dimensional (3D) simulation with three-flavor multi-energy neutrino transport. Utilizing a 70 solar mass zero metallicity progenitor, we self-consistently follow the radiation-hydrodynamics from the onset of gravitational core-collapse until the second collapse of the proto-neutron star (PNS), leading to black hole (BH) formation. Our results show that the BH formation occurs at a post-bounce time of ~300 ms for the 70 Msun star. This is significantly earlier than those in the literature where lower mass progenitors were employed. At a few ~10 ms before BH formation, we find that the stalled bounce shock is revived by intense neutrino heating from the very hot PNS, which is aided by violent convection behind the shock. In the context of 3D-GR core-collapse modeling with multi-energy neutrino transport, our numerical results present the first evidence to validate a fallback BH formation scenario of the 70 Msun star.

The existence of neutron star mergers has been supported since the discovery of the binary pulsar and the observation of its orbital energy loss, consistent with General Relativity. They are considered nucleosynthesis sites of the rapid neutron-capture process (r-process), which is responsible for creating approximately half of all heavy elements beyond Fe and is the only source of elements beyond Pb and Bi. Detailed nucleosynthesis calculations based on the decompression of neutron star matter are consistent with solar r-process abundances of heavy nuclei. Neutron star mergers have also been identified with short-duration {\gamma}-ray bursts via their IR afterglow. The high neutron densities in ejected matter permit a violent r-process, leading to fission cycling of the heaviest nuclei in regions far from (nuclear) stability. Uncertainties in several nuclear properties affect the abundance distributions. The modeling of astrophysical events also depends on the hydrodynamic treatment, the occurrence of a neutrino wind after the merger and before the possible emergence of a black hole, and the properties of black hole accretion disks. We discuss the effect of nuclear and modeling uncertainties and conclude that binary compact mergers are probably a (or the) dominant site of the production of r-process nuclei in our Galaxy.

Core-collapse supernovae are the first polluters of heavy elements in the galactic history. As such, it is important to study the nuclear compositions of their ejecta, and understand their dependence on the progenitor structure (e.g., mass, compactness, metallicity). Here, we present a detailed nucleosynthesis study based on two long-term, two-dimensional core-collapse supernova simulations of a 11.2 M$_{\odot}$ and a 17.0 M$_{\odot}$ star. We find that in both models nuclei well beyond the iron group (up to $Z \approx 44$) can be produced, and discuss in detail also the nucleosynthesis of the p-nuclei $^{92,94}$Mo and $^{96,98}$Ru. While we observe the production of $^{92}$Mo and $^{94}$Mo in slightly neutron-rich conditions in both simulations, $^{96,98}$Ru can only be produced efficiently via the $\nu$p-process. Furthermore, the production of Ru in the $\nu$p-process heavily depends on the presence of very proton-rich material in the ejecta. This disentanglement of production mechanisms has interesting consequences when comparing to the abundance ratios between these isotopes in the solar system and in presolar grains.

Atmospheric parameters and chemical compositions for ten stars with metallicities in the region of -2.2< [Fe/H] <-0.6 were precisely determined using high resolution, high signal to noise, spectra. For each star the abundances, for 14 to 27 elements, were derived using both LTE and NLTE approaches. In particular, differences by assuming LTE or NLTE are about 0.10 dex; depending on [Fe/H], Teff, gravity and element lines used in the analysis. We find that the O abundance has the largest error, ranging from 0.10 and 0.2 dex. The best measured elements are Cr, Fe, and Mn; with errors etween 0.03 and 0.11 dex. The stars in our sample were included in previous different observational work. We provide a consistent data analysis. The data dispersion introduced in the literature by different techniques and assumptions used by the different authors is within the observational errors, excepting for HD103095. We compare these results with stellar observations from different data sets and a number of theoretical galactic chemical evolution (GCE) simulations. We find a large scatter in the GCE results, used to study the origin of the elements. Within this scatter as found in previous GCE simulations, we cannot reproduce the evolution of the elemental ratios [Sc/Fe], [Ti/Fe], and [V/Fe] at different metallicities. The stellar yields from core collapse supernovae (CCSN) are likely primarily responsible for this discrepancy. Possible solutions and open problems are discussed.

We investigated r-process nucleosynthesis in magneto-rotational supernovae, based on a new explosion mechanism induced by the magneto-rotational instability. A series of axisymmetric magneto-hydrodynamical simulations with detailed microphysics including neutrino heating is performed, numerically resolving the magneto-rotational instability. Neutrino-heating dominated explosions, enhanced by magnetic fields, showed mildly neutron-rich ejecta producing nuclei up to $A \sim 130$ (i.e. the weak r-process), while explosion models with stronger magnetic fields reproduce a solar-like r-process pattern. More commonly seen abundance patterns in our models are in between the weak and regular r-process, producing lighter and intermediate mass nuclei. These {\it intermediate r-processes} exhibit a variety of abundance distributions, compatible with several abundance patterns in r-process-enhanced metal-poor stars. The amount of Eu ejecta $\sim 10^{-5} M_\odot$ in magnetically-driven jets agrees with predicted values in the chemical evolution of early galaxies. In contrast, neutrino-heating dominated explosions have a significant amount of Fe ($^{56}{\rm Ni}$) and Zn, comparable to regular supernovae and hypernovae, respectively. These results indicate magneto-rotational supernovae can produce a wide range of heavy nuclei from iron-group to r-process elements, depending on the explosion dynamics.

We present multi-dimensional core-collapse supernova simulations using the Isotropic Diffusion Source Approximation (IDSA) for the neutrino transport and a modified potential for general relativity in two different supernova codes: FLASH and ELEPHANT. Due to the complexity of the core-collapse supernova explosion mechanism, simulations require not only high-performance computers and the exploitation of GPUs, but also sophisticated approximations to capture the essential microphysics. We demonstrate that the IDSA is an elegant and efficient neutrino radiation transfer scheme, which is portable to multiple hydrodynamics codes and fast enough to investigate long-term evolutions in two and three dimensions. Simulations with a 40 solar mass progenitor are presented in both FLASH (1D and 2D) and ELEPHANT (3D) as an extreme test condition. It is found that the black hole formation time is delayed in multiple dimensions and we argue that the strong standing accretion shock instability before black hole formation will lead to strong gravitational waves.

Core-collapse supernovae (CCSNe) are the extremely energetic deaths of massive stars. They play a vital role in the synthesis and dissemination of many heavy elements in the universe. In the past, CCSN nucleosynthesis calculations have relied on artificial explosion methods that do not adequately capture the physics of the innermost layers of the star. The PUSH method, calibrated against SN1987A, utilizes the energy of heavy-flavor neutrinos emitted by the proto-neutron star (PNS) to trigger parametrized explosions. This makes it possible to follow the consistent evolution of the PNS and to ensure a more accurate treatment of the electron fraction of the ejecta. Here, we present the Iron group nucleosynthesis results for core-collapse supernovae, exploded with PUSH, for two different progenitor series. Comparisons of the calculated yields to observational metal-poor star data are also presented. Nucleosynthesis yields will be calculated for all elements and over a wide range of progenitor masses. These yields can be immensely useful for models of galactic chemical evolution.

The hadron-quark phase transition in core-collapse supernovae (CCSNe) has the potential to trigger explosions in otherwise nonexploding models. However, those hybrid supernova equations of state (EOS) shown to trigger an explosion do not support the observational 2 M$_\odot$ neutron star maximum mass constraint. In this work, we analyze cold hybrid stars by the means of a systematic parameter scan for the phase transition properties, with the aim to develop a new hybrid supernova EOS. The hadronic phase is described with the state-of-the-art supernova EOS HS(DD2), and quark matter by an EOS with a constant speed of sound (CSS) of $c_{QM}^2=1/3$. We find promising cases which meet the 2 M$_\odot$ criterion and are interesting for CCSN explosions. We show that the very simple CSS EOS is transferable into the well-known thermodynamic bag model, important for future application in CCSN simulations. In the second part, the occurrence of reconfinement and multiple phase transitions is discussed. In the last part, the influence of hyperons in our parameter scan is studied. Including hyperons no change in the general behavior is found, except for overall lower maximum masses. In both cases (with and without hyperons) we find that quark matter with $c_{QM}^2=1/3$ can increase the maximum mass only if reconfinement is suppressed or if quark matter is absolutely stable.

A phase transition to quark matter can lead to interesting phenomenological consequences in core-collapse supernovae, e.g., triggering an explosion in spherically symmetric models. However, until now this explosion mechanism was only shown to be working for equations of state that are in contradiction with recent pulsar mass measurements. Here, we identify that this explosion mechanism is related to the existence of a third family of compact stars. For the equations of state investigated, the third family is only pronounced in the hot, early stages of the protocompact star and absent or negligibly small at zero temperature and thus represents a novel kind of third family. This interesting behavior is a result of unusual thermal properties induced by the phase transition, e.g., characterized by a decrease of temperature with increasing density for isentropes, and can be related to a negative slope of the phase transition line in the temperature-pressure phase diagram.

We report on a method, PUSH, for triggering core-collapse supernova (CCSN) explosions of massive stars in spherical symmetry. This method provides a framework to study many important aspects of core collapse supernovae: the effects of the shock passage through the star, explosive supernova nucleosynthesis and the progenitor-remnant connection. Here we give an overview of the method, compare the results to multi-dimensional simulations and investigate the effects of the progenitor and the equation of state on black hole formation.

The $s$-process nucleosynthesis in Asymptotic Giant Branch (AGB) stars depends on the modeling of convective boundaries. We present models and s-process simulations that adopt a treatment of convective boundaries based on the results of hydrodynamic simulations and on the theory of mixing due to gravity waves in the vicinity of convective boundaries. Hydrodynamics simulations suggest the presence of convective boundary mixing (CBM) at the bottom of the thermal pulse-driven convective zone. Similarly, convection-induced mixing processes are proposed for the mixing below the convective envelope during third dredge-up where the 13C pocket for the s process in AGB stars forms. In this work we apply a CBM model motivated by simulations and theory to models with initial mass $M = 2$ and $M = 3M_\odot$, and with initial metal content Z = 0.01 and Z = 0.02. As reported previously, the He-intershell abundance of 12C and 16O are increased by CBM at the bottom of pulse-driven convection zone. This mixing is affecting the $^{22}Ne(\alpha,n)^{25}Mg$ activation and the s-process effciency in the 13C-pocket. In our model CBM at the bottom of the convective envelope during the third dredgeup represents gravity wave mixing. We take further into account that hydrodynamic simulations indicate a declining mixing efficiency already about a pressure scale height from the convective boundaries, compared to mixing-length theory. We obtain the formation of the 13C-pocket with a mass of $\approx 10^{-4}M_\odot$. The final $s$-process abundances are characterized by 0.36 < [s=Fe] < 0.78 and the heavy-to-light s-process ratio is 0.23 < [hs=ls] < 0.45. Finally, we compare our results with stellar observations, pre-solar grain measurements and previous work.

The neutrino mechanism of core-collapse supernova is investigated via non-relativistic, two-dimensional (2D), neutrino radiation-hydrodynamic simulations. For the transport of electron flavor neutrinos, we use the interaction rates defined by Bruenn (1985) and the isotropic diffusion source approximation (IDSA) scheme, which decomposes the transported particles into trapped particle and streaming particle components. Heavy neutrinos are described by a leakage scheme. Unlike the "ray-by-ray" approach in some other multi-dimensional supernova models, we use cylindrical coordinates and solve the trapped particle component in multiple dimensions, improving the proto-neutron star resolution and the neutrino transport in angular and temporal directions. We provide an IDSA verification by performing 1D and 2D simulations with 15 and 20 $M_\odot$ progenitors from Woosley et al.~(2007) and discuss the difference of our IDSA results with those existing in the literature. Additionally, we perform Newtonian 1D and 2D simulations from prebounce core collapse to several hundred milliseconds postbounce with 11, 15, 21, and 27 $M_\odot$ progenitors from Woosley et al.~(2002) with the HS(DD2) equation of state. General relativistic effects are neglected. We obtain robust explosions with diagnostic energies $E_{\rm dig} \gtrsim 0.1- 0.5$~B for all considered 2D models within approximately $100-300$ milliseconds after bounce and find that explosions are mostly dominated by the neutrino-driven convection, although standing accretion shock instabilities are observed as well. We also find that the level of electron deleptonization during collapse dramatically affect the postbounce evolution, e.g.~the ignorance of neutrino-electron scattering during collapse will lead to a stronger explosion.

Rotation was shown to have a strong impact on the structure and light element nucleosynthesis in massive stars. In particular, models including rotation can reproduce the primary nitrogen observed in halo extremely metal-poor (EMP) stars. Additional exploratory models showed that rotation may enhance $s$-process production at low metallicity. Here we present a large grid of massive star models including rotation and a full $s$-process network to study the impact of rotation on the weak $s$-process. We explore the possibility of producing significant amounts of elements beyond the strontium peak, which is where the weak $s$-process usually stops. We used the Geneva stellar evolution code coupled to an enlarged reaction network with 737 nuclear species up to bismuth to calculate $15-40\,\text{M}_\odot$ models at four metallicities ($Z = 0.014,10^{-3}$, $10^{-5}$, and $10^{-7}$) from the main sequence up to the end of oxygen burning. We confirm that rotation-induced mixing between the convective H-shell and He-core enables an important production of primary $^{14}$N and $^{22}$Ne and $s$-process at low metallicity. At low metallicity, even though the production is still limited by the initial number of iron seeds, rotation enhances the $s$-process production, even for isotopes heavier than strontium, by increasing the neutron to seed ratio. The increase in this ratio is a direct consequence of the primary production of $^{22}$Ne. Despite nuclear uncertainties affecting the $s$-process production and stellar uncertainties affecting the rotation-induced mixing, our results show a robust production of $s$ process at low metallicity when rotation is taken into account. Considering models with a distribution of initial rotation rates enables to reproduce the observed large range of the [Sr/Ba] ratios in (carbon-enhanced and normal) EMP stars.

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

In this work we present and discuss the observations of the Mn abundances for 247 FGK dwarfs, located in the Galactic disc with metallicity -1<Fe/H]<+0.3. The observed stars belong to the substructures of the Galaxy thick and thin discs, and to the Hercules stream. The observations were conducted using the 1.93 m telescope at Observatoire de Haute-Provence (OHP, France) equipped with the echelle type spectrographs ELODIE and SOPHIE. The abundances were derived under the LTE approximation, with an average error for the [Mn/Fe] ratio of 0.10 dex. For most of the stars in the sample Mn abundances are not available in the literature. We obtain an evolution of [Mn/Fe] ratio with the metallicity [Fe/H] consistent with previous data compilations. In particular, within the metallicity range covered by our stellar sample the [Mn/Fe] ratio is increasing with the increase of metallicity. This due to the contribution to the Galactic chemical evolution of Mn and Fe from thermonuclear supernovae. We confirm the baseline scenario where most of the Mn in the Galactic disc and in the Sun is made by thermonuclear supernovae. In particular, the effective contribution from core-collapse supernovae to the Mn in the Solar system is about 10-20%. However, present uncertainties affecting the production of Mn and Fe in thermonuclear supernovae are limiting the constraining power of the observed [Mn/Fe] trend in the Galactic discs on, e.g., the frequency of different thermonuclear supernovae populations. The different production of these two elements in different types of thermonuclear supernovae needs to be disentangled by the dependence of their relative production on the metallicity of the supernova progenitor.

For the origin of heavy r-process elements, different sources have been proposed, e.g., core-collapse supernovae or neutron star mergers. Old metal-poor stars carry the signature of the astrophysical source(s). Among the elements dominantly made by the r-process, europium (Eu) is relatively easy to observe. In this work we simulate the evolution of europium in our galaxy with the inhomogeneous chemical evolution model 'ICE', and compare our results with spectroscopic observations. We test the most important parameters affecting the chemical evolution of Eu: (a) for neutron star mergers the coalescence time scale of the merger ($t_{\mathrm{coal}}$) and the probability to experience a neutron star merger event after two supernova explosions occurred and formed a double neutron star system ($P_{\mathrm{NSM}}$) and (b) for the sub-class of magneto-rotationally driven supernovae ("Jet-SNe"), their occurrence rate compared to standard supernovae ($P_{\mathrm{Jet-SN}}$). We find that the observed [Eu/Fe] pattern in the galaxy can be reproduced by a combination of neutron star mergers and magneto-rotationally driven supernovae as r-process sources. While neutron star mergers alone seem to set in at too high metallicities, Jet-SNe provide a cure for this deficiency at low metallicities. Furthermore, we confirm that local inhomogeneities can explain the observed large spread in the europium abundances at low metallicities. We also predict the evolution of [O/Fe] to test whether the spread in $\alpha$-elements for inhomogeneous models agrees with observations and whether this provides constraints on supernova explosion models and their nucleosynthesis.

The r-process nucleosynthesis in core-collapse supernovae (CC-SNe) is studied, with a focus on the explosion scenario induced by rotation and strong magnetic fields. Nucleosynthesis calculations are conducted based on magneto-hydrodynamical explosion models with a wide range of parameters for initial rotation and magnetic fields. The explosion models are classified in two different types: i.e., prompt-magnetic-jet and delayed-magnetic-jet, for which the magnetic fields of proto-neutron stars (PNSs) during collapse and the core-bounce are strong and comparatively moderate, respectively. Following the hydrodynamical trajectories of each explosion model, we confirmed that r-processes successfully occur in the prompt-magnetic-jets, which produce heavy nuclei including actinides. On the other hand, the r-process in the delayed-magnetic-jet is suppressed, which synthesizes only nuclei up to the second peak ($A \sim 130$). Thus, the r-process in the delayed-magnetic-jets could explain only "weak r-process" patterns observed in metal-poor stars rather than the "main r-process", represented by the solar abundances. Our results imply that core-collapse supernovae are possible astronomical sources of heavy r-process elements if their magnetic fields are strong enough, while weaker magnetic explosions may produce "weak r-process" patterns ($A \lesssim 130$). We show the potential importance and necessity of magneto-rotational supernovae for explaining the galactic chemical evolution, as well as abundances of r-process enhanced metal-poor stars. We also examine the effects of the remaining uncertainties in the nature of PNSs due to weak interactions that determine the final neutron-richness of ejecta. Additionally, we briefly discuss radioactive isotope yields in primary jets (e.g., $^{56}$Ni), with relation to several optical observation of SNe and relevant high-energy astronomical phenomena.

Nowadays the number of models aimed at explaining the Type Ia supernova phenomenon is high and discriminating between them is a must-do. In this work we explore the influence of rotation in the evolution of the nuclear flame which drives the explosion in the so called gravitational confined detonation models. Assuming that the flame starts in a point-like region slightly above the center of the white dwarf (WD) and adding a moderate amount of angular velocity to the star we follow the evolution of the deflagration using a smoothed particle hydrodynamics code. We find that the results are very dependent on the angle between the rotational axis and the line connecting the initial bubble of burned material with the center of the white dwarf at the moment of the ignition. The impact of rotation is larger for angles close to 90{\deg} because the Coriolis force on a floating element of fluid is maximum, and its principal effect is to break the symmetry of the deflagration. Such symmetry breaking weakens the convergence of the nuclear flame at the antipodes of the initial ignition volume, changing the environmental conditions around the convergence region with respect to non-rotating models. These changes seem to disfavor the emergence of a detonation in the compressed volume at the antipodes, thus compromising the viability of the so called gravitational confined detonation mechanism.

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.

We report on a method, PUSH, for triggering core-collapse supernova explosions of massive stars in spherical symmetry. We explore basic explosion properties and calibrate PUSH such that the observables of SN1987A are reproduced. Our simulations are based on the general relativistic hydrodynamics code AGILE combined with the detailed neutrino transport scheme IDSA for electron neutrinos and ALS for the muon and tau neutrinos. To trigger explosions in the otherwise non-exploding simulations, we rely on the neutrino-driven mechanism. The PUSH method locally increases the energy deposition in the gain region through energy deposition by the heavy neutrino flavors. Our setup allows us to model the explosion for several seconds after core bounce. We explore the progenitor range 18-21M$_{\odot}$. Our studies reveal a distinction between high compactness (HC) and low compactness (LC) progenitor models, where LC models tend to explore earlier, with a lower explosion energy, and with a lower remnant mass. HC models are needed to obtain explosion energies around 1 Bethe, as observed for SN1987A. However, all the models with sufficiently high explosion energy overproduce $^{56}$Ni. We conclude that fallback is needed to reproduce the observed nucleosynthesis yields. The nucleosynthesis yields of $^{57-58}$Ni depend sensitively on the electron fraction and on the location of the mass cut with respect to the initial shell structure of the progenitor star. We identify a progenitor and a suitable set of PUSH parameters that fit the explosion properties of SN1987A when assuming 0.1M$_{\odot}$ of fallback. We predict a neutron star with a gravitational mass of 1.50M$_{\odot}$. We find correlations between explosion properties and the compactness of the progenitor model in the explored progenitors. However, a more complete analysis will require the exploration of a larger set of progenitors with PUSH.

Baryon density inhomogeneities in the early universe can give rise to a floor of heavy elements (up to $A\approx 270$) produced in a primordial r-process with fission cycling. A parameter study with variation of the global baryon to photon ratio $\eta$, and under inclusion of neutron diffusion effects was performed. New results concerning the dependence of the results on nuclear physics parameters are presented.

A number of possible mechanisms have been suggested to generate density inhomogeneities in the early Universe which could survive until the onset of primordial nucleosynthesis and generate neutron-rich regions. In this work we are not concerned with how the inhomogeneities were generated but we want to focus on the effect of such inhomogeneities on primordial nucleosynthesis. One of the proposed signatures of inhomogeneity, the synthesis of very heavy elements by neutron capture in a primordial r-process, was analyzed for varying baryon to photon ratios $\eta$ and fluctuation length scales $L$.

This is a White Paper in support of the mission concept of the Large Observatory for X-ray Timing (LOFT), proposed as a medium-sized ESA mission. We discuss the potential of LOFT for the study of thermonuclear X-ray bursts on accreting neutron stars. For a summary, we refer to the paper.

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.

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.

Carbon-rich grains are observed to condense in the ejecta of recent core-collapse supernovae, within a year after the explosion. Silicon carbide grains of type X are C-rich grains with isotpic signatures of explosive supernova nucleosynthesis have been found in primitive meteorites. Much rarer silicon carbide grains of type C are a special sub-group of SiC grains from supernovae. They show peculiar abundance signatures for Si and S, isotopically heavy Si and isotopically light S, which appear to to be in disagreement with model predictions. We propose that C grains are formed mostly from C-rich stellar material exposed to lower SN shock temperatures than the more common type X grains. In this scenario, extreme $^{32}$S enrichments observed in C grains may be explained by the presence of short-lived $^{32}$Si ($\tau$$_{1/2}$ = 153 years) in the ejecta, produced by neutron capture processes starting from the stable Si isotopes. No mixing from deeper Si-rich material and/or fractionation of Si from S due to molecular chemistry is needed to explain the $^{32}$S enrichments. The abundance of $^{32}$Si in the grains can provide constraints on the neutron density reached during the supernova explosion in the C-rich He shell material. The impact of the large uncertainty of the neutron capture cross sections in the $^{32}$Si region is discussed.

We compare three different statistical models for the equation of state (EOS) of stellar matter at subnuclear densities and temperatures (0.5-10 MeV) expected to occur during the collapse of massive stars and supernova explosions. The models introduce the distributions of various nuclear species in nuclear statistical equilibrium, but use somewhat different nuclear physics inputs. It is demonstrated that the basic thermodynamical quantities of stellar matter under these conditions are similar, except in the region of high densities and low temperatures. We demonstrate that mass and isotopic distributions have considerable differences related to the different assumptions of the models on properties of nuclei at these stellar conditions. Overall, the three models give similar trends, but the details reflect the uncertainties related to the modeling of medium effects, such as the temperature and density dependence of surface and bulk energies of heavy nuclei, and the nuclear shell structure effects. We discuss importance of new physics inputs for astrophysical calculations from experimental data obtained in intermediate energy heavy-ion collisions, in particular, the similarities of the conditions reached during supernova explosions and multifragmentation reactions.

About a year after core collapse supernova, dust starts to condense in the ejecta. In meteorites, a fraction of C-rich presolar grains (e.g., silicon carbide (SiC) grains of Type-X and low density graphites) are identified as relics of these events, according to the anomalous isotopic abundances. Several features of these abundances remain unexplained and challenge the understanding of core-collapse supernovae explosions and nucleosynthesis. We show, for the first time, that most of the measured C-rich grain abundances can be accounted for in the C-rich material from explosive He burning in core-collapse supernovae with high shock velocities and consequent high temperatures. The inefficiency of the $^{12}$C($\alpha$,$\gamma$)$^{16}$O reaction relative to the rest of the $\alpha$-capture chain at $T > 3.5\times10^8 \mathrm{K}$ causes the deepest He-shell material to be carbon rich and silicon rich, and depleted in oxygen. The isotopic ratio predictions in part of this material, defined here as the C/Si zone, are in agreement with the grain data. The high-temperature explosive conditions that our models reach at the bottom of the He shell, can also be representative of the nucleosynthesis in hypernovae or in the high-temperature tail of a distribution of conditions in asymmetric supernovae. Finally, our predictions are consistent with the observation of large $^{44}$Ca/$^{40}$Ca observed in the grains. This is due to the production of $^{44}$Ti together with $^{40}$Ca in the C/Si zone, and/or to the strong depletion of $^{40}$Ca by neutron captures.

The aim of this work is to present and discuss the observations of the iron peak (Fe, Ni) and neutron-capture element (Y, Zr, Ba, La, Ce, Nd, Sm, and Eu) abundances for 276 FGK dwarfs, located in the galactic disk with metallicity -1 < [Fe/H] < +0.3. Atmospheric parameters and chemical composition of the studied stars were determined from an high resolution, high signal-to-noise echelle spectra obtained with the echelle spectrograph ELODIE at the Observatoire de Haute-Provence (France). Effective temperatures were estimated by the line depth ratio method and from the H_{\alpha} line-wing fitting. Surface gravities (log g) were determined by parallaxes and the ionization balance of iron. Abundance determinations were carried out using the LTE approach, taking the hyperfine structure for Eu into account, and the abundance of Ba was computed under the NLTE approximation. We are able to assign most of the stars in our sample to the substructures of the Galaxy thick disk, thin disk, or Hercules stream according to their kinematics. The classification of 27 stars is uncertain. For most of the stars in the sample, the abundances of neutron-capture elements have not been measured earlier. For all of them, we provide the chemical composition and discuss the contribution from different nucleosynthesis processes.

We explore heavy element nucleosynthesis in the explosion of massive stars which are triggered by a quark-hadron phase transition during the early post bounce phase of core-collapse supernovae. The present study is based on general relativistic radiation hydrodynamics simulations with three-flavor Boltzmann neutrino transport in spherical symmetry, which utilize a quark-hadron hybrid equation of state based on the MIT bag model for strange quark matter. The quark-hadron phase transition inside the stellar core forms a shock wave propagating towards the surface of the proto-neutron star. The shock wave results in an explosion and ejects neutron-rich matter which is piled up or accreting on the proto-neutron star. Later, during the cooling phase, the proto-neutron star develops a proton-rich neutrino-driven wind. We present a detailed analysis of the nucleosynthesis outcome in both neutron-rich and proton-rich ejecta and compare our integrated nucleosynthesis with observations of metal poor stars.

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 examine simulations of core-collapse supernovae in spherical symmetry. Our model is based on general relativistic radiation hydrodynamics with three-flavor Boltzmann neutrino transport. We discuss the different supernova phases, including the long-term evolution up to 20 seconds after the onset of explosion during which the neutrino fluxes and mean energies decrease continuously. In addition, the spectra of all flavors become increasingly similar, indicating the change from charged- to neutral-current dominance. Furthermore, it has been shown recently by several groups independently, based on sophisticated supernova models, that collective neutrino flavor oscillations are suppressed during the early mass-accretion dominated post-bounce evolution. Here we focus on the possibility of collective flavor flips between electron and non-electron flavors during the later, on the order of seconds, evolution after the onset of an explosion with possible application for the nucleosynthesis of heavy elements.

We discuss the possible impact of strange quark matter on the evolution of core-collapse supernovae with emphasis on low critical densities for the quark-hadron phase transition. For such cases the hot proto-neutron star can collapse to a more compact hybrid star configuration hundreds of milliseconds after core-bounce. The collapse triggers the formation of a second shock wave. The latter leads to a successful supernova explosion and leaves an imprint on the neutrino signal. These dynamical features are discussed with respect to their compatibility with recent neutron star mass measurements which indicate a stiff high density nuclear matter equation of state.

Context. Rotation is known to affect the nucleosynthesis of light elements in massive stars, mainly by rotation-induced mixing. In particular, rotation boosts the primary nitrogen production. Models of rotating stars are able to reproduce the nitrogen observed in low-Z halo stars. Aims. Here we present the first grid of stellar models for rotating massive stars at low Z, where a full s-process network is used to study the impact of rotation-induced mixing on the nucleosynthesis of heavy elements. Methods. We used the Geneva stellar evolution code that includes an enlarged reaction network with nuclear species up to bismuth to calculate 25 M$_\odot$ models at three different Z and with different initial rotation rates. Results. First, we confirm that rotation-induced mixing leads to a production of primary $^{22}$Ne, which is the main neutron source for the s process in massive stars. Therefore rotation boosts the s process in massive stars at all Z. Second, the neutron-to-seed ratio increases with decreasing Z in models including rotation, which leads to the complete consumption of all iron seeds at Z < 1e-3 by the end of core He-burning. Thus at low Z, the iron seeds are the main limitation for this boosted s process. Third, as Z decreases, the production of elements up to the Ba peak increases at the expense of the elements of the Sr peak. We studied the impact of the initial rotation rate and of the uncertain $^{17}$O$(\alpha,\gamma)$ rate (which strongly affects the neutron poison strength of $^{16}$O) on our results. This study shows that rotating models can produce significant amounts of elements up to Ba over a wide range of Z. Fourth, compared to the He-core, the primary $^{22}$Ne production in the He-shell is even higher (> 1% in mass fraction at all Z), which could open the door for an explosive neutron capture nucleosynthesis in the He-shell, with a primary neutron source.

Proton-activation reactions on natural and enriched palladium samples were investigated via the activation technique in the energy range of E_p=2.75 MeV to 9 MeV, close to the upper end of the respective Gamow window of the gamma process. We have determined cross sections for 102Pd(p,gamma)103Ag, 104Pd(p,gamma)105Ag, and 105Pd(p,n)105Ag, as well as partial cross sections of 104Pd(p,n)104Ag^g, 105Pd(p,gamma)106Ag^m, 106Pd(p,n)106Ag^m, and 110Pd(p,n)110Ag^m with uncertainties between 3% and 15% for constraining theoretical Hauser-Feshbach rates and for direct use in gamma-process calculations.

We propose to perform a continuously scanning all-sky survey from 200 keV to 80 MeV achieving a sensitivity which is better by a factor of 40 or more compared to the previous missions in this energy range. The Gamma-Ray Imaging, Polarimetry and Spectroscopy (GRIPS) mission addresses fundamental questions in ESA's Cosmic Vision plan. Among the major themes of the strategic plan, GRIPS has its focus on the evolving, violent Universe, exploring a unique energy window. We propose to investigate $\gamma$-ray bursts and blazars, the mechanisms behind supernova explosions, nucleosynthesis and spallation, the enigmatic origin of positrons in our Galaxy, and the nature of radiation processes and particle acceleration in extreme cosmic sources including pulsars and magnetars. The natural energy scale for these non-thermal processes is of the order of MeV. Although they can be partially and indirectly studied using other methods, only the proposed GRIPS measurements will provide direct access to their primary photons. GRIPS will be a driver for the study of transient sources in the era of neutrino and gravitational wave observatories such as IceCUBE and LISA, establishing a new type of diagnostics in relativistic and nuclear astrophysics. This will support extrapolations to investigate star formation, galaxy evolution, and black hole formation at high redshifts.

We explore explosions of massive stars, which are triggered via the quark-hadron phase transition during the early post bounce phase of core-collapse supernovae. We construct a quark equation of state, based on the bag model for strange quark matter. The transition between the hadronic and the quark phases is constructed applying Gibbs conditions. The resulting quark-hadron hybrid equations of state are used in core-collapse supernova simulations, based on general relativistic radiation hydrodynamics and three flavor Boltzmann neutrino transport in spherical symmetry. The formation of a mixed phase reduces the adiabatic index, which induces the gravitational collapse of the central protoneutron star. The collapse halts in the pure quark phase, where the adiabatic index increases. A strong accretion shock forms, which propagates towards the protoneutron star surface. Due to the density decrease of several orders of magnitude, the accretion shock turns into a dynamic shock with matter outflow. This moment defines the onset of the explosion in supernova models that allow for a quark-hadron phase transition, where otherwise no explosions could be obtained. The shock propagation across the neutrinospheres releases a burst of neutrinos. This serves as a strong observable identification for the structural reconfiguration of the stellar core. The ejected matter expands on a short timescale and remains neutron-rich. These conditions might be suitable for the production of heavy elements via the r-process. The neutron-rich material is followed by proton-rich neutrino-driven ejecta in the later cooling phase of the protoneutron star where the vp-process might occur.

Context. Light elements and nitrogen surface abundances together can constrain the mixing efficiencies in massive stars on the main sequence, because moderate mixing at the surface leads to a depletion of light elements but only later to an enrichment in nitrogen. Aims. We want to test the rotational mixing prescriptions included in the Geneva stellar evolution code (GENEC) by following the evolution of surface abundances of light isotopes in massive stars. Methods. The GENEC is a 1D code containing sophisticated prescriptions for rotational mixing. We implemented an extended reaction network into this code including the light elements Li, Be and B, which allowed us to perform calculations testing the rotation induced mixing. Results. We followed 9, 12 and 15 solar mass models with rotation from the zero age main sequence up to the end of He burning. The calculations show the expected behaviour with faster depletion of light isotopes for faster rotating stars and more massive stars. Conclusions. We find that the mixing prescriptions used in the present rotating models for massive single stars can account for most of the observations. However the uncertainties are quite large making it hard to draw a firm conclusion on the mixing scenario.

Massive stars end their life in an explosion event with kinetic energies of the order 1 Bethe. Immediately after the explosion has been launched, a region of low density and high entropy forms behind the ejecta which is continuously subject to neutrino heating. The neutrinos emitted from the remnant at the center, the protoneutron star (PNS), heat the material above the PNS surface. This heat is partly converted into kinetic energy and the material accelerates to an outflow that is known as the neutrino driven wind. For the first time, we simulate the collapse, bounce, explosion and the neutrino driven wind phases consistently over more than 20 seconds. Our numerical model is based on spherically symmetric general relativistic radiation hydrodynamics using spectral three flavor Boltzmann neutrino transport. In simulations where no explosions are obtained naturally, we model neutrino driven explosions for low and intermediate mass Fe-core progenitor stars by enhancing the charged current reaction rates. In the case of a special progenitor star, the O-Ne-Mg-core, the explosion in spherical symmetry was obtained without enhanced opacities. The post explosion evolution is in qualitative agreement with static steady-state and parametrized dynamic models of the neutrino driven wind. On the other hand, we find generally smaller neutrino luminosities and mean neutrino energies as well as a different evolutionary behavior of the neutrino luminosities and mean neutrino energies. The neutrino driven wind is proton-rich for more than 10 seconds and the contraction of the PNS differs from the assumptions made for the conditions at the inner boundary in previous neutrino driven wind studies. Despite the moderately large entropies per baryon of about 100 and the fast expansion timescale, the conditions found in our model are unlikely to favor...

Massive stars and their supernovae are prominent sources of radioactive isotopes, the observations of which thus can help to improve our astrophysical models of those. Our understanding of stellar evolution and the final explosive endpoints such as supernovae or hypernovae or gamma-ray bursts relies on the combination of magneto-hydrodynamics, energy generation due to nuclear reactions accompanying composition changes, radiation transport, and thermodynamic properties (such as the equation of state of stellar matter). Nuclear energy production includes all nuclear reactions triggered during stellar evolution and explosive end stages, also among unstable isotopes produced on the way. Radiation transport covers atomic physics (e.g. opacities) for photon transport, but also nuclear physics and neutrino nucleon/nucleus interactions in late phases and core collapse. Here we want to focus on the astrophysical aspects, i.e. a description of the evolution of massive stars and their endpoints, with a special emphasis on the composition of their ejecta (in form of stellar winds during the evolution or of explosive ejecta). Low and intermediate mass stars end their evolution as a white dwarf with an unburned C and O composition. Massive stars evolve beyond this point and experience all stellar burning stages from H over He, C, Ne, O and Si-burning up to core collapse and explosive endstages. In this chapter we discuss the nucleosynthesis processes involved and the production of radioactive nuclei in more detail.

For finite chemical potential effective models of QCD predict a first order phase transition. In favour for the search of such a phase transition in nature, we construct an equation of state for strange quark matter based on the MIT bag model. We apply this equation of state to highly asymmetric core collapse supernova matter with finite temperatures and large baryon densities. The phase transition is constructed using the general Gibbs conditions, which results in an extended coexistence region between the pure hadronic and pure quark phases in the phase diagram, i.e. the mixed phase. The supernovae are simulated via general relativistic radiation hydrodynamics based on three flavor Boltzmann neutrino transport in spherical symmetry. During the dynamical evolution temperatures above 10 MeV, baryon densities above nuclear saturation density and a proton-to-baryon ratio below 0.2 are obtained. At these conditions the phase transition is triggered which leads to a significant softening of the EoS for matter in the mixed phase. As a direct consequence of the stiffening of the EoS again for matter in the pure quark phase, a shock wave forms at the boundary between the mixed and the pure hadronic phases. This shock is accelerated and propagates outward which releases a burst of neutrinos dominated by electron anti-neutrinos due to the lifted degeneracy of the shock-heated hadronic material. We discuss the radiation-hydrodynamic evolution of the phase transition at the example of several low and intermediate mass Fe-core progenitor stars and illustrate the expected neutrino signal from the phase transition.

Neutron-induced reaction rates, including fission, are calculated in the temperature range 1.d8 <T (K) < 1.d10 within the framework of the statistical model for targets with atomic number 83 < Z < 119 (from Po to Uuo) from the neutron to the proton drip-line. Four sets of rates have been calculated, utilizing - where possible - consistent nuclear data for neutron separation energies and fission barriers from Thomas-Fermi (TF), Extended Thomas-Fermi plus Strutinsky Integral (ETFSI), Finite-Range Droplet Model (FRDM) and Hartree-Fock-Bogolyubov (HFB) predictions. Tables of calculated values as well as analytic seven parameter fits in the standard REACLIB format are supplied. We also discuss the sensitivity of the rates to the input, aiming at a better understanding of the uncertainties introduced by the nuclear input.

Explosive astrophysical systems, such as supernovae or compact star binary mergers, provide conditions where strange quark matter can appear. The high degree of isospin asymmetry and temperatures of several MeV in such systems may cause a transition to the quark phase already around saturation density. Observable signals from the appearance of quark matter can be predicted and studied in astrophysical simulations. As input in such simulations, an equation of state with an integrated quark matter phase transition for a large temperature, density and proton fraction range is required. Additionally, restrictions from heavy ion data and pulsar observation must be considered. In this work we present such an approach. We implement a quark matter phase transition in a hadronic equation of state widely used for astrophysical simulations and discuss its compatibility with heavy ion collisions and pulsar data. Furthermore, we review the recently studied implications of the QCD phase transition during the early post-bounce evolution of core-collapse supernovae and introduce the effects from strong interactions to increase the maximum mass of hybrid stars. In the MIT bag model, together with the strange quark mass and the bag constant, the strong coupling constant $\alpha_s$ provides a parameter to set the beginning and extension of the quark phase and with this the mass and radius of hybrid stars.

The nucleosynthetic characteristics of various explosion mechanisms of Type Ia supernovae (SNe Ia) is explored based on three two-dimensional explosion simulations representing extreme cases: a pure turbulent deflagration, a delayed detonation following an approximately spherical ignition of the initial deflagration, and a delayed detonation arising from a highly asymmetric deflagration ignition. Apart from this initial condition, the deflagration stage is treated in a parameter-free approach. The detonation is initiated when the turbulent burning enters the distributed burning regime. This occurs at densities around $10^{7}$ g cm$^{-3}$ -- relatively low as compared to existing nucleosynthesis studies for one-dimensional spherically symmetric models. The burning in these multidimensional models is different from that in one-dimensional simulations as the detonation wave propagates both into unburned material in the high density region near the center of a white dwarf and into the low density region near the surface. Thus, the resulting yield is a mixture of different explosive burning products, from carbon-burning products at low densities to complete silicon-burning products at the highest densities, as well as electron-capture products synthesized at the deflagration stage. In contrast to the deflagration model, the delayed detonations produce a characteristic layered structure and the yields largely satisfy constraints from Galactic chemical evolution. In the asymmetric delayed detonation model, the region filled with electron capture species (e.g., $^{58}$Ni, $^{54}$Fe) is within a shell, showing a large off-set, above the bulk of $^{56}$Ni distribution, while species produced by the detonation are distributed more spherically (abridged).

We have investigated the (n,gamma) cross sections of p-process isotopes with the activation technique. The measurements were carried out at the Karlsruhe Van de Graaff accelerator using the 7Li(p,n)7Be source for simulating a Maxwellian neutron distribution of kT = 25 keV. Stellar cross section measurements are reported for the light p-process isotopes 102Pd, 120Te, 130,132Ba, and 156Dy. In a following paper the cross sections of 168Yb, 180W, 184Os, 190Pt, and 196Hg will be discussed. The data are extrapolated to p-process energies by including information from evaluated nuclear data libraries. The results are compared to standard Hauser-Feshbach models frequently used in astrophysics.

The implications of the formation of strange quark matter in neutron stars and in core-collapse supernovae is discussed with special emphasis on the possibility of having a strong first order QCD phase transition at high baryon densities. If strange quark matter is formed in core-collapse supernovae shortly after the bounce, it causes the launch of a second outgoing shock which is energetic enough to lead to a explosion. A signal for the formation of strange quark matter can be read off from the neutrino spectrum, as a second peak in antineutrinos is released when the second shock runs over the neutrinosphere.

We discuss the formation of stellar mass black holes via protoneutron star (PNS) collapse. In the absence of an earlier explosion, the PNS collapses to a black hole due to the continued mass accretion onto the PNS. We present an analysis of the emitted neutrino spectra of all three flavors during the PNS contraction. Special attention is given to the physical conditions which depend on the input physics, e.g. the equation of state (EoS) and the progenitor model. The PNSs are modeled as the central object in core collapse simulations using general relativistic three-flavor Boltzmann neutrino transport in spherical symmetry. The simulations are launched from several massive progenitors of 40 and 50 solar mass. We analyze the electron-neutrino luminosity dependencies and construct a simple approximation for the electron-neutrino luminosity, which depends only on the physical conditions at the electron-neutrinosphere. In addition, we analyze different mu/tau-neutrino pair-reactions separately and compare the differences during the post-bounce phase of failed core collapse supernova explosions of massive progenitors. We also investigate the connection between the increasing mu/tau-neutrino luminosity and the PNS contraction during the accretion phase before black hole formation. Comparing the different post bounce phase of the progenitor models under investigation, we find large differences in the emitted neutrino spectra. These differences and the analysis of the electron-neutrino luminosity indicate a strong progenitor model dependency of the emitted neutrino signal.

The possible role of a first order QCD phase transition at nonvanishing quark chemical potential and temperature for cold neutron stars and for supernovae is delineated. For cold neutron stars, we use the NJL model with nonvanishing color superconducting pairing gaps, which describes the phase transition to the 2SC and the CFL quark matter phases at high baryon densities. We demonstrate that these two phase transitions can both be present in the core of neutron stars and that they lead to the appearance of a third family of solution for compact stars. In particular, a core of CFL quark matter can be present in stable compact star configurations when slightly adjusting the vacuum pressure to the onset of the chiral phase transition from the hadronic model to the NJL model. We show that a strong first order phase transition can have strong impact on the dynamics of core collapse supernovae. If the QCD phase transition sets in shortly after the first bounce, a second outgoing shock wave can be generated which leads to an explosion. The presence of the QCD phase transition can be read off from the neutrino and antineutrino signal of the supernova.

We propose to advance investigations of electromagnetic radiation originating in atomic nuclei beyond its current infancy to a true astronomy. This nuclear emission is independent from conditions of gas, thus complements more traditional stronomical methods used to probe the nearby universe. Radioactive gamma-rays arise from isotopes which are made in specific locations inside massive stars, their decay in interstellar space traces an otherwise not directly observable hot and tenuous phase of the ISM, which is crucial for feedback from massive stars. Its intrinsic clocks can measure characteristic times of processes within the ISM. Frontier questions that can be addressed with studies in this field are the complex interiors of massive stars and supernovae which are key agents in galactic dynamics and chemical evolution, the history of star-forming and supernova activity affecting our solar-system environment, and explorations of occulted and inaccessible regions of young stellar nurseries in our Galaxy.

While the high-entropy wind (HEW) of Type II supernovae remains one of the more promising sites for the rapid neutron-capture (r-) process, hydrodynamic simulations have yet to reproduce the astrophysical conditions under which the latter occurs. We have performed large-scale network calculations within an extended parameter range of the HEW, seeking to identify or to constrain the necessary conditions for a full reproduction of all r-process residuals N_{r,\odot}=N_{\odot}-N_{s,\odot} by comparing the results with recent astronomical observations. A superposition of weighted entropy trajectories results in an excellent reproduction of the overall N_{r,\odot}-pattern beyond Sn. For the lighter elements, from the Fe-group via Sr-Y-Zr to Ag, our HEW calculations indicate a transition from the need for clearly different sources (conditions/sites) to a possible co-production with r-process elements, provided that a range of entropies are contributing. This explains recent halo-star observations of a clear non-correlation of Zn and Ge and a weak correlation of Sr - Zr with heavier r-process elements. Moreover, new observational data on Ru and Pd seem to confirm also a partial correlation with Sr as well as the main r-process elements (e.g. Eu).

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.

The s-process production in massive stars at very low metallicities is expected to be negligible due to the low abundance of the neutron source 22Ne, to primary neutron poisons and decreasing iron seed abundances. However, recent models of massive stars including the effects of rotation show that a strong production of 22Ne is possible in the helium core, as a consequence of the primary nitrogen production (observed in halo metal poor stars). Using the PPN post-processing code, we studied the impact of this primary 22Ne on the s process. We find a large production of s elements between strontium and barium, starting with the amount of primary 22Ne predicted by stellar models. There are several key reaction rate uncertainties influencing the s-process efficiency. Among them, 17O(alpha,gamma) may play a crucial role strongly influencing the s process efficiency, or it may play a negligible role, according to the rate used in the calculations. We also report on the development of a new parallel (MPI) post-processing code (MPPNP) designed to follow the complete nucleosynthesis in stars on highly resolved grids. We present here the first post-processing run from the ZAMS up to the end of helium burning for a 15 solar mass model.

We explore the implications of the QCD phase transition during the postbounce evolution of core-collapse supernovae. Using the MIT bag model for the description of quark matter and assuming small bag constants, we find that the phase transition occurs during the early postbounce accretion phase. This stage of the evolution can be simulated with general relativistic three-flavor Boltzmann neutrino transport. The phase transition produces a second shock wave that triggers a delayed supernova explosion. If such a phase transition happens in a future galactic supernova, its existence and properties should become observable as a second peak in the neutrino signal that is accompanied by significant changes in the energy of the emitted neutrinos. In contrast to the first neutronization burst, this second neutrino burst is dominated by the emission of anti-neutrinos because the electron-degeneracy is lifted when the second shock passes through the previously neutronized matter.

The distribution of angular momentum of the matter during X-ray bursts on neutron stars is studied by means of 3D axi-symmetric general relativistic hydrodynamics. The set of fully general relativistic Navier-Stokes equations is solved implicitly using the implicit solver GR-I-RMHD in combination with a third order spatial and second order temporal advection scheme. The viscous operators are formulated using a Kerr-like metric in the fixed background of a slowly rotating neutron star whose radius coincides with the corresponding last stable orbit. The importance of these operators and their possible simplifications are discussed as well. In the rotating case and depending on the viscosity parameter, $\alpha_\mm{tur}$, it is found that the viscously-initiated fronts at the center of bursts propagate at much faster speed than the fluid motion. These fast fronts act to decouple angular momentum from matter: angular momentum is transported outwards while matter sinks inwards into the deep gravitational well of the neutron star, thereby enhancing the compression of matter necessary for initiating ignition, that subsequently spreads over the whole surface of the neutron star on the viscous time scale. Based on the numerical simulations, we find that a viscosity parameter $\alpha_\mm{tur} = \mathcal{O}(0.1)$ is most suitable for fitting observations of neutron stars during X-ray bursts. It is argued that the spin up observed in the cooling tails of X-ray bursts is a transient phase, which eventually should be followed by a spin down phase. This delay can be attributed to a significant lengthening of the viscous time scale due to rapid cooling of matter in the outer layers.

Massive stars played a key role in the early evolution of the Universe. They formed with the first halos and started the re-ionisation. It is therefore very important to understand their evolution. In this paper, we describe the strong impact of rotation induced mixing and mass loss at very low $Z$. The strong mixing leads to a significant production of primary nitrogen 14, carbon 13 and neon 22. Mass loss during the red supergiant stage allows the production of Wolf-Rayet stars, type Ib,c supernovae and possibly gamma-ray bursts (GRBs) down to almost Z=0 for stars more massive than 60 solar masses. Galactic chemical evolution models calculated with models of rotating stars better reproduce the early evolution of N/O, C/O and C12/C13. We calculated the weak s-process production induced by the primary neon 22 and obtain overproduction factors (relative to the initial composition, Z=1.e-6) between 100-1000 in the mass range 60-90.

Since April 2005 a regularly updated stellar neutron cross section compilation is available online at http://nuclear-astrophysics.fzk.de/kadonis. This online-database is called the "Karlsruhe Astrophysical Database of Nucleosynthesis in Stars" project and is based on the previous Bao et al. compilation from the year 2000. The present version \textsc{KADoNiS} v0.2 (January 2007) includes recommended cross sections for 280 isotopes between $^{1}$H and $^{210}$Po and 75 semi-empirical estimates for isotopes without experimental information. Concerning stellar $(n,\gamma)$ cross sections of the 32 stable, proton-rich isotopes produced by the $p$ process experimental information is only available for 20 isotopes, but 9 of them have rather large uncertainties of $\geq$9%. The first part of a systematic study of stellar $(n,\gamma)$ cross sections of the $p$-process isotopes $^{74}$Se, $^{84}$Sr, $^{102}$Pd, $^{120}$Te, $^{130}$Ba, $^{132}$Ba, $^{156}$Dy, and $^{174}$Hf is presented. In another application \textsc{KADoNiS} v0.2 was used for an modification of a reaction library of Basel university. With this modified library $p$-process network calculations were carried out and compared to previous results.

The nucleosynthesis of elements beyond iron is dominated by the s and r processes. However, a small amount of stable isotopes on the proton-rich side cannot be made by neutron capture and are thought to be produced by photodisintegration reactions on existing seed nuclei in the so-called "p process". So far most of the p-process reactions are not yet accessible by experimental techniques and have to be inferred from statistical Hauser-Feshbach model calculations. The parametrization of these models has to be constrained by measurements on stable proton-rich nuclei. A series of (n,$\gamma$) activation measurements on p nuclei, related by detailed balance to the respective photodisintegrations, were carried out at the Karlsruhe Van de Graaff accelerator using the $^7$Li(p,n)$^7$Be source for simulating a Maxwellian neutron distribution of kT= 25 keV. We present here preliminary results of our extended measuring program in the mass range between A=74 and A=132, including first experimental (n,$\gamma$) cross sections of $^{74}$Se, $^{84}$Sr, $^{120}$Te and $^{132}$Ba, and an improved value for $^{130}$Ba. In all cases we find perfect agreement with the recommended MACS predictions from the Bao et al. compilation.

The "Karlsruhe Astrophysical Database of Nucleosynthesis in Stars" (KADoNiS) project is an online database for experimental cross sections relevant to the s process and p process. It is available under http://nuclear-astrophysics.fzk.de/kadonis and consists of two parts. Part 1 is an updated sequel to the previous Bao et al. compilations from 1987 and 2000 for (n,$\gamma$) cross sections relevant to the big bang and s-process nucleosynthesis. The second part will be an experimental p-process database, which is expected to be launched in winter 2005/06. The KADoNiS project started in April 2005, and a first partial update is online since August 2005. In this paper we present a short overview of the first update of the s-process database, as well as an overview of the status of stellar (n,$\gamma$) cross sections of all 32 p isotopes.

We have performed $p$-process simulations with the most recent stellar $(n,\gamma)$ cross sections from the "Karlsruhe Astrophysical Database of Nucleosynthesis in Stars" project (version v0.2, http://nuclear-astrophysics.fzk.de/kadonis). The simulations were carried out with a parametrized supernova type II shock front model (``$\gamma$ process'') of a 25 solar mass star and compared to recently published results. A decrease in the normalized overproduction factor could be attributed to lower cross sections of a significant fraction of seed nuclei located in the Bi and Pb region around the $N$=126 shell closure.

The nucleosynthesis of elements beyond iron is dominated by the s and r processes. However, a small amount of stable isotopes on the proton-rich side cannot be made by neutron capture and are thought to be produced by photodisintegration reactions on existing seed nuclei in the so-called "p process". So far most of the p-process reactions are not yet accessible by experimental techniques and have to be inferred from statistical Hauser-Feshbach model calculations. The parametrization of these models has to be constrained by measurements on stable proton-rich nuclei. A series of (n,$\gamma$) activation measurements, related by detailed balance to the respective photodisintegrations, were carried out at the Karlsruhe Van de Graaff accelerator using the $^7$Li(p,n)$^7$Be source for simulating a Maxwellian neutron distribution of kT= 25 keV. First results for the experimental (n,$\gamma$) cross sections of the light p nuclei $^{74}$Se and $^{84}$Sr are reported. These experimental values were used for an extrapolation to the Maxwellian averaged cross section at 30 keV, $<\sigma>_{30}$, yielding 271$\pm$15 mb for $^{74}$Se, and 300$\pm$17 mb for the total capture cross section of $^{84}$Sr. The partial cross section to the isomer in $^{85}$Sr was found to be 190$\pm$10 mb.

Microscopic calculations show a strong parity dependence of the nuclear level density at low excitation energy of a nucleus. Previously, this dependence has either been neglected or only implemented in the initial and final channels of Hauser-Feshbach calculations. We present an indirect way to account for a full parity dependence in all steps of a reaction, including the one of the compound nucleus formed in a reaction. To illustrate the impact on astrophysical reaction rates, we present rates for neutron captures in isotopic chains of Ni and Sn. Comparing with the standard assumption of equipartition of both parities, we find noticeable differences in the energy regime of astrophysical interest caused by the parity dependence of the nuclear level density found in the compound nucleus even at sizeable excitation energies.

Internal flows inside gravitationally stable astrophysical objects, such as the Sun, stars and compact stars are compressed and extremely subsonic. Such low Mach number flows are usually encountered when studying for example dynamo action in stars, planets, the hydro-thermodynamics of X-ray bursts on neutron stars and dwarf novae. Treating such flows is numerically complicated and challenging task. We aim to present a robust numerical tool that enables modeling the time-evolution or quasi-stationary of stratified low Mach number flows under astrophysical conditions. It is argued that astrophysical low Mach number flows cannot be considered as an asymptotic limit of incompressible flows, but rather as highly compressed flows with extremely stiff pressure terms. Unlike the pseudo-pressure in incompressible fluids, a Possion-like treatment for the pressure would smooth unnecessarily the physically induced acoustic perturbations, thereby violating the conservation character of the compressible equations. Moreover, classical dimensional splitting techniques, such as ADI or Line-Gauss-Seidel methods are found to be unsuited for modeling compressible flows with low Mach numbers. In this paper we present a nonlinear Newton-type solver that is based on the defect-correction iteration procedure and in which the Approximate Factorization Method (AFM) is used as a preconditioner. This solver is found to be sufficiently robust and is capable of capturing stationary solutions for viscous rotating flows with Mach number as small as $\mcal{M} \approx 10^{-3},$ i.e., near the incompressibility limit.

Few years ago, Boltzmann neutrino transport led to a new and reliable generation of spherically symmetric models of stellar core collapse and postbounce evolution. After the failure to prove the principles of the supernova explosion mechanism, these sophisticated models continue to illuminate the close interaction between high-density matter under extreme conditions and the transport of leptons and energy in general relativistically curved space-time. We emphasize that very different input physics is likely to be relevant for the different evolutionary phases, e.g. nuclear structure for weak rates in collapse, the equation of state of bulk nuclear matter during bounce, multidimensional plasma dynamics in the postbounce evolution, and neutrino cross sections in the explosive nucleosynthesis. We illustrate the complexity of the dynamics using preliminary 3D MHD high-resolution simulations based on parameterized deleptonization. With established spherically symmetric models we show that typical features of the different phases are reflected in the predicted neutrino signal and that a consistent neutrino flux leads to electron fractions larger than 0.5 in neutrino-driven supernova ejecta.

A generalized method to calculate the excitation-energy dependent parity ratio in the nuclear level density is presented, using the assumption of Poisson distributed independent quasi particles combined with BCS occupation numbers. It is found that it is crucial to employ a sufficiently large model space to allow excitations both from low-lying shells and to higher shells beyond a single major shell. Parity ratios are only found to equilibrate above at least 5-10 MeV of excitation energy. Furthermore, an overshooting effect close to major shells is found where the parity opposite to the ground state parity may dominate across a range of several MeV before the parity ratio finally equilibrates. The method is suited for large-scale calculations as needed, for example, in astrophysical applications. Parity distributions were computed for all nuclei from the proton dripline to the neutron dripline and from Ne up to Bi. These results were then used to recalculate astrophysical reaction rates in a Hauser-Feshbach statistical model. Although certain transitions can be considerably enhanced or suppressed, the impact on astrophysically relevant reactions remains limited, mainly due to the thermal population of target states in stellar reaction rates.

Explosive hydrogen burning in type I X-ray bursts (XRBs) comprise charged particle reactions creating isotopes with masses up to A~100. Since charged particle reactions in a stellar environment are very temperature sensitive, we use a realistic time-dependent general relativistic and self-consistent model of type I x-ray bursts to provide accurate values of the burst temperatures and densities. This allows a detailed and accurate time-dependent identification of the reaction flow from the surface layers through the convective region and the ignition region to the neutron star ocean. Using this, we determine the relative importance of specific nuclear reactions in the X-ray burst.

This manuscript reviews recent progress in our understanding of the nucleosynthesis of medium and heavy elements in supernovae. Recent hydrodynamical models of core-collapse supernovae show that a large amount of proton rich matter is ejected under strong neutrino fluxes. This matter constitutes the site of the vp-process where antineutrino absorption reactions catalyze the nucleosynthesis of nuclei with A > 64. Supernovae are also associated with the r-process responsible for the synthesis of the heaviest elements in nature. Fission during the r-process can play a major role in determining the final abundance patter and in explaining the almost universal features seen in metal-poor r-process-rich stars.

Recently, measurements of abundances in extremely metal poor (EMP) stars have brought new constraints on stellar evolution models. In an attempt to explain the origin of the abundances observed, we computed pre--supernova evolution models, explosion models and the related nucleosynthesis. In this paper, we start by presenting the pre-SN models of rotating single stars with metallicities ranging from solar metallicity down to almost metal free. We then review key processes in core-collapse and bounce, before we integrate them in a simplistic parameterization for 3D MHD models, which are well underway and allow one to follow the evolution of the magnetic fields during collapse and bounce. Finally, we present explosive nucleosynthesis results including neutrino interactions with matter, which are calculated using the outputs of the explosion models. The main results of the pre-SN models are the following. First, primary nitrogen is produced in large amount in models with an initial metallicity $Z=10^{-8}$. Second, at the same metallicity of $Z=10^{-8}$ and for models with an initial mass larger than about 60 Mo, rotating models may experience heavy mass loss (up to more than half of the initial mass of the star). The chemical composition of these winds can qualitatively reproduce the abundance patterns observed at the surface of carbon-rich EMP stars. Explosive nucleosynthesis including neutrino-matter interactions produce improved abundances for iron group elements, in particular for scandium and zinc. It also opens the way to a new neutrino and proton rich process ($\nu$p-process) able to contribute to the nucleosynthesis of elements with A > 64. (Abridged)

Core collapse supernovae are the leading actor in the story of the cosmic origin of the chemical elements. Existing models, which generally assume spherical symmetry and parameterize the explosion, have been able to broadly replicate the observed elemental pattern. However, inclusion of neutrino interactions produces noticeable improvement in the composition of the ejecta when compared to observations. Neutrino interactions may also provide a supernova source for light p-process nuclei.

With presently known input physics and computer simulations in 1D, a self-consistent treatment of core collapse supernovae does not yet lead to successful explosions, while 2D models show some promise. Thus, there are strong indications that the delayed neutrino mechanism works combined with a multi-D convection treatment for unstable layers. On the other hand there is a need to provide correct nucleosynthesis abundances for the progressing field of galactic evolution and observations of low metallicity stars. The innermost ejecta is directly affected by the explosion mechanism, i.e. most strongly the yields of Fe-group nuclei for which an induced piston or thermal bomb treatment will not provide the correct yields because the effect of neutrino interactions is not included. We apply parameterized variations to the neutrino scattering cross sections and alternatively, parameterized variations are applied to the neutrino absorption cross sections on nucleons in the ``gain region''. We find that both measures lead to similar results, causing explosions and a Ye>0.5 in the innermost ejected layers, due to the combined effect of a short weak interaction time scale and a negligible electron degeneracy, unveiling the proton-neutron mass difference. We include all weak interactions (electron and positron capture, beta-decay, neutrino and antineutrino capture on nuclei, and neutrino and antineutrino capture on nucleons) and present first nucleosynthesis results for these innermost ejected layers to discuss how they improve predictions for Fe-group nuclei. The proton-rich environment results in enhanced abundances of 45Sc, 49Ti, and 64Zn as requested by chemical evolution studies and observations of low metallicity stars as well as appreciable production of nuclei in the mass range up to A=80.

We present a new nucleosynthesis process, that we denote nu p-process, which occurs in supernovae (and possibly gamma-ray bursts) when strong neutrino fluxes create proton-rich ejecta. In this process, antineutrino absorptions in the proton-rich environment produce neutrons that are immediately captured by neutron-deficient nuclei. This allows for the nucleosynthesis of nuclei with mass numbers A >64. Making this process a possible candidate to explain the origin of the solar abundances of 92,94Mo and 96,98Ru. This process also offers a natural explanation for the large abundance of Sr seen in an hyper-metal-poor star.

Fission plays an important role in the r-process which is responsible not only for the yields of transuranium isotopes, but may have a strong influence on the formation of the majority of heavy nuclei due to fission recycling. We present calculations of beta-delayed and neutron-induced fission rates, taking into account different fission barriers predictions and mass formulae. It is shown that an increase of fission barriers results naturally in a reduction of fission rates, but that nevertheless fission leads to the termination of the r-process. Furthermore, it is discussed that the probability of triple fission could be high for $A>260$ and have an effect on the formation of the abundances of heavy nuclei. Fission after beta-delayed neutron emission is discussed as well as different aspects of the influence of fission upon r-process calculations.

This paper reports on the detailed rp-process reaction flow on an accreting neutron star and the resulting ashes of a type I X-ray burst. It is obtained by coupling a 298 isotope reaction network to a self-consistent one-dimensional model calculation with a constant accretion rate of dM/dt=1.0e17g/s (0.09 Eddington).

We ensure successful explosions (of otherwise non-explosive models) by enhancing the neutrino luminosity via reducing the neutrino scattering cross sections or by increasing the heating efficiency via enhancing the neutrino absorption cross sections in the heating region. Our investigations show that the resulting electron fraction Ye in the innermost ejecta is close to 0.5, in some areas even exceeding 0.5. We present the effects of the resulting values for Ye on the nucleosynthesis yields of the innermost zones of core collapse supernovae.

The light curves of type I X-ray bursts (XRBs) result from energy released from the atmosphere of a neutron star when accreted hydrogen and helium ignite and burn explosively via the rp-process. Since charged particle reaction rates are both density and very temperature dependent, a simulation model must provide accurate values of these variables to predict the reaction flow. This paper uses a self-consistent one-dimensional model calculation with a constant accretion rate of dM/dt=5e16g/s (0.045 Eddington) and reports on the detailed rp-process reaction flow of a given burst.

We present the results of nucleosynthesis calculations based on multidimensional (2D and 3D) hydrodynamical simulations of the thermonuclear burning phase in SNIa. The detailed nucleosynthetic yields of our explosion models are calculated by post-processing the ejecta, using passively advected tracer particles. The nuclear reaction network employed in computing the explosive nucleosynthesis contains 383 nuclear species. We analyzed two different choices of ignition conditions (centrally ignited, in which the spherical initial flame geometry is perturbated with toroidal rings, and bubbles, in which multi-point ignition conditions are simulated). We show that unburned C and O varies typically from ~40% to ~50% of the total ejected material.The main differences between all our models and standard 1D computations are, besides the higher mass fraction of unburned C and O, the C/O ratio (in our case is typically a factor of 2.5 higher than in 1D computations), and somewhat lower abundances of certain intermediate mass nuclei such as S, Cl, Ar, K, and Ca, and of 56Ni. Because explosive C and O burning may produce the iron-group elements and their isotopes in rather different proportions one can get different 56Ni-fractions (and thus supernova luminosities) without changing the kinetic energy of the explosion. Finally, we show that we need the high resolution multi-point ignition (bubbles) model to burn most of the material in the center (demonstrating that high resolution coupled with a large number of ignition spots is crucial to get rid of unburned material in a pure deflagration SNIa model).

Cross sections for the photon-induced particle-emission reactions (gamma,n), (gamma,p), and (gamma,alpha) are given for all natural isotopes from Ti to Bi. The target nuclei are assumed to be in their ground states, except for 180Ta which is naturally occurring as the isomer 180mTa. The cross sections are calculated in a statistical model (Hauser-Feshbach) approach and covering an energy range from threshold up to 7.35 MeV above the threshold (14.7 MeV above threshold for (gamma,alpha) reactions). The results are intended to aid conception and analysis of experiments which can also be used to test the methods involved in predicting astrophysical reaction rates for nucleosynthesis.

We consider the influence of potential quark matter existing at high densities in neutron star interiors on gravitational waves (GW) emitted in a binary neutron star merger event. Two types of equations of state (EoS) at zero temperatures are used, one describing pure nuclear matter, the other nuclear matter with a phase transition to quark matter at very high densities. Binary equilibrium sequences close to the innermost stable circular orbit (ISCO) are calculated to determine the GW frequencies just before merger. It is found that EoS effects begin to play a role for gravitational masses larger than $M_\infty\simeq1.5M_\odot$. The difference in the gravitational wave frequency at the ISCO grows to up to $\simeq 10%$ for the maximal allowed mass given by the EoSs used. Then, we perform 3D hydrodynamic simulations for each EoS varying the initial mass and determine the characteristic GW frequencies of the merger remnants. The implications of quark matter show up mainly in a different collapse behaviour of the merger remnant. If the collapse does not take place immediately after merger, we find a phase difference between two EoS's in the post-merger GW signal. We also compare the GW frequencies emitted by the merger remnant to values from simulations using a polytropic EoS and find an imprint of the non-constant adiabatic index of our EoSs. All calculations are based on the conformally flat (CF) approximation to general relativity and the GW signal from the merger simulation is extracted up to quadrupole order.

Type I X-ray bursts with a double peak in the bolometric luminosity have been observed from several sources. The separation between the two peaks are on the order of a few seconds. We propose a nuclear waiting point impedance in the thermonuclear reaction flow to explain these observations. Nuclear structure information suggests the potential waiting points: Mg22, Si26, S30 and Ar34, which arise in conditions, where a further reaction flow has to await a beta-decay, because the (alpha,p)-reaction is too weak to overcome the target Coulomb-barrier and the (p,gamma)-reaction is quenched by photo-disintegration at the burst temperature. The conclusion is that the effects of the experimentally unknown S30(alpha,p)Cl33 and Ar34(alpha,p)K37 might be directly visible in the observation of X-ray burst light curves.

We present an implicit finite difference representation for general relativistic radiation hydrodynamics in spherical symmetry. Our code, Agile-Boltztran, solves the Boltzmann transport equation for the angular and spectral neutrino distribution functions in self-consistent simulations of stellar core collapse and postbounce evolution. It implements a dynamically adaptive grid in comoving coordinates. Most macroscopically interesting physical quantities are defined by expectation values of the distribution function. We optimize the finite differencing of the microscopic transport equation for a consistent evolution of important expectation values. We test our code in simulations launched from progenitor stars with 13 solar masses and 40 solar masses. ~0.5 s after core collapse and bounce, the protoneutron star in the latter case reaches its maximum mass and collapses further to form a black hole. When the hydrostatic gravitational contraction sets in, we find a transient increase in electron flavor neutrino luminosities due to a change in the accretion rate. The muon- and tauon-neutrino luminosities and rms energies, however, continue to rise because previously shock-heated material with a non-degenerate electron gas starts to replace the cool degenerate material at their production site. We demonstrate this by supplementing the concept of neutrinospheres with a more detailed statistical description of the origin of escaping neutrinos. We compare the evolution of the 13 solar mass progenitor star to simulations with the MGFLD approximation, based on a recently developed flux limiter. We find similar results in the postbounce phase and validate this MGFLD approach for the spherically symmetric case with standard input physics.

The astrophysical nature of r-process sites is a long standing mystery and many probable sources have been suggested in the past, among them lower-mass core-collapse supernovae (in the range 8 - 10 Msol), higher-mass core-collapse supernovae (with masses > 20 Msol) and neutron star mergers. In this work, we present a detailed inhomogeneous chemical evolution study that considers for the first time neutron star mergers as major r-process sources, and compare this scenario to the ones in which core-collapse supernovae act as dominant r-process sites. We conclude that, due to the lack of reliable iron and r-process yields as function of progenitor mass, it is not possible at present to distinguish between the lower-mass and higher-mass supernovae scenarios within the framework of inhomogeneous chemical evolution. However, neutron-star mergers seem to be ruled out as the dominant r-process source, since their low rates of occurrence would lead to r-process enrichment that is not consistent with observations at very low metallicities. Additionally, the considerable injection of r-process material by a single neutron-star merger leads to a scatter in [r-process/Fe] ratios at later times which is much too large compared to observations.

General relativistic multi-group and multi-flavor Boltzmann neutrino transport in spherical symmetry adds a new level of detail to the numerical bridge between microscopic nuclear and weak interaction physics and the macroscopic evolution of the astrophysical object. Although no supernova explosions are obtained, we investigate the neutrino luminosities in various phases of the postbounce evolution for a wide range of progenitor stars between 13 and 40 solar masses. The signal probes the dynamics of material layered in and around the protoneutron star and is, within narrow limits, sensitive to improvements in the weak interaction physics. Only changes that dramatically exceed physical limitations allow experiments with exploding models. We discuss the differences in the neutrino signal and find the electron fraction in the innermost ejecta to exceed 0.5 as a consequence of thermal balance and weak equilibrium at the masscut.

We describe an implicit general relativistic hydrodynamics code. The evolution equations are formulated in comoving coordinates. A conservative finite differencing of the Einstein equations is outlined, and artificial viscosity and numerical diffusion are discussed. The time integration is performed with AGILE, an implicit solver for stiff algebrodifferential equations on a dynamical adaptive grid. We extend the adaptive grid technique, known from nonrelativistic hydrodynamics, to the general relativistic application and identify it with the concept of shift vectors in a 3+1 decomposition. The adaptive grid minimizes the number of required computational zones without compromising the resolution in physically important regions. Thus, the computational effort is greatly reduced when the zones are subject to computationally expensive additional processes, such as Boltzmann radiation transport or a nuclear reaction network. We present accurate results in the standard tests for supernova simulations: Sedov's point-blast explosion, the nonrelativistic and relativistic shock tube, the Oppenheimer-Snyder dust collapse, and homologous collapse.

Inhomogeneous chemical evolution models of galaxies which try to reproduce the scatter seen in element-to-iron ratios of metal-poor halo stars are heavily dependent on theoretical nucleosynthesis yields of core-collapse supernovae. Hence inhomogeneous models present themselves as a test for stellar nucleosyn- thesis calculations. Applying an inhomogeneous chemical evolution model to our Galaxy reveals a number of shortcomings of existing theoretical nucleosynthesis yields. One problem is the predicted scatter in [O/Fe] and [Mg/Fe] which is too large compared to the one observed in metal-poor halo stars. This can be either due to the O or Mg yields or due to the Fe yields (or both). However, O and Mg are alpha-elements that are produced mainly during hydrostatic burning and thus are not affected by the theoretical uncertainties afflicting the collapse and explosion of a massive star. Stellar iron yields, on the other hand, depend heavily on the choice of the mass-cut between ejecta and proto neutron star and are therefore very uncertain. We present Fe yield distributions as function of progenitor mass that are consistent with the abundance distribution of metal- poor halo stars and are in agreement with observed Ni yields of SNe II with known progenitor masses. The iron yields of lower-mass SNe II (in the range 10-20 Msol) are well constrained by those observations. Present observations, however, do not allow to determine a unique solution for higher-mass SNe. Nevertheless, the main dependence of the stellar iron yield as function of progenitor mass may be derived and can be used as constraint for future supernova/hypernova models. The results are of importance for the earliest stages of galaxy formation when the ISM is dominated by chemical inhomogeneities and the instantaneous mixing approximation is not valid.

Based on general relativistic three-flavor Boltzmann neutrino transport, we compare simulations of the postbounce evolution for five different stellar progenitors with 13 and 20 solar masses (Nomoto & Hashimoto '88), and 15, 25, and 40 solar masses (Woosley & Weaver '95).

With the exception of the Big Bang, responsible for 1,2H, 3,4He, and 7Li, stars act as sources for the composition of the interstellar medium. Cosmic rays are related to the latter and very probably due to acceleration of the mixed interstellar medium by shock waves from supernova remnants. Thus, the understanding of the abundance evolution in the interstellar medium and especially the enrichment of heavy elements, as a function of space and time, is essential. It reflects the history of star formation and the lifetimes of the diverse contributing stellar objects. Therefore, the understanding of the endpoints of stellar evolution is essential as well. These are mainly planetary nebulae and type II/Ib/Ic supernovae as evolutionary endpoints of single stars, but also events in binary systems can contribute, like e.g. supernovae of type Ia, novae and possibly X-ray bursts and neutron star or neutron star - black hole mergers. Despite many efforts, a full and self-consistent understanding of supernovae (the main contributors to nucleosynthesis in galaxies) is not existing, yet. Their fingerprints, however, seen either in spectra, lightcurves, radioactivities/decay gamma-rays or in galactic evolution, can help to constrain the composition of their ejecta and related model uncertainties.

We present a new 3D SPH code which solves the general relativistic field + hydrodynamics equations in the conformally flat approximation. Several test cases are considered to test different aspects of the code. We finally apply then the code to the coalescence of a neutron star binary system. The neutron stars are modeled by a polytropic equation of state (EoS) with adiabatic indices $\Gamma=2.0$, $\Gamma=2.6$ and $\Gamma=3.0$. We calculate the gravitational wave signals, luminosities and frequency spectra by employing the quadrupole approximation for emission and back reaction in the slow motion limit. In addition, we consider the amount of ejected mass.

We report on the stellar core collapse, bounce, and postbounce evolution of a 13 solar mass star in a self-consistent general relativistic spherically symmetric simulation based on Boltzmann neutrino transport. We conclude that approximations to exact neutrino transport and omission of general relativistic effects were not alone responsible for the failure of numerous preceding attempts to model supernova explosions in spherical symmetry. Compared to simulations in Newtonian gravity, the general relativistic simulation results in a smaller shock radius. We however argue that the higher neutrino luminosities and rms energies in the general relativistic case could lead to a larger supernova explosion energy.

The description of general relativistic radiation hydrodynamics in spherical symmetry is presented in natural coordinate choices. For hydrodynamics, comoving coordinates are chosen, and the momentum phase space for the radiation particles is described in comoving frame four momenta. We also investigate a description of the momentum phase space in terms of particle impact parameter and energy at infinity and derive a simple approximation to the general relativistic Boltzmann equation. Further developed are, however, the exact equations in comoving coordinates, because the description of the interaction between matter and radiation particles is best described in the closely related orthonormal basis comoving with the fluid elements. We achieve a conservative and concise formulation of radiation hydrodynamics that is well suited for numerical implementation by a variety of methods. The contribution of radiation to the general relativistic jump conditions at shock fronts is discussed, and artificial viscosity is consistently included in the derivations in order to support approaches relying on this option.

In a previous publication [Atomic Data and Nuclear Data Tables 75, 1 (2000)], we had given seven parameter analytical fits to theoretical reaction rates derived from nuclear cross sections calculated in the statistical model (Hauser-Feshbach formalism) for targets with 9<Z<84 (Ne to Bi) and for a mass range reaching the neutron and proton driplines. Reactions considered were (n,gamma), (n,p), (n,alpha), (p,gamma), (p,alpha), (alpha,gamma), and their inverse reactions. On public demand, we present here the theoretical nuclear cross sections and astrophysical reaction rates from which those rate fits were derived, and we provide these data as on-line electronic files. Corresponding to the fitted rates, two complete data sets are provided, one of which includes a phenomenological treatment of shell quenching for neutron-rich nuclei.

We examine the characteristics of nucleosynthesis in 'hypernovae', i.e., supernovae with very large explosion energies ($ \gsim 10^{52} $ ergs). We carry out detailed nucleosynthesis calculations for these energetic explosions and compare the yields with those of ordinary core-collapse supernovae. We find that both complete and incomplete Si-burning takes place over more extended, lower density regions, so that the alpha-rich freezeout is enhanced and produces more Ti in comparison with ordinary supernova nucleosynthesis. In addition, oxygen and carbon burning takes place in more extended, lower density regions than in ordinary supernovae. Therefore, the fuel elements O, C, Al are less abundant while a larger amount of Si, S, Ar, and Ca ("Si") are synthesized by oxygen burning; this leads to larger ratios of "Si"/O in the ejecta. Enhancement of the mass ratio between complete and incomplete Si-burning regions in the ejecta may explain the abundance ratios among iron-peak elements in metal-poor stars. Also the enhanced "Si"/O ratio may explain the abundance ratios observed in star burst galaxies. We also discuss other implications of enhanced [Ti/Fe] and [Fe/O] for Galactic chemical evolution and the abundances of low mass black hole binaries.

We calculate the rapid proton (rp) capture process of hydrogen burning on the surface of an accreting neutron star with an updated reaction network that extends up to Xe, far beyond previous work. In both steady-state nuclear burning appropriate for rapidly accreting neutron stars (such as the magnetic polar caps of accreting X-ray pulsars) and unstable burning of Type I X-ray bursts, we find that the rp process ends in a closed SnSbTe cycle. This prevents the synthesis of elements heavier than Te and has important consequences for X-ray burst profiles, the composition of accreting neutron stars, and potentially galactic nucleosynthesis of light p nuclei.

Except for 1H, 2H, 3He, 4He, and 7Li, originating from the Big Bang, all heavier elements are made in stellar evolution and stellar explosions. Nuclear physics, and in many cases nuclear structure far from stability, enters in a crucial way. Therefore, we examine in this review the role of nuclear physics in astrophysics in general and in particular how it affects stellar events and the resulting nucleosynthesis. Stellar modeling addresses four major aspects: 1. energy generation and nucleosynthesis, 2. energy transport via conduction, radiation or possibly convection, 3. hydrodynamics/hydrostatics, and finally 4. thermodynamic properties of the matter involved. Nuclear Physics enters via nuclear reaction cross sections and nuclear structure (affecting the composition changes and nuclear energy generation), neutrino-nucleon and neutrino-nucleus cross sections (affecting neutrino opacities and transport), and e.g. the equation of state at and beyond nuclear densities which creates a relation between the nuclear many-body problem and the hydrodynamic response like pressure and entropy. In the following we review these four topics by highlighting the role and impact of nuclear physics in each of these aspects of stellar modeling. The main emphasis is put on the connection to element synthesis.

We examine the characteristics of nucleosynthesis in 'hypernovae', i.e., supernovae with very large explosion energies ($\gsim 10^{52}$ ergs). Implications for the cosmic chemical evolution and the abundances in M82 are discussed.

With exact three-flavor Boltzmann neutrino transport, we simulate the stellar core collapse, bounce, and postbounce evolution of a 13 solar mass star in spherical symmetry, the Newtonian limit, without invoking convection. In the absence of convection, prior spherically symmetric models, which implemented approximations to Boltzmann transport, failed to produce explosions. We are motivated to consider exact transport to determine if these failures were due to the transport approximations made and to answer remaining fundamental questions in supernova theory. The model presented here is the first in a sequence of models beginning with different progenitors. In this model, a supernova explosion is not obtained. We discuss the ramifications of our results for the supernova mechanism.

In the framework of our investigation to explain the nucleosynthesis origin of the correlated Ca-Ti-Cr isotopic anomalies in the Ca-Al-rich ''FUN'' inclusion EK-1-4-1 of the Allende meteorite, the nuclear-physics basis in the neutron-rich N=28 region has been updated by including recent experimental data on beta-decay properties and microscopic predictions of neutron-capture cross sections. Charged-particle and subsequent r-process calculations within an entropy-based approach were performed using a complete reaction network. It is shown that there exist two astrophysical scenarios within which the observed isotopic anomalies can be reproduced simultaneously; one at low entropies (about 10) which confirms the earlier suggestrd Sn Ia mechanism, and another at high entropies (about 150) which could be compatible with the neutrino-wind scenario of a SN II.

We present the results of fully dynamical r-process network calculations for the ejecta of neutron star mergers (NSMs). The late stages of the inspiral and the final violent coalescence of a neutron star binary have been calculated in detail using a 3D hydrodynamics code (Newtonian gravity plus backreaction forces emerging from the emission of gravitational waves) and a realistic nuclear equation of state. The found trajectories for the ejecta serve as input for dynamical r-process calculations where all relevant nuclear reactions (including beta-decays depositing nuclear energy in the expanding material) are followed. We find that all the ejected material undergoes r-process. For an initial Ye close to 0.1 the abundance distributions reproduce very accurately the solar r-process pattern for nuclei with A above 130. For lighter nuclei strongly underabundant (as compared to solar) distributions are encountered. We show that this behaviour is consistent with the latest observations of very old, metal-poor stars, despite simplistic arguments that have recently been raised against the possibility of NSM as possible sources of Galactic r-process material.

The majority of nuclear reactions in astrophysics involve unstable nuclei which are not fully accessible by experiments yet. Therefore, there is high demand for reliable predictions of cross sections and reaction rates by theoretical means. The majority of reactions can be treated in the framework of the statistical model (Hauser-Feshbach). The global parametrizations of the nuclear properties needed for predictions far off stability probe our understanding of the strong force and take it to its limit. The sensitivity of astrophysical scenarios to nuclear inputs is illustrated in the framework of a detailed nucleosynthesis study in type II supernovae. Abundances resulting from calculations in the same explosion model with two different sets of reaction rates are compared. Key reactions and required nuclear information are identified.

The hydrogen-rich accreted envelopes of neutron stars in binary systems are the site for the rapid proton capture nucleosynthethic process (rp-process), which involves nuclei close to the proton dripline. An overview of the relevant reactions and nuclear properties for rp-process studies is given, along with motivations for further experimental nuclear physics studies.

Results of 3D, Newtonian hydrodynamic calculations of the last stages of the inspiral and the final coalescence of neutron star binary systems are described. The focus is on sligthly asymmetric systems with the asymmetry stemming either from different spins or different masses of the binary components. Implications for nucleosynthesis and gamma ray bursts are discussed. Test calculations for the used viscosity scheme are provided in the appendix.

The reaction rates for alpha capture processes on self-conjugate nuclei in the mass range A=20-40 have been investigated. The rates were calculated using the statistical model code NON-SMOKER taking into account isospin suppression rules. These theoretical predictions are compared with rates derived from the available experimental data about the alpha capture reactions taking also into account additional experimental information from different reaction channels populating the alpha unbound states of the self-conjugate compound nuclei.

Theoretical reaction rates in the temperature range 0.01*10^9<=T[K]<=10.*10^9 are calculated in the statistical model (Hauser-Feshbach formalism) for targets with 10<=Z<=83 (Ne to Bi) and for a mass range reaching the neutron and proton driplines. Reactions considered are (n,gamma), (n,p), (n,alpha), (p,gamma), (p,alpha),(alpha,gamma), and their inverse reactions. Reaction rates as a function of temperature for thermally populated targets are given by analytic seven parameter fits. To facilitate comparison with experimental rates, the stellar enhancement factors are also tabulated. Two complete sets of rates have been calculated, one of which includes a phenomenological treatment of shell quenching for neutron-rich nuclei. These extensive datasets are provided on-line as electronic files, while a selected subset from one calculation is given as printed tables. A summary of the theoretical inputs and advice on the use of the provided tabulations is included.

We calculate explosive nucleosynthesis in Chandrasekhar mass models for Type Ia Supernovae(SNe Ia) to obtain new constraints on the rate of matter accretion onto the progenitor white dwarf and on the ignition density of central carbon deflagration. The calculated abundance of the Fe-group neutron-rich nuclei is highly sensitive to the electron captures taking place in the central layers. The yields obtained from a slow central deflagration, and from a fast deflagration or delayed detonation in the outer layers, are combined and put to comparison with solar isotopic abundances. We found that (1) to avoid too large ratios of $^{54}$Cr/$^{56}$Fe and $^{50}$Ti/$^{56}$Fe, the ignition density should be as low as \ltsim 2 \e9 \gmc, and that (2) to avoid the overproduction of $^{58}$Ni and $^{54}$Fe, either the flame speed should not exceed a few % of the sound speed in the central low $Y_e$ layers, or the progenitor star has to be metal-poor compared with solar. Such low central densities can be realized by a rapid accretion as fast as $\dot M$ \gtsim 1 $\times$ 10$^{-7}$M$_\odot$ yr$^{-1}$. In order to reproduce the solar abundance of $^{48}$Ca, one also needs progenitor systems that undergo ignition at higher densities. We also found that the total amount of $^{56}$Ni, the Si-Ca/Fe ratio, and the abundance of elements like Mn and Cr (incomplete Si-burning ashes), depend on the density of the deflagration-detonation transition in delayed detonations. Our nucleosynthesis results favor transition densities slightly below 2.2$\times 10^7$ g cm$^{-3}$.

We introduce a stochastic halo formation model to compute the early chemical enrichment of the interstellar medium (ISM) of the halo. Contrary to 1-zone chemical evolution models, we are able to resolve local inhomogeneities in the ISM caused by single core-collapse supernovae. These inhomogeneities lead to different element abundance patterns in very metal-poor stars, which can be seen as scatter in the abundances of halo stars with metallicities [Fe/H] < 2.0. The early chemical evolution of the halo proceeds in different enrichment phases: At [Fe/H] < -3.0, the halo ISM is unmixed and dominated by local inhomogeneities caused by single core-collapse supernova (SN) events. For metallicities [Fe/H] >-2.0 the halo ISM is well mixed, showing an element abundance pattern integrated over the initial mass function. In the range -3.0 < [Fe/H] < -2.0 a continuous transition from the unmixed to the well mixed ISM occurs. For some elements (Si, Ca, Eu), the scatter in the element-to-iron ratio [El/Fe] seen in the observations of metal-poor halo stars can be reproduced. Stellar yields of other elements predict a scatter which, compared to the observations, is too large (O, Mg) or too small (Ni). Cr and Mn show a decreasing trend for lower metallicities, which can not be explained by metallicity independent yields, provided that the mixing of the ejecta with the interstellar medium does not depend on progenitor mass. This demonstrates the need for revised, self-consistent core-collapse SN yields.

The Chandrasekhar mass model for Type Ia Supernovae (SNe Ia) has received increasing support from recent comparisons of observations with light curve predictions and modeling of synthetic spectra. It explains SN Ia events via thermonuclear explosions of accreting white dwarfs in binary stellar systems, being caused by central carbon ignition when the white dwarf approaches the Chandrasekhar mass. As the electron gas in white dwarfs is degenerate, characterized by high Fermi energies for the high density regions in the center, electron capture on intermediate mass and Fe-group nuclei plays an important role in explosive burning. Electron capture affects the central electron fraction Y_e, which determines the composition of the ejecta from such explosions. Up to the present, astrophysical tabulations based on shell model matrix elements were only available for light nuclei in the sd-shell. Recently new Shell Model Monte Carlo (SMMC) and large-scale shell model diagonalization calculations have also been performed for pf-shell nuclei. These lead in general to a reduction of electron capture rates in comparison with previous, more phenomenological, approaches. Making use of these new shell model based rates, we present the first results for the composition of Fe-group nuclei produced in the central regions of SNe Ia and possible changes in the constraints on model parameters like ignition densities and burning front speeds.

The reaction rates of neutron capture reactions on light nuclei are important for reliably simulating nucleosynthesis in a variety of stellar scenarios. Neutron capture reaction rates on neutron-rich C-, N-, and O-isotopes are calculated in the framework of a hybrid compound and direct capture model. The results are tabulated and compared with the results of previous calculations as well as with experimental results.

The focus of the present review is the production of the heaviest elements in nature via the r-process. A correct understanding and modeling requires the knowledge of nuclear properties far from stability and a detailed prescription of the astrophysical environment. Experiments at CERN/ISOLDE have played a pioneering role in exploring the characteristics of nuclear structure in terms of masses and beta-decay properties. Initial examinations paid attention to far unstable nuclei with magic neutron numbers related to r-process peaks, while present activities are centered on the evolution of shell effects with the distance from the valley of stability. We first show in site-independent applications the effect of both types of nuclear properties on r-process abundances. Then, we explore the results of calculations related to two different `realistic' astrophysical sites, (i) the supernova neutrino wind and (ii) neutron star mergers. We close with a list of remaining theoretical and experimental challenges needed to overcome for a full understanding of the nature of the r-process, and the role CERN/ISOLDE can play in this process.

This review concentrates on the two principle methods used to evolve nuclear abundances within astrophysical simulations, evolution via rate equations and via equilibria. Because in general the rate equations in nucleosynthetic applications form an extraordinarily stiff system, implicit methods have proven mandatory, leading to the need to solve moderately sized matrix equations. Efforts to improve the performance of such rate equation methods are focused on efficient solution of these matrix equations, by making best use of the sparseness of these matrices. Recent work to produce hybrid schemes which use local equilibria to reduce the computational cost of the rate equations is also discussed. Such schemes offer significant improvements in the speed of reaction networks and are accurate under circumstances where calculations with complete equilibrium fail.

Rapid neutron-capture (i.e., r-process) nucleosynthesis calculations, employing internally consistent and physically realistic nuclear physics input (QRPA beta-decay rates and the ETFSI-Q nuclear mass model), have been made. These calculations are compared with ground-based and HST observations of neutron-capture elements in the metal poor halo stars CS 22892--052, HD 115444, HD 122563 and HD 126238. The elemental abundances in all four metal-poor stars are consistent with the solar r-process elemental distribution for the elements Z >/= 56. These results strongly suggest, at least for those elements, that the relative elemental r-process abundances have not changed over the history of the Galaxy. This further suggests that there is one r-process site in the Galaxy, at least for elements Z >/= 56. Employing the observed stellar abundances of stable elements, in conjunction with the solar r-process abundances to constrain the calculations, predictions for the zero decay-age abundances of the radioactive elements Th and U are made. We compare these (least-squares-fit ETFSI-Q) predictions with newly derived observational values in three very metal-poor halo stars: HD 115444, CS 22892--052 and HD 122563. Within the observational errors the observed ratio of [Th/Eu] is the same in both CS 22892--052 and HD 115444. Comparing with the theoretical ratio suggests an average age of these two very metal poor stars to be approximately 15.6 +/- 4.6 Gyr, consistent with earlier radioactive age estimates and recent globular and cosmological age estimates.

We explore the sensitivity of the nucleosynthesis of intermediate mass elements (28 < A < 80) in supernovae derived from massive stars to the nuclear reaction rates employed in the model. Two standard sources of reaction rate data (Woosley et al. 1978; and Thielemann et al. 1987) are employed in pairs of calculations that are otherwise identical. Both include as a common backbone the experimental reactions rates of Caughlan & Fowler (1988). Two stellar models are calculated for each of two main sequence masses: 15 and 25 solar masses. Each star is evolved from core hydrogen burning to a presupernova state carrying an appropriately large reaction network, then exploded using a piston near the edge of the iron core as described by Woosley & Weaver (1995). The final stellar yields from the models calculated with the two rate sets are compared and found to differ in most cases by less than a factor of two over the entire range of nuclei studied. Reasons for the major discrepancies are discussed in detail along with the physics underlying the two reaction rate sets employed. The nucleosynthesis results are relatively robust and less sensitive than might be expected to uncertainties in nuclear reaction rates, though they are sensitive to the stellar model employed.

Proc. IX. Workshop on Nuclear Astrophysics, Ringberg Castle, Tegernsee, Bavaria; March 1998

We explore the effects on nucleosynthesis in Type II supernovae of various parameters (mass cut, neutron excess, explosion energy, progenitor mass) in order to explain the observed trends of the iron-peak element abundance ratios ([Cr/Fe], [Mn/Fe], [Co/Fe] and [Ni/Fe]) in halo stars as a function of metallicity for the range $ -4 \le$ [Fe/H] $\le -2.5$. [Cr/Fe] and [Mn/Fe] decrease with decreasing [Fe/H], while [Co/Fe] behaves the opposite way and increases. We show that such a behavior can be explained by a variation of mass cuts in Type II supernovae as a function of progenitor mass, which provides a changing mix of nucleosynthesis from an alpha-rich freeze-out of Si-burning and incomplete Si-burning. This explanation is consistent with the amount of ejected $^{56}$Ni determined from modeling the early light curves of individual supernovae. We also suggest that the ratio [H/Fe] of halo stars is mainly determined by the mass of interstellar hydrogen mixed with the ejecta of a single supernova which is larger for larger explosion energy and the larger Str\"omgren radius of the progenitor.

We present the results of 3D Newtonian SPH simulations of the merger of a neutron star binary. The microscopic properties of matter are described by the physical equation of state of Lattimer and Swesty (LS-EOS). To test for the robustness of our results we check the sensitivity to the approximations of our model as well as to the binary system parameters. The main and new result is that for the realistic LS-EOS, depending on the initial spin, between 4e-3 and 4e-2 solar masses of material become unbound. If, as suggested, large parts of this matter consist of r-process nuclei, neutron star mergers could account for the whole observed r-process material in the Galaxy.

The influence of the initial composition and structure of the exploding white dwarf on the nucleosynthesis and structure of Type Ia Supernovae has been studied. The progenitor structures are based on detailed stellar evolutionary tracks for stars in the mass range between 1 to 9 $M_\odot$ using the state of the art code FRANEC. The calculations of the thermonuclear explosions are based on a set of delayed detonation models which give a good account of the optical and infrared light curves and of the spectral evolution. Our code solves the hydrodynamical equations explicitly by the piecewise parabolic method. Nuclear burning is taken into account using an extended network of 218 nuclei. In principle, our calculations allow the observed spectra and light curve to be linked to the progenitor. Moreover, our study is relevant to estimate potential evolution in the progenitor population at cosmological time scales.

Having examined the application of quasi-equilibrium to hydrostatic silicon burning in Paper I of this series, Hix & Thielemann (1996), we now turn our attention to explosive silicon burning. Previous authors have shown that for material which is heated to high temperature by a passing shock and then cooled by adiabatic expansion, the results can be divided into three broad categories; \emph{incomplete burning}, \emph{normal freezeout} and \emph{$\alpha$-rich freezeout}, with the outcome depending on the temperature, density and cooling timescale. In all three cases, we find that the important abundances obey quasi-equilibrium for temperatures greater than approximately 3 GK, with relatively little nucleosynthesis occurring following the breakdown of quasi-equilibrium. We will show that quasi-equilibrium provides better abundance estimates than global nuclear statistical equilibrium, even for normal freezeout and particularly for $\alpha$-rich freezeout. We will also examine the accuracy with which the final nuclear abundances can be estimated from quasi-equilibrium.

Recent calculations place questions on the ability of supernovae to produce r-process nuclei in the correct amounts. We present results from 3D Newtonian SPH calculations of mergers of equal mass neutron star binaries. We find the amounts of ejected mass to be very promising for an explanation of the observed r-process abundances. Preliminary abundance calculations show that practically all the material is subject to r-process conditions. The calculated abundance patterns reproduce the basic features of the solar r-process abundances very well. Especially the peak around A= 195 is easily reproduced without any tuning of the initial entropies. Thus, eventually all the observed r-process material could be explained by mass ejection during neutron star mergers.

Examination of the process of silicon burning, the burning stage that leads to the production of the iron peak nuclei, reveals that the nuclear evolution is dominated by large groups of nuclei in mutual equilibrium. These quasi-equilibrium (QSE) groups form well in advance of the global Nuclear Statistical Equilibrium (NSE). We present an improved ``minimal'' nuclear network, which takes advantage of quasi-equilibrium in order to calculate the nuclear evolution and energy generation while further reducing the computational cost compared to a conventional \alpha-chain network. During silicon burning, the resultant \emph{QSE-reduced} \alpha network is twice as fast as the full \alpha network it replaces and requires the tracking of only half as many abundance variables, without significant loss of accuracy. When the QSE-reduced \alpha network is used in combination with a conventional \alpha network stretching from He to Si, the combined \alpha7 network provides an accurate approximation for all of the burning stages from He burning to NSE, while tracking only 7 abundances. These reductions in computational cost and the number of species evolved make the \alpha7 network well suited for inclusion within hydrodynamic simulations, particularly those in multi-dimension.

An improved code for the calculation of astrophysical reaction rates in the statistical model is presented. It includes the possibility to study isospin effects. Such effects heavily affect capture rates involving self-conjugate nuclei and may also be found in reactions on other intermediate and heavy targets.

This review concentrates on nucleosynthesis processes in general and their applications to massive stars and supernovae. A brief initial introduction is given to the physics in astrophysical plasmas which governs composition changes. We present the basic equations for thermonuclear reaction rates and nuclear reaction networks. The required nuclear physics input for reaction rates is discussed, i.e. cross sections for nuclear reactions, photodisintegrations, electron and positron captures, neutrino captures, inelastic neutrino scattering, and beta-decay half-lives. We examine especially the present state of uncertainties in predicting thermonuclear reaction rates, while the status of experiments is discussed by others in this volume (see M. Wiescher). It follows a brief review of hydrostatic burning stages in stellar evolution before discussing the fate of massive stars, i.e. the nucleosynthesis in type II supernova explosions (SNe II). Except for SNe Ia, which are explained by exploding white dwarfs in binary stellar systems (which will not be discussed here), all other supernova types seem to be linked to the gravitational collapse of massive stars (M$>$8M$_\odot$) at the end of their hydrostatic evolution. SN1987A, the first type II supernova for which the progenitor star was known, is used as an example for nucleosynthesis calculations. Finally, we discuss the production of heavy elements in the r-process up to Th and U and its possible connection to supernovae.

The influence of the initial composition of the exploding white dwarf on the nucleosynthesis, light curves and spectra of Type Ia supernovae has been studied in order to evaluate the size of evolutionary effects on cosmological time scales, how the effects can be recognized and how one may be able to correct for them. The calculations are based on a set of delayed detonation models which give a good account of the optical and infrared light curves and of the spectral evolution. The explosions and light curves are calculated using a one- dimensional Lagrangian radiation-hydro code including a nuclear network. NLTE- spectra are computed for various epochs using the structure resulting from the light curve code. The following questions are addressed : What do we learn about the progenitor evolution and its metallicity? What are the systematic effects for the determination of the cosmological parameters $\Omega_M$ and $\Lambda$ and how can we recognize this potential 'pitfalls' and correct for evolutionary effects?

The prediction of cross sections for nuclei far off stability is crucial in the field of nuclear astrophysics. We discuss the model mostly employed for such calculations: the statistical model (Hauser-Feshbach). Special emphasis is put on the uncertainties arising from nuclear level density descriptions and an improved global description is presented. Furthermore, criteria for the applicability of the statistical model are investigated and a "map" for the applicability of the model to reactions of stable and unstable nuclei with neutral and charged particles is given.

Among the major uncertainties involved in the Chandrasekhar mass models for Type Ia supernovae are the companion star of the accreting white dwarf (or the accretion rate that determines the carbon ignition density) and the flame speed after ignition. We present nucleosynthesis results from relatively slow deflagration (1.5 - 3 % of the sound speed) to constrain the rate of accretion from the companion star. Because of electron capture, a significant amount of neutron-rich species such as ^{54}Cr, ^{50}Ti, ^{58}Fe, ^{62}Ni, etc. are synthesized in the central region. To avoid the too large ratios of ^{54}Cr/^{56}Fe and ^{50}Ti/^{56}Fe, the central density of the white dwarf at thermonuclear runaway must be as low as \ltsim 2 \e9 \gmc. Such a low central density can be realized by the accretion as fast as $\dot M \gtsim 1 \times 10^{-7} M_\odot yr^{-1}$. These rapidly accreting white dwarfs might correspond to the super-soft X-ray sources.

Presupernova evolution and explosive nucleosynthesis in massive stars for main-sequence masses from 13 $M_\odot$ to 70 $M_\odot$ are calculated. We examine the dependence of the supernova yields on the stellar mass, $^{12}C(\alpha, \gamma) ^{16}O}$ rate, and explosion energy. The supernova yields integrated over the initial mass function are compared with the solar abundances.

Oxygen to iron abundance ratios of metal-poor stars provide information on nucleosynthesis yields from massive stars which end in Type II supernova explosions. Using a standard model of chemical evolution of the Galaxy we have reproduced the solar neighborhood abundance data and estimated the oxygen and iron yields of genuine SN II origin. The estimated yields are compared with the theoretical yields to derive the relation between the lower and upper mass limits in each generation of stars and the IMF slope. Independently of this relation, we furthermore derive the relation between the lower mass limit and the IMF slope from the stellar mass to light ratio in the solar neighborhood. These independent relations unambiguously determine the upper mass limit of $m_u=50 \pm 10 M_sun$ and the IMF slope index of 1.3 - 1.6 above 1 M_sun. This upper mass limit corresponds to the mass beyond which stars end as black holes without ejecting processed matter into the interstellar medium. We also find that the IMF slope index below 0.5 M_sun cannot be much shallower than 0.8.

Hot (explosive) hydrogen burning or the Rapid Proton Capture Process (rp-process) occurs in a number of astrophysical environments. Novae and X-ray bursts are the most prominent ones, but accretion disks around black holes and other sites are candidates as well. The expensive and often multidimensional hydro calculations for such events require an accurate prediction of the thermonuclear energy generation, while avoiding full nucleosynthesis network calculations. In the present investigation we present an approximation scheme applicable in a temperature range which covers the whole range of all presently known astrophysical sites. It is based on the concept of slowly varying hydrogen and helium abundances and assumes a kind of local steady flow by requiring that all reactions entering and leaving a nucleus add up to a zero flux. This scheme can adapt itself automatically and covers situations at low temperatures, characterized by a steady flow of reactions, as well as high temperature regimes where a $(p,\gamma)-(\gamma,p)$-equilibrium is established. In addition to a gain of a factor of 15 in computational speed over a full network calculation, and an energy generation accurate to more than 15 %, this scheme also allows to predict correctly individual isotopic abundances. Thus, it delivers all features of a full network at a highly reduced cost and can easily be implemented in hydro calculations.

Nuclear Astrophysics requires the knowledge of reaction rates over a wide range of nuclei and temperatures. In recent calculations the nuclear level density - as an important ingredient to the statistical model (Hauser-Feshbach) - has shown the highest uncertainties. In a back-shifted Fermi-gas formalism utilizing an energy-dependent level density parameter and employing microscopic corrections from a recent FRDM mass formula, we obtain a highly improved fit to experimental level densities. The resulting level density is used for determining criteria for the applicability of the statistical model on neutron-induced reactions.

The prediction of cross sections for nuclei far off stability is crucial in the field of nuclear astrophysics. In recent calculations the nuclear level density -- as an important ingredient to the statistical model (Hauser-Feshbach) -- has shown the highest uncertainties. We present a global parametrization of nuclear level densities based on the back-shifted Fermi-Gas formalism. Employment of an energy-dependent level density parameter $a$ and microscopic corrections from a recent FRDM mass formula by M\"oller et al.\ leads to a highly improved fit of level densities at the neutron-separation energy in the mass range $20\le A \le 245$. The importance of using proper microscopic corrections from mass formulae is emphasized. The resulting level description is well suited for astrophysical applications. The level density can also provide clues to the applicability of the statistical model which is only correct for a high density of excited states. Using the above description one can derive a ``map'' for the applicability of the model for reactions of stable and unstable nuclei with neutral and charged particles.

From the spectra and light curves it is clear that SNIa are thermonuclear explosions of white dwarfs. However, details of the explosion are highly under debate. Here, we present detailed models which are consistent with respect to the explosion mechanism, the optical and infrared light curves (LC), and the spectral evolution. This leaves the description of the burning front and the structure of the white dwarf as the only free parameters. The explosions are calculated using one-dimensional Lagrangian codes including nuclear networks. Subsequently, optical and IR-LCs are constructed. Detailed NLTE-spectra are computed for several instants of time using the density, chemical and luminosity structure resulting from the LCs. The general methods and critical tests are presented (sect. 2). Different models for the thermonuclear explosion are discussed including detonations, deflagrations, delayed detonations, pulsating delayed detonations (PDD) and helium detonations (sect.3). Comparisons between theoretical and observed LCs and spectra provide an insight into details of the explosion and nature of the progenitor stars (sect. 4/5). We try to answer several related questions. Are subluminous SNe Ia a group different from `normal' SN Ia (sect. 5)? Can we understand observed properties of the LCs and spectra (sect. 4)? What do we learn about the progenitor evolution and its metallicity (sect. 3, Figs. 4,5)? Do successful SN~Ia models depend on the type of the host galaxy (Table 2)? Using both the spectral and LC information, theoretical models allow for a determination of the Hubble constant independent from `local' distance indicators such as delta-Cephei stars. Ho is found to be 67 +- 9km/s/Mpc and, from SN1988U, qo equals 0.7 +- 1. within 95 percent confidence levels.

As the ultimate stage of stellar nucleosynthesis, and the source of the iron peak nuclei, silicon burning is important to our understanding of the evolution of massive stars and supernovae. Our reexamination of silicon burning, using results gleaned from simulation work done with a large nuclear network (299 nuclei and more than 3000 reactions) and from independent calculations of equilibrium abundance distributions, offers new insights into the quasi-equilibrium mechanism and the approach to nuclear statistical equilibrium. We find that the degree to which the matter has been neutronized is of great importance, not only to the final products but also to the rate of energy generation and the membership of the quasi-equilibrium groups. A small increase in the global neutronization results in much larger free neutron fluences, increasing the abundances of more neutron-rich nuclei. As a result, incomplete silicon burning results in neutron richness among the isotopes of the iron peak much larger than the global neutronization would indicate. Finally, we briefly discuss the limitations and pitfalls of models for silicon burning currently employed within hydrodynamic models. In a forthcoming paper we will present a new approximation to the full nuclear network which preserves the most important features of the large nuclear network calculations at a significant improvement in computational speed. Such improved methods are ideally suited for hydrodynamic calculations which involve the production of iron peak nuclei, where the larger network calculation proves unmanageable.

We present 3D Newtonian simulations of the coalescence of two neutron stars, using a Smoothed Particle Hydrodynamics (SPH) code. We begin the simulations with the two stars in a hard, circular binary, and have them spiral together as angular momentum is lost through gravitational radiation at the rate predicted by modeling the system as two point masses. We model the neutron stars as hard polytropes ($\gamma=2.4$) of equal mass, and investigate the effect of the initial spin of the two stars on the coalescence. The process of coalescence, from initial contact to the formation of an axially symmetric object, takes only a few orbital periods. Some of the material from the two neutron stars is shed, forming a thick disk around the central, coalesced object. The mass of this disk is dependent on the initial neutron star spins; higher spin rates resulting in greater mass loss, and thus more massive disks. For spin rates that are most likely to be applicable to real systems, the central coalesced object has a mass of $2.4M_\odot$, which is tantalizingly close to the maximum mass allowed by any neutron star equation of state for an object that is supported in part by rotation. Using a realistic nuclear equation of state, we estimate the temperature of the material after the coalescence. We find that the central object is at a temperature of $\sim 10$MeV, whilst the disk is heated by shocks to a temperature of 2-4MeV.