Sigma Single-Particle Potential in Nuclear Matter

• Sigma single-particle potential is defined as the effective potential experienced by a Sigma hyperon in nuclear matter, derived via self-consistent BHF or G-matrix methods.
• It encapsulates the impact of two-body and three-body forces along with channel couplings (e.g., ΛN–ΣN) that influence the depth and sign of the potential.
• Quantitative models show that potential values vary with density and momentum, with increased repulsion at high density playing a key role in addressing the hyperon puzzle in neutron stars.

The Sigma ($\Sigma$) single-particle potential is a central quantity in nuclear and hypernuclear many-body theory, entering both as a non-relativistic mean-field for $\Sigma$ hyperons in nuclear matter and as the analog of the nucleonic mass operator when extended to hyperonic degrees of freedom. In the context of hyperonic and baryonic matter calculations, $U_\Sigma(\rho, k)$ characterizes the effective potential experienced by a $\Sigma$ baryon of momentum $k$ in a nuclear medium of density $\rho$. It encapsulates the effects of two-body and, when included, three-body baryonic forces, as well as channel couplings such as $\Lambda N \leftrightarrow \Sigma N$, and is essential for the understanding of hypernuclear structure, heavy-ion observables, and the equation of state of dense matter.

1. Formal Definition and Theoretical Framework

The $\Sigma$ single-particle (s.p.) potential, $U_\Sigma(k, \rho)$, is obtained self-consistently via the Brueckner–Hartree–Fock (BHF) or $G$-matrix methods in nuclear matter, employing modern hyperon-nucleon ($YN$) interactions derived from chiral effective field theory (EFT) or from lattice QCD approaches.

In the BHF approach, the in-medium $G$-matrix, $G_{YN,Y'N'}(\omega)$, satisfies the Bethe–Goldstone equation: $$G_{YN,Y'N'}(\omega) = V_{YN,Y'N'} + \sum_{Y''N''} V_{YN,Y''N''} \frac{Q_{Y''N''}}{\omega - E_{Y''}(k_{Y''}) - E_{N''}(k_{N''}) + i\eta} G_{Y''N'',Y'N'}(\omega)$$ where $V_{YN,Y'N'}$ is the bare baryon-baryon interaction (e.g., from chiral EFT), $Q_{Y''N''}$ is the Pauli projection operator, $\omega$ is the starting energy, and $E_{B}(k) = m_B + k^2/(2m_B) + U_B(k)$ includes the self-consistent potential for all baryons.

The $\Sigma$ single-particle potential in symmetric nuclear matter (SNM) at density $\rho$ and momentum $k$ is given by: $$U_\Sigma(k, \rho) = \sum_{N \leq k_F} \Re \langle \Sigma k, N k' | G_{\Sigma N, \Sigma N}(\omega) | \Sigma k, N k' \rangle$$ with $k_F$ the nucleon Fermi momentum and the sum restricted to occupied nucleon states.

2. Key Contributions and Channel Couplings

Hyperon-nucleon interactions, particularly the $\Sigma N$ and $\Lambda N$ sectors, involve strong coupling, notably between $\Lambda N$ and $\Sigma N$ via one-pion exchange. In chiral EFT, S-wave and P-wave contact terms, as well as meson-exchange components, are determined from low-energy scattering data, recently constrained by experiments such as the J-PARC E40 $\Sigma^+p$ scattering measurements. Channel couplings modify the effective $\Sigma N$ interaction and, by extension, the depth and sign of $U_\Sigma$.

For realistic nuclear matter, three-body forces (3BF) among nucleons and hyperons are included by normal-ordering the 3BF with respect to the nucleon Fermi sea, resulting in density-dependent effective two-body terms in $V_{YN}$:

• $\Lambda NN-\Lambda NN$ and $\Lambda NN-\Sigma NN$ 3BFs are prominent.
• For the $\Sigma$ potential, 3BF effects appear indirectly via the modification of the $\Lambda$ potential in the $G$-matrix denominator. Direct $\Sigma NN-\Sigma NN$ 3BFs are often omitted.

In finite nuclei, the situation is handled via mass operators and mean fields, as in the Theory of Finite Fermi Systems, but in hyperonic sectors, the self-consistent nuclear matter approach using the $G$-matrix is standard.

3. Quantitative Results and Model Dependence

Tabulated below are key results for $U_\Sigma(0)$ at saturation density ($\rho_0 \simeq 0.16$ fm$^{-3}$), highlighting variations among interaction models and computational approaches:

Approach / Model	$U_\Sigma(\rho_0, k=0)$ [MeV]	Reference
Lattice QCD + BHF (HAL QCD method)	+11	(Inoue et al., 2016)
Chiral EFT NLO13 (old, pre-E40 data)	+3.7 to +11.2	(Jinno et al., 4 Feb 2026)
Chiral EFT SMS NLO, N$^2$LO, E40-constrained	–10 to –11	(Jinno et al., 4 Feb 2026)
Chiral EFT NLO + 3BF (most-repulsive)	$+30\pm20$	(Jinno et al., 27 Aug 2025)
Chiral EFT NLO + 3BF, G-matrix (BHF)	+15	(Kohno, 2018)

At higher density ($2\rho_0$), the $\Sigma$ potential generally increases in repulsion, reaching values of $+40$ MeV (Kohno, 2018) or even $+80$–$+100$ MeV depending on the implementation of 3BFs and the underlying $YN$ interaction (Jinno et al., 27 Aug 2025).

4. Physical Interpretation and Empirical Constraints

The $\Sigma$ single-particle potential determines the possibility of bound $\Sigma$ states in nuclei and plays a crucial role in the appearance of hyperons in neutron star matter. Key observations and physical consequences include:

• Sign and Magnitude: Early models and lattice QCD indicate a moderately repulsive $U_\Sigma(\rho_0,0)=+11$ MeV, in fair qualitative agreement with empirical expectations from $\Sigma^-$-atoms and quasifree production, which suggest $+20$ MeV with $\pm5$–10 MeV uncertainty (Inoue et al., 2016).
• Chiral EFT Evolution: Modern chiral EFT constrained by the J-PARC E40 $\Sigma^+p$ scattering data yields a shift to weakly attractive values, $U_\Sigma(\rho_0,0) \approx -10$ MeV, distinct from the earlier consensus on repulsion (Jinno et al., 4 Feb 2026).
• Momentum and Density Dependence: The potential decreases in magnitude with increasing $k$, remaining weakly attractive or nearly vanishing by $k\approx1.5$ fm$^{-1}$. With increasing density, $U_\Sigma$ grows more repulsive in most models that include hard three-body repulsion (Jinno et al., 27 Aug 2025).
• Empirical Significance: The potential determines the depth of possible $\Sigma$ hypernuclei and regulates the onset of $\Sigma^-$ hyperons in neutron stars. The increasing repulsion at high density provides a mechanism to stiffen the equation of state and alleviate the "hyperon puzzle" in compact star physics (Jinno et al., 27 Aug 2025, Kohno, 2018).

5. Computational Techniques and Parametrizations

For practical applications in simulations spanning from heavy-ion collisions to astrophysical modeling, $U_\Sigma(\rho, k)$ is parametrized:

$$U_\Sigma(\rho, k) = U_\rho\!\left(u = \frac{\rho}{\rho_0}\right) + U_m^0(\rho, k)$$

with density ($U_\rho$) and momentum ($U_m^0$) dependent parts,

$U_\rho(u) = a\,u + b\,u^{4/3} + c\,u^{5/3}$

$$U_m^0(\rho, k) = \frac{C}{\rho_0} \int d^3k' \frac{f(k')}{1 + [(\mathbf{k} - \mathbf{k}')/\mu]^2}$$

where coefficients $\{a, b, c, C, \mu\}$ are tuned to match $G$-matrix results up to $3\rho_0$ and $k\lesssim2.5$ fm$^{-1}$ (Jinno et al., 27 Aug 2025). This facilitates embedding in transport codes (e.g., RQMDv2) for macroscopic observables such as hyperon flow in heavy-ion collisions.

Uncertainty estimates draw on cutoff variations and chiral expansion systematics, with typical errors of $\pm10$–20 MeV at $\rho_0$ (Jinno et al., 4 Feb 2026, Inoue et al., 2016), increasing at higher density.

6. Comparison with $\Lambda$ Single-Particle Potential and Isospin Effects

The $\Sigma$ potential is systematically less attractive (or more repulsive) than the $\Lambda$ potential, which is empirically set at $U_\Lambda(\rho_0,0)\simeq-30$ MeV to match hypernuclear data (Inoue et al., 2016, Jinno et al., 27 Aug 2025).

Isospin effects are non-negligible: in pure neutron matter (PNM), the $\Sigma^-$ potential remains repulsive and insensitive to $\Lambda N$–$\Sigma N$ coupling, while $\Sigma^0$ and $\Sigma^+$ potentials are modified by coupling and renormalization effects (Kohno, 2018).

Isospin Channel	$U_\Sigma(\rho_0,0)$ [MeV]	Reference
$\Sigma$ (SNM, Lattice)	+11	(Inoue et al., 2016)
$\Sigma^-$ (PNM, Lattice)	+12	(Inoue et al., 2016)
$\Sigma^0$ (PNM, Lattice)	+10	(Inoue et al., 2016)
$\Sigma^-$ (PNM, NLO)	+20...+45	(Kohno, 2018)

This structure reflects the interplay between SU(3) irreducible channels and the empirical constraint from $\Sigma^-$–atom data.

7. Implications for Hypernuclear and Astrophysical Systems

The depth and sign of $U_\Sigma$ dictate the existence of bound $\Sigma$ hypernuclei and the threshold for hyperonization in neutron stars:

• A strongly repulsive $U_\Sigma$ suppresses bound $\Sigma$ states and delays the appearance of $\Sigma^-$ in neutron star matter, helping maintain a stiff equation of state and thereby supporting massive $\sim2M_\odot$ neutron stars (Kohno, 2018, Jinno et al., 27 Aug 2025).
• A weakly attractive or small $U_\Sigma$ would permit earlier $\Sigma$ onset and potentially soften the EOS, reinstating the "hyperon puzzle" (Jinno et al., 4 Feb 2026).
• Heavy-ion collision observables, such as directed flow of $\Sigma^0$, are sensitive to the magnitude of $U_\Sigma$, providing indirect empirical validation (Jinno et al., 27 Aug 2025).

The current experimental constraints do not yet fully resolve the width of theoretical predictions for $U_\Sigma$, and ongoing measurements, including those of $\Sigma$ hypernuclei and hyperon-nucleon scattering, remain critical for reducing uncertainties.

---

References

• (Inoue et al., 2016): Lattice QCD + BHF extraction of $U_\Sigma$ and comparison with empirical data.
• (Jinno et al., 27 Aug 2025): Chiral EFT G-matrix three-body-force constrained $\Sigma$ potential for astrophysical and heavy-ion applications.
• (Jinno et al., 4 Feb 2026): Chiral EFT NLO/N$^2$LO analysis including J-PARC E40 constraints; qualitative shift in sign of $U_\Sigma$.
• (Kohno, 2018): Chiral EFT NLO $G$-matrix results for $U_\Sigma$ with and without $\Lambda N$–$\Sigma N$ coupling and empirical discussion.