$\Lambda NN$ Input to Neutron Stars from Hypernuclear Data

Basic Information

• Paper ID: 2411.11751
• Title: $\Lambda NN$ input to neutron stars from hypernuclear data
• Authors: E. Friedman, A. Gal (Hebrew University, Jerusalem)
• Classification: nucl-th nucl-ex
• Conference: Proceedings of International Conference on Exotic Atoms and Related Topics and Conference on Low Energy Antiprotons (EXA-LEAP2024), 26-30 August 2024, Vienna
• Paper Link: https://arxiv.org/abs/2411.11751

Abstract

This study extends the authors' two important works from 2023 by establishing a well-defined $\Lambda$-nucleus optical potential through fitting single-particle binding energy data of $\Lambda$ particles in hypernuclei. The potential comprises a Pauli-corrected linear density term ($\Lambda N$) and quadratic density term ($\Lambda NN$). This work extends the analysis to 21 $\Lambda$ single-particle data points, including 1$d_\Lambda$ and 1$f_\Lambda$ states in medium-heavy and heavy nuclei. The potential depths obtained at nuclear matter density $\rho_0=0.17$ fm$^{-3}$ are: $D^{(2)}_\Lambda=-37.5\mp 0.7$ MeV and $D^{(3)}_\Lambda=+9.8\pm 1.2$ MeV, with total potential depth $D_\Lambda=-27.7\pm 0.5$ MeV. These results provide critical input for resolving the "hyperon puzzle" in neutron star cores.

Research Background and Motivation

Scientific Problem

The core issue addressed by this research is the "hyperon puzzle": observed neutron star masses exceed two solar masses, yet theoretical predictions suggest that $\Lambda$ hyperons in neutron star cores would soften the equation of state, limiting neutron star masses to approximately 1.5 solar masses.

Research Significance

1. Astrophysical Importance: Understanding neutron star internal structure and equation of state is crucial for modern astrophysics
2. Fundamental Physics Value: The $\Lambda NN$ three-body interaction is an essential component of strong interaction theory
3. Experimental Verification: Provides theoretical predictions and testing standards for future JLab experiments

Limitations of Existing Methods

1. Skyrme-Hartree-Fock (SHF) Method: Cannot produce reliable values for $D^{(2)}_\Lambda$ and $D^{(3)}_\Lambda$
2. High-Order Density Terms: Requires consideration of higher-order powers of nuclear density, but their necessity remains unproven
3. Pauli Correlation Effects: Most SHF studies lack explicit treatment of Pauli correlations

Core Contributions

1. Extended Dataset: Expansion from 14 to 21 $\Lambda$ single-particle binding energy data points
2. Precise Potential Parameters: Determination of two-body and three-body potential depths at nuclear matter density
3. Isospin Dependence: Discovery and quantification of isospin suppression effects in $\Lambda NN$ interactions
4. Astrophysical Application: Provision of critical repulsive $\Lambda NN$ interaction strength for resolving the neutron star "hyperon puzzle"
5. Experimental Predictions: Testable predictions for the forthcoming JLab E12-15-008 experiment

Methodology Details

Optical Potential Model

A density-dependent $\Lambda$-nucleus optical potential is employed: $V^{\text{OPT}}_\Lambda(\rho) = V^{(2)}_\Lambda(\rho) + V^{(3)}_\Lambda(\rho)$

Two-Body $\Lambda N$ Interaction Term

$V^{(2)}_\Lambda(\rho) = -\frac{4\pi}{2\mu_\Lambda} \frac{f^{(2)}_A b_0}{1 + \frac{3k_F}{2\pi} f^{(2)}_A b_0} \rho$

This includes Pauli correlation effects, manifested through the Fermi momentum $k_F = (3\pi^2\rho/2)^{1/3}$.

Three-Body $\Lambda NN$ Interaction Term

$V^{(3)}_\Lambda(\rho) = +\frac{4\pi}{2\mu_\Lambda} f^{(3)}_A B_0 \frac{\rho^2}{\rho_0}$

Technical Innovations

1. Pauli Correlation Treatment: Explicit inclusion of Fermi momentum dependence, which is neglected in most SHF studies
2. Isospin Suppression Factor: $F = \frac{(2Z)^2 + (N-Z)^2}{A^2}$ Used to suppress $\Lambda NN$ interactions involving excess neutrons and symmetric nuclear cores
3. Density Decomposition: Separation of nuclear density into symmetric core density and excess neutron density

Experimental Setup

Dataset

• 21 $\Lambda$ single-particle binding energy data points
• Including 1$s_\Lambda$, 1$p_\Lambda$, 1$d_\Lambda$, 1$f_\Lambda$ states
• Coverage from $^{16}_\Lambda$N to $^{208}_\Lambda$Pb mass range
• Exclusion of $^{12}_\Lambda$B and $^{13}_\Lambda$C (poor fit for light nuclei)

Nuclear Density Treatment

1. Light Nuclei: Harmonic oscillator-type density distribution
2. Heavy Nuclei: Two-parameter or three-parameter Fermi distribution
3. Proton Density: Based on experimental charge distribution data
4. Neutron Density: Accounting for neutron skin effects

Fitting Method

• Least-squares $\chi^2$ minimization
• 100% correlation between parameters $b_0$ and $B_0$
• Inclusion of statistical and systematic error analysis

Experimental Results

Main Fitting Results

21 Data Point Fit:

• $b_0 = 1.32 \pm 0.072$ fm (attractive)
• $B_0 = 0.162 \pm 0.020$ fm (repulsive)
• $\chi^2$/degrees of freedom = 0.95

Potential Depths (at $\rho_0 = 0.17$ fm$^{-3}$):

• $D^{(2)}_\Lambda = -37.5 \mp 0.7$ MeV
• $D^{(3)}_\Lambda = +9.8 \pm 1.2$ MeV
• $D_\Lambda = -27.7 \pm 0.5$ MeV

Model Comparison Verification

Extrapolation predictions for $^{16}_\Lambda$N show:

• Model X (complete $\rho^2$ term): Insufficient binding in heavy nuclei
• Model Y (suppression factor F): Good agreement with experiment
• Suppression factor reaches F=0.67 for Pb

JLab Experiment Predictions

For $^{40,48}$Ca$(e,e'K^+)^{40,48}_\Lambda$K reactions:

• 1$s_\Lambda$ state: Approximately 2.5 MeV difference
• 1$p_\Lambda$ state: Exceeding 2 MeV difference
• Within JLab E12-15-008 experimental precision

Theoretical Model Comparison

All models' $D^{(2)}_\Lambda$ values are over-bound relative to the empirical value of approximately -30 MeV, suggesting the need for repulsive $\Lambda NN$ contributions.

Physical Significance and Discussion

Hyperon Puzzle Resolution Mechanism

1. Repulsive $\Lambda NN$ Interaction: The repulsive $D^{(3)}_\Lambda \approx +10$ MeV is sufficient to suppress hyperon appearance in neutron star cores
2. Stiffened Equation of State: Prevents hyperon softening effects, allowing neutron star masses exceeding two solar masses

Physical Origin of Isospin Dependence

• Originates from $\Sigma$ and $\Sigma^*$ intermediate states in single-pion exchange diagrams of $\Lambda NN$
• Isospin dependence $\vec{\tau}_1 \cdot \vec{\tau}_2$ leads to different couplings for symmetric cores and excess neutrons

Consistency with Scattering Data

The over-binding of $D^{(2)}_\Lambda$ in this work is quantitatively consistent with recent $\Lambda p$ scattering data and joint analysis of correlation functions PLB 850 (2024) 138550.

Conclusions and Prospects

Main Conclusions

1. Precise Potential Parameters: First precise determination of $\Lambda N$ and $\Lambda NN$ potential depths from hypernuclear data
2. Repulsive Three-Body Force: $D^{(3)}_\Lambda = +9.8 \pm 1.2$ MeV provides critical input for resolving the hyperon puzzle
3. Isospin Effects: Quantification of isospin suppression in $\Lambda NN$ interactions, testable in JLab experiments

Future Directions

1. Microscopic Model Testing: Provides testing standards for microscopic models in astrophysical calculations
2. JLab Experimental Verification: E12-15-008 experiment will directly test isospin dependence predictions
3. Neutron Star Applications: Application of results to neutron star equation of state calculations

In-Depth Evaluation

Strengths

1. Methodological Innovation: Explicit treatment of Pauli correlation effects, addressing SHF method deficiencies
2. Comprehensive Data Coverage: From light to heavy nuclei, including multiple angular momentum states
3. Clear Physical Picture: Isospin suppression mechanism has explicit theoretical foundation
4. High Practical Value: Directly serves the frontier problem of neutron star physics

Limitations

1. Light Nucleus Deviation: $^{12}_\Lambda$B and $^{13}_\Lambda$C cannot be well fitted
2. Model Simplification: Neglects non-local terms and other secondary effects
3. Systematic Errors: Nuclear density uncertainties introduce approximately 1 MeV systematic error

Scientific Impact

1. Theoretical Contribution: Provides standard optical potential model for $\Lambda$ hypernuclear physics
2. Experimental Guidance: Offers theoretical predictions for large-scale experimental facilities like JLab
3. Astrophysical Application: Directly serves neutron star physics research

This study, through precise phenomenological analysis, provides important experimental constraints for understanding the many-body nature of strong interactions and resolving the neutron star "hyperon puzzle," representing a significant contribution to the interdisciplinary field of nuclear and astrophysics.