Neutrons in the cores of neutron stars are expected to form a spin-triplet 3P2^{3}P_{2}–3F2^{3}F_{2} superfluid, and the magnitude, density dependence, and order-parameter structure of this pairing gap remain among the largest open uncertainties in the physics of dense matter. The gap controls the core heat capacity, neutrino emissivity, and cooling behavior of thermally emitting neutron stars, yet present calculations disagree qualitatively: BCS studies with realistic two-nucleon interactions plus chiral three-nucleon forces give gaps of at most a few tenths of an MeV with strong regulator dependence, while the most comprehensive diagrammatic treatment to date suggests the gap may be radically suppressed by many-body correlations. The problem is to determine the 3P2^{3}P_{2}–3F2^{3}F_{2} gap function over neutron-star-core densities with controlled theoretical uncertainties — combining consistent two- and three-nucleon interactions (including induced three-body terms), medium corrections beyond the static first-order self-energy with demonstrated convergence, the magnetic-sublevel gap structure below the critical temperature, and an uncertainty budget spanning effective-field-theory truncation, regulators, and many-body truncation.

Background

A neutron star's interior contains, outside a solid crust, a liquid core of neutron-rich matter at and above nuclear saturation density. Soon after the discovery of pulsars it was proposed that the neutrons in this core become superfluid: like electrons in a superconductor, they form Cooper pairs bound by an attractive component of the nuclear interaction and condense. In the notation of quantum mechanics, pairing states are labeled by spin, orbital angular momentum, and total angular momentum; at core densities the dominant attraction is expected in the spin-triplet (S=1S=1), p-wave (L=1L=1) channel with total angular momentum J=2J=2, denoted 3P2^{3}P_{2}, which is coupled by the tensor force to the 3F2^{3}F_{2} (L=3L=3) wave and must be treated as a coupled two-channel problem. Within Bardeen–Cooper–Schrieffer (BCS) theory — pairing of fermions near a Fermi surface, characterized by a gap Δ\Delta that separates paired quasiparticle excitations — this 3P2^{3}P_{2}–3F2^{3}F_{2} gap opens at densities above roughly one nucleon per several cubic femtometers (Fermi momentum kF≳1 fm−1k_F \gtrsim 1\ \mathrm{fm}^{-1}) and, at the bare two-body level, survives to high density with gaps up to about an MeV, with strong dependence on which realistic nucleon–nucleon (NN) interaction is used.

Modern calculations start from chiral effective field theory (EFT), in which nuclear forces are organized as a systematic expansion in powers of momentum over a breakdown scale Λχ≈500\Lambda_\chi \approx 500–600 MeV600\ \mathrm{MeV}; beyond leading orders this expansion generates three-nucleon (3N) forces, and the truncation order bounds the theoretical error. Interactions are commonly softened with the Similarity Renormalization Group (SRG), a unitary flow that decouples high- and low-momentum modes while preserving two-body scattering data; the flow, however, induces many-body (3N and higher) forces that must be evolved consistently, and any residual dependence of an observable on the flow scale λ\lambda signals missing many-body corrections. In the medium, the pairing interaction and the single-particle spectrum are also modified by self-energy effects (effective mass m∗m^*, quasiparticle residue ZZ) and by medium polarization (screening and antiscreening of the pairing interaction by particle-hole excitations). For the 3P2^{3}P_{2} channel these corrections are individually large: effective-mass and ZZ-factor corrections, self-consistent Green's function studies, and polarization calculations all reduce the gap and shift its closure (the density at which pairing disappears) downward. A further technical limitation is the angle-averaged approximation, which collapses the coupled, angle-dependent gap equation to a scalar one and discards the dependence of the gap on the magnetic quantum number mm of the pair — information needed to decide whether the gap is isotropic or has nodes, which controls neutrino emission and heat transport below the transition temperature TcT_c.

Srinivas and Ramanan (2016) computed angle-averaged 3P2^{3}P_{2}–3F2^{3}F_{2} BCS gaps in pure neutron matter using SRG-evolved chiral NN interactions augmented by the leading N2LO (next-to-next-to-leading order) 3N force added as a density-dependent effective two-body force, plus a static first-order self-energy. It used the λ\lambda-dependence of its results as a diagnostic of missing many-body corrections and found: strong dependence on the 3N regulator cutoff, signaling missing short-range and higher-order 3N effects; that self-energy corrections reduce the gap and shift closure to lower density; and that the trustworthy density window is bounded by the chiral cutoff (conservatively kF≲2.0 fm−1k_F \lesssim 2.0\ \mathrm{fm}^{-1}). It explicitly left open: higher-order and SRG-induced 3N forces; energy-dependent self-energies and higher-order many-body/medium effects beyond first order; and the mm- and jj-resolved gap structure below TcT_c; while its use of resolution-scale and 3N-cutoff dependence as diagnostics of the missing many-body physics motivates, but does not itself constitute, a combined uncertainty budget. Subsequent work has advanced individual pieces but not the whole: Drischler, Krüger, Hebeler, and Schwenk (2017) included subleading N3LO 3N forces in the gap equation for the first time with order-by-order uncertainty bands, but found non-overlapping bands at low chiral orders, strong regulator dependence, and concluded that no final quantitative statement on the 3P2^{3}P_{2}–3F2^{3}F_{2} gap size was possible, while many-body corrections beyond BCS remained unassessed. Papakonstantinou and Clark (2017) surveyed the literature and found extreme variation in the net effect of 3N forces, coupled channels, and mass renormalization on the gap, in some inputs even pointing to absence of a triplet gap. Most significantly, Krotscheck, Papakonstantinou, and Wang (2023, 2024, 2026), summing parquet Feynman diagrams (a simultaneous resummation of ladder and ring diagrams plus important contributions outside the parquet class) with realistic v8v_8-type interactions (phenomenological two-nucleon potentials built from the eight leading spin–isospin operator channels, fitted to nucleon–nucleon scattering data rather than derived from chiral EFT) up to kF=1.8 fm−1k_F = 1.8\ \mathrm{fm}^{-1}, found that many-body correlations strongly reduce the in-medium spin–orbit interaction, leading to a radical — near-complete — suppression of the 3P2^{3}P_{2}–3F2^{3}F_{2} gap and an enhancement of 3P0^{3}P_{0} pairing; the authors note that extending to higher densities requires techniques still to be developed. In parallel, the structure of the assumed 3P2^{3}P_{2} condensate below TcT_c is studied phenomenologically through Ginzburg–Landau theory, spin-polarized phases under magnetic fields (Mizushima, Yasui, Inotani, and Nitta 2021), and vortex physics (Hattori and Sekizawa 2026). Reviews (Sedrakian and Clark 2019; Ramanan and Urban 2021) continue to list the triplet gap among the most uncertain inputs to neutron-star cooling theory. What is still missing is a single determination that combines consistent NN+3N inputs (including induced three-body terms), converged medium corrections, the resolved gap structure below TcT_c, and a quantified uncertainty budget over the relevant density range.

Problem Statement

Determine, with quantified theoretical uncertainties, the neutron 3P2^{3}P_{2}–3F2^{3}F_{2} pairing gap in pure neutron matter (the standard approximation to neutron-star-core matter) as a function of density — its magnitude Δ(kF)\Delta(k_F), its closure density (the Fermi momentum at which pairing disappears), and its magnetic-sublevel structure below the transition temperature TcT_c (the dependence of the gap on the magnetic quantum number mm of the pair — the projection of the pair's total angular momentum J=2J=2 — and on the total angular momentum jj of the individual neutrons, rather than an angle-averaged scalar gap) — from microscopic nuclear interactions, in a calculation that jointly provides: (i) consistent two- and three-nucleon inputs, in the chiral EFT framework including higher-order 3N forces and three-body terms induced by any interaction softening (e.g. SRG evolution), with regulator dependence under control; (ii) medium corrections beyond the static first-order self-energy — energy-dependent self-energies, effective-mass and quasiparticle-strength effects, and screening/polarization (vertex) corrections — with demonstrated convergence of these corrections in the triplet channel at core densities; and (iii) an uncertainty budget that combines EFT truncation errors, NN and 3N regulator dependence, resolution-scale dependence, and many-body truncation errors, over the density range where the underlying interactions are valid and with an honest statement of the reach beyond it. These components jointly serve one objective — a quantitative, uncertainty-quantified triplet-gap determination usable for neutron-star cooling and transport modeling — and answers by any microscopic many-body route (chiral EFT with BCS-type or self-consistent treatments, diagrammatic resummations such as the parquet approach, or ab initio methods) are admissible, provided the full objective — requirements (i)–(iii) together with the resolved gap structure, not a subset of its elements — is discharged; in particular, a route based on interactions that omit the consistent chiral NN+3N inputs of requirement (i) does not by itself discharge the objective.

References

  1. Sarath Srinivas and S. Ramanan, "Triplet pairing in pure neutron matter", Physical Review C 94, 064303 (2016). DOI: 10.1103/PhysRevC.94.064303; arXiv:1606.09053; https://arxiv.org/abs/1606.09053
  2. C. Drischler, T. Krüger, K. Hebeler, and A. Schwenk, "Pairing in neutron matter: New uncertainty estimates and three-body forces", Physical Review C 95, 024302 (2017). DOI: 10.1103/PhysRevC.95.024302; arXiv:1610.05213; https://arxiv.org/abs/1610.05213
  3. P. Papakonstantinou and J. W. Clark, "Three-Nucleon Forces and Triplet Pairing in Neutron Matter", Journal of Low Temperature Physics 189, 361–382 (2017). DOI: 10.1007/s10909-017-1808-9; arXiv:1705.10463 (preprint titled "Three-nucleon forces and superfluidity in neutron matter"); https://arxiv.org/abs/1705.10463
  4. Armen Sedrakian and John W. Clark, "Superfluidity in nuclear systems and neutron stars", European Physical Journal A 55, 167 (2019). DOI: 10.1140/epja/i2019-12863-6; arXiv:1802.00017; https://arxiv.org/abs/1802.00017
  5. M. Urban and S. Ramanan, "Neutron pairing with medium polarization beyond the Landau approximation", Physical Review C 101, 035803 (2020). DOI: 10.1103/PhysRevC.101.035803; arXiv:1911.10772; https://arxiv.org/abs/1911.10772
  6. S. Ramanan and M. Urban, "Pairing in pure neutron matter", European Physical Journal Special Topics 230, 567–577 (2021). DOI: 10.1140/epjs/s11734-021-00008-0; arXiv:2010.09632; https://arxiv.org/abs/2010.09632
  7. Eckhard Krotscheck, Panagiota Papakonstantinou, and Jiawei Wang, "Triplet Pairing in Neutron Matter", The Astrophysical Journal 955, 76 (2023). DOI: 10.3847/1538-4357/acee7c; arXiv:2305.07096; https://arxiv.org/abs/2305.07096
  8. Eckhard Krotscheck, Panagiota Papakonstantinou, and Jiawei Wang, "Variational and parquet-diagram calculations for neutron matter. V. Triplet pairing", Physical Review C 109, 015803 (2024). DOI: 10.1103/PhysRevC.109.015803; arXiv:2308.00873; https://arxiv.org/abs/2308.00873
  9. Panagiota Papakonstantinou, Eckhard Krotscheck, and Jiawei Wang, "Triplet pairing in neutron matter in a comprehensive diagrammatic approach", EPJ Web of Conferences 378, 03002 (2026). DOI: 10.1051/epjconf/202637803002; arXiv:2512.00099; https://arxiv.org/abs/2512.00099
  10. Takeshi Mizushima, Shigehiro Yasui, Daisuke Inotani, and Muneto Nitta, "Spin-polarized phases of 3P2^{3}P_{2} superfluids in neutron stars", Physical Review C 104, 045803 (2021). DOI: 10.1103/PhysRevC.104.045803; arXiv:2108.01256; https://arxiv.org/abs/2108.01256
  11. Tatsuhiro Hattori and Kazuyuki Sekizawa, "Exploring Interplays Between 3P2^{3}P_{2} Neutron Superfluid Vortices and 1S0^{1}S_{0} Proton Fluxtubes in the Outer Core of Neutron Stars", EPJ Web of Conferences 378, 06021 (2026). DOI: 10.1051/epjconf/202637806021; arXiv:2512.22577; https://arxiv.org/abs/2512.22577