We generate a new complete equation of state (EOS) of nuclear matter for a wide range of temperatures, densities, and proton fractions ready for use in astrophysical simulations of supernovae and neutron star mergers. Our previous two papers tabulated the EOS at over 180,000 grid points in the temperature range $T$ = 0 to 80 MeV, the density range $n_B$ = 10$^{-8}$ to 1.6 fm$^{-3}$, and the proton fraction range $Y_P$ = 0 to 0.56. In this paper we combine these data points using a suitable interpolation scheme to generate a single equation of state table on a finer grid. This table is thermodynamically consistent and conserves entropy during adiabatic compression tests. We present various thermodynamic quantities and the composition of matter in the new EOS, along with several comparisons with existing EOS tables. Our equation of state table is available for download.

We study the equation of state (EOS) of nuclear matter at subnuclear density in a Virial expansion for a nonideal gas. The gas consists of neutrons, protons, alpha particles, and 8980 species of nuclei with $A \ge 12$ and masses from the finite-range droplet model (FRDM). At very low density, the Virial expansion reduces to nuclear statistical equilibrium. At higher density, the Virial results match smoothly to the relativistic mean field results discussed in our previous paper. We tabulate the resulting EOS at over 73,000 grid points in the temperature range $T$ = 0.158 to 15.8 MeV, the density range $n_B$ = 10$^{-8}$ to 0.1 fm$^{-3}$, and the proton fraction range $Y_P$ = 0.05 to 0.56. In the future we plan to match these low density results to our earlier high density mean field results, and generate a full EOS table for use in supernova and neutron star merger simulations. This Virial EOS is exact in the low density limit.

We calculate the equation of state (EoS) of dense matter, using a relativistic mean field (RMF) model with a density dependent coupling that is a slightly modified form of the original NL3 interaction. For nonuniform nuclear matter we approximate the unit lattice as a spherical Wigner-Seitz cell, wherein the meson mean fields and nucleon Dirac wave functions are solved fully self-consistently. We also calculate uniform nuclear matter for a wide range of temperatures, densities, and proton fractions, and match them to non-uniform matter as the density decreases. The calculations took over 6,000 CPU days in Indiana University's supercomputer clusters. We tabulate the resulting EoS at over 107,000 grid points in the proton fraction range $Y_P$ = 0 to 0.56. For the temperature range $T$ = 0.16 to 15.8 MeV we cover the density range $n_B$ = 10$^{-4}$ to 1.6 fm$^{-3}$; and for the higher temperature range $T$ = 15.8 to 80 MeV we cover the larger density range $n_B$ = 10$^{-8}$ to 1.6 fm$^{-3}$. In the future we plan to study low density, low temperature (T$<$15.8 MeV), nuclear matter using a Virial expansion, and we will match the low density and high density results to generate a complete EoS table for use in astrophysical simulations of supernova and neutron star mergers.