Gorkov algebraic diagrammatic construction for infinite nuclear matter

This paper introduces a novel many-body truncation for Gorkov self-consistent Green's function theory that combines first-order pairing correlations with third-order particle-number-conserving dynamical correlations to accurately predict the equation of state and spectral properties of infinite nuclear matter using modern chiral effective field theory Hamiltonians.

The Gist

Imagine you are trying to understand the behavior of a massive, invisible crowd of people (nucleons) packed tightly together in a giant, invisible box. This "crowd" is nuclear matter, the stuff that makes up the cores of neutron stars and the hearts of heavy atoms.

Physicists have been trying to write a "rulebook" for how this crowd behaves for decades. The problem is that these particles don't just walk around; they are constantly grabbing, pushing, and dancing with each other in complex, chaotic ways. To predict how this crowd acts, you need a super-accurate mathematical model.

This paper introduces a new, smarter way to build that model. Here is the breakdown using simple analogies:

1. The Problem: The "Two-Step Dance"

In the world of nuclear physics, particles do two main things:

• The General Crowd: They move around, bump into each other, and create pressure (like a busy dance floor).
• The Pairs: Sometimes, they pair up and dance in perfect sync (this is called pairing or superfluidity).

Old methods tried to handle the "General Crowd" and the "Pairs" separately, or they tried to do everything at once but got stuck in a mathematical mess. It was like trying to film a chaotic mosh pit and a synchronized ballet at the same time with a shaky camera. The result was often blurry or unstable.

2. The Solution: A Hybrid Camera System

The authors (Marino, Barbieri, and Colò) built a hybrid camera system. They didn't try to film the whole chaotic mess at once. Instead, they split the job:

• The "Pairing" Lens (First Order): They used a simple, stable lens to focus just on the dancing pairs. They treated this part simply but accurately. Think of this as taking a clear, steady photo of the couples dancing.
• The "Crowd" Lens (Third Order): For the chaotic mosh pit (the rest of the particles), they used a high-tech, super-detailed lens called ADC(3). This lens captures the complex, long-range interactions and the "ripples" in the crowd.

The Magic Trick: The paper's biggest innovation is a "translator" that allows these two different lenses to talk to each other. They created a way to describe the complex, pairing-heavy system using a simpler, "particle-number-conserving" reference state. It's like describing a complex jazz improvisation by first writing down the basic sheet music, then adding the improvisation on top, rather than trying to write the whole jazz song from scratch every time.

3. The "Optimized Reference" (The Perfect Map)

To make this work, they had to create a "perfect map" of the crowd's average behavior, which they call the Optimized Reference State (OpRS).

Imagine you are trying to predict traffic in a city. You could guess the average speed, but that's often wrong. Instead, this method calculates the best possible average speed based on the current traffic conditions, then uses that as a baseline to calculate the specific jams and accidents (the complex correlations).

They tested four different ways to draw this "map" and found one that was the most stable and accurate, ensuring their predictions didn't wobble or crash.

4. The Results: Predicting the Universe

Using this new method, they simulated Infinite Nuclear Matter (a theoretical, endless block of neutrons and protons) at zero temperature.

• The Equation of State (EOS): This is the "pressure vs. density" rulebook. It tells us how hard it is to squeeze the matter. Their results matched other top-tier methods perfectly, giving them confidence that their "rulebook" is correct.
• The "Ghost" Effects: They didn't just calculate the total weight of the crowd; they looked at the momentum distribution. This is like asking, "How many people are walking slowly vs. running fast?"
  • In a perfect, boring world (Hartree-Fock theory), everyone walks at a set speed.
  • In their realistic simulation, they found that because of the complex interactions, some "slow walkers" are actually running fast, and some "runners" are slowing down. This "smearing" of speeds is a sign of strong quantum correlations.

5. Why This Matters

Why do we care about a theoretical box of infinite particles?

• Neutron Stars: These are the densest objects in the universe. To understand how they cool down, how they spin, and how they crash into each other (creating gravitational waves), we need to know exactly how nuclear matter behaves under extreme pressure.
• Better Tools: This method is faster and more stable than previous attempts. It allows scientists to use the most modern, accurate forces (derived from Chiral Effective Field Theory) without the math breaking down.

The Bottom Line

The authors built a hybrid mathematical engine that combines a simple, stable view of particle pairing with a high-definition view of particle chaos. By using a clever "translator" (the optimized reference state), they managed to simulate infinite nuclear matter with unprecedented accuracy.

This isn't just about numbers; it's about giving astrophysicists a better map to understand the most extreme environments in our universe, from the cooling of neutron stars to the explosions of supernovas. They didn't just solve a math problem; they polished the lens through which we view the cosmos.

Technical Summary

1. Problem Statement

The accurate ab initio description of infinite nuclear matter, particularly in the presence of pairing correlations (superfluidity) and at low densities, remains a significant challenge in theoretical nuclear physics.

• The Pairing Instability: Standard particle-number-conserving methods (like Dyson Self-Consistent Green's Function or Coupled Cluster) struggle at low densities where the system undergoes a superfluid phase transition. The particle-hole excitation gap vanishes, leading to numerical instabilities.
• The Gorkov Challenge: While Gorkov Green's Function theory naturally handles pairing by breaking particle number symmetry in the reference state, extending high-order many-body approximations (like third-order Algebraic Diagrammatic Construction, ADC(3)) within the full Gorkov framework is computationally expensive and complex due to the proliferation of anomalous diagrams.
• The Gap: Existing approaches often rely on finite-temperature formalisms to avoid pairing instabilities (extrapolating to $T=0$) or use lower-order truncations (ADC(2)) that lack the precision required for modern chiral effective field theory ($\chi$EFT) interactions. There is a need for a method that combines the stability of Gorkov theory for pairing with the high-order accuracy of Dyson-ADC(3) for dynamical correlations.

2. Methodology

The authors propose a novel hybrid Gorkov-Dyson scheme that treats pairing and dynamical correlations at different levels of approximation within a unified framework.

• Hybrid Truncation:

  • Pairing (Static): Pairing correlations are handled at the first order (static level) using the full Gorkov formalism. This includes a static pairing field ($\Delta$) and a static mean field ($\Lambda$) derived from the dressed propagator.
  • Dynamical Correlations: The dynamic part of the self-energy is treated using the Dyson-ADC(3) scheme. This assumes that the dynamical dependence of the anomalous (pairing) self-energy is negligible, while the normal (particle-number conserving) self-energy is computed up to third order.
  • Justification: Pairing is predominantly a Fermi-surface phenomenon, while strong many-body dynamical correlations (fragmentation of single-particle strength) occur over a broader energy range. Neglecting the frequency dependence of the anomalous self-energy is a well-justified approximation for ground-state energies.

• Optimized Reference State (OpRS):

  • To bridge the Gorkov propagator (which breaks particle number) and the Dyson-ADC(3) calculation (which requires particle-number conservation), the authors introduce an Optimized Reference State.
  • The dressed Gorkov propagator is approximated by a particle-number-conserving propagator ($g^{OpRS}$) with the same number of poles as a mean-field propagator.
  • The poles of this reference state are determined by matching moments of the spectral function (specifically the "Centroid" prescription, $p=+1$, combined with a Fermi momentum constraint, $k_F$). This ensures the reference state preserves the chemical potential and particle number parity while capturing essential correlation effects.

• Corrections and Extensions:

  • Non-Skeleton Diagrams: The method includes non-skeleton diagrams to correct for the partial self-consistency introduced by the OpRS. These diagrams effectively revert the interaction energies in the ADC matrices to Hartree-Fock (HF) values, reducing dependence on the specific reference state choice.
  • ADC(3)-D Truncation: The authors implement the ADC(3)-D scheme, which incorporates Coupled-Cluster (CC) amplitudes (specifically from CCD calculations) into the ADC coupling vertices. This effectively includes infinite classes of higher-order diagrams (beyond standard ADC(3)) without increasing the computational scaling, improving accuracy in strongly correlated regimes.

• Implementation Details:

  • Calculations are performed for Symmetric Nuclear Matter (SNM) and Pure Neutron Matter (PNM).
  • Interactions are derived from Chiral Effective Field Theory ($\chi$EFT) at Next-to-Next-to-Leading Order (NNLO), including 2-body and 3-body forces (via the Normal-Ordered 2-Body approximation).
  • Boundary conditions include both Periodic Boundary Conditions (PBC) and Special-Point Twisted-Angle Boundary Conditions (sp-TABC) to minimize finite-size effects.

3. Key Contributions

1. Novel Formalism: Development of a stable, high-order Gorkov-SCGF method that successfully integrates pairing at the static level with third-order dynamical correlations via a hybrid Dyson-Gorkov approach.
2. OpRS Scheme: Introduction and validation of the Optimized Reference State scheme, which allows the use of particle-number-conserving ADC(3) machinery within a superfluid context.
3. ADC(3)-D for Nuclear Matter: First application of the ADC(3)-D truncation (combining ADC and CC) to infinite nuclear matter, demonstrating its superiority over standard ADC(3) and ADC(2) in capturing correlation energy.
4. Comprehensive Validation: Extensive testing of model-space convergence, boundary condition effects (PBC vs. TABC), and sensitivity to different $\chi$EFT Hamiltonians (NNLOsat, $\Delta$NNLO$_{go}$).

4. Results

• Equation of State (EOS):

  • The method provides state-of-the-art predictions for the EOS of SNM and PNM.
  • ADC(3)-D yields the most bound results, correcting earlier estimates and showing excellent agreement with Coupled Cluster (CC) and Many-Body Perturbation Theory (MBPT) benchmarks.
  • For SNM, the results with NNLOsat(450) and $\Delta$NNLO$_{go}$(394) are compatible with the empirical saturation point ($\rho_0 \approx 0.16$ fm$^{-3}$, $E_0 \approx -16$ MeV).
  • The symmetry energy $S(\rho)$ and its slope $L$ are calculated, with NNLOsat(450) underestimating the density dependence compared to the GO models.

• Single-Particle Properties:

  • Momentum Distributions ($\rho(k)$): The method predicts a depletion of occupation numbers below the Fermi surface ($\rho(k_F^-) < 1$) and a non-zero tail above $k_F$. SNM shows stronger depletion than PNM due to stronger correlations.
  • Spectral Functions: The spectral functions reveal the fragmentation of single-particle strength. Near the Fermi surface, well-defined quasiparticle peaks emerge (narrow resonances), while states far from the Fermi surface show significant incoherent background (satellite peaks).
  • Finite-Size Effects: Comparisons between PBC and sp-TABC show that while PBCs are sufficient for the EOS, TABCs provide a much finer momentum mesh, offering a more accurate description of momentum distributions and spectral functions.

• Stability: The hybrid approach remains numerically stable across a wide range of densities, including the low-density regime where superfluidity sets in, avoiding the instabilities that plague standard Dyson-SCGF or CC methods in pure neutron matter.

5. Significance

• Astrophysical Relevance: The robust predictions for the EOS and pairing gaps in neutron matter are crucial for modeling neutron star structure, cooling mechanisms, and glitches. The ability to handle low-density superfluidity ab initio is a major step forward.
• Bridging Theories: The work successfully bridges the gap between Gorkov theory (pairing) and high-order Dyson/CC methods (dynamical correlations), offering a computationally efficient alternative to full Gorkov-ADC(3) or Bogoliubov CC.
• Foundation for EDFs: The results provide a microscopic, ab initio foundation for Nuclear Energy Density Functionals (EDFs). By grounding the static and dynamic response of nuclear matter in fundamental interactions, the authors pave the way for constraining the isovector and pairing sectors of EDFs, improving their predictive power for exotic nuclei and astrophysical scenarios.
• Methodological Benchmark: The paper establishes ADC-SCGF as a reliable tool for infinite matter, validating the hybrid truncation and the ADC(3)-D extension against other state-of-the-art many-body methods.