---
title: Sigma Single-Particle Potential in Nuclear Matter
url: https://www.emergentmind.com/topics/sigma-single-particle-potential
type: topic
---

# Sigma Single-Particle Potential in Nuclear Matter

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                           | [1612.08399]      |
| Chiral EFT NLO13 (old, pre-E40 data)         | +3.7 to +11.2                 | [2602.04281]      |
| Chiral EFT SMS NLO, N$^2$LO, E40-constrained | –10 to –11                    | [2602.04281]      |
| Chiral EFT NLO + 3BF (most-repulsive)        | $+30\pm20$                    | [2508.19560]      |
| Chiral EFT NLO + 3BF, G-matrix (BHF)         | +15                           | [1802.05388]      |

At higher density ($2\rho_0$), the $\Sigma$ potential generally increases in repulsion, reaching values of $+40$ MeV ([1802.05388]) or even $+80$–$+100$ MeV depending on the implementation of 3BFs and the underlying $YN$ interaction ([2508.19560]).

## 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 ([1612.08399]).
- **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 ([2602.04281]).
- **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 ([2508.19560]).
- **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 ([2508.19560], [1802.05388]).

## 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}$ ([2508.19560]). 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$ ([2602.04281], [1612.08399]), 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 ([1612.08399], [2508.19560]).

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 ([1802.05388]).

| Isospin Channel           | $U_\Sigma(\rho_0,0)$ [MeV] | Reference    |
|--------------------------|----------------------------|--------------|
| $\Sigma$ (SNM, Lattice)  | +11                        | [1612.08399] |
| $\Sigma^-$ (PNM, Lattice)| +12                        | [1612.08399] |
| $\Sigma^0$ (PNM, Lattice)| +10                        | [1612.08399] |
| $\Sigma^-$ (PNM, NLO)    | +20...+45                  | [1802.05388] |

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 ([1802.05388], [2508.19560]).
- A weakly attractive or small $U_\Sigma$ would permit earlier $\Sigma$ onset and potentially soften the EOS, reinstating the "hyperon puzzle" ([2602.04281]).
- Heavy-ion collision observables, such as directed flow of $\Sigma^0$, are sensitive to the magnitude of $U_\Sigma$, providing indirect empirical validation ([2508.19560]).

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**

- [1612.08399]: Lattice QCD + BHF extraction of $U_\Sigma$ and comparison with empirical data.
- [2508.19560]: Chiral EFT G-matrix three-body-force constrained $\Sigma$ potential for astrophysical and heavy-ion applications.
- [2602.04281]: Chiral EFT NLO/N$^2$LO analysis including J-PARC E40 constraints; qualitative shift in sign of $U_\Sigma$.
- [1802.05388]: Chiral EFT NLO $G$-matrix results for $U_\Sigma$ with and without $\Lambda N$–$\Sigma N$ coupling and empirical discussion.

Source: https://www.emergentmind.com/topics/sigma-single-particle-potential