Mean-Field Quark Model of Confinement

• Mean-Field Quark Model of Confinement is a relativistic effective-field framework that implements quark confinement via scalar–vector harmonic oscillator potentials and bag boundaries.
• It unifies quark-level corrections with meson-exchange interactions to accurately reproduce baryonic spectra, nuclear saturation properties, and neutron star observables.
• The model enables exploration of hypernuclear structure, dense matter equations of state, and QCD phase transitions through self-consistent mean-field and chiral extensions.

The mean-field quark model of confinement encompasses a class of relativistic effective-field approaches in which the fundamental color-confinement mechanism is implemented at the level of Dirac quarks subjected to one or more explicit confining potentials, most commonly of scalar–vector harmonic oscillator type, or, in hybrid models, an additional bag boundary. In these frameworks, the collective properties of nuclei, hypernuclei, and stellar matter are modeled by embedding the quark-level structure of baryons into a hadronic mean-field theory, in which baryon–baryon interactions arise from meson exchange and/or mean fields. This synthesis enables a microscopic, unified parametric connection between baryonic spectra, saturation properties of nuclear matter, hyperon physics, and the behavior of dense matter in neutron stars.

1. Quark-Level Confinement Mechanisms

In mean-field quark models such as the quark mean-field (QMF) model, confinement is imposed on each constituent quark (e.g., $u$, $d$, $s$, $c$) within a baryon via a phenomenological scalar–vector mixed harmonic oscillator potential,

$U_q(r) = \frac{1}{2} (1+\gamma^0)(a_q r^2 + V_q),$

where $a_q$ and $V_q$ are flavor-dependent potential parameters fit to baryon masses and charge radii (Xing et al., 2016, Wu et al., 2020, Xing et al., 2017). The quark Dirac Hamiltonian, with this confining term, reads

$$H_q = \boldsymbol{\alpha} \cdot \mathbf{p} + \beta m_q + U_q(r).$$

For baryons embedded in nuclear matter, mean fields ($\sigma$, $\omega$, $\rho$) couple to the quarks via density-dependent shifts in mass and chemical potential,

$$m_q^* = m_q - g_{\sigma q} \sigma, \qquad \epsilon_q^* = \epsilon_q - g_{\omega q} \omega - \tau_3 g_{\rho q} \rho.$$

For each baryon, the ground-state energies of the confined quarks are combined and corrected for center-of-mass motion, pion-loop and one-gluon-exchange effects to reproduce physical baryon masses. In the QMF bag (QMFB) variant, an infinite wall at radius $R$ is imposed, enforcing vanishing quark flux at the MIT bag boundary and dominating over the soft confining potential as the primary confinement mechanism (Zhu et al., 2018).

2. Determination and Calibration of Model Parameters

The confining-potential strengths $(a_q, V_q)$ are fixed by demanding that the resultant baryon mass spectrum and root-mean-square radii match empirical data. For $u,d$ quarks, couplings are typically constrained to yield $M_N = 939$ MeV and $r_N \approx 0.87$ fm, while $s$ and $c$ quark parameters are calibrated using the masses of strange and charm baryons ($\Lambda$, $\Sigma$, $\Xi$, $\Lambda^+_c$, $\Sigma^+_c$, $\Xi^{++}_{cc}$) (Wu et al., 2020, Xing et al., 2017, Hu et al., 2013). The energy spectrum for a confined quark is found by solving the radial Dirac equation with proper boundary conditions, leading to quantization constraints that connect potential parameters to the quark masses and energies. For bag-type models, the bag constant $B$ and zero-point correction $Z_0$ are additional parameters, determined to maintain stability of baryons and the saturation properties of nuclear matter (Zhu et al., 2018).

Pion and gluon corrections are evaluated perturbatively. The one-pion self-energy for a baryon is

$$\delta M_B^\pi = -C_B f_{NN\pi}^2 I_\pi, \qquad I_\pi = \frac{1}{\pi m_\pi^2} \int_0^\infty \! dk \, \frac{k^4 u^2(k)}{k^2 + m_\pi^2},$$

where $u(k)$ is the quark wavefunction form factor and $C_B$ an SU(6) spin-flavor coefficient. One-gluon-exchange corrections are computed separately for color-electric and color-magnetic channels using quark current overlaps. Center-of-mass motion is subtracted using explicit formulas (Xing et al., 2016, Xing et al., 2017).

3. Embedding Quark-Level Structure in Hadronic Mean-Field Theory

The effective baryon masses $M_B^*(\sigma)$, derived from their constituent quark structure and media-induced shifts, serve as input to a relativistic hadronic Lagrangian coupling baryons to mesonic mean fields:

$$\mathcal{L} = \sum_B \bar{\psi}_B \left[i\gamma^\mu\partial_\mu - M_B^*(\sigma) - g_{\omega B} \gamma^0 \omega_0 - g_{\rho B}\gamma^0\tau_3 \rho_0 \right] \psi_B - \mathcal{V}_{\text{meson}},$$

where $\mathcal{V}_{\text{meson}}$ encapsulates $\sigma$, $\omega$, and $\rho$ kinetic and self-interaction terms. The baryon-meson couplings $g_{M B}$ are connected to quark-meson couplings $g_{Mq}$ via counting rules ($g_{\omega N} = 3g_{\omega q}$, $g_{\omega \Lambda} = 2g_{\omega q}$, etc.) or further adjusted for SU(3) breaking and hyperon potentials (Hu et al., 2013, Xing et al., 2017).

Mean-field equations for the meson fields and baryon Dirac equations are solved self-consistently, with source terms (scalar, vector, isovector, and charge densities) determined from the baryon ground states. The coupled system ensures that environmental effects (nuclear medium, isospin asymmetry, hyperonic composition) feed back to the quark substructure via the mean fields, producing density-dependent baryon properties and nuclear matter EOS (Xing et al., 2016, Xing et al., 2017).

4. Self-Consistency Loop and Numerical Implementation

Numerical solution proceeds via iterative self-consistency. Starting from initial guesses for the meson fields, quark-level Dirac equations in the baryons yield updated effective baryon masses $M_B^*(\sigma)$. These are then used to populate the baryonic ground states, generating new densities and currents for the meson field equations. The meson solutions are fed back into the quark Dirac Hamiltonians, and the procedure is repeated until convergence in all fields, single-particle energies, and the total binding energy is achieved. The total energy in a static nucleus or infinite matter is

$$E_{\text{tot}} = \sum_B \sum_{n\ell j} n_B \varepsilon_{n\ell j}^{(B)} - \frac{1}{2} \int d^3r \, [g_{\sigma B} \sigma \rho_s^B - g_{\omega B} \omega \rho_v^B - g_{\rho N} \rho_3 \rho_3^v - eA\rho_c^v] + \text{nonlinear meson terms},$$

where the mean-field contributions to the EOS are computed directly from the output fields (Wu et al., 2020).

5. Physical Implications and Applications

The mean-field quark model of confinement underpins a computational framework that can simultaneously and consistently address:

• Nucleon and hyperon spectra, radii, and spin–orbit splittings;
• Systematics of single-$\Lambda$ and single-$\Xi^0$ levels in hypernuclei, with typical agreement to measured $\Lambda$ energies at the $\sim2{-}3\%$ level (Xing et al., 2017, Wu et al., 2020);
• Properties of bulk nuclear matter, including effective nucleon mass $M^*_N(\sigma)$, saturation energy and density, symmetry energy $S(\rho)$, and its slope $L$ (Xing et al., 2016);
• Neutron star structure: equations of state with and without hyperonic degrees of freedom, maximum masses up to $2.1\,M_\odot$ (for QMF-NK1S-3S parameter sets with hyperons), radii $R_{1.4} \sim 13\, \text{km}$, and the softening effects of hyperons and bag confinement on the EOS (Xing et al., 2017, Zhu et al., 2018, Hu et al., 2013);
• Sensitivity of observable properties (hyperon thresholds, neutron star maximum mass) to quark model details such as SU(3) breaking or bag boundary conditions.

The QMF and QMFB models differ in that the latter, incorporating a hard-wall bag boundary, shifts the burden of confinement from the soft harmonic potential to the bag. This sharpens the density dependence of effective baryon masses, increases the scalar response, enhances high-density symmetry energy, and suppresses the usual $R$–$L$ correlation in neutron star radii (Zhu et al., 2018).

6. Extensions and Alternative Mean-Field Confinement Schemes

Beyond the canonical QMF/QMFB implementations, further developments have incorporated additional dynamical confinement mechanisms:

• Chiral mean-field models with Polyakov loop variables realize statistical confinement and deconfinement transitions via an effective background field $\Phi$, reproducing QCD phase diagram features and lattice QCD thermodynamics (Rau et al., 2013).
• Extended Nambu–Jona-Lasinio (NJL) models supplement the traditional four-quark interaction with a six-quark color-singlet channel for baryons and momentum-dependent "confining" couplings $G_3(|\tilde{P}|)$ which diverge in the infrared, effectively preventing free three-quark propagation at low momentum, and dynamically generating nucleons as RPA poles beyond mean field (Cao, 2024).
• Confinement via positivity violation and nonlocal quark mass functions: Models inspired by Gribov–Zwanziger theory yield a quark propagator with complex conjugate poles and a spectral function violating reflection positivity, signifying absence of physical asymptotic single-quark states (Dudal et al., 2013).

These alternative approaches emphasize that mean-field quark confinement can be implemented by a variety of mechanisms—potential well, bag boundary, nonlocal propagators, or collective interaction vertices—all enforcing unobservability of colored states and the dynamical generation of baryons.

7. Limitations and Outlook

The mean-field quark model of confinement, while successful in connecting hadronic spectra to nuclear and astrophysical observables within a unified framework, remains dependent on phenomenological parameterizations of the confining potential and on empirical calibration at the level of baryonic masses and nuclear properties. The absence of first-principles derivation of confinement from QCD and the variable treatment of quark mass function, bag boundary, and medium effects limit strict predictive power. Nevertheless, QMF and its extensions provide testable bridges between QCD-scale dynamics and the macroscopic behavior of strongly interacting matter, enabling systematic exploration of hypernuclear structure, neutron star phenomenology, and the properties of matter at extreme densities (Wu et al., 2020, Hu et al., 2013, Xing et al., 2017, Xing et al., 2016, Zhu et al., 2018, Rau et al., 2013, Cao, 2024, Dudal et al., 2013).