Exploring Hyperon Skyrme Forces in Multi-$Λ$ Hypernuclei and Neutron Star Matter

This study employs a comprehensive Bayesian analysis within the Skyrme Hartree-Fock framework to constrain $\Lambda\Lambda$ and $\Lambda\Lambda N$ interaction parameters using both hypernuclear data and astrophysical observations, revealing that repulsive components in these interactions are essential for reconciling hyperon-rich equations of state with the existence of $\sim2\,M_{\odot}$ neutron stars.

The Gist

Imagine a neutron star as the universe's ultimate pressure cooker. It is a city-sized ball of matter so dense that a single teaspoon would weigh a billion tons. Inside this cosmic pressure cooker, the rules of physics get weird. Usually, these stars are made of neutrons, but under such extreme pressure, some neutrons might transform into heavier, stranger cousins called hyperons (specifically, the $\Lambda$ hyperon).

For a long time, scientists have had a major headache trying to understand these stars, known as the "Hyperon Puzzle." Here's the problem: When you add hyperons to the mix, they act like a soft pillow in a mattress. They make the star's internal structure "squishy" (softening the equation of state). If the star gets too squishy, it collapses under its own gravity. But we know from telescopes that some neutron stars are incredibly heavy (about twice the mass of our Sun). If hyperons make them squishy, how do they stay so heavy without collapsing?

This paper is like a team of detectives using a massive amount of evidence to solve the mystery of how these hyperons behave.

The Detective Work: Mixing Two Worlds

The researchers used a method called Bayesian analysis, which is like a super-smart guessing game. They combined two very different types of clues:

1. Lab Clues (Nuclear Data): Experiments on Earth where scientists create tiny "hypernuclei" (atoms with a hyperon inside). This tells them how hyperons behave at low densities, like in a calm room.
2. Space Clues (Astrophysical Data): Observations of real neutron stars, including their mass, size, and how they wobble when they crash into each other (gravitational waves). This tells them how hyperons behave in the extreme pressure of a star.

The Toolkit: The "Skyrme" Force

To model this, the team used a mathematical toolkit called Skyrme forces. Think of this as a recipe book for how particles talk to each other. The recipe has five main ingredients (parameters) that control the interaction between hyperons:

• The "Hug" ($\lambda_0$): A local, short-range attraction.
• The "Push" ($\lambda_1, \lambda_2$): Momentum-dependent forces that act like a repulsive push when particles move fast or get crowded.
• The "Crowd Control" ($\lambda_3, \alpha$): Three-body forces that kick in when there are many particles together, acting like a strong repulsive barrier at high densities.

The Big Discovery: The "Spring" Effect

The paper found that the behavior of hyperons isn't just one thing; it changes depending on how crowded the star is. They discovered a crucial switch:

1. At Low Density (The Hug): When the star isn't too dense yet, the hyperons like to stick together. The "Hug" parameter is strong and attractive. This makes the star a bit softer, just like the old puzzle suggested.
2. At High Density (The Spring): As the star gets squeezed tighter and tighter, the "Push" and "Crowd Control" ingredients take over. The interaction flips from a hug to a repulsive spring.

The Analogy: Imagine a crowd of people in a room.

• Low Density: They are friendly and might even hold hands (attraction).
• High Density: As the room gets packed, they start elbowing each other and pushing back hard to make space (repulsion).

This "spring" effect is the key to solving the puzzle. Even though hyperons try to make the star squishy at first, the repulsive force at high densities acts like a stiffening agent. It stops the star from collapsing, allowing it to support the massive weight of 2 Suns.

What the Numbers Say

The researchers didn't just guess; they calculated the exact "recipe" that fits all the data:

• The Two-Body Force: They found that the direct interaction between two hyperons is tightly constrained. It starts attractive but becomes repulsive at high speeds/densities.
• The Three-Body Force: They found that interactions involving three particles (two hyperons and a nucleon) are essential. These forces act like a final safety net, adding extra stiffness to the star's core.
• The Result: By including these repulsive forces, the maximum weight a neutron star can hold increases by up to 22%. With the extra help of three-body forces, the star can gain another 0.1 solar masses, easily explaining how we see stars that are twice as heavy as our Sun.

The Bottom Line

This paper doesn't just say "hyperons exist." It provides a detailed, experimentally grounded map of how they behave. It shows that nature has a clever trick: hyperons start out friendly but turn into a rigid, repulsive force when the pressure gets too high. This repulsion is what allows the universe's densest stars to remain stable giants rather than collapsing into black holes.

The study is a major step forward, bridging the gap between tiny experiments in a lab and the massive, invisible giants floating in space, finally giving us a coherent picture of what happens inside the heart of a neutron star.

Technical Summary

1. Problem Statement

The internal structure of neutron stars (NS) containing hyperons (strange baryons) presents a significant challenge in nuclear astrophysics known as the hyperon puzzle.

• The Conflict: The inclusion of hyperons (specifically $\Lambda$ hyperons) in the equation of state (EOS) of dense matter typically introduces strong attractive interactions that soften the EOS. This softening reduces the maximum mass of neutron stars, often predicting values below the observed $\sim 2 M_\odot$ (solar masses), which contradicts observations of massive pulsars like PSR J0740+6620.
• The Uncertainty: The interactions between hyperons ($\Lambda\Lambda$) and between hyperons and nucleons ($\Lambda\Lambda N$ three-body forces) are poorly constrained. Experimental data is limited to a few double-$\Lambda$ hypernuclei events, leading to large uncertainties in the $\Lambda\Lambda$ potential depth and the nature of three-body repulsion.
• The Goal: To systematically constrain the parameters of $\Lambda\Lambda$ two-body and $\Lambda\Lambda N$ three-body interactions within the Skyrme Hartree-Fock (SHF) framework by combining terrestrial hypernuclear data with multi-messenger astrophysical observations.

2. Methodology

The authors employ a Bayesian inference approach to constrain the interaction parameters, integrating diverse datasets.

• Theoretical Framework:

  • Model: Skyrme Hartree-Fock (SHF) density functional theory.
  • Interactions:
    • $NN$ (Nucleon-Nucleon): Three standard parameter sets (SLy4, SKI3, SGI).
    • $\Lambda N$: Three parameter sets (YMR, HP$\Lambda$2, SLL4).
    • $\Lambda\Lambda$: Five parameter sets (SLL1, SLL2, SLL3, SLL1', SLL3') based on fitting double-$\Lambda$ hypernuclei.
  • Hamiltonian: The energy density includes terms for local/momentum-independent ($\lambda_0$), momentum-dependent ($\lambda_1, \lambda_2$), and density-dependent three-body ($\lambda_3, \alpha$) interactions.

• Bayesian Analysis:

  • Posterior Distribution: Calculated using $P(\theta|D) \propto P(D|\theta)P(\theta)$, where $\theta$ represents the interaction parameters and $D$ is the dataset.
  • Datasets (Likelihoods):
    1. Astrophysical: Mass-radius measurements from NICER (PSR J0030+0451, J0740+6620, J0437-4715) and tidal deformability constraints from the gravitational wave event GW170817.
    2. Nuclear: The potential depth of the $\Lambda\Lambda$ interaction in pure $\Lambda$ matter at $\rho \approx \rho_0/5$, constrained by experimental data from double-$\Lambda$ hypernuclei ($^{10}_{\Lambda\Lambda}\text{Be}$, $^{13}_{\Lambda\Lambda}\text{B}$, $^{6}_{\Lambda\Lambda}\text{He}$).
  • Sampling: Performed using bilby and pymultinest packages.

• Inference Strategy: Three scenarios were tested: (i) Astrophysical data only (+Astro), (ii) Nuclear data only (+Nucl), and (iii) Combined data (+Astro+Nucl).

3. Key Contributions

• First Comprehensive Bayesian Constraint: This is the first study to simultaneously constrain $\Lambda\Lambda$ and $\Lambda\Lambda N$ Skyrme parameters using both hypernuclear experimental data and modern multi-messenger astrophysical observations.
• Decoupling Interaction Roles: The study successfully isolates the roles of two-body vs. three-body forces. It demonstrates that while the two-body $\Lambda\Lambda$ interaction is tightly constrained by nuclear data, the three-body $\Lambda\Lambda N$ parameters are primarily driven by astrophysical requirements to support massive stars.
• Mechanism Identification: The paper elucidates the microscopic mechanism by which repulsive hyperon interactions resolve the hyperon puzzle, specifically through the density-dependent sign change of the $\Lambda\Lambda$ potential and the stiffening effect of three-body forces on effective masses.

4. Key Results

A. Parameter Constraints

• $\lambda_0$ (Local, momentum-independent): Tightly constrained to be attractive (peak at $\approx -722$ MeV fm$^3$). This is driven by nuclear hypernuclear data.
• $\lambda_1, \lambda_2$ (Momentum-dependent): Prefer positive (repulsive) values. These terms provide necessary repulsion at high densities.
• $\lambda_3, \alpha$ (Three-body $\Lambda\Lambda N$): Constrained almost exclusively by astrophysical data. $\lambda_3$ is large and positive (repulsive), while $\alpha$ favors smaller values, indicating a moderate density dependence.
• Bayes Factors: The combination of the SLy4 ($NN$) and SLL4 ($\Lambda N$) interaction sets yielded the highest Bayes factor, indicating the best simultaneous fit to all data.

B. The $\Lambda\Lambda$ Potential Depth ($U_{\Lambda\Lambda}$)

• The posterior distribution reveals a density-dependent crossover:
  • At low densities ($\rho_\Lambda \lesssim 0.27 \rho_0$), the interaction is attractive, consistent with hypernuclear data.
  • At higher densities, the interaction becomes repulsive. This transition is crucial for stiffening the EOS at high densities.

C. Impact on Neutron Star Properties

• Solving the Hyperon Puzzle:
  • Without hyperon interactions, the inclusion of $\Lambda$ hyperons reduces the maximum mass to $\sim 1.5 M_\odot$ (the "puzzle").
  • Including the constrained $\Lambda\Lambda$ interactions increases the maximum mass by up to 22% (to $\sim 1.98 M_\odot$).
  • Adding $\Lambda\Lambda N$ three-body forces further stiffens the EOS, raising the maximum mass to $\sim 2.03 M_\odot$, consistent with observed massive pulsars.
• Particle Fractions: The repulsive components of $\Lambda\Lambda$ and $\Lambda\Lambda N$ interactions suppress the $\Lambda$ hyperon fraction at high densities, preventing the EOS from becoming too soft.
• Effective Masses: The interactions significantly alter the density evolution of baryon effective masses. The inclusion of three-body forces increases the effective mass of $\Lambda$ hyperons at high density, delaying their accumulation and further stiffening the EOS.
• Radius and Tidal Deformability: While the maximum mass is significantly affected, the radius and tidal deformability of $1.4 M_\odot$ stars remain largely unchanged by the inclusion of these hyperon forces, staying consistent with NICER and GW170817 constraints.

5. Significance

This work provides a unified, experimentally grounded description of dense matter containing strangeness.

• It demonstrates that the "hyperon puzzle" can be resolved without invoking exotic phases of matter (like quark matter) or ad-hoc repulsion, provided that $\Lambda\Lambda$ and $\Lambda\Lambda N$ interactions are correctly constrained.
• It establishes a robust link between terrestrial hypernuclear experiments (J-PARC, FAIR) and astrophysical observations (NICER, LIGO/Virgo).
• The Bayesian framework offers a statistical pathway to refine these constraints as new experimental data on double-$\Lambda$ hypernuclei and gravitational waves become available.
• The study highlights the critical necessity of three-body hyperonic forces in nuclear theory to explain the existence of $2 M_\odot$ neutron stars.

Limitations & Future Work: The current study focuses only on $\Lambda$ hyperons. Future extensions must include $\Sigma$ and $\Xi$ hyperons, whose interactions are even less constrained but may further soften the EOS, requiring even stronger repulsive mechanisms.