We present an error band on neutron matter properties at finite temperature (finite-T) which comprehends uncertainties on the nuclear interaction, the many-body method convergence, and the thermodynamical consistency of the approach. This study provides nonperturbative predictions for finite-T neutron matter employing chiral interactions which are selected on the basis of their performance in both finite nuclei and infinite matter at zero temperature. Since proper theoretical uncertainties at finite-T are still generally lacking, the band provided here represents a first step towards setting first-principles constraints on thermal aspects of the nuclear matter equation of state.

We exploit the many-body self-consistent Green's function method to analyze finite-temperature properties of infinite nuclear matter and to explore the behavior of the thermal index used to simulate thermal effects in equations of state for astrophysical applications. We show how the thermal index is both density and temperature dependent, unlike often considered, and we provide an error estimate based on our ${\it ab~initio}$ calculations. The inclusion of many-body forces is found to be critical for the density dependence of the thermal index. We also compare our results to a parametrization in terms of the density dependence of the nucleon effective mass. Our study questions the validity of predictions made for the gravitational-wave signal from neutron-star merger simulations with a constant thermal index.

Neutron matter is an ideal laboratory for nuclear interactions derived from chiral effective field theory since all contributions are predicted up to next-to-next-to-next-to-leading order (N$^3$LO) in the chiral expansion. By making use of recent advances in the partial-wave decomposition of three- nucleon (3N) forces, we include for the first time N$^3$LO 3N interactions in many-body perturbation theory (MBPT) up to third order and in self-consistent Green's function theory (SCGF). Using these two complementary many-body frameworks we provide improved predictions for the equation of state of neutron matter at zero temperature and also analyze systematically the many-body convergence for different chiral EFT interactions. Furthermore, we present an extension of the normal-ordering framework to finite temperatures. These developments open the way to improved calculations of neutron-rich matter including estimates of theoretical uncertainties for astrophysical applications.