Monte Carlo Study of Supernova Neutrino Spectra Formation

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.Comment: 33 pages, 16 eps-figs, submitted to ApJ. Sections added: weak magnetism, discussion of different analytic fits to the spectra and detailed spectral shap

Testing Approximations of Thermal Effects in Neutron Star Merger Simulations

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.Comment: 10 pages, 6 figures, 9 eps files; revised with minor additions due to referee comments; accepted by Phys.Rev.

Ledoux-Convection in Protoneutron Stars --- a Clue to Supernova Nucleosynthesis?

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.Comment: 17 pages, LaTeX, 8 postscript figures, uses epsf.sty. To appear in ApJ 473 (Letters), 1996 December 1

Growth of Perturbation in Gravitational Collapse and Accretion

When a self-gravitating spherical gas cloud collapses or accretes onto a central mass, the inner region of the cloud develops a density profile $\rho\propto r^{-3/2}$ and the velocity approaches free-fall. We show that in this region, nonspherical perturbations grow with decreasing radius. In the linear regime, the tangential velocity perturbation increases as $r^{-1}$, while the Lagrangian density perturbation, $\Delta\rho/\rho$, grows as $r^{-1/2}$. Faster growth occurs if the central collapsed object maintains a finite multiple moment, in which case $\Delta\rho/\rho$ increases as $r^{-l}$, where $l$ specifies the angular degree of the perturbation. These scaling relations are different from those obtained for the collapse of a homogeneous cloud. Our numerical calculations indicate that nonspherical perturbations are damped in the subsonic region, and that they grow and approach the asymptotic scalings in the supersonic region. The implications of our results to asymmetric supernova collapse and to black hole accretion are briefly discussed.Comment: 23 pages including 6 ps figures; Minor changes and update; To appear in ApJ, 200

Electron Neutrino Pair Annihilation: A New Source for Muon and Tau Neutrinos in 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.Comment: 11 pages, 9 eps-figs, submitted to Ap

Core-collapse supernova simulations: Variations of the input physics

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.Comment: 6 pages, 6 ps figures (in color), to appear in W. Hillebrandt and E. Mueller, eds., Proceedings of the 11th Workshop on "Nuclear Astrophysics" held at Ringberg Castle, February 11-16, 200

Two-Dimensional Hydrodynamic Core-Collapse Supernova Simulations with Spectral Neutrino Transport II. Models for Different Progenitor Stars

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)Comment: 42 pages, 44 figures; revised according to referee comments; accepted to Astronomy & Astrophysic

Core Collapse and Then? The Route to Massive Star Explosions

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.Comment: 14 pages; 3 figures. Invited Review, in: ``From Twilight to Highlight: The Physics of Supernovae'', Eds. W. Hillebrandt and B. Leibundgut, Springer Series ``ESO Astrophysics Symposia'', Berli

Exploiting the neutronization burst of a galactic supernova

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.Comment: 14 pages, 17 eps figures, revtex4 style, minor comments adde

Short Gamma-Ray Bursts from Binary Neutron Star Mergers

We present the results from new relativistic hydrodynamic simulations of binary neutron star mergers using realistic non-zero temperature equations of state. We vary several unknown parameters in the system such as the neutron star (NS) masses, their spins and the nuclear equation of state. The results are then investigated with special focus on the post-merger torus-remnant system. Observational implications on the Gamma-ray burst (GRB) energetics are discussed and compared with recent observations