Constraining the $ΛΛ$ interaction with terrestrial and astronomical data

This study utilizes a Skyrme energy density functional framework constrained by terrestrial double-$\Lambda$ hypernuclear data and pseudodata to demonstrate that including repulsive $p$-wave and $N\Lambda\Lambda$ three-body forces yields equations of state consistent with neutron star observations, thereby highlighting the necessity of future experimental data on heavier hypernuclei.

The Gist

Imagine the universe as a giant, cosmic kitchen. Most of the ingredients we know—like the protons and neutrons that make up the atoms in your body—are like flour and sugar. They are stable and easy to work with. But in the extreme, high-pressure environments of exploding stars or the dense cores of neutron stars, the kitchen gets so hot and crowded that weird, exotic ingredients start to appear. One of these exotic ingredients is the Lambda particle (or $\Lambda$).

This paper is a recipe book for understanding how these exotic Lambda particles behave when they hang out together in the densest matter in the universe.

Here is the story of the research, broken down into simple concepts:

1. The Mystery of the "Double-Lambda"

In normal matter, protons and neutrons stick together to form nuclei. Sometimes, a Lambda particle sneaks in. But the real mystery happens when two Lambda particles try to hang out in the same nucleus. This is called a "double-Lambda hypernucleus."

Scientists want to know: How do these two Lambdas interact? Do they like each other? Do they push each other away?

• The Problem: We don't have enough data. We've only seen a few tiny, light "double-Lambda" atoms in experiments on Earth. It's like trying to figure out how a whole crowd of people behaves in a stadium by only watching two people in a small room. You can't see the big picture.

2. The "Virtual" Experiment (The Magic Trick)

Since we can't build heavy double-Lambda atoms in a lab yet, the researchers used a clever trick. They built a computer model (a "three-body model") that acts like a virtual laboratory.

• The Analogy: Imagine you want to know how a heavy truck handles a bumpy road, but you only have a toy car. Instead of just guessing, you use a super-accurate simulation to predict how the toy car would behave if it were a truck.
• The team used this simulation to generate "pseudodata" (fake but scientifically sound data) for heavier atoms. They combined the real, tiny data from Earth with this virtual data for bigger atoms.

3. Tuning the "Dial" (The Parameters)

The researchers used a mathematical framework called the Skyrme Energy Density Functional. Think of this as a complex recipe with four main "dials" (parameters) that control how the Lambda particles interact:

• Dial 1 & 2 (The S-Wave): These control how the Lambdas interact when they are close and calm. The team found that using only the tiny, real atoms didn't give them enough information to turn these dials correctly. It was like trying to tune a radio with a broken antenna.
• The Breakthrough: Once they added their "virtual" data for heavier atoms, the dials snapped into place. They could finally figure out exactly how these particles stick together in different-sized systems.

4. The Neutron Star Puzzle (The Heavyweight Champion)

Now, let's zoom out to Neutron Stars. These are the dead cores of massive stars, crushed so tightly that a teaspoon of their material weighs a billion tons.

• The Puzzle: For a long time, scientists thought that if Lambda particles appeared inside a neutron star, they would make the star "squishy" (soft). If the star gets too squishy, gravity wins, and the star collapses into a black hole.
• The Conflict: But, we have observed neutron stars that are twice as heavy as our Sun. If they were too squishy, they would have collapsed long ago. This is known as the "Hyperon Puzzle." How can these stars be so heavy if the Lambdas inside them make them weak?

5. The Solution: Adding "Repulsion"

The researchers tested their new recipe against the neutron stars. They found that the Lambda particles need a little bit of repulsion (a "push") to keep the star from collapsing.

• The Analogy: Imagine a crowded elevator. If everyone just stands there, the elevator might feel heavy and unstable. But if everyone leans back slightly (creating a little space), the elevator becomes more stable and can hold more weight.
• The study showed that by adding a specific "repulsive" force between the Lambdas (specifically the "p-wave" interaction), the neutron star becomes stiff enough to support its own massive weight.

The Big Takeaway

This paper is a bridge between two worlds:

1. The Micro World: Tiny atoms in a lab (and virtual simulations).
2. The Macro World: Giant, heavy neutron stars in space.

By combining data from both, the team created a new, more accurate "rulebook" for how exotic matter behaves. They proved that:

• You need to look at heavy systems (even if they are just virtual) to understand the rules of light systems.
• The "push" between Lambda particles is the secret sauce that allows massive neutron stars to exist without turning into black holes.

In short: They used a mix of real experiments and smart computer simulations to solve a cosmic mystery, showing us exactly how the universe's densest objects stay together. And they are now calling for future experiments to build those heavy double-Lambda atoms in real life to double-check their work!

Technical Summary

1. Problem Statement

The interaction between two $\Lambda$ hyperons ($\Lambda\Lambda$) is poorly constrained compared to nucleon-nucleon ($NN$) or nucleon-hyperon ($N\Lambda$) interactions. This uncertainty leads to significant ambiguities in the Equation of State (EoS) for dense matter containing hyperons, directly impacting our understanding of neutron star (NS) structure.

• The Hyperon Puzzle: The appearance of hyperons in neutron star cores typically softens the EoS, reducing the maximum mass of neutron stars below the observed $2M_\odot$ limit. Resolving this requires precise knowledge of repulsive $\Lambda\Lambda$ and three-body $N\Lambda\Lambda$ forces.
• Data Limitations: Existing experimental data on double-$\Lambda$ hypernuclei are limited to very light systems (e.g., $^9_\Lambda\text{Li}$, $^9_\Lambda\text{Be}$, $^{11}_\Lambda\text{Be}$). Previous studies found that these light systems are insufficient to simultaneously constrain the bulk (volume) and surface parameters of the $\Lambda\Lambda$ interaction, leading to large uncertainties in the EoS.

2. Methodology

The authors employ a Skyrme-type Energy Density Functional (EDF) framework based on the KIDS (Korea-IBS-Daegu-SKKU) models.

• Theoretical Framework:
  • The total energy density is decomposed into $E_{NN}$ (nucleonic), $E_{N\Lambda}$ (nucleon-$\Lambda$), and $E_{\Lambda\Lambda}$ ($\Lambda$-$\Lambda$) contributions.
  • The $\Lambda\Lambda$ interaction is modeled with a Skyrme force containing five parameters:
    • $\lambda_0, \lambda_1$: $s$-wave interaction terms.
    • $\lambda_2$: $p$-wave interaction term.
    • $\lambda_3, \alpha$: Density-dependent term representing an effective repulsive $N\Lambda\Lambda$ three-body force (with $\alpha$ fixed to 1).
• Constraint Strategy (Two-Step Approach):
  1. Terrestrial Constraints (Hypernuclei):
     • The authors fit the $s$-wave parameters ($\lambda_0, \lambda_1$) using experimental binding energies of light double-$\Lambda$ hypernuclei.
     • Crucial Innovation: To break the degeneracy between $\lambda_0$ and $\lambda_1$, they supplement experimental data with pseudodata for heavier systems ($^{18}_\Lambda\text{O}$, $^{34}_\Lambda\text{S}$, $^{42}_\Lambda\text{Ca}$, $^{58}_\Lambda\text{Ni}$). These pseudodata are generated using a core + $2\Lambda$ three-body model (solved via the Gaussian Expansion Method) calibrated to single-$\Lambda$ data and the $^6_\Lambda\text{He}$ double-$\Lambda$ hypernucleus.
     • The $p$-wave ($\lambda_2$) and density-dependent ($\lambda_3$) terms are left unconstrained in this step because they have negligible effects on the ground states of light double-$\Lambda$ systems (where $\Lambda$s occupy $1s$ orbitals).
  2. Astronomical Constraints (Neutron Stars):
     • The constrained $s$-wave parameters are used as a baseline.
     • The remaining parameters ($\lambda_2$ and $\lambda_3$) are explored by calculating Neutron Star properties using the Tolman-Oppenheimer-Volkov (TOV) equations.
     • Observational Benchmarks: The models are tested against:
       • Maximum mass constraint: $M_{max} \geq 2M_\odot$.
       • Mass-Radius (M-R) relation: Consistency with NICER data for PSR J0740+6620.
       • Causality: The speed of sound must not exceed the speed of light.

3. Key Contributions

• Resolution of Parameter Degeneracy: The study demonstrates that light hypernuclear data alone cannot uniquely determine the $s$-wave parameters ($\lambda_0, \lambda_1$). By incorporating pseudodata for medium-mass hypernuclei, the authors successfully constrain these parameters to a narrow, well-defined region.
• Complementary Constraints: The paper establishes a clear division of labor between terrestrial and astronomical data:
  • Hypernuclear data (augmented by three-body models) tightly constrain the $s$-wave interaction ($\lambda_0, \lambda_1$).
  • Neutron star observables are required to constrain the $p$-wave ($\lambda_2$) and density-dependent ($\lambda_3$) terms, which are insensitive to finite nuclei ground states but critical for dense matter.
• Phenomenological EoS: The authors provide a set of phenomenologically acceptable EoSs that satisfy both hypernuclear binding energies and massive neutron star observations.

4. Key Results

• $s$-wave Parameters: The inclusion of pseudodata for $A \geq 18$ reduced the uncertainty in $\lambda_0$ and $\lambda_1$. The resulting parameter sets (labeled LL2) reproduce experimental and pseudodata binding energies with a mean absolute relative deviation (MD) of $< 3\%$.
• Role of Repulsion in Neutron Stars:
  • Without repulsive $\Lambda\Lambda$ terms ($\lambda_2=0, \lambda_3=0$), the EoS becomes too soft, and the maximum neutron star mass drops below $2M_\odot$ (the hyperon puzzle).
  • $p$-wave repulsion ($\lambda_2$): This is the dominant factor in stiffening the EoS. The study finds that $\lambda_2 \gtrsim 300 \text{ MeV fm}^5$ (for KIDS-A) or $\gtrsim 100 \text{ MeV fm}^5$ (for KIDS-D) is required to support $2M_\odot$ stars.
  • Density-dependent term ($\lambda_3$): While less critical than $\lambda_2$, a positive $\lambda_3$ (representing $N\Lambda\Lambda$ repulsion) helps shift M-R curves upward, slightly lowering the threshold for $\lambda_2$.
• Allowed Parameter Ranges: The study identifies the following phenomenologically acceptable ranges:
  • $\lambda_0 \in [-500, -200] \text{ MeV fm}^3$
  • $\lambda_1 \in [250, 500] \text{ MeV fm}^5$
  • $\lambda_2 \in [100, 600] \text{ MeV fm}^5$
  • $\lambda_3 \in [0, 1000] \text{ MeV fm}^6$
• Hyperon Fraction: In massive neutron stars within the acceptable parameter region, the central fraction of $\Lambda$ hyperons ($Y_{\Lambda,c}$) is found to be approximately 0.2–0.3, significantly lower than the $\sim 0.6–0.7$ found in models without repulsive interactions.

5. Significance

• Solving the Hyperon Puzzle: The work provides a concrete mechanism (repulsive $p$-wave and three-body forces) to resolve the hyperon puzzle within a Skyrme EDF framework, allowing for the existence of $2M_\odot$ neutron stars with hyperon cores.
• Future Experimental Guidance: The results highlight the critical need for experimental data on heavier double-$\Lambda$ hypernuclei. The study shows that without such data, the $s$-wave parameters remain unconstrained, limiting the predictive power of the EoS.
• Methodological Benchmark: The paper validates the use of three-body model pseudodata to bridge the gap between light experimental data and the requirements of dense matter physics. It offers a robust starting point for future studies involving multi-$\Lambda$ hypernuclei and relativistic mean-field models with $SU(3)$ flavor symmetry.

In conclusion, this paper successfully integrates terrestrial nuclear physics with astrophysical observations to place rigorous constraints on the $\Lambda\Lambda$ interaction, offering a coherent description of matter from finite hypernuclei to the cores of massive neutron stars.