Eigen-microstate Signatures of Criticality in Relativistic Heavy-Ion Collisions

This paper introduces an eigen-microstate framework that analyzes final-state particle fluctuations in relativistic heavy-ion collisions to identify the QCD critical point by extracting dominant critical modes as a robust order parameter, offering a background-insensitive tool directly applicable to RHIC Beam Energy Scan data without relying on equilibrium assumptions.

The Gist

Imagine trying to understand a massive, chaotic crowd at a concert. Usually, if you look at the crowd, you see a jumble of people moving randomly. But sometimes, if a specific event happens (like a famous singer stepping on stage), the crowd might suddenly start moving in a synchronized wave, or forming distinct clusters.

This paper is about a new way to spot those "synchronized waves" or "clusters" in the chaotic aftermath of relativistic heavy-ion collisions. These are experiments where scientists smash heavy atoms together at near-light speeds to recreate the conditions of the early universe (a soup of quarks and gluons).

Here is the breakdown of their discovery using simple analogies:

1. The Problem: Finding a Needle in a Haystack

Scientists want to find a specific "Critical Point" in these collisions—a moment where the matter changes phase, similar to water turning into steam.

• The Old Way: Scientists used to look for specific signals (like counting how many particles of a certain type are made). But this is like trying to hear a whisper in a hurricane. The "noise" of the collision (random statistical fluctuations, particles decaying, etc.) is so loud that it drowns out the signal. You need millions of events to see anything, and even then, it's hard to be sure.
• The New Idea: Instead of listening for a whisper, the authors propose looking at the entire pattern of the crowd at once.

2. The Solution: The "Eigen-Microstate" Framework

The authors developed a new mathematical tool called the Eigen-Microstate Approach (EMA). Here is how it works, step-by-step:

• Step 1: Taking a Snapshot (The Microstate)
  Imagine every single collision is a unique photograph. In this photo, we don't just count the people; we look at exactly where they are standing and how they are moving. The authors call this a "microstate." It captures the unique "personality" of that one specific crash.

• Step 2: The Group Photo (The Ensemble)
  They take thousands of these snapshots and stack them together. They ask: "If we look at all these photos together, is there a common theme that keeps appearing?"

• Step 3: Finding the "Main Character" (The Eigen-Microstate)
  Using a method similar to how Netflix analyzes your viewing habits to find your "top genre," this math breaks down the thousands of chaotic photos into a few "Main Characters" (called Eigen-Microstates).

  • The "Background" Character: Most of the time, the "Main Character" is just random noise or standard physics (like the crowd shuffling randomly).
  • The "Critical" Character: If a critical point exists, a new Main Character emerges. This character represents a synchronized pattern (like the crowd suddenly forming a perfect circle or a wave).

3. The "Volume Knob" (The Order Parameter)

The most important part of their discovery is a number they call the largest eigenvalue.

• Think of this number as a volume knob for the "Critical Character."
• If the knob is turned down (low number), the system is chaotic and disordered (just random noise).
• If the knob is turned up (high number), it means the "Critical Character" has taken over. The system has become ordered, and a specific, large-scale pattern has formed.
• The authors found that as they added more "critical signal" to their computer simulations, this volume knob turned up, and the patterns became clearer and more organized (like the crowd forming distinct patches or clusters).

4. Why This is a Game-Changer

The paper highlights four main advantages of this new method:

1. It ignores the noise: It doesn't need to subtract out the "background noise" manually. It naturally separates the interesting pattern from the random chaos.
2. It doesn't need "perfect" conditions: Traditional methods assume the system settles down and becomes calm (thermal equilibrium) before you can measure it. This new method works even if the system is still chaotic and evolving rapidly (which is exactly what happens in these collisions).
3. It finds the "hidden" signal: It can detect the critical pattern even when it's mixed with a lot of non-critical data.
4. It's efficient: You don't need billions of events to see the result; a few thousand are enough to see the pattern emerge.

The Bottom Line

The authors tested this on computer simulations (mixing "normal" collision data with "critical" data). They found that their method successfully spotted the critical patterns, identified them as distinct "shapes" (like patches or rings), and measured their strength.

They conclude that this tool is ready to be applied to real data from the RHIC Beam Energy Scan (a major experiment at Brookhaven National Laboratory). It offers a fresh, powerful way to hunt for the elusive "Critical Point" in the universe's building blocks without getting lost in the noise.

Technical Summary

Technical Summary: Eigen-microstate Signatures of Criticality in Relativistic Heavy-Ion Collisions

Problem Statement
Locating the QCD critical point (CP) and mapping the phase boundary between hadronic matter and quark-gluon plasma (QGP) are central goals of relativistic heavy-ion programs (e.g., RHIC BES, CBM, NICA). While theoretical models predict a first-order phase transition at high baryon density terminating at a second-order CP, experimental identification remains elusive. Conventional observables, such as higher-order moments of conserved charges or factorial moments of multiplicity distributions, face significant challenges: they require extremely high statistics, necessitate careful subtraction of non-critical backgrounds (e.g., statistical fluctuations, resonance decays, conservation laws), and rely on assumptions of thermal equilibrium. Furthermore, the degree of equilibration in heavy-ion collisions is intrinsically uncertain, and the relevant order parameter may not be directly accessible through final-state observables.

Methodology
The authors apply the Eigen-Microstate Approach (EMA), a framework originally developed for systems without equilibrium assumptions, to relativistic heavy-ion collisions. The methodology proceeds as follows:

1. Original Microstate (OM) Construction: Each collision event is treated as a microstate defined by fluctuations in the final-state charged particle multiplicity within transverse-momentum ($p_T$) space. The $p_T$ plane is partitioned into an $L \times L$ lattice (restricted to $p_T \in [-2.0, 2.0]$ GeV/c). For the $I$-th event, the fluctuation in the $i$-th cell is defined as $\Delta N_{ch,i}^I = N_{ch,i}^I - \langle N_{ch,i} \rangle$, where $\langle N_{ch,i} \rangle$ is the event-averaged multiplicity. These fluctuations form a normalized vector $A^I$.
2. Correlation Matrix and Eigen-decomposition: A temporal-spatial correlation matrix $C_{IJ} = [A^I]^T \cdot A^J$ is constructed from an ensemble of $M$ events. Solving the eigenvalue equation $C \mathbf{b}^I = \lambda_I \mathbf{b}^I$ yields ordered eigenvalues $\lambda_1 \ge \lambda_2 \ge \dots \ge \lambda_M$ and corresponding eigenvectors.
3. Eigen-Microstates (EMs) and Order Parameter: The eigenvectors define the eigen-microstates $E^I$. The square of the amplitude of an EM corresponds to its eigenvalue ($|E^I|^2 = \lambda_I$), which serves as a statistical weight $w_I$. The largest weight, $w_1$, is proposed as a robust order parameter. A finite $w_1$ in the limit $M \to \infty$ indicates the "condensation" of a dominant collective mode, analogous to Bose-Einstein condensation or magnetization in the Ising model.

Simulation Setup
To validate the method, the authors generated event samples using two models:

• UrQMD (Ultra-relativistic Quantum Molecular Dynamics): Simulates non-critical heavy-ion collisions (0–5% central Au+Au at $\sqrt{s_{NN}} = 19.6$ GeV), incorporating collective flow and standard dynamics.
• CMC (Critical Monte Carlo): Generates momentum-space fluctuations with fractal criticality via Lévy random walks.
• Hybrid Samples: Critical signals were introduced into UrQMD events by replacing a fraction $\alpha_p$ of particles with CMC particles, subject to kinematic cuts ($|\eta| < 0.5$, specific $p_T$ windows) and momentum conservation constraints.

Key Results

• Extraction of Critical Modes: In the baseline UrQMD sample, the first eigen-microstate (EM) exhibits a ring-like structure with uniform positive fluctuations, while higher EMs show random, rotationally symmetric distributions. In hybrid samples with increasing $\alpha_p$, the EM patterns evolve. At $\alpha_p = 6\%$, nascent two-patch structures appear in the second and third EMs. As $\alpha_p$ increases to 12%, the first two EMs display pronounced two-patch patterns, and the third shows a four-patch pattern. These multi-patch structures reflect enhanced correlation lengths and critical clustering.
• Order Parameter Behavior: The weight $w_1$ (and $w_2$) increases with the critical signal fraction $\alpha_p$, peaking near $\alpha_p \approx 70\%$ before saturating. The weight cumulant $c(m) = \sum_{I=1}^m w_I$ rises more rapidly and saturates earlier for samples with higher $\alpha_p$, indicating that fewer leading EMs dominate the ensemble, a signature of an ordered phase.
• Fractal and Scaling Properties: The patch structures in the EMs remain qualitatively unchanged across different lattice resolutions ($L=10, 40, 100$), confirming the self-similar fractal nature of the critical mode. Furthermore, the largest eigenvalue $w_1$ exhibits a linear dependence on $L$ in a double-log scale, and the ratios $w_2/w_1$ and $w_3/w_1$ approach fixed-point values for $\alpha_p > 9\%$, confirming finite-size scaling (FSS) behavior.

Significance and Claims
The paper claims that the EMA offers a powerful new tool for the search for the QCD critical point with four distinct advantages over conventional methods:

1. Background Insensitivity: It extracts collective patterns directly from the event ensemble without requiring explicit subtraction of non-critical backgrounds (e.g., resonance decays, statistical fluctuations).
2. No Equilibrium Assumption: The framework is formulated for arbitrary event ensembles, making it applicable even when thermal equilibration is incomplete, a common feature of relativistic heavy-ion collisions.
3. Effective Order Parameter: It provides a practical indicator ($w_1$) of critical mode strength when the true QCD order parameter is experimentally inaccessible.
4. Moderate Statistical Requirements: Stable extraction is achievable with $\sim 20,000$ events, which is computationally efficient compared to higher-order cumulant analyses.

The authors conclude that while further work is needed to assess sensitivity under realistic detector effects, the EMA provides a complementary perspective for ongoing searches in RHIC BES-II and future heavy-ion programs, offering a systematic way to analyze critical patterns in non-equilibrium systems.