⚛️ nuclear theory

Ab initio Green's function theory of superfluid neutron matter

This paper presents a unified ab initio description of superfluid neutron matter across a wide density range using the nonperturbative Gorkov algebraic diagrammatic construction method, successfully quantifying pairing gaps and identifying microscopic signatures of superfluidity while assessing the sensitivity of these results to many-body truncations.

Original authors: F. Marino, W. G. Jiang, C. Barbieri, G. Colò

Published 2026-10-05

📖 5 min read🧠 Deep dive

CC BY 4.0

Original authors: F. Marino, W. G. Jiang, C. Barbieri, G. Colò

Original paper licensed under CC BY 4.0 (http://creativecommons.org/licenses/by/4.0/). ✨ This is an AI-generated explanation of the paper below. It is not written or endorsed by the authors. For technical accuracy, refer to the original paper. Read full disclaimer

Deep inside the crushing heart of a neutron star, where matter is squeezed to densities found nowhere else in the universe, a strange transformation occurs. These stars are not solid balls of rock or gas, but vast, exotic fluids made almost entirely of neutrons. In the outer layers, these neutrons behave like a normal liquid, but as you move deeper into the inner crust, the environment changes. The density drops, and the neutrons begin to pair up, forming a superfluid. This is a state of matter where particles flow without any friction at all, a phenomenon that dictates how these stars cool down over millions of years and how they spin, occasionally causing sudden glitches in their rotation. Understanding this superfluid state is crucial for decoding the life cycle of these stellar remnants, yet for decades, scientists have struggled to describe it with precision. The difficulty lies in the sheer complexity of the interactions between these neutrons; they influence one another in ways that are too tangled for simple, standard calculations to untangle.

A team of researchers has now stepped into this complex landscape with a new, more powerful approach to map the behavior of this superfluid neutron matter. Instead of relying on older methods that make simplifying guesses about how these particles interact, they employed a sophisticated computational framework known as the Gorkov algebraic diagrammatic construction. Think of this method as a way to track the energy and movement of every single neutron in the system, accounting for the fact that they are constantly bumping into and influencing one another in a crowded, chaotic environment. By running these simulations across a wide range of densities, from the thin gas of the inner crust to the denser core, the team was able to generate a unified picture of how the neutrons behave, bridging the gap between the sparse outer layers and the dense interior.

The results of these simulations reveal a clear and distinct signature of superfluidity. In the normal, non-superfluid state, neutrons can be excited to higher energy levels with very little effort, creating a smooth transition in their energy spectrum. However, when the system enters the superfluid phase, a gap opens up in this spectrum. It is as if a barrier forms, requiring a specific, minimum amount of energy to break the pairs of neutrons that have formed. The researchers were able to measure the size of this energy gap, finding that it reaches a maximum of about 2 million electron volts in the density range where superfluidity is strongest. This is a significant finding because it confirms that the pairing is robust, but it also shows that the gap is smaller than what older, simpler theories had predicted.

One of the most important contributions of this work is the rigorous assessment of how much confidence we can place in these numbers. Previous studies often used approximations that were difficult to control, leaving scientists unsure if their results were accurate or just artifacts of the math used. This new study explicitly tested different levels of complexity in their calculations. They found that while the basic predictions hold up, the finer details of the interaction between neutrons do matter. Specifically, they discovered that the dynamic, fleeting interactions between particles tend to suppress the size of the pairing gap compared to the simpler, static models. By comparing their most advanced simulations against slightly less complex versions, they were able to estimate the uncertainty in their results, showing that the predictions are stable and reliable up to a certain density.

The team also explored how different theoretical models of the forces between neutrons affect the outcome. They tested several variations of the nuclear force, all derived from modern theories of particle physics, and found that while the details of the force change the exact value of the gap at higher densities, the predictions remain remarkably consistent in the low-density region where the superfluid forms. This consistency suggests that the fundamental physics of the pairing is well-understood in this regime. Furthermore, the study highlighted a critical limitation in other popular methods used to study neutron stars. When the researchers tried to use a different, widely used technique to simulate the same low-density environment, the calculations failed to converge, essentially breaking down. In contrast, their new method continued to work smoothly, capturing the subtle changes in the energy of the system that signal the onset of superfluidity.

Ultimately, this work provides a solid, quantitative foundation for understanding the interior of neutron stars. By pinning down the energy required to break neutron pairs and quantifying the uncertainties in the calculation, the researchers have removed a major source of guesswork from the field. This clarity is essential for interpreting observations of neutron stars, such as how quickly they cool or how they respond to sudden changes in rotation. With a more accurate map of the superfluid state, scientists can now better connect the microscopic world of individual neutrons to the macroscopic behavior of these cosmic giants, paving the way for a deeper understanding of the most extreme matter in the universe.

Technical Summary: Ab initio Green's function theory of superfluid neutron matter

Problem Statement
The inner crust of neutron stars is characterized by a regime where neutron-rich nuclei coexist with a dilute neutron gas at nearly zero temperature. In this low-density environment, neutron-neutron interactions in the attractive $^1S_0$ channel drive the formation of Cooper pairs, leading to superfluidity. This phenomenon significantly influences macroscopic astrophysical observables, including star cooling rates and rotational dynamics (e.g., pulsar glitches). However, a precise quantitative understanding of this system has been historically hindered by uncontrolled many-body approximations. Standard descriptions rely on BCS mean-field theory (or Hartree-Fock-Bogoliubov, HFB), which assumes weak interactions—a premise questionable in strongly correlated nuclear systems. While beyond-mean-field corrections are known to suppress pairing gaps relative to BCS predictions, previous calculations lacked rigorous assessments of theoretical uncertainties, particularly regarding the sensitivity of results to many-body truncations and the specific choice of chiral Hamiltonians.

Methodology
The authors employ a unified ab initio framework combining the nonperturbative Gorkov algebraic diagrammatic construction (ADC) method with the Gorkov Green's function formalism. This approach allows for the calculation of the one-body spectral function $S_{11}(k, \omega)$ , which resolves both momentum ( $k$ ) and energy ( $\omega$ ) dependencies.

Formalism: The study utilizes the Gorkov Self-consistent Green's function (SCGF) formalism at zero temperature. The spectral function is defined via quasiparticle amplitudes ( $U^q_k, V^q_k$ ) and energies ( $\hbar\omega_q$ ), determined by solving a self-consistent Gorkov eigenvalue problem involving the normal self-energy $\Sigma_{11}$ and the anomalous self-energy $\Sigma_{12}$ .
Approximations: The dynamical self-energy $\tilde{\Sigma}_{11}(k, \omega)$ is computed within the ADC hierarchy. The primary results utilize the third-order ADC approximation [ADC(3)], which includes dominant dynamical correlations (both ladder and ring diagrams) in the normal channel. For comparison, the authors also employ the Hartree-Fock-Bogoliubov (HFB) approximation (static self-energies only), the second-order ADC [ADC(2)], and the ADC(ld,2) truncation (an all-order resummation of ladder diagrams).
Computational Setup: Simulations are performed for homogeneous neutron matter using $N=66$ particles in a cubic box with special-point twisted-angle boundary conditions (sp-TABCs). This boundary condition is critical for resolving momentum-dependent observables like the pairing gap, which standard periodic boundary conditions fail to capture accurately due to coarse momentum grids.
Interactions: The study tests three different chiral potentials at next-to-next-to-leading order (NNLO): $\Delta$ NNLO $_{GO}$ (394), $\Delta$ NNLO $_{GO}$ (450), and NNLO $_{sat}$ (450).

Key Results

Microscopic Signatures of Superfluidity: The spectral functions reveal a distinct qualitative difference between the normal and superfluid phases. In the normal phase, quasiparticle bands cross the chemical potential near the Fermi momentum ( $k_F$ ). In the superfluid phase, these bands repel each other, creating a finite gap in the spectral strength around the chemical potential. This gap is quantified by fitting the dispersion branches to a Bogoliubov-like formula.
Suppression of Pairing Gaps: Comparing ADC(3) results with HFB reveals that dynamical correlations significantly suppress the pairing gap. For the interactions studied, HFB predicts a maximum gap of approximately 2.5 MeV, whereas ADC(3) reduces this to about 2.0 MeV. This quenching is attributed to the energy-dependent contribution to the effective mass (the quasiparticle weight $Z < 1$ ) arising from the dynamical self-energy.
Convergence and Uncertainty: The study provides the first quantitative estimate of uncertainties associated with many-body truncation.
- At low Fermi momenta ( $k_F \lesssim 0.7$ fm $^{-1}$ ), ADC(3) and ADC(ld,2) results overlap closely, indicating that ladder diagrams are the crucial ingredient for convergence in this regime.
- Theoretical uncertainties (estimated as the deviation between ADC(3) and ADC(ld,2)) are negligible up to $k_F \approx 0.8$ fm $^{-1}$ , where predictions across different chiral potentials show remarkable consistency.
- Above the gap maximum, interaction dependence becomes pronounced due to finite-range effects and three-nucleon contributions, leading to larger theoretical uncertainties.
Equation of State (EOS) and Method Breakdown: The paper investigates the Equation of State (EOS) in the superfluid regime. While HFB and Hartree-Fock (HF) results diverge below $k_F \approx 0.8$ fm $^{-1}$ (with HFB being energetically favored), correlated methods show a more complex behavior. Coupled-cluster calculations at the doubles truncation level (CCD) fail to converge below $k_F \approx 0.7$ fm $^{-1}$ . Attempts to extrapolate CCD results using a density-continuation emulator fail to reproduce the subtle change in EOS concavity observed in the ADC(3) results. The discrepancy between ADC(3) and normal-state CCD reaches nearly 20% of the total energy in the region $k_F \approx 0.4–0.6$ fm $^{-1}$ , coinciding with the region of maximal pairing gaps.

Significance and Claims
The authors claim that this work provides a robust, unified ab initio description of neutron matter spanning from the dilute regime to saturation density. The primary significance lies in:

Uncontrolled Approximations: Moving beyond uncontrolled mean-field theories by employing a nonperturbative Green's function method that systematically accounts for dynamical correlations.
Uncertainty Quantification: Explicitly estimating the sensitivity of pairing gaps to many-body truncation, demonstrating that gaps are well-constrained up to $k_F \approx 0.8$ fm $^{-1}$ once ladder diagrams are included.
EOS Concavity: Identifying a distinct qualitative change in the concavity of the EOS induced by superfluidity, a feature that normal-state frameworks (like standard CCD) cannot capture, even with extrapolation techniques.
Future Utility: Establishing a framework capable of exploring electroweak processes relevant to neutron stars, specifically neutrino emission rates and thermal evolution, which depend critically on the accurate description of pairing and dynamical correlations.

The paper concludes that the Gorkov-ADC method offers a single, robust framework for understanding neutron matter across density regimes, paving the way for more accurate modeling of neutron star crusts and their thermal evolution.

Technical Summary: Ab initio Green's function theory of superfluid neutron matter

More like this