Many of the currently available equations of state for core-collapse supernova simulations give large neutron star radii and do not provide large enough neutron star masses, both of which are inconsistent with some recent neutron star observations. In addition, one of the critical uncertainties in the nucleon-nucleon interaction, the nuclear symmetry energy, is not fully explored by the currently available equations of state. In this article, we construct two new equations of state which match recent neutron star observations and provide more flexibility in studying the dependence on nuclear matter properties. The equations of state are also provided in tabular form, covering a wide range in density, temperature, and asymmetry, suitable for astrophysical simulations. These new equations of state are implemented into our spherically symmetric core-collapse supernova model, which is based on general relativistic radiation hydrodynamics with three-flavor Boltzmann neutrino transport. The results are compared with commonly used equations of state in supernova simulations of 11.2 and 40 M progenitors. We consider only equations of state which are fitted to nuclear binding energies and other experimental and observational constraints. We find that central densities at bounce are weakly correlated with L and that there is a moderate influence of the symmetry energy on the evolution of the electron fraction. The new models also obey the previously observed correlation between the time to black hole formation and the maximum mass of an s = 4 neutron star.

We revise the bound from the supernova SN1987A on the coupling of ultralight axion-like particles (ALPs) to photons. In a core-collapse supernova, ALPs would be emitted via the Primakoff process, and eventually convert into gamma rays in the magnetic field of the Milky Way. The lack of a gamma-ray signal in the GRS instrument of the SMM satellite in coincidence with the observation of the neutrinos emitted from SN1987A therefore provides a strong bound on their coupling to photons. Due to the large uncertainty associated with the current bound, we revise this argument, based on state-of-the-art physical inputs both for the supernova models and for the Milky-Way magnetic field. Furthermore, we provide major amendments, such as the consistent treatment of nucleon-degeneracy effects and of the reduction of the nuclear masses in the hot and dense nuclear medium of the supernova. With these improvements, we obtain a new upper limit on the photon-ALP coupling:g_aγ≲ 5.3 × 10^-12GeV^-1, form_a≲ 4.4 × 10^-10eV,and we also give its dependence at larger ALP massesm_a. Moreover, we discuss how much the Fermi-LAT satellite experiment could improve this bound, should a close-enough supernova explode in the near future.

Massive stars end their lives in explosions with kinetic energies on the order of 1051 erg. Immediately after the explosion has been launched, a region of low density and high entropy forms behind the ejecta, which is continuously subject to neutrino heating. The neutrinos emitted from the remnant at the center, the protoneutron star (PNS), heat the material above the PNS surface. This heat is partly converted into kinetic energy, and the material accelerates to an outflow that is known as the neutrino-driven wind. For the first time we simulate the collapse, bounce, explosion, and the neutrino-driven wind phases consistently over more than 20 s. Our numerical model is based on spherically symmetric general relativistic radiation hydrodynamics using spectral three-flavor Boltzmann neutrino transport. In simulations where no explosions are obtained naturally, we model neutrino-driven explosions for low- and intermediate-mass Fe-core progenitor stars by enhancing the charged current reaction rates. In the case of a special progenitor star, the 8.8 O-Ne-Mg-core, the explosion in spherical symmetry was obtained without enhanced opacities. The post-explosion evolution is in qualitative agreement with static steady-state and parametrized dynamic models of the neutrino-driven wind. On the other hand, we generally find lower neutrino luminosities and mean neutrino energies, as well as a different evolutionary behavior of the neutrino luminosities and mean neutrino energies. The neutrino-driven wind is proton-rich for more than 10 s and the contraction of the PNS differs from the assumptions made for the conditions at the inner boundary in previous neutrino-driven wind studies. Despite the moderately high entropies of about 100 k B/baryon and the fast expansion timescales, the conditions found in our models are unlikely to favor r-process nucleosynthesis. The simulations are carried out until the neutrino-driven wind settles down to a quasi-stationary state. About 5 s after the bounce, the peak temperature inside the PNS already starts to decrease because of the continued deleptonization. This moment determines the beginning of a cooling phase dominated by the emission of neutrinos. We discuss the physical conditions of the quasi-static PNS evolution and take the effects of deleptonization and mass accretion from early fallback into account.

We identify an observable imprint of a first-order hadron-quark phase transition at supranuclear densities on the gravitational-wave (GW) emission of neutron-star mergers. Specifically, we show that the dominant postmerger GW frequency $f_{\mathrm{peak}}$ may exhibit a significant deviation from an empirical relation between $f_{\mathrm{peak}}$ and the tidal deformability if a strong first-order phase transition leads to the formation of a gravitationally stable extended quark matter core in the postmerger remnant. A comparison of the GW signatures from a large, representative sample of microphysical, purely hadronic equations of state indicates that this imprint is only observed in those systems which undergo a strong first-order phase transition. Such a shift of the dominant postmerger GW frequency can be revealed by future GW observations, which would provide evidence for the existence of a strong first-order phase transition in the interior of neutron-stars.