Superfluid neutron matter with a twist

TL;DR

This paper introduces a novel application of twisted boundary conditions in the context of neutron superfluidity, demonstrating that twist-averaging significantly reduces finite-size effects in calculations relevant for neutron star physics.

Contribution

It pioneers the use of twisted boundary conditions for pairing calculations in superfluid neutron matter, improving the accuracy of finite system extrapolations.

Findings

Twist-averaging reduces finite-size effects to within 2% for systems with N≥50.

First application of twisted boundary conditions in pairing calculations.

Enhanced accuracy in modeling neutron star superfluid properties.

Abstract

Superfluid neutron matter is a key ingredient in the composition of neutron stars. The physics of the inner crust is largely dependent on that of its $S$-wave neutron superfluid which has made its presence known through pulsar glitches and modifications on the neutron star cooling. Moreover, with recent gravitational-wave observations of neutron star mergers, the need for an equation of state for the matter of these compact stars is further accentuated and a model-independent treatment of neutron superfluidity is important. \textit{Ab initio} techniques developed for finite systems can be guided to perform extrapolations to the thermodynamic limit and attain this model-independent extraction of various quantities of infinite superfluid neutron matter. To inform such an extrapolation scheme, we performed calculations of the neutron $^1S_0$ pairing gap using the model-independent odd-even staggering in the context of the particle-conserving, projected BCS theory under twisted boundary conditions. While the practice of twisted boundary conditions is standard in solid state physics and has been used repeatedly in the past to reduce finite-size effects, this is the first time it is employed in the context of pairing. We find that a twist-averaging approach results in a substantial reduction of the finite-size effects, bringing systems with $N\gtrapprox 50$ within a $2\%$ error margin from the infinite system. This can significantly reduce extrapolation-related errors in the extraction of superfluid neutron matter quantities.