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.

An equation of state (EOS) for uniform nuclear matter is constructed at zero and finite temperatures with the variational method starting from the realistic nuclear Hamiltonian composed of the Argonne V18 and UIX potentials. The energy is evaluated in the two-body cluster approximation with the three-body-force contribution treated phenomenologically so as to reproduce the empirical saturation conditions. The obtained energies for symmetric nuclear matter and neutron matter at zero temperature are in fair agreement with those by Akmal, Pandharipande and Ravenhall, and the maximum mass of the neutron star is 2.2 Msolar. At finite temperatures, a variational method by Schmidt and Pandharipande is employed to evaluate the free energy, which is used to derive various thermodynamic quantities of nuclear matter necessary for supernova simulations. The result of this variational method at finite temperatures is found to be self-consistent.

A new scenario for neutron-star cooling is proposed, based on the correspondence between pion condensation, occurring in neutron matter due to critical spin-isospin fluctuations, and the metal-insulator phase transition in a two-dimensional electron gas. Beyond the threshold density for pion condensation, where neutron-star matter loses its spatial homogeneity, the neutron single-particle spectrum acquires an insulating gap that quenches neutron contributions to neutrino-production reactions and to the star's specific heat. In the liquid phase at densities below the transition point, spin-isospin fluctuations are found to play dual roles. On the one hand, they lead to a multi-sheeted neutron Fermi surface that extends to low momenta, thereby activating the normally forbidden direct-Urca cooling mechanism; on the other, they amplify the nodeless $P$-wave neutron superfluid gap while suppressing $S$-wave pairing. In this picture, lighter stars without a pion-condensed core experience slow cooling, while enhanced cooling occurs in heavier stars through direct-Urca emission from a narrow shell of the interior.