Model Study of Eigen-Microstate Signatures of Criticality in Relativistic Heavy-Ion Collisions

This study demonstrates that the eigen-microstate approach (EMA) serves as a robust, background-independent method for identifying critical fluctuations in relativistic heavy-ion collisions by effectively filtering non-critical correlations and capturing the fractal nature of criticality through characteristic eigenvalue patterns and finite-size scaling behaviors.

The Gist

The Big Picture: Finding a Needle in a Haystack

Imagine physicists are trying to find a specific, rare type of weather pattern (a "critical point") inside a massive, chaotic storm (a heavy-ion collision). The problem is that the storm is full of normal wind, rain, and thunder (background noise) that looks very similar to the rare pattern they are hunting for.

For decades, scientists have tried to spot this rare pattern by measuring specific things, like how loud the thunder is or how fast the wind blows. But the paper argues that these methods get confused by all the normal noise.

Instead, the authors propose a new detective tool called Eigen-Microstate Approach (EMA). Think of this not as measuring the wind speed, but as looking at the entire shape of the storm cloud to see if it has a hidden, repeating structure.

How the New Tool Works: The "Group Photo" Analogy

Imagine you take 20,000 photos of a crowd of people at a concert.

• The Old Way: You count how many people are wearing red shirts, or how many are jumping. This is like looking at individual particles.
• The New Way (EMA): You lay all 20,000 photos on a table and ask a super-smart computer to find the "common themes" that connect them.

The computer breaks the crowd down into "modes" or "patterns":

1. The Main Pattern (The "Condensation"): If everyone is just standing there, the main pattern is just "a crowd."
2. The Critical Pattern: If a secret group of people starts dancing in a specific, synchronized, fractal way (like a fractal snowflake), the computer spots this as a distinct, dominant shape that stands out from the noise.

The paper claims that if a "critical point" exists in the collision, it will create a specific, organized "dance" (a cluster-like pattern) that the computer can see clearly, even if it's mixed with millions of other random movements.

The Experiments: Testing the Tool

The authors tested this tool using three different "crowds" (simulations):

1. The "Normal" Crowds (UrQMD and Stochastic Models)
They simulated heavy-ion collisions that don't have a critical point.

• The Result: The computer looked at the data and said, "I see a crowd, and I see some random noise." It found no organized "dance."
• The Lesson: The tool is very good at ignoring normal physics (like particles bouncing off each other or conservation laws). It filters out the background noise so it doesn't get fooled.

2. The "Fake" Critical Crowds (Hybrid Models)
They took the normal simulations and secretly swapped in some "critical" events (using a model called CMC that mimics the fractal nature of a critical point). They did this in two ways:

• Scenario A (Event-Level): They replaced entire photos of the crowd with photos of the "dancing" group.
  • Result: The computer spotted the dance immediately, even if only 1% of the photos were swapped.
• Scenario B (Particle-Level): They took a normal photo and swapped out just a few people in the crowd with "dancers."
  • Result: The computer needed to swap out a much larger percentage (around 9-12%) before it could clearly see the dance pattern.

The Takeaway: The tool is much better at spotting a "critical point" if the whole event is critical, rather than just a few particles. However, it can still find the signal even if it's hidden in a small fraction of the data.

The "Magic Number" and the "Fixed Point"

The paper introduces two key ways to confirm they found the real thing:

1. The "Leader" (The Largest Eigenvalue):
   Think of the computer finding a "leader" of the patterns. In a normal crowd, the leader is weak. But when the critical "dance" appears, this leader suddenly becomes very strong and dominant. The paper suggests this "strength" acts like a thermometer: as you get closer to the critical point, this number goes up and stabilizes.

2. The "Zoom Test" (Finite-Size Scaling):
   Imagine looking at the "dance" pattern through a microscope.

   • If you zoom in (look at a small area) or zoom out (look at the whole room), does the pattern look the same?
   • Real critical phenomena are fractal, meaning they look the same at every scale (like a fern leaf or a coastline).
   • The authors tested their tool at different "zoom levels" (different grid sizes). They found that when the critical signal is strong, the ratio of the "second strongest pattern" to the "strongest pattern" becomes the same regardless of the zoom level. This "fixed point" is a strong fingerprint that the signal is genuine criticality, not just random noise.

Summary

This paper is a "model study," meaning they tested their new method on computer simulations, not real experimental data yet.

They concluded that:

• The Eigen-Microstate Approach is a robust way to find critical fluctuations.
• It successfully filters out the "noise" of normal particle collisions.
• It can detect critical signals even when they are a tiny fraction of the total data.
• It identifies the critical point by looking for organized, fractal patterns and a dominant "leader" pattern that behaves consistently no matter how you scale the data.

The authors suggest that this method should be used on real data from the RHIC (Relativistic Heavy Ion Collider) and future experiments to finally locate the elusive QCD critical point.

Technical Summary

Technical Summary: Model Study of Eigen-Microstate Signatures of Criticality in Relativistic Heavy-Ion Collisions

Problem Statement
The identification of the Quantum Chromodynamics (QCD) critical point (CP) remains a primary objective in contemporary nuclear physics. While theoretical models predict a first-order phase transition at high baryon density terminating at a second-order CP, experimental confirmation via the Relativistic Heavy Ion Collider (RHIC) Beam Energy Scan (BES) and future facilities (FAIR, NICA, HIRFL-CSR) has been elusive. A significant challenge is that conventional observables (e.g., higher-order cumulants, factorial moments) are heavily contaminated by non-critical backgrounds, including resonance decays, global conservation laws, hadronic rescattering, and collective flow. Furthermore, the short-lived, non-equilibrium nature of heavy-ion collisions complicates the extraction of genuine critical fluctuations from these large backgrounds. There is a need for methods that can isolate critical modes without relying heavily on equilibrium assumptions or explicit background subtraction.

Methodology
The authors employ the Eigen-Microstate Approach (EMA), a framework that treats each collision event as a microstate and analyzes the ensemble via the diagonalization of an event-event correlation matrix. The study utilizes the following models and procedures:

• Baseline Models:
  • UrQMD: A realistic cascade model for non-critical heavy-ion collisions (Au+Au at $\sqrt{s_{NN}} = 19.6$ GeV), incorporating standard dynamics like conservation laws and collective flow.
  • Stochastic Models: Two reference ensembles constructed to share UrQMD's multiplicity distribution but differ in kinetic constraints. Stochastic Model I has no kinetic constraints (uniform momentum assignment), while Stochastic Model II imposes global kinetic constraints (Gamma-distributed $p_T$ and elliptic flow modulation) but lacks intrinsic particle-particle correlations.
• Critical Embedding:
  • Critical Monte Carlo (CMC): A model generating scale-invariant fractal fluctuations via Lévy random walks in momentum space, characteristic of systems near a CP.
  • Hybrid Construction: Critical fluctuations are embedded into the UrQMD baseline via two scenarios:
    1. Event-level replacement: A fraction $\alpha_e$ of UrQMD events is replaced entirely by CMC events.
    2. Particle-level replacement: A fraction $\alpha_p$ of particles within each UrQMD event is replaced by particles from a CMC event, subject to a $p_T$-matching condition to preserve the baseline spectrum.
• Analysis: The correlation matrix $C_{ij}$ is constructed from transverse momentum fluctuations. Diagonalization yields eigenvalues ($\lambda_n$) and eigenvectors, which define the Eigen-Microstates (EMs). The study analyzes the spatial patterns of the leading EMs, the hierarchy of eigenvalue weights ($w_n = \lambda_n$), cumulative weights $C(m)$, and finite-size scaling behaviors across different binning scales ($L$).

Key Contributions and Results

1. Filtering of Non-Critical Backgrounds:

   • In non-critical ensembles (UrQMD and Stochastic Model II), the leading eigen-microstate (EM1) represents an average global mode reflecting the overall radial momentum distribution, while EM2 and EM3 exhibit random, noise-like fluctuations without stable cluster structures.
   • The cumulative weight $C(m)$ for UrQMD and Stochastic Model II shows rapid initial saturation, indicating that the eigenvalue spectrum is dominated by global kinetic constraints (average $p_T$ and $\phi$ distributions) rather than specific dynamical correlations.
   • Crucially, EMA is shown to be largely insensitive to conventional short-range correlations (resonance decays, rescattering) present in UrQMD, effectively filtering them out to leave only the average global shape.

2. Emergence of Critical Signatures in Hybrid Models:

   • Pattern Formation: As the critical fraction increases, the EMs develop distinct, coherent cluster-like (patch) structures that break the rotational symmetry of the baseline.
   • Sensitivity Comparison: EMA is significantly more sensitive to event-level criticality than particle-level criticality. A very small fraction of fully critical events ($\alpha_e \sim 0.1\% - 1\%$) is sufficient to reshape EM2 and EM3 and generate visible cluster patterns. In contrast, particle-level replacement requires much larger fractions ($\alpha_p \gtrsim 9\%$) to produce comparable cluster structures, reflecting the weaker coherence of localized critical droplets.
   • Order Parameter Behavior: The largest eigenvalue weight ($w_1$) behaves as an order-parameter-like quantity. It increases smoothly with the critical fraction and eventually saturates. The cumulative weight $C(m)$ saturates earlier in the presence of criticality, indicating that a small number of leading EMs dominate the ensemble.

3. Finite-Size Scaling and Fractal Nature:

   • Scale Invariance: The spatial patterns of the EMs exhibit self-similarity across different binning scales ($L = 10, 40, 100$), consistent with the fractal nature of critical fluctuations.
   • Fixed-Point Behavior: The ratio of eigenvalue weights ($w_3/w_1$) exhibits fixed-point behavior. At low critical fractions, the ratio depends on the system size $L$. However, as the critical fraction increases ($\alpha_e \gtrsim 4\%$, $\alpha_p \gtrsim 9\%$), the curves for different $L$ converge and overlap, confirming that EMA captures genuine critical scaling.

Significance and Claims
The paper claims that the Eigen-Microstate Approach offers a robust, background-independent framework for identifying critical fluctuations in relativistic heavy-ion collisions. By demonstrating that EMA effectively filters out conventional non-critical correlations while remaining highly sensitive to coherent critical modes, the study suggests that EMA can serve as a powerful tool for critical-point searches.

The authors posit that the methodology provides direct guidance for analyzing data from RHIC BES II and future experiments. Specifically, they propose that scanning collision energy and centrality for the onset of cluster-like EM patterns and fixed-point behavior in eigenvalue ratios constitutes a promising strategy for locating the QCD critical point. The study emphasizes that the largest eigenvalue acts as an effective order parameter, and the fractal nature of the EMs provides a distinct signature of criticality that is less susceptible to the systematic uncertainties plaguing traditional observables.