Proton High-Order Cumulants in Au+Au Collisions at High Baryon Density from JAM with a Centrality-Independent Framework

This study utilizes the JAM model and a novel Centrality-Independent Genuine Cumulant Analysis (CIGAR) framework to systematically analyze higher-order proton cumulants in Au+Au collisions at high baryon densities, providing a dynamic non-critical baseline essential for the search for the QCD critical point.

The Gist

Imagine the universe as a giant, chaotic kitchen. Physicists are trying to figure out the "recipe" of matter by smashing heavy atoms (like gold) together at incredibly high speeds. When these atoms crash, they create a tiny, super-hot soup called Quark-Gluon Plasma (QGP). Scientists believe that if they smash these atoms at just the right speed, they might find a "Critical Point"—a specific moment where the soup changes its state in a dramatic, non-linear way, similar to how water suddenly turns into steam.

To find this Critical Point, scientists look at the protons (tiny particles inside the atoms) that fly out of the crash. They don't just count them; they look at the patterns of how many protons appear in each crash. This is where Cumulants come in.

The Problem: The "Bumpy Table" Effect

Imagine you are trying to measure the average height of people in a room. But, every time you take a measurement, the room itself changes size slightly. Sometimes it's a small closet; sometimes it's a giant ballroom.

• If you count people in a closet, you get a small number.
• If you count people in a ballroom, you get a big number.

If you mix these measurements together without accounting for the room size, your data looks messy and misleading. In physics, this is called volume fluctuation. In heavy-ion collisions, the "room size" (the volume of the collision) changes slightly from one crash to the next. Traditional methods tried to fix this by grouping crashes into "bins" (like sorting them by how hard they hit), but at low energies, this method is like trying to sort a messy pile of sand with a sieve that has holes too big—it lets too much noise through.

The Solution: The "CIGAR" Tool

The authors of this paper developed a new, smarter way to measure these patterns called CIGAR (Centrality-Independent Genuine Cumulant Analysis Framework).

Think of CIGAR as a high-tech noise-canceling headphone for data. Instead of just sorting the crashes into bins, CIGAR uses a sophisticated mathematical recipe (combining "differential evolution" and "Bayesian inference") to reconstruct the true shape of the proton distribution. It effectively "subtracts" the background noise caused by the changing room sizes, leaving behind only the genuine signal of the physics happening inside the crash.

What They Found

Using a computer simulation called the JAM model (which acts like a virtual physics lab), the team smashed gold atoms together at specific low energies (3.2 to 4.5 GeV) and applied their new CIGAR tool.

1. CIGAR Works Better: When they compared CIGAR to the old method (CBWC), CIGAR produced results that were much cleaner and more consistent. It successfully removed the "bumpy table" noise, proving it's a reliable tool for finding the Critical Point in the future.
2. The "Spectator" Effect: In a collision, not every part of the atom hits the other one. Some parts just graze by, like a car that misses a crash but still spins out. These are called spectators. The paper found that these "grazing" parts mess up the data, especially at lower energies and when looking at a wide area of the crash. It's like trying to hear a whisper in a room where people are shouting in the hallway; the further you are from the center, the more the hallway noise interferes.
3. The Pattern: They found that the patterns of protons change depending on how "central" (head-on) the crash is. In the most central crashes, the patterns hit a peak and then drop off, likely because the rules of conservation (you can't create or destroy protons out of thin air) start to limit how much the numbers can fluctuate.

Why This Matters

This paper doesn't claim to have found the Critical Point yet. Instead, it provides a clean, reliable baseline.

Imagine you are looking for a rare, specific bird in a forest. If the forest is full of noise and wind (volume fluctuations), you might mistake a rustling leaf for a bird. This paper builds a better pair of binoculars (CIGAR) and maps out exactly how the wind behaves (the spectator effect) so that when real experiments (like those at RHIC or the future CBM experiment) look for the Critical Point, they know exactly what a "false alarm" looks like.

In short: The authors built a better mathematical filter to clean up the noise in heavy-ion collision data, allowing scientists to see the true signals of matter's phase changes more clearly, which is a crucial step toward finding the elusive QCD Critical Point.

Technical Summary

Technical Summary: Proton High-Order Cumulants in Au+Au Collisions at High Baryon Density from JAM with a Centrality-Independent Framework

Problem Statement
The primary objective of heavy-ion collision experiments is to map the QCD phase diagram and locate the QCD critical point (CP). High-order cumulants of conserved charges (net-baryon, net-electric charge, net-strangeness) are identified as sensitive observables for this search due to their dependence on the correlation length $\xi$, which diverges at the critical point. While the RHIC-STAR collaboration has measured net-proton cumulants across a broad energy range ($\sqrt{s_{NN}} = 3 - 200$ GeV), a significant challenge persists in the high-baryon-density regime (low energy, $\sqrt{s_{NN}} < 7.7$ GeV). In this regime, initial volume fluctuations significantly influence cumulant measurements derived using the traditional Centrality Bin Width Correction (CBWC) method. These fluctuations can obscure non-monotonic signals indicative of the critical point. Furthermore, the specific impact of spectator nucleons on cumulant measurements in this energy range requires systematic investigation to establish a reliable non-critical baseline.

Methodology
This study utilizes the Jet AA Microscopic Transport Model (JAM) to simulate Au+Au collisions at center-of-mass energies of $\sqrt{s_{NN}} = 3.2, 3.5, 3.9,$ and $4.5$ GeV. The simulations cover impact parameters from 0 to 14 fm, with particle acceptance matching the RHIC-STAR fixed-target detector configuration ($0.4 < p_T < 2.0$ GeV/c and varying rapidity windows).

To address the volume fluctuation issue, the authors employ a novel method termed the Centrality-Independent Genuine Cumulant Analysis Framework (CIGAR). Unlike the traditional CBWC approach, CIGAR reconstructs the particle number distribution by optimizing Edgeworth expansion parameters. This is achieved through a hybrid approach combining differential evolution and Bayesian inference. The framework fits the total distribution as a sum of sub-distributions corresponding to consecutive centrality bins (2.5% intervals), effectively eliminating initial volume fluctuations without relying on centrality binning corrections.

The study calculates:

1. Ordinary cumulants ($C_1$ to $C_4$) and their ratios ($C_2/C_1, C_3/C_1, C_4/C_2$).
2. Factorial cumulants ($\kappa_n$) and their ratios, which isolate genuine $n$-particle correlations.
3. The effects of spectator nucleons by comparing simulations with and without spectator contributions.
4. Dependence on rapidity window size ($\Delta y$) and centrality.

Key Contributions and Results

• Validation of CIGAR: The study demonstrates that the CIGAR method effectively suppresses initial volume fluctuations. When comparing CIGAR results against the $N_{part}$-based CBWC baseline (which serves as a reference with minimal volume fluctuations), the CIGAR results closely overlap across all four energies. In contrast, the traditional RefMult3-based CBWC method shows systematic deviations from the $N_{part}$ baseline, particularly at $\sqrt{s_{NN}} = 3.2$ GeV. This confirms that CBWC suffers from reduced centrality resolution at lower energies due to smaller reference multiplicities, whereas CIGAR maintains consistent performance.

• Cumulant Behavior:

  • Centrality Dependence: The first-order cumulant ($C_1$) increases approximately linearly with $N_{part}$, with slopes decreasing at higher energies, reflecting weakened baryon stopping. The fourth-order cumulant ($C_4$) exhibits a non-monotonic behavior, rising with $N_{part}$ to a maximum around $N_{part} \sim 150-200$ before decreasing in the most central collisions. This suppression is attributed to baryon number conservation effects when the acceptance window constitutes a small fraction of the total system volume.
  • Ratio Trends: All cumulant ratios ($C_2/C_1, C_3/C_1, C_4/C_2$) show a monotonic decrease from peripheral to central collisions. The magnitude of these ratios decreases with increasing collision energy, and the centrality dependence flattens at higher energies.

• Spectator Nucleon Effects: The systematic investigation of spectator contributions reveals a centrality- and energy-dependent effect.

  • For $C_2/C_1$, spectator removal enhances the ratio in central collisions but reduces it in peripheral collisions.
  • For $C_4/C_2$, spectator exclusion consistently lowers the ratio in both centrality classes, with the reduction most pronounced at 3.2 GeV and becoming negligible at 4.5 GeV.
  • Spectator effects are found to be most significant at low collision energies and for wide rapidity windows, diminishing at higher energies and in narrow windows.

• Factorial Cumulants: The factorial cumulant ratios ($\kappa_2/\kappa_1, \kappa_3/\kappa_1, \kappa_4/\kappa_1$) are approximately an order of magnitude smaller than ordinary cumulant ratios and remain close to zero. Spectator effects on these ratios are generally negligible, except for a modest enhancement in $\kappa_3/\kappa_1$ in central collisions.

Significance and Claims
The paper establishes a dynamic non-critical baseline for higher-order cumulants, factorial cumulants, and their ratios in the high-baryon-density regime. By validating the CIGAR method as a superior tool for eliminating volume fluctuations compared to traditional CBWC, the study provides a robust framework for interpreting experimental data.

The authors claim that their results are crucial for the search for the QCD critical point. Specifically, the established baselines allow for more accurate comparisons with experimental data from the RHIC-STAR fixed-target program ($\sqrt{s_{NN}} = 3.2 - 7.7$ GeV) and future experiments such as the Compressed Baryonic Matter (CBM) experiment at FAIR ($\sqrt{s_{NN}} = 2.4 - 4.9$ GeV). The study emphasizes that the spectator effect, which is pronounced at low energies and large rapidity windows, must be carefully accounted for when comparing theoretical baselines with experimental measurements to avoid misinterpreting non-critical background effects as signals of the critical point.