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.

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

We introduce a new algorithm for the calculation of multidimensional optical depths in approximate radiative transport schemes, equally applicable to neutrinos and photons. Motivated by (but not limited to) neutrino transport in three-dimensional simulations of core-collapse supernovae and neutron star mergers, our method makes no assumptions about the geometry of the matter distribution, apart from expecting optically transparent boundaries. Based on local information about opacities, the algorithm figures out an escape route that tends to minimize the optical depth without assuming any pre-defined paths for radiation. Its adaptivity makes it suitable for a variety of astrophysical settings with complicated geometry (e.g., core-collapse supernovae, compact binary mergers, tidal disruptions, star formation, etc.). We implement the MODA algorithm into both a Eulerian hydrodynamics code with a fixed, uniform grid and into an SPH code where we make use a tree structure that is otherwise used for searching neighbours and calculating gravity. In a series of numerical experiments, we compare the MODA results with analytically known solutions. We also use snapshots from actual 3D simulations and compare the results of MODA with those obtained with other methods such as the global and local ray-by-ray method. It turns out that MODA achieves excellent accuracy at a moderate computational cost. In an appendix we also discuss implementation details and parallelization strategies.

By implementing widely-used equations of state (EOS) from Lattimer & Swesty (LS) and H. Shen et al. (SHEN) in core-collapse supernova simulations, we explore possible impacts of these EOS on the post-bounce dynamics prior to the onset of neutrino-driven explosions. Our spherically symmetric (1D) and axially symmetric (2D) models are based on neutrino radiation hydrodynamics including spectral transport, which is solved by the isotropic diffusion source approximation. We confirm that in 1D simulations neutrino-driven explosions cannot be obtained for any of the employed EOS. Impacts of the EOS on the post-bounce hydrodynamics are more clearly visible in 2D simulations. In 2D models of a 15 M_sun progenitor using the LS EOS, the stalled bounce shock expands to increasingly larger radii, which is not the case using the SHEN EOS. Keeping in mind that the omission of the energy drain by heavy-lepton neutrinos in the present scheme could facilitate explosions, we find that 2D models of an 11.2 M_sun progenitor produce neutrino-driven explosions for all the EOS under investigation. Models using the LS EOS are slightly more energetic compared to those with the SHEN EOS. The more efficient neutrino heating in the LS models coincides with a higher electron antineutrino luminosity and a larger mass that is enclosed within the gain region. The models based on the LS EOS also show a more vigorous and aspherical downflow of accreting matter to the surface of the protoneutron star (PNS). The accretion pattern is essential for the production and strength of outgoing pressure waves, that can push in turn the shock to larger radii and provide more favorable conditions for the explosion. [abbreviated]

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.

By performing a series of one- and two-dimensional (1-, 2D) hydrodynamic simulations with spectral neutrino transport, we study possible impacts of collective neutrino oscillations on the dynamics of core-collapse supernovae. To model the spectral swapping which is one of the possible outcome of the collective neutrino oscillations, we parametrize the onset time when the spectral swap begins, the radius where the spectral swap occurs, and the threshold energy above which the spectral interchange between heavy-lepton neutrinos and electron/anti-electron neutrinos takes place, respectively. By doing so, we systematically study how the neutrino heating enhanced by the spectral swapping could affect the shock evolution as well as the matter ejection. We also investigate the progenitor dependence using a suite of progenitor models (13, 15, 20, and 25 $M_\odot$). We find that there is a critical heating rate induced by the spectral swapping to trigger explosions, which significantly differs between the progenitors. The critical heating rate is generally smaller for 2D than 1D due to the multidimensionality that enhances the neutrino heating efficiency. For the progenitors employed in this paper, the final remnant masses are estimated to range in 1.1-1.5$M_\odot$. For our 2D model of the $15M_\odot$ progenitor, we find a set of the oscillation parameters that could account for strong supernova explosions ($\sim 10^{51}$ erg), simultaneously leaving behind the remnant mass close to $\sim 1.4 M_\odot$.

By performing axisymmetric hydrodynamic simulations of core-collapse supernovae with spectral neutrino transport based on the isotropic diffusion source approximation scheme, we support the assumption that the neutrino-heating mechanism aided by the standing accretion shock instability and convection can initiate an explosion of a 13 $M_{\odot}$ star. Our results show that bipolar explosions are more likely to be associated with models which include rotation. We point out that models, which form a north-south symmetric bipolar explosion, can lead to larger explosion energies than for the corresponding unipolar explosions.

Predictions of the thermodynamic conditions for phase transitions at high baryon densities and large chemical potentials are currently uncertain and largely phenomenological. Neutrino observations of core-collapse supernovae can be used to constrain the situation. Recent simulations of stellar core collapse that include a description of quark matter predict a sharp burst of anti \nu_e several hundred milliseconds after the prompt \nu_e neutronization burst. We study the observational signatures of that anti \nu_e burst at current neutrino detectors - IceCube and Super-Kamiokande. For a Galactic core-collapse supernova, we find that signatures of the QCD phase transition can be detected, regardless of the neutrino oscillation scenario. The detection would constitute strong evidence of a phase transition in the stellar core, with implications for the equation of state at high matter density and the supernova explosion mechanism.

We present the gravitational wave (GW) analysis of an extensive series of 3D MHD core-collapse simulations. Our 25 models are launched from a 15 solar mass progenitor, a spherically symmetric effective general relativistic potential, the Lattimer-Swesty or the Shen equation of state (EoS), and a neutrino parametrisation scheme which is accurate until about 5ms postbounce. For 3 representative models, we also include long-term neutrino physics by means of a leakage scheme. Non- or only slowly rotating models show GW emission due to prompt and proto-neutron star convection, allowing the distinction between the two different nuclear EoS. For moderately or fast rotation rates models, we find, in agreement with recent results, only a type I GW signature at core bounce. Models which are set up with an initial central angular velocity of >~ 2pi rad/s emit GWs due to the low T/|W| dynamical instability during the postbounce phase. Weak B-fields do not notably influence the dynamical evolution of the core and thus the GW emission. However, for strong initial poloidal B-fields (~1e12 G),flux-freezing and field winding leads to conditions where P_{mag}/P_{mat} ~ 1, causing the onset of a jet-like supernova explosion and hence the emission of a type IV GW signal. In contradiction to axisymmetric simulations, we find evidence that nonaxisymmetric fluid modes can counteract or even suppress jet formation for models with strong initial toroidal B-fields. We point out that the inclusion of the deleptonisation during the postbounce phase is an indispensable issue for the quantitative prediction of GWs from core-collapse supernovae, as it can alter the GW amplitude up to a factor of 10 compared to a pure hydrodynamical treatment.

We have performed a set of 11 three-dimensional magnetohydrodynamical core collapse supernova simulations in order to investigate the dependencies of the gravitational wave signal on the progenitor's initial conditions. We study the effects of the initial central angular velocity and different variants of neutrino transport. Our models are started up from a 15 solar mass progenitor and incorporate an effective general relativistic gravitational potential and a finite temperature nuclear equation of state. Furthermore, the electron flavour neutrino transport is tracked by efficient algorithms for the radiative transfer of massless fermions. We find that non- and slowly rotating models show gravitational wave emission due to prompt- and lepton driven convection that reveals details about the hydrodynamical state of the fluid inside the protoneutron stars. Furthermore we show that protoneutron stars can become dynamically unstable to rotational instabilities at T/|W| values as low as ~2 % at core bounce. We point out that the inclusion of deleptonization during the postbounce phase is very important for the quantitative GW prediction, as it enhances the absolute values of the gravitational wave trains up to a factor of ten with respect to a lepton-conserving treatment.

FISH is a fast and simple ideal magneto-hydrodynamics code that scales to ~10 000 processes for a Cartesian computational domain of ~1000^3 cells. The simplicity of FISH has been achieved by the rigorous application of the operator splitting technique, while second order accuracy is maintained by the symmetric ordering of the operators. Between directional sweeps, the three-dimensional data is rotated in memory so that the sweep is always performed in a cache-efficient way along the direction of contiguous memory. Hence, the code only requires a one-dimensional description of the conservation equations to be solved. This approach also enable an elegant novel parallelisation of the code that is based on persistent communications with MPI for cubic domain decomposition on machines with distributed memory. This scheme is then combined with an additional OpenMP parallelisation of different sweeps that can take advantage of clusters of shared memory. We document the detailed implementation of a second order TVD advection scheme based on flux reconstruction. The magnetic fields are evolved by a constrained transport scheme. We show that the subtraction of a simple estimate of the hydrostatic gradient from the total gradients can significantly reduce the dissipation of the advection scheme in simulations of gravitationally bound hydrostatic objects. Through its simplicity and efficiency, FISH is as well-suited for hydrodynamics classes as for large-scale astrophysical simulations on high-performance computer clusters. In preparation for the release of a public version, we demonstrate the performance of FISH in a suite of astrophysically orientated test cases.

Astrophysical observations originate from matter that interacts with radiation or transported particles. We develop a pragmatic approximation in order to enable multi-dimensional simulations with basic spectral radiative transfer when the computational resources are not sufficient to solve the complete Boltzmann transport equation. The distribution function of the transported particles is decomposed into trapped and streaming particle components. Their separate evolution equations are coupled by a source term that converts trapped particles into streaming particles. We determine this source term by requiring the correct diffusion limit. For a smooth transition to the free streaming regime, this 'diffusion source' is limited by the matter emissivity. The resulting streaming particle emission rates are integrated over space to obtain the streaming particle flux. A geometric estimate of the flux factor is used to convert the particle flux to the streaming particle density. The efficiency of the scheme results from the freedom to use different approximations for each particle component. In supernovae, reactions with trapped particles on fast time scales establish equilibria that reduce the number of primitive variables required to evolve the trapped particle component. On the other hand, a stationary-state approximation facilitates the treatment of the streaming particle component. Different approximations may apply in applications to stellar atmospheres, star formation, or cosmological radiative transfer. We compare the isotropic diffusion source approximation with Boltzmann neutrino transport of electron flavour neutrinos in spherically symmetric supernova models and find good agreement. An extension of the scheme to the multi-dimensional case is also discussed.

We present the gravitational wave analysis from rotating (model s15g) and nearly non-rotating (model s15h) 3D MHD core collapse supernova simulations at bounce and the first couple of ten milliseconds afterwards. The simulations are launched from 15M_{\odot} progenitor models stemming from stellar evolution calculations. Gravity is implemented by a spherically symmetric effective general relativistic potential. The input physics uses the Lattimer-Swesty equation of state for hot, dense matter and a neutrino parametrisation scheme that is accurate until the first few ms after bounce. The 3D simulations allow us to study features already known from 2D simulations as well as nonaxisymmetric effects. In agreement with recent results we find only type I gravitational wave signals at core bounce. In the later stage of the simulations, one of our models (s15g) shows nonaxisymmetric gravitational wave emission caused by a low T/|W| dynamical instability, while the other model radiates gravitational waves due to a convective instability in the protoneutron star. The total energy released in gravitational waves within the considered time intervals is 1.52\times10^{-7}M_{\odot} (s15g) and 4.72\times10^{-10}M_{\odot} (s15h). Both core collapse simulations indicate that corresponding events in our Galaxy would be detectable either by the LIGO or Advanced LIGO detector.

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.

Neutrino transport in spherically symmetric models of stellar core collapse and bounce has achieved a technically complete level, rewarded by the agreement among independent groups that a multi-dimensional treatment of the fluid-instabilities in the post-bounce phase is indispensable to model supernova explosions. While much effort is required to develop a reliable neutrino transport technique in axisymmetry, we explore neutrino physics approximations and parameterizations for an efficient three-dimensional simulation of the fluid-instabilities in the shock-heated matter that accumulates between the accretion shock and the protoneutron star. We demonstrate the reliability of a simple parameterization scheme in the collapse phase and extend our 3D magneto-hydrodynamical collapse simulations to a preliminary postbounce evolution. The growth of magnetic fields is investigated.

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.

A simple and computationally efficient parameterization of the deleptonization, the entropy changes, and the neutrino stress is presented for numerical simulations of stellar core collapse. The parameterization of the neutrino physics is based on a bounce profile of the electron fraction as it results from state-of-the-art collapse simulations with multi-group Boltzmann neutrino transport in spherical symmetry. Two additional parameters include an average neutrino escape energy and a neutrino trapping density. The parameterized simulations reproduce the consequences of the delicate neutrino thermalization/diffusion process during the collapse phase and provide a by far more realistic alternative to the adiabatic approximation, which has often been used in the investigation of the emission of gravitational waves, of multi-dimensional general relativistic effects, of the evolution of magnetic fields, or even of the nucleosynthesis in simulations of core collapse and bounce. For supernova codes that are specifically designed for the postbounce phase, the parameterization builds a convenient bridge between the point where the applicability of a stellar evolution code ends and the point where the postbounce evolution begins.

The role of ion-ion correlations in suppressing neutrino-nucleus elastic scattering during stellar core collapse is reinvestigated, using two different equations of state. We test the improved description by Itoh et al. against the treatment suggested by Horowitz and find that the stronger cross section reduction for small momentum transfer in the former case does not lead to noticeable changes of the core deleptonization and entropy increase during collapse, because the improvements are relevant below neutrino trapping conditions only for very low neutrino energies, corresponding to a very small phase space volume. Treating screening effects for ionic mixtures by the linear mixing rule, applied to the collection of representative heavy nucleus, alpha particles, and free nucleons, which is assumed to characterize the composition in nuclear statistical equilibrium, we cannot determine mentionable differences during stellar collapse, because alpha particles are not sufficiently abundant and their coherent scattering opacity is too small.

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.

Accurate neutrino transport has been built into spherically symmetric simulations of stellar core collapse and postbounce evolution. The results of such simulations agree that spherically symmetric models with standard microphysical input fail to explode by the delayed, neutrino-driven mechanism. Independent groups implemented fundamentally different numerical methods to tackle the Boltzmann neutrino transport equation. Here we present a direct and detailed comparison of such neutrino radiation-hydrodynamical simulations for two codes, Agile-Boltztran of the Oak Ridge-Basel group and Vertex of the Garching group. The former solves the Boltzmann equation directly by an implicit, general relativistic discrete angle method on the adaptive grid of a conservative implicit hydrodynamics code with second-order TVD advection. In contrast, the latter couples a variable Eddington factor technique with an explicit, moving-grid, conservative high-order Riemann solver with important relativistic effects treated by an effective gravitational potential. The presented study is meant to test both neutrino radiation-hydrodynamics implementations and to provide a data basis for comparisons and verifications of supernova codes to be developed in the future. Results are discussed for simulations of the core collapse and post-bounce evolution of a 13 solar mass star with Newtonian gravity and a 15 solar mass star with relativistic gravity.

All recent numerical simulations agree that stars in the main sequence mass range of 9-40 solar masses do not produce a prompt hydrodynamic ejection of the outer layers after core collapse and bounce. Rather they suggest that stellar core collapse and supernova explosion are dynamically distinct astrophysical events, separated by an unspectacular accretion phase of at least ~40 ms duration. As long as the neutrinospheres remain convectively stable, the explosion dynamics is determined by the neutrons, protons, electrons and neutrinos in the layer of impact-heated matter piling up on the protoneutron star. The crucial role of neutrino transport in this regime has been emphasized in many previous investigations. Here, we search for efficient means to address the role of magnetic fields and fluid instabilities in stellar core collapse and the postbounce phase.

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.

The number in the title is a birthday-inspired overestimate by about an order of magnitude. Nevertheless, this contribution reminds of well-known reasons that prevented straightforward explosions in many numerical supernova models. It emphasizes a separation between the two events of core collapse and supernova explosion.

The most important weak nuclear interaction to the dynamics of stellar core collapse is electron capture, primarily on nuclei with masses larger than 60. In prior simulations of core collapse, electron capture on these nuclei has been treated in a highly parameterized fashion, if not ignored. With realistic treatment of electron capture on heavy nuclei come significant changes in the hydrodynamics of core collapse and bounce. We discuss these as well as the ramifications for the post-bounce evolution in core collapse supernovae.

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.

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.

The remnant resulting from the merger of two neutron stars produces neutrinos in copious amounts. In this paper we present the neutrino emission results obtained via Newtonian, high-resolution simulations of the coalescence event. These simulations use three-dimensional smoothed particle hydrodynamics together with a nuclear, temperature dependent equation of state and a multi-flavour neutrino leakage scheme. We present the details of our scheme, discuss the neutrino emission results from a neutron star coalescence and compare them to the core-collapse supernova case where neutrino emission has been studied for several decades. The average neutrino energies are similar to those in the supernova case, but contrary to the latter, the luminosities are dominated by electron-type antineutrinos which are produced in the hot, neutron-rich, thick disk of the merger remnant. The cooler parts of this disk contain substantial fractions of heavy nuclei, which, however, do not influence the overall neutrino emission results significantly. Our total neutrino luminosities from the merger event are considerably lower than those found in previous investigations. This has serious consequences for the ability to produce a gamma-ray burst via neutrino annihilation. The neutrinos are emitted preferentially along the initial binary rotation axis, an event seen ``pole-on'' would appear much brighter in neutrinos than a similar event seen ``edge-on''.

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.

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 propose an approximation to general relativity that captures the main gravitational effects of dynamical importance in supernovae. The conceptual link between this formalism and the Newtonian limit is such that it could likely be implemented in existing multidimensional Newtonian gravitational hydrodynamics codes employing a Poisson solver. As a test of the formalism's utility, we display results for rapidly rotating (and therefore highly deformed) neutron stars.

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.

We present self-consistent general relativistic simulations of stellar core collapse, bounce, and postbounce evolution for 13, 15, and 20 solar mass progenitors in spherical symmetry. Our simulations implement three-flavor Boltzmann neutrino transport and standard nuclear physics. The results are compared to our corresponding simulations with Newtonian hydrodynamics and O(v/c) Boltzmann transport.

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.

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.

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.