Uncertainties in our knowledge of the properties of dense matter near and above nuclear saturation density are among the main sources of variations in multimessenger signatures predicted for core-collapse supernovae (CCSNe) and the properties of neutron stars (NSs). We construct 97 new finite-temperature equations of state (EOSs) of dense matter that obey current experimental, observational, and theoretical constraints and discuss how systematic variations in the EOS parameters affect the properties of cold nonrotating NSs and the core collapse of a 20–M ⊙ progenitor star. The core collapse of the 20–M ⊙ progenitor star is simulated in spherical symmetry using the general-relativistic radiation-hydrodynamics code GR1D where neutrino interactions are computed for each EOS using the NuLib library. We conclude that the effective mass of nucleons at densities above nuclear saturation density is the largest source of uncertainty in the CCSN neutrino signal and dynamics even though it plays a subdominant role in most properties of cold NS matter. Meanwhile, changes in other observables affect the properties of cold NSs, while having little effect in CCSNe. To strengthen our conclusions, we perform six octant three-dimensional CCSN simulations varying the effective mass of nucleons at nuclear saturation density. Here, we conclude that neutrino heating and, thus, the likelihood of explosion is significantly increased for EOSs where the effective mass of nucleons at nuclear saturation density is large.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

DOE OSTI2017

Neutrino-nucleon scattering in supernova matter from the virial expansion

Here we extend our virial approach to study the neutral-current neutrino response of nuclear matter at low densities. In the long-wavelength limit, the virial expansion makes model-independent predictions for neutrino-nucleon scattering rates and the density S V and spin S A responses. We find that S A is significantly reduced from one even at low densities. We provide a simple fit $S^f_A$(n,T,Y p ) of the axial response as a function of density n, temperature T, and proton fraction Y p , which can be incorporated into supernova simulations in a straightforward manner. This fit reproduces our virial results at low densities and the Burrows and Sawyer random-phase approximation (RPA) model calculations at high densities. Preliminary one-dimensional supernova simulations suggest that the virial reduction in the axial response may enhance neutrino heating rates in the gain region during the accretion phase of a core-collapse supernovae.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

DOE OSTI2013

BLACK HOLE-NEUTRON STAR MERGERS WITH A HOT NUCLEAR EQUATION OF STATE: OUTFLOW AND NEUTRINO-COOLED DISK FOR A LOW-MASS, HIGH-SPIN CASE

Neutrino emission significantly affects the evolution of the accretion tori formed in black hole-neutron star mergers. It removes energy from the disk, alters its composition, and provides a potential power source for a gamma-ray burst. To study these effects, simulations in general relativity with a hot microphysical equation of state (EOS) and neutrino feedback are needed. We present the first such simulation, using a neutrino leakage scheme for cooling to capture the most essential effects and considering a moderate mass (1.4 M{sub ☉} neutron star, 5.6 M{sub ☉} black hole), high-spin (black hole J/M {sup 2} = 0.9) system with the K{sub 0} = 220 MeV Lattimer-Swesty EOS. We find that about 0.08 M{sub ☉} of nuclear matter is ejected from the system, while another 0.3 M{sub ☉} forms a hot, compact accretion disk. The primary effects of the escaping neutrinos are (1) to make the disk much denser and more compact, (2) to cause the average electron fraction Y{sub e} of the disk to rise to about 0.2 and then gradually decrease again, and (3) to gradually cool the disk. The disk is initially hot (T ∼ 6 MeV) and luminous in neutrinos (L{sub ν} ∼ 10{sup 54} erg s{sup –1}), but the neutrino luminosity decreases by an order of magnitude over 50 ms of post-merger evolution.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2013

THE PROGENITOR DEPENDENCE OF THE PRE-EXPLOSION NEUTRINO EMISSION IN CORE-COLLAPSE SUPERNOVAE

We perform spherically symmetric general-relativistic simulations of core collapse and the postbounce pre-explosion phase in 32 presupernova stellar models of solar metallicity with zero-age main-sequence masses of 12-120 M {sub Sun }. Using energy-dependent three-species neutrino transport in the two-moment approximation with an analytic closure, we show that the emitted neutrino luminosities and spectra follow very systematic trends that are correlated with the compactness ({approx}M/R) of the progenitor star's inner regions via the accretion rate in the pre-explosion phase. We find that these qualitative trends depend only weakly on the nuclear equation of state (EOS), but quantitative observational statements will require independent constraints on the EOS and the rotation rate of the core as well as a more complete understanding of neutrino oscillations. We investigate the simulated response of water Cherenkov detectors to the electron antineutrino fluxes from our models and find that the large statistics of a galactic core collapse event may allow robust conclusions on the inner structure of the progenitor star.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2012

A NEW MONTE CARLO METHOD FOR TIME-DEPENDENT NEUTRINO RADIATION TRANSPORT

Monte Carlo approaches to radiation transport have several attractive properties such as simplicity of implementation, high accuracy, and good parallel scaling. Moreover, Monte Carlo methods can handle complicated geometries and are relatively easy to extend to multiple spatial dimensions, which makes them potentially interesting in modeling complex multi-dimensional astrophysical phenomena such as core-collapse supernovae. The aim of this paper is to explore Monte Carlo methods for modeling neutrino transport in core-collapse supernovae. We generalize the Implicit Monte Carlo photon transport scheme of Fleck and Cummings and gray discrete-diffusion scheme of Densmore et al. to energy-, time-, and velocity-dependent neutrino transport. Using our 1D spherically-symmetric implementation, we show that, similar to the photon transport case, the implicit scheme enables significantly larger timesteps compared with explicit time discretization, without sacrificing accuracy, while the discrete-diffusion method leads to significant speed-ups at high optical depth. Our results suggest that a combination of spectral, velocity-dependent, Implicit Monte Carlo and discrete-diffusion Monte Carlo methods represents a robust approach for use in neutrino transport calculations in core-collapse supernovae. Our velocity-dependent scheme can easily be adapted to photon transport.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2012

THE ARDUOUS JOURNEY TO BLACK HOLE FORMATION IN POTENTIAL GAMMA-RAY BURST PROGENITORS

We present a quantitative study on the properties at death of fast-rotating massive stars evolved at low-metallicity-objects that are proposed as likely progenitors of long-duration {gamma}-ray bursts (LGRBs). We perform one-dimensional+rotation stellar-collapse simulations on the progenitor models of Woosley and Heger, and critically assess their potential for the formation of a black hole and a Keplerian disk (namely, a collapsar) or a proto-magnetar. We note that theoretical uncertainties in the treatment of magnetic fields and the approximate handling of rotation compromise the accuracy of stellar-evolution models. We find that only the fastest rotating progenitors achieve sufficient compactness for black hole formation while the bulk of models possess a core density structure typical of garden-variety core-collapse supernova (SN) progenitors evolved without rotation and at solar metallicity. Of the models that do have sufficient compactness for black hole formation, most of them also retain a large amount of angular momentum in the core, making them prone to a magneto-rotational explosion, therefore preferentially leaving behind a proto-magnetar. A large progenitor angular-momentum budget is often the sole criterion invoked in the community today to assess the suitability for producing a collapsar. This simplification ignores equally important considerations such as the core compactness, which conditions black hole formation, the core angular momentum, which may foster a magneto-rotational explosion preventing black hole formation, or the metallicity and the residual envelope mass which must be compatible with inferences from observed LGRB/SNe. Our study suggests that black hole formation is non-trivial, that there is room for accommodating both collapsars and proto-magnetars as LGRB progenitors, although proto-magnetars seem much more easily produced by current stellar-evolutionary models.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2011

BLACK HOLE FORMATION IN FAILING CORE-COLLAPSE SUPERNOVAE

We present results of a systematic study of failing core-collapse supernovae and the formation of stellar-mass black holes (BHs). Using our open-source general-relativistic 1.5D code GR1D equipped with a three-species neutrino leakage/heating scheme and over 100 presupernova models, we study the effects of the choice of nuclear equation of state (EOS), zero-age main sequence (ZAMS) mass and metallicity, rotation, and mass-loss prescription on BH formation. We find that the outcome, for a given EOS, can be estimated, to first order, by a single parameter, the compactness of the stellar core at bounce. By comparing protoneutron star (PNS) structure at the onset of gravitational instability with solutions of the Tolman-Oppenheimer-Volkof equations, we find that thermal pressure support in the outer PNS core is responsible for raising the maximum PNS mass by up to 25% above the cold NS value. By artificially increasing neutrino heating, we find the critical neutrino heating efficiency required for exploding a given progenitor structure and connect these findings with ZAMS conditions, establishing, albeit approximately, for the first time based on actual collapse simulations, the mapping between ZAMS parameters and the outcome of core collapse. We also study the effect of progenitor rotation and find that the dimensionless spin of nascent BHs may be robustly limited below a* = Jc/GM{sup 2} = 1 by the appearance of nonaxisymmetric rotational instabilities.

79 ASTRONOMY AND ASTROPHYSICS↗