Search for the QCD Critical Point in High Energy Nuclear Collisions: A Status Report

This paper reviews recent STAR experiment results on net-proton multiplicity fluctuations from RHIC BES-II collisions to search for the QCD critical point by comparing fourth-order cumulant and factorial cumulant ratios with non-critical theoretical models while addressing initial volume fluctuations and outlining future research directions.

The Gist

Imagine the universe as a giant, cosmic kitchen. For the first few microseconds after the Big Bang, the ingredients of matter (quarks and gluons) were swimming freely in a super-hot, super-dense soup called the Quark-Gluon Plasma (QGP). As the universe cooled, this soup froze into solid "ice cubes" called protons and neutrons, which make up the atoms in our world today.

Physicists believe that if you could heat up this "ice" again and squeeze it with enough pressure, it wouldn't just melt smoothly back into soup. Instead, there might be a specific spot on the map of temperature and pressure where the transition gets messy and chaotic. This spot is called the QCD Critical Point. Finding it is like finding the exact moment water turns to steam, but for the fundamental building blocks of the universe.

This paper is a status report from the STAR experiment at the Relativistic Heavy Ion Collider (RHIC). Here's what they are doing, explained simply:

1. The Experiment: Smashing Gold Balls

Think of the RHIC as a giant particle accelerator, essentially a very fast racetrack. The scientists smash gold atoms (Au+Au) together at nearly the speed of light.

• The Goal: By changing the speed of the crash, they can control how "hot" and how "dense" the resulting fireball is.
• The Range: They have been smashing these atoms at energies ranging from very high (200 GeV) down to very low (3 GeV). This covers a huge range of the "map" where the Critical Point might be hiding.
• The New Gear: Recently, they upgraded their detectors (like adding better cameras and sensors) to look at the lowest energy crashes (the "fixed-target" mode), which creates the densest, most pressure-heavy conditions.

2. The Clue: Counting the "Net-Protons"

When the gold balls smash, they create a shower of new particles. The scientists can't see neutrons easily, so they focus on protons (which are positively charged).

• The Analogy: Imagine you are at a crowded party. You want to know if the crowd is behaving normally or if something strange is happening. You count the number of people wearing red hats (protons) versus blue hats (anti-protons).
• The Fluctuation: In a normal, calm party, the number of red hats varies slightly from one room to another, but it's predictable. However, if the party is near a "Critical Point" (like a dance floor about to collapse), the number of red hats might suddenly jump up and down wildly.
• The Math: The scientists use complex math (called "cumulants") to measure these wild swings. They look at the average, the spread, and the "skewness" of the counts.

3. The Results: A Bumpy Road

The paper reports on what they found when they plotted these numbers against the collision energy:

• High Energy (The Smooth Part): At high speeds, the data matches the "boring" models. The fluctuations are calm and predictable, just like a smooth highway.
• Low Energy (The Bump): As they slowed down the collisions to create higher density, the data started to act weird.
  • At around 20 GeV, the data showed a "bump" or a deviation from the smooth models. This is a potential "smoking gun" for the Critical Point. The significance of this bump is about 2 to 5 times the expected noise level (a 2-5 sigma signal), which is exciting but not yet a confirmed discovery.
  • At the lowest energies (below 10 GeV), the data started rising again in ways the standard models couldn't explain. This suggests that at high densities, protons start "attracting" each other in a way they don't at high speeds.

4. The Challenge: The "Volume" Problem

One of the biggest headaches in this experiment is the Volume Fluctuation.

• The Metaphor: Imagine trying to measure the average height of people in a room. If you accidentally measure a room that is slightly bigger in one shot and smaller in another, your average height calculation will be wrong, even if the people didn't change.
• The Fix: In heavy-ion collisions, the "size" of the fireball changes slightly from crash to crash. The paper discusses a new, clever mathematical method to "subtract" this size effect so they can see the true signal of the particles themselves. They tested this on computer simulations (UrQMD) and it works well.

5. What's Next?

The story isn't over.

• More Data: The STAR experiment is currently gathering final results for the lowest energy crashes (3 to 4.5 GeV). This is the "high pressure" zone where the Critical Point is most likely to be found.
• New Competitors: Other labs around the world (like FAIR in Germany and NICA in Russia) are building their own "gold-smashing" machines. They will cover the same low-energy, high-density territory, allowing scientists to cross-check their findings.
• China's HIAF: A new facility in Huizhou, China, is being built to do similar physics with even higher intensity beams.

The Bottom Line

This paper is a "work in progress" report. The scientists have found some very interesting "bumps" in the data that suggest the QCD Critical Point might be real and located in the high-density, low-energy region. However, they are being careful. They know that other effects (like the size of the collision or how particles move) can fake these signals.

They are currently refining their tools and gathering more data to confirm if they have truly found the "Holy Grail" of nuclear physics: the point where matter changes its state in the most dramatic way possible.

Technical Summary

1. Problem Statement

The primary goal of heavy-ion collision physics is to map the phase diagram of Quantum Chromodynamics (QCD), specifically locating the QCD Critical Point (CP).

• Theoretical Context: Lattice QCD predicts a smooth crossover transition between Quark-Gluon Plasma (QGP) and hadronic matter at low baryon chemical potential ($\mu_B$). However, at high $\mu_B$, the nature of the transition is unknown due to the fermion sign problem in lattice calculations. Many theoretical models predict a first-order phase transition at high $\mu_B$ terminating at a second-order critical point.
• Experimental Challenge: The STAR experiment at the Relativistic Heavy Ion Collider (RHIC) conducts the Beam Energy Scan (BES) to locate this point by measuring fluctuations of conserved quantities (net-baryon number).
• Key Obstacles:
  • Interpretation: Fluctuation observables are sensitive to critical phenomena but are also heavily influenced by non-critical effects such as finite-size effects, critical slowing down, baryon transport, and initial volume fluctuations (IVF).
  • Low-Energy Limitations: In fixed-target mode (low collision energies, high $\mu_B$), limited centrality resolution due to smaller reference multiplicities introduces significant volume fluctuation effects, complicating the extraction of critical signals.
  • Model Discrepancies: Experimental data often deviate from non-critical baseline models (e.g., UrQMD, HRG), but distinguishing these deviations from critical signatures requires rigorous statistical analysis and improved methodology.

2. Methodology

The paper reviews data from the STAR experiment's BES-II program, covering collider mode ($\sqrt{s_{NN}} = 7.7 - 27$ GeV) and fixed-target mode ($\sqrt{s_{NN}} = 3.0 - 3.9$ GeV).

• Observables: The study focuses on net-proton multiplicity fluctuations (used as a proxy for net-baryon number) and proton factorial cumulants.
  • Cumulants ($C_n$): Calculated up to the 4th order ($C_1$ to $C_4$).
  • Factorial Cumulants ($\kappa_n$): Calculated up to the 4th order.
  • Ratios: Ratios like $C_2/\langle N \rangle$, $C_3/C_1$, $C_4/C_2$, and factorial equivalents are used to cancel volume dependence and compare with theoretical susceptibilities.
• Data Analysis Techniques:
  • Centrality Selection: Events are classified into centrality bins (e.g., 0-5%, 70-80%) based on charged hadron multiplicity.
  • Corrections:
    • Centrality Bin Width Correction (CBWC): Applied to suppress volume fluctuations within broad centrality bins.
    • Efficiency Correction: Analytical corrections assuming a binomial detector response are applied to account for tracking and particle identification efficiency.
  • Baseline Comparison: Experimental ratios are compared against non-critical theoretical models:
    • Hydrodynamic models with excluded-volume corrections (Hydro EV).
    • Hadron Resonance Gas with canonical ensemble (HRG CE).
    • UrQMD transport model.
    • Lattice QCD calculations.
• Novel Approach (Sec 3.2): To address the severe volume fluctuation issues in low-energy fixed-target data, the authors propose a centrality-independent method:
  1. Parameterize cumulants as polynomial functions of $C_1$ (mean multiplicity).
  2. Ensure weights of proton distributions across centrality bins are similar.
  3. Reconstruct proton distributions using the Edgeworth expansion based on cumulant information.
  4. Optimize parameters using a differential evolution algorithm integrated with Bayesian optimization to match reconstructed distributions with experimental data.

3. Key Contributions

1. Comprehensive Review of BES-II Data: The paper synthesizes the latest results from both collider and fixed-target modes, extending the energy coverage down to $\sqrt{s_{NN}} = 3.0$ GeV.
2. Introduction of Factorial Cumulants: It highlights the utility of proton factorial cumulant ratios ($\kappa_n$) as complementary probes to standard cumulants, particularly for disentangling volume effects.
3. Development of IVF Suppression Method: A significant technical contribution is the proposal and validation of a centrality-independent method to suppress Initial Volume Fluctuations. This method was tested in UrQMD simulations and shown to reproduce results comparable to $N_{part}$-based centrality (which is difficult to measure experimentally).
4. Efficiency Correction Validation: The authors demonstrated the validity of their efficiency correction method in UrQMD simulations, showing that corrected results align closely with true values.

4. Key Results

• High Energy ($\sqrt{s_{NN}} > 11$ GeV):
  • Net-proton cumulant ratios ($C_2/\langle N \rangle$, $C_3/C_1$, $C_4/C_2$) show a monotonic decrease as energy decreases from 200 GeV, consistent with non-critical models.
  • Critical Deviation: A maximum deviation of the $C_4/C_2$ ratio from non-critical baselines is observed at $\sqrt{s_{NN}} \approx 19.6$ GeV in central (0-5%) collisions. The significance of this deviation reaches 2 to 5$\sigma$, consistent with the predicted signature of the QCD critical point.
• Low Energy ($\sqrt{s_{NN}} < 11$ GeV):
  • Rising Trends: Below 10 GeV, lower-order ratios ($C_2/\langle N \rangle$, $C_3/C_1$, $\kappa_2/\kappa_1$, $\kappa_3/\kappa_1$) exhibit a clear increase as energy decreases (moving to higher $\mu_B$).
  • Model Failure: Standard non-critical models (UrQMD, HRG) fail to reproduce this rising trend. This suggests the necessity of attractive interactions in the high baryon density region, a condition required for the formation of a critical point (e.g., in Van der Waals interactions).
  • Fixed-Target Data: New preliminary results from $\sqrt{s_{NN}} = 3.0 - 3.9$ GeV confirm the monotonic decrease of $C_4/C_2$ but show the rising behavior in lower-order ratios.
• Volume Fluctuation Impact: The study confirms that without proper suppression (like the new centrality-independent method), volume fluctuations significantly distort higher-order cumulants, particularly in mid-central and central collisions at low energies.

5. Significance and Outlook

• Evidence for Critical Point: The non-monotonic behavior of $C_4/C_2$ around 20 GeV and the rising trends of lower-order ratios at low energies provide the strongest experimental evidence to date for the existence of the QCD critical point in the high $\mu_B$ region.
• Methodological Advancement: The proposed centrality-independent method is crucial for future analyses of fixed-target data, where traditional centrality definitions are prone to large volume fluctuation errors.
• Future Directions:
  • Final Results: Final results from the STAR fixed-target program ($\sqrt{s_{NN}} = 3 - 4.5$ GeV) are needed to confirm the trend.
  • Global Effort: The paper outlines the upcoming capabilities of the CBM experiment at FAIR and the MPD experiment at NICA, which will cover the overlapping energy range ($2.0 \le \sqrt{s_{NN}} \le 11.2$ GeV) with higher statistics.
  • HIAF: The new High Intensity Heavy-ion Facility in China will offer unique capabilities for studying baryon polarizations and hyperon-nucleon interactions relevant to critical fluctuations.

In conclusion, the paper establishes that while deviations from non-critical baselines are observed, a robust search for the QCD critical point requires the integration of advanced dynamical modeling, finite-size scaling, and improved experimental techniques to suppress volume fluctuations, particularly in the high baryon density regime.