Some massive stars end their lives as failed core-collapse supernovae (CCSNe) and become black holes (BHs). Although in this class of phenomena the stalled supernova (SN) shock is not revived, the outer stellar envelope can still be partially ejected. This occurs because the hydrodynamic equilibrium of the star is disrupted by the gravitational mass loss of the protoneutron star (PNS) due to neutrino emission. We develop a simple parameterized model that emulates PNS evolution and its neutrino emission and use it to simulate failed CCSNe in spherical symmetry for a wide range of progenitor stars. Our model allows us to study mass ejection of failed CCSNe where the PNS collapses into a BH within ∼100 ms and up to ∼10⁶ s. We perform failed CCSNe simulations for 262 different pre-SN progenitors and determine how the energy and mass of the ejecta depend on progenitor properties and the equation of state (EOS) of dense matter. In the case of a future failed CCSN observation, the trends obtained in our simulations can be used to place constraints on the pre-SN progenitor characteristics, the EOS, and on PNS properties at BH formation time.

We determine the threshold mass for prompt (no bounce) black hole formation in equal-mass neutron star (NS) mergers using a new set of 227 numerical relativity simulations. We consider 23 phenomenological and microphysical finite-temperature equations of state (EOS), including models with hyperons and first-order phase transitions to deconfined quarks. We confirm the existence of EOS-insensitive relations between the threshold mass, binary tidal parameter at the threshold (Λth), maximum mass of nonrotating NSs, and radii of reference mass NSs. We combine the EOS-insensitive relations, phenomenological constraints on NS properties, and observational data from GW170817 to derive an improved lower limit on radii of maximum mass and a 1.6  M⊙ NS of 9.81 and 10.90 km, respectively. We also constrain the radius and quadrupolar tidal deformability (Λ) of a 1.4  M⊙ NS to be larger than 10.74 km and 172, respectively. We consider uncertainties in all independent parameters—fitting coefficients as well as GW170817 masses while reporting the range of radii constraints. We discuss an approach to constrain the upper as well as lower limit of NS maximum mass using future binary NS detections and their identification as prompt or delayed collapse. With future observations, it will be possible to derive even tighter constraints on the properties of matter at and above nuclear density using the method proposed in this work.

Gravitational waves (GWs) provide unobscured insight into the birthplace of neutron stars and black holes in core-collapse supernovae (CCSNe). The nuclear equation of state (EOS) describing these dense environments is yet uncertain, and variations in its prescription affect the proto−neutron star (PNS) and the post-bounce dynamics in CCSN simulations, subsequently impacting the GW emission. We perform axisymmetric simulations of CCSNe with Skyrme-type EOSs to study how the GW signal and PNS convection zone are impacted by two experimentally accessible EOS parameters, (1) the effective mass of nucleons, m⋆, which is crucial in setting the thermal dependence of the EOS, and (2) the isoscalar incompressibility modulus, Ksat. While Ksat shows little impact, the peak frequency of the GWs has a strong effective mass dependence due to faster contraction of the PNS for higher values of m⋆ owing to a decreased thermal pressure. These more compact PNSs also exhibit more neutrino heating, which drives earlier explosions and correlates with the GW amplitude via accretion plumes striking the PNS, exciting the oscillations. We investigate the spatial origin of the GWs and show the agreement between a frequency-radial distribution of the GW emission and a perturbation analysis. We do not rule out overshoot from below via PNS convection as another moderately strong excitation mechanism in our simulations. We also study the combined effect of effective mass and rotation. In all our simulations we find evidence for a power gap near ∼1250 Hz; we investigate its origin and report its EOS dependence.

We study the consequences of a hadron-quark phase transition (PT) in failing core-collapse supernovae (CCSNe) that give birth to stellar-mass black holes (BH). We perform a suite of neutrino-transport general-relativistic hydrodynamic simulations in spherical symmetry with 21 progenitor models and a hybrid equation of state (EoS) including hadrons and quarks. We find that the effect of the PT on the CCSN postbounce dynamics is a function of the bounce compactness parameter xi(2.2). For xi(2.2) greater than or similar to 0.24, the PT leads to a second dynamical collapse of the protocompact star (PCS). While BH formation starts immediately after this second collapse for models with xi(2.2) greater than or similar to 0.51, the PCS experiences a second bounce and oscillations for models with 0.24 less than or similar to x xi(2.2) less than or similar to 0.51. These models emit potent oscillatory neutrino signals with a period of similar to 1 ms for tens of milliseconds after the second bounce, which can be a strong indicator of the PT in failing CCSNe if detected in the future. However, no shock revival occurs and BH formation inevitably takes place in our spherically symmetric simulations. Furthermore, via a diagram of mass-specific entropy evolution of the PCS, the progenitor dependence can be understood through the appearance of a third family of compact stars emerging at large entropy induced by the PT.

Despite their astrophysical relevance, nuclear pasta phases are relatively unstudied at high temperatures. We present molecular dynamics simulations of symmetric nuclear matter with several topologies of lasagna at a range of temperatures to study the pasta-uniform transition. Using the Minkowski functionals, we quantify trends in the occupied volume, surface area, mean breadth, and Euler characteristic. The amplitude of surface displacements of the pasta increase with temperature which produce short-lived topological defects such as holes and filaments near melting, resulting in power laws for increasing surface curvature with temperature. We calculate the static structure factor and report the shear viscosity and thermal conductivity of pasta, finding that the shear viscosity is minimized at the melting temperature. These results may have implications for the thermoelastic properties of nuclear pasta and finite-temperature corrections to the equation of state at pasta densities.

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 pre-supernova 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-specific entropy phase space M-grav(PNS) - (s) over bar is well correlated with the progenitor compactness and is almost EOS independent, apart from the final end point. Furthermore, BH formation occurs, to a very good approximation, soon after the PNS overcomes the maximum gravitational mass supported by a hot NS with constant specific entropy equal to (s) over bar. These results show a path to constraining the temperature dependence of the EOS through the detection of neutrinos from a failed galactic supernova.

Differences in the equation of state (EOS) of dense matter translate into differences in astrophysical simulations and their multimessenger signatures. Thus, extending the number of EOSs for astrophysical simulations allows us to probe the effect of different aspects of the EOS in astrophysical phenomena. In this work, we construct the EOS of hot and dense matter based on the Akmal, Pandharipande, and Ravenhall (APR) model and thereby extend the open-source SROEOS code which computes EOSs of hot dense matter for Skyrme-type parametrizations of the nuclear forces. Unlike Skrme-type models, in which parameters of the interaction are fit to reproduce the energy density of nuclear matter and/or properties of heavy nuclei, the EOS of APR is obtained from potentials resulting from fits to nucleon-nucleon scattering and properties of light nuclei. In addition, this EOS features a phase transition to a spin-isospin ordered state of nucleons, termed a neutral pion condensate, at supranuclear densities. We show that differences in the effective masses between EOSs have consequences for the properties of nuclei in the subnuclear inhomogeneous phase of matter. We also test the new EOS of APR in spherically symmetric core-collapse of massive stars with 15 M-circle dot and 40 M-circle dot, respectively. We find that the phase transition in the EOS of APR speeds up the collapse of the star. However, this phase transition does not generate a second shock wave or another neutrino burst as reported for the hadron-to-quark phase transition. The reason for this difference is that the width of the coexistence region and the latent heat in the EOS of APR are substantially smaller than in the quark-to-hadron transition employed earlier, which results in a significantly smaller softening of the high density EOS.

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-circle dot progenitor star. The core collapse of the 20-M-circle dot progenitor star is simulated in spherical symmetry using the general-relativistic radiation-hydrodynamics code GRID 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.