Blue-supergiant stars develop into core-collapse supernovae --- one of the most energetic outbursts in the universe --- when all nuclear burning fuel is exhausted in the stellar core. Previous attempts failed to explain observed explosions of such stars which have a zero-age main sequence mass of 50~M$_\odot$ or more. Here we exploit the largely uncertain state of matter at high density, and connect the modeling of such stellar explosions with a first-order phase transition from nuclear matter to the quark-gluon plasma. The resulting energetic supernova explosions can account for a large variety of lightcurves, from peculiar type II to super-luminous events. The remnants are neutron stars with quark matter core, known as hybrid stars, of about 2~M$_\odot$ at birth. A galactic event of this kind could be observable due to the release of a second neutrino burst. Its observation would confirm such a first-order phase transition at densities relevant for astrophysics.

A commonly applied quark matter model in astrophysics is the thermodynamic bag model (tdBAG). The original MIT bag model approximates the effect of quark confinement, but does not explicitly account for the breaking of chiral symmetry, an important property of Quantum Chromodynamics (QCD). It further ignores vector repulsion. The vector-interaction-enhanced bag model (vBag) improves the tdBAG approach by accounting for both dynamical chiral symmetry breaking and repulsive vector interactions. The latter is of particular importance to studies of dense matter in beta-equilibriumto explain the two solar mass maximum mass constraint for neutron stars. The model is motivated by analyses of QCD based Dyson-Schwinger equations (DSE), assuming a simple quark-quark contact interaction. Here, we focus on the study of hybrid neutron star properties resulting from the application of vBag and will discuss possible extensions.

In this review article we discuss selected developments regarding the role of the equation of state (EOS) in simulations of core-collapse supernovae. There are no first-principle calculations of the state of matter under supernova conditions since a wide range of conditions is covered, in terms of density, temperature and isospin asymmetry. Instead, model EOS are commonly employed in supernova studies. These can be divided into regimes with intrinsically different degrees of freedom: heavy nuclei at low temperatures, inhomogeneous nuclear matter where light and heavy nuclei coexist together with unbound nucleons, and the transition to homogeneous matter at high densities and temperatures. In this article we discuss each of these phases with particular view on their role in supernova simulations.

We give an overview about equations of state (EOS) which are currently available for simulations of core-collapse supernovae and neutron star mergers. A few selected important aspects of the EOS, such as the symmetry energy, the maximum mass of neutron stars, and cluster formation, are confronted with constraints from experiments and astrophysical observations. There are just very few models which are compatible even with this very restricted set of constraints. These remaining models illustrate the uncertainty of the uniform nuclear matter EOS at high densities. In addition, at finite temperatures the medium modifications of nuclear clusters represent a conceptual challenge. In conclusion, there has been significant development in the recent years, but there is still need for further improved general purpose EOS tables.

We review various theoretical approaches for the equation of state (EoS) of dense matter, relevant for the description of core-collapse supernovae, compact stars and compact star mergers. The emphasis is put on models that are applicable to all of these scenarios. Such EoS models have to cover large ranges in baryon number density, temperature and isospin asymmetry. The characteristics of matter change dramatically within these ranges, from a mixture of nucleons, nuclei, and electrons to uniform, strongly interacting matter containing nucleons, and possibly other particles such as hyperons or quarks. As the development of an EoS requires joint efforts from many directions we consider different theoretical approaches and discuss relevant experimental and observational constraints which provide insights for future research. Finally, results from applications of the discussed EoS models are summarized.

We compare two classes of hybrid equations of state with a hadron-to-quark matter phase transition in their application to core collapse supernova simulations. The first one uses the quark bag model and describes the transition to three-flavor quark matter at low critical densities. The second one employs a Polyakov-loop extended Nambu-Jona-Lasinio (PNJL) model with parameters describing a phase transition to two-flavor quark matter at higher critical densities. These models possess a distinctly different temperature dependence of their transition densities which turns out to be crucial for the possible appearance of quark matter in supernova cores. During the early post bounce accretion phase quark matter is found only if the phase transition takes place at sufficiently low densities as in the study based on the bag model. The increase critical density with increasing temperature, as obtained for our PNJL parametrization, prevents the formation of quark matter. The further evolution of the core collapse supernova as obtained applying the quark bag model leads to a structural reconfiguration of the central proto-neutron star where, in addition to a massive pure quark matter core, a strong hydrodynamic shock wave forms and a second neutrino burst is released during the shock propagation across the neutrinospheres. We discuss the severe constraints in the freedom of choice of quark matter models and their parametrization due to the recently observed 2 solar mass pulsar and their implications for further studies of core collapse supernovae in the QCD phase diagram.

We investigate nuclear matter at finite temperature and density, including the formation of light clusters up to the alpha particle The novel feature of this work is to include the formation of clusters as well as their dissolution due to medium effects in a systematic way using two many-body theories: a microscopic quantum statistical (QS) approach and a generalized relativistic mean field (RMF) model. Nucleons and clusters are modified by medium effects. Both approaches reproduce the limiting cases of nuclear statistical equilibrium (NSE) at low densities and cluster-free nuclear matter at high densities. The treatment of the cluster dissociation is based on the Mott effect due to Pauli blocking, implemented in slightly different ways in the QS and the generalized RMF approaches. We compare the numerical results of these models for cluster abundances and thermodynamics in the region of medium excitation energies with temperatures T <= 20 MeV and baryon number densities from zero to a few times saturation density. The effect of cluster formation on the liquid-gas phase transition and on the density dependence of the symmetry energy is studied. Comparison is made with other theoretical approaches, in particular those, which are commonly used in astrophysical calculations. The results are relevant for heavy-ion collisions and astrophysical applications.