We present a much improved equation of state for neutron star matter, QHC19, with a smooth crossover from the hadronic regime at lower densities to the quark regime at higher densities. We now use the Togashi et al.~equation of state (Togashi:2017), a generalization of the Akmal-Pandharipande-Ravenhall equation of state of uniform nuclear matter, in the entire hadronic regime; the Togashi equation of state consistently describes non-uniform as well as uniform matter, and matter at beta equilibrium without the need for an interpolation between pure neutron and symmetric nuclear matter. We describe the quark matter regime at higher densities with the Nambu--Jona--Lasinio model, now identifying tight constraints on the phenomenological universal vector repulsion between quarks and the pairing interaction between quarks arising from the requirements of thermodynamic stability and causal propagation of sound. The resultant neutron star properties agree very well with the inferences of the LIGO/Virgo collaboration, from GW170817, of the pressure vs. baryon density, neutron star radii, and tidal deformabilities. The maximum neutron star mass allowed by QHC19 is 2.35 $M_\odot$, consistent with all neutron star mass determinations.

We perform an extensive study of the influence of nuclear weak interactions on core-collapse supernovae (CCSNe), paying particular attention to consistency between nuclear abundances in the equation of state (EOS) and nuclear weak interactions. We compute properties of uniform matter based on the variational method. For inhomogeneous nuclear matter, we take a full ensemble of nuclei into account with various finite-density and thermal effects and directly use the nuclear abundances to compute nuclear weak interaction rates. To quantify the impact of a consistent treatment of nuclear abundances on CCSN dynamics, we carry out spherically symmetric CCSN simulations with full Boltzmann neutrino transport, systematically changing the treatment of weak interactions, EOSs, and progenitor models. We find that the inconsistent treatment of nuclear abundances between the EOS and weak interaction rates weakens the EOS dependence of both the dynamics and neutrino signals. We also test the validity of two artificial prescriptions for weak interactions of light nuclei and find that both prescriptions affect the dynamics. Furthermore, there are differences in neutrino luminosities by ~10% and in average neutrino energies by 0.25-1 MeV from those of the fiducial model. We also find that the neutronization burst neutrino signal depends on the progenitor more strongly than on the EOS, preventing a detection of this signal from constraining the EOS.

A proto-neutron star (PNS) is a newly formed compact object in a core collapse supernova. In this Letter, the neutrino emission from the cooling process of a PNS is investigated using two types of nuclear equation of state (EOS). It is found that the neutrino signal is mainly determined by the high-density EOS. The neutrino luminosity and mean energy are higher and the cooling time scale is longer for the softer EOS. Meanwhile, the neutrino mean energy and the cooling time scale are also affected by the low-density EOS because of the difference in the population of heavy nuclei. Heavy nuclei have a large scattering cross section with neutrinos owing to the coherent effects and act as thermal insulation near the surface of a PNS. The neutrino mean energy is higher and the cooling time scale is longer for an EOS with a large symmetry energy at low densities, namely a small density derivative coefficient of the symmetry energy, $L$.

The existence of neutron stars with $2M_\odot$ requires the strong stiffness of the equation of state (EoS) of neutron-star matter. We introduce a multi-pomeron exchange potential (MPP) working universally among 3- and 4-baryons to stiffen the EoS. Its strength is restricted by analyzing the nucleus-nucleus scattering with the G-matrix folding model. The EoSs are derived using the Brueckner-Hartree-Fock (BHF) and the cluster variational method (CVM) with the nuclear interactions ESC and AV18. The mass-radius relations are derived by solving the Tolmann-Oppenheimer-Volkoff (TOV) equation, where the maximum masses over $2M_\odot$ are obtained on the basis of the terrestrial data. Neutron-star radii $R$ at a typical mass $1.5M_\odot$ are predicted to be $12.3\!\sim\!13.0$ km. The uncertainty of calculated radii is mainly from the ratio of 3- and 4-pomeron coupling constants, which cannot be fixed by any terrestrial experiment. Though values of $R(1.5M_\odot)$ are not influenced by hyperon-mixing effects, finely-observed values for them indicate degrees of EoS softening by hyperon mixing in the region of $M\!\sim\!2M_\odot$. If $R(1.5M_\odot)$ is less than about 12.4 km, the softening of EoS by hyperon mixing has to be weak. Useful information can be expected by the space mission NICER offering precise measurements for neutron-star radii within $\pm 5\%$.

We have constructed a nuclear equation of state (EOS) that includes a full nuclear ensemble for use in core-collapse supernova simulations. It is based on the EOS for uniform nuclear matter that two of the authors derived recently, applying a variational method to realistic two- and there-body nuclear forces. We utilize an extended liquid drop model of heavy nuclei and a quantum-theoretical mass evaluation for light nuclei. In addition to realistic nuclear forces, the inclusion of in-medium effects on the full ensemble of nuclei makes the new EOS one of the most realistic EOS's for supernova simulations. We make comparisons with the FYSS EOS, which is based on the same formulation for the nuclear ensemble but adopts the relativistic mean field (RMF) theory with the TM1 parameter set for uniform nuclear matter. The new EOS is softer than the FYSS EOS around and above nuclear saturation densities. We find that neutron-rich nuclei with small mass numbers are more abundant in the new EOS than in the FYSS EOS because of the larger saturation densities and smaller symmetry energy of nuclei in the former. We apply the two EOS's to 1D supernova simulations and find that the new EOS gives lower electron fractions and higher temperatures in the collapse phase owing to the smaller symmetry energy. As a result, the inner core has smaller masses for the new EOS. It is more compact, on the other hand, due to the softness of the new EOS and bounces at higher densities. The ensuing outward propagations of the shock wave in the outer core are very similar in the two simulations, which may be an artifact, though, caused by the use of the same tabulated electron capture rates for heavy nuclei ignoring differences in the nuclear composition between the two EOS's in these computations.

A new table of the nuclear equation of state (EOS) based on realistic nuclear potentials is constructed for core-collapse supernova numerical simulations. Adopting the EOS of uniform nuclear matter constructed by two of the present authors with the cluster variational method starting from the Argonne v18 and Urbana IX nuclear potentials, the Thomas-Fermi calculation is performed to obtain the minimized free energy of a Wigner-Seitz cell in non-uniform nuclear matter. As a preparation for the Thomas-Fermi calculation, the EOS of uniform nuclear matter is modified so as to remove the effects of deuteron cluster formation in uniform matter at low densities. Mixing of alpha particles is also taken into account following the procedure used by Shen et al. (1998, 2011). The critical densities with respect to the phase transition from non-uniform to uniform phase with the present EOS are slightly higher than those with the Shen EOS at small proton fractions. The critical temperature with respect to the liquid-gas phase transition decreases with the proton fraction in a more gradual manner than in the Shen EOS. Furthermore, the mass and proton numbers of nuclides appearing in non-uniform nuclear matter with small proton fractions are larger than those of the Shen EOS. These results are consequences of the fact that the density derivative coefficient of the symmetry energy of our EOS is smaller than that of the Shen EOS.

We investigate the effects of the odd-state part of bare $\Lambda \Lambda$ interactions on the structure of neutron stars (NSs) by constructing equations of state (EOSs) for uniform nuclear matter containing $\Lambda$ and $\Sigma^-$ hyperons with use of the cluster variational method. The isoscalar part of the Argonne v18 two-nucleon potential and the Urbana IX three-nucleon potential are employed as the interactions between nucleons, whereas, as the bare $\Lambda N$ and even-state $\Lambda \Lambda$ interactions, two-body central potentials that are determined so as to reproduce the experimental data on single- and double-$\Lambda$ hypernuclei are adopted. In addition, the $\Sigma^- N$ interaction is constructed so as to reproduce the empirical single-particle potential of $\Sigma^-$ in symmetric nuclear matter. Since the odd-state part of the $\Lambda \Lambda$ interaction is not known owing to lack of experimental data, we construct four EOSs of hyperonic nuclear matter, each with a different odd-state part of the $\Lambda \Lambda$ interaction. The EOS obtained for NS matter becomes stiffer as the odd-state $\Lambda \Lambda$ interaction becomes more repulsive, and correspondingly the maximum mass of NSs increases. It is interesting that the onset density of $\Sigma^-$ depends strongly on the repulsion of the odd-state $\Lambda \Lambda$ interaction. Furthermore, we take into account the three-baryon repulsive force to obtain results that are consistent with observational data on heavy NSs.

The equation of state (EOS) for hot asymmetric nuclear matter which is constructed with the variational method starting from the Argonne v18 and Urbana IX nuclear forces is applied to spherically symmetric core-collapse supernovae (SNe). We first investigate the EOS of isentropic beta-stable SN matter, and find that the matter with the variational EOS is more neutron-rich than that with the Shen EOS. Using the variational EOS for uniform matter supplemented by the Shen EOS of non-uniform matter at low densities, we perform general-relativistic spherically symmetric simulations of core-collapse SNe with and without neutrino transfer, starting from a presupernova model of 15 solar mass. In the adiabatic simulation without neutrino transfer, the explosion is successful, and the explosion energy with the variational EOS is larger than that with the Shen EOS. In the case of the simulation with neutrino transfer, the shock wave stalls and then the explosion fails, as in other spherically symmetric simulations. The inner core with the variational EOS is more compact than that with the Shen EOS, due to the relative softness of the variational EOS. This implies that the variational EOS is more advantageous for SN explosions than the Shen EOS.

An equation of state (EOS) for uniform asymmetric nuclear matter (ANM) is constructed at zero and finite temperatures by the variational method starting from the nuclear Hamiltonian that is composed of the Argonne v18 and Urbana IX potentials. At zero temperature, the two-body energy is calculated with the Jastrow wave function in the two-body cluster approximation which is supplemented by Mayer's condition and the healing-distance condition so as to reproduce the result by Akmal, Pandharipande and Ravenhall. The energy caused by the three-body force is treated somewhat phenomenologically so that the total energy reproduces the empirical saturation conditions. The masses and radii of neutron stars obtained with the EOS are consistent with recent observational data. At finite temperatures, thermodynamic quantities such as free energy, internal energy, entropy, pressure and chemical potentials are calculated with an extension of the method by Schmidt and Pandharipande. The validity of the frozen-correlation approximation employed in this work is confirmed as compared with the result of the fully minimized calculation. The quadratic proton-fraction-dependence of the energy of ANM is confirmed at zero temperature, whereas the free energy of ANM deviates from the quadratic proton-fraction-dependence markedly at finite temperatures. The obtained EOS of ANM will be an important ingredient of a new nuclear EOS for supernova numerical simulations.