We develop the parquet-diagram summation method for neutron matter interacting via potentials that include spin, tensor, and spin-orbit components. For that purpose, we derive an exact expression for the sum of all ring diagrams in terms of effective local particle-hole interactions involving the above four operators. The parquet equations are closed by deriving the spin-orbit contribution to those particle-hole interactions. We show that many-body correlations screen the bare spin-orbit potential considerably, and the corrections of that screened spin-orbit potential to the other three interaction channels are quite small. We apply our method to the calculation of the response of neutron matter to density and both longitudinal and transverse spin-dependent external fields.
Recent calculations with an effective isospin-dependent contact interaction show the possibility of the crossover from superfluidity of neutron Cooper pairs in the
1S0 pairing channel to Bose-Einstein condensation (BEC) of di-neutron bound states in dilute nuclear matter. The density and spin correlation functions are calculated for a di-neutron condensate in asymmetric nuclear matter with the aim of finding the possible features of the BCS-BEC crossover. It is shown that the zero-momentum transfer spin correlation function satisfies the sum rule at zero temperature. In symmetric nuclear matter, the density correlation function changes sign at low momentum transfer across the BCS-BEC transition, and this feature can be considered as a signature of the crossover. At finite isospin asymmetry, this criterion gives too large a value for the critical asymmetry ${\ensuremath{\alpha}}_{c}^{d}~0.9$, at which the BEC state is quenched. Therefore, it can be trusted for the description of the density-driven BCS-BEC crossover of neutron pairs only at small isospin asymmetry. This result generalizes the conclusion of the study in Phys. Rev. Lett. 95, 090402 (2005), where the change of sign of the density correlation function at low momentum transfer in two-component quantum fermionic atomic gas with the balanced populations of fermions of different species was considered as an unambiguous signature of the BCS-BEC transition.
Low-density neutron matter is relevant to the study of neutron-rich nuclei and neutron star crusts. Unpolarized neutron matter has been studied extensively over a number of decades, while experimental guidance has recently started to emerge from the field of ultracold atomic gases. In this work, we study population-imbalanced neutron matter (possibly relevant to magnetars and to density functionals of nuclei) by applying a quantum Monte Carlo method that has proven successful in the field of cold atoms. We report on the first ab initiosimulations of superfluid low-density polarized neutron matter. For systems with small imbalances, we find a linear dependence of the energy on the polarization, the proportionality coefficient being dependent on the density. We also present results for the momentum and pair distributions of the two fermionic components.