Fluctuations and correlations of conserved charges in the Polyakov chiral SU(3) quark mean field model

This paper computes generalized susceptibilities and correlations of conserved charges up to eighth order within the Polyakov chiral SU(3) quark mean field model, comparing vacuum-inclusive and no-sea variants to reveal distinct behaviors in chiral-deconfinement splitting, off-diagonal correlator amplitudes, and higher-order kurtosis ratios as functions of temperature and baryon chemical potential.

The Gist

Imagine the universe as a giant, cosmic soup. In the very early moments after the Big Bang, or inside the heart of colliding heavy atoms in a particle accelerator, this soup is so hot and dense that the basic building blocks of matter—protons and neutrons—melt apart into a "quark-gluon plasma." It's like ice melting into water, but instead of water, you have a swirling sea of tiny, free-floating particles called quarks.

This paper is a detailed recipe book for understanding how this cosmic soup behaves when you change the temperature or the "pressure" (specifically, the density of matter) inside it. The authors, Dhananjay Singh and Arvind Kumar, use a sophisticated mathematical model called the Polyakov chiral SU(3) quark mean field (PCQMF) model to predict how the soup reacts.

Here is a breakdown of their work using simple analogies:

1. The Two Main Transitions: Unfreezing and Un-gluing

In this cosmic soup, there are two major changes happening as things cool down:

• Chiral Symmetry Breaking (Unfreezing): Think of quarks as dancers. At high temperatures, they are free to dance anywhere. As it cools, they pair up and get "stuck" in a specific formation (forming protons and neutrons). This is like the soup freezing into a solid block.
• Deconfinement (Un-gluing): This is when the "glue" holding the quarks together breaks. At high heat, the glue snaps, and quarks roam free. At lower heat, the glue holds them tight.

The authors wanted to see if these two events happen at the exact same time or if they are slightly separated, like two doors opening one after the other.

2. The Secret Ingredient: The "Vacuum" Term

The most important part of this study is testing two different versions of their recipe:

• Version A (vac=1): Includes the "fermion vacuum term." Imagine this as accounting for the "background noise" or the invisible energy of empty space that still affects the particles. It's like realizing that even when a room is empty, the air pressure and temperature still exist and affect how a balloon behaves.
• Version B (vac=0): Ignores this background energy. It's a simpler recipe that assumes empty space is truly nothing.

The authors found that including this "background noise" (Version A) changes the results significantly. It makes the transition between the "stuck" and "free" states sharper and creates a clearer separation between the two "doors" (the chiral and deconfinement transitions).

3. Measuring the "Fluctuations" (The Soup's Jitters)

To understand the soup, the scientists didn't just look at the average temperature; they looked at the fluctuations or "jitters."

• Imagine a crowd of people. If everyone is calm, the crowd is still. If they are excited, they jostle and bump into each other.
• The authors calculated how much the "charge" (like electric charge or the number of baryons) jitters around. They looked at these jitters up to the eighth order.
  • Analogy: If the "first order" is just the average number of people in a room, the "second order" is how much that number wiggles. The "eighth order" is looking at incredibly complex, subtle patterns in how the crowd moves—like detecting a specific rhythm in the jostling that only happens right before the crowd breaks into a dance.

4. Key Findings: What the "Vacuum" Changed

• Splitting the Transitions: When they included the "vacuum" term, they saw a clear gap between the two transitions. The "unfreezing" happened at a slightly different temperature than the "un-gluing." Without the vacuum term, these two events looked more like they were happening at the same time.
• Twin Peaks: When they looked at the complex "jitters" (higher-order fluctuations), the version with the vacuum term showed twin peaks (two distinct humps) in the data. This is like hearing two distinct drumbeats instead of one long thud. This proves the two transitions are separate events.
• Strange Quarks: They also looked at "strange" particles (a heavier type of quark). They found that the "vacuum" version was better at describing the behavior of light particles, while the "no-vacuum" version surprisingly did a better job describing the behavior of the heavy "strange" particles when they were melting.

5. Comparing to Reality (Lattice QCD)

The authors compared their mathematical soup to data from Lattice QCD, which is like a super-computer simulation of the universe that acts as the "gold standard" or the actual measurement.

• Their model generally matched the trends seen in the super-computer data.
• However, like any model, it had some limitations. For example, it underestimated the "jitters" of electric charge at low temperatures because the model treats pions (light particles) as frozen statues rather than wiggly, active particles.

6. Pushing the Limits (High Density)

Finally, they tested what happens if you squeeze the soup even harder (increasing the density of matter, or $\mu_B$).

• They found that as the density increases, the "jitters" become wilder and more complex.
• One specific ratio they measured (related to how "spiky" the distribution of particles is) turned negative in the version with the vacuum term, but stayed positive in the version without it. This is a crucial difference that could help experimentalists at facilities like the RHIC (Relativistic Heavy Ion Collider) figure out which version of the physics is correct.

Summary

In short, this paper is a deep dive into the "recipe" for the early universe's soup. The authors discovered that including the "background energy" of empty space (the vacuum term) makes the model more realistic. It reveals that the transition from free quarks to bound matter happens in two distinct steps, and it creates unique, complex patterns in how particles fluctuate. These patterns serve as a fingerprint that scientists can look for in real-world experiments to understand the fundamental nature of matter.

Technical Summary

Technical Summary: Fluctuations and Correlations of Conserved Charges in the Polyakov Chiral SU(3) Quark Mean Field Model

Problem and Motivation
Quantum Chromodynamics (QCD) at finite temperature ($T$) and baryon chemical potential ($\mu_B$) describes the phase structure of matter relevant to the early universe and relativistic heavy-ion collisions. While lattice QCD simulations at $\mu_B = 0$ have established that chiral symmetry restoration and deconfinement occur as a smooth crossover, the behavior at finite $\mu_B$ remains constrained by the fermion sign problem. Experimental efforts, such as the RHIC Beam Energy Scan (BES), measure higher-order cumulants of net-proton, net-kaon, and net-charge distributions to probe critical dynamics and locate a potential critical endpoint (CEP).

Theoretical effective models are essential for interpreting these data. However, existing mean-field chiral models often suffer from deficiencies, such as underestimating second-order charge susceptibilities ($\chi_Q^2$) due to the "freezing" of pion fluctuations, and lacking comprehensive calculations of the full set of generalized baryon-charge-strangeness (BQS) susceptibilities. Furthermore, the role of the fermion vacuum (Dirac-sea) term in the Polyakov chiral SU(3) quark mean field (PCQMF) model has been largely unexplored regarding thermodynamic susceptibilities, despite its known importance in restoring universality classes and smoothing crossovers in related models like PNJL and PQM.

Methodology
The authors compute generalized susceptibilities of conserved charges ($B, Q, S$) within the PCQMF model. The model incorporates non-strange ($\sigma$) and strange ($\zeta$) scalar fields, an isovector scalar ($\delta$), a dilaton field ($\chi$) for the trace anomaly, and the Polyakov loop ($\Phi, \bar{\Phi}$) for deconfinement. Vector mesons ($\omega, \rho, \phi$) provide repulsion, though the vector coupling is set to zero ($g_v=0$) in the primary analysis.

The study compares two variants:

1. vac=1: Includes the fermion vacuum term (one-loop quark contribution) in the thermodynamic potential.
2. vac=0: A "no-sea" variant where the vacuum term is omitted, with model parameters independently refitted to the same vacuum observables (pion/kaon masses, $\sigma$-meson mass, vacuum stationarity).

The generalized susceptibilities $\chi_{ijk}^{BQS}$ are defined as derivatives of the scaled pressure with respect to reduced chemical potentials. The calculation involves solving coupled gap equations for mean fields and computing derivatives up to the eighth order for diagonal susceptibilities and all twelve independent fourth-order off-diagonal correlators.

Key Contributions
This work presents the first comprehensive calculation of the full set of generalized BQS susceptibilities in the PCQMF model including the fermion vacuum term. Specifically:

• At $\mu_B = 0$, it covers diagonal susceptibilities $\chi_{n}^{B,Q,S}$ up to $n=8$ and all twelve independent fourth-order off-diagonal correlators.
• At finite $\mu_B$ (up to 500 MeV, with $\mu_Q=\mu_S=0$), it computes diagonal susceptibilities up to eighth order for baryons and fourth order for charge/strangeness, along with second-order off-diagonals and odd-order baryon susceptibilities ($\chi_1^B, \chi_3^B, \chi_5^B$).
• It provides a direct comparison between the vacuum-included and no-sea variants to isolate the impact of the Dirac-sea term on fluctuation observables.

Results

• Order Parameters and Thermodynamics:

  • The chiral pseudocritical temperature ($T_{pc}$) is found to be 170.5 MeV (vac=1) and 166.4 MeV (vac=0) at $\mu_B=0$. The deconfinement temperature ($T_{dec}$) is 144.4 MeV (vac=1) and 146.6 MeV (vac=0).
  • The inclusion of the vacuum term (vac=1) leads to a steeper drop in the scalar field $\sigma$ and a distinct inflection in the subtracted chiral condensate derivative ($-d\Delta_{l,s}/dT$) near $T_{dec}$, a feature absent in vac=0.
  • The model underpredicts bulk thermodynamic quantities (pressure, energy density) at low $T$ due to the absence of hadronic degrees of freedom, a known limitation of mean-field quark models.

• Diagonal Susceptibilities:

  • Second Order: All diagonal susceptibilities ($\chi_2$) lie below lattice data at low $T$, with the deficit most severe for $\chi_Q^2$ due to the lack of pion fluctuations.
  • Higher Orders: The chiral-deconfinement splitting, induced by the glue-to-Yang-Mills temperature mapping, is resolved in higher-order derivatives. In vac=1, this manifests as twin maxima in $\chi_4^B$ and $\chi_6^Q$, twin minima in $\chi_8^B$ and $\chi_8^Q$, and multiple zero crossings in $\chi_6^B$. The amplitudes of these structures are significantly larger in vac=1 than in vac=0.
  • Charge/Strangeness: For charge ($\chi_4^Q$), vac=1 amplitudes exceed vac=0. For strangeness ($\chi_4^S, \chi_6^S, \chi_8^S$), the hierarchy reverses in the strange-melting region, where vac=0 dominates.

• Off-Diagonal Correlators:

  • Among the twelve fourth-order off-diagonal correlators, the BQ (baryon-charge) components show vac=1 dominance across the chiral crossover.
  • Conversely, the BS (baryon-strangeness), QS (charge-strangeness), and mixed BQS components peak near the strange-melting temperature, where vac=0 amplitudes exceed vac=1.
  • The $\chi_{11}^{BQ}$ correlator tracks lattice data through the peak, while $\chi_{11}^{BS}$ and $\chi_{11}^{QS}$ are generally overestimated by the model compared to lattice results.

• Finite $\mu_B$ Behavior:

  • As $\mu_B$ increases, the peak amplitudes of susceptibilities grow, with fourth-order moments growing much faster than second-order ones.
  • The chiral-deconfinement splitting becomes more pronounced, with twin peaks merging into single sharp peaks in some channels as $T_{pc}$ drops faster than $T_{dec}$.
  • Odd-order baryon susceptibilities ($\chi_1^B, \chi_3^B, \chi_5^B$) emerge at finite $\mu_B$, exhibiting rich oscillatory structures in vac=1.

• Cumulant Ratios:

  • The kurtosis ratio $R_4^B \equiv \chi_4^B / \chi_2^B$ along the pseudocritical line starts positive but crosses zero at $\mu_B/T_{pc} \approx 2.15$ in the vac=1 variant, remaining positive for vac=0.
  • Higher-order ratios $R_5^B$ and $R_6^B$ start negative in vac=1 and become increasingly negative with $\mu_B$, whereas they start marginally positive in vac=0 before turning negative at higher $\mu_B$.
  • The vac=1 amplitudes for these ratios significantly exceed those of vac=0.

Significance and Claims
The paper claims to provide the first complete set of generalized BQS susceptibilities in the PCQMF model with the fermion vacuum term. The authors assert that the vacuum term is crucial for capturing the correct universality class in the chiral limit and for generating specific structural features in higher-order susceptibilities, such as the inflection in the condensate derivative and the enhanced amplitude of oscillatory structures.

The study highlights that the vacuum term does not uniformly enhance all observables; rather, it creates a channel-dependent hierarchy where vac=1 dominates chiral crossover features (BQ sector) while vac=0 dominates strange-melting features (BS, QS, BQS sectors). The authors note that while the model reproduces qualitative crossover behaviors seen in lattice QCD, quantitative discrepancies (such as the underestimation of $\chi_Q^2$ and the overshoot of higher-order ratios) persist, attributing these to the mean-field approximation's neglect of mesonic fluctuations. The work establishes a baseline for future comparisons with experimental data along the heavy-ion freeze-out trajectory, a task reserved for subsequent studies.