Nonperturbative Extraction of the Effective Mass in Neutron Matter

Structure Factors for Hot Neutron Matter from Ab Initio Lattice Simulations with High-Fidelity Chiral Interactions

We present the first ab initio lattice calculations of spin and density correlations in hot neutron matter using high-fidelity interactions at next-to-next-to-next-to-leading order in chiral effective field theory. These correlations have a large impact on neutrino heating and shock revival in core-collapse supernovae and are encapsulated in functions called structure factors. Unfortunately, calculations of structure factors using high-fidelity chiral interactions were well out of reach using existing computational methods. In this Letter, we solve the problem using a computational approach called the rank-one operator (RO) method. The RO method is a general technique with broad applications to simulations of fermionic many-body systems. It solves the problem of exponential scaling of computational effort when using perturbation theory for higher-body operators and higher-order corrections. Using the RO method, we compute the vector and axial static structure factors for hot neutron matter as a function of temperature and density. The ab initio lattice results are in good agreement with virial expansion calculations at low densities but are more reliable at higher densities. Random phase approximation codes used to estimate neutrino opacity in core-collapse supernovae simulations can now be calibrated with ab initio lattice calculations.

BCS-BEC Crossover Effects and Pseudogap in Neutron Matter

Due to the large neutron–neutron scattering length, dilute neutron matter resembles the unitary Fermi gas, which lies half-way in the crossover from the BCS phase of weakly coupled Cooper pairs to the Bose–Einstein condensate of dimers. We discuss crossover effects in analogy with the T-matrix theory used in the physics of ultracold atoms, which we generalize to the case of a non-separable finite-range interaction. A problem of the standard Nozières–Schmitt-Rink approach and different ways to solve it are discussed. It is shown that in the strong-coupling regime, the spectral function exhibits a pseudo-gap at temperatures above the critical temperature Tc. The effect of the correlated density on the density dependence of Tc is found to be rather weak, but a possibly important effect due to the reduced quasiparticle weight is identified.