A first-order quantum chromodynamics (QCD) phase transition (PT) may take place in the protocompact star (PCS) produced by a core-collapse supernova (CCSN). In this work, we study the consequences of such a PT in a non-rotating CCSN with axisymmetric hydrodynamic simulations. We find that the PT leads to the collapse of the PCS and results in a loud burst of gravitational waves (GWs). The amplitude of this GW burst is $\sim30$ times larger than the post-bounce GW signal normally found for non-rotating CCSN. It shows a broad peak at high frequencies ($\sim2500-4000$ Hz) in the spectrum, has a duration of $\lesssim5 {\rm ms}$, and carries $\sim3$ orders of magnitude more energy than the other episodes. Also, the peak frequency of the PCS oscillation increases dramatically after the PT-induced collapse. In addition to a second neutrino burst, the GW signal, if detected by the ground-based GW detectors, is decisive evidence of the first-order QCD PT inside CCSNe and provides key information about the structure and dynamics of the PCS.

With the advent of modern neutrino and gravitational wave detectors, the promise of multi-messenger detections of the next galactic core-collapse supernova has become very real. Such detections will give insight into the core-collapse supernova mechanism, the structure of the progenitor star, and may resolve longstanding questions in fundamental physics. In order to properly interpret these detections, a thorough understanding of the landscape of possible core-collapse supernova events, and their multi-messenger signals, is needed. We present detailed predictions of neutrino and gravitational wave signals from 1D simulations of stellar core collapse, spanning the landscape of core-collapse progenitors from $9-120\,\mathrm{M}_{\odot}$. In order to achieve explosions in 1D, we use the STIR model, which includes the effects of turbulence and convection in 1D supernova simulations to mimic the 3D explosion mechanism. We study the gravitational wave emission from the 1D simulations using an astroseismology analysis of the proto-neutron star. We find that the neutrino and gravitational wave signals are strongly correlated with the structure of the progenitor star and remnant compact object. Using these correlations, future detections of the first few seconds of neutrino and gravitational wave emission from a galactic core-collapse supernova may be able to provide constraints on stellar evolution independent of pre-explosion imaging and the mass of the compact object remnant prior to fallback accretion.

The core collapse of a massive star results in the formation of a proto-neutron star (PNS). If enough material is accreted onto a PNS it will become gravitationally unstable and further collapse into a black-hole (BH). We perform a systematic study of failing core-collapse supernovae in spherical symmetry for a wide range of presupernova progenitor stars and equations of state (EOSs) of nuclear matter. We analyze how variations in progenitor structure and the EOS of dense matter above nuclear saturation density affect the PNS evolution and subsequent BH formation. Comparisons of core-collapse for a given progenitor star and different EOSs show that the path traced by the PNS in mass-entropy phase space $M_{\mathrm{grav}}^{\mathrm{PNS}}-\tilde{s}$ is well correlated with the progenitor compactness and almost EOS independent, apart from the final endpoint. Furthermore, BH formation occurs, to a very good approximation, soon after the PNS overcomes the maximum \textit{gravitational} mass supported by a hot NS with constant entropy equal to $\tilde{s}$. These results show a path to constraining the temperature dependence of the EOS through the detection of neutrinos from a failed galactic supernova.

The core-collapse supernova (CCSN) mechanism is fundamentally three-dimensional with instabilities, convection, and turbulence playing crucial roles in aiding neutrino-driven explosions. Simulations of CCNSe including accurate treatments of neutrino transport and sufficient resolution to capture key instabilities remain amongst the most expensive numerical simulations in astrophysics, prohibiting large parameter studies in 2D and 3D. Studies spanning a large swath of the incredibly varied initial conditions of CCSNe are possible in 1D, though such simulations must be artificially driven to explode. We present a new method for including the most important effects of convection and turbulence in 1D simulations of neutrino-driven CCSNe, called Supernova Turbulence In Reduced-dimensionality, or STIR. Our new approach includes crucial terms resulting from the turbulent and convective motions of the flow. We estimate the strength of convection and turbulence using a modified mixing length theory (MLT) approach introducing a few free parameters to the model which are fit to the results of 3D simulations. For sufficiently large values of the mixing length parameter, turbulence-aided neutrino-driven explosions are obtained. We compare the results of STIR to high-fidelity 3D simulations and perform a parameter study of CCSN explosion using 200 solar-metallicity progenitor models from 9 to 120 $M_\odot$. We find that STIR is a better predictor of which models will explode in multidimensional simulations than other methods of driving explosions in 1D. We also present a preliminary investigation of predicted observable characteristics of the CCSN population from STIR, such as the distributions of explosion energies and remnant masses.

We investigate correlated gravitational wave and neutrino signals from rotating core-collapse supernovae with simulations. Using an improved mode identification procedure based on mode function matching, we show that a linear quadrupolar mode of the core produces a dual imprint on gravitational waves and neutrinos in the early post-bounce phase of the supernova. The angular harmonics of the neutrino emission are consistent with the mode energy around the neutrinospheres, which points to a mechanism for the imprint on neutrinos. Thus, neutrinos carry information about the mode amplitude in the outer region of the core, whereas gravitational waves probe deeper in. We also find that the best-fit mode function has a frequency bounded above by $\sim 420$ Hz, and yet the mode's frequency in our simulations is $\sim 15\%$ higher, due to the use of Newtonian hydrodynamics and a widely used pseudo-Newtonian gravity approximation. This overestimation is particularly important for the analysis of gravitational wave detectability and asteroseismology, pointing to limitations of pseudo-Newtonian approaches for these purposes, possibly even resulting in excitation of incorrect modes. In addition, mode frequency matching (as opposed to mode function matching) could be resulting in mode misidentification in recent work. Lastly, we evaluate the prospects of a multimessenger detection of the mode using current technology. The detection of the imprint on neutrinos is most challenging, with a maximum detection distance of $\sim\!1$ kpc using the IceCube Neutrino Observatory. The maximum distance for detecting the complementary gravitational wave imprint is $\sim\!5$ kpc using Advanced LIGO at design sensitivity.

We present a new, multi-dimensional implementation of the Advanced Spectral Leakage (ASL) scheme with the purpose of modelling neutrino-matter interactions in neutron star mergers. A major challenge is the neutrino absorption in the semi-transparent regime, which is responsible for driving winds from the merger remnant. The composition of such winds is crucial in the understanding of the electromagnetic emission in the recently observed macronova following GW170817. Compared to the original version, we introduce an optical-depth-dependent flux factor to model the average angle of neutrino propagation, and a modulation that accounts for flux anisotropies in non-spherical geometries. We scrutinise our approach by first comparing the new scheme against the original one for a spherically symmetric core-collapse supernova snapshot, both in 1D and in 3D, and additionally against a two-moment (M1) scheme as implemented in 1D into the code GR1D. The luminosities and mean energies agree to a few percents in most tests. Finally, for the case of a binary merger remnant snapshot we compare the new ASL scheme with the M1 scheme that is implemented in the Eulerian adaptive mesh refinement code FLASH. We find that the neutrino absorption distribution in the semi-transparent regime is overall well reproduced. Both approaches agree to within $\lesssim 15\%$ for the average energies and to better than $\sim 35 \%$ in the total luminosities.

We explore the influence of progenitor mass and rotation on the gravitational-wave (GW) emission from core-collapse supernovae, during the postbounce, preexplosion, accretion-phase. We present the results from 15 two-dimensional (2D) neutrino radiation-hydrodynamic simulations from initial stellar collapse to $\sim$300 ms after core bounce. We examine the features of the GW signals for four zero-age main sequence (ZAMS) progenitor masses ranging from 12 $M_\odot$ to 60 $M_\odot$ and four core rotation rates from 0 to 3 rad s$^{-1}$. We find that GW strain immediately around core bounce is fairly independent of ZAMS mass and---consistent with previous findings---that it is more heavily dependent on the core angular momentum. At later times, all nonrotating progenitors exhibit loud GW emission, which we attribute to vibrational g-modes of the protoneutron star (PNS) excited by convection in the postshock layer and the standing accretion shock instability (SASI). We find that increasing rotation rates results in muting of the accretion-phase GW signal due to centrifugal effects that inhibit convection in the postshock region, quench the SASI, and slow the rate at which the PNS peak vibrational frequency increases. Additionally, we verify the efficacy of our approximate general relativistic (GR) effective potential treatment of gravity by comparing our core bounce GW strains with the recent 2D GR results of other groups.

Uncertainties in our knowledge of the properties of dense matter near and above nuclear saturation density are among the main sources of variations in multi-messenger 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_\odot$ progenitor star. The core collapse of the $20\,M_\odot$ 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. 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.

As of today, we have directly detected exactly one source in both gravitational waves (GWs) and electromagnetic (EM) radiation, the binary neutron star merger GW170817, its associated gamma-ray burst GRB170817A, and the subsequent kilonova SSS17a/AT 2017gfo. Within ten years, we will detect hundreds of events, including new classes of events such as neutron-star-black-hole mergers, core-collapse supernovae, and almost certainly something completely unexpected. As we build this sample, we will explore exotic astrophysical topics ranging from nucleosynthesis, stellar evolution, general relativity, high-energy astrophysics, nuclear matter, to cosmology. The discovery potential is extraordinary, and investments in this area will yield major scientific breakthroughs. Here we outline some of the most exciting scientific questions that can be answered by combining GW and EM observations.

We perform hydrodynamical simulations of the interaction between supernova (SN) ejecta and circumstellar medium (CSM) for SN 1993J and SN 2011dh, and calculate the radio and X-ray emissions expected from the shocked gas at late epochs ($t$). Considering the ejecta structure from multi-group radiation hydrodynamics simulation, we find that the observed rapid drop in radio and X-ray light curves of SN 1993J at $t>$3000 days can be due to a change in the mass-loss rate ($\dot M$) around $\sim$6500 years prior to the explosion of the SN. The exact epoch scales inversely with the assumed wind velocity of $v_{\rm w}=10~ km~s^{-1}$. The progenitor of this SN very likely belonged to a binary system, where, during its evolution, the primary had transferred material to the secondary. It is argued in the paper that the change in $\dot M$ can happen because of a change in the mass accretion efficiency ($\eta$) of the companion star. It is possible that before $\sim6500~(v_{\rm w}/10~km~s^{-1})^{-1}$years prior to the explosion, $\eta$ was high, thus the CSM was tenuous, which causes the late time downturn in fluxes. In the case of SN 2011dh, the late time evolution is found to be consistent with a wind medium with $\dot M/v_{\rm w}=4\times10^{-6 }~M_{\odot}~ yr^{-1}/10 ~{km ~s^{-1}}$. It is difficult from our analysis to predict whether the progenitor of this SN had a binary companion, however, if future observations show similar decrease in radio and X-ray fluxes, then that would give strong support to a scenario where both SNe had undergone similar kind of binary evolution before explosion.

Neutrinos play a crucial role in the core-collapse supernova (CCSN) explosion mechanism. The requirement of accurately calculating the transport of neutrinos makes simulations of the CCSN mechanism extremely challenging and computationally expensive. Historically, this stiff challenge has been met by making approximations to the full transport equation. In this work, we compare CCSN simulations in one- and two-dimensions with three approximate neutrino transport schemes, each implemented in the {\tt FLASH} simulation framework. We compare a two moment M1 scheme with an analytic closure (M1), the Isotropic Diffusion Source Approximation (IDSA), and the Advanced Spectral Leakage (ASL) method. We identify and discuss the advantages and disadvantages of each scheme. For each approximate transport scheme, we use identical grid setups, hydrodynamics, and gravity solvers to investigate the transport effects on supernova shock dynamics and neutrino quantities. We find that the transport scheme has a small effect on the evolution of protoneutron star (PNS) radius, PNS mass, and the mass accretion rate. The neutrino luminosities, mean energies, and shock radii have a $\sim$ 10-20\% quantitative difference but the overall qualitative trends are fairly consistent between all three approximations. We find larger differences in the gain region properties, including the gain region mass and the net heating rate in the gain region, as well as the strength of PNS convection in the core. We investigate the progenitor, nuclear equation of state, and stochastic perturbation dependence of our simulations and find similar magnitudes of impact on key quantities. We also compare the computational expense of the various approximations.

In July 2018 an FRIB Theory Alliance program was held on the implications of GW170817 and its associated kilonova for r-process nucleosynthesis. Topics of discussion included the astrophysical and nuclear physics uncertainties in the interpretation of the GW170817 kilonova, what we can learn about the astrophysical site or sites of the r process from this event, and the advances in nuclear experiment and theory most crucial to pursue in light of the new data. Here we compile a selection of scientific contributions to the workshop, broadly representative of progress in r-process studies since the GW170817 event.

The details of the physical mechanism that drives core-collapse supernovae (CCSNe) remain uncertain. While there is an emerging consensus on the qualitative outcome of detailed CCSN mechanism simulations in 2D, only recently have high-fidelity 3D simulations become possible. Here we present the results of an extensive set of 3D CCSN simulations using high-fidelity multidimensional neutrino transport, high-resolution hydrodynamics, and approximate general relativistic gravity. We employ a state-of-the-art 20 solar mass progenitor generated using the Modules for Experiments in Stellar Astrophysics (MESA; Farmer et al. (2016) Paxton et al. (2011, 2013, 2015, 2018) and the SFHo equation of state of Steiner et al. (2013). While none of our 3D CCSN simulations explode within ~500ms after core bounce, we find that the presence of large scale aspherical motion in the Si and O shells surrounding the collapsing iron core aid shock expansion and bring the models closer to the threshold of explosion. We also find some dependence on resolution and geometry (octant vs. full 4$\pi$). As has been noted in other recent works, we find that the post-shock turbulence plays an important role in determining the overall dynamical evolution of our simulations. We find a strong standing accretion shock instability (SASI) that develops at late times during the shock recession epoch. The SASI aids in transient shock expansion phases, but is not enough to result in shock revival. We also report that for a subset of our simulations, we find conclusive evidence for the LESA first reported in Tamborra et al. (2014), but until now, not confirmed by other simulation codes. Both the progenitor asphericities and the SASI-induced transient shock expansion phases generate transient gravitational waves and neutrino signal modulations via perturbations of the protoneutron star by turbulent motions. (abridged)

We present a covariant ray tracing algorithm for computing high-resolution neutrino distributions in general relativistic numerical spacetimes with hydrodynamical sources. Our formulation treats the very important effect of elastic scattering of neutrinos off of nuclei and nucleons (changing the neutrino's direction but not energy) by incorporating estimates of the background neutrino fields. Background fields provide information about the spectra and intensities of the neutrinos scattered into each ray. These background fields may be taken from a low-order moment simulation or be ignored, in which case the method reduces to a standard state-of-the-art ray tracing formulation. The method handles radiation in regimes spanning optically thick to optically thin. We test the new code, highlight its strengths and weaknesses, and apply it to a simulation of a neutron star merger to compute neutrino fluxes and spectra, and to demonstrate a neutrino flavor oscillation calculation. In that environment, we find qualitatively different fluxes, spectra, and oscillation behaviors when elastic scattering is included.

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.

Each of the potential signals from a black hole-neutron star merger should contain an imprint of the neutron star equation of state: gravitational waves via its effect on tidal disruption, the kilonova via its effect on the ejecta, and the gamma ray burst via its effect on the remnant disk. These effects have been studied by numerical simulations and quantified by semi-analytic formulae. However, most of the simulations on which these formulae are based use equations of state without finite temperature and composition-dependent nuclear physics. In this paper, we simulate black hole-neutron star mergers varying both the neutron star mass and the equation of state, using three finite-temperature nuclear models of varying stiffness. Our simulations largely vindicate formulae for ejecta properties but do not find the expected dependence of disk mass on neutron star compaction. We track the early evolution of the accretion disk, largely driven by shocking and fallback inflow, and do find notable equation of state effects on the structure of this early-time, neutrino-bright disk.

We present results from simulations of core-collapse supernovae in FLASH using a newly-implemented multidimensional neutrino transport scheme and a newly-implemented general relativistic (GR) treatment of gravity. We use a two-moment method with an analytic closure (so-called M1 transport) for the neutrino transport. This transport is multienergy, multispecies, velocity-dependent and truly multidimensional, i.e., we do not assume the commonly used "ray-by-ray" approximation. Our GR gravity is implemented in our Newtonian hydrodynamics simulations via an effective relativistic potential (GREP) that closely reproduces the GR structure of neutron stars and has been shown to match GR simulations of core collapse quite well. In axisymmetry, we simulate core-collapse supernovae with four different progenitor models in both Newtonian and GR gravity. We find that the more compact protoneutron star structure realized in simulations with GR gravity gives higher neutrino luminosities and higher neutrino energies. These differences in turn give higher neutrino heating rates (upwards of $\sim$20-30% over the corresponding Newtonian gravity simulations) that increase the efficacy of the neutrino mechanism. Three of the four models successfully explode in the simulations assuming GREP gravity. In our Newtonian gravity simulations, two of the four models explode, but at times much later than observed in our GR gravity simulations. Our results, both in Newtonian and GR gravity, compare well with several other studies in the literature. These results conclusively show that the approximation of Newtonian gravity for simulating the core-collapse supernova central engine is not acceptable. We also simulate four additional models in GR gravity to highlight the growing disparity between parameterized 1D models of core-collapse supernovae and the current generation of 2D models.

Core-collapse supernova explosions are driven by a central engine that converts a small fraction of the gravitational binding energy released during core collapse to outgoing kinetic energy. The suspected mode for this energy conversion is the neutrino mechanism, where a fraction of the neutrinos emitted from the newly formed protoneutron star are absorbed by and heat the matter behind the supernova shock. Accurate neutrino-matter interaction terms are crucial for simulating these explosions. In this proceedings for IAUS 331, SN 1987A, 30 years later, we explore several corrections to the neutrino-nucleon scattering opacity and demonstrate the effect on the dynamics of the core-collapse supernova central engine via two dimensional neutrino-radiation-hydrodynamics simulations. Our results reveal that the explosion properties are sensitive to corrections to the neutral-current scattering cross section at the 10-20% level, but only for densities at or above $\sim 10^{12}$ g cm$^{-3}$

We report on a set of long-term general-relativistic three-dimensional (3D) multi-group (energy-dependent) neutrino-radiation hydrodynamics simulations of core-collapse supernovae. We employ a full 3D two-moment scheme with the local M1 closure, three neutrino species, and 12 energy groups per species. With this, we follow the post-core-bounce evolution of the core of a nonrotating $27$-$M_\odot$ progenitor in full unconstrained 3D and in octant symmetry for $\gtrsim$$ 380\,\mathrm{ms}$. We find the development of an asymmetric runaway explosion in our unconstrained simulation. We test the resolution dependence of our results and, in agreement with previous work, find that low resolution artificially aids explosion and leads to an earlier runaway expansion of the shock. At low resolution, the octant and full 3D dynamics are qualitatively very similar, but at high resolution, only the full 3D simulation exhibits the onset of explosion.

Gravitational waves (GWs) generated by axisymmetric rotating collapse, bounce, and early postbounce phases of a galactic core-collapse supernova will be detectable by current-generation gravitational wave observatories. Since these GWs are emitted from the quadrupole-deformed nuclear-density core, they may encode information on the uncertain nuclear equation of state (EOS). We examine the effects of the nuclear EOS on GWs from rotating core collapse and carry out 1824 axisymmetric general-relativistic hydrodynamic simulations that cover a parameter space of 98 different rotation profiles and 18 different EOS. We show that the bounce GW signal is largely independent of the EOS and sensitive primarily to the ratio of rotational to gravitational energy, and at high rotation rates, to the degree of differential rotation. The GW frequency of postbounce core oscillations shows stronger EOS dependence that can be parameterized by the core's EOS-dependent dynamical frequency $\sqrt{G\bar{\rho}_c}$. We find that the ratio of the peak frequency to the dynamical frequency follows a universal trend that is obeyed by all EOS and rotation profiles and that indicates that the nature of the core oscillations changes when the rotation rate exceeds the dynamical frequency. We find that differences in the treatments of low-density nonuniform nuclear matter, of the transition from nonuniform to uniform nuclear matter, and in the description of nuclear matter up to around twice saturation density can mildly affect the GW signal. We find that approximations and uncertainties in electron capture rates can lead to variations in the GW signal that are of comparable magnitude to those due to different nuclear EOS. This emphasizes the need for reliable nuclear electron capture rates and for self-consistent multi-dimensional neutrino radiation-hydrodynamic simulations of rotating core collapse.

Binary neutron star mergers are promising sources of gravitational waves for ground-based detectors such as Advanced LIGO. Neutron-rich material ejected by these mergers may also be the main source of r-process elements in the Universe, while radioactive decays in the ejecta can power bright electromagnetic post-merger signals. Neutrino-matter interactions play a critical role in the evolution of the composition of the ejected material, which significantly impacts the outcome of nucleosynthesis and the properties of the associated electromagnetic signal. In this work, we present a simulation of a binary neutron star merger using an improved method for estimating the average neutrino energies in our energy-integrated neutrino transport scheme. These energy estimates are obtained by evolving the neutrino number density in addition to the neutrino energy and flux densities. We show that significant changes are observed in the composition of the polar ejecta when comparing our new results with earlier simulations in which the neutrino spectrum was assumed to be the same everywhere in optically thin regions. In particular, we find that material ejected in the polar regions is less neutron rich than previously estimated. Our new estimates of the composition of the polar ejecta make it more likely that the color and timescale of the electromagnetic signal depend on the orientation of the binary with respect to an observer's line-of-sight. These results also indicate that important observable properties of neutron star mergers are sensitive to the neutrino energy spectrum, and may need to be studied through simulations including a more accurate, energy-dependent neutrino transport scheme.

Nuclear pasta, with nucleons arranged into tubes, sheets, or other complex shapes, is expected in core collapse supernovae (SNe) at just below nuclear density. We calculate the additional opacity from neutrino-pasta coherent scattering using molecular dynamics simulations. We approximately include this opacity in simulations of SNe. We find that pasta slows neutrino diffusion and greatly increases the neutrino signal at late times of 10 or more seconds after stellar core collapse. This signal, for a galactic SN, should be clearly visible in large detectors such as Super-Kamiokande.

We generalize 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 S_A is significantly reduced from one even at low densities. We provide a simple fit S_A^f(n,T,Y_p) of the axial response as a function of density n, temperature T and proton fraction Y_p. This fit reproduces our model independent virial results at low densities and reproduces the Burrows and Sawyer random phase approximation (RPA) results at high densities. Our fit can be incorporated into supernova simulations in a straight forward manner. Preliminary one dimensional supernova simulations suggest that the reduction in the axial response may enhance neutrino heating rates in the gain region during the accretion phase of a core-collapse supernovae.

The DUNE/LBNF program aims to address key questions in neutrino physics and astroparticle physics. Realizing DUNE's potential to reconstruct low-energy particles in the 10-100 MeV energy range will bring significant benefits for all DUNE's science goals. In neutrino physics, low-energy sensitivity will improve neutrino energy reconstruction in the GeV range relevant for the kinematics of DUNE's long-baseline oscillation program. In astroparticle physics, low-energy capabilities will make DUNE's far detectors the world's best apparatus for studying the electron-neutrino flux from a supernova. This will open a new window to unrivaled studies of the dynamics and neutronization of a star's central core in real time, the potential discovery of the neutrino mass hierarchy, provide new sensitivity to physics beyond the Standard Model, and evidence of neutrino quantum-coherence effects. The same capabilities will also provide new sensitivity to `boosted dark matter' models that are not observable in traditional direct dark matter detectors.

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 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_{\odot} neutron star, 5.6 M_{\odot} black hole), high spin (black hole J/M^2=0.9) system with the K_0=220 MeV Lattimer-Swesty equation of state. We find that about 0.08 M_{\odot} of nuclear matter is ejected from the system, while another 0.3 M_{\odot} forms a hot, compact accretion disk. The primary effects of the escaping neutrinos are (i) to make the disk much denser and more compact, (ii) to cause the average electron fraction Y_e of the disk to rise to about 0.2 and then gradually decrease again, and (iii) to gradually cool the disk. The disk is initially hot (T~6 MeV) and luminous in neutrinos (L_{\nu}~10^{54} erg s^{-1}), but the neutrino luminosity decreases by an order of magnitude over 50 ms of post-merger evolution.

We study the merger of binary neutron stars with different mass ratios adopting three different realistic, microphysical nuclear equations of state, as well as incorporating neutrino cooling effects. In particular, we concentrate on the influence of the equation of state on the gravitational wave signature and also on its role, in combination with neutrino cooling, in determining the properties of the resulting hypermassive neutron star, of the neutrinos produced, and of the ejected material. The ejecta we find are consistent with other recent studies that find that small mass ratios produce more ejecta than equal mass cases (up to some limit) and this ejecta is more neutron rich. This trend indicates the importance with future kilonovae observations of measuring the individual masses of an associated binary neutron star system, presumably from concurrent gravitational wave observations, in order to be able to extract information about the nuclear equation of state

Neutron star mergers are among the most promising sources of gravitational waves for advanced ground-based detectors. These mergers are also expected to power bright electromagnetic signals, in the form of short gamma-ray bursts, infrared/optical transients, and radio emission. Simulations of these mergers with fully general relativistic codes are critical to understand the merger and post-merger gravitational wave signals and their neutrinos and electromagnetic counterparts. In this paper, we employ the SpEC code to simulate the merger of low-mass neutron star binaries (two $1.2M_\odot$ neutron stars) for a set of three nuclear-theory based, finite temperature equations of state. We show that the frequency peaks of the post-merger gravitational wave signal are in good agreement with predictions obtained from simulations using a simpler treatment of gravity. We find, however, that only the fundamental mode of the remnant is excited for long periods of time: emission at the secondary peaks is damped on a millisecond timescale in the simulated binaries. For such low-mass systems, the remnant is a massive neutron star which, depending on the equation of state, is either permanently stable or long-lived. We observe strong excitations of l=2, m=2 modes, both in the massive neutron star and in the form of hot, shocked tidal arms in the surrounding accretion torus. We estimate the neutrino emission of the remnant using a neutrino leakage scheme and, in one case, compare these results with a gray two-moment neutrino transport scheme. We confirm the complex geometry of the neutrino emission, also observed in previous simulations with neutrino leakage, and show explicitly the presence of important differences in the neutrino luminosity, disk composition, and outflow properties between the neutrino leakage and transport schemes.

We present Sedonu, a new open source, steady-state, special relativistic Monte Carlo (MC) neutrino transport code, available at bitbucket.org/srichers/sedonu. The code calculates the energy- and angle-dependent neutrino distribution function on fluid backgrounds of any number of spatial dimensions, calculates the rates of change of fluid internal energy and electron fraction, and solves for the equilibrium fluid temperature and electron fraction. We apply this method to snapshots from two-dimensional simulations of accretion disks left behind by binary neutron star mergers, varying the input physics and comparing to the results obtained with a leakage scheme for the case of a central black hole and a central hypermassive neutron star. Neutrinos are guided away from the densest regions of the disk and escape preferentially around 45 degrees from the equatorial plane. Neutrino heating is strengthened by MC transport a few scale heights above the disk midplane near the innermost stable circular orbit, potentially leading to a stronger neutrino-driven wind. Neutrino cooling in the dense midplane of the disk is stronger when using MC transport, leading to a globally higher cooling rate by a factor of a few and a larger leptonization rate by an order of magnitude. We calculate neutrino pair annihilation rates and estimate that an energy of 2.8e46 erg is deposited within 45 degrees of the symmetry axis over 300 ms when a central BH is present. Similarly, 1.9e48 erg is deposited over 3 s when an HMNS sits at the center, but neither estimate is likely to be sufficient to drive a GRB jet.

We present an open-source update to the spherically-symmetric, general-relativistic hydrodynamics, core-collapse supernova (CCSN) code GR1D (O'Connor & Ott, 2010, CQG, 27, 114103). The source code is available at http://www.GR1Dcode.org. We extend its capabilities to include a general relativistic treatment of neutrino transport based on the moment formalisms of Shibata et al., 2011, PTP, 125, 1255 and Cardall et al., 2013, PRD, 87 103004. We pay special attention to implementing and testing numerical methods and approximations that lessen the computational demand of the transport scheme by removing the need to invert large matrices. This is especially important for the implementation and development of moment-like transport methods in two and three dimensions. A critical component of neutrino transport calculations are the neutrino-matter interaction coefficients that describe the production, absorption, scattering, and annihilation of neutrinos. In this article we also describe our open-source, neutrino interaction library NuLib (available at http://www.nulib.org). We believe that an open-source approach to describing these interactions is one of the major steps needed to progress towards robust models of CCSNe and robust predictions of the neutrino signal. We show, via comparisons to full Boltzmann neutrino transport simulations of CCSNe, that our neutrino transport code performs remarkably well. Furthermore, we show that the methods and approximations we employ to increase efficiency do not decrease the fidelity of our results. We also test the ability of our general relativistic transport code to model failed CCSN by evolving a 40 solar-mass progenitor to the onset of collapse to a black hole.

A weak-rate library aimed at investigating the sensitivity of astrophysical environments to variations of electron-capture rates on medium-heavy nuclei has been developed. With this library, the sensitivity of the core-collapse and early post-bounce phases of core-collapse supernovae to nuclear electron-capture is examined by systematically and statistically varying electron-capture rates of individual nuclei. The rates are adjusted by factors consistent with uncertainties indicated by comparing theoretical rates to those deduced from charge-exchange and $\beta$-decay measurements. To ensure a model independent assessment, sensitivity studies across a comprehensive set of progenitors and equations of state are performed. In our systematic study, we find a +16/-4 % range in the mass of the inner-core at the time of shock formation and a $\pm$20% range of peak {\nu}e-luminosity during the deleptonization burst. These ranges are each five times as large as those seen from a separate progenitor study in which we evaluated the sensitivity of these parameters to 32 presupernova stellar models. It is also found that the simulations are more sensitive to a reduction in the electron-capture rates than an enhancement, and in particular to the reduction in the rates for neutron-rich nuclei near the N = 50 closed neutron-shell. As measurements for medium-heavy (A > 65) and neutron-rich nuclei are sparse, and because accurate theoretical models which account for nuclear structure considerations on the individual nucleus level are not readily available, rates for these nuclei may be overestimated. If more accurate estimates confirm this, results from this study indicate that significant changes to the core-collapse trajectory can be expected. For this reason, experimental and theoretical efforts should focus in this region of the nuclear chart.

We study the merger of binary neutron stars using different realistic, microphysical nuclear equations of state, as well as incorporating magnetic field and neutrino cooling effects. In particular, we concentrate on the influence of the equation of state on the gravitational wave signature and also on its role, in combination with cooling and electromagnetic effects, in determining the properties of the hypermassive neutron star resulting from the merger, the production of neutrinos, and the characteristics of ejecta from the system. The ejecta we find are consistent with other recent studies that find soft equations of state produce more ejecta than stiffer equations of state. Moreover, the degree of neutron richness increases for softer equations of state. In light of reported kilonova observations (associated to GRB~130603B and GRB~060614) and the discovery of relatively low abundances of heavy, radioactive elements in deep sea deposits (with respect to possible production via supernovae), we speculate that a soft EoS might be preferred---because of its significant production of sufficiently neutron rich ejecta---if such events are driven by binary neutron star mergers. We also find that realistic magnetic field strengths, obtained with a sub-grid model tuned to capture magnetic amplification via the Kelvin-Helmholtz instability at merger, are generally too weak to affect the gravitational wave signature post-merger within a time scale of $\approx 10$~ms but can have subtle effects on the post-merger dynamics.

We present a first simulation of the post-merger evolution of a black hole-neutron star binary in full general relativity using an energy-integrated general relativistic truncated moment formalism for neutrino transport. We describe our implementation of the moment formalism and important tests of our code, before studying the formation phase of a disk after a black hole-neutron star merger. We use as initial data an existing general relativistic simulation of the merger of a neutron star of 1.4 solar mass with a black hole of 7 solar mass and dimensionless spin a/M=0.8. Comparing with a simpler leakage scheme for the treatment of the neutrinos, we find noticeable differences in the neutron to proton ratio in and around the disk, and in the neutrino luminosity. We find that the electron neutrino luminosity is much lower in the transport simulations, and that the remnant is less neutron-rich. The spatial distribution of the neutrinos is significantly affected by relativistic effects. Over the short timescale evolved, we do not observe purely neutrino-driven outflows. However, a small amount of material (3e-4Msun) is ejected in the polar region during the circularization of the disk. Most of that material is ejected early in the formation of the disk, and is fairly neutron rich. Through r-process nucleosynthesis, that material should produce high-opacity lanthanides in the polar region, and could thus affect the lightcurve of radioactively powered electromagnetic transients. We also show that by the end of the simulation, while the bulk of the disk is neutron-rich, its outer layers have a higher electron fraction. As that material would be the first to be unbound by disk outflows on longer timescales, the changes in Ye experienced during the formation of the disk could have an impact on the nucleosynthesis outputs from neutrino-driven and viscously-driven outflows. [Abridged]

We present a first exploration of the results of neutron star-black hole mergers using black hole masses in the most likely range of $7M_\odot-10M_\odot$, a neutrino leakage scheme, and a modeling of the neutron star material through a finite-temperature nuclear-theory based equation of state. In the range of black hole spins in which the neutron star is tidally disrupted ($\chi_{\rm BH}\gtrsim 0.7$), we show that the merger consistently produces large amounts of cool ($T\lesssim 1\,{\rm MeV}$), unbound, neutron-rich material ($M_{\rm ej}\sim 0.05M_\odot-0.20M_\odot$). A comparable amount of bound matter is initially divided between a hot disk ($T_{\rm max}\sim 15\,{\rm MeV}$) with typical neutrino luminosity $L_\nu\sim 10^{53}\,{\rm erg/s}$, and a cooler tidal tail. After a short period of rapid protonization of the disk lasting $\sim 10\,{\rm ms}$, the accretion disk cools down under the combined effects of the fall-back of cool material from the tail, continued accretion of the hottest material onto the black hole, and neutrino emission. As the temperature decreases, the disk progressively becomes more neutron-rich, with dimmer neutrino emission. This cooling process should stop once the viscous heating in the disk (not included in our simulations) balances the cooling. These mergers of neutron star-black hole binaries with black hole masses $M_{\rm BH}\sim 7M_\odot-10M_\odot$ and black hole spins high enough for the neutron star to disrupt provide promising candidates for the production of short gamma-ray bursts, of bright infrared post-merger signals due to the radioactive decay of unbound material, and of large amounts of r-process nuclei.

The merger of two neutron stars leaves behind a rapidly spinning hypermassive object whose survival is believed to depend on the maximum mass supported by the nuclear equation of state, angular momentum redistribution by (magneto-)rotational instabilities, and spindown by gravitational waves. The high temperatures (~5-40 MeV) prevailing in the merger remnant may provide thermal pressure support that could increase its maximum mass and, thus, its life on a neutrino-cooling timescale. We investigate the role of thermal pressure support in hypermassive merger remnants by computing sequences of spherically-symmetric and axisymmetric uniformly and differentially rotating equilibrium solutions to the general-relativistic stellar structure equations. Using a set of finite-temperature nuclear equations of state, we find that hot maximum-mass critically spinning configurations generally do not support larger baryonic masses than their cold counterparts. However, subcritically spinning configurations with mean density of less than a few times nuclear saturation density yield a significantly thermally enhanced mass. Even without decreasing the maximum mass, cooling and other forms of energy loss can drive the remnant to an unstable state. We infer secular instability by identifying approximate energy turning points in equilibrium sequences of constant baryonic mass parametrized by maximum density. Energy loss carries the remnant along the direction of decreasing gravitational mass and higher density until instability triggers collapse. Since configurations with more thermal pressure support are less compact and thus begin their evolution at a lower maximum density, they remain stable for longer periods after merger.

We incorporate realistic, tabulated equations of state into fully relativistic simulations of magnetized neutron stars along with a neutrino leakage scheme which accounts for cooling via neutrino emission. Both these improvements utilize open-source code (GR1D) and tables from http://stellarcollapse.org. Our implementation makes use of a novel method for the calculation of the optical depth which simplifies its use with distributed adaptive mesh refinement. We present various tests with and without magnetization and preliminary results both from single stars and from the merger of a binary system.

Three-dimensional simulations of core-collapse supernovae are granting new insight into the as-yet uncertain mechanism that drives successful explosions. While there is still debate about whether explosions are obtained more easily in 3D than in 2D, it is undeniable that there exist qualitative and quantitative differences between the results of 3D and 2D simulations. We present an extensive set of high-resolution one-, two-, and three-dimensional core-collapse supernova simulations with multispecies neutrino leakage carried out in two different progenitors. Our simulations confirm the results of Couch (2013) indicating that 2D explodes more readily than 3D. We argue that this is due to the inadequacies of 2D to accurately capture important aspects of the three-dimensional dynamics. We find that without artificially enhancing the neutrino heating rate we do not obtain explosions in 3D. We examine the development of neutrino-driven convection and the standing accretion shock instability and find that, in separate regimes, either instability can dominate. We find evidence for growth of the standing accretion shock instability for both 15-$M_\odot$ and 27-$M_\odot$ progenitors, however, it is weaker in 3D exploding models. The growth rate of both instabilities is artificially enhanced along the symmetry axis in 2D as compared with our axis-free 3D Cartesian simulations. Our work highlights the growing consensus that core-collapse supernovae must be studied in 3D if we hope to solve the mystery of how the explosions are powered.

We study the three-dimensional (3D) hydrodynamics of the post-core-bounce phase of the collapse of a 27-solar-mass star and pay special attention to the development of the standing accretion shock instability (SASI) and neutrino-driven convection. To this end, we perform 3D general-relativistic simulations with a 3-species neutrino leakage scheme. The leakage scheme captures the essential aspects of neutrino cooling, heating, and lepton number exchange as predicted by radiation-hydrodynamics simulations. The 27-solar-mass progenitor was studied in 2D by B. Mueller et al. (ApJ 761:72, 2012), who observed strong growth of the SASI while neutrino-driven convection was suppressed. In our 3D simulations, neutrino-driven convection grows from numerical perturbations imposed by our Cartesian grid. It becomes the dominant instability and leads to large-scale non-oscillatory deformations of the shock front. These will result in strongly aspherical explosions without the need for large-scale SASI shock oscillations. Low-l-mode SASI oscillations are present in our models, but saturate at small amplitudes that decrease with increasing neutrino heating and vigor of convection. Our results, in agreement with simpler 3D Newtonian simulations, suggest that once neutrino-driven convection is started, it is likely to become the dominant instability in 3D. Whether it is the primary instability after bounce will ultimately depend on the physical seed perturbations present in the cores of massive stars. The gravitational wave signal, which we extract and analyze for the first time from 3D general-relativistic models, will serve as an observational probe of the postbounce dynamics and, in combination with neutrinos, may allow us to determine the primary hydrodynamic instability.

We present results from a new set of 3D general-relativistic hydrodynamic simulations of rotating iron core collapse. We assume octant symmetry and focus on axisymmetric collapse, bounce, the early postbounce evolution, and the associated gravitational wave (GW) and neutrino signals. We employ a finite-temperature nuclear equation of state, parameterized electron capture in the collapse phase, and a multi-species neutrino leakage scheme after bounce. The latter captures the important effects of deleptonization, neutrino cooling and heating and enables approximate predictions for the neutrino luminosities in the early evolution after core bounce. We consider 12-solar-mass and 40-solar-mass presupernova models and systematically study the effects of (i) rotation, (ii) progenitor structure, and (iii) postbounce neutrino leakage on dynamics, GW, and, neutrino signals. We demonstrate, that the GW signal of rapidly rotating core collapse is practically independent of progenitor mass and precollapse structure. Moreover, we show that the effects of neutrino leakage on the GW signal are strong only in nonrotating or slowly rotating models in which GW emission is not dominated by inner core dynamics. In rapidly rotating cores, core bounce of the centrifugally-deformed inner core excites the fundamental quadrupole pulsation mode of the nascent protoneutron star. The ensuing global oscillations (f~700-800 Hz) lead to pronounced oscillations in the GW signal and correlated strong variations in the rising luminosities of antineutrino and heavy-lepton neutrinos. We find these features in cores that collapse to protoneutron stars with spin periods <~ 2.5 ms and rotational energies sufficient to drive hyper-energetic core-collapse supernova explosions. Hence, joint GW + neutrino observations of a core collapse event could deliver strong evidence for or against rapid core rotation. [abridged]

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 1D+rotation stellar-collapse simulations on the progenitor models of Woosley & Heger (2006) 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 compromises 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.

We perform spherically-symmetric general-relativistic simulations of core collapse and the postbounce preexplosion phase in 32 presupernova stellar models of solar metallicity with zero-age-main-sequence masses of 12 M_{sun} to 120 M_{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 (~M/R) of the progenitor star's inner regions via the accretion rate in the preexplosion phase. We find that these qualitative trends depend only weakly on the nuclear equation of state, but quantitative observational statements will require independent constraints on the equation of state 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.

Core-collapse supernovae are among the most energetic cosmic cataclysms. They are prodigious emitters of neutrinos and quite likely strong galactic sources of gravitational waves. Observation of both neutrinos and gravitational waves from the next galactic or near extragalactic core-collapse supernova will yield a wealth of information on the explosion mechanism, but also on the structure and angular momentum of the progenitor star, and on aspects of fundamental physics such as the equation of state of nuclear matter at high densities and low entropies. In this contribution to the proceedings of the Neutrino 2012 conference, we summarize recent progress made in the theoretical understanding and modeling of core-collapse supernovae. In this, our emphasis is on multi-dimensional processes involved in the explosion mechanism such as neutrino-driven convection and the standing accretion shock instability. As an example of how supernova neutrinos can be used to probe fundamental physics, we discuss how the rise time of the electron antineutrino flux observed in detectors can be used to probe the neutrino mass hierarchy. Finally, we lay out aspects of the neutrino and gravitational-wave signature of core-collapse supernovae and discuss the power of combined analysis of neutrino and gravitational wave data from the next galactic core-collapse supernova.

Core-collapse supernovae may depend sensitively on charged current neutrino interactions in warm, low density neutron rich matter. A proton in neutron rich matter is more tightly bound than is a neutron. This energy shift \Delta U increases the electron energy in \nu_e + n --> p + e, increasing the available phase space and absorption cross section. Likewise \Delta U decreases the positron energy in \bar \nu_e + p --> n + e^+, decreasing the phase space and cross section. We have calculated \Delta U using a model independent virial expansion and we find \Delta U is much larger, at low densities, than the predictions of many mean field models. Therefore \Delta U could have a significant impact on charged current neutrino interactions in supernovae. Preliminary simulations of the accretion phase of core-collapse supernovae find that \Delta U increases \bar \nu_e energies and decreases the \nu_e luminosity.

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

We explore the effects of collective neutrino flavor oscillations due to neutrino-neutrino interactions on the neutrino heating behind a stalled core-collapse supernova shock. We carry out axisymmetric (2D) radiation-hydrodynamic core-collapse supernova simulations, tracking the first 400 ms of the post-core-bounce evolution in 11.2 solar mass and 15 solar mass progenitor stars. Using inputs from these 2D simulations, we perform neutrino flavor oscillation calculations in multi-energy single-angle and multi-angle single-energy approximations. Our results show that flavor conversions do not set in until close to or outside the stalled shock, enhancing heating by not more than a few percent in the most optimistic case. Consequently, we conclude that the postbounce pre-explosion dynamics of standard core-collapse supernovae remains unaffected by neutrino oscillations. Multi-angle effects in regions of high electron density can further inhibit collective oscillations, strengthening our conclusion.

Core-collapse supernovae are among Nature's grandest explosions. They are powered by the energy released in gravitational collapse and include a rich set of physical phenomena involving all fundamental forces and many branches of physics and astrophysics. We summarize the current state of core-collapse supernova theory and discuss the current set of candidate explosion mechanisms under scrutiny as core-collapse supernova modeling is moving towards self-consistent three-dimensional simulations. Recent work in nuclear theory and neutron star mass and radius measurements are providing new constraints for the nuclear equation of state. We discuss these new developments and their impact on core-collapse supernova modeling. Neutrino-neutrino forward scattering in the central regions of core-collapse supernovae can lead to collective neutrino flavor oscillations that result in swaps of electron and heavy-lepton neutrino spectra. We review the rapid progress that is being made in understanding these collective oscillations and their potential impact on the core-collapse supernova explosion mechanism.

We generate a second equation of state (EOS) of nuclear matter for a wide range of temperatures, densities, and proton fractions for use in supernovae, neutron star mergers, and black hole formation simulations. We employ full relativistic mean field (RMF) calculations for matter at intermediate density and high density, and the Virial expansion of a non-ideal gas for matter at low density. For this EOS we use the RMF effective interaction FSUGold, whereas our earlier EOS was based on the RMF effective interaction NL3. The FSUGold interaction has a lower pressure at high densities compared to the NL3 interaction. We calculate the resulting EOS at over 100,000 grid points in the temperature range $T$ = 0 to 80 MeV, the density range $n_B$ = 10$^{-8}$ to 1.6 fm$^{-3}$, and the proton fraction range $Y_p$ = 0 to 0.56. We then interpolate these data points using a suitable scheme to generate a thermodynamically consistent equation of state table on a finer grid. We discuss differences between this EOS, our NL3 based EOS, and previous EOSs by Lattimer-Swesty and H. Shen et al for the thermodynamic properties, composition, and neutron star structure. The original FSUGold interaction produces an EOS, that we call FSU1.7, that has a maximum neutron star mass of 1.7 solar masses. A modification in the high density EOS is introduced to increase the maximum neutron star mass to 2.1 solar masses and results in a slightly different EOS that we call FSU2.1. The EOS tables for FSU1.7 and FSU2.1 are available for download.

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^2 = 1 by the appearance of nonaxisymmetric rotational instabilities.

We investigate several aspects of black hole formation in failing core-collapse supernovae using 1D general-relativistic hydrodynamic simulations. We use the open-source code GR1D and incorporate into it nucleon-nucleon Bremsstrahlung, a crucial neutrino pair-production channel. We focus on how various thermal effects can influence the postbounce supernova evolution towards black hole formation. By performing simulations with and without nucleon-nucleon Bremsstrahlung, we investigate the sensitivity of black hole formation to thermal support in the protoneutron star. We also investigate delayed black hole formation by artificially driving explosions in an extreme model where the protoneutron star is initially thermally supported above the maximum baryonic cold neutron star mass but then collapses to a black hole after neutrino cooling removes sufficient thermal support.

We perform 3+1 general relativistic simulations of rotating core collapse in the context of the collapsar model for long gamma-ray bursts. We employ a realistic progenitor, rotation based on results of stellar evolution calculations, and a simplified equation of state. Our simulations track self-consistently collapse, bounce, the postbounce phase, black hole formation, and the subsequent early hyperaccretion phase. We extract gravitational waves from the spacetime curvature and identify a unique gravitational wave signature associated with the early phase of collapsar formation.

We discuss results from simulations of black hole formation in failing core-collapse supernovae performed with the code GR1D, a new open-source Eulerian spherically-symmetric general-relativistic hydrodynamics code. GR1D includes rotation in an approximate way (1.5D), comes with multiple finite-temperature nuclear equations of state (EOS), and treats neutrinos in the post-core-bounce phase via a 3-flavor leakage scheme and a heating prescription. We chose the favored K_0=220 MeV-variant of the Lattimer & Swesty (1990) EOS and present collapse calculations using the progenitor models of Limongi & Chieffi (2006). We show that there is no direct (or ``prompt'') black hole formation in the collapse of ordinary massive stars (8 M_Sun ~< M_ZAMS ~< 100 M_Sun) and present first results from black hole formation simulations that include rotation.

We present the new open-source spherically-symmetric general-relativistic (GR) hydrodynamics code GR1D. It is based on the Eulerian formulation of GR hydrodynamics (GRHD) put forth by Romero-Ibanez-Gourgoulhon and employs radial-gauge, polar-slicing coordinates in which the 3+1 equations simplify substantially. We discretize the GRHD equations with a finite-volume scheme, employing piecewise-parabolic reconstruction and an approximate Riemann solver. GR1D is intended for the simulation of stellar collapse to neutron stars and black holes and will also serve as a testbed for modeling technology to be incorporated in multi-D GR codes. Its GRHD part is coupled to various finite-temperature microphysical equations of state in tabulated form that we make available with GR1D. An approximate deleptonization scheme for the collapse phase and a neutrino-leakage/heating scheme for the postbounce epoch are included and described. We also derive the equations for effective rotation in 1D and implement them in GR1D. We present an array of standard test calculations and also show how simple analytic equations of state in combination with presupernova models from stellar evolutionary calculations can be used to study qualitative aspects of black hole formation in failing rotating core-collapse supernovae. In addition, we present a simulation with microphysical EOS and neutrino leakage/heating of a failing core-collapse supernova and black hole formation in a presupernova model of a 40 solar mass zero-age main-sequence star. We find good agreement on the time of black hole formation (within 20%) and last stable protoneutron star mass (within 10%) with predictions from simulations with full Boltzmann neutrino radiation hydrodynamics.

Core-collapse supernovae are among Nature's most energetic events. They mark the end of massive star evolution and pollute the interstellar medium with the life-enabling ashes of thermonuclear burning. Despite their importance for the evolution of galaxies and life in the universe, the details of the core-collapse supernova explosion mechanism remain in the dark and pose a daunting computational challenge. We outline the multi-dimensional, multi-scale, and multi-physics nature of the core-collapse supernova problem and discuss computational strategies and requirements for its solution. Specifically, we highlight the axisymmetric (2D) radiation-MHD code VULCAN/2D and present results obtained from the first full-2D angle-dependent neutrino radiation-hydrodynamics simulations of the post-core-bounce supernova evolution. We then go on to discuss the new code Zelmani which is based on the open-source HPC Cactus framework and provides a scalable AMR approach for 3D fully general-relativistic modeling of stellar collapse, core-collapse supernovae and black hole formation on current and future massively-parallel HPC systems. We show Zelmani's scaling properties to more than 16,000 compute cores and discuss first 3D general-relativistic core-collapse results.

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.

We extend the virial equation of state to include 3H and 3He nuclei, and predict significant mass-three fractions near the neutrinosphere in supernovae. While alpha particles are often more abundant, we demonstrate that energy transfer cross-sections for muon and tau neutrinos at low densities are dominated by breakup of the loosely-bound 3H and 3He nuclei. The virial coefficients involving A=3 nuclei are calculated directly from the corresponding nucleon-3H and nucleon-3He scattering phase shifts. For the neutral-current inelastic cross-sections and the energy transfer cross sections, we perform ab-initio calculations based on microscopic two- and three-nucleon interactions and meson-exchange currents.

We have developed radiative transfer models of the radio emission from colliding-wind binaries (CWB) based on a hydrodynamical treatment of the wind-collision region (WCR). The archetype of CWB systems is the 7.9-yr period binary WR140, which exhibits dramatic variations at radio wavelengths. High-resolution radio observations of WR140 permit a determination of several system parameters, particularly orbit inclination and distance, that are essential for any models of this system. A model fit to data at orbital phase 0.9 is shown, and some short comings of our model described.

Radio emission models of colliding wind binaries (CWBs) have been discussed by Dougherty et al. (2003). We extend these models by considering the temporal and spatial evolution of the energy distribution of relativistic electrons as they advect downstream from their shock acceleration site. The energy spectrum evolves significantly due to the strength of inverse-Compton (IC) cooling in these systems, and a full numerical evaluation of the synchrotron emission and absorption coefficients is made. We have demonstrated that the geometry of the WCR and the streamlines of the flow within it lead to a spatially dependent break frequency in the synchrotron emission. We therefore do not observe a single, sharp break in the synchrotron spectrum integrated over the WCR, but rather a steepening of the synchrotron spectrum towards higher frequencies. We also observe that emission from the wind-collision region (WCR) may appear brightest near the shocks, since the impact of IC cooling on the non-thermal electron distribution is greatest near the contact discontinuity (CD). We show that the thermal emission from a CWB can mimic a thermal plus non-thermal composite spectrum if the thermal emission from the WCR becomes comparable to that from the unshocked winds, and obtain a number of insights relevant to models of closer systems such as WR140. Finally, we apply our new models to the very wide system WR147. The acceleration of non-thermal electrons appears to be very efficient in our models of WR147, and we suggest that the shock structure may be modified by feedback from the accelerated particles. (Abridged).

WR146 is a WC6+O8 colliding-wind binary (CWB) system with thermal emission from the stellar winds of the two stars, and bright non-thermal emission from the wind-collision region (WCR) where the winds collide. We present high resolution radio observations from 1.4 to 43 GHz that give one of the best quality radio spectra of any CWB to date. Observations at 22 GHz now span 8 years, and reveal the proper motion of the system, allowing comparison of multi-epoch data. VLBI observations show the location of the WCR relative to the stellar components, from which the wind momentum ratio can be shown to be 0.06+/-0.15. The radio spectrum and the spatial distribution of emission are modelled, and we determine the contribution of both stellar winds and the WCR to the observed emission. We show that our current models fail to account for the high frequency spectrum of WR146, and also produce too much emission far from the stagnation point of the wind collision.