Neutron matter is interesting both as an extension of terrestrial nuclear physics and due to its significance for the study of neutron stars. In this work, after some introductory comments on nuclear forces, nuclear ab initio theory, and nuclear phenomenology, we employ two techniques, Quantum Monte Carlo (QMC) and Energy Density Functionals, to practically handle an extended system composed of strongly interacting neutrons. We start by summarizing work on the static response of neutron matter, which considers the impact of external influences on the time-independent system. We then proceed to discuss new results of the energy of quasiparticle excitations in neutron matter, including QMC calculations with chiral or phenomenological nucleon-nucleon interactions. As part of this study, we carefully study the approach of our finite-number computations toward the infinite-system limit.

We carry out non-perturbative calculations of the single-particle excitation spectrum in strongly interacting neutron matter. These are microscopic quantum Monte Carlo computations of many-neutron energies at different densities as well as several distinct excited states. As input, we employ both phenomenological and chiral two- and three-nucleon interactions. We use the single-particle spectrum to extract the effective mass in neutron matter. With a view to systematizing the error involved in this extraction, we carefully assess the impact of finite-size effects on the quasiparticle dispersion relation. We find an effective-mass ratio that drops from 1 as the density is increased. We conclude by connecting our results with the physics of ultracold gases as well as with energy-density functional theories of nuclei and neutron-star matter.

We investigate the problem of periodically modulated strongly interacting neutron matter. We carry out ab initio non-perturbative auxiliary-field diffusion Monte Carlo calculations using an external sinusoidal potential in addition to phenomenological two- and three-nucleon interactions. Several choices for the wave function ansatz are explored and special care is taken to extrapolate finite-sized results to the thermodynamic limit. We perform calculations at various densities as well as at different strengths and periodicities of the one-body potential. Our microscopic results are then used to constrain the isovector term from energy-density functional theories of nuclei at many different densities, while making sure to separate isovector contributions from bulk properties. Lastly, we use our results to extract the static density-density linear response function of neutron matter at different densities. Our findings provide insights into inhomogeneous neutron matter and are related to the physics of neutron star crusts and neutron-rich nuclei.

We generalize the problem of strongly interacting neutron matter by adding a periodic external modulation. This allows us to study from first principles a neutron system that is extended and inhomogeneous, with connections to the physics of both neutron-star crusts and neutron-rich nuclei. We carry out fully non-perturbative microscopic Quantum Monte Carlo calculations of the energy of neutron matter at different densities, as well as different strengths and periodicities of the external potential. In order to remove systematic errors, we examine finite-size effects and the impact of the wave function ansatz. We also make contact with energy-density functional theories of nuclei and disentangle isovector gradient contributions from bulk properties. Finally, we calculate the static density-density linear response function of neutron matter and compare it with the response of other physical systems.

We survey the current status of understanding of pairing and superfluidity of neutrons and protons in neutron stars from a theoretical perspective, with emphasis on basic physical properties. During the past two decades, the blossoming of the field of ultracold atomic gases and the development of quantum Monte Carlo methods for solving the many-body problem have been two important sources of inspiration, and we shall describe how these have given insight into neutron pairing gaps. The equilibrium properties and collective oscillations of the inner crust of neutron stars, where neutrons paired in a $^1$S$_0$ state coexist with a lattice of neutron-rich nuclei, are also described. While pairing gaps are well understood at densities less than one tenth of the nuclear saturation density, significant uncertainties exist at higher densities due to the complicated nature of nucleon-nucleon interactions, the difficulty of solving the many-body problem under these conditions, and the increasing importance of many-nucleon interactions. We also touch more briefly on the subject of pairing of neutrons in other angular momentum states, specifically the $^3$P$_2$ state, as well as pairing of protons.

Neutron matter is an intriguing nuclear system with multiple connections to other areas of physics. Considerable progress has been made over the last two decades in exploring the properties of pure neutron fluids. Here we begin by reviewing work done to explore the behavior of very low density neutron matter, which forms a strongly paired superfluid and is thus similar to cold Fermi atoms, though at energy scales differing by many orders of magnitude. We then increase the density, discussing work that ties the study of neutron matter with the determination of the properties of neutron-rich nuclei and neutron-star crusts. After this, we review the impact neutron matter at even higher densities has on the mass-radius relation of neutron stars, thereby making contact with astrophysical observations.

Neutron matter presents a unique system in chiral effective field theory (EFT), because all many-body forces among neutrons are predicted to next-to-next-to-next-to-leading order (N3LO). We discuss perturbative and first Quantum Monte Carlo (QMC) calculations of neutron matter with chiral EFT interactions and their astrophysical impact for the equation of state and neutron stars.