Canonical statistical hadronization with local baryon conservation for higher-order cumulants

This paper establishes that local baryon conservation within a canonical statistical hadronization framework can drive higher-order net-proton cumulant ratios to small or negative values in restricted rapidity acceptances, necessitating careful accounting of this baseline to avoid misinterpreting upcoming LHC measurements as signals of chiral criticality.

The Gist

Imagine a giant, chaotic party where thousands of guests (particles) are dancing in a massive hall (the collision zone). In high-energy physics, scientists crash two heavy nuclei together to create this "fireball" of particles. One of the most important rules of this party is Conservation of Baryon Number. Think of "baryons" as the VIP guests (like protons and neutrons). The rule says: The total number of VIPs minus the number of anti-VIPs must always stay the same. You can't just create a VIP out of thin air, and you can't make them disappear without a trace.

This paper is about understanding how this strict "VIP rule" changes the way we count the guests, especially when we only look at a small corner of the dance floor.

The Problem: The "Global" vs. "Local" View

Imagine you are a security guard trying to count how many VIPs are in a specific room.

• The Old Way (Global Conservation): The guard assumes that if a VIP enters the room, an anti-VIP must have left the entire building somewhere, even if that building is huge and the exit is on the other side of the world. This assumes the whole party is one giant, connected unit.
• The New Way (Local Conservation): The guard realizes that in reality, if a VIP enters the room, the anti-VIP that balances them out is likely standing right next to them, or at least in the same hallway. They are "locally" balanced.

The authors of this paper argue that for high-energy collisions (like those at the Large Hadron Collider, or LHC), the "Local" view is much more accurate. If you assume the balancing act happens instantly across the whole universe, you get the wrong math. If you assume the balancing happens in a small neighborhood (a few meters in "rapidity space"), the math changes significantly.

The Analogy: The Gaussian "Balancing Act"

The authors use a clever mathematical tool called a Gaussian Kernel. Think of this as a "blur" or a "smear."

• If you have a VIP at point A, the "anti-VIP" isn't just at point A. It's smeared out in a bell-curve shape around A.
• The width of this bell curve is called $\sigma_\eta$.
  • Narrow curve: The anti-VIP is very close (Ultra-local).
  • Wide curve: The anti-VIP can be further away (approaching the global view).

The paper calculates what happens when you count guests in a specific window (the "acceptance") while this "smearing" effect is happening.

The Big Surprise: The "Negative" Number

The most exciting discovery in the paper is about higher-order cumulants.

• Simple Analogy: Imagine you are measuring the "wiggliness" of the crowd.
  • 2nd Order: How much does the crowd size vary? (Standard deviation).
  • 4th Order & 6th Order: How "spiky" or "bumpy" is the distribution? Are there extreme outliers?

Scientists have been looking for a specific signal in these "spikiness" measurements. They believe that if the matter created in the collision undergoes a special phase change (called Chiral Criticality, related to how particles get their mass), the 6th-order measurement ($\kappa_6$) should turn negative.

The Paper's Warning:
The authors found that you don't need a special phase change to get a negative number.
Even if the party is just a boring, normal gas of particles (an "Ideal Gas"), the simple act of Local Baryon Conservation can naturally drive that 6th-order number down to zero or even negative values, if you are only looking at a small section of the dance floor.

Why this matters:
If scientists see a negative number in their data, they might scream, "We found the Chiral Critical Point!" But this paper says, "Wait! It might just be because of the local VIP rule. You have to subtract this 'boring' effect first before you can claim you found something new."

The Tools and Results

1. Better Math: They generalized the math to handle up to the 6th order (previously, most people only looked at the 2nd or 4th). They proved their math matches perfectly with a different method called the "Diffusion Master Equation" (which models how particles slowly drift apart over time).
2. The "Box" vs. The "Gaussian": Previous models used a "Box" approach (assuming the balancing happens perfectly within a sharp, hard-edged box). The authors show that a "Gaussian" (smooth, bell-curve) approach is more realistic and gives different results, especially when you look at larger areas of the fireball.
3. Predictions for O-O and Pb-Pb Collisions: They made specific predictions for upcoming experiments with Oxygen-Oxygen (O-O) and Lead-Lead (Pb-Pb) collisions at the LHC.
   • They provide a "Baseline": A set of numbers that represents what we expect to see if only conservation laws are at play, with no exotic physics.
   • They show that for the 6th-order ratio, the "Local" baseline can be negative, while the "Global" baseline stays positive.

The Takeaway

This paper is a "reality check" for experimental physicists. It says: "Before you celebrate finding a new state of matter, make sure you've correctly accounted for the fact that VIPs and Anti-VIPs tend to stick together locally."

If you ignore this local balancing act, you might mistake a simple mathematical consequence of conservation for a revolutionary discovery. The authors have provided the precise "correction factor" (the baseline) that future experiments need to use to ensure their discoveries are real.

Technical Summary

Technical Summary: Canonical Statistical Hadronization with Local Baryon Conservation for Higher-Order Cumulants

Problem Statement
Fluctuations of conserved charges, particularly net-proton cumulants, serve as critical probes for the QCD phase structure, including the search for the chiral critical point and the chiral crossover transition. While lattice QCD predicts a smooth crossover at the LHC ($\mu_B \approx 0$) characterized by specific signatures in higher-order susceptibilities (e.g., a drop in $\chi_4/\chi_2$ and a negative $\chi_6$), experimental measurements are complicated by non-critical effects. The most significant of these is the suppression of fluctuations due to exact baryon number conservation.

Previous approaches to modeling this conservation have limitations. The global canonical ensemble assumes conservation over the entire fireball, creating unphysical correlations between causally disconnected rapidity regions. Conversely, the "correlation volume" ($V_c$) approach, which enforces conservation within a sharp rapidity window, fails to describe intermediate cases where correlations spread over a finite range and neglects hadrons outside the window. Furthermore, existing models for local conservation have largely been restricted to second-order cumulants or coordinate-space analyses, lacking a systematic framework for higher-order moments (up to 6th order) and realistic kinematic acceptance.

Methodology
The authors generalize the density correlator formalism for local baryon conservation to higher orders ($n \le 6$). The core of the methodology involves:

1. Density Correlator Formalism: The $n$-point density correlators $C_n$ are decomposed into local (self-correlation) terms and non-local balancing terms required to satisfy global conservation constraints. The authors derive explicit expressions for $C_n$ up to $n=6$, generalizing previous work to include arbitrary local conservation kernels $\kappa_n$.
2. Gaussian Local Conservation: Instead of a sharp $V_c$ window, the authors model local conservation using a Gaussian kernel in spatial rapidity space. The two-point kernel $\kappa_2$ is defined as a Gaussian with width $\sigma_\eta$.
3. Finite-Volume Boundary Conditions: To address the limitations of periodic boundary conditions (which create unphysical "echoes" at fireball edges), the authors implement reflecting boundary conditions. This is achieved via image summation using the Poisson summation formula, ensuring that correlations decay physically toward the fireball boundaries.
4. Higher-Order Kernels: For $n \ge 3$, a symmetric Gaussian-cluster ansatz is adopted. This assumes $n$ coordinates are mutually correlated via a Gaussian weight depending on all pairwise distances, controlled by a single width parameter $\sigma_\eta$ shared across all orders.
5. Acceptance Integration: The formalism derives closed-form acceptance formulas for net-baryon and species-decomposed cumulants. By utilizing a "common-source representation," the $n$-dimensional integrals over the acceptance function are reduced to one-dimensional quadratures, facilitating efficient calculation of cumulants up to the 6th order.
6. Phenomenological Application: The coordinate-space results are mapped to momentum space using the blast-wave model to simulate kinematic cuts (pseudorapidity and transverse momentum) relevant for ALICE measurements in O–O and Pb–Pb collisions.

Key Contributions

• Analytical Derivation: The paper provides the first closed-form expressions for 5th and 6th order density correlators with local conservation kernels, extending the formalism beyond the previously available 4th order.
• Refined Boundary Conditions: The implementation of reflecting boundary conditions for Gaussian kernels resolves unphysical artifacts found in periodic or minimum-image conventions, particularly for higher-order cumulants involving multiple coordinates.
• Equivalence to Diffusion: The authors demonstrate that their static Gaussian local-conservation model yields results in exact agreement with the diffusion master equation approach (Ref. [21]) for all cumulant ratios up to $\kappa_6/\kappa_2$, provided the kernel width is matched to the diffusion length ($\sigma_\eta = \sqrt{2}d$). This establishes the Gaussian model as an analytically tractable static representation of diffusive broadening.
• Baseline Predictions: The study establishes a non-critical baseline for net-proton cumulants in O–O and Pb–Pb collisions at the LHC, quantifying the ideal-gas effects of local baryon conservation.

Results

• Cumulant Ratios: In the coordinate-space limit, the ratios $\kappa_4/\kappa_2$ and $\kappa_6/\kappa_2$ exhibit a characteristic double-minimum structure as a function of acceptance fraction $\alpha$. In the limit of very local conservation ($\sigma_\eta \to 0$), these ratios approach constant plateau values independent of $\alpha$.
• Negative $\kappa_6$: A significant finding is that local baryon conservation alone can drive the ratio $\kappa_6/\kappa_2$ to small or even negative values in restricted acceptance. Specifically, in the small-$\sigma_\eta$ limit, the analytical plateau value is calculated to be $\approx -0.025$.
• Comparison with $V_c$: The Gaussian approach agrees with the $V_c$ approach at small acceptance fractions ($\alpha \lesssim 0.1$), relevant for midrapidity measurements. However, for larger acceptances, the Gaussian model predicts a smooth interpolation between local and global regimes, whereas the $V_c$ model drops abruptly to zero once the acceptance exceeds the conservation window.
• Species Decomposition: The study shows that while net-baryon fluctuations are suppressed, the total baryon-plus-antibaryon variance ($\kappa_2[B+\bar{B}]$) remains Poissonian (unity) in the ideal gas limit, regardless of acceptance. This provides a distinct signature to test conservation mechanisms.
• Experimental Predictions: For upcoming LHC Run 3 measurements (O–O and Pb–Pb), the model predicts that $\kappa_6/\kappa_2$ can range from $\approx 0.70$ (global conservation) to $\approx -0.01$ (ultra-local) at full ALICE acceptance. The separation between global and local scenarios is most pronounced in higher-order ratios.

Significance
The paper argues that the observation of negative $\kappa_6$ or specific suppression patterns in net-proton cumulants cannot be immediately attributed to chiral criticality without a rigorous accounting of the conservation baseline. The authors emphasize that local baryon conservation alone is sufficient to drive $\kappa_6/\kappa_2$ to negative values in restricted acceptance. Therefore, precise theoretical baselines that incorporate local conservation effects are essential for interpreting upcoming high-precision LHC data. The developed formalism provides a systematic, analytically tractable framework to establish these baselines and disentangle conservation effects from potential critical phenomena.