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.

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.