We investigate the phase transition from hadron to quark matter in the general case without the assumption of chemical equilibrium with respect to weak decays. The effects of net strangeness on charge and isospin fractions, chemical potentials, and temperature are studied in the context of the Chiral Mean Field (CMF) model that incorporates chiral symmetry restoration and deconfinement. The extent to which these quantities are probed during deconfinement in conditions expected to exist in protoneutron stars, binary neutron-star mergers, and heavy-ion collisions is analyzed quantitatively via the construction of 3-dimensional phase diagrams.

Neutron star interiors provide the opportunity to probe properties of cold dense matter in the QCD phase diagram. Utilizing models of dense matter in accord with nuclear systematics at nuclear densities, we investigate the compatibility of deconfined quark cores with current observational constraints on the maximum mass and tidal deformability of neutron stars. We explore various methods of implementing the hadron-to-quark phase transition, specifically, first-order transitions with sharp (Maxwell construction) and soft (Gibbs construction) interfaces, and smooth crossover transitions. We find that within the models we apply, hadronic matter has to be stiff for a first-order phase transition and soft for a crossover transition. In both scenarios and for the equations of state we employed, quarks appear at the center of pre-merger neutron stars in the mass range $\approx 1.0-1.6\,{\rm M}_{\odot}$, with a squared speed of sound $c^2_{\rm QM}\gtrsim 0.4$ characteristic of strong repulsive interactions required to support the recently discovered neutron star masses $\geq 2\,{\rm M}_{\odot}$. We also identify equations of state and phase transition scenarios that are consistent with the bounds placed on tidal deformations of neutron stars in the recent binary merger event GW170817. We emphasize that distinguishing hybrid stars with quark cores from normal hadronic stars is very difficult from the knowledge of masses and radii alone, unless drastic sharp transitions induce distinctive disconnected hybrid branches in the mass-radius relation.

In this work, we discuss the dense matter equation of state (EOS) for the extreme range of conditions encountered in neutron stars and their mergers. The calculation of the properties of such an EOS involves modeling different degrees of freedom (such as nuclei, nucleons, hyperons, and quarks), taking into account different symmetries, and including finite density and temperature effects in a thermodynamically consistent manner. We begin by addressing subnuclear matter consisting of nucleons and a small admixture of light nuclei in the context of the excluded volume approach. We then turn our attention to supranuclear homogeneous matter as described by the Chiral Mean Field (CMF) formalism. Finally, we present results from realistic neutron-star-merger simulations performed using the CMF model that predict signatures for deconfinement to quark matter in gravitational wave signals.

In this work, we study matter in the cores of proto-neutron stars, focusing on the impact of their composition on the stellar structure. We begin by examining the effects of finite temperature (through a fixed entropy per baryon) and lepton fraction on purely nucleonic matter by making use of the DSH model . We then turn our attention to a relativistic mean-field model containing exotic degrees of freedom, the Chiral Mean Field (CMF) model, again, under the conditions of finite temperature and trapped neutrinos. In the latter, since both hyperons and quarks are found in the cores of large-mass stars, their interplay and the possibility of mixtures of phases is taken into account and analyzed. Finally, we discuss how stellar rotation can affect our results.

Differences in the equation of state (EOS) of dense matter translate into differences in astrophysical simulations and their multi-messenger 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 neutral pion condensate at supra-nuclear densities. We show that differences in the effective masses between EOSs have consequences for the properties of nuclei in the sub-nuclear inhomogeneous phase of matter. We also test the new EOS of APR in spherically symmetric core-collapse of massive stars with $15M_\odot$ and $40M_\odot$, 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 onset of the phase transition in the EOS of APR occurs at larger densities than for the quark-to-hadron transition employed earlier which results in a significantly smaller softening of the high density EOS.

In simulations of binary neutron star mergers, the dense matter equation of state (EOS) is required over wide ranges of density and temperature as well as under conditions in which neutrinos are trapped, and the effects of magnetic fields and rotation prevail. Here we assess the status of dense matter theory and point out the successes and limitations of approaches currently in use. A comparative study of the excluded volume (EV) and virial approaches for the $np\alpha$ system using the equation of state of Akmal, Pandharipande and Ravenhall for interacting nucleons is presented in the sub-nuclear density regime. Owing to the excluded volume of the $\alpha$-particles, their mass fraction vanishes in the EV approach below the baryon density 0.1 fm$^{-3}$, whereas it continues to rise due to the predominantly attractive interactions in the virial approach. The EV approach of Lattimer et al. is extended here to include clusters of light nuclei such as d, $^3$H and $^3$He in addition to $\alpha$-particles. Results of the relevant state variables from this development are presented and enable comparisons with related but slightly different approaches in the literature. We also comment on some of the sweet and sour aspects of the supra-nuclear EOS. The extent to which the neutron star gravitational and baryon masses vary due to thermal effects, neutrino trapping, magnetic fields and rotation are summarized from earlier studies in which the effects from each of these sources were considered separately. Increases of about $20\% (\gtrsim 50\%)$ occur for rigid (differential) rotation with comparable increases occurring in the presence of magnetic fields only for fields in excess of $10^{18}$ Gauss. Comparatively smaller changes occur due to thermal effects and neutrino trapping. Some future studies to gain further insight into the outcome of dynamical simulations are suggested.

We present a thermodynamically consistent method by which equations of state based on nonrelativistic potential models can be modified so that they respect causality at high densities, both at zero and finite temperature (entropy). We illustrate the application of the method using the high density phase parametrization of the well known APR model in its pure neutron matter configuration as an example. We also show that, for models with only contact interactions, the adiabatic speed of sound is independent of the temperature in the limit of very large temperature. This feature is approximately valid for models with finite-range interactions as well, insofar as the temperature dependence they introduce to the Landau effective mass is weak. In addition, our study reveals that in first principle nonrelativistic models of hot and dense matter, contributions from higher than two-body interactions must be screened at high density to preserve causality.

We explore the thermal properties of hot and dense matter using a model that reproduces the empirical properties of isospin symmetric and asymmetric bulk nuclear matter, optical model fits to nucleon-nucleus scattering data, heavy-ion flow data in the energy range 0.5-2 GeV/A, and the largest well-measured neutron star mass of 2 $\rm{M}_\odot$. Results of this model which incorporates finite range interactions through Yukawa type forces are contrasted with those of a zero-range Skyrme model that yields nearly identical zero-temperature properties at all densities for symmetric and asymmetric nucleonic matter and the maximum neutron star mass, but fails to account for heavy-ion flow data due to the lack of an appropriate momentum dependence in its mean field. Similarities and differences in the thermal state variables and the specific heats between the two models are highlighted. Checks of our exact numerical calculations are performed from formulas derived in the strongly degenerate and non-degenerate limits. Our studies of the thermal and adiabatic indices, and the speed of sound in hot and dense matter for conditions of relevance to core-collapse supernovae, the thermal evolution of neutron stars from their birth and mergers of compact binary stars reveal that substantial variations begin to occur at sub-saturation densities before asymptotic values are reached at supra-nuclear densities.

We investigate the thermal properties of the potential model equation of state of Akmal, Pandharipande and Ravenhall. This equation of state approximates the microscopic model calculations of Akmal and Pandharipande, which feature a neutral pion condensate. We treat the bulk homogeneous phase for isospin asymmetries ranging from symmetric nuclear matter to pure neutron matter and for temperatures and densities relevant for simulations of core-collapse supernovae, proto-neutron stars, and neutron star mergers. Numerical results of the state variables are compared with those of a typical Skyrme energy density functional with similar properties at nuclear densities, but which differs substantially at supra-nuclear densities. Analytical formulas, which are applicable to non-relativistic potential models such as the equations of state we are considering, are derived for all state variables and their thermodynamic derivatives. A highlight of our work is its focus on thermal response functions in the degenerate and non-degenerate situations, which allow checks of the numerical calculations for arbitrary degeneracy. These functions are sensitive to the density dependent effective masses of neutrons and protons, which determine the thermal properties in all regimes of degeneracy. We develop the "thermal asymmetry free energy" and establish its relation to the more commonly used nuclear symmetry energy. We also explore the role of the pion condensate at supra-nuclear densities and temperatures. Tables of matter properties as functions of baryon density, composition (i.e., proton fraction) and temperature are being produced which are suitable for use in astrophysical simulations of supernovae and neutron stars.