The high temperature and electron degeneracy attained during a supernova allow for the formation of a large muon abundance within the core of the resulting proto-neutron star. If new pseudoscalar degrees of freedom have large couplings to the muon, they can be produced by this muon abundance and contribute to the cooling of the star. By generating the largest collection of supernova simulations with muons to date, we show that observations of the cooling rate of SN 1987A place strong constraints on the coupling of axion-like particles to muons, limiting the coupling to $g_{a\mu} < 10^{-8.1}~\text{GeV}^{-1}$.

We present 3D full-sphere supernova simulations of non-rotating low-mass (~9 Msun) progenitors, covering the entire evolution from core collapse through bounce and shock revival, through shock breakout from the stellar surface, until fallback is completed several days later. We obtain low-energy explosions [~(0.5-1.0)x 10^{50} erg] of iron-core progenitors at the low-mass end of the core-collapse supernova (LMCCSN) domain and compare to a super-AGB (sAGB) progenitor with an oxygen-neon-magnesium core that collapses and explodes as electron-capture supernova (ECSN). The onset of the explosion in the LMCCSN models is modelled self-consistently using the Vertex-Prometheus code, whereas the ECSN explosion is modelled using parametric neutrino transport in the Prometheus-HOTB code, choosing different explosion energies in the range of previous self-consistent models. The sAGB and LMCCSN progenitors that share structural similarities have almost spherical explosions with little metal mixing into the hydrogen envelope. A LMCCSN with less 2nd dredge-up results in a highly asymmetric explosion. It shows efficient mixing and dramatic shock deceleration in the extended hydrogen envelope. Both properties allow fast nickel plumes to catch up with the shock, leading to extreme shock deformation and aspherical shock breakout. Fallback masses of <~5x10^{-3} Msun have no significant effects on the neutron star (NS) masses and kicks. The anisotropic fallback carries considerable angular momentum, however, and determines the spin of the newly-born NS. The LMCCSNe model with less 2nd dredge-up results in a hydrodynamic and neutrino-induced NS kick of >40 km/s and a NS spin period of ~30 ms, both not largely different from those of the Crab pulsar at birth.

Neutrinos are unique probes of core-collapse supernova dynamics, especially in the case of black hole (BH) forming stellar collapses, where the electromagnetic emission may be faint or absent. By investigating two 3D hydrodynamical simulations of BH-forming stellar collapses of mass 40 and 75 M_sun, we identify the physical processes preceding BH formation through neutrinos, and forecast the neutrino signal expected in the existing IceCube and Super-Kamiokande detectors, as well as in the future generation DUNE facility. Prior to the abrupt termination of the neutrino signal corresponding to BH formation, both models develop episodes of strong and long-lasting activity by the spiral standing accretion shock instability (SASI). We find that the spiral SASI peak in the Fourier power spectrum of the neutrino event rate will be distinguishable at 3 sigma above the detector noise for distances up to O(30) kpc in the most optimistic scenario, with IceCube having the highest sensitivity. Interestingly, given the long duration of the spiral SASI episodes, the spectrograms of the expected neutrino event rate carry clear signs of the evolution of the blue spiral SASI frequency as a function of time, as the shock radius and post-shock fluid velocity evolve. Due to the high accretion luminosity and its large-amplitude SASI-induced modulations, any contribution from asymmetric (dipolar or quadrupolar) neutrino emission associated with the lepton emission self-sustained asymmetry (LESA) is far subdominant in the neutrino signal.

The possible detection of a compact object in the remnant of SN 1987A presents an unprecedented opportunity to follow its early evolution. The suspected detection stems from an excess of infrared emission from a dust blob near the compact object's predicted position. The infrared excess could be due to the decay of isotopes like 44Ti, accretion luminosity from a neutron star or black hole, magnetospheric emission or a wind originating from the spindown of a pulsar, or thermal emission from an embedded, cooling neutron star (NS 1987A). It is shown that the last possibility is the most plausible as the other explanations are disfavored by other observations and/or require fine-tuning of parameters. Not only are there indications the dust blob overlaps the predicted location of a kicked compact remnant, but its excess luminosity also matches the expected thermal power of a 30 year old neutron star. Furthermore, models of cooling neutron stars within the Minimal Cooling paradigm readily fit both NS 1987A and Cas A, the next-youngest known neutron star. If correct, a long heat transport timescale in the crust and a large effective stellar temperature are favored, implying relatively limited crustal n-1S0 superfluidity and an envelope with a thick layer of light elements, respectively. If the locations don't overlap, then pulsar spindown or accretion might be more likely, but the pulsar's period and magnetic field or the accretion rate must be rather finely tuned. In this case, NS 1987A may have enhanced cooling and/or a heavy-element envelope.

Comparison of theoretical line profiles to observations provides important tests for supernova explosion models. We study the shapes of radioactive decay lines predicted by current 3D core-collapse explosion simulations, and compare these to observations of SN 1987A and Cas A. Both the widths and shifts of decay lines vary by several thousand kilometers per second depending on viewing angle. The line profiles can be complex with multiple peaks. By combining observational constraints from 56Co decay lines, 44Ti decay lines, and Fe IR lines, we delineate a picture of the morphology of the explosive burning ashes in SN 1987A. For M_ZAMS=15-20 Msun progenitors exploding with ~1.5 *10^51 erg, ejecta structures suitable to reproduce the observations involve a bulk asymmetry of the 56Ni of at least ~400 km/s and a bulk velocity of at least ~1500 km/s. By adding constraints to reproduce the UVOIR bolometric light curve of SN 1987A up to 600d, an ejecta mass around 14 Msun is favoured. We also investigate whether observed decay lines can constrain the neutron star (NS) kick velocity. The model grid provides a constraint V_NS > V_redshift, and applying this to SN 1987A gives a NS kick of at least 500 km/s. For Cas A, our single model provides a satisfactory fit to the NuSTAR observations and reinforces the result that current neutrino-driven core-collapse SN models can achieve enough bulk asymmetry in the explosive burning material. Finally, we investigate the internal gamma-ray field and energy deposition, and compare the 3D models to 1D approximations.

The mass function for black holes and neutron stars at birth is explored for mass-losing helium stars. These should resemble, more closely than similar studies of single hydrogen-rich stars, the results of evolution in close binary systems. The effects of varying the mass-loss rate and metallicity are calculated using a simple semi-analytic approach to stellar evolution that is tuned to reproduce detailed numerical calculations. Though the total fraction of black holes made in stellar collapse events varies considerably with metallicity, mass-loss rate, and mass cutoff, from 5$\%$ to 30$\%$, the shapes of their birth functions are very similar for all reasonable variations in these quantities. Median neutron star masses are in the range 1.32 - 1.37 $M_\odot$ regardless of metallicity. The median black hole mass for solar metallicity is typically 8 to 9 $M_\odot$ if only initial helium cores below 40 $M_\odot$ (ZAMS mass less than 80 $M_\odot$) are counted, and 9 - 13 $M_\odot$, in most cases, if helium cores with initial masses up to 150 $M_\odot$ (ZAMS mass less than 300 $M_\odot$) contribute. As long as the mass-loss rate as a function of mass exhibits no strong non-linearities, the black hole birth function from 15 to 35 $M_\odot$ has a slope that depends mostly on the initial mass function for main sequence stars. These findings imply the possibility of constraining the initial mass function and the properties of mass loss in close binaries using ongoing measurements of gravitational wave radiation. The expected rotation rates of the black holes are briefly discussed.

The release of spin-down energy by a magnetar is a promising scenario to power several classes of extreme explosive transients. However, it lacks a firm basis because magnetar formation still represents a theoretical challenge. Using the first three-dimensional simulations of a convective dynamo based on a protoneutron star interior model, we demonstrate that the required dipolar magnetic field can be consistently generated for sufficiently fast rotation rates. The dynamo instability saturates in the magnetostrophic regime with the magnetic energy exceeding the kinetic energy by a factor of up to 10. Our results are compatible with the observational constraints on galactic magnetar field strength and provide strong theoretical support for millisecond protomagnetar models of gamma-ray burst and superluminous supernova central engines.

We present a carefully designed, systematic study of the angular resolution dependence of simulations with the Prometheus-Vertex neutrino-hydrodynamics code. Employing a simplified neutrino heating-cooling scheme in the Prometheus hydrodynamics module allows us to sample the angular resolution between 4 degrees and 0.5 degrees. With a newly-implemented static mesh refinement (SMR) technique on the Yin-Yang grid, the angular coordinates can be refined in concentric shells, compensating for the diverging structure of the spherical grid. In contrast to previous studies with Prometheus and other codes, we find that higher angular resolution and therefore lower numerical viscosity provides more favorable explosion conditions and faster shock expansion. We discuss the possible reasons for the discrepant results. The overall dynamics seem to converge at a resolution of about 1 degree. Applying the SMR setup to marginally exploding progenitors is disadvantageous for the shock expansion, however, because kinetic energy of downflows is dissipated to internal energy at resolution interfaces, leading to a loss of turbulent pressure support and a steeper temperature gradient. We also present a way to estimate the numerical viscosity on grounds of the measured turbulent kinetic-energy spectrum, leading to smaller values that are better compatible with the flow behavior witnessed in our simulations than results following calculations in previous literature. Interestingly, the numerical Reynolds numbers in the turbulent, neutrino-heated postshock layer (some 10 to several 100) are in the ballpark of expected neutrino-drag effects on the relevant length scales in the turbulent postshock layer. We provide a formal derivation and quantitative assessment of the neutrino drag terms in an appendix.

Neutrinos from a supernova (SN) might undergo fast flavor conversions near the collapsed stellar core. We perform a detailed study of this intriguing possibility, analyzing time-dependent state-of-the-art 3D SN models of 9 and 20 Msun. Both models were computed with multi-D three-flavor neutrino transport based on a two-moment solver, and both exhibit the presence of the lepton-number emission self-sustained asymmetry (LESA). The transport solution does not provide the angular distributions of the neutrino fluxes, which are crucial to track the fast flavor instability. To overcome this limitation, we use a recently proposed approach based on the angular moments of the energy-integrated electron lepton-number distribution. With this method we find the possibility of fast neutrino flavor instability at radii <~20 km, which is well interior to the neutrinosphere. Our results confirm recent observations in a 2D SN model and in 2D/3D models with fixed matter background, which were computed with Boltzmann neutrino transport. However, the flavor unstable locations are not isolated points as discussed previously, but thin skins surrounding volumes where electron antineutrinos are more abundant than electron neutrinos. These volumes grow with time and appear first in the convective layer of the proto-neutron star (PNS), where a decreasing electron fraction (Ye) and high temperatures favor the occurrence of regions with negative neutrino chemical potential. Since Ye remains higher in the LESA dipole direction, where convective lepton-number transport out from the nonconvective PNS core slows down the deleptonization, flavor unstable conditions become more widespread in the opposite hemisphere. This interesting phenomenon deserves further investigation, since its impact on SN modeling and possible consequences for SN dynamics and neutrino observations are presently unclear. (abridged)

Light curves, explosion energies, and remnant masses are calculated for a grid of supernovae resulting from massive helium stars that have been evolved including mass loss. These presupernova stars should approximate the results of binary evolution for stars in interacting systems that lose their envelopes close to the time of helium core ignition. Initial helium star masses are in the range 2.5 to 40\,\Msun, which correspond to main sequence masses of about 13 to 90\,\Msun. Common Type Ib and Ic supernovae result from stars whose final masses are approximately 2.5 to 5.6\,\Msun. For heavier stars, a large fraction of collapses lead to black holes, though there is an island of explodability for presupernova masses near 10\,\Msun. The median neutron star mass in binaries is 1.35--1.38\,\Msun \ and the median black hole mass is between 9 and 11\,\Msun. Even though black holes less massive than 5 \Msun\ are rare, they are predicted down to the maximum neutron star mass. There is no empty ``gap'', only a less populated mass range. For standard assumptions regarding the explosions and nucleosynthesis, the models predict light curves that are fainter than the brighter common Type Ib and Ic supernovae. Even with a very liberal, but physically plausible increase in $^{56}$Ni production, the highest energy models are fainter, at peak, than 10$^{42.6}$\,erg\,s$^{-1}$, and very few approach that limit. The median peak luminosity ranges from 10$^{42.0}$ to 10$^{42.3}$\,erg\,s$^{-1}$. Possible alternatives to the standard neutrino-powered and radioactive-illuminated models are explored. Magnetars are a promising alternative. Several other unusual varieties of Type I supernovae at both high and low mass are explored.

We present a seven-minute long $4\pi$-3D simulation of a shell merger event in a non-rotating $18.88\, M_\odot$ supernova progenitor before the onset of gravitational collapse. The key motivation is to capture the large-scale mixing and asymmetries in the wake of the shell merger before collapse using a self-consistent approach. The $4\pi$ geometry is crucial as it allows us to follow the growth and evolution of convective modes on the largest possible scales. We find significant differences between the kinematic, thermodynamic and chemical evolution of the 3D and the 1D model. The 3D model shows vigorous convection leading to more efficient mixing of nuclear species. In the 3D case the entire oxygen shell attains convective Mach numbers of $\mathord{\approx}\, 0.1$, whereas in the 1D model, the convective velocities are much lower and there is negligible overshooting across convective boundaries. In the 3D case, the convective eddies entrain nuclear species from the neon (and carbon) layers into the deeper part of the oxygen burning shell, where they burn and power a violent convection phase with outflows. This is a prototypical model of a convective-reactive system. Due to the strong convection and the resulting efficient mixing, the interface between the neon layer and the silicon-enriched oxygen layer disappears during the evolution, and silicon is mixed far out into merged oxygen/neon shell. Neon entrained inwards by convective downdrafts burns, resulting in lower neon mass in the 3D model compared to the 1D model at time of collapse. In addition, the 3D model develops remarkable large-scale, large-amplitude asymmetries, which may have important implications for the impending gravitational collapse and the subsequent explosion.

The compact binary radio pulsar system J0453+1559 (Martinez et al. 2015) consists of a recycled pulsar as primary component of $1.559(5)\,M_\odot$ and an unseen companion star of $1.174(4)\,M_\odot$. Because of the relatively large orbital eccentricity of $e = 0.1125$, it was argued that the companion is a neutron star, making it the neutron star with the lowest accurately determined mass to date. However, a direct observational determination of the nature of the companion is currently not feasible. Moreover, state-of-the-art stellar evolution and supernova modeling are contradictive concerning the possibility to produce such a low-mass neutron star remnant. Here we challenge the neutron star interpretation by reasoning that the lower-mass component could instead be a white dwarf born in a thermonuclear electron-capture supernova (tECSN) event, in which oxygen-neon deflagration in the degenerate stellar core of an ultra-stripped progenitor ejects several $0.1\,M_\odot$ of matter and leaves a bound ONeFe white dwarf as the second-formed compact remnant. We determine the ejecta mass and remnant kick needed in this scenario to explain the properties of PSR~J0453+1559 by a neutron star--white dwarf system. More work on tECSNe is needed to assess the viability of this scenario.

We present high angular resolution (~80 mas) ALMA continuum images of the SN 1987A system, together with CO $J$=2 $\!\rightarrow\!$ 1, $J$=6 $\!\rightarrow\!$ 5, and SiO $J$=5 $\!\rightarrow\!$ 4 to $J$=7 $\!\rightarrow\!$ 6 images, which clearly resolve the ejecta (dust continuum and molecules) and ring (synchrotron continuum) components. Dust in the ejecta is asymmetric and clumpy, and overall the dust fills the spatial void seen in H$\alpha$ images, filling that region with material from heavier elements. The dust clumps generally fill the space where CO $J$=6 $\!\rightarrow\!$ 5 is fainter, tentatively indicating that these dust clumps and CO are locationally and chemically linked. In these regions, carbonaceous dust grains might have formed after dissociation of CO. The dust grains would have cooled by radiation, and subsequent collisions of grains with gas would also cool the gas, suppressing the CO $J$=6 $\!\rightarrow\!$ 5 intensity. The data show a dust peak spatially coincident with the molecular hole seen in previous ALMA CO $J$=2 $\!\rightarrow\!$ 1 and SiO $J$=5 $\!\rightarrow\!$ 4 images. That dust peak, combined with CO and SiO line spectra, suggests that the dust and gas could be at higher temperatures than the surrounding material, though higher density cannot be totally excluded. One of the possibilities is that a compact source provides additional heat at that location. Fits to the far-infrared--millimeter spectral energy distribution give ejecta dust temperatures of 18--23K. We revise the ejecta dust mass to $\mathrm{M_{dust}} = 0.2-0.4$M$_\odot$ for carbon or silicate grains, or a maximum of $<0.7$M$_\odot$ for a mixture of grain species, using the predicted nucleosynthesis yields as an upper limit.

We present the new code NADA-FLD to solve multi-dimensional neutrino-hydrodynamics in full general relativity (GR) in spherical polar coordinates. The energy-dependent neutrino transport assumes the flux-limited diffusion (FLD) approximation and evolves the neutrino energy densities measured in the frame comoving with the fluid. Operator splitting is used to avoid multi-dimensional coupling of grid cells in implicit integration steps involving matrix inversions. Terms describing lateral diffusion and advection are integrated explicitly using the Allen-Cheng or the Runge-Kutta-Legendre method, which remain stable even in the optically thin regime. We discuss several toy-model problems in one and two dimensions to test the basic functionality and individual components of the transport scheme. We also perform fully dynamic core-collapse supernova (CCSN) simulations in spherical symmetry. For a Newtonian model we find good agreement with the M1 code ALCAR, and for a GR model we reproduce the main effects of GR in CCSNe already found by previous works.

Rotation of core-collapse supernovae (SNe) affects the neutrino emission characteristics. By comparing the neutrino properties of three three-dimensional SN simulations of a 15 M_sun progenitor (one non-rotating model and two models rotating at different velocities), we investigate how the neutrino emission varies with the flow dynamics in the SN core depending on the different degrees of rotation. The large-amplitude sinusoidal modulations due to the standing accretion-shock instability (SASI) are weaker in both the rotating models than in the non-rotating case. The SN progenitor rotation reduces the radial velocities and radial component of the kinetic energy associated with convection interior to the proto-neutron star. This effect seems to disfavor the growth of the hemispheric neutrino-emission asymmetries associated with the lepton-emission self-sustained asymmetry (LESA). An investigation of the multipole expansion of the neutrino luminosity and the electron neutrino lepton number flux shows a dominant quadrupolar mode in rotating SN models. Our findings highlight the power of using neutrinos as probes of SN hydrodynamics.

During the first few hundred days after the explosion, core-collapse supernovae (SNe) emit down-scattered X-rays and gamma-rays originating from radioactive line emissions, primarily from the $^{56}$Ni $\rightarrow$ $^{56}$Co $\rightarrow$ $^{56}$Fe chain. We use SN models based on three-dimensional neutrino-driven explosion simulations of single stars and mergers to compute this emission and compare the predictions with observations of SN 1987A. A number of models are clearly excluded, showing that high-energy emission is a powerful way of discriminating between models. The best models are almost consistent with the observations, but differences that cannot be matched by a suitable choice of viewing angle are evident. Therefore, our self-consistent models suggest that neutrino-driven explosions are able to produce, in principle, sufficient mixing, although remaining discrepancies may require small changes to the progenitor structures. The soft X-ray cutoff is primarily determined by the metallicity of the progenitor envelope. The main effect of asymmetries is to vary the flux level by a factor of ${\sim}$3. For the more asymmetric models, the shapes of the light curves also change. In addition to the models of SN 1987A, we investigate two models of Type II-P SNe and one model of a stripped-envelope Type IIb SN. The Type II-P models have similar observables as the models of SN 1987A, but the stripped-envelope SN model is significantly more luminous and evolves faster. Finally, we make simple predictions for future observations of nearby SNe.

The neutrino process that occurs in the outer stellar shells during a supernova explosion and involves neutrino-nucleus reactions produces a range of rare, stable and radioactive isotopes. We improve previous $\nu$-process studies by using, for the first time, the time-dependent neutrino emission spectra, as predicted from supernova simulations, rather than a simplified parametric description modeled after the neutron-star cooling phase. In particular, our calculations use time-dependent neutrino spectra for all neutrino species, consider their deviation from a Fermi-Dirac distribution and account for the neutrino emission from the neutrino burst and accretion phases. We find that the time-dependent treatment of the neutrino emission spectra results in higher yields for the selected nuclei produced by the $\nu$~process as compared to previous studies and also compared to the approximation of assuming constant neutrino energies corresponding to the time-averaged mean energy radiated in each species. The effect is largest for nuclides produced by charged-current reactions. Our results reflect the dynamical competition between neutrino-induced reactions and the effect of the shock passage through the star. By varying the neutrino burst luminosity and the duration of the accretion phase, we study the impact of these early emission phases and their uncertainties on the $\nu$-process nucleosynthesis. We find that the deviation of the neutrino spectra from a Fermi-Dirac distribution calculated in supernova simulations has a negligible effect on the $\nu$-process yields.

A set of eight self-consistent, time-dependent supernova (SN) simulations in three spatial dimensions (3D) for 9 solar-mass and 20 solar-mass progenitors is evaluated for the presence of dipolar asymmetries of the electron lepton-number emission as discovered by Tamborra et al. and termed lepton-number emission self-sustained asymmetry (LESA). The simulations were performed with the Aenus-Alcar neutrino/hydrodynamics code, which treats the energy- and velocity-dependent transport of neutrinos of all flavors by a two-moment scheme with algebraic M1 closure. For each of the progenitors, results with fully multi-dimensional (FMD) neutrino transport and with ray-by-ray-plus (RbR+) approximation are considered for two different grid resolutions. While the 9 solar-mass models develop explosions, the 20 solar-mass progenitor does not explode with the employed version of simplified neutrino opacities. In all 3D models we observe the growth of substantial dipole amplitudes of the lepton-number (electron neutrino minus antineutrino) flux with stable or slowly time-evolving direction and overall properties fully consistent with the LESA phenomenon. Models with RbR+ transport develop LESA dipoles somewhat faster and with temporarily higher amplitudes, but the FMD calculations exhibit cleaner hemispheric asymmetries with a far more dominant dipole. In contrast, the RbR+ results display much wider multipole spectra of the neutrino-emission anisotropies with significant power also in the quadrupole and higher-order modes. Our results disprove speculations that LESA is a numerical artifact of RbR+ transport. We also discuss LESA as consequence of a dipolar convection flow inside of the nascent neutron star and establish, tentatively, a connection to Chandrasekhar's linear theory of thermal instability in spherical shells.

We present predictions for the gravitational-wave (GW) emission of three-dimensional supernova (SN) simulations performed for a 15 solar-mass progenitor with the Prometheus-Vertex code using energy-dependent, three-flavor neutrino transport. The progenitor adopted from stellar evolution calculations including magnetic fields had a fairly low specific angular momentum (j_Fe <~ 10^{15} cm^2/s) in the iron core (central angular velocity ~0.2 rad/s), which we compared to simulations without rotation and with artificially enhanced rotation (j_Fe <~ 2*10^{16} cm^2/s; central angular velocity ~0.5 rad/s). Our results confirm that the time-domain GW signals of SNe are stochastic, but possess deterministic components with characteristic patterns at low frequencies (<~200 Hz), caused by mass motions due to the standing accretion shock instability (SASI), and at high frequencies, associated with gravity-mode oscillations in the surface layer of the proto-neutron star (PNS). Non-radial mass motions in the post-shock layer as well as PNS convection are important triggers of GW emission, whose amplitude scales with the power of the hydrodynamic flows. There is no monotonic increase of the GW amplitude with rotation, but a clear correlation with the strength of SASI activity. Our slowly rotating model is a fainter GW emitter than the non-rotating model because of weaker SASI activity and damped convection in the post-shock layer and PNS. In contrast, the faster rotating model exhibits a powerful SASI spiral mode during its transition to explosion, producing the highest GW amplitudes with a distinctive drift of the low-frequency emission peak from ~80-100 Hz to ~40-50 Hz. This migration signifies shock expansion, whereas non-exploding models are discriminated by the opposite trend.

We describe a multi-messenger interpretation of GW170817, which yields a robust lower limit on NS radii. This excludes NSs with radii smaller than about 10.7 km and thus rules out very soft nuclear matter. We stress the potential of this type of constraints when future detections become available. A very similar argumentation may yield an upper bound on the maximum mass of nonrotating NSs. We also discuss simulations of NS mergers, which undergo a first-order phase transition to quark matter. We point out a different dynamical behavior. Considering the gravitational-wave signal, we identify an unambiguous signature of the QCD phase transition in NS mergers. The occurrence of quark matter through a strong first-order phase transition during merging leads to a characteristic shift of the dominant postmerger frequency. The frequency shift is indicative for a phase transition if it is compared to the postmerger frequency which is expected for purely hadronic EoS models. A very strong deviation of several 100 Hz is observed for hybrid EoSs in an otherwise tight relation between the tidal deformability and the postmerger frequency. We address the potential impact of a first-order phase transition on the electromagnetic counterpart of NS mergers. Our simulations suggest that there would be no significant qualitative differences between a system undergoing a phase transition to quark matter and purely hadronic mergers. The quantitative differences are within the spread which is found between different hadronic EoS models. This implies on the one hand that GW170817 is compatible with a possible transition to quark matter. On the other hand these considerations show that it may not be easy to identify quantitative differences between purely hadronic mergers and events in which quark matter occurs considering solely their electromagnetic counterpart or their nucleosynthesis products. (abridged)

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.

With the same method as used previously, we investigate neutrino-driven explosions of a larger sample of blue supergiant models. The larger sample includes three new presupernova stars. The results are compared with light-curve observations of the peculiar type IIP SN 1987A. The explosions were modeled in 3D with the neutrino-hydrodynamics code PROMETHEUS-HOTB, and light-curve calculations were performed in spherical symmetry with the radiation-hydrodynamics code CRAB. Our results confirm the basic findings of the previous work: 3D neutrino-driven explosions with SN 1987A-like energies synthesize an amount of Ni-56 that is consistent with the radioactive tail of the light curve. Moreover, the models mix hydrogen inward to minimum velocities below 400 km/s as required by spectral observations. Hydrodynamic simulations with the new progenitor models, which possess smaller radii than the older ones, show much better agreement between calculated and observed light curves in the initial luminosity peak and during the first 20 days. A set of explosions with similar energies demonstrated that a high growth factor of Rayleigh-Taylor instabilities at the (C+O)/He composition interface combined with a weak interaction of fast Rayleigh-Taylor plumes, where the reverse shock occurs below the He/H interface, provides a sufficient condition for efficient outward mixing of Ni-56 into the hydrogen envelope. This condition is realized to the required extent only in one of the older stellar models, which yielded a maximum velocity of around 3000 km/s for the bulk of ejected Ni-56, but failed to reproduce the helium-core mass of 6 Msun inferred from the absolute luminosity of the presupernova star. We conclude that none of the single-star progenitor models proposed for SN 1987A to date satisfies all constraints set by observations. (Abridged)

We present a new, computationally efficient, energy-integrated approximation for neutrino effects in hot and dense astrophysical environments such as supernova cores and compact binary mergers and their remnants. Our new method, termed ILEAS for Improved Leakage-Equilibration-Absorption Scheme, improves the lepton-number and energy losses of traditional leakage descriptions by a novel prescription of the diffusion time-scale based on a detailed energy integral of the flux-limited diffusion equation. The leakage module is supplemented by a neutrino-equilibration treatment that ensures the proper evolution of the total lepton number and medium plus neutrino energies as well as neutrino-pressure effects in the neutrino-trapping domain. Moreover, we employ a simple and straightforwardly applicable ray-tracing algorithm for including re-absorption of escaping neutrinos especially in the decoupling layer and during the transition to semi-transparent conditions. ILEAS is implemented on a three-dimensional (3D) Cartesian grid with a minimum of free and potentially case-dependent parameters and exploits the basic physics constraints that should be fulfilled in the neutrino-opaque and free-streaming limits. We discuss a suite of tests for stationary and time-dependent proto-neutron star models and post-merger black-hole-torus configurations, for which 3D ILEAS results are demonstrated to agree with energy-dependent 1D and 2D two-moment (M1) neutrino transport on the level of 10--15 percent in basic neutrino properties. This also holds for the radial profiles of the neutrino luminosities and of the electron fraction. Even neutrino absorption maps around torus-like neutrino sources are qualitatively similar without any fine-tuning, confirming that ILEAS can satisfactorily reproduce local losses and re-absorption of neutrinos as found in sophisticated transport calculations.

Self-consistent, time-dependent supernova (SN) simulations in three spatial dimensions (3D) are conducted with the Aenus-Alcar code, comparing, for the first time, calculations with fully multi-dimensional (FMD) neutrino transport and the ray-by-ray-plus (RbR+) approximation, both based on a two-moment solver with algebraic M1 closure. We find good agreement between 3D results with FMD and RbR+ transport for both tested grid resolutions in the cases of a 20 solar-mass progenitor, which does not explode with the employed simplified set of neutrino opacities, and of an exploding 9 solar-mass model. This is in stark contrast to corresponding axisymmetric (2D) simulations, which confirm previous claims that the RbR+ approximation can foster explosions in 2D in particular in models with powerful axial sloshing of the stalled shock due to the standing accretion shock instability (SASI). However, while local and instantaneous variations of neutrino fluxes and heating rates can still be considerably higher with RbR+ transport in 3D, the time-averaged quantities are very similar to FMD results because of the absence of a fixed, artificial symmetry axis that channels the flow. Therefore, except for stochastic fluctuations, the neutrino signals and the post-bounce evolution of 3D simulations with FMD and RbR+ transport are also very similar, in particular for our calculations with the better grid resolution. Higher spatial resolution has clearly a more important impact than the differences by the two transport treatments. Our results back up the use of the RbR+ approximation for neutrino transport in 3D SN modeling.

Measuring the rotation of core-collapse supernovae (SN) and of their progenitor stars is extremely challenging. Here it is demonstrated that neutrinos may potentially be employed as stellar gyroscopes, if phases of activity by the standing accretion-shock instability (SASI) affect the neutrino emission prior to the onset of the SN explosion. This is shown by comparing the neutrino emission properties of self-consistent, three-dimensional (3D) SN simulations of a 15 M_sun progenitor without rotation as well as slow and fast rotation compatible with observational constraints. The explosion of the fast rotating model gives rise to long-lasting, massive polar accretion downflows with stochastic time-variability, detectable e.g. by the IceCube Neutrino Observatory for any observer direction. While spectrograms of the neutrino event rate of non-rotating SNe feature a well-known sharp peak due to SASI for observers located in the proximity of the SASI plane, the corresponding spectrograms of rotating models show activity over a wide range of frequencies, most notably above 200 Hz for rapid rotation. In addition, the Fourier power spectra of the event rate for rotating models exhibit a SASI peak with lower power than in non-rotating models. The spectra for the rotating models also show secondary peaks at higher frequencies with greater relative heights compared to the main SASI peak than for non-rotating cases. These rotational imprints will be detectable for SNe at 10 kpc or closer.

We present spherically symmetric (1D) and axisymmetric (2D) supernova simulations for a convection-dominated 9 Msun and a 20 Msun progenitor that develops violent activity by the standing-accretion-shock instability (SASI). We compare in detail the Aenus-Alcar code, which uses fully multidimensional two-moment neutrino transport with an M1 closure, with a ray-by-ray-plus (RbR+) version of this code and with the Prometheus-Vertex code that employs RbR+ two-moment transport with a Boltzmann closure. Besides testing consequences of ignored non-radial neutrino-flux components in the RbR+ approximation, we also discuss the influence of various transport ingredients applied or not applied in recent literature, namely simplified neutrino-pair processes, neutrino-electron scattering, velocity-dependent and gravitational-redshift terms, and strangeness and many-body corrections for neutrino-nucleon scattering. Alcar and Vertex show excellent agreement in 1D and 2D despite a slightly but systematically smaller radius (~1km) and stronger convection of the proto-neutron star with Alcar. As found previously, the RbR+ approximation is conducive to explosions, but much less severely in the convection-dominated 9 Msun case than in the marginally exploding 20 Msun model, where the onset time of explosion also exhibits big stochastic variations, and the RbR+ approximation has no distinctly stronger supportive effect than simplified pair processes or strangeness and many-body corrections. Neglecting neutrino-electron scattering has clearly unfavorable effects for explosions, while ignoring velocity and gravitational-redshift effects can both promote or delay the explosion. The ratio of advection timescale to neutrino-heating timescale in 1D simulations is a sensitive indicator of the influence of physics ingredients on explosions also in multidimensional simulations.

Neutron stars (NSs) obtain kicks of typically several 100 km/s at birth. The gravitational tug-boat mechanism can explain these kicks as consequences of asymmetric mass ejection during the supernova (SN) explosion. Support for this hydrodynamic explanation is provided by observations of SN remnants with associated NSs, which confirm the prediction that the bulk of the explosion ejecta, in particular chemical elements between silicon and the iron group, are dominantly expelled in the hemisphere opposite to the direction of the NS kick. Here, we present a large set of two- and three-dimensional explosion simulations of electron-capture SNe, considering explosion energies between ~3x10^49 erg and ~1.6x10^50 erg. We find that the fast acceleration of the SN shock in the steep density gradient delimiting the O-Ne-Mg core of the progenitor enables such a rapid expansion of neutrino-heated matter that the growth of neutrino-driven convection freezes out quickly in a high-mode spherical harmonics pattern. Since the corresponding momentum asymmetry of the ejecta is very small and the gravitational acceleration by the fast-expanding ejecta abates rapidly, the NS kick velocities are at most a few km/s. The extremely low core compactness of O-Ne-Mg-core progenitors therefore favors hydrodynamic NS kicks much below the ~160 km/s measured for the Crab pulsar. This suggests either that the Crab Nebula is not the remnant of an electron-capture SN, but of a low-mass iron-core progenitor, or that the Crab pulsar was not accelerated by the gravitational tug-boat mechanism but received its kick by a non-hydrodynamic mechanism such as, e.g., anisotropic neutrino emission.

Despite more than 30 years of searches, the compact object in Supernova (SN) 1987A has not yet been detected. We present new limits on the compact object in SN 1987A using millimeter, near-infrared, optical, ultraviolet, and X-ray observations from ALMA, VLT, HST, and Chandra. The limits are approximately 0.1 mJy ($0.1\times 10^{-26}$ erg s$^{-1}$ cm$^{-2}$ Hz$^{-1}$) at 213 GHz, 1 Lsun ($6\times 10^{-29}$ erg s$^{-1}$ cm$^{-2}$ Hz$^{-1}$) in optical if our line-of-sight is free of ejecta dust, and $10^{36}$ erg s$^{-1}$ ($2\times 10^{-30}$ erg s$^{-1}$ cm$^{-2}$ Hz$^{-1}$) in 2-10 keV X-rays. Our X-ray limits are an order of magnitude less constraining than previous limits because we use a more realistic ejecta absorption model based on three-dimensional neutrino-driven SN explosion models (presented in an accompanying article). The allowed bolometric luminosity of the compact object is 22 Lsun if our line-of-sight is free of ejecta dust, or 138 Lsun if dust-obscured. Depending on assumptions, these values limit the effective temperature of a neutron star to <4-8 MK and do not exclude models, which typically are in the range 3-4 MK. For the simplest accretion model, the accretion rate for an efficiency $\eta$ is limited to $< 10^{-11} \eta^{-1}$ Msun yr$^{-1}$, which excludes most predictions. For pulsar activity modeled by a rotating magnetic dipole in vacuum, the limit on the magnetic field strength ($B$) for a given spin period ($P$) is $B < 10^{14} P^2$ G s$^{-2}$. By combining information about radiation reprocessing and geometry, it is likely that the compact object is a dust-obscured thermally-emitting neutron star, which may appear as a region of higher-temperature ejecta dust emission.

We present a new method for solving the three-dimensional gravitational potential of a density field on the Yin-Yang grid. Our algorithm is based on a multipole decomposition and completely symmetric with respect to the two Yin-Yang grid patches. It is particularly efficient on distributed-memory machines with a large number of compute tasks, because the amount of data being explicitly communicated is minimized. All operations are performed on the original grid without the need for interpolating data onto an auxiliary spherical mesh.

The material expelled by core-collapse supernova (SN) explosions absorbs X-rays from the central regions. We use SN models based on three-dimensional neutrino-driven explosions to estimate optical depths to the center of the explosion, compare different progenitor models, and investigate the effects of explosion asymmetries. The optical depths below 2 keV for progenitors with a remaining hydrogen envelope are expected to be high during the first century after the explosion due to photoabsorption. A typical optical depth is $100 t_4^{-2} E^{-2}$, where $t_4$ is the time since the explosion in units of 10 000 days (${\sim}$27 years) and $E$ the energy in units of keV. Compton scattering dominates above 50 keV, but the scattering depth is lower and reaches unity already at ${\sim}$1000 days at 1 MeV. The optical depths are approximately an order of magnitude lower for hydrogen-stripped progenitors. The metallicity of the SN ejecta is much higher than in the interstellar medium, which enhances photoabsorption and makes absorption edges stronger. These results are applicable to young SN remnants in general, but we explore the effects on observations of SN 1987A and the compact object in Cas A in detail. For SN 1987A, the absorption is high and the X-ray upper limits of ${\sim}$100 Lsun on a compact object are approximately an order of magnitude less constraining than previous estimates using other absorption models. The details are presented in an accompanying paper. For the central compact object in Cas A, we find no significant effects of our more detailed absorption model on the inferred surface temperature.

We present a comparison between several simulation codes designed to study the core-collapse supernova mechanism. We pay close attention to controlling the initial conditions and input physics in order to ensure a meaningful and informative comparison. Our goal is three-fold. First, we aim to demonstrate the current level of agreement between various groups studying the core-collapse supernova central engine. Second, we desire to form a strong basis for future simulation codes and methods to compare to. Lastly, we want this work to be a stepping stone for future work exploring more complex simulations of core-collapse supernovae, i.e., simulations in multiple dimensions and simulations with modern neutrino and nuclear physics. We compare the early (first ~500ms after core bounce) spherically-symmetric evolution of a 20 solar mass progenitor star from six different core-collapse supernovae codes: 3DnSNe-IDSA, AGILE-BOLTZTRAN, FLASH, F{\sc{ornax}}, GR1D, and PROMETHEUS-VERTEX. Given the diversity of neutrino transport and hydrodynamic methods employed, we find excellent agreement in many critical quantities, including the shock radius evolution and the amount of neutrino heating. Our results provide an excellent starting point from which to extend this comparison to higher dimensions and compare the development of hydrodynamic instabilities that are crucial to the supernova explosion mechanism, such as turbulence and convection.

The birth properties of neutron stars yield important information on the still debated physical processes that trigger the explosion and on intrinsic neutron-star physics. These properties include the high space velocities of young neutron stars with average values of several 100 km/s, whose underlying "kick" mechanism is not finally clarified. There are two competing possibilities that could accelerate NSs during their birth: anisotropic ejection of either stellar debris or neutrinos. We here present new evidence from X-ray measurements that chemical elements between silicon and calcium in six young gaseous supernova remnants are preferentially expelled opposite to the direction of neutron star motion. There is no correlation between the kick velocities and magnetic field strengths of these neutron stars. Our results support a hydrodynamic origin of neutron-star kicks connected to asymmetric explosive mass ejection, and they conflict with neutron-star acceleration scenarios that invoke anisotropic neutrino emission caused by particle and nuclear physics in combination with very strong neutron-star magnetic fields.

We revisit the diffuse supernova neutrino background in light of recent systematic studies of stellar core collapse that reveal the quantitative impacts of the progenitor conditions on the collapse process. In general, the dependence of the progenitor on the core-collapse neutrino emission is not monotonic in progenitor initial mass, but we show that it can, at first order, be characterized by the core compactness. For the first time, we incorporate the detailed variations in the neutrino emission over the entire mass range $8$-$100 {\rm M}_\odot$, based on (i) a long-term simulation of the core collapse of a $8.8 {\rm M}_\odot$ O-Ne-Mg core progenitor, (ii) over 100 simulations of iron core collapse to neutron stars, and (iii) half a dozen simulations of core collapse to black holes (the "failed channel"). The fraction of massive stars that undergo the failed channel remains uncertain, but in view of recent simulations which reveal high compactness to be conducive to collapse to black holes, we characterize the failed fraction by considering a threshold compactness above which massive stars collapse to black holes and below which the final remnant is a neutron star. We predict that future detections of the diffuse supernova neutrino background may have the power to reveal this threshold compactness, if its value is relatively small as suggested by interpretations of several recent astronomical observations.

The remnant of neutron star mergers is dense in neutrinos. By employing inputs from one hydrodynamical simulation of a binary neutron star merger remnant with a black hole of $3\ M_\odot$ in the center, dimensionless spin parameter $0.8$ and an accretion torus of $0.3\ M_\odot$, the neutrino emission properties are investigated as the merger remnant evolves. Initially, the local number density of $\bar{\nu}_e$ is larger than that of $\nu_e$ everywhere above the remnant. Then, as the torus approaches self-regulated equilibrium, the local abundance of neutrinos overcomes that of antineutrinos in a funnel around the polar region. The region where the fast pairwise flavor conversions can occur shrinks accordingly as time evolves. Still, we find that fast flavor conversions do affect most of the neutrino-driven ejecta. Assuming that fast flavor conversions lead to flavor equilibration, a significant enhancement of nuclei with mass numbers $A>130$ is found as well as a change of the lanthanide mass fraction by more than a factor of a thousand. Our findings hint towards a potentially relevant role of neutrino flavor oscillations for the prediction of the kilonova (macronova) lightcurves and motivate further work in this direction.

We examine the nucleosynthesis in the innermost, neutrino-processed ejecta (a few $10^{-3}\,M_\odot$) of self-consistent, two-dimensional explosion models of core-collapse supernovae for six progenitor stars with different initial masses. Three models have initial masses near the low-mass end of the supernova range, $8.8\,M_\odot$ (e8.8; electron-capture supernova), $9.6\,M_\odot$ (z9.6), and $8.1\,M_\odot$ (u8.1), with initial metallicities of 1, 0, and $10^{-4}$ times the solar metallicity, respectively. The other three are solar-metallicity models with initial masses of $11.2\,M_\odot$ (s11), $15\,M_\odot$ (s15), and $27\,M_\odot$ (s27). The low-mass models e8.8, z9.6, and u8.1 exhibit high production factors (nucleosynthetic abundances relative to the solar ones) of 100--200 for light trans-iron elements from Zn to Zr. This is associated with appreciable ejection of neutron-rich matter in these models. Remarkably, the nucleosynthetic outcomes for progenitors e8.8 and z9.6 are almost identical, including interesting productions of $^{48}$Ca and $^{60}$Fe, irrespective of their quite different (O-Ne-Mg and Fe) cores prior to collapse. In the more massive models s11, s15, and s27, several proton-rich isotopes of light trans-iron elements, including the $p$-isotope $^{92}$Mo (for s27) are made, up to production factors of $\sim$30. Both electron-capture and core-collapse supernovae near the low-mass end can therefore be dominant contributors to the Galactic inventory of light trans-iron elements from Zn to Zr and probably $^{48}$Ca and live $^{60}$Fe. The innermost ejecta of more massive supernovae may have only sub-dominant contributions to the chemical enrichment of the Galaxy except for $^{92}$Mo.

We present the first self-consistent, three-dimensional (3D) core-collapse supernova simulations performed with the Prometheus-Vertex code for a rotating progenitor star. Besides using the angular momentum of the 15 solar-mass model as obtained in the stellar evolution calculation with an angular frequency of about 0.001 rad/s (spin period of more than 6000 s) at the Si/Si-O interface, we also computed 2D and 3D cases with no rotation and with a ~300 times shorter rotation period and different angular resolutions. In 2D, only the nonrotating and slowly rotating models explode, while rapid rotation prevents an explosion within 500 ms after bounce because of lower radiated neutrino luminosities and mean energies and thus reduced neutrino heating. In contrast, only the fast rotating model develops an explosion in 3D when the Si/Si-O interface collapses through the shock. The explosion becomes possible by the support of a powerful SASI spiral mode, which compensates for the reduced neutrino heating and pushes strong shock expansion in the equatorial plane. Fast rotation in 3D leads to a "two-dimensionalization" of the turbulent energy spectrum (yielding roughly a -3 instead of a -5/3 power-law slope at intermediate wavelengths) with enhanced kinetic energy on the largest spatial scales. We also introduce a generalization of the "universal critical luminosity condition" of Summa et al. (2016) to account for the effects of rotation, and demonstrate its viability for a set of more than 40 core-collapse simulations including 9 and 20 solar-mass progenitors as well as black-hole forming cases of 40 and 75 solar-mass stars to be discussed in forthcoming papers.

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.

We introduce a new, powerful method to constrain properties of neutron stars (NSs). We show that the total mass of GW170817 provides a reliable constraint on the stellar radius if the merger did not result in a prompt collapse as suggested by the interpretation of associated electromagnetic emission. The radius R_1.6 of nonrotating NSs with a mass of 1.6 M_sun can be constrained to be larger than 10.68_{-0.04}^{+0.15} km, and the radius R_max of the nonrotating maximum mass configuration must be larger than 9.60_{-0.03}^{+0.14} km. We point out that detections of future events will further improve these constraints. Moreover, we show that a future event with a signature of a prompt collapse of the merger remnant will establish even stronger constraints on the NS radius from above and the maximum mass M_max of NSs from above. These constraints are particularly robust because they only require a measurement of the chirp mass and a distinction between prompt and delayed collapse of the merger remnant, which may be inferred from the electromagnetic signal or even from the presence/absence of a ringdown gravitational-wave (GW) signal. This prospect strengthens the case of our novel method of constraining NS properties, which is directly applicable to future GW events with accompanying electromagnetic counterpart observations. We emphasize that this procedure is a new way of constraining NS radii from GW detections independent of existing efforts to infer radius information from the late inspiral phase or postmerger oscillations, and it does not require particularly loud GW events.

A large fraction of core-collapse supernovae (CCSNe), 30-50%, are expected to originate from the low-mass end of progenitors with $M_{\rm ZAMS}~= 8-12~M_\odot$. However, degeneracy effects make stellar evolution modelling of such stars challenging, and few predictions for their supernova light curves and spectra have been presented. Here we calculate synthetic nebular spectra of a 9 $M_\odot$ Fe CCSN model exploded with the neutrino mechanism. The model predicts emission lines with FWHM$\sim$1000 km/s, including signatures from each deep layer in the metal core. We compare this model to observations of the three subluminous IIP SNe with published nebular spectra; SN 1997D, SN 2005cs, and SN 2008bk. The prediction of both line profiles and luminosities are in good agreement with SN 1997D and SN 2008bk. The close fit of a model with no tuning parameters provides strong evidence for an association of these objects with low-mass Fe CCSNe. For SN 2005cs, the interpretation is less clear, as the observational coverage ended before key diagnostic lines from the core had emerged. We perform a parameterised study of the amount of explosively made stable nickel, and find that none of these three SNe show the high $^{58}$Ni/$^{56}$Ni ratio predicted by current models of electron capture SNe (ECSNe) and ECSN-like explosions. Combined with clear detection of lines from O and He shell material, these SNe rather originate from Fe core progenitors. We argue that the outcome of self-consistent explosion simulations of low-mass stars, which gives fits to many key observables, strongly suggests that the class of subluminous Type IIP SNe is the observational counterpart of the lowest mass CCSNe.

The evolutionary pathways of core-collapse supernova progenitors at the low-mass end of the spectrum are beset with major uncertainties. In recent years, a variety of evolutionary channels has been discovered in addition to the classical electron capture supernova channel of super-AGB stars. The few available progenitor models at the low-mass end have been studied with great success in supernova simulations as the peculiar density structure makes for robust neutrino-driven explosions in this mass range. Detailed nucleosynthesis calculations have been conducted both for models of electron capture supernovae and low-mass iron core supernovae and revealed an interesting production of the lighter trans-iron elements (such as Zn, Sr, Y, Zr) as well as rare isotopes like Ca-48 and Fe-60. We stress the need to explore the low-mass end of the supernova spectrum further and link various observables to understand the diversity of explosions in this regime.

Double neutron star (DNS) systems represent extreme physical objects and the endpoint of an exotic journey of stellar evolution and binary interactions. Large numbers of DNS systems and their mergers are anticipated to be discovered using the Square-Kilometre-Array searching for radio pulsars and high-frequency gravitational wave detectors (LIGO/VIRGO), respectively. Here we discuss all key properties of DNS systems, as well as selection effects, and combine the latest observational data with new theoretical progress on various physical processes with the aim of advancing our knowledge on their formation. We examine key interactions of their progenitor systems and evaluate their accretion history during the high-mass X-ray binary stage, the common envelope phase and the subsequent Case BB mass transfer, and argue that the first-formed NSs have accreted at most $\sim 0.02\;M_{\odot}$. We investigate DNS masses, spins and velocities, and in particular correlations between spin period, orbital period and eccentricity. Numerous Monte Carlo simulations of the second supernova (SN) events are performed to extrapolate pre-SN stellar properties and probe the explosions. All known close-orbit DNS systems are consistent with ultra-stripped exploding stars. Although their resulting NS kicks are often small, we demonstrate a large spread in kick magnitudes which may, in general, depend on the past interaction history of the exploding star and thus correlate with the NS mass. We analyze and discuss NS kick directions based on our SN simulations. Finally, we discuss the terminal evolution of close-orbit DNS systems until they merge and possibly produce a short $\gamma$-ray burst.

Almost since the beginning, massive stars and their resultant supernovae have played a crucial role in the Universe. These objects produce tremendous amounts of energy and new, heavy elements that enrich galaxies, encourage new stars to form and sculpt the shapes of galaxies we see today. The end of millions of years of massive star evolution and the beginning of hundreds or thousands of years of supernova evolution are separated by a matter of a few seconds, in which some of the most extreme physics found in the Universe causes the explosive and terminal disruption of the star. Key questions remain unanswered in both the studies of how massive stars evolve and the behaviour of supernovae, and it appears the solutions may not lie on just one side of the explosion or the other or in just the domain of the stellar evolution or the supernova astrophysics communities. The need to view massive star evolution and supernovae as continuous phases in a single narrative motivated the Theo Murphy international scientific meeting "Bridging the gap: from massive stars to supernovae" at Chicheley Hall in June 2016, with the specific purpose to simultaneously address the scientific connections between theoretical and observational studies of massive stars and their supernovae, through engaging astronomers from both communities.

Type II-plateau supernovae (SNe IIP) are the most numerous subclass of core-collapse SNe originating from massive stars. In the framework of the neutrino-driven explosion mechanism, we study the SN outburst properties for a red supergiant progenitor model and compare the corresponding light curves with observations of the ordinary Type IIP SN 1999em. Three-dimensional (3D) simulations of (parametrically triggered) neutrino-driven explosions are performed with the (explicit, finite-volume, Eulerian, multifluid hydrodynamics) code PROMETHEUS, using a presupernova model of a 15 Msun star as initial data. At approaching homologous expansion, the hydrodynamical and composition variables of the 3D models are mapped to a spherically symmetric configuration, and the simulations are continued with the (implicit, Lagrangian radiation-hydrodynamics) code CRAB to follow the blast-wave evolution during the SN outburst. Our 3D neutrino-driven explosion model with an explosion energy of about 0.5x10^51 erg produces Ni-56 in rough agreement with the amount deduced from fitting the radioactively powered light-curve tail of SN 1999em. The considered presupernova model, 3D explosion simulations, and light-curve calculations can explain the basic observational features of SN 1999em, except for those connected to the presupernova structure of the outer stellar layers. Our 3D simulations show that the distribution of Ni-rich matter in velocity space is asymmetric with a strong dipole component that is consistent with the observations of SN 1999em. The monotonic luminosity decline from the plateau to the radioactive tail in ordinary SNe IIP is a manifestation of the intense turbulent mixing at the He/H composition interface.

We study the impact of large-scale perturbations from convective shell burning on the core-collapse supernova explosion mechanism using three-dimensional (3D) multi-group neutrino hydrodynamics simulations of an 18 solar mass progenitor. Seed asphericities in the O shell, obtained from a recent 3D model of O shell burning, help trigger a neutrino-driven explosion 330ms after bounce whereas the shock is not revived in a model based on a spherically symmetric progenitor for at least another 300ms. We tentatively infer a reduction of the critical luminosity for shock revival by ~20% due to pre-collapse perturbations. This indicates that convective seed perturbations play an important role in the explosion mechanism in some progenitors. We follow the evolution of the 18 solar mass model into the explosion phase for more than 2s and find that the cycle of accretion and mass ejection is still ongoing at this stage. With a preliminary value of 0.77 Bethe for the diagnostic explosion energy, a baryonic neutron star mass of 1.85 solar masses, a neutron star kick of ~600km/s and a neutron star spin period of ~20ms at the end of the simulation, the explosion and remnant properties are slightly atypical, but still lie comfortably within the observed distribution. Although more refined simulations and a larger survey of progenitors are still called for, this suggests that a solution to the problem of shock revival and explosion energies in the ballpark of observations are within reach for neutrino-driven explosions in 3D.

The merger of two neutron stars may give birth to a long-lived hypermassive neutron star. If it harbours a strong magnetic field of magnetar strength, its spin-down could explain several features of short gamma-ray burst afterglows. The magnetorotational instability (MRI) has been proposed as a mechanism to amplify the magnetic field to the required strength. Previous studies have, however, neglected neutrinos, which may have an important impact on the MRI by inducing a viscosity and drag. We investigate the impact of these neutrino effects on the linear growth of the MRI by applying a local stability analysis to snapshots of a neutron star merger simulation. We find that neutrinos have a significant impact inside the hypermassive neutron star, but have at most a marginal effect in the torus surrounding it. Inside the hypermassive neutron star, the MRI grows in different regimes depending on the radius and on the initial magnetic field strength. For magnetic fields weaker than $10^{13}-10^{14}\,{\rm G}$, the growth rate of the MRI is significantly reduced due to the presence of neutrinos. We conclude that neutrinos should be taken into account when studying the growth of the MRI from realistic initial magnetic fields. Current numerical simulations, which neglect neutrino viscosity, are only consistent, i.e. in the adopted ideal regime, if they start from artificially strong initial magnetic fields above $\sim10^{14}\,{\rm G}$. One should be careful when extrapolating these results to lower initial magnetic fields, where the MRI growth is strongly affected by neutrino viscosity or drag.

Extremely strong magnetic fields of the order of $10^{15}\,{\rm G}$ are required to explain the properties of magnetars, the most magnetic neutron stars. Such a strong magnetic field is expected to play an important role for the dynamics of core-collapse supernovae, and in the presence of rapid rotation may power superluminous supernovae and hypernovae associated to long gamma-ray bursts. The origin of these strong magnetic fields remains, however, obscure and most likely requires an amplification over many orders of magnitude in the protoneutron star. One of the most promising agents is the magnetorotational instability (MRI), which can in principle amplify exponentially fast a weak initial magnetic field to a dynamically relevant strength. We describe our current understanding of the MRI in protoneutron stars and show recent results on its dependence on physical conditions specific to protoneutron stars such as neutrino radiation, strong buoyancy effects and large magnetic Prandtl number.

Most massive stars end their lives in core-collapse supernova explosions and enrich the interstellar medium with explosively nucleosynthesized elements. Following core collapse, the explosion is subject to instabilities as the shock propagates outwards through the progenitor star. Observations of the composition and structure of the innermost regions of a core-collapse supernova provide a direct probe of the instabilities and nucleosynthetic products. SN 1987A in the Large Magellanic Cloud (LMC) is one of very few supernovae for which the inner ejecta can be spatially resolved but are not yet strongly affected by interaction with the surroundings. Our observations of SN 1987A with the Atacama Large Millimeter/submillimeter Array (ALMA) are of the highest resolution to date and reveal the detailed morphology of cold molecular gas in the innermost regions of the remnant. The 3D distributions of carbon and silicon monoxide (CO and SiO) emission differ, but both have a central deficit, or torus-like distribution, possibly a result of radioactive heating during the first weeks ("nickel heating"). The size scales of the clumpy distribution are compared quantitatively to models, demonstrating how progenitor and explosion physics can be constrained.

The spatial and velocity distributions of nuclear species synthesized in the innermost regions of core-collapse supernovae can yield important clues about explosion asymmetries and the operation of the still disputed explosion mechanism. Recent observations of radioactive $^{44}$Ti with high-energy satellite telescopes (Nuclear Spectroscopic Telescope Array [NuSTAR], INTEGRAL) have measured gamma-ray line details, which provide direct evidence of large-scale explosion asymmetries in SN 1987A and in Cassiopeia A (Cas A) even by mapping of the spatial brightness distribution (NuSTAR). Here we discuss a 3D simulation of a neutrino-driven explosion, using a parameterized neutrino engine, whose $^{44}$Ti distribution is mostly concentrated in one hemisphere pointing opposite to the neutron star (NS) kick velocity. Both exhibit intriguing resemblance to the observed morphology of the Cas A remnant, although neither the progenitor nor the explosion was fine-tuned for a perfect match. Our results demonstrate that the asymmetries observed in this remnant can, in principle, be accounted for by a neutrino-driven explosion, and that the high $^{44}$Ti abundance in Cas A may be explained without invoking rapid rotation or a jet-driven explosion, because neutrino-driven explosions generically eject large amounts of high-entropy matter. The recoil acceleration of the NS is connected to mass ejection asymmetries and is opposite to the direction of the stronger explosion, fully compatible with the gravitational tugboat mechanism. Our results also imply that Cas A and SN 1987A could possess similarly "one-sided" Ti and Fe asymmetries, with the difference that Cas A is viewed from a direction with large inclination angle to the NS motion, whereas the NS in SN 1987A should have a dominant velocity component pointing toward us.

Fostered by the possibilities of multi-dimensional computational modeling, in particular the advent of three-dimensional (3D) simulations, our understanding of the neutrino-driven explosion mechanism of core-collapse supernovae (SNe) has experienced remarkable progress over the past decade. First self-consistent, first-principle models have shown successful explosions in 3D, and even failed cases may be cured by moderate changes of the microphysics inside the neutron star (NS), better grid resolution, or more detailed progenitor conditions at the onset of core collapse, in particular large-scale perturbations in the convective Si and O burning shells. 3D simulations have also achieved to follow neutrino-driven explosions continuously from the initiation of the blast wave, through the shock breakout from the progenitor surface, into the radioactively powered evolution of the SN, and towards the free expansion phase of the emerging remnant. Here we present results from such simulations, which form the basis for direct comparisons with observations of SNe and SN remnants in order to derive constraints on the still disputed explosion mechanism. It is shown that predictions based on hydrodynamic instabilities and mixing processes associated with neutrino-driven explosions yield good agreement with measured NS kicks, light-curve properties of SN 1987A, and asymmetries of iron and 44Ti distributions observed in SN 1987A and Cassiopeia A.

We present gravitational wave (GW) signal predictions from four 3D multi-group neutrino hydrodynamics simulations of core-collapse supernovae of progenitors with 11.2 Msun, 20 Msun, and 27 Msun. GW emission in the pre-explosion phase strongly depends on whether the post-shock flow is dominated by the standing accretion shock instability (SASI) or convection and differs considerably from 2D models. SASI activity produces a strong signal component below 250 Hz through asymmetric mass motions in the gain layer and a non-resonant coupling to the proto-neutron star (PNS). Both convection- and SASI-dominated models show GW emission above 250 Hz, but with considerably lower amplitudes than in 2D. This is due to a different excitation mechanism for high-frequency l=2 motions in the PNS surface, which are predominantly excited by PNS convection in 3D. Resonant excitation of high-frequency surface g-modes in 3D by mass motions in the gain layer is suppressed compared to 2D because of smaller downflow velocities and a lack of high-frequency variability in the downflows. In the exploding 20 Msun model, shock revival results in enhanced low-frequency emission due to a change of the preferred scale of the convective eddies in the PNS convection zone. Estimates of the expected excess power in two frequency bands suggests that second-generation detectors will only be able to detect very nearby events, but that third-generation detectors could distinguish SASI- and convection-dominated models at distances of ~10 kpc.

According to recent studies, the collective flavor evolution of neutrinos in core-collapse supernovae depends strongly on the flavor-dependent angular distribution of the local neutrino radiation field, notably on the angular intensity of the electron-lepton number carried by neutrinos. To facilitate further investigations of this subject, we study the energy and angle distributions of the neutrino radiation field computed with the Vertex neutrino-transport code for several spherically symmetric (1D) supernova simulations (of progenitor masses 11.2, 15 and 25 M_sun) and explain how to extract this information from additional models of the Garching group. Beginning in the decoupling region ("neutrino sphere"), the distributions are more and more forward peaked in the radial direction with an angular spread that is largest for $\nu_e$, smaller for $\bar\nu_e$, and smallest for $\nu_x$, where $x=\mu$ or $\tau$. While the energy-integrated $\nu_e$ minus $\bar\nu_e$ angle distribution has a dip in the forward direction, it does not turn negative in any of our investigated cases.

Supernovae are the most powerful cosmic sources of MeV neutrinos. These elementary particles play a crucial role when the evolution of a massive star is terminated by the collapse of its core to a neutron star or a black hole and the star explodes as supernova. The release of electron neutrinos, which are abundantly produced by electron captures, accelerates the catastrophic infall and causes a gradual neutronization of the stellar plasma by converting protons to neutrons as dominant constituents of neutron star matter. The emission of neutrinos and antineutrinos of all flavors carries away the gravitational binding energy of the compact remnant and drives its evolution from the hot initial to the cold final state. The absorption of electron neutrinos and antineutrinos in the surroundings of the newly formed neutron star can power the supernova explosion and determines the conditions in the innermost supernova ejecta, making them an interesting site for the nucleosynthesis of iron-group elements and trans-iron nuclei. In this Chapter the basic neutrino physics in supernova cores and nascent neutron stars will be discussed. This includes the most relevant neutrino production, absorption, and scattering processes, elementary aspects of neutrino transport in dense environments, the characteristic neutrino emission phases with their typical signal features, and the perspectives connected to a measurement of the neutrino signal from a future galactic supernova.

The question why and how core-collapse supernovae (SNe) explode is one of the central and most long-standing riddles of stellar astrophysics. A solution is crucial for deciphering the SN phenomenon, for predicting observable signals such as light curves and spectra, nucleosynthesis, neutrinos, and gravitational waves, for defining the role of SNe in the evolution of galaxies, and for explaining the birth conditions and properties of neutron stars (NSs) and stellar-mass black holes. Since the formation of such compact remnants releases over hundred times more energy in neutrinos than the SN in the explosion, neutrinos can be the decisive agents for powering the SN outburst. According to the standard paradigm of the neutrino-driven mechanism, the energy transfer by the intense neutrino flux to the medium behind the stagnating core-bounce shock, assisted by violent hydrodynamic mass motions (sometimes subsumed by the term "turbulence"), revives the outward shock motion and thus initiates the SN blast. Because of the weak coupling of neutrinos in the region of this energy deposition, detailed, multidimensional hydrodynamic models including neutrino transport and a wide variety of physics are needed to assess the viability of the mechanism. Owing to advanced numerical codes and increasing supercomputer power, considerable progress has been achieved in our understanding of the physical processes that have to act in concert for the success of neutrino-driven explosions. First studies begin to reveal observational implications and avenues to test the theoretical picture by data from individual SNe and SN remnants but also from population-integrated observables. While models will be further refined, a real breakthrough is expected through the next Galactic core-collapse SN, when neutrinos and gravitational waves can be used to probe the conditions deep inside the dying star. (abridged)

Asymmetric mass ejection in the early phase of supernova (SN) explosions can impart a kick velocity to the new-born neutron star (NS). For neutrino-driven explosions the NS acceleration was shown to be mainly caused by the gravitational attraction of the anisotropically expelled inner ejecta, while hydrodynamic forces contribute on a subdominant level, and asymmetric neutrino emission plays only a secondary role. Two- and three-dimensional hydrodynamic simulations demonstrated that this gravitational tug-boat mechanism can explain the observed space velocities of young NSs up to more than 1000 km/s. Here, we discuss how the NS kick depends on the energy, ejecta mass, and asymmetry of the SN explosion, and which role the compactness of the pre-collapse stellar core plays for the momentum transfer to the NS. We also provide simple analytic expressions for the NS velocity in terms of these quantities. Referring to results of hydrodynamic simulations in the literature, we argue why within the discussed scenario of NS acceleration, electron-capture SNe, low-mass Fe-core SNe and ultra-stripped SNe can be expected to have considerably lower intrinsic NS kicks than core-collapse SNe of massive stellar cores. Our basic arguments remain valid also if progenitor stars possess large-scale asymmetries in their convective silicon and oxygen burning layers. Possible scenarios for spin-kick alignment are sketched. Much of our discussion remains on a conceptual and qualitative level and more work is needed on the numerical modeling side to determine the dependences of involved parameters, whose prescriptions will be needed for recipes that can be used to better describe NS kicks in binary evolution and population synthesis studies.

We present the first 3D simulation of the last minutes of oxygen shell burning in an 18 solar mass supernova progenitor up to the onset of core collapse. A moving inner boundary is used to accurately model the contraction of the silicon and iron core according to a 1D stellar evolution model with a self-consistent treatment of core deleptonization and nuclear quasi-equilibrium. The simulation covers the full solid angle to allow the emergence of large-scale convective modes. Due to core contraction and the concomitant acceleration of nuclear burning, the convective Mach number increases to ~0.1 at collapse, and an l=2 mode emerges shortly before the end of the simulation. Aside from a growth of the oxygen shell from 0.51 to 0.56 solar masses due to entrainment from the carbon shell, the convective flow is reasonably well described by mixing length theory, and the dominant scales are compatible with estimates from linear stability analysis. We deduce that artificial changes in the physics, such as accelerated core contraction, can have precarious consequences for the state of convection at collapse. We argue that scaling laws for the convective velocities and eddy sizes furnish good estimates for the state of shell convection at collapse and develop a simple analytic theory for the impact of convective seed perturbations on shock revival in the ensuing supernova. We predict a reduction of the critical luminosity for explosion by 12--24% due to seed asphericities for our 3D progenitor model relative to the case without large seed perturbations.

Neutrino-pair production and annihilation through nucleon-nucleon bremsstrahlung is included in current supernova simulations by rates that are based on the one-pion-exchange approximation. Here we explore the consequences of bremsstrahlung rates based on a modern nuclear interactions for proto-neutron star cooling and the corresponding neutrino emission. We find that despite a reduction of the bremsstrahlung emission by a factor of 2-5 in the neutrinospheric region, models with the improved treatment exhibit only $\lesssim$5% changes of the neutrino luminosities and an increase of $\lesssim$0.7 MeV of the average energies of the radiated neutrino spectra, with the largest effects for the antineutrinos of all flavors and at late times. Overall, the proto-neutron star cooling evolution is slowed down modestly by $\lesssim$0.5-1 s.

We present self-consistent, axisymmetric core-collapse supernova simulations performed with the Prometheus-Vertex code for 18 pre-supernova models in the range of 11-28 solar masses, including progenitors recently investigated by other groups. All models develop explosions, but depending on the progenitor structure, they can be divided into two classes. With a steep density decline at the Si/Si-O interface, the arrival of this interface at the shock front leads to a sudden drop of the mass-accretion rate, triggering a rapid approach to explosion. With a more gradually decreasing accretion rate, it takes longer for the neutrino heating to overcome the accretion ram pressure and explosions set in later. Early explosions are facilitated by high mass-accretion rates after bounce and correspondingly high neutrino luminosities combined with a pronounced drop of the accretion rate and ram pressure at the Si/Si-O interface. Because of rapidly shrinking neutron star radii and receding shock fronts after the passage through their maxima, our models exhibit short advection time scales, which favor the efficient growth of the standing accretion-shock instability. The latter plays a supportive role at least for the initiation of the re-expansion of the stalled shock before runaway. Taking into account the effects of turbulent pressure in the gain layer, we derive a generalized condition for the critical neutrino luminosity that captures the explosion behavior of all models very well. We validate the robustness of our findings by testing the influence of stochasticity, numerical resolution, and approximations in some aspects of the microphysics.

Long-lived radioactive nuclei play an important role as nucleo-cosmochronometers and as cosmic tracers of nucleosynthetic source activity. In particular nuclei in the actinide region like thorium, uranium, and plutonium can testify to the enrichment of an environment by the still enigmatic astrophysical sources that are responsible for the production of neutron-rich nuclei by the rapid neutron-capture process (r-process). Supernovae and merging neutron-star (NS) or NS-black hole binaries are considered as most likely sources of the r-nuclei. But arguments in favour of one or the other or both are indirect and make use of assumptions; they are based on theoretical models with remaining simplifications and shortcomings. An unambiguous observational determination of a production event is still missing. In order to facilitate searches in this direction, e.g.\ by looking for radioactive tracers in stellar envelopes, the interstellar medium or terrestrial reservoirs, we provide improved theoretical estimates and corresponding uncertainty ranges for the actinide production ($^{232}$Th, $^{235,236,238}$U, $^{237}$Np, $^{244}$Pu, and $^{247}$Cm) in neutrino-driven winds of core-collapse supernovae. Since state-of-the-art supernova models do not yield r-process viable conditions ---but still lack, for example, the effects of strong magnetic fields--- we base our investigation on a simple analytical, Newtonian, adiabatic and steady-state wind model and consider the superposition of a large number of contributing components, whose nucleosynthesis-relevant parameters (mass weight, entropy, expansion time scale, and neutron excess) are constrained by the assumption that the integrated wind nucleosynthesis closely reproduces the solar system distribution of r-process elements. We also test the influence of uncertain nuclear physics.

Nonspherical mass motions are a generic feature of core-collapse supernovae, and hydrodynamic instabilities play a crucial role for the explosion mechanism. First successful neutrino-driven explosions could be obtained with self-consistent, first-principle simulations in three spatial dimensions (3D). But 3D models tend to be less prone to explosion than corresponding axisymmetric (2D) ones. This has been explained by 3D turbulence leading to energy cascading from large to small spatial scales, inversely to the 2D case, thus disfavoring the growth of buoyant plumes on the largest scales. Unless the inertia to explode simply reflects a lack of sufficient resolution in relevant regions, it suggests that some important aspect may still be missing for robust and sufficiently energetic neutrino-powered explosions. Such deficits could be associated with progenitor properties like rotation, magnetic fields or pre-collapse perturbations, or with microphysics that could lead to an enhancement of neutrino heating behind the shock. 3D simulations have also revealed new phenomena that are not present in 2D, for example spiral modes of the standing accretion shock instability (SASI) and a stunning dipolar lepton-emission self-sustained asymmetry (LESA). Both impose time- and direction-dependent variations on the detectable neutrino signal. The understanding of these effects and of their consequences is still in its infancy.

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

Nucleosynthesis, light curves, explosion energies, and remnant masses are calculated for a grid of supernovae resulting from massive stars with solar metallicity and masses from 9.0 to 120 solar masses. The full evolution is followed using an adaptive reaction network of up to 2000 nuclei. A novel aspect of the survey is the use of a one-dimensional neutrino transport model for the explosion. This explosion model has been calibrated to give the observed energy for SN 1987A, using several standard progenitors, and for the Crab supernova using a 9.6 solar mass progenitor. As a result of using a calibrated central engine, the final kinetic energy of the supernova is variable and sensitive to the structure of the presupernova star. Many progenitors with extended core structures do not explode, but become black holes, and the masses of exploding stars do not form a simply connected set. The resulting nucleosynthesis agrees reasonably well with the sun provided that a reasonable contribution from Type Ia supernovae is also allowed, but with a deficiency of light s-process isotopes. The resulting neutron star IMF has a mean gravitational mass near 1.4 solar masses. The average black hole mass is about 9 solar masses if only the helium core implodes, and 14 solar masses if the entire presupernova star collapses. Only ~10% of supernovae come from stars over 20 solar masses and some of these are Type Ib or Ic. Some useful systematics of Type IIp light curves are explored.

Neutrinos play a crucial role in the collapse and explosion of massive stars, governing the infall dynamics of the stellar core, triggering and fueling the explosion and driving the cooling and deleptonization of the newly formed neutron star. Due to their role neutrinos carry information from the heart of the explosion and, due to their weakly interacting nature, offer the only direct probe of the dynamics and thermodynamics at the center of a supernova. In this paper, we review the present status of modelling the neutrino physics and signal formation in collapsing and exploding stars. We assess the capability of current and planned large underground neutrino detectors to yield faithful information of the time and flavor dependent neutrino signal from a future Galactic supernova. We show how the observable neutrino burst would provide a benchmark for fundamental supernova physics with unprecedented richness of detail. Exploiting the treasure of the measured neutrino events requires a careful discrimination of source-generated properties from signal features that originate on the way to the detector. As for the latter, we discuss self-induced flavor conversions associated with neutrino-neutrino interactions that occur in the deepest stellar regions; matter effects that modify the pattern of flavor conversions in the dynamical stellar envelope; neutrino-oscillation signatures that result from structural features associated with the shock-wave propagation as well as turbulent mass motions in post-shock layers. Finally, we highlight our current understanding of the formation of the diffuse supernova neutrino background and we analyse the perspectives for a detection of this relic signal that integrates the contributions from all past core-collapse supernovae in the Universe.

The coalescence of a neutron-star binary is likely to result in the formation of a neutron-star merger remnant for a large range of binary mass configurations. The massive merger remnant shows strong oscillations, which are excited by the merging process, and emits gravitational waves. Here we discuss possibilities and prospects of inferring unknown stellar properties of neutron stars by the detection of postmerger gravitational-wave emission, which thus leads to constraints of the equation of state of high-density matter. In particular, the dominant oscillation frequency of the postmerger remnant provides tight limits to neutron-star radii. We mention first steps towards a practical implementation of future gravitational-wave searches for the postmerger emission. Moreover, we outline possibilities to estimate the unknown maximum mass of nonrotating neutron stars from such types of measurements. Finally, we review the origin and scientific implications of secondary peaks in the gravitational-wave spectrum of neutron-star mergers and differences in the dynamical behavior of the postmerger remnant depending on the binary configuration. These considerations lead to a unified picture of the post-merger gravitational-wave emission and the post-merger dynamics.

We present the first special relativistic, axisymmetric hydrodynamic simulations of black hole-torus systems (approximating general relativistic gravity) as remnants of binary-neutron star (NS-NS) and neutron star-black hole (NS-BH) mergers, in which the viscously driven evolution of the accretion torus is followed with self-consistent energy-dependent neutrino transport and the interaction with the cloud of dynamical ejecta expelled during the NS-NS merging is taken into account. The modeled torus masses, BH masses and spins, and the ejecta masses, velocities, and spatial distributions are adopted from relativistic merger simulations. We find that energy deposition by neutrino annihilation can accelerate outflows with initially high Lorentz factors along polar low-density funnels, but only in mergers with extremely low baryon pollution in the polar regions. NS-BH mergers, where polar mass ejection during the merging phase is absent, provide sufficiently baryon-poor environments to enable neutrino-powered, ultrarelativistic jets with terminal Lorentz factors above 100 and considerable dynamical collimation, favoring short gamma-ray bursts (sGRBs), although their typical energies and durations might be too small to explain the majority of events. In the case of NS-NS mergers, however, neutrino emission of the accreting and viscously spreading torus is too short and too weak to yield enough energy for the outflows to break out from the surrounding ejecta shell as highly relativistic jets. We conclude that neutrino annihilation alone cannot power sGRBs from NS-NS mergers.

Thus far, judging the fate of a massive star (either a neutron star (NS) or a black hole) solely by its structure prior to core collapse has been ambiguous. Our work and previous attempts find a non-monotonic variation of successful and failed supernovae with zero-age main-sequence mass, for which no single structural parameter can serve as a good predictive measure. However, we identify two parameters computed from the pre-collapse structure of the progenitor, which in combination allow for a clear separation of exploding and non-exploding cases with only few exceptions (~1-2.5%) in our set of 621 investigated stellar models. One parameter is M4, defining the normalized enclosed mass for a dimensionless entropy per nucleon of s=4, and the other is mu4 = d(m/M_sun)/d(r/1000 km) at s=4, being the normalized mass-derivative at this location. The two parameters mu4 and M4*mu4 can be directly linked to the mass-infall rate, Mdot, of the collapsing star and the electron-type neutrino luminosity of the accreting proto-NS, L_nue ~ M_ns*Mdot, which play a crucial role in the "critical luminosity" concept for the theoretical description of neutrino-driven explosions as runaway phenomenon of the stalled accretion shock. All models were evolved employing the approach of Ugliano et al. for simulating neutrino-driven explosions in spherical symmetry. The neutrino emission of the accretion layer is approximated by a gray transport solver, while the uncertain neutrino emission of the 1.1 M_sun proto-NS core is parametrized by an analytic model. The free parameters connected to the core-boundary prescription are calibrated to reproduce the observables of Supernova 1987A for five different progenitor models.

Remnants of neutron-star mergers are essentially massive, hot, differentially rotating neutron stars, which are initially strongly oscillating. They represent a unique probe for high-density matter because the oscillations are detectable via gravitational-wave measurements and are strongly dependent on the equation of state. The impact of the equation of state is apparent in the frequency of the dominant oscillation mode of the remnant. For a fixed total binary mass a tight relation between the dominant postmerger frequency and the radii of nonrotating neutron stars exists. Inferring observationally the dominant postmerger frequency thus determines neutron star radii with high accuracy of the order of a few hundred meters. By considering symmetric and asymmetric binaries of the same chirp mass, we show that the knowledge of the binary mass ratio is not critical for this kind of radius measurements. We summarize different possibilities to deduce the maximum mass of nonrotating neutron stars. We clarify the nature of the three most prominent features of the postmerger gravitational-wave spectrum and argue that the merger remnant can be considered to be a single, isolated, self-gravitating object that can be described by concepts of asteroseismology. The understanding of the different mechanisms shaping the gravitational-wave signal yields a physically motivated analytic model of the gravitational-wave emission, which may form the basis for template-based gravitational-wave data analysis. We explore the observational consequences of a scenario of two families of compact stars including hadronic and quark stars. We find that this scenario leaves a distinctive imprint on the postmerger gravitational-wave signal. In particular, a strong discontinuity in the dominant postmerger frequency as function of the total mass will be a strong indication for two families of compact stars. (abridged)

We present a comprehensive study of r-process element nucleosynthesis in the ejecta of compact binary mergers (CBMs) and their relic black-hole (BH)-torus systems. The evolution of the BH-accretion tori is simulated for seconds with a Newtonian hydrodynamics code including viscosity effects, pseudo-Newtonian gravity for rotating BHs, and an energy-dependent two-moment closure scheme for the transport of electron neutrinos and antineutrinos. The investigated cases are guided by relativistic double neutron star (NS-NS) and NS-BH merger models, producing ~3-6 Msun BHs with rotation parameters of A~0.8 and tori of 0.03-0.3 Msun. Our nucleosynthesis analysis includes the dynamical (prompt) ejecta expelled during the CBM phase and the neutrino and viscously driven outflows of the relic BH-torus systems. While typically ~20-25% of the initial accretion-torus mass are lost by viscously driven outflows, neutrino-powered winds contribute at most another ~1%, but neutrino heating enhances the viscous ejecta significantly. Since BH-torus ejecta possess a wide distribution of electron fractions (0.1-0.6) and entropies, they produce heavy elements from A~80 up to the actinides, with relative contributions of A>130 nuclei being subdominant and sensitively dependent on BH and torus masses and the exact treatment of shear viscosity. The combined ejecta of CBM and BH-torus phases can reproduce the solar abundances amazingly well for A>90. Varying contributions of the torus ejecta might account for observed variations of lighter elements with 40<Z<56 relative to heavier ones, and a considerable reduction of the prompt ejecta compared to the torus ejecta, e.g. in highly asymmetric NS-BH mergers, might explain the composition of heavy-element deficient stars.

We have performed r-process calculations for matter ejected dynamically in neutron star mergers based on a complete set of trajectories from a three-dimensional relativistic smoothed particle hydrodynamic simulation. Our calculations consider an extended nuclear network, including spontaneous, $\beta$- and neutron-induced fission and adopting fission yield distributions from the ABLA code. We have studied the sensitivity of the r-process abundances to nuclear masses by using different models. Most of the trajectories, corresponding to 90% of the ejected mass, follow a relatively slow expansion allowing for all neutrons to be captured. The resulting abundances are very similar to each other and reproduce the general features of the observed r-process abundance (the second and third peaks, the rare-earth peak and the lead peak) for all mass models as they are mainly determined by the fission yields. We find distinct differences in the abundance yields at and just above the third peak, which can be traced back to different predictions of neutron separation energies for r-process nuclei around neutron number $N=130$. The remaining trajectories, which contribute 10% by mass to the total integrated abundances, follow such a fast expansion that the r process does not use all the neutrons. This also leads to a larger variation of abundances among trajectories as fission does not dominate the r-process dynamics. The total integrated abundances are dominated by contributions from the slow abundances and hence reproduce the general features of the observed r-process abundances. We find that at timescales of weeks relevant for kilonova light curve calculations, the abundance of actinides is larger than the one of lanthanides. Hence actinides can be even more important than lanthanides to determine the photon opacities under kilonova conditions. (Abridged)

We present the new code ALCAR developed to model multidimensional, multi energy-group neutrino transport in the context of supernovae and neutron-star mergers. The algorithm solves the evolution equations of the 0th- and 1st-order angular moments of the specific intensity, supplemented by an algebraic relation for the 2nd-moment tensor to close the system. The scheme takes into account frame-dependent effects of order O(v/c) as well as the most important types of neutrino interactions. The transport scheme is significantly more efficient than a multidimensional solver of the Boltzmann equation, while it is more accurate and consistent than the flux-limited diffusion method. The finite-volume discretization of the essentially hyperbolic system of moment equations employs methods well-known from hydrodynamics. For the time integration of the potentially stiff moment equations we employ a scheme in which only the local source terms are treated implicitly, while the advection terms are kept explicit, thereby allowing for an efficient computational parallelization of the algorithm. We investigate various problem setups in one and two dimensions to verify the implementation and to test the quality of the algebraic closure scheme. In our most detailed test, we compare a fully dynamic, one-dimensional core-collapse simulation with two published calculations performed with well-known Boltzmann-type neutrino-hydrodynamics codes and we find very satisfactory agreement.

During the accretion phase of a core-collapse supernova (SN), the deleptonization flux has recently been found to develop a global dipole pattern (LESA---Lepton Emission Self-sustained Asymmetry). The $\nu_e$ minus $\bar\nu_e$ flux essentially vanishes in one direction, potentially facilitating self-induced flavor conversion. On the other hand, below the stalled shock wave, self-induced flavor conversion is typically suppressed by multi-angle matter effects, preventing any impact of flavor conversion on SN explosion dynamics. In a schematic model of SN neutrino fluxes, we study the impact of modified $\bar\nu_e$-$\nu_e$ flux asymmetries on collective flavor conversion. In the parameter space consisting of matter density and effective neutrino density, the region of instability with regard to self-induced flavor conversion is much larger for a vanishing lepton number flux, yet this modification does not intersect a realistic SN profile. Therefore, it appears that, even in the presence of LESA, self-induced flavor conversion remains suppressed below the shock front.

Interactions with neutrons and protons play a crucial role for the neutrino opacity of matter in the supernova core. Their current implementation in many simulation codes, however, is rather schematic and ignores not only modifications for the correlated nuclear medium of the nascent neutron star, but also free-space corrections from nucleon recoil, weak magnetism or strange quarks, which can easily add up to changes of several 10% for neutrino energies in the spectral peak. In the Garching supernova simulations with the Prometheus-Vertex code, such sophistications have been included for a long time except for the strange-quark contributions to the nucleon spin, which affect neutral-current neutrino scattering. We demonstrate on the basis of a 20 M_sun progenitor star that a moderate strangeness-dependent contribution of g_a^s = -0.2 to the axial-vector coupling constant g_a = 1.26 can turn an unsuccessful three-dimensional (3D) model into a successful explosion. Such a modification is in the direction of current experimental results and reduces the neutral-current scattering opacity of neutrons, which dominate in the medium around and above the neutrinosphere. This leads to increased luminosities and mean energies of all neutrino species and strengthens the neutrino-energy deposition in the heating layer. Higher nonradial kinetic energy in the gain layer signals enhanced buoyancy activity that enables the onset of the explosion at ~300 ms after bounce, in contrast to the model with vanishing strangeness contributions to neutrino-nucleon scattering. Our results demonstrate the close proximity to explosion of the previously published, unsuccessful 3D models of the Garching group.

The well-studied type IIP SN 1987A, produced by the explosion of a blue supergiant (BSG) star, is a touchstone for massive-star evolution, simulations of neutrino-driven explosions, and modeling of light curves and spectra. In the framework of the neutrino-driven mechanism, we study the dependence of explosion properties on the structure of four different BSGs and compare the corresponding light curves with observations of SN 1987A. We perform 3D simulations with the PROMETHEUS code until about one day and map the results to the 1D code CRAB for the light curve calculations. All of our 3D models with explosion energies compatible with SN 1987A produce 56Ni in rough agreement with the amount deduced from fitting the radioactively powered light-curve tail. One of the progenitors yields maximum velocities of ~3000 km/s for the bulk of ejected 56Ni, consistent with observations. In all of our models inward mixing of hydrogen during the 3D evolution leads to minimum H-velocities below 100 km/s, in good agreement with spectral observations. The considered BSG models, 3D explosion simulations, and light-curve calculations can thus explain basic observational features of SN 1987A. However, all progenitors have too large pre-SN radii to reproduce the narrow initial luminosity peak, and the structure of their outer layers is not suitable to match the observed light curve during the first 30-40 days. Only one stellar model has a structure of the He core and the He/H composition interface that enables sufficient outward mixing of 56Ni and inward mixing of hydrogen to produce a good match of the dome-like shape of the observed light-curve maximum. But this model falls short of the He-core mass of 6 Msun inferred from the absolute luminosity of the pre-SN star. The lack of an adequate pre-SN model for SN 1987A is a pressing challenge for the theory of massive-star evolution. (Abridged)

Neutrino oscillations, especially to light sterile states, can affect the nucleosynthesis yields because of their possible feedback effect on the electron fraction (Ye). For the first time, we perform nucleosynthesis calculations for neutrino-driven wind trajectories from the neutrino-cooling phase of an 8.8 Msun electron-capture supernova, whose hydrodynamic evolution was computed in spherical symmetry with sophisticated neutrino transport and whose Ye evolution was post-processed by including neutrino oscillations both between active and active-sterile flavors. We also take into account the alpha-effect as well as weak magnetism and recoil corrections in the neutrino absorption and emission processes. We observe effects on the Ye evolution which depend in a subtle way on the relative radial positions of the sterile MSW resonances, of collective flavor transformations, and on the formation of alpha particles. For the adopted supernova progenitor, we find that neutrino oscillations, also to a sterile state with eV-mass, do not significantly affect the element formation and in particular cannot make the post-explosion wind outflow neutron rich enough to activate a strong r-process. Our conclusions become even more robust when, in order to mimic equation-of-state dependent corrections due to nucleon potential effects in the dense-medium neutrino opacities, six cases with reduced Ye in the wind are considered. In these cases, despite the conversion of active neutrinos to sterile neutrinos, Ye increases or is not significantly lowered compared to the values obtained without oscillations and active flavor transformations. This is a consequence of a complicated interplay between sterile-neutrino production, neutrino-neutrino interactions, and alpha-effect.

We investigate beta-interactions of free nucleons and their impact on the electron fraction (Y_e) and r-process nucleosynthesis in ejecta characteristic of binary neutron star mergers (BNSMs). For that we employ trajectories from a relativistic BNSM model to represent the density-temperature evolutions in our parametric study. In the high-density environment, positron captures decrease the neutron richness at the high temperatures predicted by the hydrodynamic simulation. Circumventing the complexities of modelling three-dimensional neutrino transport, (anti)neutrino captures are parameterized in terms of prescribed neutrino luminosities and mean energies, guided by published results and assumed as constant in time. Depending sensitively on the adopted neutrino-antineutrino luminosity ratio, neutrino processes increase Y_e to values between 0.25 and 0.40, still allowing for a successful r-process compatible with the observed solar abundance distribution and a significant fraction of the ejecta consisting of r-process nuclei. If the electron neutrino luminosities and mean energies are relatively large compared to the antineutrino properties, the mean Y_e might reach values >0.40 so that neutrino captures seriously compromise the success of the r-process. In this case, the r-abundances remain compatible with the solar distribution, but the total amount of ejected r-material is reduced to a few percent, because the production of iron-peak elements is favored. Proper neutrino physics, in particular also neutrino absorption, have to be included in BNSM simulations before final conclusions can be drawn concerning r-processing in this environment and concerning observational consequences like kilonovae, whose peak brightness and color temperature are sensitive to the composition-dependent opacity of the ejecta.

The oscillations of a merger remnant forming after the coalescence of two neutron stars are very characteristic for the high-density equation of state. The dominant oscillation frequency occurs as a pronounced peak in the kHz range of the gravitational-wave spectrum. We describe how the dominant oscillation frequency of the remnant can be employed to infer the radii of non-rotating neutron stars.

Since core-collapse supernova simulations still struggle to produce robust neutrino-driven explosions in 3D, it has been proposed that asphericities caused by convection in the progenitor might facilitate shock revival by boosting the activity of non-radial hydrodynamic instabilities in the post-shock region. We investigate this scenario in depth using 42 relativistic 2D simulations with multi-group neutrino transport to examine the effects of velocity and density perturbations in the progenitor for different perturbation geometries that obey fundamental physical constraints (like the anelastic condition). As a framework for analysing our results, we introduce semi-empirical scaling laws relating neutrino heating, average turbulent velocities in the gain region, and the shock deformation in the saturation limit of non-radial instabilities. The squared turbulent Mach number, <Ma^2>, reflects the violence of aspherical motions in the gain layer, and explosive runaway occurs for <Ma^2>~0.3, corresponding to a reduction of the critical neutrino luminosity by ~25% compared to 1D. In the light of this theory, progenitor asphericities aid shock revival mainly by creating anisotropic mass flux onto the shock: Differential infall efficiently converts velocity perturbations in the progenitor into density perturbations (Delta rho/rho) at the shock of the order of the initial convective Mach number Ma. The anisotropic mass flux and ram pressure deform the shock and thereby amplify post-shock turbulence. Large-scale (l=2,l=1) modes prove most conducive to shock revival, whereas small-scale perturbations require unrealistically high convective Mach numbers. Initial density perturbations in the progenitor are only of order Ma^2 and therefore play a subdominant role.

We present the first successful simulation of a neutrino-driven supernova explosion in three dimensions (3D), using the Prometheus-Vertex code with an axis-free Yin-Yang grid and a sophisticated treatment of three-flavor, energy-dependent neutrino transport. The progenitor is a nonrotating, zero-metallicity 9.6 Msun star with an iron core. While in spherical symmetry outward shock acceleration sets in later than 300 ms after bounce, a successful explosion starts at ~130 ms postbounce in two dimensions (2D). The 3D model explodes at about the same time but with faster shock expansion than in 2D and a more quickly increasing and roughly 10 percent higher explosion energy of >10^50 erg. The more favorable explosion conditions in 3D are explained by lower temperatures and thus reduced neutrino emission in the cooling layer below the gain radius. This moves the gain radius inward and leads to a bigger mass in the gain layer, whose larger recombination energy boosts the explosion energy in 3D. These differences are caused by less coherent, less massive, and less rapid convective downdrafts associated with postshock convection in 3D. The less violent impact of these accretion downflows in the cooling layer produces less shock heating and therefore diminishes energy losses by neutrino emission. We thus have, for the first time, identified a reduced mass accretion rate, lower infall velocities, and a smaller surface filling factor of convective downdrafts as consequences of 3D postshock turbulence that facilitate neutrino-driven explosions and strengthen them compared to the 2D case.

We present 3D simulations of core-collapse supernovae from blast-wave initiation by the neutrino-driven mechanism to shock breakout from the stellar surface, considering two 15 Msun red supergiants (RSG) and two blue supergiants (BSG) of 15 Msun and 20 Msun. We demonstrate that the metal-rich ejecta in homologous expansion still carry fingerprints of asymmetries at the beginning of the explosion, but the final metal distribution is massively affected by the detailed progenitor structure. The most extended and fastest metal fingers and clumps are correlated with the biggest and fastest-rising plumes of neutrino-heated matter, because these plumes most effectively seed the growth of Rayleigh-Taylor (RT) instabilities at the C+O/He and He/H composition-shell interfaces after the passage of the SN shock. The extent of radial mixing, global asymmetry of the metal-rich ejecta, RT-induced fragmentation of initial plumes to smaller-scale fingers, and maximal Ni and minimal H velocities do not only depend on the initial asphericity and explosion energy (which determine the shock and initial Ni velocities) but also on the density profiles and widths of C+O core and He shell and on the density gradient at the He/H transition, which lead to unsteady shock propagation and the formation of reverse shocks. Both RSG explosions retain a great global metal asymmetry with pronounced clumpiness and substructure, deep penetration of Ni fingers into the H-envelope (with maximum velocities of 4000-5000 km/s for an explosion energy around 1.5 bethe) and efficient inward H-mixing. While the 15 Msun BSG shares these properties (maximum Ni speeds up to ~3500 km/s), the 20 Msun BSG develops a much more roundish geometry without pronounced metal fingers (maximum Ni velocities only ~2200 km/s) because of reverse-shock deceleration and insufficient time for strong RT growth and fragmentation at the He/H interface.

The magnetorotational instability (MRI) is a promising mechanism to amplify the magnetic field in fast rotating protoneutron stars. The diffusion of neutrinos trapped in the PNS induces a transport of momentum, which can be modelled as a viscosity on length-scales longer than the neutrino mean free path. This neutrino-viscosity can slow down the growth of MRI modes to such an extent that a minimum initial magnetic field strength of $\gtrsim 10^{12}\, {\rm G}$ is needed for the MRI to grow on a sufficiently short time-scale to potentially affect the explosion. It is uncertain whether the magnetic field of fast rotating progenitor cores is strong enough to yield such an initial magnetic field in PNS. At MRI wavelengths shorter than the neutrino mean free path, on the other hand, neutrino radiation does not act as a viscosity but rather induces a drag on the velocity with a damping rate independent of the wavelength. We perform a linear analysis of the MRI in this regime, and apply our analytical results to the PNS structure from a one-dimensional numerical simulation. We show that in the outer layers of the PNS, the MRI can grow from weak magnetic fields at wavelengths shorter than the neutrino mean free path, while deeper in the PNS MRI growth takes place in the viscous regime and requires a minimum magnetic field strength.

We present the first comprehensive study of r-process element nucleosynthesis in the ejecta of compact binary mergers (CBMs) and their relic black-hole (BH)-torus systems. The evolution of the BH-accretion tori is simulated for seconds with a Newtonian hydrodynamics code including viscosity effects, pseudo-Newtonian gravity for rotating BHs, and an energy-dependent two-moment closure scheme for the transport of electron neutrinos and antineutrinos. The investigated cases are guided by relativistic double neutron star (NS-NS) and NS-BH merger models, producing ~3-6 Msun BHs with rotation parameters of A~0.8 and tori of 0.03-0.3 Msun. Our nucleosynthesis analysis includes the dynamical (prompt) ejecta expelled during the CBM phase and the neutrino and viscously driven outflows of the relic BH-torus systems. While typically ~20-25% of the initial accretion-torus mass are lost by viscously driven outflows, neutrino-powered winds contribute at most another ~1%, but neutrino heating enhances the viscous ejecta significantly. Since BH-torus ejecta possess a wide distribution of electron fractions (0.1-0.6) and entropies, they produce heavy elements from A~80 up to the actinides, with relative contributions of A>130 nuclei being subdominant and sensitively dependent on BH and torus masses and the exact treatment of shear viscosity. The combined ejecta of CBM and BH-torus phases can reproduce the solar abundances amazingly well for A>90. Varying contributions of the torus ejecta might account for observed variations of lighter elements with 40<Z<56 relative to heavier ones, and a considerable reduction of the prompt ejecta compared to the torus ejecta, e.g. in highly asymmetric NS-BH mergers, might explain the composition of heavy-element deficient stars.

During the shock-wave propagation in a core-collapse supernova (SN), matter turbulence may affect neutrino flavor conversion probabilities. Such effects have been usually studied by adding parametrized small-scale random fluctuations (with arbitrary amplitude) on top of coarse, spherically symmetric matter density profiles. Recently, however, two-dimensional (2D) SN models have reached a space resolution high enough to directly trace anisotropic density profiles, down to scales smaller than the typical neutrino oscillation length. In this context, we analyze the statistical properties of a large set of SN matter density profiles obtained in a high-resolution 2D simulation, focusing on a post-bounce time (2 s) suited to study shock-wave effects on neutrino propagation on scales as small as O(100) km and possibly below. We clearly find the imprint of a broken (Kolmogorov-Kraichnan) power-law structure, as generically expected in 2D turbulence spectra. We then compute the flavor evolution of SN neutrinos along representative realizations of the turbulent matter density profiles, and observe no or modest damping of the neutrino crossing probabilities on their way through the shock wave. In order to check the effect of possibly unresolved fluctuations at scales below O(100) km, we also apply a randomization procedure anchored to the power spectrum calculated from the simulation, and find consistent results within \pm 1 sigma fluctuations. These results show the importance of anchoring turbulence effects on SN neutrinos to realistic, fine-grained SN models.

We derive constraints on the time-averaged event rate of neutron star-black hole (NS-BH) mergers by using estimates of the population-integrated production of heavy rapid neutron-capture (r-process) elements with nuclear mass numbers A > 140 by such events in comparison to the Galactic repository of these chemical species. Our estimates are based on relativistic hydrodynamical simulations convolved with theoretical predictions of the binary population. This allows us to determine a strict upper limit of the average NS-BH merger rate of ~6*10^-5 per year. We quantify the uncertainties of this estimate to be within factors of a few mostly because of the unknown BH spin distribution of such systems, the uncertain equation of state of NS matter, and possible errors in the Galactic content of r-process material. Our approach implies a correlation between the merger rates of NS-BH binaries and of double NS systems. Predictions of the detection rate of gravitational-wave signals from such compact-object binaries by Advanced LIGO and Advanced Virgo on the optimistic side are incompatible with the constraints set by our analysis.

The neutrino emission characteristics of the first full-scale three-dimensional supernova simulations with sophisticated three-flavor neutrino transport for three models with masses 11.2, 20 and 27 M_sun are evaluated in detail. All the studied progenitors show the expected hydrodynamical instabilities in the form of large-scale convective overturn. In addition, the recently identified LESA phenomenon (lepton-number emission self-sustained asymmetry) is generic for all our cases. Pronounced SASI (standing accretion-shock instability) activity appears in the 20 and 27 M_sun cases, partly in the form of a spiral mode, inducing large but direction and flavor-dependent modulations of neutrino emission. These modulations can be clearly identified in the existing IceCube and future Hyper-Kamiokande detectors, depending on distance and detector location relative to the main SASI sloshing direction.

During the stalled-shock phase of our 3D hydrodynamical core-collapse simulations with energy-dependent, 3-flavor neutrino transport, the lepton-number flux (nue minus antinue) emerges predominantly in one hemisphere. This novel, spherical-symmetry breaking neutrino-hydrodynamical instability is termed LESA for "Lepton-number Emission Self-sustained Asymmetry." While the individual nue and antinue fluxes show a pronounced dipole pattern, the heavy-flavor neutrino fluxes and the overall luminosity are almost spherically symmetric. Initially, LESA seems to develop stochastically from convective fluctuations, it exists for hundreds of milliseconds or more, and it persists during violent shock sloshing associated with the standing accretion shock instability. The nue minus antinue flux asymmetry originates mainly below the neutrinosphere in a region of pronounced proto-neutron star (PNS) convection, which is stronger in the hemisphere of enhanced lepton-number flux. On this side of the PNS, the mass-accretion rate of lepton-rich matter is larger, amplifying the lepton-emission asymmetry, because the spherical stellar infall deflects on a dipolar deformation of the stalled shock. The increased shock radius in the hemisphere of less mass accretion and minimal lepton-number flux (antinue flux maximum) is sustained by stronger convection on this side, which is boosted by stronger neutrino heating because the average antinue energy is higher than the average nue energy. Asymmetric heating thus supports the global deformation despite extremely nonstationary convective overturn behind the shock. While these different elements of LESA form a consistent picture, a full understanding remains elusive at present. There may be important implications for neutrino-flavor oscillations, the neutron-to-proton ratio in the neutrino-heated supernova ejecta, and neutron-star kicks, which remain to be explored.

Considering general relativistic, two-dimensional (2D) supernova (SN) explosion models of progenitor stars between 8.1 and 27 solar masses, we systematically analyze the properties of the neutrino emission from core collapse and bounce to the post-explosion phase. The models were computed with the Vertex-CoCoNuT code, using three-flavor, energy-dependent neutrino transport in the ray-by-ray-plus approximation. Our results confirm the close similarity of the mean energies of electron antineutrinos and heavy-lepton neutrinos and even their crossing during the accretion phase for stars with M>10 M_sun as observed in previous 1D and 2D simulations with state-of-the-art neutrino transport. We establish a roughly linear scaling of the electron antineutrino mean energy with the proto-neutron star (PNS) mass, which holds in time as well as for different progenitors. Convection inside the PNS affects the neutrino emission on the 10-20% level, and accretion continuing beyond the onset of the explosion prevents the abrupt drop of the neutrino luminosities seen in artificially exploded 1D models. We demonstrate that a wavelet-based time-frequency analysis of SN neutrino signals in IceCube will offer sensitive diagnostics for the SN core dynamics up to at least ~10kpc distance. Strong, narrow-band signal modulations indicate quasi-periodic shock sloshing motions due to the standing accretion shock instability (SASI), and the frequency evolution of such "SASI neutrino chirps" reveals shock expansion or contraction. The onset of the explosion is accompanied by a shift of the modulation frequency below 40-50Hz, and post-explosion, episodic accretion downflows will be signaled by activity intervals stretching over an extended frequency range in the wavelet spectrogram.

The equation of state plays a critical role in the physics of the merger of two neutron stars. Recent numerical simulations with microphysical equation of state suggest the outcome of such events depends on the mass of the neutron stars. For less massive systems, simulations favor the formation of a hypermassive, quasi-stable neutron star, whose oscillations produce a short, high frequency burst of gravitational radiation. Its dominant frequency content is tightly correlated with the radius of the neutron star, and its measurement can be used to constrain the supranuclear equation of state. In contrast, the merger of higher mass systems results in prompt gravitational collapse to a black hole. We have developed an algorithm which combines waveform reconstruction from a morphology-independent search for gravitational wave transients with Bayesian model selection, to discriminate between post-merger scenarios and accurately measure the dominant oscillation frequency. We demonstrate the efficacy of the method using a catalogue of simulated binary merger signals in data from LIGO and Virgo, and we discuss the prospects for this analysis in advanced ground-based gravitational wave detectors. From the waveforms considered in this work and assuming an optimally oriented source, we find that the post-merger neutron star signal may be detectable by this technique to $\sim 10\text{--}25$\,Mpc. We also find that we successfully discriminate between the post-merger scenarios with $\sim 95\%$ accuracy and determine the dominant oscillation frequency of surviving post-merger neutron stars to within $\sim 10$\,Hz, averaged over all detected signals. This leads to an uncertainty in the estimated radius of a non-rotating 1.6\,M$_{\odot}$ reference neutron star of $\sim 100\,$m.

We present a novel method for revealing the equation of state of high-density neutron star matter through gravitational waves emitted during the postmerger phase of a binary neutron star system. The method relies on a small number of detections of the peak frequency in the postmerger phase for binaries of different (relatively low) masses, in the most likely range of expected detections. From such observations, one can construct the derivative of the peak frequency versus the binary mass, in this mass range. Through a detailed study of binary neutron star mergers for a large sample of equations of state, we show that one can extrapolate the above information to the highest possible mass (the threshold mass for black hole formation in a binary neutron star merger). In turn, this allows for an empirical determination of the maximum mass of cold, nonrotating neutron stars to within 0.1 M_sun, while the corresponding radius is determined to within a few percent. Combining this with the determination of the radius of cold, nonrotating neutron stars of 1.6 M_sun (to within a few percent, as was demonstrated in Bauswein et al., PRD, 86, 063001, 2012), allows for a clear distinction of a particular candidate equation of state among a large set of other candidates. Our method is particularly appealing because it reveals simultaneously the moderate and very high-density parts of the equation of state, enabling the distinction of mass-radius relations even if they are similar at typical neutron star masses. Furthermore, our method also allows to deduce the maximum central energy density and maximum central rest-mass density of cold, nonrotating neutron stars with an accuracy of a few per cent.

The VERTEX code is employed for multi-dimensional neutrino-radiation hydrodynamics simulations of core-collapse supernova explosions from first principles. The code is considered state-of-the-art in supernova research and it has been used for modeling for more than a decade, resulting in numerous scientific publications. The computational performance of the code, which is currently deployed on several high-performance computing (HPC) systems up to the Tier-0 class (e.g. in the framework of the European PRACE initiative and the German GAUSS program), however, has so far not been extensively documented. This paper presents a high-level overview of the relevant algorithms and parallelization strategies and outlines the technical challenges and achievements encountered along the evolution of the code from the gigaflops scale with the first, serial simulations in 2000, up to almost petaflops capabilities, as demonstrated lately on the SuperMUC system of the Leibniz Supercomputing Centre (LRZ). In particular, we shall document the parallel scalability and computational efficiency of VERTEX at the large scale and on the major, contemporary HPC platforms. We will outline upcoming scientific requirements and discuss the resulting challenges for the future development and operation of the code.

The success of the neutrino mechanism of core-collapse supernovae relies on the supporting action of two hydrodynamic instabilities: neutrino-driven convection and the Standing Accretion Shock Instability (SASI). Depending on the structure of the stellar progenitor, each of these instabilities can dominate the evolution of the gain region prior to the onset of explosion, with implications for the ensuing asymmetries. Here we examine the flow dynamics in the neighborhood of explosion by means of parametric two-dimensional, time-dependent hydrodynamic simulations for which the linear stability properties are well understood. We find that systems for which the convection parameter is sub-critical (SASI-dominated) develop explosions once large-scale, high-entropy bubbles are able to survive for several SASI oscillation cycles. These long-lived structures are seeded by the SASI during shock expansions. Finite-amplitude initial perturbations do not alter this outcome qualitatively, though they can lead to significant differences in explosion times. Supercritical systems (convection-dominated) also explode by developing large-scale bubbles, though the formation of these structures is due to buoyant activity. Non-exploding systems achieve a quasi-steady state in which the time-averaged flow adjusts itself to be convectively sub-critical. We characterize the turbulent flow using a spherical Fourier-Bessel decomposition, identifying the relevant scalings and connecting temporal and spatial components. Finally, we verify the applicability of these principles on the general relativistic, radiation-hydrodynamic simulations of Mueller, Janka, & Heger (2012), and discuss implications for the three-dimensional case.

The likely outcome of a compact object merger event is a central black hole surrounded by a rapidly accreting torus of debris. This disk of debris is a rich source of element synthesis, the outcome of which is needed to predict electromagnetic counterparts of individual events and to understand the contribution of mergers to galactic chemical evolution. Here we study disk outflow nucleosynthesis in the context of a two-dimensional, time-dependent black hole-neutron star merger accretion disk model. We use two time snapshots from this model to examine the impact of the evolution of the neutrino fluxes from the disk on the element synthesis. While the neutrino fluxes from the early-time disk snapshot appear to favor neutron-rich outflows, by the late-time snapshot the situation is reversed. As a result we find copious production of Nickel-56 in the outflows.

Neutron star (NS) merger ejecta offer a viable site for the production of heavy r-process elements with nuclear mass numbers A > 140. The crucial role of fission recycling is responsible for the robustness of this site against many astrophysical uncertainties, but calculations sensitively depend on nuclear physics. In particular the fission fragment yields determine the creation of 110 < A < 170 nuclei. Here we apply a new scission-point model, called SPY, to derive the fission fragment distribution (FFD) of all relevant neutron-rich, fissioning nuclei. The model predicts a doubly asymmetric FFD in the abundant A ~ 278 mass region that is responsible for the final recycling of the fissioning material. Using ejecta conditions based on relativistic NS merger calculations we show that this specific FFD leads to a production of the A ~ 165 rare-earth peak that is nicely compatible with the abundance patterns in the Sun and metal-poor stars. This new finding further strengthens the case of NS mergers as possible dominant origin of r-nuclei with A > 140.

The first full-scale three-dimensional (3D) core-collapse supernova (SN) simulations with sophisticated neutrino transport show pronounced effects of the standing accretion shock instability (SASI) for two high-mass progenitors (20 and 27 M_sun). In a low-mass progenitor (11.2 M_sun), large-scale convection is the dominant nonradial hydrodynamic instability in the postshock accretion layer. The SASI-associated modulation of the neutrino signal (80 Hz in our two examples) will be clearly detectable in IceCube or the future Hyper-Kamiokande detector, depending on progenitor properties, distance, and observer location relative to the main SASI sloshing direction. The neutrino signal from the next galactic SN can therefore diagnose the nature of the hydrodynamic instability.

We perform hydrodynamical simulations of neutron-star mergers for a large sample of temperature-dependent, nuclear equations of state, and determine the threshold mass above which the merger remnant promptly collapses to form a black hole. We find that, depending on the equation of state, the threshold mass is larger than the maximum mass of a non-rotating star in isolation by between 30 and 70 per cent. Our simulations also show that the ratio between the threshold mass and maximum mass is tightly correlated with the compactness of the non-rotating maximum-mass configuration. We speculate on how this relation can be used to derive constraints on neutron-star properties from future observations.

We investigate the nucleosynthesis of the radionuclide 60Fe in electron-capture supernovae (ECSNe). The nucleosynthetic results are based on a self-consistent, two-dimensional simulation of an ECSN as well as models in which the densities are systematically increased by some factors (low-entropy models). 60Fe is found to be appreciably made in neutron-rich ejecta during the nuclear quasi-equilibrium phase with greater amounts being produced in the lower-entropy models. Our results, combining them with the yields of core-collapse supernovae (CCSNe) in the literature, suggest that ECSNe account for at least 4-30% of live 60Fe in the Milky Way. ECSNe co-produce neutron-rich isotopes, 48Ca, 50Ti, 54Cr, some light trans-iron elements, and possibly weak r-process elements including some radionuclides such as 93Zr, 99Tc, and 107Pd, whose association with 60Fe might have been imprinted in primitive meteorites or in the deep ocean crust on the Earth.

We investigate systematically the dynamical mass ejection, r-process nucleosynthesis, and properties of electromagnetic counterparts of neutron-star (NS) mergers in dependence on the uncertain properties of the nuclear equation of state (EoS) by employing 40 representative, microphysical high-density EoSs in relativistic, hydrodynamical simulations. The crucial parameter determining the ejecta mass is the radius R_1.35 of a 1.35 M_sun NS. NSs with smaller R_1.35 ("soft" EoS) eject systematically higher masses. These range from ~10^-3 M_sun to ~10^-2 M_sun for 1.35-1.35 M_sun binaries and from ~5*10^-3 M_sun to ~2*10^-2 M_sun for 1.2-1.5 M_sun systems (with kinetic energies between ~5*10^49 erg and 10^51 erg). Correspondingly, the bolometric peak luminosities of the optical transients of symmetric (asymmetric) mergers vary between 3*10^41 erg/s and 14*10^41 erg/s (9*10^41 erg/s and 14.5*10^41 erg/s) on timescales between ~2 h and ~12 h. If these signals with absolute bolometric magnitudes from -15.0 to -16.7 are measured, the tight correlation of their properties with those of the merging NSs might provide valuable constraints on the high-density EoS. The r-process nucleosynthesis exhibits a remarkable robustness independent of the EoS, producing a nearly solar abundance pattern above mass number 130. By the r-process content of the Galaxy and the average production per event the Galactic merger rate is limited to 4*10^-5/yr (4*10^-4/yr) for a soft (stiff) NS EoS, if NS mergers are the main source of heavy r-nuclei. The production ratio of radioactive 232Th to 238U attains a stable value of 1.64-1.67, which does not exclude NS mergers as potential sources of heavy r-material in the most metal-poor stars.

The relevance of the standing accretion shock instability (SASI) compared to neutrino-driven convection in three-dimensional (3D) supernova-core environments is still highly controversial. Studying a 27 Msun progenitor, we demonstrate, for the first time, that violent SASI activity can develop in 3D simulations with detailed neutrino transport despite the presence of convection. This result was obtained with the Prometheus-Vertex code with the same sophisticated neutrino treatment so far used only in 1D and 2D models. While buoyant plumes initially determine the nonradial mass motions in the postshock layer, bipolar shock sloshing with growing amplitude sets in during a phase of shock retraction and turns into a violent spiral mode whose growth is only quenched when the infall of the Si/SiO interface leads to strong shock expansion in response to a dramatic decrease of the mass accretion rate. In the phase of large-amplitude SASI sloshing and spiral motions, the postshock layer exhibits nonradial deformation dominated by the lowest-order spherical harmonics (l=1, m=0,-1,+1) in distinct contrast to the higher multipole structures associated with neutrino-driven convection. We find that the SASI amplitudes, shock asymmetry, and nonradial kinetic energy in 3D can exceed those of the corresponding 2D case during extended periods of the evolution. We also perform parametrized 3D simulations of a 25 Msun progenitor, using a simplified, gray neutrino transport scheme, an axis-free Yin-Yang grid, and different amplitudes of random seed perturbations. They confirm the importance of the SASI for another progenitor, its independence of the choice of spherical grid, and its preferred growth for fast accretion flows connected to small shock radii and compact proto-neutron stars as previously found in 2D setups.

Integrating trajectories of low-mass X-ray binaries containing black holes within the Galactic potential, Repetto, Davies & Sigurdsson recently showed that the large distances of some systems above the Galactic plane can only be explained if black holes receive appreciable natal kicks. Surprisingly, they found that the distribution of black hole kick velocities (rather than that of the momenta) should be similar to that of neutron stars. Here I argue that this result can be understood if neutron star and black hole kicks are a consequence of large-scale asymmetries created in the supernova ejecta by the explosion mechanism. The corresponding anisotropic gravitational attraction of the asymmetrically expelled matter does not only accelerate new-born neutron stars by the "gravitational tug-boat mechanism". It can also lead to delayed black-hole formation by asymmetric fallback of the slowest parts of the initial ejecta onto the transiently existing neutron star, in course of which the momentum of the black hole can grow with the fallback mass. Black hole kick velocities will therefore not be reduced by the ratio of neutron star to black hole mass as would be expected for kicks caused by anisotropic neutrino emission of the nascent neutron star.

We report that electron-capture supernovae (ECSNe), arising from collapsing oxygen-neon-magnesium cores, are a possible source of 48Ca, whose origin has remained a long-standing puzzle. Our two-dimensional, self-consistent explosion model of an ECSN predicts ejection of neutron-rich matter with electron fractions Ye = 0.40-0.42 and relatively low entropies, s = 13-15 kB per nucleon (kB is the Boltzmann constant). Post-processing nucleosynthesis calculations result in appreciable production of 48Ca in such neutron-rich and low-entropy matter during the quasi-nuclear equilibrium and subsequent freezeout phases. The amount of ejected 48Ca can account for that in the solar inventory when we consider possible uncertainties in the entropies or ejecta-mass distribution. ECSNe could thus be a site of 48Ca production in addition to a hypothetical, rare class of high-density Type Ia supernovae.

We present a detailed theoretical analysis of the gravitational-wave (GW) signal of the post-bounce evolution of core-collapse supernovae (SNe), employing for the first time relativistic, two-dimensional (2D) explosion models with multi-group, three-flavor neutrino transport based on the ray-by-ray-plus approximation. The waveforms reflect the accelerated mass motions associated with the characteristic evolutionary stages that were also identified in previous works: A quasi-periodic modulation by prompt postshock convection is followed by a phase of relative quiescence before growing amplitudes signal violent hydrodynamical activity due to convection and the standing accretion shock instability during the accretion period of the stalled shock. Finally, a high-frequency, low-amplitude variation from proto-neutron star (PNS) convection below the neutrinosphere appears superimposed on the low-frequency trend associated with the aspherical expansion of the SN shock after the onset of the explosion. Relativistic effects in combination with detailed neutrino transport are shown to be essential for quantitative predictions of the GW frequency evolution and energy spectrum, because they determine the structure of the PNS surface layer and its characteristic g-mode frequency. Burst-like high-frequency activity phases, correlated with sudden luminosity increase and spectral hardening of electron (anti-)neutrino emission for some 10ms, are discovered as new features after the onset of the explosion. They correspond to intermittent episodes of anisotropic accretion by the PNS in the case of fallback SNe. We find stronger signals for more massive progenitors with large accretion rates. The typical frequencies are higher for massive PNSs, though the time-integrated spectrum also strongly depends on the model dynamics.

We present 3D simulations of supernova (SN) explosions of nonrotating stars, triggered by the neutrino-heating mechanism with a suitable choice of the core-neutrino luminosity. Our results show that asymmetric mass ejection caused by hydrodynamic instabilities can accelerate the neutron star (NS) up to recoil velocities of more than 700 km/s by the "gravitational tug-boat mechanism", which is enough to explain most observed pulsar velocities. The associated NS spin periods are about 100 ms to 8 s without any correlation between spin and kick magnitudes or directions. This suggests that faster spins and a possible spin-kick alignment might require angular momentum in the progenitor core prior to collapse. Our simulations for the first time demonstrate a clear correlation between the size of the NS kick and anisotropic ejection of heavy elements created by explosive burning behind the shock. In the case of large NS kicks the explosion is significantly stronger opposite to the kick vector. Therefore the bulk of the Fe-group elements, in particular nickel, is ejected mostly in large clumps against the kick direction. This contrasts with the case of low recoil velocity, where the Ni-rich lumps are more isotropically distributed. Intermediate-mass nuclei heavier than Si (like Ca and Ti) also exhibit a significant enhancement in the hemisphere opposite to the direction of fast NS motion, while the distribution of C, O, and Ne is not affected, and that of Mg only marginally. Mapping the spatial distribution of the heavy elements in SN remnants with identified pulsar motion may offer an important diagnostic test of the kick mechanism. Different from kick scenarios based on anisotropic neutrino emission, our hydrodynamical acceleration model predicts enhanced ejection of Fe-group elements and of their nuclear precursors in the direction opposite to the NS recoil. (abridged)

To study the capabilities of supernova neutrino detectors, the instantaneous spectra are often represented by a quasi-thermal distribution of the form f(E) = E^alpha e^{-(alpha+1)E/E_{av}} where E_{av} is the average energy and alpha a numerical parameter. Based on a spherically symmetric supernova model with full Boltzmann neutrino transport we have, at a few representative post-bounce times, re-converged the models with vastly increased energy resolution to test the fit quality. For our examples, the spectra are well represented by such a fit in the sense that the counting rates for a broad range of target nuclei, sensitive to different parts of the spectrum, are reproduced very well. Therefore, the mean energy and root-mean-square energy of numerical spectra hold enough information to provide the correct alpha and to forecast the response of multi-channel supernova neutrino detection.

Core-collapse supernovae are among the most fascinating phenomena in astrophysics and provide a formidable challenge for theoretical investigation. They mark the spectacular end of the lives of massive stars and, in an explosive eruption, release as much energy as the sun produces during its whole life. A better understanding of the astrophysical role of supernovae as birth sites of neutron stars, black holes, and heavy chemical elements, and more reliable predictions of the observable signals from stellar death events are tightly linked to the solution of the long-standing puzzle how collapsing stars achieve to explode. In this article our current knowledge of the processes that contribute to the success of the explosion mechanism are concisely reviewed. After a short overview of the sequence of stages of stellar core-collapse events, the general properties of the progenitor-dependent neutrino emission will be briefly described. Applying sophisticated neutrino transport in axisymmetric (2D) simulations with general relativity as well as in simulations with an approximate treatment of relativistic effects, we could find successful neutrino-driven explosions for a growing set of progenitor stars. First results of three-dimensional (3D) models have been obtained, and magnetohydrodynamic simulations demonstrate that strong initial magnetic fields in the pre-collapse core can foster the onset of neutrino-powered supernova explosions even in nonrotating stars. These results are discussed in the context of the present controversy about the value of 2D simulations for exploring the supernova mechanism in realistic 3D environments, and they are interpreted against the background of the current disagreement on the question whether the standing accretion shock instability (SASI) or neutrino-driven convection is the crucial agency that supports the onset of the explosion.

The neutrino-driven explosion mechanism for core-collapse supernovae in its modern flavor relies on the additional support of hydrodynamical instabilities in achieving shock revival. Two possible candidates, convection and the so-called standing accretion shock instability (SASI), have been proposed for this role. In this paper, we discuss new successful simulations of supernova explosions that shed light on the relative importance of these two instabilities. While convection has so far been observed to grow first in self-consistent hydrodynamical models with multi-group neutrino transport, we here present the first such simulation in which the SASI grows faster while the development of convection is initially inhibited. We illustrate the features of this SASI-dominated regime using an explosion model of a 27 solar mass progenitor, which is contrasted with a convectively-dominated model of an 8.1 solar mass progenitor with subsolar metallicity, whose early post-bounce behavior is more in line with previous 11.2 and 15 solar mass explosion models. We analyze the conditions discriminating between the two different regimes, showing that a high mass-accretion rate and a short advection time-scale are conducive for strong SASI activity. We also briefly discuss some important factors for capturing the SASI-driven regime, such as general relativity, the progenitor structure, a nuclear equation of state leading to a compact proto-neutron star, and the neutrino treatment. Finally, we evaluate possible implications of our findings for 2D and 3D supernova simulations. Our results show that a better understanding of the SASI and convection in the non-linear regime is required.

We study supernova neutrino flux variations in the IceCube detector, using 3D models based on a simplified neutrino transport scheme. The hemispherically integrated neutrino emission shows significantly smaller variations compared with our previous study of 2D models, largely because of the reduced SASI activity in this set of 3D models which we interpret as a pessimistic extreme. For the studied cases, intrinsic flux variations up to about 100 Hz frequencies could still be detected in a supernova closer than about 2 kpc.

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

Neutron-star (NS) merger simulations are conducted for 38 representative microphysical descriptions of high-density matter in order to explore the equation-of-state dependence of the postmerger ring-down phase. The formation of a deformed, oscillating, differentially rotating very massive NS is the typical outcome of the coalescence of two stars with 1.35 $M_{\odot}$ for most candidate EoSs. The oscillations of this object imprint a pronounced peak in the gravitational-wave (GW) spectra, which is used to characterize the emission for a given model. The peak frequency of this postmerger GW signal correlates very well with the radii of nonrotating NSs, and thus allows to constrain the high-density EoS by a GW detection. In the case of 1.35-1.35 $M_{\odot}$ mergers the peak frequency scales particularly well with the radius of a NS with 1.6 $M_{\odot}$, where the maximum deviation from this correlation is only 60 meters for fully microphysical EoSs which are compatible with NS observations. Combined with the uncertainty in the determination of the peak frequency it appears likely that a GW detection can measure the radius of a 1.6 $M_{\odot}$ NS with an accuracy of about 100 to 200 meters. We also uncover relations of the peak frequency with the radii of nonrotating NSs with 1.35 $M_{\odot}$ or 1.8 $M_{\odot}$, with the radius or the central energy density of the maximum-mass Tolman-Oppenheimer-Volkoff configuration, and with the pressure or sound speed at a fiducial rest-mass density of about twice nuclear saturation density. Furthermore, it is found that a determination of the dominant postmerger GW frequency can provide an upper limit for the maximum mass of nonrotating NSs. The prospects for a detection of the postmerger GW signal and a determination of the dominant GW frequency are estimated to be in the range of 0.015 to 1.2 events per year with the upcoming Advanced LIGO detector.

We present the first two-dimensional general relativistic (GR) simulations of stellar core collapse and explosion with the CoCoNuT hydrodynamics code in combination with the VERTEX solver for energy-dependent, three-flavor neutrino transport, using the extended conformal flatness condition for approximating the spacetime metric and a ray-by-ray-plus ansatz to tackle the multi-dimensionality of the transport. For both of the investigated 11.2 and 15 solar mass progenitors we obtain successful, though seemingly marginal, neutrino-driven supernova explosions. This outcome and the time evolution of the models basically agree with results previously obtained with the PROMETHEUS hydro solver including an approximative treatment of relativistic effects by a modified Newtonian potential. However, GR models exhibit subtle differences in the neutrinospheric conditions compared to Newtonian and pseudo-Newtonian simulations. These differences lead to significantly higher luminosities and mean energies of the radiated electron neutrinos and antineutrinos and therefore to larger energy-deposition rates and heating efficiencies in the gain layer with favorable consequences for strong non-radial mass motions and ultimately for an explosion. Moreover, energy transfer to the stellar medium around the neutrinospheres through nucleon recoil in scattering reactions of heavy-lepton neutrinos also enhances the mentioned effects. Together with previous pseudo-Newtonian models the presented relativistic calculations suggest that the treatment of gravity and energy-exchanging neutrino interactions can make differences of even 50-100% in some quantities and is likely to contribute to a finally successful explosion mechanism on no minor level than hydrodynamical differences between different dimensions.

Following a simulation approach of recent publications we explore the viability of the neutrino-heating explosion mechanism in dependence on the spatial dimension. Our results disagree with previous findings. While we also observe that two-dimensional (2D) models explode for lower driving neutrino luminosity than spherically symmetric (1D) models, we do not find that explosions in 3D occur easier and earlier than in 2D. Moreover, we find that the average entropy of matter in the gain layer hardly depends on the dimension and thus is no good diagnostic quantity for the readiness to explode. Instead, mass, integrated entropy, total neutrino-heating rate, and nonradial kinetic energy in the gain layer are higher when models are closer to explosion. Coherent, large-scale mass motions as typically associated with the standing accretion-shock instability (SASI) are observed to be supportive for explosions because they drive strong shock expansion and thus enlarge the gain layer. While 2D models with better angular resolution explode clearly more easily, the opposite trend is seen in 3D. We interpret this as a consequence of the turbulent energy cascade, which transports energy from small to large spatial scales in 2D, thus fostering SASI activity. In contrast, the energy flow in 3D is in the opposite direction, feeding fragmentation and vortex motions on smaller scales and thus making the 3D evolution with finer grid resolution more similar to 1D. More favorable conditions for explosions in 3D may therefore be tightly linked to efficient growth of low-order SASI modes including nonaxisymmetric ones.

Neutrinos streaming from a supernova (SN) core occasionally scatter in the envelope, producing a small "neutrino halo" with a much broader angle distribution than the primary flux originating directly from the core. Cherry et al. (2012) have recently pointed out that, during the accretion phase, the halo actually dominates neutrino-neutrino refraction at distances exceeding some 100 km. However, the multiangle matter effect (which increases if the angle distribution is broader) still appears to suppress self-induced flavor conversion during the accretion phase.

Supernova theory, numerical and analytic, has made remarkable progress in the past decade. This progress was made possible by more sophisticated simulation tools, especially for neutrino transport, improved microphysics, and deeper insights into the role of hydrodynamic instabilities. Violent, large-scale nonradial mass motions are generic in supernova cores. The neutrino-heating mechanism, aided by nonradial flows, drives explosions, albeit low-energy ones, of ONeMg-core and some Fe-core progenitors. The characteristics of the neutrino emission from new-born neutron stars were revised, new features of the gravitational-wave signals were discovered, our notion of supernova nucleosynthesis was shattered, and our understanding of pulsar kicks and explosion asymmetries was significantly improved. But simulations also suggest that neutrino-powered explosions might not explain the most energetic supernovae and hypernovae, which seem to demand magnetorotational driving. Now that modeling is being advanced from two to three dimensions, more realism, new perspectives, and hopefully answers to long-standing questions are coming into reach.

The rise time of a Galactic supernova (SN) bar-nue lightcurve, observable at a high-statistics experiment such as the IceCube Cherenkov detector, can provide a diagnostic tool for the neutrino mass hierarchy at "large" 1-3 leptonic mixing angle theta_13. Thanks to the combination of matter suppression of collective effects at early postbounce times on one hand and the presence of the ordinary Mikheyev-Smirnov-Wolfenstein effect in the outer layers of the SN on the other hand, a sufficiently fast rise time on O(100) ms scale is indicative of an inverted mass hierarchy. We investigate results from an extensive set of stellar core-collapse simulations, providing a first exploration of the astrophysical robustness of these features. We find that for all the models analyzed (sharing the same weak interaction microphysics) the rise times for the same hierarchy are similar not only qualitatively, but also quantitatively, with the signals for the two classes of hierarchies significantly separated. We show via Monte Carlo simulations that the two cases should be distinguishable at IceCube for SNe at a typical Galactic distance 99% of the times. Finally, a preliminary survey seems to show that the faster rise time for inverted hierarchy as compared to normal hierarchy is a qualitatively robust feature predicted by several simulation groups. Since the viability of this signature ultimately depends on the quantitative assessment of theoretical/numerical uncertainties, our results motivate an extensive campaign of comparison of different code predictions at early accretion times with implementation of microphysics of comparable sophistication, including effects such like nucleon recoils in weak interactions.

We discuss the prospects for improved upper limits on neutrino masses that may be provided by a core-collapse supernova explosion in our galaxy, if it exhibits time variations in the neutrino emissions on the scale of a few milliseconds as suggested by recent two-dimensional simulations. Analyzing simulations of such neutrino emissions using the wavelet technique adopted in [1], we find that an upper limit m_nu ~ 0.14 eV could be established at the 95% confidence level if the time variations in emissions were to be preserved during neutrino propagation to the Earth.

We propose the liquid-scintillator detector LENA (Low Energy Neutrino Astronomy) as a next-generation neutrino observatory on the scale of 50 kt. The outstanding successes of the Borexino and KamLAND experiments demonstrate the large potential of liquid-scintillator detectors in low-energy neutrino physics. LENA's physics objectives comprise the observation of astrophysical and terrestrial neutrino sources as well as the investigation of neutrino oscillations. In the GeV energy range, the search for proton decay and long-baseline neutrino oscillation experiments complement the low-energy program. Based on the considerable expertise present in European and international research groups, the technical design is sufficiently mature to allow for an early start of detector realization.

Self-induced flavor conversions of supernova (SN) neutrinos can strongly modify the flavor dependent fluxes. We perform a linearized flavor stability analysis with accretion-phase matter profiles of a 15 M_sun spherically symmetric model and corresponding neutrino fluxes. We use realistic energy and angle distributions, the latter deviating strongly from quasi-isotropic emission, thus accounting for both multi-angle and multi-energy effects. For our matter and neutrino density profile we always find stable conditions: flavor conversions are limited to the usual MSW effect. In this case one may distinguish the neutrino mass hierarchy in a SN neutrino signal if the mixing angle theta_13 is as large as suggested by recent experiments.

Supernova explosions provide the most sensitive probes of neutrino propagation, such as the possibility that neutrino velocities might be affected by the foamy structure of space-time thought to be generated by quantum-gravitational (QG) effects. Recent two-dimensional simulations of the neutrino emissions from core-collapse supernovae suggest that they might exhibit variations in time on the scale of a few milliseconds. We analyze simulations of such neutrino emissions using a wavelet technique, and consider the limits that might be set on a linear or quadratic violation of Lorentz invariance in the group velocities of neutrinos of different energies, v/c = [1 \pm (E/M_{nuLV1})] or [1 \pm (E/M_{\nuLV2})^2], if variations on such short time scales were to be observed, where the mass scales M_{nuLVi} might appear in models of quantum gravity. We find prospective sensitivities to M_{nuLV1} ~ 2 X 10^{13} GeV and M_{nuLV2} ~ 10^6 GeV at the 95% confidence level, up to two orders of magnitude beyond estimates made using previous one-dimensional simulations of core-collapse supernovae. We also analyze the prospective sensitivities to scenarios in which the propagation times of neutrinos of fixed energies are subject to stochastic fluctuations.

We present results of general relativistic simulations of collapsing supermassive stars with and without rotation using the two-dimensional general relativistic numerical code Nada, which solves the Einstein equations written in the BSSN formalism and the general relativistic hydrodynamics equations with high resolution shock capturing schemes. These numerical simulations use an equation of state which includes effects of gas pressure, and in a tabulated form those associated with radiation and the electron-positron pairs. We also take into account the effect of thermonuclear energy released by hydrogen and helium burning. We find that objects with a mass of 5x10^{5} solar mass and an initial metallicity greater than Z_{CNO}~0.007 do explode if non-rotating, while the threshold metallicity for an explosion is reduced to Z_{CNO}~0.001 for objects uniformly rotating. The critical initial metallicity for a thermonuclear explosion increases for stars with mass ~10^{6} solar mass. For those stars that do not explode we follow the evolution beyond the phase of black hole formation. We compute the neutrino energy loss rates due to several processes that may be relevant during the gravitational collapse of these objects. The peak luminosities of neutrinos and antineutrinos of all flavors for models collapsing to a BH are ~10^{55} erg/s. The total radiated energy in neutrinos varies between ~10^{56} ergs for models collapsing to a BH, and ~10^{45}-10^{46} ergs for models exploding.

Motivated by recent hints for sterile neutrinos from the reactor anomaly, we study active-sterile conversions in a three-flavor scenario (2 active + 1 sterile families) for three different representative times during the neutrino-cooling evolution of the proto-neutron star born in an electron-capture supernova. In our "early model" (0.5 s post bounce), the nu_e-nu_s MSW effect driven by Delta m^2=2.35 eV^2 is dominated by ordinary matter and leads to a complete nu_e-nu_s swap with little or no trace of collective flavor oscillations. In our "intermediate" (2.9 s p.b.) and "late models" (6.5 s p.b.), neutrinos themselves significantly modify the nu_e-nu_s matter effect, and, in particular in the late model, nu-nu refraction strongly reduces the matter effect, largely suppressing the overall nu_e-nu_s MSW conversion. This phenomenon has not been reported in previous studies of active-sterile supernova neutrino oscillations. We always include the feedback effect on the electron fraction Y_e due to neutrino oscillations. In all examples, Y_e is reduced and therefore the presence of sterile neutrinos can affect the conditions for heavy-element formation in the supernova ejecta, even if probably not enabling the r-process in the investigated outflows of an electron-capture supernova. The impact of neutrino-neutrino refraction is strong but complicated, leaving open the possibility that with a more complete treatment, or for other supernova models, active-sterile neutrino oscillations could generate conditions suitable for the r-process.

We present results from the first generation of multi-dimensional general relativistic neutrino hydrodynamics simulations of core-collapse supernovae. A comparison with models computed using either the purely Newtonian approximation or the "effective gravitational potential" approach reveals appreciable quantitative differences in the heating conditions and the gravitational wave spectra. Our results underscore the important role of general relativity in the supernova problem (which appears to be on par with other important factors such as the dimensionality and the equation of state) both for our understanding of the explosion dynamics as well as for predictions of observable signatures.

The quest for the supernova explosion mechanism has been one of the outstanding challenges in computational astrophysics for several decades. Simulations have now progressed to a stage at which the solution appears close and neutrino and gravitational wave signals from self-consistent explosion models are becoming available. Here we focus one of the recent advances in supernova modeling, the inclusion of general relativity in multi-dimensional neutrino hydrodynamics simulations, and present the latest simulation results for an 11.2 and a 15 solar mass progenitor. We also mention 3D effects as another aspect in supernova physics awaiting further, more thorough investigation.

We examine r-process nucleosynthesis in the neutrino-driven wind from the thick accretion disk (or "torus") around a black hole. Such systems are expected as emnants of binary neutron star or neutron star -- black hole mergers. We consider a simplified, analytic, time-dependent evolution model of a 3M_sun central black hole surrounded by a neutrino emitting accretion torus with 90km radius, which serves as basis for computing spherically symmetric neutrino-driven wind solutions. We find that ejecta with modest entropies (~30 per nucleon in units of the Boltzmann constant) and moderate expansion timescales (~100ms) dominate in the mass outflow. The mass-integrated nucleosynthetic abundances are in good agreement with the solar system r-process abundance distribution if a minimal value of the electron fraction at the charged-particle freezeout, Ye,min~0.2, is achieved. In the case of Ye,min~0.3, the production of r-elements beyond A~130 does not reach to the third peak but could be still important for an explanation of the abundance signatures in r-process deficient stars in the early Galaxy. The total mass of the ejected r-process nuclei is estimated to be ~1x10^-3 M_sun. If our model was representative, this demands a Galactic event rate of ~2x10^-4 yr^-1 for black-hole-torus winds from merger remnants to be the dominant source of the r-process elements. Our result thus suggests that black-hole-torus winds from compact binary mergers have the potential to be a major but probably not the dominant production site of r-process elements.

Time-dependent and direction-dependent neutrino and gravitational-wave (GW) signatures are presented for a set of 3D hydrodynamic models of parametrized, neutrino-driven supernova explosions of non-rotating 15 and 20 solar mass stars. We employ an approximate treatment of neutrino transport. Due to the excision of the high-density core of the proto-neutron star and the use of an axis-free overset grid, the models can be followed from the post-bounce accretion phase for more than one second without imposing any symmetry restrictions. GW and neutrino emission exhibit the generic time-dependent features known from 2D models. Non-radial hydrodynamic mass motions in the accretion layer and their interaction with the outer layers of the proto-neutron star together with anisotropic neutrino emission give rise to a GW signal with an amplitude of ~5-20 cm and frequencies 100--500 Hz. The GW emission from mass motions reaches a maximum before the explosion sets in. Afterwards the GW signal exhibits a low-frequency modulation, in some cases describing a quasi-monotonic growth, associated with the non-spherical expansion of the explosion shock wave and the large-scale anisotropy of the escaping neutrino flow. Variations of the mass-quadrupole moment due to convective activity inside the nascent neutron star contribute a high-frequency component to the GW signal during the post-explosion phase. The GW signals exhibit strong variability between the two polarizations, different explosion simulations and different observer directions, and does not possess any template character. The neutrino emission properties show fluctuations over the neutron star surface on spatial and temporal scales that reflect the different types of non-spherical mass motions. The modulation amplitudes of the measurable neutrino luminosities and mean energies are significantly smaller than predicted by 2D simulations.

We demonstrate by a large set of merger simulations for symmetric binary neutron stars (NSs) that there is a tight correlation between the frequency peak of the postmerger gravitational-wave (GW) emission and the physical properties of the nuclear equation of state (EoS), e.g. expressed by the radius of the maximum-mass Tolman-Oppenheimer-Volkhoff configuration. Therefore, a single measurement of the peak frequency of the postmerger GW signal will constrain the NS EoS significantly. For plausible optimistic merger-rate estimates a corresponding detection with Advanced LIGO is likely to happen within an operation time of roughly a year.

Long Gamma-Ray Bursts (GRBs) are the most dramatic examples of massive stellar deaths, usually as- sociated with supernovae (Woosley et al. 2006). They release ultra-relativistic jets producing non-thermal emission through synchrotron radiation as they interact with the surrounding medium (Zhang et al. 2004). Here we report observations of the peculiar GRB 101225A (the "Christmas burst"). Its gamma-ray emission was exceptionally long and followed by a bright X-ray transient with a hot thermal component and an unusual optical counterpart. During the first 10 days, the optical emission evolved as an expanding, cooling blackbody after which an additional component, consistent with a faint supernova, emerged. We determine its distance to 1.6 Gpc by fitting the spectral-energy distribution and light curve of the optical emission with a GRB-supernova template. Deep optical observations may have revealed a faint, unresolved host galaxy. Our proposed progenitor is a helium star-neutron star merger that underwent a common envelope phase expelling its hydrogen envelope. The resulting explosion created a GRB-like jet which gets thermalized by interacting with the dense, previously ejected material and thus creating the observed black-body, until finally the emission from the supernova dominated. An alternative explanation is a minor body falling onto a neutron star in the Galaxy (Campana et al. 2011).

Although the rapid neutron-capture process, or r-process, is fundamentally important for explaining the origin of approximately half of the stable nuclei with A > 60, the astrophysical site of this process has not been identified yet. Here we study r-process nucleosynthesis in material that is dynamically ejected by tidal and pressure forces during the merging of binary neutron stars (NSs) and within milliseconds afterwards. For the first time we make use of relativistic hydrodynamical simulations of such events, defining consistently the conditions that determine the nucleosynthesis, i.e., neutron enrichment, entropy, early density evolution and thus expansion timescale, and ejecta mass. We find that 10^{-3}-10^{-2} solar masses are ejected, which is enough for mergers to be the main source of heavy (A > 140) galactic r-nuclei for merger rates of some 10^{-5} per year. While asymmetric mergers eject 2-3 times more mass than symmetric ones, the exact amount depends weakly on whether the NSs have radii of ~15 km for a "stiff" nuclear equation of state (EOS) or ~12 km for a "soft" EOS. R-process nucleosynthesis during the decompression becomes largely insensitive to the detailed conditions because of efficient fission recycling, producing a composition that closely follows the solar r-abundance distribution for nuclei with mass numbers A > 140. Estimating the light curve powered by the radioactive decay heating of r-process nuclei with an approximative model, we expect high emission in the B-V-R bands for 1-2 days with potentially observable longer duration in the case of asymmetric mergers because of the larger ejecta mass.

The rapid neutron-capture process, or r-process, is known to be fundamental for explaining the origin of approximately half of the A>60 stable nuclei observed in nature. In recent years nuclear astrophysicists have developed more and more sophisticated r-process models, by adding new astrophysical or nuclear physics ingredients to explain the solar system composition in a satisfactory way. Despite these efforts, the astrophysical site of the r-process remains unidentified. The composition of the neutron star outer crust material is investigated after the decompression that follows its possible ejection. The composition of the outer crust of a neutron star is estimated before and after decompression. Two different possible initial conditions are considered, namely an idealized crust composed of cold catalyzed matter and a crust initially in nuclear statistical equilibrium at temperatures around 10 GK. We show that in this second case before decompression and at temperatures typically corresponding to 8 GK, the Coulomb effect due to the high densities in the crust leads to an overall composition of the outer crust in neutron-rich nuclei with a mass distribution close to the solar system r-abundance distribution. Such distributions differ, however, from the solar one due to a systematic shift in the second peak to lower values. After decompression, the capture of the few neutrons per seed nucleus available in the hot outer crust leads to a final distribution of stable neutron-rich nuclei with a mass distribution of 80 < A < 140 nuclei in excellent agreement with the solar distribution, provided the outer crust is initially at temperatures around 8 GK and all layers of the outer crust are ejected. The decompression of the neutron star matter from the outer crust provides suitable conditions for a robust r-processing of the light species, i.e., r-nuclei with A < 140.

We study the excitation of nonaxisymmetric modes in the post-merger phase of binary compact object mergers and the associated gravitational wave emission. Our analysis is based on general-relativistic simulations, in the spatial conformal flatness approximation, using smoothed-particle-hydrodynamics for the evolution of matter, and we use a set of equal and unequal mass models, described by two nonzero-temperature hadronic equations of state and by one strange star equation of state. Through Fourier transforms of the evolution of matter variables, we can identify a number of oscillation modes, as well as several nonlinear components (combination frequencies). We focus on the dominant m=2 mode, which forms a triplet with two nonlinear components that are the result of coupling to the quasiradial mode. A corresponding triplet of frequencies is identified in the gravitational wave spectrum, when the individual masses of the compact objects are in the most likely range of 1.2 to 1.35 $M_\odot$. We can thus associate, through direct analysis of the dynamics of the fluid, a specific frequency peak in the gravitational wave spectrum with the nonlinear component resulting from the difference between the m=2 mode and the quasiradial mode. Once such observations become available, both the m=2 and quasiradial mode frequencies could be extracted, allowing for the application of gravitational-wave asteroseismology to the post-merger remnant and leading to tight constraints on the equation of state of high-density matter.

An 8.8 solar mass electron-capture supernova (SN) was simulated in spherical symmetry consistently from collapse through explosion to nearly complete deleptonization of the forming neutron star. The evolution time of about 9 s is short because of nucleon-nucleon correlations in the neutrino opacities. After a brief phase of accretion-enhanced luminosities (~200 ms), luminosity equipartition among all species becomes almost perfect and the spectra of electron antineutrinos and muon/tau antineutrinos very similar. We discuss consequences for the neutrino-driven wind as a nucleosynthesis site and for flavor oscillations of SN neutrinos.

Context. The influence of magnetic fields on stellar core collapse and explosion is not well explored. It depends on the possibility to amplify the pre-collapse fields. Without rotation this can happen by compression, convection, the standing accretion shock instability, and the accumulation and growth of Alfven waves in the accretion flows. Aims. We investigate such amplification mechanisms of the magnetic field during the collapse and post-bounce evolution of the core of a non-rotating 15 solar mass star with varied initial field strengths, taking into account the microphysical equation of state and neutrino physics crucial for supernova cores. Methods. We perform simulations of ideal MHD with neutrino transport in axisymmetry. The transport of electron neutrinos and antineutrinos is treated with a new scheme that solves the energy-dependent set of radiation energy and momentum equations in 2d by using an analytic closure relation. Results. The magnetic field undergoes amplification by turbulent flows. We find indications for amplification by interacting waves in accretion streams. The fields can reach up to equipartition with the velocity field. Very high magnetic field strengths require very strong pre-collapse fields and are able to shape the post-bounce flow, leading to a pattern dominated by low-order multipoles. Such models are closest to a successful explosion. Conclusions. Magnetic fields can build up to interesting strengths even in non-rotating collapsing stellar cores. Starting with fields in the pre-collapse core as predicted by present stellar evolution models, typical neutron star fields emerge, whereas stronger progenitor fields lead to fields of magnetar strength. Only in the latter case the fields have dynamical effects on the flows in the supernova core. However, in none of our 2D simulations, we find an explosion until 500 ms post-bounce.

We examine how the uncertainties involved in supernova dynamics as well as in nuclear data inputs affect the nup-process in the neutrino-driven winds. For the supernova dynamics, we find that the wind termination by the preceding dense ejecta shell, as well as the electron fraction (Y_{e, 3}; at 3 10^9 K) play a crucial role. A wind termination within the temperature range of (1.5-3) 10^9 K greatly enhances the efficiency of the nu p-process. This implies that the early wind phase, when the innermost layer of the preceding supernova ejecta is still 200-1000 km from the center, is most relevant to the nup-process. The outflows with Y_{e, 3} = 0.52-0.60 result in the production of the p-nuclei up to A=108 with interesting amounts. Furthermore, the p-nuclei up to A=152 can be produced if Y_{e, 3} = 0.65 is achieved. For the nuclear data inputs, we test the sensitivity to the rates relevant to the breakout from the pp-chain region (A < 12), to the (n, p) rates on heavy nuclei, and to the nuclear masses along the nup-process pathway. We find that a small variation of the rates of triple-alpha and of the (n, p) reaction on 56Ni leads to a substantial change in the p-nuclei production. We also find that 96Pd (N=50) on the nup-process path plays a role as a second seed nucleus for the production of heavier p-nuclei. The uncertainty in the nuclear mass of 82Zr can lead to a factor of two reduction in the abundance of the p-isotope 84Sr.

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

We examine electron-capture supernovae (ECSNe) as sources of elements heavier than iron in the solar system and in Galactic halo stars. Nucleosynthesis calculations are performed on the basis of thermodynamic histories of mass elements from a fully self-consistent, two-dimensional (2D) hydrodynamic explosion model of an ECSN. We find that neutron-rich convective lumps with an electron fraction down to Ye,min=0.40, which are absent in the one-dimensional (1D) counterpart, allow for interesting production of elements between the iron group and N=50 nuclei (from Zn to Zr, with little Ga) in nuclear (quasi-)equilibrium. Our models yield very good agreement with the Ge, Sr, Y, and Zr abundances of r-process deficient Galactic halo stars and constrain the occurrence of ECSNe to ~4% of all stellar core-collapse events. If tiny amounts of additional material with slightly lower Ye,min down to ~0.30-0.35 were also ejected - which presently cannot be excluded because of the limitations of resolution and two-dimensionality of the model -, a weak r-process can yield elements beyond N=50 up to Pd, Ag, and Cd as observed in the r-process deficient stars.

Using three-dimensional (3D) simulations of neutrino-powered supernova explosions we show that the hydrodynamical kick scenario proposed by Scheck et al. on the basis of two-dimensional (2D) models can yield large neutron star (NS) recoil velocities also in 3D. Although the shock stays relatively spherical, standing accretion-shock and convective instabilities lead to a globally asymmetric mass and energy distribution in the postshock layer. An anisotropic momentum distribution of the ejecta is built up only after the explosion sets in. Total momentum conservation implies the acceleration of the NS on a timescale of 1-3 seconds, mediated mainly by long-lasting, asymmetric accretion downdrafts and the anisotropic gravitational pull of large inhomogeneities in the ejecta. In a limited set of 15 solar-mass models with an explosion energy of about 10^51 erg this stochastic mechanism is found to produce kicks from <100 km/s to >500 km/s, and >1000 km/s seem possible. Strong rotational flows around the accreting NS do not develop in our collapsing, non-rotating progenitors. The NS spins therefore remain low with estimated periods of about 500-1000 ms and no alignment with the kicks.

We perform three-dimensional relativistic hydrodynamical calculations of neutron star mergers to assess the reliability of an approximate treatment of thermal effects in such simulations by combining an ideal-gas component with zero-temperature, micro-physical equations of state. To this end we compare the results of simulations that make this approximation to the outcome of models with a consistent treatment of thermal effects in the equation of state. In particular we focus on the implications for observable consequences of merger events like the gravitational-wave signal. It is found that the characteristic gravitational-wave oscillation frequencies of the post-merger remnant differ by about 50 to 250 Hz (corresponding to frequency shifts of 2 to 8 per cent) depending on the equation of state and the choice of the characteristic index of the ideal-gas component. In addition, the delay time to black hole collapse of the merger remnant as well as the amount of matter remaining outside the black hole after its formation are sensitive to the description of thermal effects.

In the delayed explosion scenario of core-collapse supernovae (SNe), the accretion phase shows pronounced convective overturns and a low-multipole hydrodynamic instability, the standing accretion shock instability (SASI). These effects imprint detectable fast time variations on the emerging neutrino flux. Among existing detectors, IceCube is best suited to this task, providing an event rate of ~1000 events per ms during the accretion phase for a fiducial SN distance of 10 kpc, comparable to what could be achieved with a megaton water Cherenkov detector. If the SASI activity lasts for several hundred ms, a Fourier component with an amplitude of 1% of the average signal clearly sticks out from the shot noise. We analyze in detail the output of axially symmetric hydrodynamical simulations that predict much larger amplitudes up to frequencies of a few hundred Hz. If these models are roughly representative for realistic SNe, fast time variations of the neutrino signal are easily detectable in IceCube or future megaton-class instruments. We also discuss the information that could be deduced from such a measurement about the physics in the SN core and the explosion mechanism of the SN.

All the current r-process scenarios relevant to core-collapse supernovae are facing severe difficulties. In particular, recent core-collapse simulations with neutrino transport show no sign of a neutron-rich wind from the proto-neutron star. In this paper, we discuss nucleosynthesis of the r-process in an alternative astrophysical site, "black hole winds", which are the neutrino-driven outflow from the accretion torus around a black hole. This condition is assumed to be realized in double neutron star mergers, neutron star - black hole mergers, or hypernovae.

We present a new general relativistic (GR) code for hydrodynamic supernova simulations with neutrino transport in spherical and azimuthal symmetry (1D/2D). The code is a combination of the CoCoNuT hydro module, which is a Riemann-solver based, high-resolution shock-capturing method, and the three-flavor, energy-dependent neutrino transport scheme VERTEX. VERTEX integrates the neutrino moment equations with a variable Eddington factor closure computed from a model Boltzmann equation and uses the ray-by-ray plus approximation in 2D, assuming the neutrino distribution to be axially symmetric around the radial direction, and thus the neutrino flux to be radial. Our spacetime treatment employs the ADM 3+1 formalism with the conformal flatness condition for the spatial three-metric. This approach is exact in 1D and has been shown to yield very accurate results also for rotational stellar collapse. We introduce new formulations of the energy equation to improve total energy conservation in relativistic and Newtonian hydro simulations with Eulerian finite-volume codes. Moreover, a modified version of the VERTEX scheme is developed that simultaneously conserves energy and lepton number with better accuracy and higher numerical stability. To verify our code, we conduct a series of tests, including a detailed comparison with published 1D results for stellar core collapse. Long-time simulations of proto-neutron star cooling over several seconds both demonstrate the robustness of the new CoCoNuT-VERTEX code and show the approximate treatment of GR effects by means of an effective gravitational potential as in PROMETHEUS-VERTEX to be remarkably accurate in 1D. (abridged)

We present the first three-dimensional (3D) simulations of the large-scale mixing that takes place in the shock-heated stellar layers ejected in the explosion of a 15.5 solar-mass blue supergiant star. The outgoing supernova shock is followed from its launch by neutrino heating until it breaks out from the stellar surface more than two hours after the core collapse. Violent convective overturn in the post-shock layer causes the explosion to start with significant asphericity, which triggers the growth of Rayleigh-Taylor (RT) instabilities at the composition interfaces of the exploding star. Deep inward mixing of hydrogen (H) is found as well as fast-moving, metal-rich clumps penetrating with high velocities far into the H-envelope of the star as observed, e.g., in the case of SN 1987A. Also individual clumps containing a sizeable fraction of the ejected iron-group elements (up to several 0.001 solar masses) are obtained in some models. The metal core of the progenitor is partially turned over with Ni-dominated fingers overtaking oxygen-rich bullets and both Ni and O moving well ahead of the material from the carbon layer. Comparing with corresponding 2D (axially symmetric) calculations, we determine the growth of the RT fingers to be faster, the deceleration of the dense metal-carrying clumps in the He and H layers to be reduced, the asymptotic clump velocities in the H-shell to be higher (up to ~4500 km/s for the considered progenitor and an explosion energy of 10^{51} ergs, instead of <2000 km/s in 2D), and the outward radial mixing of heavy elements and inward mixing of hydrogen to be more efficient in 3D than in 2D. We present a simple argument that explains these results as a consequence of the different action of drag forces on moving objects in the two geometries. (abridged)

Context: Mergers of neutron stars (NS) and black holes (BH) are among the strongest sources of gravitational waves and are potential central engines for short gamma-ray bursts. Aims: We aim to compare the general relativistic (GR) results by other groups with Newtonian calculations of models with equivalent parameters. We vary the mass ratios between NS and BH and the compactness of the NS. The mass of the NS is 1.4 M_sol. We compare the dynamics in the parameter-space regions where the NS is expected to reach the innermost stable circular orbit (ISCO) before being tidally disrupted (mass shedding, MS) and vice versa. Methods: The hydrodynamics is evolved by a Newtonian PPM scheme with four levels of nested grids. We use a polytropic EoS (Gamma=2), as was done in the GR simulations. However, instead of full GR we use a Newtonian potential supplemented by a Paczynski-Wiita-Artemova potential for the BH, both disregarding and including rotation of the BH. Results: If the NS is compact (C=0.18) it is accreted by the BH more quickly, and only a small amount of mass remains outside the BH. If the mass ratio is small (Q=2 or 3) or the NS is less compact (C=0.16 or less) the NS is tidally torn apart before being accreted. Although most of the mass is absorbed by the BH, some 0.1 M_sol remain in a tidal arm. For small mass ratios the tidal arm can wrap around the BH to form a thick disk. When including the effects of BH spin-up or spin-down by the accreted matter, more mass remains in the surroundings (0.2-0.3 M_sol). Conclusions: Although details and quantitative results differ, the general trends of our Newtonian calculations are similar to the GR calculations. A clear delimiting line that separates ISCO from the MS cases is not found. Inclusion of BH rotation as well as sufficient numerical resolution are extremely important.

We perform three-dimensional relativistic hydrodynamical simulations of the coalescence of strange stars (SSs) and explore the possibility to decide on the strange matter hypothesis by means of gravitational-wave (GW) measurements. Selfbinding of strange quark matter (SQM) and the generally more compact stars yield features that clearly distinguish SS from neutron star (NS) mergers, e.g. hampering tidal disruption during the plunge of quark stars. Furthermore, instead of forming dilute halos around the remnant as in the case of NS mergers, the coalescence of SSs results in a differentially rotating hypermassive object with a sharp surface layer surrounded by a geometrically thin, clumpy high-density SQM disk. We also investigate the importance of including non-zero temperature equations of state (EoSs) in NS and SS merger simulations. In both cases we find a crucial sensitivity of the dynamics and outcome of the coalescence to thermal effects, which, e.g., determine the outer remnant structure and the delay time of the dense remnant core to black hole collapse. For comparing and classifying the GW signals, we use a number of characteristic quantities like the maximum frequency during inspiral or the dominant frequency of oscillations of the postmerger remnant. In general, these frequencies are higher for SS mergers. If not, additional features of the GW luminosity spectrum may help to discriminate coalescence events of the different types. Future GW measurements may thus help to decide on the existence of SQM stars. (abridged)

(Abridged.) The accretion-induced collapse (AIC) of a white dwarf (WD) may lead to the formation of a protoneutron star and a collapse-driven supernova explosion. This process represents a path alternative to thermonuclear disruption of accreting white dwarfs in Type Ia supernovae. Neutrino and gravitational-wave (GW) observations may provide crucial information necessary to reveal a potential AIC. Motivated by the need for systematic predictions of the GW signature of AIC, we present results from an extensive set of general-relativistic AIC simulations using a microphysical finite-temperature equation of state and an approximate treatment of deleptonization during collapse. Investigating a set of 114 progenitor models in rotational equilibrium, with a wide range of rotational configurations, temperatures and central densities, we extend previous Newtonian studies and find that the GW signal has a generic shape akin to what is known as a "Type III" signal in the literature. We discuss the detectability of the emitted GWs, showing that the signal-to-noise ratio for current or next-generation interferometer detectors could be high enough to detect such events in our Galaxy. Some of our AIC models form massive quasi-Keplerian accretion disks after bounce. In rapidly differentially rotating models, the disk mass can be as large as ~0.8-Msun. Slowly and/or uniformly rotating models produce much smaller disks. Finally, we find that the postbounce cores of rapidly spinning white dwarfs can reach sufficiently rapid rotation to develop a nonaxisymmetric rotational instability.

We determine the Galactic production rate of strangelets as a canonical input to calculations of the measurable cosmic ray flux of strangelets by performing simulations of strange star mergers and combining the results with recent estimates of stellar binary populations. We find that the flux depends sensitively on the bag constant of the MIT bag model of QCD and disappears for high values of the bag constant and thus more compact strange stars. In the latter case strange stars could coexist with ordinary neutron stars as they are not converted by the capture of cosmic ray strangelets. An unambiguous detection of an ordinary neutron star would then not rule out the strange matter hypothesis.

We present 2D hydrodynamic simulations of the long-time accretion phase of a 15 solar mass star after core bounce and before the launch of a supernova explosion. Our simulations are performed with the Prometheus-Vertex code, employing multi-flavor, energy-dependent neutrino transport and an effective relativistic gravitational potential. Testing the influence of a stiff and a soft equation of state for hot neutron star matter, we find that the non-radial mass motions in the supernova core due to the standing accretion shock instability (SASI) and convection impose a time variability on the neutrino and gravitational-wave signals. These variations have larger amplitudes as well as higher frequencies in the case of a more compact nascent neutron star. After the prompt shock-breakout burst of electron neutrinos, a more compact accreting remnant radiates neutrinos with higher luminosities and larger mean energies. The observable neutrino emission in the direction of SASI shock oscillations exhibits a modulation of several 10% in the luminosities and ~1 MeV in the mean energies with most power at typical SASI frequencies of 20-100 Hz. At times later than 50-100 ms after bounce the gravitational-wave amplitude is dominated by the growing low-frequency (<200 Hz) signal associated with anisotropic neutrino emission. A high-frequency wave signal is caused by nonradial gas flows in the outer neutron star layers, which are stirred by anisotropic accretion from the SASI and convective regions. The gravitational-wave power then peaks at about 300-800 Hz with distinctively higher spectral frequencies originating from the more compact and more rapidly contracting neutron star. The detectability of the SASI effects in the neutrino and gravitational-wave signals is briefly discussed. (abridged)

We examine nucleosynthesis in the electron capture supernovae of progenitor AGB stars with an O-Ne-Mg core (with the initial stellar mass of 8.8 M_\odot). Thermodynamic trajectories for the first 810 ms after core bounce are taken from a recent state-of-the-art hydrodynamic simulation. The presented nucleosynthesis results are characterized by a number of distinct features that are not shared with those of other supernovae from the collapse of stars with iron core (with initial stellar masses of more than 10 M_\odot). First is the small amount of 56Ni (= 0.002-0.004 M_\odot) in the ejecta, which can be an explanation for observed properties of faint supernovae such as SNe 2008S and 1997D. In addition, the large Ni/Fe ratio is in reasonable agreement with the spectroscopic result of the Crab nebula (the relic of SN 1054). Second is the large production of 64Zn, 70Ge, light p-nuclei (74Se, 78Kr, 84Sr, and 92Mo), and in particular, 90Zr, which originates from the low Y_e (= 0.46-0.49, the number of electrons per nucleon) ejecta. We find, however, that only a 1-2% increase of the minimum Y_e moderates the overproduction of 90Zr. In contrast, the production of 64Zn is fairly robust against a small variation of Y_e. This provides the upper limit of the occurrence of this type of events to be about 30% of all core-collapse supernovae.

We study here the formation of heavy r-process nuclei in the high-entropy environment of rapidly expanding neutrino-driven winds from compact objects. In particular, we explore the sensitivity of the element creation in the A>130 region to the low-temperature behavior of the outflows. For this purpose we employ a simplified model of the dynamics and thermodynamical evolution for radiation dominated, adiabatic outflows. It consists of a first stage of fast, exponential cooling, followed by a second phase of slower evolution, either assuming constant density and temperature or a power-law decay of these quantities. These cases are supposed to capture the most relevant effects of a strong deceleration or decreasing acceleration of the transsonic outflows, respectively, e.g. in a wind termination shock caused by the collision with the slower, preceding supernova ejecta. We find that not only the transition temperature between the two expansion phases can make a big difference in the formation of the platinum peak, but also the detailed cooling law during the later phase. Unless the transition temperature and corresponding (free neutron) density become too small (T < 2*10^8 K), a lower temperature or faster temperature decline during this phase allow for a stronger appearance of the third abundance peak. Since the nuclear photodisintegration rates between ~2*10^8 K and ~10^9 K are more sensitive to the temperature than the n-capture rates are to the free neutron density, a faster cooling in this temperature regime shifts the r-process path closer to the n-drip line. With low (gamma,n)- but high beta-decay rates, the r-processing then does not proceed through a (gamma,n)-(n,gamma) equilibrium but through a quasi-equilibrium of (n,gamma)-reactions and beta-decays, as recently also pointed out by Wanajo.

We present results of 2D hydrodynamic simulations of stellar core collapse, which confirm that the neutrino-heating mechanism remains viable for the explosion of a wider mass range of supernova progenitors with iron cores. We used an energy-dependent treatment of the neutrino transport based on the "ray-by-ray plus" approximation, in which the number, energy, and momentum equations are closed with a variable Eddington factor obtained by iteratively solving a model Boltzmann equation. We focus on the evolution of a 15 Msun progenitor and show that shock revival and the explosion are initiated at about 600 ms post bounce, powered by neutrino energy deposition. Similar to previous findings for an 11.2 Msun star, but significantly later, the onset of the explosion is fostered by the standing accretion shock instability (SASI). This instability exhibits highest growth rates for the dipole and quadrupole modes, which lead to large-amplitude bipolar shock oscillations and push the shock to larger radii, thus increasing the time accreted matter is exposed to neutrino heating in the gain layer. Therefore also convective overturn behind the shock is strengthened. A "soft" nuclear equation of state that causes a rapid contraction and a smaller radius of the forming neutron star and thus a fast release of gravitational binding energy, seems to be more favorable for an explosion. Rotation has the opposite effect because it leads to a more extended and cooler neutron star and thus lower neutrino luminosities and mean energies and overall less neutrino heating. Neutron star g-mode oscillations and the acoustic mechanism play no important role in our simulations. (abridged)

We present results from detailed general relativistic simulations of stellar core collapse to a proto-neutron star, using two different microphysical nonzero-temperature nuclear equations of state as well as an approximate description of deleptonization during the collapse phase. Investigating a wide variety of rotation rates and profiles as well as masses of the progenitor stars and both equations of state, we confirm in this very general setup the recent finding that a generic gravitational wave burst signal is associated with core bounce, already known as type I in the literature. The previously suggested type II (or "multiple-bounce") waveform morphology does not occur. Despite this reduction to a single waveform type, we demonstrate that it is still possible to constrain the progenitor and postbounce rotation based on a combination of the maximum signal amplitude and the peak frequency of the emitted gravitational wave burst. Our models include to sufficient accuracy the currently known necessary physics for the collapse and bounce phase of core-collapse supernovae, yielding accurate and reliable gravitational wave signal templates for gravitational wave data analysis. In addition, we assess the possiblity of nonaxisymmetric instabilities in rotating nascent proto-neutron stars. We find strong evidence that in an iron core-collapse event the postbounce core cannot reach sufficiently rapid rotation to become subject to a classical bar-mode instability. However, many of our postbounce core models exhibit sufficiently rapid and differential rotation to become subject to the recently discovered dynamical instability at low rotation rates.

We consider hot accretion disk outflows from black hole - neutron star mergers in the context of the nucleosynthesis they produce. We begin with a three dimensional numerical model of a black hole - neutron star merger and calculate the neutrino and antineutrino fluxes emitted from the resulting accretion disk. We then follow the element synthesis in material outflowing the disk along parameterized trajectories. We find that at least a weak r-process is produced, and in some cases a main r-process as well. The neutron-rich conditions required for this production of r-process nuclei stem directly from the interactions of the neutrinos emitted by the disk with the free neutrons and protons in the outflow.

We discuss the flavor conversion of neutrinos from core collapse supernovae that have oxygen-neon-magnesium (ONeMg) cores. Using the numerically calculated evolution of the star up to 650 ms post bounce, we find that, for the normal mass hierarchy, the electron neutrino flux in a detector shows signatures of two typical features of an ONeMg-core supernova: a sharp step in the density profile at the base of the He shell and a faster shock wave propagation compared to iron core supernovae. Before the shock hits the density step (t ~ 150 ms), the survival probability of electron neutrinos is about 0.68, in contrast to values of 0.32 or less for an iron core supernova. The passage of the shock through the step and its subsequent propagation cause a decrease of the survival probability and a decrease of the amplitude of oscillations in the Earth, reflecting the transition to a more adiabatic propagation inside the star. These changes affect the lower energy neutrinos first; they are faster and more sizable for larger theta_13. They are unique of ONeMg-core supernovae, and give the possibility to test the speed of the shock wave. The time modulation of the Earth effect and its negative sign at the neutronization peak are the most robust signatures in a detector.

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.

It has been recently proposed that the shocked surface layers of exploding O-Ne-Mg cores provide the conditions for r-process nucleosynthesis, because their rapid expansion and high entropies enable heavy r-process isotopes to form even in an environment with very low initial neutron excess of the matter. We show here that the most sophisticated available hydrodynamic simulations (in spherical and axial symmetry) do not support this new r-process scenario because they fail to provide the necessary conditions of temperature, entropy, and expansion timescale by significant factors. This suggests that, either the formation of r-process elements works differently than suggested by Ning et al. (2007, NQM07), or that some essential core properties with influence on the explosion dynamics might be different from those predicted by Nomoto's progenitor model.

Nucleosynthesis in early neutrino winds is investigated. Presented is a brief overview of two recent problems of supernova nucleosynthesis. In the first part we investigate the effect of nuclear parameters on the synthesis of Mo92 and Mo94. Based on recent experimental results, we find that the proton rich winds of the model investigated here can not be the source of solar Mo92 and Mo94. In the second part we investigate the nucleosynthesis from neutron rich bubbles and show that they do not contribute to the overall nucleosynthesis.

We have studied detailed nucleosynthesis in the shocked surface layers of an Oxygen-Neon-Magnesium core collapse supernova with an eye to determining if the conditions are suitable for r process nucleosynthesis. We find no such conditions in an unmodified model, but do find overproduction of N=50 nuclei (previously seen in early neutron-rich neutrino winds) in amounts that, if ejected, would pose serious problems for galactic chemical evolution.

We report here on recent progress in understanding the birth conditions of neutron stars and the way how supernovae explode. More sophisticated numerical models have led to the discovery of new phenomena in the supernova core, for example a generic hydrodynamic instability of the stagnant supernova shock against low-mode nonradial deformation and the excitation of gravity-wave activity in the surface and core of the nascent neutron star. Both can have supportive or decisive influence on the inauguration of the explosion, the former by improving the conditions for energy deposition by neutrino heating in the postshock gas, the latter by supplying the developing blast with a flux of acoustic power that adds to the energy transfer by neutrinos. While recent two-dimensional models suggest that the neutrino-driven mechanism may be viable for stars from about 8 solar masses to at least 15 solar masses, acoustic energy input has been advocated as an alternative if neutrino heating fails. Magnetohydrodynamic effects constitute another way to trigger explosions in connection with the collapse of sufficiently rapidly rotating stellar cores, perhaps linked to the birth of magnetars. The global explosion asymmetries seen in the recent simulations offer an explanation of even the highest measured kick velocities of young neutron stars.

Recent nucleosynthesis calculations of Type II supernovae using advanced neutrino transport determine that the early neutrino winds are proton-rich. However, a fraction of the ejecta emitted at the same time is composed of neutron-rich pockets. In this paper we calculate the nucleosynthesis contribution from the neutron-rich pockets in the hot convective bubbles of a core-collapse supernova and show that they do not contribute significantly to the total nucleosynthesis.

By 2D hydrodynamic simulations including a detailed equation of state and neutrino transport, we investigate the interplay between different non-radial hydrodynamic instabilities that play a role during the postbounce accretion phase of collapsing stellar cores. The convective mode of instability, which is driven by negative entropy gradients caused by neutrino heating or by time variations of the shock strength, can be identified clearly by the development of typical Rayleigh-Taylor mushrooms. However, in cases where the gas in the postshock region is rapidly advected towards the gain radius, the growth of such a buoyancy instability can be suppressed. In such a situation the shocked flow nevertheless can develop non-radial asymmetry with an oscillatory growth of the amplitude. This phenomenon has been termed ``standing accretion shock instability'' (SASI). It is shown here that the SASI oscillations can trigger convective instability and like the latter they lead to an increase of the average shock radius and of the mass in the gain layer. Both hydrodynamic instabilities in combination stretch the advection time of matter through the neutrino-heating layer and thus enhance the neutrino energy deposition in support of the neutrino-driven explosion mechanism. A rapidly contracting and more compact nascent NS turns out to be favorable for explosions, because the accretion luminosity and neutrino heating are larger and the growth rate of the SASI is higher. Moreover, we show that the oscillation period of the SASI and a variety of other features in our simulations agree with estimates for the advective-acoustic cycle (AAC), in which perturbations are carried by the accretion flow from the shock to the neutron star and pressure waves close an amplifying global feedback loop. (abridged)

We analyze the gravitational wave (GW) emission from our recently published set of relativistic neutron star (NS) merger simulations and determine characteristic signal features that allow one to link GW measurements to the properties of the merging binary stars. We find that the distinct peak in the GW energy spectrum that is associated with the formation of a hypermassive merger remnant has a frequency that depends strongly on the properties of the nuclear equation of state (EoS) and on the total mass of the binary system, whereas the mass ratio and the NS spins have a weak influence. If the total mass can be determined from the inspiral chirp signal, the peak frequency of the postmerger signal is a sensitive indicator of the EoS.

The neutrino-heating mechanism remains a viable possibility for the cause of the explosion in a wide mass range of supernova progenitors. This is demonstrated by recent two-dimensional hydrodynamic simulations with detailed, energy-dependent neutrino transport. Neutrino-driven explosions were not only found for stars in the range of 8-10 solar masses with ONeMg cores and in case of the iron core collapse of a progenitor with 11 solar masses, but also for a ``typical'' progenitor model of 15 solar masses. For such more massive stars, however, the explosion occurs significantly later than so far thought, and is crucially supported by large-amplitude bipolar oscillations due to the nonradial standing accretion shock instability (SASI), whose low (dipole and quadrupole) modes can develop large growth rates in conditions where convective instability is damped or even suppressed. The dominance of low-mode deformation at the time of shock revival has been recognized as a possible explanation of large pulsar kicks and of large-scale mixing phenomena observed in supernovae like SN 1987A.

We perform general relativistic simulations of stellar core collapse to a proto-neutron star, using a microphysical equation of state as well as an approximate description of deleptonization. We show that for a wide variety of rotation rates and profiles the gravitational wave burst signals from the core bounce are of a generic type, known as Type I in the literature. In our systematic study, using both general relativity and Newtonian gravity, we identify and individually quantify the micro- and macrophysical mechanisms leading to this result, i.e. the effects of rotation, the equation of state, and deleptonization. Such a generic type of signal templates will likely facilitate a more efficient search in current and future gravitational wave detectors of both interferometric and resonant type.

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

We present detailed results from performing general relativistic (GR) simulations of stellar core collapse to a proto-neutron star, using a microphysical equation of state (EoS) as well as an approximate description of deleptonization during the collapse phase. We show that for a wide variety of rotation rates and profiles the gravitational wave (GW) burst signals from the core bounce are of a generic type, already known as Type I in the literature. In addition, for most models the characteristic frequency of the GW burst signal lies in a narrow range around approximately 718 Hz. In our systematic study, using both GR and Newtonian gravity, we identify, individually quantify, and discuss in detail the micro- and macrophysical mechanisms leading to this result, i.e. the effects of rotation, the EoS, and deleptonization. We also discuss the detectability prospects of such GW burst signals by GW detectors, and infer that such a generic type of signal templates will likely facilitate a more efficient search in current and future detectors of both interferometric and resonant type.

We present results from the first 2D and 3D simulations of the collapse of rotating stellar iron cores in general relativity employing a finite-temperature equation of state and an approximate treatment of deleptonization during collapse. We compare fully nonlinear and conformally flat spacetime evolution methods and find that the conformally flat treatment is sufficiently accurate for the core-collapse supernova problem. We focus on the gravitational wave (GW) emission from rotating collapse, core bounce, and early postbounce phases. Our results indicate that the GW signature of these phases is much more generic than previously estimated. In addition, we track the growth of a nonaxisymmetric instability of dominant m = 1 character in one of our models that leads to prolonged narrow-band GW emission at ~930 Hz over several tens of milliseconds.

We present results from the first 2+1 and 3+1 simulations of the collapse of rotating stellar iron cores in general relativity employing a finite-temperature equation of state and an approximate treatment of deleptonization during collapse. We compare full 3+1 and conformally-flat spacetime evolution methods and find that the conformally-flat treatment is sufficiently accurate for the core-collapse supernova problem. We focus on the gravitational wave (GW) emission from rotating collapse, core bounce, and early postbounce phases. Our results indicate that the GW signature of these phases is much more generic than previously estimated. In addition, we track the growth of a nonaxisymmetric instability of dominant m = 1 character in two of our models that leads to prolonged narrow-band GW emission at ~ 920-930 Hz over several tens of milliseconds.

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

An extended set of binary neutron star (NS) merger simulations is performed with an approximative conformally flat treatment of general relativity to systematically investigate the influence of the nuclear equation of state (EoS), the neutron star masses, and the NS spin states prior to merging. We employ the two non-zero temperature EoSs of Shen et al. (1998a,b) and Lattimer & Swesty (1991). In addition, we use the cold EoS of Akmal et al. (1998) with a simple ideal-gas-like extension according to Shibata & Taniguchi (2006), and an ideal-gas EoS with parameters fitted to the supernuclear part of the Shen-EoS. We estimate the mass sitting in a dilute high-angular momentum ``torus'' around the future black hole (BH). The dynamics and outcome of the models is found to depend strongly on the EoS and on the binary parameters. Larger torus masses are found for asymmetric systems (up to ~0.3 M_sun for a mass ratio of 0.55), for large initial NSs, and for a NS spin state which corresponds to a larger total angular momentum. We find that the postmerger remnant collapses either immediately or after a short time when employing the soft EoS of Lattimer& Swesty, whereas no sign of post-merging collapse is found within tens of dynamical timescales for all other EoSs used. The typical temperatures in the torus are found to be about 3-10 MeV depending on the strength of the shear motion at the collision interface between the NSs and thus depending on the initial NS spins. About 10^{-3}-10^{-2} M_sun of NS matter become gravitationally unbound during or right after the merging process. This matter consists of a hot/high-entropy component from the collision interface and (only in case of asymmetric systems) of a cool/low-entropy component from the spiral arm tips. (abridged)

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

We investigate the energy-momentum deposition due to neutrino-antineutrino annihilation in the vicinity of axisymmetric, accreting black holes (BHs) by numerically ray-tracing neutrino trajectories in a Kerr space-time. Hyperaccreting stellar-mass BHs are widely considered as energy sources that can drive ultrarelativistic outflows with the potential to produce gamma-ray bursts. In contrast to earlier works, we provide an extensive and detailed parameter study of the influence of general relativistic (GR) effects and of different neutrinosphere geometries. These include idealized thin disks, tori, and spheres, or are constructed as non-selfgravitating equilibrium matter distributions for varied BH rotation. Considering isothermal neutrinospheres with the same temperature and surface area, we confirm previous results that compared to Newtonian calculations, GR effects increase the annihilation rate measured by an observer at infinity by a factor of 2 when the neutrinosphere is a disk. However, in case of a torus and a sphere the influence of GR effects is globally only ~25%, although locally it can be significantly larger. Independent of whether GR effects are included, disks yield the highest energy deposition rates, followed by tori and, with the lowest rates, spheres. For disks and tori, increasing the angular momentum of the BH from 0 to 1 enhances the energy deposition rate measured by an observer at infinity by roughly a factor of 2 due to the shrinking inner radius of the neutrinosphere. (abridged)

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.

We analyze the linear stability of a stalled accretion shock in a perfect gas with a parametrized cooling function L ~ rho^{beta-alpha} P^alpha. The instability is dominated by the l=1 mode if the shock radius exceeds 2-3 times the accretor radius, depending on the parameters of the cooling function. The growth rate and oscillation period are comparable to those observed in the numerical simulations of Blondin & Mezzacappa (2006). The instability mechanism is analyzed by separately measuring the efficiencies of the purely acoustic cycle and the advective-acoustic cycle. These efficiencies are estimated directly from the eigenspectrum, and also through a WKB analysis in the high frequency limit. Both methods prove that the advective-acoustic cycle is unstable, and that the purely acoustic cycle is stable. Extrapolating these results to low frequency leads us to interpret the dominant mode as an advective-acoustic instability, different from the purely acoustic interpretation of Blondin & Mezzacappa (2006). A simplified characterization of the instability is proposed, based on an advective-acoustic cycle between the shock and the radius r_nabla where the velocity gradients of the stationary flow are strongest. The importance of the coupling region in this mechanism calls for a better understanding of the conditions for an efficient advective-acoustic coupling in a decelerated, nonadiabatic flow, in order to extend these results to core-collapse supernovae.

A toy model is analyzed in order to evaluate the linear stability of the gain region immediately behind a stalled accretion shock, after core bounce. This model demonstrates that a negative entropy gradient is not sufficient to warrant linear instability. The stability criterion is governed by the ratio \chi of the advection time through the gain region divided by the local timescale of buoyancy. The gain region is linearly stable if \chi< 3. The classical convective instability is recovered in the limit \chi\gg 3. For \chi>3, perturbations are unstable in a limited range of horizontal wavelengths centered around twice the vertical size H of the gain region. The threshold horizontal wavenumbers k_{min} and k_{max} follow simple scaling laws such that Hk_{min}\propto 1/{\chi} and Hk_{max}\propto \chi. The convective stability of the l=1 mode in spherical accretion is discussed, in relation with the asymmetric explosion of core collapse supernovae. The advective stabilization of long wavelength perturbations weakens the possible influence of convection alone on a global l=1 mode.

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.

(abridged) Three-dimensional hydrodynamic simulations with an Eulerian PPM code are presented for the time-dependent evolution of accretion tori around nonrotating and rotating stellar-mass black holes (BHs), using a pseudo-Newtonian gravitational potential to approximate the effects of general relativity. The initial configurations are assumed to be remnants of binary neutron star (NS) or NS+BH mergers and consist of a 4 solar mass BH with varied spin, girded by a torus with a mass between 0.01 and 0.2 solar masses. The evolution of tori without and with physical shear viscosity is simulated, using a realistic equation of state and following the energy loss and lepton number change due to neutrino emission by a neutrino-trapping scheme. The time-dependent efficiency of converting rest-mass energy to neutrinos is found to reach 10 percent, the efficiency of converting neutrino energy into a pair-photon fireball by neutrino annihilation can reach several percent. The rate of the latter process declines with time much less steeply than the total neutrino luminosity, because the ongoing protonization of the torus ensures a rather stable product of neutrino and antineutrino luminosities. The neutrino emission increases steeply with higher viscosity, larger torus mass, and larger BH spin in corotation with the torus. For rotation rates as expected for post-merger BHs (a > 0.5) and reasonable values of the alpha viscosity (alpha ~ 0.1), the considered tori release sufficient energy in neutrinos to account for the energetics of the well-localized short gamma-ray bursts recently detected by Hete and Swift, if collimation is invoked as predicted by hydrodynamic jet simulations.

1D and 2D supernova simulations for stars between 11 and 25 solar masses are presented, making use of the Prometheus/Vertex neutrino-hydrodynamics code, which employs a full spectral treatment of the neutrino transport. Multi-dimensional transport aspects are treated by the ``ray-by-ray plus'' approximation described in Paper I. Our set of models includes a 2D calculation for a 15 solar mass star whose iron core is assumed to rotate rigidly with an angular frequency of 0.5 rad/s before collapse. No important differences were found depending on whether random seed perturbations for triggering convection are included already during core collapse, or whether they are imposed on a 1D collapse model shortly after bounce. Convection below the neutrinosphere sets in about 40 ms p.b. at a density above 10**12 g/cm^3 in all 2D models, and encompasses a layer of growing mass as time goes on. It leads to a more extended proto-neutron star structure with accelerated lepton number and energy loss and significantly higher muon and tau neutrino luminosities, but reduced mean energies of the radiated neutrinos, at times later than ~100 ms p.b. In case of an 11.2 solar mass star we find that low (l = 1,2) convective modes cause a probably rather weak explosion by the convectively supported neutrino-heating mechanism after ~150 ms p.b. when the 2D simulation is performed with a full 180 degree grid, whereas the same simulation with 90 degree wedge fails to explode like all other models. This sensitivity demonstrates the proximity of our 2D models to the borderline between success and failure, and stresses the need of simulations in 3D, ultimately without the axis singularity of a polar grid. (abridged)

By means of two-dimensional (2D) simulations we study hydrodynamic instabilities during the first seconds of neutrino-driven supernova explosions, using a PPM hydrodynamics code, supplemented with a gray, non-equilibrium approximation of radial neutrino transport. We consider three 15 solar mass progenitors with different structures and one rotating model, in which we replace the dense core of the newly formed neutron star (NS) by a contracting inner grid boundary, and trigger neutrino-driven explosions by systematically varying the neutrino fluxes emitted at this boundary. Confirming more idealized studies as well as supernova simulations with spectral transport, we find that random seed perturbations can grow by hydrodynamic instabilities to a globally asymmetric mass distribution, leading to a dominance of dipole (l=1) and quadrupole (l=2) modes in the explosion ejecta. Anisotropic gravitational and hydrodynamic forces are found to accelerate the NS on a timescale of 2-3 seconds. Since the explosion anisotropies develop chaotically, the magnitude of the corresponding kick varies stochastically in response to small differences in the fluid flow. Our more than 70 models separate into two groups, one with high and the other with low NS velocities and accelerations after 1s of post-bounce evolution, depending on whether the l=1 mode is dominant in the ejecta or not. This leads to a bimodality of the distribution when the NS velocities are extrapolated to their terminal values. The fast group has an average velocity of about 500 km/s and peak values in excess of 1000 km/s. Establishing a link to the measured distribution of pulsar velocities, however, requires a much larger set of calculations and ultimately 3D modeling. (abridged)

Based on recent models of relativistic jet formation by thermal energy deposition around black hole-torus systems, the relation between the on- and off-axis appearance of short, hard gamma-ray bursts (GRBs) is discussed in terms of energetics, duration, average Lorentz factor, and probability of observation, assuming that the central engines are remnants of binary neutron star or neutron star-black hole mergers. As a consequence of the interaction with the torus matter at the jet basis and the subsequent expansion of the jets into an extremely low-density environment, the collimated ultrarelativistic outflows possess flat core profiles with only little variation of radially-averaged properties, and are bounded by very steep lateral edges. Owing to the rapid decrease of the isotropic-equivalent energy near the jet edges, the probability of observing the lateral, lower Lorentz factor wings is significantly reduced and most short GRBs should be seen with on-axis-like properties. Taking into account cosmological and viewing angle effects, theoretical predictions are made for the short-GRB distributions with redshift z, fluence, and isotropic-equivalent energy. The observational data for short bursts with determined redshifts are found to be compatible with the predictions only if either the intrinsic GRB rate density drops rapidly at z>1, or a large number of events at z>1 are missed, implying that the subenergetic GRB 050509b was an extremely rare low-fluence event with detectable photon flux only because of its proximity and shortness. It appears unlikely that GRB 050509b can be explained as an off-axis event. The detection of short GRBs with small Lorentz factors is statistically disfavored, suggesting a possible reason for the absence of soft short bursts in the duration-hardness diagram.

Two-dimensional simulations of strongly anisotropic supernova explosions of a nonrotating 15 solar mass blue supergiant progenitor are presented, which follow the hydrodynamic evolution from times shortly after shock formation until hours later. It is shown that explosions which around the time of shock revival are dominated by low-order unstable modes (i.e. by a superposition of the l=2 and l=1 modes, in which the former is strongest), are consistent with all major observational features of SN 1987A, in contrast to models which show high-order mode perturbations only and were published in earlier work. Among other items, the low-mode models exhibit final iron-group velocities of up to 3300 km/s, strong mixing at the He/H composition interface, with hydrogen being mixed downward in velocity space to only 500 km/s, and a final prolate anisotropy of the ejecta with a major to minor axis ratio of about 1.6. The success of low-mode explosions with an energy of about 2x10**51 erg to reproduce these observed features is based on two effects: the (by 40%) larger initial maximum velocities of metal-rich clumps compared to our high-mode models, and the initial global deformation of the shock. The latter triggers the growth of a strong Richtmyer-Meshkov instability at the He/H interface that results in a global anisotropy of the inner ejecta at late times (i.e. t > 10000 s), although the shock itself has long become spherical by then. The simulations suggest a coherent picture, which explains the observational data of SN 1987A within the framework of the neutrino-driven explosion mechanism using a minimal set of assumptions. It is therefore argued that other paradigms, which are based on (more) controversial physics, may not be required to explain this event. (abbreviated)

Merging neutron stars (NSs) are hot candidates for the still enigmatic sources of short gamma-ray bursts (GRBs). If the central engines of the huge energy release are accreting relic black holes (BHs) of such mergers, it is important to understand how the properties of the BH-torus systems, in particular disc masses and mass and rotation rate of the compact remnant, are linked to the characterizing parameters of the NS binaries. For this purpose we present relativistic smoothed particle hydrodynamics simulations with conformally flat approximation of the Einstein field equations and a physical, non-zero temperature equation of state. Thick disc formation is highlighted as a dynamical process caused by angular momentum transfer through tidal torques during the merging process of asymmetric systems or in the rapidly spinning triaxial post-merger object. Our simulations support the possibility that the first well-localized short and hard GRBs 050509b, 050709, 050724, 050813 have originated from NS merger events and are powered by neutrino-antineutrino annihilation around a relic BH-torus system. Using model parameters based on this assumption, we show that the measured GRB energies and durations lead to estimates for the accreted masses and BH mass accretion rates which are compatible with theoretical expectations. In particular, the low energy output and short duration of GRB 050509b set a very strict upper limit of less than 100 ms for the time interval after the merging until the merger remnant has collapsed to a BH, leaving an accretion torus with a small mass of only about 0.01 solar masses. This favors a (nearly) symmetric NS+NS binary with a typical mass as progenitor system.

Supernovae are nature's grandest explosions and an astrophysical laboratory in which unique conditions exist that are not achievable on Earth. They are also the furnaces in which most of the elements heavier than carbon have been forged. Scientists have argued for decades about the physical mechanism responsible for these explosions. It is clear that the ultimate energy source is gravity, but the relative roles of neutrinos, fluid instabilities, rotation and magnetic fields continue to be debated.

We present results of simulations of stellar collapse and explosions in spherical symmetry for progenitor stars in the 8-10 solar mass range with an O-Ne-Mg core. The simulations were continued until nearly one second after core bounce and were performed with the Prometheus/Vertex code with a variable Eddington factor solver for the neutrino transport, including a state-of-the-art treatment of neutrino-matter interactions. Particular effort was made to implement nuclear burning and electron capture rates with sufficient accuracy to ensure a smooth continuation, without transients, from the progenitor evolution to core collapse. Using two different nuclear equations of state (EoSs), a soft version of the Lattimer & Swesty EoS and the significantly stiffer Wolff & Hillebrandt EoS, we found no prompt explosions, but instead delayed explosions, powered by neutrino heating and the neutrino-driven baryonic wind which sets in about 200 ms after bounce. The models eject little nickel (< 0.015 solar masses), explode with an energy of about or slightly more than 10**50 erg, and leave behind neutron stars (NSs) with a baryonic mass near 1.36 solar masses. Different from previous models of such explosions, the ejecta during the first second have a proton-to-baryon ratio of Ye > 0.46, which suggests a chemical composition that is not in conflict with galactic abundances. No low-entropy matter with Ye << 0.5 is ejected. This excludes such explosions as sites of a low-entropy r-process. The low explosion energy and nucleosynthetic implications are compatible with the observed properties of the Crab supernova, and the small nickel mass supports the possibility that our models explain some subluminous Type II-P supernovae.

One of the outstanding unsolved riddles of nuclear astrophysics is the origin of the so called ``p-process'' nuclei from A = 92 to 126. Both the lighter and heavier p-process nuclei are adequately produced in the neon and oxygen shells of ordinary Type II supernovae, but the origin of these intermediate isotopes, especially 92,94Mo and 96,98Ru, has long been mysterious. Here we explore the production of these nuclei in the neutrino-driven wind from a young neutron star. We consider such early times that the wind still contains a proton excess because the rates for electron neutrino and positron captures on neutrons are faster than those for the inverse captures on protons. Following a suggestion by Frohlich et al. 2005, we also include the possibility that, in addition to the protons, alpha-particles, and heavy seed, a small flux of neutrons is maintained by the reaction p(bar(nu_e),e+)n. This flux of neutrons is critical in bridging the long waiting points along the path of the rp-process by (n,p) and (n,gamma) reactions. Using the unmodified ejecta histories from a recent two-dimensional supernova model by Janka et al. 2003, we find synthesis of p-rich nuclei up to 102Pd. However, if the entropy of these ejecta is increased by a factor of two, the synthesis extends to 120Te. Still larger increases in entropy, that might reflect the role of magnetic fields or vibrational energy input neglected in the hydrodynamical model, result in the production of numerous r-, s-, and p-process nuclei up to A approximately 170, even in winds that are proton-rich.

Supernova models with a full spectral treatment of the neutrino transport are presented, employing the Prometheus/Vertex neutrino-hydrodynamics code with a ``ray-by-ray plus'' approximation for treating two- (or three-) dimensional problems. The method is described in detail and critically assessed with respect to its capabilities, limitations, and inaccuracies in the context of supernova simulations. In this first paper of a series, 1D and 2D core-collapse calculations for a (nonrotating) 15 M_sun star are discussed, uncertainties in the treatment of the equation of state -- numerical and physical -- are tested, Newtonian results are compared with simulations using a general relativistic potential, bremsstrahlung and interactions of neutrinos of different flavors are investigated, and the standard approximation in neutrino-nucleon interactions with zero energy transfer is replaced by rates that include corrections due to nucleon recoil, thermal motions, weak magnetism, and nucleon correlations. Models with the full implementation of the ``ray-by-ray plus'' spectral transport were found not to explode, neither in spherical symmetry nor in 2D with a 90 degree lateral wedge. The success of previous 2D simulations with grey, flux-limited neutrino diffusion can therefore not be confirmed. Omitting the radial velocity terms in the neutrino momentum equation leads to ``artificial'' explosions by increasing the neutrino energy density in the convective gain layer by about 20--30% and thus the integral neutrino energy deposition in this region by about a factor of two. (abbreviated)

Based on new developments in the understanding of supernovae (SNe) as gravitational-wave (GW) sources we estimate the GW background from all cosmic SNe. For a broad range of frequencies around 1 Hz, this background is crudely comparable to the GW background expected from standard inflationary models. While our estimate remains uncertain within several orders of magnitude, the SN GW background may become detectable by second-generation space-based interferometers such as the proposed Big-Bang Observatory (BBO). By the same token, the SN GWs may become a foreground for searches of the inflationary GWs, in particular for sub-Hz frequencies where the SN background is Gaussian and where the BBO will be most sensitive. SN simulations lasting far beyond the usual cutoff of about 1 second are needed for more robust predictions in the sub-Hz frequency band. An even larger GW background can arise from a hypothetical early population of massive stars, although their GW source strength as well as their abundance are currently poorly understood.

We investigate the possibility to approximate relativistic effects in hydrodynamical simulations of stellar core collapse and post-bounce evolution by using a modified gravitational potential in an otherwise standard Newtonian hydrodynamic code. Different modifications of the effective relativistic potential introduced by Rampp & Janka (2002) are discussed. Corresponding hydrostatic solutions are compared with solutions of the TOV equations, and hydrodynamic simulations with two different codes are compared with fully relativistic results. One code is applied for one- and two-dimensional calculations with a simple equation of state, and employs either the modified effective relativistic potential in a Newtonian framework or solves the general relativistic field equations under the assumption of the conformal flatness condition (CFC) for the three-metric. The second code allows for full-scale supernova runs including a microphysical equation of state and neutrino transport based on the solution of the Boltzmann equation and its moments equations. We present prescriptions for the effective relativistic potential for self-gravitating fluids to be used in Newtonian codes, which produce excellent agreement with fully relativistic solutions in spherical symmetry, leading to significant improvements compared to previously published approximations. Moreover, they also approximate qualitatively well relativistic solutions for models with rotation.

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.

One of the robust features found in simulations of core-collapse supernovae (SNe) is the prompt neutronization burst, i.e. the first $\sim 25$ milliseconds after bounce when the SN emits with very high luminosity mainly $\nu_e$ neutrinos. We examine the dependence of this burst on variations in the input of current SN models and find that recent improvements of the electron capture rates as well as uncertainties in the nuclear equation of state or a variation of the progenitor mass have only little effect on the signature of the neutronization peak in a megaton water Cherenkov detector for different neutrino mixing schemes. We show that exploiting the time-structure of the neutronization peak allows one to identify the case of a normal mass hierarchy and large 13-mixing angle $\theta_{13}$, where the peak is absent. The robustness of the predicted total event number in the neutronization burst makes a measurement of the distance to the SN feasible with a precision of about 5%, even in the likely case that the SN is optically obscured.

New simulations demonstrate that low-mode, nonradial hydrodynamic instabilities of the accretion shock help starting hot-bubble convection in supernovae and thus support explosions by the neutrino-heating mechanism. The prevailing conditions depend on the high-density equation of state which governs stellar core collapse, core bounce, and neutron star formation. Tests of this sensitivity to nuclear physics variations are shown for spherically symmetric models. Implications of current explosion models for r-process nucleosynthesis are addressed.

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.

The identification of the astrophysical site and the specific conditions in which r-process nucleosynthesis takes place remain unsolved mysteries of astrophysics. The present paper emphasizes some important future challenges faced by nuclear physics in this problem, particularly in the determination of the radiative neutron capture rates by exotic nuclei close to the neutron drip line and the fission probabilities of heavy neutron-rich nuclei. These quantities are particularly relevant to determine the composition of the matter resulting from the decompression of initially cold neutron star matter. New detailed r-process calculations are performed and the final composition of ejected inner and outer neutron star crust material is estimated. We discuss the impact of the many uncertainties in the astrophysics and nuclear physics on the final composition of the ejected matter. The similarity between the predicted and the solar abundance pattern for A > 140 nuclei as well as the robustness of the prediction with varied input parameters makes this scenario one of the most promising that deserves further exploration.

A few seconds after bounce in a core-collapse supernova, the shock wave passes the density region corresponding to resonant neutrino oscillations with the ``atmospheric'' neutrino mass difference. The transient violation of the adiabaticity condition manifests itself in an observable modulation of the neutrino signal from a future galactic supernova. In addition to the shock wave propagation effects that were previously studied, a reverse shock forms when the supersonically expanding neutrino-driven wind collides with the slower earlier supernova ejecta. This implies that for some period the neutrinos pass two subsequent density discontinuities, giving rise to a ``double dip'' feature in the average neutrino energy as a function of time. We study this effect both analytically and numerically and find that it allows one to trace the positions of the forward and reverse shocks. We show that the energy dependent neutrino conversion probabilities allow one to detect oscillations even if the energy spectra of different neutrino flavors are the same as long as the fluxes differ. These features are observable in the \bar\nu_e signal for an inverted and in the \nu_e signal for a normal neutrino mass hierarchy, provided the 13-mixing angle is ``large'' (sin^2\theta_{13}\gg 10^{-5}).

As an explosion develops in the collapsed core of a massive star, neutrino emission drives convection in a hot bubble of radiation, nucleons, and pairs just outside a proto-neutron star. Shortly thereafter, neutrinos drive a wind-like outflow from the neutron star. In both the convective bubble and the early wind, weak interactions temporarily cause a proton excess (Ye > 0.5) to develop in the ejected matter. This situation lasts for at least the first second, and the approximately 0.05 - 0.1 solar masses that is ejected has an unusual composition that may be important for nucleosynthesis. Using tracer particles to follow the conditions in a two-dimensional model of a successful supernova explosion calculated by Janka, Buras & Rampp (2003), we determine the composition of this material. Most of it is helium and 56Ni. The rest is relatively rare species produced by the decay of proton-rich isotopes unstable to positron emission. In the absence of pronounced charged-current neutrino capture, nuclear flow will be held up by long-lived waiting point nuclei in the vicinity of 64Ge. The resulting abundance pattern can be modestly rich in a few interesting rare isotopes like 45Sc, 49Ti, and 64Zn. The present calculations imply yields that, when compared with the production of major species in the rest of the supernova, are about those needed to account for the solar abundance of 45Sc and 49Ti. Since the synthesis will be nearly the same in stars of high and low metallicity, the primary production of these species may have discernible signatures in the abundances of low metallicity stars. We also discuss uncertainties in the nuclear physics and early supernova evolution to which abundances of interesting nuclei are sensitive.

Two- and three-dimensional simulations demonstrate that hydrodynamic instabilities can lead to low-mode (l=1,2) asymmetries of the fluid flow in the neutrino-heated layer behind the supernova shock. This provides a natural explanation for aspherical mass ejection and for pulsar recoil velocities even in excess of 1000 km/s. We propose that the bimodality of the pulsar velocity distribution might be a consequence of a dominant l=1 mode in case of the fast component, while higher-mode anisotropy characterizes the postshock flow and SN ejecta during the birth of the slow neutron stars. We argue that the observed large asymmetries of supernovae and the measured high velocities of young pulsars therefore do not imply rapid rotation of the iron core of the progenitor star, nor do they require strong magnetic fields to play a crucial role in the explosion. Anisotropic neutrino emission from accretion contributes to the neutron star acceleration on a minor level, and pulsar kicks do not make a good case for non-standard neutrino physics in the nascent neutron star.

We present the first general relativistic hydrodynamic models of the launch and evolution of relativistic jets and winds, driven by thermal energy deposition, possibly due to neutrino-antineutrino annihilation, in the close vicinity of black hole-accretion torus systems. The latter are considered to be the remnants of compact object mergers. Our two-dimensional simulations establish the link between such mergers and future observations of short gamma-ray bursts (GRBs) by the SWIFT satellite. They show that ultrarelativistic outflow with maximum terminal Lorentz factors (Gamma) around 1000 develops for polar energy deposition rates above some 1e48 erg/s per steradian, provided the merger environment has a sufficiently low baryon density. Due to the collimation by the dense accretion torus the typical semi-opening angles of the Gamma > 100 cone are 5-10 degrees, corresponding to about 0.4-1.5% of the hemisphere and apparent isotropized energies (kinetic plus internal) up to ~1e51 erg. 10-30% of the deposited energy are transferred to the outflow with Gamma > 100. Our models confirm the viability of post-merger BH-torus systems as engines of short, hard GRBs and can explain the durations of all observed short GRBs, because different propagation velocities of the front and rear ends lead to a radial stretching of the ultrarelativistic fireball before transparency is reached. The ultrarelativistic flow reveals a highly non-uniform structure with Lorentz factor variations up to factors of a few, caused by the action of Kelvin-Helmholtz instabilities that originate at the fireball-torus interface (abbreviated).

We show by two-dimensional and first three-dimensional simulations of neutrino-driven supernova explosions that low (l=1,2) modes can dominate the flow pattern in the convective postshock region on timescales of hundreds of milliseconds after core bounce. This can lead to large global anisotropy of the supernova explosion and pulsar kicks in excess of 500 km/s.

We briefly summarize recent efforts in Garching for modeling stellar core collapse and post-bounce evolution in one and two dimensions. The transport of neutrinos of all flavors is treated by iteratively solving the coupled system of frequency-dependent moment equations together with a model Boltzmann equation which provides the closure. A variety of progenitor stars, different nuclear equations of state, stellar rotation, and global asymmetries due to large-mode hydrodynamic instabilities have been investigated to ascertain the road to finally successful, convectively supported neutrino-driven explosions.

We present three-dimensional hydrodynamic simulations of the evolution of selfgravitating, thick accretion discs around hyperaccreting stellar-mass black holes. The black hole-torus systems are considered to be remnants of compact object mergers, in which case the disc is not fed by an external mass reservoir and the accretion is non-stationary. Our models take into account viscous dissipation, described by an alpha-law, a detailed equation of state for the disc gas, and an approximate treatment of general relativistic effects on the disc structure by using a pseudo-Newtonian potential for the black hole including its possible rotation and spin-up during accretion. Magnetic fields are ignored. The neutrino emission of the hot disc is treated by a neutrino-trapping scheme, and the neutrino-antineutrino annihilation near the disc is evaluated in a post-processing step. Our simulations show that the neutrino emission and energy deposition by neutrino-antineutrino annihilation increase sensitively with the disc mass, with the black hole spin in case of a disc in corotation, and in particular with the alpha-viscosity. We find that for sufficiently large alpha-viscosity neutrino-antineutrino annihilation can be a viable energy source for gamma-ray bursts.

Our current knowledge of neutron star formation, progenitors, and natal masses, spins, magnetic fields, and space velocities is briefly reviewed from a theorist's perspective. More observational information is badly needed to constrain theoretical possibilities.

Recent progress in modeling core-collapse supernovae is summarized and set in perspective. Two-dimensional simulations with state-of-the-art treatment of neutrino transport still fail to produce powerful explosions, but evidence is presented that they are very close to success.

Assuming that the neutrino luminosity from the neutron star core is sufficiently high to drive supernova explosions by the neutrino-heating mechanism, we show that low-mode (l = 1, 2) convection can develop from random seed perturbations behind the shock. A slow onset of the explosion is crucial, requiring the core luminosity to vary slowly with time, in contrast to the burst-like exponential decay assumed in previous work. Gravitational and hydrodynamic forces by the globally asymmetric supernova ejecta were found to accelerate the remnant neutron star on a timescale of more than a second to velocities above 500 km/s, in agreement with observed pulsar proper motions.

We have computed the gravitational wave signal from supernova core collapse using the presently most realistic input physics available. We start from state-of-the-art progenitor models of rotating and non-rotating massive stars, and simulate the dynamics of their core collapse by integrating the equations of axisymmetric hydrodynamics together with the Boltzmann equation for the neutrino transport including an elaborate description of neutrino interactions, and a realistic equation of state. We compute the quadrupole wave amplitudes, the Fourier wave spectra, the amount of energy radiated in form of gravitational waves, and the S/N ratios for the LIGO and the tuned Advanced LIGO interferometers resulting both from non-radial mass motion and anisotropic neutrino emission. The simulations demonstrate that the dominant contribution to the gravitational wave signal is produced by neutrino-driven convection behind the supernova shock. For stellar cores rotating at the extreme of current stellar evolution predictions, the core-bounce signal is detectable with advanced LIGO up to a distance of 5kpc, whereas the signal from post-shock convection is observable up to a distance of about 100kpc. If the core is non-rotating its gravitational wave emission can be measured up to a distance of 15kpc, while the signal from the Ledoux convection in the deleptonizing, nascent neutron star can be detected up to a distance of 10kpc. Both kinds of signals are generically produced by convection in any core collapse supernova.

Two-dimensional simulations of a Type II and a Type Ib-like supernova explosion are presented that encompass shock revival by neutrino heating, neutrino-driven convection, explosive nucleosynthesis, the growth of Rayleigh-Taylor instabilities, and the propagation of newly formed metal clumps through the exploding star. In both cases we find very high Ni56 velocities of 17000 km/s shortly after shock-revival, and a complete fragmentation of the progenitor's metal core within the first few minutes after core bounce, due to the growth of Rayleigh-Taylor instabilities at the Si/O and (C+O)/He composition interfaces. This leads to the formation of high-velocity, metal-rich clumps which eventually decouple from the flow and move ballistically through the ejecta. Maximum final metal velocities of 3500-5500 km/s and 1200 km/s are obtained for the Type Ib model and the Type II model, respectively. The low maximum metal velocities in the Type II model, which are significantly smaller than those observed in SN 1987A, are due to the massive hydrogen envelope of the progenitor. The envelope forces the supernova shock to decelerate strongly, leaving behind a reverse shock below the He/H interface, which interacts with the clumps and slows them down significantly. This reverse shock is absent in the Type Ib-like model. The latter is in fairly good agreement with observations of Type Ib supernovae. In addition, in this case the pattern of clump formation in the ejecta is correlated with the convective pattern prevailing during the shock-revival phase. This might be used to deduce observational constraints for the dynamics during this early phase of the evolution, and the role of neutrino heating in initiating the explosion.

Transporting nu_mu and nu_tau in a supernova (SN) core involves several processes that have been neglected in traditional simulations. Based on a Monte Carlo study we find that the flavor-dependent spectral differences are much smaller than is often stated in the literature. A full-scale SN simulation using a Boltzmann solver and including all relevant neutrino reactions confirms these results. The flavor-dependent flux differences are largest during the initial accretion phase.

Two-dimensional hydrodynamic simulations of stellar core-collapse with and without rotation are presented which for the first time were performed by solving the Boltzmann equation for the neutrino transport including a state-of-the-art description of neutrino interactions. Although convection develops below the neutrinosphere and in the neutrino-heated region behind the supernova shock, the models do not explode. This suggests missing physics, possibly with respect to the nuclear equation of state and weak interactions in the subnuclear regime. However, it might also indicate a fundamental problem of the neutrino-driven explosion mechanism.

The neutrino flux and spectra formation in a supernova core is studied by using a Monte Carlo code. The dominant opacity contribution for nu_mu and nu_tau is elastic scattering on nucleons. In addition we switch on or off a variety of processes which allow for the exchange of energy or the creation and destruction of neutrino pairs, notably nucleon bremsstrahlung, the e^+ e^- pair annihilation process and nu_e-bar nu_e -> nu_{mu,tau} nu_{mu,tau}-bar, recoil and weak magnetism in elastic nucleon scattering, elastic scattering on electrons and positrons and elastic scattering on electron neutrinos and anti-neutrinos. The least important processes are neutrino-neutrino scattering and e^+ e^- annihilation. The formation of the spectra and fluxes of nu_mu is dominated by the nucleonic processes, i.e. bremsstrahlung and elastic scattering with recoil, but also nu_e nu_e-bar annihilation and nu_mu e^\pm scattering contribute significantly. When all processes are included, the spectral shape of the emitted neutrino flux is always ``pinched,'' i.e. the width of the spectrum is smaller than that of a thermal spectrum with the same average energy. In all of our cases we find that the average nu_mu-bar energy exceeds the average nu_e-bar energy by only a small amount, 10% being a typical number. Weak magnetism effects cause the opacity of nu_mu to differ slightly from that of nu_mu-bar, translating into differences of the luminosities and average energies of a few percent. Depending on the density, temperature, and composition profile, the flavor-dependent luminosities L_{nu_e}$, L_{nu_e-bar}, and L_{nu_mu} can mutually differ from each other by up to a factor of two in either direction.

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.

One of the central problems in supernova theory is the question how massive stars explode. Understanding the physical processes that drive the explosion is crucial for linking the stellar progenitors to the final remnants and for predicting observable properties like explosion energies, neutron star and black hole masses, nucleosynthetic yields, explosion anisotropies, and pulsar kicks. Here we review different suggestions for the explosion mechanism and discuss the constraints that can or cannot be deduced from observations. The most refined and detailed hydrodynamical models are available for neutrino-driven supernova explosions, although the viability of this mechanism has yet to be demonstrated convincingly. Since spherical models do not explode, the hope rests on the helpful effects of convection inside the nascent neutron star and in the neutrino-heated region behind the stalled shock. We present the first two-dimensional simulations of these processes which have been performed with a Boltzmann solver for the neutrino transport and a state-of-the-art description of neutrino-matter interactions. Our most complete models fail, but convection brings them encouragingly close to a success. An explosion could be obtained by just a minor modification of the neutrino transport, in which case the exploding model fulfills important requirements from observations (abridged).

The rapidly growing base of observational data for supernova explosions of massive stars demands theoretical explanations. Central of these is a self-consistent model for the physical mechanism that provides the energy to start and drive the disruption of the star. We give arguments why the delayed neutrino-heating mechanism should still be regarded as the standard paradigm to explain most explosions of massive stars and show how large-scale and even global asymmetries can result as a natural consequence of convective overturn in the neutrino-heating region behind the supernova shock. Since the explosion is a threshold phenomenon and depends sensitively on the efficiency of the energy transfer by neutrinos, even relatively minor differences in numerical simulations can matter on the secular timescale of the delayed mechanism. To enhance this point, we present some results of recent one- and two-dimensional computations, which we have performed with a Boltzmann solver for the neutrino transport and a state-of-the-art description of neutrino-matter interactions. Although our most complete models fail to explode, the simulations demonstrate that one is encouragingly close to the critical threshold because a modest variation of the neutrino transport in combination with postshock convection leads to a weak neutrino-driven explosion with properties that fulfill important requirements from observations.

We present here the first results of two-dimensional hydrodynamical simulations of the neutrino-heating phase in the collapsed core of a 15 solar mass star, where the neutrino transport is treated with a variable Eddington factor method for solving the Boltzmann transport equation, and the neutrino interactions include nucleon-nucleon bremsstrahlung, nucleon recoils and correlations, and weak-magnetism effects as well as direct interactions between neutrinos of different flavors. With the given input physics (neutrino reactions and nuclear equation of state), our best simulations do not develop strong convection in the neutrino-heating layer behind the shock and do not yield explosions. With about 30% higher neutrino-energy deposition behind the shock, however, an explosion occurs on a timescale of 150 ms after core bounce. It leaves behind a neutron star with an initial baryonic mass of 1.4 solar masses and ejects the N = 50 isotopes of Sr, Y and Zr in amounts consistent with Galactic abundances.

Neutrino transport and neutrino interactions in dense matter play a crucial role in stellar core collapse, supernova explosions and neutron star formation. Here we present a detailed description of a new numerical code for treating the time and energy dependent neutrino transport in hydrodynamical simulations of such events. The code is based on a variable Eddington factor method to deal with the integro-differential character of the Boltzmann equation. The moments of the neutrino distribution function and the energy and lepton number exchange with the stellar medium are determined by iteratively solving the zeroth and first order moment equations in combination with a model Boltzmann equation. The latter is discretized on a grid of tangent rays. The integration of the transport equations and the neutrino source terms is performed in a time-implicit way. In the present version of the program, the transport part is coupled to an explicit hydrodynamics code which follows the evolution of the stellar plasma by a finite-volume method with piecewise parabolic interpolation, using a Riemann solver for calculating the hydrodynamic states. The neutrino source terms are implemented in an operator-split step. Neutrino transport and hydrodynamics can be calculated with different spatial grids and different time steppings. The structure of the described code is modular and offers a high degree of flexibility for an application to relativistic and multi-dimensional problems at different levels of refinement and accuracy. We critically evaluate results for a number of test cases, including neutrino transport in rapidly moving stellar media and approximate relativistic core collapse, and suggest a path for generalizing the code to be used in multi-dimensional simulations of convection in neutron stars and supernovae.

We show that in a supernova core the annihilation process nu_e nu_e-bar -> nu_{mu,tau} nu_{mu,tau}-bar is always more important than the traditional reaction e^+ e^- -> nu_{mu,tau} nu_{mu,tau}-bar as a source for muon and tau neutrino pairs. We study the impact of the new process by means of a Monte Carlo transport code with a static stellar background model and by means of a self-consistent hydrodynamical simulation with Boltzmann neutrino transport. Nucleon bremsstrahlung NN -> NN nu_{mu,tau} nu_{mu,tau}-bar is also included as another important source term. Taking into account nu_e nu_e-bar -> nu_{mu,tau} nu_{mu,tau}-bar increases the nu_mu and nu_tau luminosities by as much as 20% while the spectra remain almost unaffected. In our hydrodynamical simulation the shock was somewhat weakened. Elastic nu_{mu,tau} nu_e and nu_{mu,tau} nu_e scattering is not negligible but less important than nu_{mu,tau} e^+ or e^- scattering. Its influence on the nu_{mu,tau} fluxes and spectra is small after all other processes have been included.

Spherically symmetric simulations of stellar core collapse and post-bounce evolution are used to test the sensitivity of the supernova dynamics to different variations of the input physics. We consider a state-of-the-art description of the neutrino-nucleon interactions, possible lepton-number changing neutrino reactions in the neutron star, and the potential impact of hydrodynamic mixing behind the supernova shock.

Supermassive black holes are now realized to exist in the centers of most galaxies. The recent discoveries of luminous quasars at redshifts higher than 6 require that these black holes were assembled already when the Universe was less than a billion years old. They might originate from the collapse of supermassive stars, a scenario which could ensure a sufficiently rapid formation. Supermassive stars are dominated by photon pressure and radiate at their Eddington limit, which drives their quasi-static evolution to a final relativistic instability. Above some critical value of the metallicity, their collapse can lead to a gigantic explosion, powered by the energy release due to hydrogen burning, but below this critical metallicity their collapse inevitably ends in the formation of a black hole, accompanied by the emission of huge amounts of energy in the form of neutrinos. Although collapsing supermassive stars are the most powerful known burst sources of neutrinos, the associated conditions do not appear favorable for producing highly relativistic outflows that can explain cosmic gamma-ray bursts.

(Abridged) In this paper we present a compilation of results from our most advanced neutron star merger simulations, including a description of the employed numerical procedures and a more complete overview over a large number of computed models. The three-dimensional hydrodynamic simulations were done with a code based on the Piecewise Parabolic Method with up to five levels of nested Cartesian grids. The simulations are basically Newtonian, but gravitational-wave emission and the corresponding back-reaction are taken into account. The use of a physical nuclear equation of state allows us to follow the thermodynamic history of the stellar medium and to compute the energy and lepton number loss due to the emission of neutrinos. The computed models differ concerning the neutron star masses and mass ratios, the neutron star spins, the numerical resolution expressed by the cell size of the finest grid and the number of grid levels, and the calculation of the temperature from the solution of the entropy equation instead of the energy equation. Our simulations show that the details of the gravitational-wave emission are still sensitive to the numerical resolution, even in our highest-quality calculations. The amount of mass which can be ejected from neutron star mergers depends strongly on the angular momentum of the system. Our results do not support the initial conditions of temperature and proton-to-nucleon ratio assumed in recent work for producing a solar r-process pattern for nuclei around and above the A approx 130 peak. The improved models confirm our previous conclusion that gamma-ray bursts are not powered by neutrino emission during the dynamical phase of the merging of two neutron stars.

We present the results of numerical simulations of the spherically symmetric gravitational collapse of supermassive stars (SMS). The collapse is studied using a general relativistic hydrodynamics code. The coupled system of Einstein and fluid equations is solved employing observer time coordinates, by foliating the spacetime by means of outgoing null hypersurfaces. The code contains an equation of state which includes effects due to radiation, electrons and baryons, and detailed microphysics to account for electron-positron pairs. In addition energy losses by thermal neutrino emission are included. We are able to follow the collapse of SMS from the onset of instability up to the point of black hole formation. Several SMS with masses in the range $5\times 10^5 M_{\odot}- 10^9 M_{\odot}$ are simulated. In all models an apparent horizon forms initially, enclosing the innermost 25% of the stellar mass. From the computed neutrino luminosities, estimates of the energy deposition by $\nu\bar{\nu}$-annihilation are obtained. Only a small fraction of this energy is deposited near the surface of the star, where, as proposed recently by Fuller & Shi (1998), it could cause the ultrarelativistic flow believed to be responsible for $\gamma$-ray bursts. Our simulations show that for collapsing SMS with masses larger than $5\times 10^5 M_{\odot}$ the energy deposition is at least two orders of magnitude too small to explain the energetics of observed long-duration bursts at cosmological redshifts. In addition, in the absence of rotational effects the energy is deposited in a region containing most of the stellar mass. Therefore relativistic ejection of matter is impossible.

The current picture of the collapse and explosion of massive stars and the formation of neutron stars is reviewed. According to the favored scenario, however by no means proven and undisputed, neutrinos deposit the energy of the explosion in the stellar medium which surrounds the nascent neutron star. Observations, in particular of Supernova~1987A, suggest that mixing processes play an important role in the expanding star, and multi-dimensional simulations show that these are linked to convective instabilities in the immediate vicinity of the neutron star. Convectively enhanced energy transport inside the neutron star can have important consequences for the neutrino emission and thus the neutrino-heating mechanism. This also holds for a suppression of the neutrino interactions at nuclear densities. Multi-dimensional hydrodynamics, general relativity, and a better understanding of the neutrino interactions in neutron star matter may be crucial to resolve the problem that state-of-the-art spherical models do not yield explosions even with a very accurate treatment of neutrino transport by solving the Boltzmann equation.

Merging neutron stars and neutron star-black holes binaries are powerful sources of gravitational waves. They have also been suggested as possible sources of cosmic gamma-ray bursts and are discussed as sites for the formation of r-process elements. Whereas the first aspect is undoubtable, the latter two are very uncertain. The current status of our knowledge and of the numerical modeling is briefly reviewed and the results of simulations are critically discussed concerning their significance and implications for potentially observable signals.

Energy deposition by neutrinos can rejuvenate the stalled bounce shock and can provide the energy for the supernova explosion of a massive star. This neutrino-heating mechanism, however, is not finally accepted or proven as the trigger of the explosion. Part of the problem is that the complexity of the hydrodynamic models often hampers a clear and simple interpretation of the results. This demands a deeper theoretical understanding of the requirements of a successful shock revival. A toy model is presented here for discussing the neutrino heating phase analytically by a time-dependent treatment, which allows one to calculate the radius and velocity of the supernova shock from global properties of the gain layer as solutions of an initial value problem. A criterion is derived for the requirements of shock revival. It confirms the existence of a minimum neutrino luminosity needed for shock expansion, but also demonstrates the importance of a sufficiently large mass infall rate to the shock. The possibility of very energetic neutrino-driven explosions seems excluded because the total specific energy transferred to nucleons is limited by about 1e52 ergs per solar mass (about 5 MeV per nucleon) and the total mass in the gain layer is typically only around 0.1 solar masses. Energy transport by convection from the region of maximum heating to radii closer behind the shock is found to support the explosion by reducing the energy loss associated with the inward advection of neutrino-heated matter through the gain radius. (abridged)

We study if the neutrino mixing parameters suggested by the atmospheric neutrino anomaly imply chemical equilibrium between mu- and tau-flavored leptons in a supernova (SN) core. The initial flavor-conversion rate would indeed be fast if the nu_mu-nu_tau-mixing angle were not suppressed by second-order refractive effects. The neutrino diffusion coefficients are different for nu_mu, anti-nu_mu, nu_tau and anti-nu_tau so that neutrino transport will create a net mu and tau lepton number density. This will typically lead to a situation where the usual first-order refractive effects dominate, further suppressing the rate of flavor conversion. Altogether, neutrino refraction has the nontrivial consequence of guaranteeing the separate conservation of e, mu, and tau lepton number in a SN core on the infall and cooling time scales, even when neutrino mixing angles are large.

New two-dimensional, high-resolution calculations of a core collapse supernova in a 15 Msol blue supergiant are presented, which cover the entire evolution from shock revival until the first few hours of the explosion. Explosive nucleosynthesis, its dependence upon convective mixing during the first second of the evolution and the growth of Rayleigh-Taylor instabilities at the composition interfaces of the progenitor star are all modeled consistently and allow for a comparison with observational data. We confirm our earlier findings, that the perturbations induced by neutrino driven convection are sufficiently strong to seed large-scale Rayleigh-Taylor mixing and to destroy the onion-shell structure of the stellar He-core. As in our earlier calculations, the strong deceleration of the nickel clumps in the layers adjacent to the He/H interface suggests that the high velocities of iron-group elements observed in SN 1987A cannot be explained on the basis of currently favored progenitor models. Possible solutions to this dilemma and the implications of the mixing for type Ib explosions are briefly discussed.

We present a spherically symmetric, Newtonian core-collapse simulation of a 15 solar mass star with a 1.28 solar mass iron core. The time-, energy-, and angle-dependent transport of electron neutrinos and antineutrinos was treated with a new code which iteratively solves the Boltzmann equation and the equations for neutrino number, energy and momentum to order O(v/c) in the velocity v of the stellar medium. The supernova shock expands to a maximum radius of 350 km instead of only about 240 km as in a comparable calculation with multi-group flux-limited diffusion (MGFLD) by Bruenn, Mezzacappa, & Dineva (1995). This may be explained by stronger neutrino heating due to the more accurate transport in our model. Nevertheless, after 180 ms of expansion the shock finally recedes to a radius around 250 km (compared to about 170 km in the MGFLD run). The effect of an accurate neutrino transport is helpful, but not large enough to cause an explosion of the considered 15 solar mass star. Therefore postshock convection and/or an enhancement of the core neutrino luminosity by convection or reduced neutrino opacities in the neutron star seem necessary for neutrino-driven explosions of such stars. We find an electron fraction Y_e > 0.5 in the neutrino-heated matter, which suggests that the overproduction problem of neutron-rich nuclei with mass numbers around A = 90 in exploding models may be absent when a Boltzmann solver is used for the electron neutrino and antineutrino transport.

High-resolution two-dimensional simulations were performed for the first five minutes of the evolution of a core collapse supernova explosion in a 15 solar mass blue supergiant progenitor. The computations start shortly after bounce and include neutrino-matter interactions by using a light-bulb approximation for the neutrinos, and a treatment of the nucleosynthesis due to explosive silicon and oxygen burning. We find that newly formed iron-group elements are distributed throughout the inner half of the helium core by Rayleigh-Taylor instabilities at the Ni+Si/O and C+O/He interfaces, seeded by convective overturn during the early stages of the explosion. Fast moving nickel mushrooms with velocities up to about 4000 km/s are observed. This offers a natural explanation for the mixing required in light curve and spectral synthesis studies of Type Ib explosions. A continuation of the calculations to later times, however, indicates that the iron velocities observed in SN 1987 A cannot be reproduced because of a strong deceleration of the clumps in the dense shell left behind by the shock at the He/H interface.

We present results of high-resolution two-dimensional simulations which follow the first five minutes of a core collapse supernova explosion in a 15 solar mass blue supergiant progenitor. The computations start shortly after core bounce and include neutrino-matter interactions by using a light-bulb approximation for the neutrinos, and a treatment of the nucleosynthesis due to explosive silicon and oxygen burning. We find that newly formed iron-group elements are distributed throughout a significant fraction of the stellar helium core by the concerted action of convective and Rayleigh-Taylor instabilities. Fast moving nickel mushrooms with velocities up to 4000 km/s are observed. A continuation of the calculations to later times, however, indicates, that the iron velocities observed in SN 1987 A cannot be reproduced due to a strong deceleration of the clumps during their interaction with the dense shell left behind by the shock at the He/H interface. Therefore, we cannot confirm the claim that convective "premixing" of the nickel in the early phases of the explosion solves the problem of the high iron velocities.

Hydrodynamic simulations of the merger of stellar mass black hole - neutron star binaries (BH/NS) are compared with mergers of binary neutron stars (NS/NS). The simulations are Newtonian, but take into account the emission and backreaction of gravitational waves. The use of a physical nuclear equation of state allows us to include the effects of neutrino emission. For low neutron star to black hole mass ratios the neutron star transfers mass to the black hole during a few cycles of orbital decay and subsequent widening before finally being disrupted, whereas for ratios near unity the neutron star is already distroyed during its first approach. A gas mass between about 0.3 and about 0.7 solar masses is left in an accretion torus around the black hole and radiates neutrinos at a luminosity of several 10^{53} erg/s during an estimated accretion time scale of about 0.1 s. The emitted neutrinos and antineutrinos annihilate into electron-positron pairs with efficiencies of 1-3% percent and rates of up to 2*10^{52} erg/s, thus depositing an energy of up to 10^{51} erg above the poles of the black hole in a region which contains less than 10^{-5} solar masses of baryonic matter. This could allow for relativistic expansion with Lorentz factors around 100 and is sufficient to explain apparent burst luminosities of up to several 10^{53} erg/s for burst durations of approximately 0.1-1 s, if the gamma emission is collimated in two moderately focussed jets in a fraction of about 1/100-1/10 of the sky.

Radioactive isotopes such as 44Ti offer probably the most direct probe into nucleosynthesis environments in the interior of exploding stars, when the associated gamma-ray activities in the explosion remnant are detected and translated back to the isotopic abundances at the time of the explosion. In this paper, we assert that the procedure may not necessarily be straightforward, at least in the case of 44Ti, an orbital-electron capture decay isotope. Using the analytic model of McKee and Truelove (1995) for young supernova remnants, and assuming the existence of overdense 56Fe-dominated clumps that contain also 44Ti, we show that a high degree of ionization may be caused by the reverse shock so that the electron-capture rate of 44Ti could be significantly reduced from its laboratory value. When applied to Cas A, this increases under certain conditions the current 44Ti-activity by a factor 1.5 - 2.5, which yields a better compatibility between the COMPTEL observation of the 1.16 MeV line activity associated with the 44Ti decay and the SN model predictions of the initial 44Ti abundance. This possibility is, however, subject to various uncertainties, and in particular to the unknown properties and radial distribution of the the clumps in the ejecta.

In a supernova core, magnetic fields cause a directional variation of the neutrino refractive index so that resonant flavor oscillations would lead to a deformation of the "neutrinosphere" for, say, tau neutrinos. The associated anisotropic neutrino emission was proposed as a possible origin of the observed pulsar proper motions. We argue that this effect was vastly overestimated because the variation of the temperature over the deformed neutrinosphere is not an adequate measure for the anisotropy of neutrino emission. The neutrino flux is generated inside the neutron star core and is transported through the atmosphere at a constant luminosity, forcing the temperature gradient in the atmosphere to adjust to the inflow of energy from below. Therefore, no emission anisotropy is caused by a deformation of the neutrinosphere to lowest order. An estimate of the higher-order corrections must take into account the modified atmospheric temperature profile in response to the deformation of the neutrinosphere and the corresponding feedback on the core. We go through this exercise in the framework of a simplified model which can be solved analytically.

Multi-dimensional hydrodynamic simulations of the post-bounce evolution of collapsed stellar iron cores have demonstrated that convective overturn between the stalled shock and the neutrinosphere can have an important effect on the neutrino-driven explosion mechanism. Whether a model yields a successful explosion or not, however, still depends on the power of neutrino energy deposition behind the stalled shock. The neutrino interaction with the stellar gas in the ``hot bubble'' also determines the duration of the shock stagnation phase, the explosion energy, and the composition of the neutrino-heated supernova ejecta. More accurate models require a more precise calculation of the neutrino luminosities and spectra and of the angular distributions of the neutrinos in the heating region. Therefore it is necessary to improve the numerical treatment of the neutrino transport, to develop a better understanding of the neutrino opacities of the dense nuclear medium, and to take into account convective processes {\it inside} the newly formed neutron star.

Hydrodynamic simulations were performed of the dynamical phase of the merging of binary neutron stars (NS-NS) and of neutron star black hole binaries (NS-BH), using the physical nuclear equation of state of Lattimer and Swesty (1991) and taking into account the emission of gravitational waves and neutrinos.

By means of three-dimensional hydrodynamic simulations with a Eulerian PPM code we investigate the formation and the properties of the accretion torus around the stellar mass black hole which originates from the merging of two neutron stars. The simulations are performed with four nested cartesian grids which allow for both a good resolution near the central black hole and a large computational volume. They include the use of a physical equation of state as well as the neutrino emission from the hot matter of the torus. The gravity of the black hole is described with a Newtonian and alternatively with a Paczynski-Wiita potential. In a post-processing step, we evaluate our models for the energy deposition by nu-nubar annihilation around the accretion torus. Our models show that nu-nubar annihilation can yield the energy to account for weak, short gamma-ray bursts, if moderate beaming is involved. In fact, the barrier of the dense baryonic gas of the torus suggests that the low-density pair-photon-plasma is beamed as axial jets into a fraction 2 delta Omega/ (4 pi) between 1/100 and 1/10 of the sky, corresponding to opening half-angles of roughly ten to several tens of degrees. Thus gamma-burst energies of 10^{50}--10^{51} erg seem within the reach of our models (if the source is interpreted as radiating isotropically), corresponding to luminosities around 10^{51} erg/s for typical burst durations of 0.1--1 s. Gravitational capture of radiation by the black hole, redshift and ray bending do not reduce the jet energy significantly. Effects associated with the Kerr character of the rapidly rotating black hole, however, could increase the gamma-burst energy considerably, and effects due to magnetic fields might even be required to get the energies for long complex gamma-ray bursts.

We have coded a Boltzmann solver based on a finite difference scheme (S_N method) aiming at calculations of neutrino transport in type II supernovae. Close comparison between the Boltzmann solver and a Monte Carlo transport code has been made for realistic atmospheres of post bounce core models under the assumption of a static background. We have also investigated in detail the dependence of the results on the numbers of radial, angular, and energy grid points and the way to discretize the spatial advection term which is used in the Boltzmann solver. A general relativistic calculation has been done for one of the models. We find overall good agreement between the two methods. However, because of a relatively small number of angular grid points (which is inevitable due to limitations of the computation time) the Boltzmann solver tends to underestimate the flux factor and the Eddington factor outside the (mean) ``neutrinosphere'' where the angular distribution of the neutrinos becomes highly anisotropic. This fact suggests that one has to be cautious in applying the Boltzmann solver to a calculation of the neutrino heating in the hot-bubble region because it might tend to overestimate the local energy deposition rate. A comparison shows that this trend is opposite to the results obtained with a multi-group flux-limited diffusion approximation of neutrino transport. The accuracy of the Boltzmann solver can be considerably improved by using a variable angular mesh to increase the angular resolution in the semi-transparent regime.

Three-dimensional hydrodynamical simulations are presented for the direct head-on or off-center collision of two neutron stars, employing a basically Newtonian PPM code but including the emission of gravitational waves and their back-reaction on the hydrodynamical flow. A physical nuclear equation of state is used that allows us to follow the thermodynamical evolution of the stellar matter and to compute the emission of neutrinos. Predicted gravitational wave signals, luminosities and waveforms, are presented. The models are evaluated for their implications for gamma-ray burst scenarios. We find an extremely luminous outburst of neutrinos with a peak luminosity of more than 4E54 erg/s for several milliseconds. This leads to an efficiency of about 1% for the annihilation of neutrinos with antineutrinos, corresponding to an average energy deposition rate of more than 1E52 erg/s and a total energy of about 1E50 erg deposited in electron-positron pairs around the collision site within 10ms. Although these numbers seem very favorable for gamma-ray burst scenarios, the pollution of the $e^\pm$ pair-plasma cloud with nearly 0.1$M_{\odot}$ of dynamically ejected baryons is 5 orders of magnitude too large. Therefore the formation of a relativistically expanding fireball that leads to a gamma-ray burst powered by neutrino emission from colliding neutron stars is definitely ruled out.

Collection of "Extended Abstracts" of talks presented at the SFB-375 workshop NEUTRINO ASTROPHYSICS, Ringberg Castle, Tegernsee, Germany, 20-24 Oct. 1997. Topics include Solar Neutrinos, Supernovae, Gamma-Ray Bursts, High-Energy Neutrinos (atmospheric, astrophysical), Cosmology, and Future Prospects. The focus is on astrophysical and experimental/observational aspects of astrophysical neutrinos while particle theory and neutrino laboratory experiments are not covered.

The observation of neutrinos from Supernova 1987A has confirmed the theoretical conjecture that these particles play a crucial role during the collapse of the core of a massive star. Only one per cent of the energy they carry away from the newly formed neutron star may account for all the kinetic and electromagnetic energy responsible for the spectacular display of the supernova explosion. However, the neutrinos emitted from the collapsed stellar core at the center of the explosion couple so weakly to the surrounding matter that convective processes behind the supernova shock and/or inside the nascent neutron star might be required to increase the efficiency of the energy transfer to the stellar mantle and envelope. The conditions for a successful explosion by the neutrino-heating mechanism and the possible importance of convection in and around the neutron star are shortly discussed. Neutrino-driven explosions turn out to be very sensitive to the parameters describing the neutrino emission of the proto-neutron star and to the details of the dynamical processes in the collapsed stellar core. Therefore uniform explosions with a well defined energy seem unlikely and type-II supernova explosions do not offer promising perspectives for being useful as standard candles.

The observation of neutrinos from Supernova~1987A has confirmed the theoretical conjecture that these particles play a crucial role during the collapse of the core of a massive star. Only one per cent of the energy they carry away from the newly formed neutron star may account for all the kinetic and electromagnetic energy responsible for the spectacular display of the supernova explosion. However, the neutrinos emitted from the collapsed stellar core at the center of the explosion couple so weakly to the surrounding matter that convective processes behind the supernova shock and/or inside the nascent neutron star might be required to increase the efficiency of the energy transfer to the stellar mantle and envelope. The conditions for a successful explosion by the neutrino-heating mechanism and the possible importance of convection in and around the neutron star are shortly reviewed.

We present the results of numerical studies of supernova remnant evolution and their effects on galactic and globular cluster evolution. We show that parameters such as the density and the metallicity of the environment significantly influence the evolution of the remnant, and thus change its effects on the global environment (e.g., globular clusters, galaxies) as a source of thermal and kinetic energy. We conducted our studies using a one-dimensional hydrodynamics code, in which we implemented a metallicity dependent cooling function. Global time-dependent quantities such as the total kinetic and thermal energies and the radial extent are calculated for a grid of parameter sets. The quantities calculated are the total energy, the kinetic energy, the thermal energy, the radial extent, and the mass. We distinguished between the hot, rarefied bubble and the cold, dense shell, as those two phases are distinct in their roles in a gas-stellar system. We also present power-law fits to those quantities as a function of environmental parameters after the extensive cooling has ceased. The power-law fits enable simple incorporation of improved supernova energy input and matter redistribution (including the effect of the local conditions) in galactic/globular cluster models. Our results for the energetics of supernova remnants in the late stages of their expansion give total energies ranging from 9e49 to 3e50 ergs, with a typical case being 1e50 erg, depending on the surrounding environment. About 8.5e49 erg of this energy can be found in the form of kinetic energy. Supernovae play an important role in the evolution of the interstellar medium

We re-examine the very stringent limits on the axion mass based on the strength and duration of the neutrino signal from SN 1987A, in the light of new measurements of the axial-vector coupling strength of nucleons, possible suppression of axion emission due to many-body effects, and additional emission processes involving pions. The suppression of axion emission due to nucleon spin fluctuations induced by many-body effects degrades previous limits by a factor of about 2. Emission processes involving thermal pions can strengthen the limits by a factor of 3-4 within a perturbative treatment that neglects saturation of nucleon spin fluctuations. Inclusion of saturation effects, however, tends to make the limits less dependent on pion abundances. The resulting axion mass limit also depends on the precise couplings of the axion and ranges from 0.5x10**(-3) eV to 6x10**(-3) eV.

The dynamics, time evolution of the mass distribution, and gravitational wave signature of coalescing neutron stars described by polytropes are compared with three simulations published previously: (a) ``Run 2'' of Zhuge et al. (1994), (b) ``Model III'' of Shibata et al. (1992), and (c) ``Model A64'' of Ruffert et al. (1996). We aim at studying the differences due to the use of different numerical methods, different implementations of the gravitational wave backreaction, and different equations of state. Comparison (a) confronts the results of our grid-based PPM scheme with those from an SPH code. We found that due to the lower numerical viscosity of the PPM code, the post-merging oscillations and pulsations could be followed for a longer time and lead to larger secondary and tertiary maxima of the gravitational wave luminosity. In case (b) two grid based codes with the same backreaction formalism but differing hydrodynamic integrators and different numerical resolution are compared. Satisfactory agreement of the amplitude of the gravitational wave luminosity is established, although due to the different initial conditions a small time delay develops in the onset of the dynamical instability. In (c) we find that using a polytropic equation of state instead of the high-density equation of state of Lattimer & Swesty (1991) does not change the overall dynamical evolution of the merger and yields agreement of the gravitational wave signature to within 20% accuracy. However, differences of the structure and evolution of the outer layers of the neutron stars are present, which has important implications for questions like mass loss and disk formation during the merging of binary neutron stars.

The results of recent multi-dimensional simulations of type-II supernovae are reviewed. They show that convective instabilities in the collapsed stellar core might play an important role already during the first second after the formation of the supernova shock. Convectively unstable situations occur below and near the neutrinosphere as well as in the neutrino-heated region between the nascent neutron star and the supernova shock after the latter has stalled at a radius of typically 100-200 km. While convective overturn in the layer of neutrino energy deposition clearly helps the explosion to develop and potentially provides an explanation of strong mantle and envelope mixing, asphericities, and non-uniform 56Ni distribution observed in supernova SN 1987A, its presence and importance depends on the strength of the neutrino heating and thus on the size of the neutrino fluxes from the neutron star. Convection in the hot-bubble region can only develop if the growth timescale of the instabilities and the heating timescale are both shorter than the accretion timescale of the matter advected through the stagnant shock. For too small neutrino luminosities this requirement is not fulfilled and convective activity cannot develop, leading to very weak explosions or even fizzling models, just as in the one-dimensional situation. Convectively enhanced neutrino luminosities from the protoneutron star can therefore provide an essential condition for the explosion of the star. Very recent two-dimensional, self-consistent, general relativistic simulations of the cooling of a newly-formed neutron star demonstrate and confirm the possibility that Ledoux convection, driven by negative lepton number and

Two-dimensional hydrodynamical simulations of the deleptonization of a newly formed neutron star were performed. Driven by negative lepton fraction and entropy gradients, convection starts near the neutrinosphere about 20-30 ms after core bounce, but moves deeper into the protoneutron star, and after about one second the whole protoneutron star is convective. The deleptonization of the star proceeds much faster than in the corresponding spherically symmetrical model because the lepton flux and the neutrino luminosities increase by up to a factor of two. The convection below the neutrinosphere raises the neutrinospheric temperatures and mean energies of the emitted neutrinos by 10-20%. This can have important implications for the supernova explosion mechanism and changes the detectable neutrino signal from the Kelvin-Helmholtz cooling of the protoneutron star. In particular, the enhanced electron neutrino flux relative to the electron antineutrino flux during the early post-bounce evolution might solve the overproduction problem of certain elements in the neutrino-heated ejecta in models of type-II supernova explosions.

Three-dimensional hydrodynamical, Newtonian calculations of the coalescence of equal-mass binary neutron stars are performed, including a physical high-density equation of state and a treatment of the neutrino emission of the heated matter. The total neutrino luminosity climbs to a maximum value of 1--$1.5\cdot 10^{53}$~erg/s of which 90--95\% originate from the toroidal gas cloud surrounding the very dense core formed after the merging. When the neutrino luminosities are highest, $\nu\bar\nu$-annihilation deposits about 0.2--0.3\% of the emitted neutrino energy in the immediate neighborhood of the merger, and the maximum integral energy deposition rate is 3--$4\cdot 10^{50}$~erg/s. Since the $3\,M_{\odot}$ core of the merged object will most likely collapse into a black hole within milliseconds, the energy that can be pumped into a pair-photon fireball is insufficient by a factor of about 1000 to explain $\gamma$-ray bursts at cosmological distances with an energy of the order of $10^{51}/(4\pi)$~erg/steradian. Analytical estimates show that the additional energy provided by the annihilation of $\nu\bar\nu$ pairs emitted from a possible accretion torus of $\sim 0.1\,M_{\odot}$ around the central black hole is still more than a factor of 10 too small, unless focussing of the fireball into a jet-like expansion plays an important role. About $10^{-4}$--$10^{-3}$~$M_\odot$ of material lost during the neutron star merging and swept out from the system in a neutrino-driven wind might be a site for nucleosythesis. Aspects of a possible r-processing in these ejecta are discussed.

In the hot and dense medium of a supernova (SN) core, the nucleon spins fluctuate so fast that the axial-vector neutrino opacity and the axion emissivity are expected to be significantly modified. Axions with $m_a\alt10^{-2}\,{\rm eV}$ are not excluded by SN~1987A. A substantial transfer of energy in neutrino-nucleon ($\nu N$) collisions is enabled which may alter the spectra of SN neutrinos relative to calculations where energy-conserving $\nu N$ collisions had been assumed near the neutrinosphere.

We perform three-dimensional hydrodynamical simulations of the coalescence of binary neutron stars. We include the emission and backreaction of gravitational waves into the Newtonian ``Piecewise Parabolic Method''. The use of the physical equation of state of Lattimer & Swesty (1991) allows us to take into account the production of neutrinos. We evaluate our models for the efficiency of neutrino-antineutrino annihilation in the surroundings of the coalescing neutron stars. The corresponding energy deposition prior to and during merging turns out to be 2--3 orders of magnitude too small to power a typical gamma-ray burst with an energy output of about 10^{51}/(4 pi) erg/sterad at cosmological distances. Analytical estimates of the subsequent evolution of the disk which possibly surrounds the central black hole show that even under the most favorable conditions the energy in a fireball of electron-positron pairs and photons that is created by neutrino-antineutrino annihilation falls short by at least an order of magnitude.

We investigate the dynamics and evolution of coalescing neutron stars. Although the code (Piecewise Parabolic Method) is purely Newtonian, we do include the emission of gravitational waves and their backreaction on the hydrodynamic flow. The properties of neutron star matter are described by the physical equation of state of Lattimer \& Swesty (1991). Energy loss by all types of neutrinos and changes of the electron fraction due to the emission of electron neutrinos and antineutrinos are taken into account by an elaborate ``neutrino leakage scheme''. We simulate the coalescence of two identical, cool neutron stars with a baryonic mass of $\approx\!1.6\,M_\odot$ and a radius of $\approx\!15$~km and with an initial center-to-center distance of 42~km. The initial distributions of density and electron concentration are given from a model of a cold neutron star in hydrostatic equilibrium (central temperature about $8\,{\rm MeV}$). We investigate three cases which differ by the initial velocity distribution in the neutron stars, representing different cases of the neutron star spins relative to the direction of the orbital angular momentum vector. Within about 1~ms the neutron stars merge into a rapidly spinning ($P_{\rm spin}\approx 1$~ms), high-density body ($\rho\approx 10^{14}$~g/cm$^3$) with a surrounding thick disk of material with densities $\rho\approx 10^{10}-10^{12}$~g/cm$^3$ and orbital velocities of~0.3--0.5~c. In this work we evaluate the models in detail with respect to the gravitational wave emission using the quadrupole approximation. In a forthcoming paper we will concentrate on the neutrino emission and implications for gamma-ray bursters. A maximum luminosity in excess of $10^{55}$~erg/s is reached for about 1~ms.

The neutrino signal from SN~1987A provides an excellent opportunity to constrain physical theories for matter at extreme conditions and properties of particles that are produced in supernova (SN) cores. Phase transitions in the supranuclear equation of state (EOS) may change the cooling history of the collapsed stellar core and may even lead to a collapse of the protoneutron star to a black hole on a time scale of several seconds. The implications of such a scenario for the characteristics of the neutrino emission are tested. Moreover, the consequences of a reduction of the neutrino opacities by a possible suppression of the axial-vector interaction in a dense medium (Raffelt and Seckel 1994) are investigated. Since the strongest constraints from SN~1987A for the considered theories rely on the signal duration marked by the last three events in Kamiokande, we also discuss attempts to explain these events different from the standard protoneutron star cooling scenario. In particular, attention is focussed on the scenario of a rapid cooling of the protoneutron star within the time interval of the first 8 Kamiokande events or the IMB signal which is only about 2--$5\,{\rm s}$. Such a short neutrino emission could be the consequence of low neutrino opacities or convectively enhanced neutrino transport. Since the driving force of neutrino heating around the protoneutron star diminishes, the pressure support of the overlying material decreases and some matter might fall back to the neutron star. However, analytical estimates and numerical simulations indicate that the neutrino emission associated with the accretion of this material is hardly energetic enough to account for the last bunch of neutrinos at Kamiokande.

Muon and tau neutrinos ($\nu_x$) interact with protoneutron star matter only via neutral currents and exchange energy with the stellar gas predominantly by neutrino-electron scattering and neutrino-pair processes. In contrast, electron neutrinos and antineutrinos ($\nu_e$ and $\bar\nu_e$) are frequently absorbed and produced in charged-current mediated reactions with nucleons. Therefore the emergent $\nu_e$ and $\bar\nu_e$ originate from layers with lower temperatures further out in the star and are emitted with much lower characteristic spectral temperatures. In addition, a major contribution to the $\nu_e$ and $\bar\nu_e$ opacities is due to absorptions, while the opacity of $\nu_x$ is strongly dominated by scattering reactions with nucleons and nuclei in which the $\nu_x$ energy is (essentially) conserved. Therefore the $\nu_x$ distribution is nearly isotropic when $\nu_x$ decouple energetically and their outward diffusion is slowed down. In a generalized form to include this effect, the Stefan-Boltzmann Radiation Law can account for both the facts that $\nu_e$ ($\bar\nu_e$) and $\nu_x$ emerge from the star with similar luminosities but with very different characteristic spectral temperatures. Simple analytical expressions to estimate the effect are given. If, as recently argued, even at densities significantly below nuclear matter density neutral-current scatterings were associated with considerable energy transfer between neutrino and target particle, one might expect spectral temperatures of $\nu_x$ much closer to those of $\nu_e$ and $\bar\nu_e$. This is of relevance for the detection of neutrino signals from supernovae.

We investigate the dynamics and evolution of coalescing neutron stars. The three-dimensional Newtonian equations of hydrodynamics are integrated by the `Piecewise Parabolic Method' However, we do include the effects of the emission of gravitational waves on the hydrodynamics. The properties of neutron star matter are described by the equation of state of Lattimer & Swesty. In addition to the fundamental hydrodynamic quantities, density, momentum, and energy, we follow the time evolution of the electron density in the stellar gas. Energy losses and changes of the electron abundance due to the emission of neutrinos are taken into account by an elaborate ``neutrino leakage scheme'', which employs a careful calculation of the lepton number and energy source terms of all neutrino types. The grid is Cartesian and equidistant with a resolution of 64**3 or 128**3, which allows us to calculate the self-gravity via fast Fourier transforms.

One- and two-dimensional hydrodynamical simulations of neutrino-driven supernova explosions are discussed. The simulations cover the phase between the stagnation of the prompt shock and about one second after core bounce. Systematic variation of the neutrino fluxes from the neutrino sphere shows that the explosion energy, explosion time scale, initial mass of the protoneutron star, and explosive nucleosynthesis of iron-group elements depend sensitively on the strength of the neutrino heating during the first few 100 ms after shock formation. Convective overturn in the neutrino-heated region behind the shock is a crucial help for the explosion only in a narrow window of neutrino luminosities. Here powerful explosions can be obtained only in the multi-dimensional case. For higher core-neutrino fluxes also spherically symmetrical models yield energetic explosions, while for lower luminosities even with convection no strong explosions occur.

Axial-vector interactions of neutrinos and axions with a medium are based on structure functions which cannot be calculated in simple perturbative terms for nuclear densities. We use the SN~1987A neutrino signal duration to estimate the range of allowed neutrino opacities of a supernova (SN) core. We perform numerical simulations of the evolution from the hot, lepton-rich post-collapse stage to the cold, deleptonized neutron star and compare the predicted neutrino signal characteristics with the SN~1987A data. Besides the neutrino opacities we vary the mass, temperature profile, and equation of state of the initial model. Within plausible limits of these quantities the neutrino opacity cannot be much smaller than its ``standard'' value. This limits the width of the spin-density structure function to much smaller values than indicated by a perturbative calculation, implying that bremsstrahlung processes such as axion emission probably saturate at around 10\% nuclear density. A substantial suppression at yet larger densities that might have been expected does not seem to occur.