Parameter adjustment of nuclear leading-order local pairing energy density functionals

This study benchmarks a protocol for adjusting parameters of a local leading-order T=1 pairing energy density functional by matching infinite nuclear matter pairing gaps at the chemical potential, demonstrating that this approach yields consistent predictions for nuclear masses and moments of inertia while highlighting critical pitfalls such as spurious di-nucleon condensation and the significant influence of spin-gradient terms and mean-field contributions on nuclear observables.

The Gist

The Big Picture: Tuning the "Glue" of the Atom

Imagine the atomic nucleus as a crowded dance floor filled with two types of dancers: protons and neutrons. These dancers are constantly moving, but they also have a special habit: they love to pair up. When they pair up, they move in perfect sync, which makes the whole dance floor more stable. In physics, we call this pairing.

Scientists use complex mathematical recipes called Energy Density Functionals (EDFs) to predict how these nuclei behave. Think of these recipes like a cookbook.

• The Main Dish (Particle-Hole Interaction): This describes how the dancers move individually and how they push or pull on each other to form the shape of the nucleus.
• The Secret Sauce (Pairing Interaction): This describes the "glue" that makes the dancers pair up.

The Problem:
For a long time, scientists have had a problem. The "Main Dish" recipes (which describe the individual dancers) are never perfect. They get the general shape right, but they often miss the tiny details of where each dancer stands. Because the "Secret Sauce" (pairing) depends heavily on exactly where the dancers are standing, if you use a slightly different Main Dish recipe, you have to completely rewrite the Secret Sauce recipe to get the right results. It's like trying to bake a cake: if you change the flour brand, you can't just use the same amount of sugar; you have to adjust the whole recipe.

This paper asks: Is there a way to tune the "Secret Sauce" so it works perfectly, no matter which "Main Dish" recipe we use?

The Solution: The "Infinite Dance Floor"

The authors realized that trying to tune the sauce by looking at specific, finite nuclei (like a specific dance floor with 50 dancers) is messy. The details of the dance floor (the shape, the specific dancers) get in the way.

Instead, they decided to look at a theoretical, infinite dance floor (Infinite Nuclear Matter).

• The Analogy: Imagine an endless, perfectly flat dance floor with no walls and no specific corners. The dancers are spread out evenly.
• Why this helps: On this infinite floor, there are no messy edges or weird shapes. The physics becomes smooth and predictable. It's like testing a car engine on a perfect, endless highway rather than on a bumpy, winding mountain road.

The authors built a new computer program to simulate this infinite dance floor. They used it as a calibration tool.

The Experiment: Tuning the Recipe

They took several different "Main Dish" recipes (Skyrme parameter sets) that were known to be good but had different properties (like different "effective masses," which is a fancy way of saying how heavy the dancers feel when they move).

1. The Reference: They picked one specific, trusted "Secret Sauce" recipe (called ULB) that worked well with one specific Main Dish.
2. The Calibration: They used their infinite dance floor simulator to say, "Okay, we want the pairing glue to feel exactly the same on this infinite floor, no matter which Main Dish recipe we use."
3. The Adjustment: They tweaked the numbers in the Secret Sauce recipe (specifically how the glue changes with density) until the infinite dance floor looked identical for all the different Main Dish recipes.

The Results: Does it Work?

They tested these new, adjusted recipes on real, finite nuclei (the actual dance floors with walls).

• Success: They found that by tuning the sauce on the "infinite floor," they could get consistent, accurate results for real nuclei, regardless of which Main Dish recipe they started with.
• The "Odd-Even" Staggering: One specific thing they checked was the "Odd-Even Staggering." Imagine a line of dancers. If you have an even number, they pair up perfectly. If you have an odd number, one dancer is left out. This makes the odd-numbered nuclei slightly less stable. The new method predicted this stability difference very accurately across the board.
• Rotational Inertia: They also checked how the nucleus spins. The new method predicted how fast they spin and how much energy it takes to spin them up, matching real-world data very well.

The Warning: Don't Copy the Wrong Chef

The paper also tested a different approach. Some scientists try to copy the "Secret Sauce" from a different type of recipe entirely (the Gogny force, which is like a different cuisine).

• The Mistake: They tried to tune their local recipe to match the gaps (the pairing strength) found in the Gogny recipe on the infinite floor.
• The Result: It failed. When they applied this to real nuclei, the results were way off.
• The Lesson: Just because two recipes look the same on a smooth, infinite highway doesn't mean they will drive the same car on a bumpy mountain road. The Gogny recipe has "hidden features" (momentum dependence) that a simple local recipe can't mimic just by matching the average gap. You can't just copy-paste the sauce; you have to understand the underlying physics.

The "Bose-Einstein Condensate" Bug

The authors also found a weird glitch. If you tune the sauce parameters too aggressively (specifically, if you make the density dependence too sharp), the simulation predicts that at very low densities, the neutrons stop acting like individual dancers and suddenly collapse into a giant, clumpy blob called a Bose-Einstein Condensate (a state of matter where particles act as a single wave).

• The Reality Check: In real life, neutrons don't do this in the way the math predicted. This is a "spurious instability"—a bug in the math. The paper warns that if you use certain parameter settings, your computer might predict that a nucleus dissolves into a cloud of di-neutron blobs, which is physically impossible.

Summary

What did they do?
They created a new, robust way to tune the "pairing glue" in nuclear physics models. Instead of guessing based on messy real-world data, they used a clean, theoretical "infinite universe" to calibrate the glue.

Why does it matter?
It allows scientists to mix and match different nuclear models without breaking the pairing physics. It makes predictions more reliable for things like the stability of heavy elements and the behavior of neutron stars.

The Takeaway:
To get the best cake, you don't just look at the frosting; you have to understand how the frosting interacts with the specific cake batter you are using. By testing on a "perfect" infinite cake, they found a way to make the frosting work on any cake batter, as long as you don't use the wrong recipe entirely.

Technical Summary

1. Problem Statement

In nuclear structure physics, Energy Density Functionals (EDFs) are used to model the effective interaction between nucleons. While the particle-hole (ph) channel (mean field) is often described by Skyrme functionals, the particle-particle (pp) channel (pairing) is frequently treated with separate phenomenological functionals.

• The Consistency Issue: The strength of pairing correlations depends heavily on the single-particle level density near the chemical potential, which is determined by the ph EDF (specifically the effective mass). Consequently, pairing parameters adjusted for one Skyrme parameter set (e.g., SLy4) often fail to produce consistent results when applied to other Skyrme sets with different effective masses or functional forms.
• The Limitation of Finite Nuclei Fitting: Adjusting pairing parameters directly to observables in finite nuclei (like odd-even mass staggering) is problematic because the single-particle spectra predicted by EDFs often contain systematic errors (shell structure deficiencies). Parameters adjusted to these flawed spectra may "absorb" the errors, leading to poor predictive power when extrapolating to other nuclei or phenomena.
• The Gap in Microscopic Mapping: Attempts to map finite-range interactions (like Gogny) or microscopic calculations onto local Leading-Order (LO) pairing EDFs by matching pairing gaps in infinite nuclear matter (INM) have yielded mixed results, sometimes leading to unrealistic predictions in finite nuclei.

2. Methodology

The authors propose a protocol to adjust the parameters of a local, density-dependent, leading-order (LO) $T=1$ (like-particle) pairing EDF based on the properties of homogeneous infinite nuclear matter (INM) rather than finite nuclei.

• The Pairing EDF: The study utilizes the standard LO pairing functional:
  $$E_{pair} = \int d^3r \frac{V_0}{4} \left[ 1 - \eta \left( \frac{D_0(r)}{\rho_{ref}} \right)^\sigma \right] \sum_{q=n,p} \left( |\tilde{D}_q(r)|^2 + |\tilde{C}_q(r)|^2 \right)$$
  Parameters to be adjusted:

  • $V_0$: Overall strength.
  • $\eta$: Controls the volume ($\eta=0$) vs. surface ($\eta=1$) character.
  • $\sigma$: Controls the width of the density dependence peak.

• Reference Framework:

  • INM Solver: The authors implemented a self-consistent HFB solver for homogeneous, non-polarized INM at arbitrary isospin asymmetry. This solver handles both NLO (Next-to-Leading Order) and N2LO (Next-to-Next-to-Leading Order) Skyrme functionals.
  • Reference Data: Two reference pairing gaps in symmetric INM were used:
    1. ULB Reference: Gaps calculated with the standard SLy4 ph EDF + the established "ULB" pairing parameters.
    2. Gogny D1S Reference: Gaps calculated with the finite-range Gogny D1S force.
  • Adjustment Protocol: The parameters ($V_0, \eta, \sigma$) for various Skyrme ph EDFs (SLy4, 1T2T series with varying effective masses, and SN2LO1) were adjusted to reproduce the density dependence of the pairing gap $\Delta(k_{\lambda})$ at the chemical potential in symmetric INM, matching the reference curves.

• Validation: The resulting parameter sets were tested against three key observables in finite nuclei (Sn, Pb, and Yb isotopic chains):

  1. Odd-even staggering of binding energies ($\Delta^{(5)}_n$).
  2. Odd-even staggering of charge radii ($\Delta^{(3)}_r$).
  3. Rotational moments of inertia (via cranked HFB calculations).

3. Key Contributions and Findings

A. Successful Transfer via INM Gaps

The study demonstrates that adjusting LO pairing parameters to reproduce the pairing gap in symmetric INM allows for the transfer of pairing properties across different Skyrme ph EDFs with varying effective masses.

• When parameters are readjusted to match the INM gaps of a reference (e.g., SLy4+ULB), the resulting combinations (e.g., 1T2T(0.80)+LOULB) yield consistent predictions for odd-even mass staggering and rotational moments of inertia, regardless of the underlying effective mass.
• This confirms that the global trend of pairing correlations is governed by the effective mass and level density, which can be captured by the INM gap.

B. Failure of Finite-Range to Local Mapping

The paper highlights a critical limitation: Matching the gap at the chemical potential ($\Delta(k_\lambda)$) is insufficient to replicate finite-range interactions.

• When the authors adjusted a local LO EDF to reproduce the INM gaps of the finite-range Gogny D1S force (LOD1S parameters), the resulting functionals overestimated pairing gaps in finite nuclei significantly.
• Reason: Finite-range interactions have a strong momentum dependence in the pairing matrix elements ($\Delta_k$) that changes sign at high momenta. A local EDF with a cutoff mimics this only partially. The INM gap $\Delta(k_\lambda)$ does not capture the full momentum structure required to reproduce the matrix elements in the complex wave functions of finite nuclei.

C. Spurious Bose-Einstein Condensate (BEC) Instability

The authors identified a non-physical instability in certain regions of the parameter space.

• For pairing EDFs with very small values of the density-dependence power $\sigma$ (e.g., $\sigma < 0.4$), the HFB equations in INM at low densities yield solutions corresponding to a Bose-Einstein condensate (BEC) of di-nucleons.
• This manifests as a spurious transition where the chemical potential drops below the bottom of the single-particle spectrum, leading to unphysical long-range tails in the pair density of weakly-bound nuclei. This instability must be avoided when choosing parameters.

D. Impact of Time-Odd and Spin-Gradient Terms

• Spin-Gradient Terms: The study reveals that the inclusion or exclusion of spin-gradient terms in the ph EDF (present in 1T2T/SN2LO1 but omitted in SLy4) significantly affects the odd-even staggering of masses. These terms introduce an odd-even effect in the binding energy that is independent of pairing correlations.
• Pairing Contribution to Mean Field: Including the contribution of the pairing EDF to the normal mean field (rearrangement terms) is crucial for reproducing the correct magnitude of the odd-even staggering of charge radii. Omitting this contribution leads to systematically smaller or inverted staggering.

4. Results Summary

Observable	Performance of LOULB (INM-adjusted)	Performance of LOD1S (Gogny-mapped)
Odd-Even Mass Staggering	Excellent agreement across different Skyrme sets (SLy4, 1T2T, SN2LO1).	Systematic overestimation of gaps.
Rotational Moments of Inertia	Consistent description of backbending and $E_\gamma$ trends.	Slightly overestimated (due to stronger pairing).
Charge Radii Staggering	Reasonable size if pairing contribution to mean field is included.	Poor description; fails to capture pair density variations.
INM Stability	Stable for $\sigma \approx 1$.	Unstable (BEC transition) for low $\sigma$.

5. Significance and Outlook

• Robust Parameterization Strategy: The paper establishes a robust protocol for generating families of pairing EDFs. By anchoring the adjustment to the global trends of INM gaps rather than specific finite-nuclei data, one can create consistent pairing functionals for any Skyrme ph EDF, facilitating the study of effective mass dependence and higher-order gradient terms.
• Limitations of Local Approximations: The work provides a strong warning against simply mapping finite-range interactions (like Gogny) or microscopic gaps onto local EDFs solely by matching the INM gap at the chemical potential. The momentum dependence of the interaction is lost in this mapping, leading to failures in finite nuclei.
• Future Directions: The authors suggest that to fully replicate finite-range or microscopic results, one might need to introduce gradient terms in the pairing EDF (active only in inhomogeneous systems), though this would require re-adjusting to finite nuclei, partially negating the benefit of the INM protocol.
• Physical Insight: The identification of the BEC instability in the LO pairing EDF serves as a critical diagnostic tool for the validity of parameter sets, particularly for weakly-bound nuclei and neutron-rich matter.

In conclusion, this study provides a systematic framework for decoupling the adjustment of pairing parameters from the specific deficiencies of single-particle spectra in finite nuclei, while simultaneously highlighting the subtle but critical role of momentum dependence and time-odd terms in achieving predictive accuracy.