The most efficient axion production mechanism in a supernova (SN) core is the nucleon-nucleon bremsstrahlung. This process has been often modeled at the level of the vacuum one-pion exchange (OPE) approximation. Starting from this naive recipe, we revise the calculation including systematically different effects, namely a non-vanishing mass for the exchanged pion, the contribution from the two-pions exchange, effective in-medium nucleon masses and multiple nucleon scatterings. Moreover, we allow for an arbitrary degree of nucleon degeneracy. A self consistent treatment of the axion emission rate including all these effects is currently missing. The aim of this work is to provide such an analysis. Furthermore, we demonstrate that the OPE potential with all the previous corrections gives rise to similar results as the on-shell T-matrix, and is therefore well justified for our and similar studies. We find that the axion emissivity is reduced by over an order of magnitude with respect to the basic OPE calculation, after all these effects are accounted for. The implications for the axion mass bound and the impact for the next generation experimental axion searches is also discussed.

Neutrinos produced during a supernova explosion induce reactions on abundant nuclei in the outer stellar shells and contribute in this way to the synthesis of the elements in the Universe. This neutrino nucleosynthesis process has been identified as an important contributor to the origin of $^7$Li, $^{11}$B,$^{19}$F, $^{138}$La, and $^{180}$Ta, but also to the long-lived radionuclides $^{22}$Na and $^{26}$Al, which are both key isotopes for $\gamma$-ray astronomy. The manuscript summarizes the recent progress achieved in simulations of neutrino nucleosynthesis.

The kilonova emission observed following the binary neutron star merger event GW170817 provided the first direct evidence for the synthesis of heavy nuclei through the r-process. The late-time transition in the spectral energy distribution to near-infrared wavelengths was interpreted as indicating the production of lanthanide nuclei, with atomic mass number A>140. However, compelling evidence for the presence of heavier third-peak r-process elements (e.g., gold, platinum) or translead nuclei remains elusive. At early times (~days) most of the r-process heating arises from a large statistical ensemble of beta-decays, which thermalize efficiently while the ejecta is still dense, generating a heating rate that is reasonably approximated by a single power-law. However, at later times (weeks to months), the decay energy input can possibly be dominated by a discrete number of alpha-decays, 223Ra (half-life t_{1/2}=11.43d), 225Ac (t_{1/2}=10.0d, following the beta-decay of 225Ra with t_{1/2}=14.9d), and the fissioning isotope 254Cf (t_{1/2}=60.5d), which liberate more energy per decay and thermalize with greater efficiency than beta-decay products. Late-time nebular observations of kilonovae which constrain the radioactive power provide the potential to identify signatures of these individual isotopes, thus confirming the production of heavy nuclei. In order to constrain the bolometric light to the required accuracy, multi-epoch & wide-band observations are required with sensitive instruments like the James Webb Space Telescope. In addition, by comparing the nuclear heating rate obtained with an abundance distribution following the Solar r abundance pattern, to the bolometric lightcurve of AT2017gfo, we find that the yet-uncertain r abundance of 72Ge plays a decisive role in powering the lightcurve, if one assumes that GW170817 has produced a full range of the Solar r abundances down to A~70.

Muons can be created in nascent neutron stars (NSs) due to the high electron chemical potentials and the high temperatures. Because of their relatively lower abundance compared to electrons, their role has so far been ignored in numerical simulations of stellar core collapse and NS formation. However, the appearance of muons softens the NS equation of state, triggers faster NS contraction and thus leads to higher luminosities and mean energies of the emitted neutrinos. This strengthens the postshock heating by neutrinos and can facilitate explosions by the neutrino-driven mechanism.

The ultimate understanding of Type Ia Supernovae diversity is one of the most urgent issues to exploit thermonuclear explosions of accreted White Dwarfs (WDs) as cosmological yardsticks. In particular, we investigate the impact of the progenitor system metallicity on the physical and chemical properties of the WD at the explosion epoch. We analyze the evolution of CO WDs through the accretion and simmering phases by using evolutionary models based on time-dependent convective mixing and an extended nuclear network including the most important electron captures, beta decays and URCA processes. We find that, due to URCA processes and electron-captures, the neutron excess and density at which the thermal runaway occurs are substantially larger than previously claimed. Moreover, we find that the higher the progenitor metallicity, the larger the neutron excess variation during the accretion and simmering phases and the higher the central density and the convective velocity at the explosion. Hence, the simmering phase acts as an amplifier of the differences existing in SNe Ia progenitors. When applying our results to the neutron excess estimated for the Tycho and Kepler young Supernova remnants, we derive that the metallicity of the progenitors should be in the range Z=0.030-0.032, close to the average metallicity value of the thin disk of the Milky Way. As the amount of ${^{56}}$Ni produced in the explosion depends on the neutron excess and central density at the thermal runaway, our results suggest that the light curve properties depend on the progenitor metallicity.

We study the optical and near-infrared luminosities and detectability of radioactively powered electromagnetic transients ('macronovae') occuring in the aftermath of binary neutron star and neutron star black hole mergers. We explore the transients that result from the dynamic ejecta and those from different types of wind outflows. Based on full nuclear network simulations we calculate the resulting light curves in different wavelength bands. We scrutinize the robustness of the results by comparing a) two different nuclear reaction networks and b) two macronova models. We explore in particular how sensitive the results are to the production of alpha-decaying trans-lead nuclei. We compare two frequently used mass models: the Finite-Range Droplet Model (FRDM) and the nuclear mass model of Duflo and Zuker (DZ31). We find that the abundance of alpha-decaying trans-lead nuclei has a significant impact on the observability of the resulting macronovae. For example, the DZ31 model yields considerably larger abundances resulting in larger heating rates and thermalization efficiencies and therefore predicts substantially brighter macronova transients. We find that the dynamic ejecta from NSNS models can reach peak K-band magnitudes in excess of -15 while those from NSBH cases can reach beyond -16. Similar values can be reached by some of our wind models. Several of our models (both wind and dynamic ejecta) yield properties that are similar to the transient that was observed in the aftermath of the short GRB 130603B. We further explore the expected macronova detection frequencies for current and future instruments such as VISTA, ZTF and LSST.

Neutrino interactions with matter play an important role in determining the nucleosynthesis outcome in explosive astrophysical environments such as core-collapse supernovae or mergers of compact objects. In this article, we first discuss our recent work on the importance of studying the time evolution of collective neutrino oscillations among active flavors in determining their effects on nucleosynthesis. We then consider the possible active-sterile neutrino mixing and demonstrate the need of a consistent approach to evolve neutrino flavor oscillations, matter composition, and the hydrodynamics when flavor oscillations can happen very deep inside the supernovae.

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.

In this paper, we explore the effects of neutrino flavor oscillations on supernova nucleosynthesis and on the neutrino signals. Our study is based on detailed information about the neutrino spectra and their time evolution from a spherically-symmetric supernova model for an 18 M progenitor. We find that collective neutrino oscillations are not only sensitive to the detailed neutrino energy and angular distributions at emission, but also to the time evolution of both the neutrino spectra and the electron density profile. We apply the results of neutrino oscillations to study the impact on supernova nucleosynthesis and on the neutrino signals from a Galactic supernova. We show that in our supernova model, collective neutrino oscillations enhance the production of rare isotopes 138La and 180Ta but have little impact on the nu p-process nucleosynthesis. In addition, the adiabatic MSW flavor transformation, which occurs in the C/O and He shells of the supernova, may affect the production of light nuclei such as 7Li and 11B. For the neutrino signals, we calculate the rate of neutrino events in the Super-Kamiokande detector and in a hypothetical liquid argon detector. Our results suggest the possibility of using the time profiles of the events in both detectors, along with the spectral information of the detected neutrinos, to infer the neutrino mass hierarchy.

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

We study the impact of neutrino-pair production from the de-excitation of highly excited heavy nuclei on core-collapse supernova simulations, following the evolution up to several 100 ms after core bounce. Our study is based on the AGILE-Boltztran supernova code, which features general relativistic radiation hydrodynamics and accurate three-flavor Boltzmann neutrino transport in spherical symmetry. In our simulations the nuclear de-excitation process is described in two different ways. At first we follow the approach proposed by Fuller and Meyer [Astrophys. J. 376,701 (1991)], which is based on strength functions derived in the framework of the nuclear Fermi-gas model of non-interacting nucleons. Secondly, we parametrize the allowed and forbidden strength distributions in accordance with measurements for selected nuclear ground states. We determine the de-excitation strength by applying the Brink hypothesis and detailed balance. For both approaches, we find that nuclear de-excitation has no effect on the supernova dynamics. However, we find that nuclear de-excitation is the leading source for the production of electron antineutrinos as well as heavy-lepton flavor (anti)neutrinos during the collapse phase. At sufficiently high densities, the associated neutrino spectra are influenced by interactions with the surrounding matter, making proper simulations of neutrino transport important for the determination of the neutrino-energy loss rate. We find that even including nuclear de-excitations, the energy loss during the collapse phase is overwhelmingly dominated by electron neutrinos produced by electron captures.

The stellar mass range 8<M/Mo<12 corresponds to the most massive AGB stars and the most numerous massive stars. It is host to a variety of supernova progenitors and is therefore very important for galactic chemical evolution and stellar population studies. In this paper, we study the transition from super-AGB star to massive star and find that a propagating neon-oxygen burning shell is common to both the most massive electron capture supernova (EC-SN) progenitors and the lowest mass iron-core collapse supernova (FeCCSN) progenitors. Of the models that ignite neon burning off-center, the 9.5Mo model would evolve to an FeCCSN after the neon-burning shell propagates to the center, as in previous studies. The neon-burning shell in the 8.8Mo model, however, fails to reach the center as the URCA process and an extended (0.6 Mo) region of low Ye (0.48) in the outer part of the core begin to dominate the late evolution; the model evolves to an EC-SN. This is the first study to follow the most massive EC-SN progenitors to collapse, representing an evolutionary path to EC-SN in addition to that from SAGB stars undergoing thermal pulses. We also present models of an 8.75Mo super-AGB star through its entire thermal pulse phase until electron captures on 20Ne begin at its center and of a 12Mo star up to the iron core collapse. We discuss key uncertainties and how the different pathways to collapse affect the pre-supernova structure. Finally, we compare our results to the observed neutron star mass distribution.

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

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

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

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

We study proton rich nucleosynthesis in windlike outflows from gamma-ray bursts accretion disks with the aim to determine if such outflows are a site of the vp-process. The efficacy of this vp-process depends on thermodynamic and hydrodynamic factors. We discuss the importance of the entropy of the material, the outflow rate, the initial ejection point and accretion rate of the disk. In some cases the vp-process pushes the nucleosynthesis out to A~100 and produces light p-nuclei. However, even when these nuclei are not produced, neutrino induced interactions can significantly alter the abundance pattern and cannot be neglected.

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

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

We have performed shell-model calculations of the half-lives and neutron-branching probabilities of the r-process waiting point nuclei at the magic neutron number N=82. These new calculations use a larger model space than previous shell model studies and an improved residual interaction which is adjusted to recent spectroscopic data around A=130. Our shell-model results give a good account of all experimentally known half-lives and $Q_\beta$-values for the N=82 r-process waiting point nuclei. Our half-life predictions for the N=82 nuclei with Z=42--46 agree well with recent estimates based in the energy-density functional method.

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

Based on the shell model for Gamow-Teller and the Random Phase Approximation for forbidden transitions, we have calculated reaction rates for inelastic neutrino-nucleus scattering (INNS) under supernova (SN) conditions, assuming a matter composition given by Nuclear Statistical Equilibrium. The rates have been incorporated into state-of-the-art stellar core-collapse simulations with detailed energy-dependent neutrino transport. While no significant effect on the SN dynamics is observed, INNS increases the neutrino opacities noticeably and strongly reduces the high-energy tail of the neutrino spectrum emitted in the neutrino burst at shock breakout. Relatedly the expected event rates for the observation of such neutrinos by earthbound detectors are reduced by up to about 60%.

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.

Advances in our understanding and the modeling of stellar core-collapse and supernova explosions over the past 15 years are reviewed, concentrating on the evolution of hydrodynamical simulations, the description of weak interactions and nuclear equation of state effects, and new insights into the nucleosynthesis occurring in the early phases of the explosion, in particular the neutrino-p process. The latter is enabled by the proton-richness of the early ejecta, which was discovered because of significant progress has been made in the treatment of neutrino transport and weak interactions. This progress has led to a new generation of sophisticated Newtonian and relativistic hydrodynamics simulations in spherical symmetry. Based on these, it is now clear that the prompt bounce-shock mechanism is not the driver of supernova explosions, and that the delayed neutrino-heating mechanism can produce explosions without the aid of multi-dimensional processes only if the progenitor star has an ONeMg core inside a very dilute He-core, i.e., has a mass in the 8--10 solar mass range. Hydrodynamic instabilities of various kinds have indeed been recognized to occur in the supernova core and to be of potential importance for the explosion. Neutrino-driven explosions, however, have been seen in two-dimensional simulations with sophisticated neutrino transport so far only when the star has a small iron core and low density in the surrounding shells as being found in stars near 10--11 solar masses. The explosion mechanism of more massive progenitors is still a puzzle. It might involve effects of three-dimensional hydrodynamics or might point to the relevance of rapid rotation and magnetohydrodynamics, or to still incompletely explored properties of neutrinos and the high-density equation of state.

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 review the nuclear physics input necessary for the study of the collapse of massive stars as precursor to supernova explosions. Recent theoretical advances for the calculation of the relevant weak-interaction processes and their influence in the collapse and postbounce dynamics are discussed. Future improvements are expected to come with the advent of radioactive ion-beam facilities.

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.

We study neutrino process nucleosynthesis in massive stars using newly calculated cross sections, an expanded reaction network, and complete and self-consistent models of the progenitor star. We reevaluate the production of light isotopes from abundant progenitors as well as that of rare, heavy, proton-rich isotopes. In particular, new results are given for B11, F19, La138, and Ta180. The production of these isotopes places limits on neutrino spectrum and oscialltions.

Supernova simulations to date have assumed that during core collapse electron captures occur dominantly on free protons, while captures on heavy nuclei are Pauli-blocked and are ignored. We have calculated rates for electron capture on nuclei with mass numbers A=65-112 for the temperatures and densities appropriate for core collapse. We find that these rates are large enough so that, in contrast to previous assumptions, electron capture on nuclei dominates over capture on free protons. This leads to significant changes in core collapse simulations.

A simple description for obtaining the parity distribution of nuclear levels in the pf + g9/2 shell as a function of excitation energy was recently derived. We implement this in a global nuclear level density model. In the framework of the statistical model, cross sections and astrophysical reaction rates are calculated in the Fe region and compared to rates obtained with the common assumption of an equal distribution of parities. We find considerable differences, especially for reactions involving particles in the exit channel.

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.

Recent experimental data and progress in nuclear structure modeling have lead to improved descriptions of astrophysically important weak-interaction processes. The review discusses these advances and their applications to hydrostatic solar and stellar burning, to the slow and rapid neutron-capture processes, to neutrino nucleosynthesis, and to explosive hydrogen burning. Special emphasis is given to the weak-interaction processes associated with core-collapse supernovae. Despite some significant progress, important improvements are still warranted. Such improvements are expected to come from future radioactive ion-beam facilities.

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.

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.

Using the recent shell model evaluation of stellar weak interaction rates we have calculated the neutrino spectra arising from electron capture on pf-shell nuclei under presupernova conditions. We present a simple parametrization of the spectra which allows for an easy implementation into collapse simulations. We discuss that the explicit consideration of thermal ensembles in the parent nucleus broadens the neutrino spectra and results in larger average neutrino energies. The capture rates and neutrino spectra can be easily modified to account for phase space blocking by neutrinos which becomes increasingly important during the final stellar collapse.

Recent shell-model calculations of weak-interaction rates for nuclei in the mass range A = 45 - 65 have resulted in substantial revisions to the hitherto standard set of Fuller, Fowler, & Newman (FFN). In particular, key electron-capture rates, such as that for Co60 are much smaller. We consider here the effects of these revised rates on the presupernova (post-oxygen burning) evolution of massive stars in the mass range 11 to 40 M_sun. Moreover, we include, for the first time in models by our group, the effects of modern rates for beta-decay in addition to electron capture and positron emission. Values for the central electron mole number at the time of iron core collapse in the new models are typically larger, by delta Y_e = 0.005 to 0.015, than those of Woosley & Weaver 1995, with a tendency for the more massive models to display larger differences. About half of this change is a consequence of including beta-decay; the other half, result of the smaller rates for electron capture. Unlike what might be expected solely on basis of the larger Y_e's, the new iron core masses are systematically smaller owing to a decrease in the entropy in the outer iron core. The changes in iron core mass range from zero to 0.1 M_sun. We also observe, as predicted by Aufderheide et al. (1994), a tendency towards beta-equilibrium just prior to the collapse of the core, and the subsequent loss of that equilibrium as core collapse proceeds. We discuss the key weak reaction rates, both beta-decay and electron-capture, responsible for the evolution of Y_e and make suggestions for future measurements.

Improved values for stellar weak interaction rates have been recently calculated based upon a large shell model diagonalization. Using these new rates (for both beta decay and electron capture), we have examined the presupernova evolution of massive stars in the range 15-40 Msun. Comparing our new models with a standard set of presupernova models by Woosley and Weaver, we find significantly larger values for the electron-to-baryon ratio Ye at the onset of collapse and iron core masses reduced by approximately 0.1 Msun. The inclusion of beta-decay accounts for roughly half of the revisions, while the other half is a consequence of the improved nuclear physics. These changes will have important consequences for nucleosynthesis and the supernova explosion mechanism.

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.

Based on large-scale shell model calculations we have determined the electron capture, positron capture and beta-decay rates on more than 100 nuclei in the mass range A=45-65. The rates are given for densities \rho Y_e =10^7-10^{10} mol/cm^3 and temperatures T=10^9-10^{10} K and hence are relevant for both types of supernovae (Type Ia and Type II). The shell model electron capture rates are significantly smaller than currently assumed. For proton-to-baryon ratios Y_e=0.42-0.46 mol/g, the beta-decay rates are faster than the electron capture rates during the core collapse of a massive star.

We have performed large-scale shell-model calculations of the half-lives and neutron-branching probabilties of the r-process waiting point nuclei at the magic neutron number N=82. We find good agreement with the measured half-lives of 129Ag and 130Cd. Our shell-model half-lives are noticeably shorter than those currently adopted in r-process simulations. Our calculation suggests that 130Cd is not produced in beta-flow equilibrium with the other N=82 isotones on the r-process path.

A measurement of the 56Ni cosmic ray abundance has been discussed as a possible tool to determine the acceleration time scale of relativistic particles in cosmic rays. This conjecture will depend on the halflife of totally ionized 56Ni which can only decay by higher-order forbidden transitions. We have calculated this halflife within large-scale shell model calculations and find t_{1/2} \approx 4 \times 10^4 years, only slightly larger than the currently available experimental lower limit, but too short for 56Ni to serve as a cosmic ray chronometer.

We calculate supernova electron capture and beta decay rates for various pf-shell nuclei using large-scale shell model techniques. We show that the centroid of the Gamow-Teller strength distribution has been systematically misplaced in previous rate estimates. Our total electron capture rates are significantly smaller than currently adopted in core collapse calculations, while the total beta decay rates change less. Our calculation shows that for electron-to-baryon ratios Ye=0.42-0.46 beta decay rates are larger than the competing electron capture rates.

At densities between 10^8 and 10^10} g/cm^3 electron capture in a presupernova collapse is believed to mainly occur on odd-odd nuclei. We have derived the rates for six of the most important electron capturing nuclei, {54,56,58}Mn and {56,58,60}Co, based on calculations of the Gamow-Teller strength distributions for the ground states and first excited states. These calculations have been performed by shell model diagonalization in the pf shell using a recently modified version of the KB3 interaction. The shell model rates are noticeably smaller than the presently adopted rates as the latter have been derived by placing the Gamow-Teller (GT) resonance at too low excitation energies.

We have calculated the Gamow-Teller strength distributions for the ground states and first excited states in 55Co and 56Ni. These calculations have been performed by shell model diagonalization in the pf shell using the KB3 interaction. The Gamow-Teller distributions are used to calculate the electron capture rates for typical presupernova conditions. Our 55Co rate is noticeably smaller than the presently adopted rate as it is dominated by weak low-lying transitions rather than the strong Gamow-Teller (GT) resonance which is located at a higher excitation energy in the daughter than usually parametrized. Although our 56Ni rate agrees with the presently adopted rate, we do not confirm the conventional parametrization of the GT centroid. Our results support general trends suggested on the basis of shell model Monte Carlo calculations.

The nucleus 54Mn has been observed in cosmic rays. In astrophysical environments it is fully stripped of its atomic electrons and its decay is dominated by the beta- branch to the 54Fe ground state. Application of 54Mn based chronometer to study the confinement of the iron group cosmic rays requires knowledge of the corresponding halflife, but its measurement is impossible at the present time. However, the branching ratio for the related beta+ decay of 54Mn was determined recently. We use the shell model with only a minimal truncation and calculate both beta+ and beta- decay rates of 54Mn. Good agreement for the beta+ branch suggests that the calculated partial halflife of the beta- decay, (4.94 \pm 0.06) x 10^5 years, should be reliable. However, this halflife is noticeably shorter than the range 1-2 x 10^6 y indicated by the fit based on the 54Mn abundance in cosmic rays. We also evaluate other known unique second forbidden beta decays from the nuclear p and sd shells (10Be, 22Na, and two decay branches of 26Al) and show that the shell model can describe them with reasonable accuracy as well.