Net Quark Number Gain in QCD

• Net quark number gain is an observable in finite-temperature, finite-density QCD that quantifies the medium’s change in net quark content when a static quark or antiquark probe is inserted.
• It distinguishes between meson-like screening at low chemical potential and baryon-like screening at high chemical potential, linking confinement and deconfinement phases.
• Defined via derivatives of the Polyakov loop with respect to the quark chemical potential, it serves as a complementary probe alongside generalized susceptibilities for QCD phase-diagram analysis.

Searching arXiv for papers on net quark number gain and related QCD fluctuation observables. Net quark number gain is an observable in finite-temperature, finite-density QCD that measures how the medium’s total net quark number responds to the insertion of a static quark or antiquark probe, excluding the explicit conserved charge carried by the probe itself. In its modern formulation, it is tied to derivatives of the Polyakov loop with respect to the quark chemical potential, and it probes how the thermal bath screens a color source. At high temperature the observable is insignificant, whereas at low temperature it reveals whether the medium prefers meson-like or baryon-like screening configurations; in that sense it connects conserved-charge response, confinement, and the organization of the QCD phase diagram (Surkau et al., 30 Sep 2025).

1. Definition as a probe-induced conserved-charge response

For a static quark probe, the net quark number gain is defined by

$$\Delta Q_q + 1 = 1 + T \frac{\partial}{\partial \mu}\ln \ell,$$

with the Polyakov loop related to the free-energy cost of inserting the probe through

$\ell = e^{-\beta \Delta F_q}.$

An analogous definition applies to an antiquark probe through the anti-Polyakov loop $\bar \ell$. The observable measures the change in the medium’s net quark content, not the explicit quark number of the probe itself. In the detailed formulation summarized for “Using net quark number gain to probe the phases of QCD,” it is described as a gauge-invariant and renormalization-group-invariant observable (Surkau et al., 30 Sep 2025).

This definition is distinct from the usual net quark number density and from generalized susceptibilities. In lattice and effective-theory studies, generalized quark number susceptibilities are defined as

$$\chi_n^q = \frac{1}{V\beta}\frac{\partial^n \ln Z_\mu}{\partial \mu^n},$$

with $\chi_1^q$ giving the net quark number density and higher orders encoding higher responses to the chemical potential (Gattringer et al., 2014). The distinction is structural: the susceptibilities differentiate the partition function of the unprobed medium, whereas net quark number gain is a global difference between the system with and without a static color source.

That distinction is physically consequential. Ordinary net quark number density quantifies how a medium loads conserved charge under a change of $\mu$, while net quark number gain asks what additional quark or antiquark content the bath itself must supply to make the insertion of a single color source thermodynamically consistent.

2. Polyakov loops, confinement, and screening channels

The central physical interpretation of net quark number gain is screening. In the deconfined regime, the Polyakov loops approach unity, and the medium’s net quark content changes negligibly when a static quark or antiquark is inserted. In the confined regime, by contrast, isolated color sources are not supported as asymptotic objects; the bath must provide additional quarks or antiquarks so that the inserted source is embedded into a color-neutral hadron-like configuration (Surkau et al., 30 Sep 2025).

For temperatures much smaller than the constituent quark masses, the low-temperature asymptotics give

$$\Delta Q_q + 1 \simeq \frac{3}{1 + e^{-3\beta\mu} f_{\beta M}/f_{2\beta M}}.$$

In the $T \to 0$ limit, this tends to a step function. For $\mu < M/3$, one has meson-like screening, with

$\Delta Q_q + 1 \to 0,$

so that the medium supplies an antiquark and the total screened configuration has vanishing net quark number. For $\ell = e^{-\beta \Delta F_q}.$0, one has baryon-like screening, with

$\ell = e^{-\beta \Delta F_q}.$1

so that the medium supplies two additional quarks and the total screened configuration carries baryon-like quantum numbers (Surkau et al., 30 Sep 2025).

This step-like behavior is the sharpest manifestation of the observable’s content. It does not merely register the presence of quarks in the bath; it identifies the preferred screening cloud around a static source. It also clarifies how a single quark or antiquark can be added to a confining medium in the first place: the medium furnishes the missing quarks or antiquarks needed to complete a hadron-like state. At finite temperature the transition between the meson-like and baryon-like regimes is smooth rather than discontinuous, but the same screening logic remains operative.

3. Relation to quark-number fluctuations and complex distributions

Although net quark number gain is defined through a static probe, it belongs to a broader family of observables characterizing conserved-charge response in QCD. At finite chemical potential, the quark number operator

$\ell = e^{-\beta \Delta F_q}.$2

is non-hermitian, so its value on individual gauge configurations can be complex. The corresponding distribution

$\ell = e^{-\beta \Delta F_q}.$3

therefore lives in the complex plane, and because the fermion determinant is complex, the measure is not real and positive. In leading-order chiral perturbation theory and for $\ell = e^{-\beta \Delta F_q}.$4, the mean net quark number vanishes while the variance is nonzero; the physical result arises from delicate cancellations between large real and imaginary contributions (Lombardo et al., 2010).

This point is conceptually important for any discussion of “gain.” A vanishing average conserved charge does not imply an absence of large underlying fluctuations. In the small-$\ell = e^{-\beta \Delta F_q}.$5 regime, the distribution is strongly oscillatory and complex-valued in the imaginary part, with amplitude growing exponentially with the volume, and all values of the phase of the fermion determinant are important for the ensemble average (Lombardo et al., 2010).

In effective chiral models, the same theme appears in a different language. Near the chiral crossover, the net quark number probability distribution

$\ell = e^{-\beta \Delta F_q}.$6

is substantially modified by critical dynamics. In the quark-meson model treated with the functional renormalization group, the distribution near the $\ell = e^{-\beta \Delta F_q}.$7 chiral crossover is narrower than the corresponding mean-field and Skellam distributions, and higher cumulants such as the sixth-order one encode the sensitivity of the tails to critical fluctuations (Morita et al., 2013). In the Polyakov loop-extended quark-meson model, higher moments of the net-quark number density exhibit a peculiar structure near the phase transition, and ratios such as $\ell = e^{-\beta \Delta F_q}.$8 act as probes of deconfinement and chiral restoration (Skokov et al., 2010).

The broader implication is that probe-induced charge gain, ordinary density response, and the full probability distribution of net quark number are not interchangeable quantities, but they are tightly coupled aspects of the same thermodynamic sector.

4. Critical behavior and phase-diagram diagnostics

Generalized susceptibilities of the net-baryon number,

$\ell = e^{-\beta \Delta F_q}.$9

are widely used to diagnose QCD criticality. In the small-quark-mass mapping of the three-dimensional Ising model onto the QCD $\bar \ell$0 plane, the sixth-, eighth-, and tenth-order susceptibilities display a robust positive peak near the critical point, whereas the frequently discussed negative dip is not universal for small mapping angle $\bar \ell$1; its presence depends on the mapping parameters and on the distance to the phase-transition line, and sub-leading critical contributions enhance it only modestly (Pan et al., 1 Apr 2025).

This result constrains how net quark number gain and related observables should be interpreted in critical-point searches. A positive peak in high-order susceptibilities is a stable critical signature in the summarized analysis, while a negative dip is not. Claims of universality for the latter are therefore not supported in that setting (Pan et al., 1 Apr 2025).

A complementary route uses imaginary chemical potential. The Fourier coefficients of the net-baryon number density,

$\bar \ell$2

encode the singularity structure of the QCD partition function. Their large-$\bar \ell$3 behavior is governed by Lee-Yang edge singularities, and scaling fits in a narrow temperature interval below the Roberge-Weiss transition temperature yield a Lee-Yang edge singularity at $\bar \ell$4 from the quoted lattice analysis (Schmidt, 2023). In this framework, the net-baryon response to chemical potential is read through the analytic structure of the partition function itself.

Model studies further separate chiral and deconfinement sensitivity. In the Hybrid Quark-Meson-Nucleon model, higher cumulants of the net-baryon number are substantially enhanced around the chiral phase transition but are not as sensitive to the deconfinement transition, because the auxiliary bag field governing quark-hadron switching is massive and does not act as a critical mode (Marczenko et al., 2017). In FRG studies of the quark-meson model, the pion mass $\bar \ell$5 emerges as a natural soft momentum scale at which cumulants saturate at their critical values, while for momentum scales larger than $\bar \ell$6 the characteristic $\bar \ell$7 structure of higher-order cumulants is lost (Morita et al., 2014).

Taken together, these results show that critical information in the net-quark sector is distributed across several observables: probe-induced gain, cumulants, Fourier coefficients, and momentum-resolved responses. They do not all diagnose the same physics with the same robustness.

5. External fields, anisotropy, and lattice computation

External conditions can qualitatively reorganize the conserved-charge sector. In a strongly magnetized QCD medium, rotational symmetry is broken, the pressure becomes anisotropic, and the second-order quark number susceptibility becomes anisotropic as well. The longitudinal and transverse susceptibilities are defined by

$\bar \ell$8

and the summarized result is that fluctuations are larger along the magnetic-field direction (Karmakar et al., 2021). In such a medium, any interpretation of net quark number response must account for directional dependence.

On the lattice, fugacity expansion provides a direct route to generalized susceptibilities at finite chemical potential. In this approach,

$\bar \ell$9

and quark-number observables are reconstructed from moments of the canonical determinants. Results quoted for $$\chi_n^q = \frac{1}{V\beta}\frac{\partial^n \ln Z_\mu}{\partial \mu^n},$$0 Wilson fermions give generalized susceptibilities up to fourth order, with ratios agreeing with the Hadron Resonance Gas model below the crossover and with the Stefan–Boltzmann limit above it, up to $$\chi_n^q = \frac{1}{V\beta}\frac{\partial^n \ln Z_\mu}{\partial \mu^n},$$1 (Gattringer et al., 2014).

These formulations do not define net quark number gain in the strict Polyakov-loop sense, but they provide the computational environment in which conserved-charge response is quantified. A plausible implication is that a complete phase-diagram analysis should treat probe observables and bulk susceptibilities as complementary, not competing, diagnostics.

6. Nonequilibrium production, hadronization, and heavy-ion phenomenology

In nonequilibrium settings, the phraseology of “gain” must be handled carefully. Real-time lattice simulations of a longitudinally expanding QCD plasma starting from a highly occupied gluonic state with vacuum quark sector show that the total quark number, after an initial rapid increase, grows almost linearly with time, and that the growth rate is consistent with a kinetic-theory estimate based on $$\chi_n^q = \frac{1}{V\beta}\frac{\partial^n \ln Z_\mu}{\partial \mu^n},$$2 scattering in the small-angle approximation (Tanji et al., 2017). However, in that setup the produced quarks and antiquarks come in pairs, so the net baryon number remains zero. This is quark production, not net quark number gain in the probe-induced sense.

By contrast, in non-equilibrium hadronization models the net quark number is conserved throughout expansion, constituent-mass generation, recombination, and freeze-out:

$$\chi_n^q = \frac{1}{V\beta}\frac{\partial^n \ln Z_\mu}{\partial \mu^n},$$3

The dynamics redistribute fixed net quark content among baryons, antibaryons, and mesons, but no net quark number is gained or lost (Zschocke et al., 2011). In quark combination models, generating-function methods make this mapping explicit and allow the calculation of final-state net-proton cumulants from initial quark and antiquark content (Yang et al., 2019).

Heavy-ion phenomenology adds yet another layer. As $$\chi_n^q = \frac{1}{V\beta}\frac{\partial^n \ln Z_\mu}{\partial \mu^n},$$4 decreases, net baryon number at midrapidity increases because baryon stopping transports entrance-channel quarks to midrapidity. In a two-component coalescence picture, this can violate simple constituent quark number scaling and generate systematic differences such as $$\chi_n^q = \frac{1}{V\beta}\frac{\partial^n \ln Z_\mu}{\partial \mu^n},$$5 without requiring the absence of partonic degrees of freedom (Dunlop et al., 2011). Experimentally, net-proton cumulant ratios are used as proxies for net-baryon and net-quark fluctuations; one recent model study describes BES-I high-energy data with a strongly coupled QGP description and interprets the absence of the expected critical variation in recent low-energy BES-II data as evidence for a weakly coupled, valence-quark-dominated regime at lower collision energies (Mamo, 2024).

The literature therefore juxtaposes several distinct quantities: probe-induced net quark number gain, bulk susceptibilities, total nonequilibrium quark production, and conserved-charge transport to midrapidity. This suggests that the term “net quark number gain” is most precise when reserved for the Polyakov-loop-based response to a static source, while related fluctuation and transport observables should be identified by their own operational definitions.