Superfluid dynamics in neutron star crusts

A simple description of superfluid hydrodynamics in the inner crust of a neutron star is given. Particular attention is paid to the effect of the lattice of nuclei on the properties of the superfluid neutrons, and the effects of entrainment, the fact that some fraction of the neutrons are locked to the motion of the protons in nuclei.

💡 Research Summary

The paper presents a concise yet comprehensive formulation of superfluid hydrodynamics in the inner crust of a neutron star, with particular emphasis on how the crystalline lattice of nuclei modifies the behavior of the superfluid neutrons. The authors begin by delineating the physical environment of the inner crust: a periodic array of neutron‑rich nuclei embedded in a sea of dripped neutrons that have undergone Cooper pairing and thus form a superfluid. In this mixed phase, the neutrons are not completely free; the periodic nuclear potential gives rise to a band structure, which in turn modifies the neutrons’ effective mass (m*). This “band‑mass” effect directly influences the inertia of the superfluid component and alters the speed of sound for superfluid phonons.

A central novelty of the work is the explicit inclusion of entrainment – the phenomenon whereby a fraction of the superfluid neutrons is dynamically locked to the motion of the protons bound in the nuclei. The authors introduce an entrainment coefficient, often denoted fₑ, which quantifies the locked‑neutron fraction as a function of density, temperature, and lattice geometry. By employing Skyrme‑Hartree‑Fock calculations of the neutron band structure, they provide realistic estimates of fₑ(ρ) and the corresponding effective mass m*(ρ). These quantities are then incorporated into a two‑fluid hydrodynamic model: one fluid represents the superfluid neutrons, the other the combined lattice‑proton system. The resulting equations of motion contain cross‑terms proportional to fₑ, embodying the momentum exchange between the two components.

Linearizing the coupled equations yields two distinct acoustic modes. The first is a pure superfluid phonon, whose dispersion relation is governed by the superfluid density and the renormalized effective mass. The second is a crustal acoustic mode that involves collective motion of the lattice and the entrained neutrons; its speed depends on the lattice shear modulus, the entrainment coefficient, and the effective mass. The authors show that mode mixing becomes significant near the critical temperature for superfluidity, potentially leaving observable imprints on pulsar timing noise and quasi‑periodic oscillations seen in magnetar flares.

The paper then connects these theoretical results to the long‑standing puzzle of pulsar glitches. In traditional glitch models, the entire superfluid component is assumed to rotate independently of the crust, providing a reservoir of angular momentum that can be suddenly transferred to the solid lattice. However, entrainment reduces the effective angular momentum reservoir to (1 – fₑ) times the nominal superfluid moment of inertia. By inserting realistic values of fₑ (ranging from 0.2 to 0.5 depending on depth), the authors demonstrate that the observed fractional spin‑up ΔΩ/Ω (10⁻⁶–10⁻⁸) can be reproduced with a much smaller superfluid fraction than previously thought. This resolves a tension between glitch amplitudes and constraints from neutron‑star mass–radius measurements.

Another important implication concerns the critical temperature T_c for neutron pairing. The band structure induced by the lattice reduces the pairing gap, thereby lowering T_c. Consequently, at typical neutron‑star interior temperatures (∼10⁸ K), the region where neutrons remain superfluid may be thinner than estimates that ignore lattice effects. This has ramifications for thermal evolution models, as the specific heat and neutrino emissivity of the crust are strongly dependent on the presence of a superfluid.

In the concluding section, the authors acknowledge the limitations of their current framework. The analysis assumes linear perturbations, uniform temperature, and a perfect crystalline lattice without defects. Non‑linear vortex dynamics, temperature gradients, and lattice imperfections (e.g., dislocations, impurities) are not treated but are expected to play a role in realistic glitch events and in the damping of crustal oscillations. The paper calls for future work employing large‑scale quantum Monte‑Carlo simulations, fully relativistic two‑fluid hydrodynamics, and the integration of upcoming observational data from gravitational‑wave detectors, NICER, and next‑generation X‑ray telescopes. Such efforts will refine the entrainment parameters, improve our understanding of neutron‑star interior composition, and ultimately provide a more unified picture of how superfluid dynamics, nuclear lattice structure, and macroscopic astrophysical phenomena intertwine.