⚛️ nuclear experiments

Locating the QCD critical point with finite-size scaling of proton cumulants

This paper analyzes net-proton cumulants in Au+Au collisions using finite-size scaling to identify a potential QCD critical point near $\mu_B \approx 700$ MeV, but concludes that while the observed scaling patterns are consistent with criticality, they are not yet definitive evidence due to similar behaviors arising from non-critical dynamical simulations and experimental acceptance effects.

Original authors: Agnieszka Sorensen, Paul Sorensen

Published 2026-09-01

📖 5 min read🧠 Deep dive

CC BY 4.0

Original authors: Agnieszka Sorensen, Paul Sorensen

Original paper licensed under CC BY 4.0 (http://creativecommons.org/licenses/by/4.0/). ✨ This is an AI-generated explanation of the paper below. It is not written or endorsed by the authors. For technical accuracy, refer to the original paper. Read full disclaimer

The universe as we know it is built from matter that is far more complex than the solid objects we touch every day. Deep inside the cores of neutron stars and in the first fractions of a second after the Big Bang, this matter exists in a state where the fundamental building blocks of protons and neutrons melt together into a hot, dense soup called the quark-gluon plasma. Physicists have long suspected that as this soup cools and changes back into ordinary matter, it does not always do so smoothly. Instead, under certain conditions of extreme density and temperature, it might undergo a dramatic shift, similar to how water suddenly turns to ice, but with a twist: there may be a specific point where this change becomes a sharp, explosive transition rather than a gradual one. Finding this specific point, known as the critical point, is one of the great quests in modern physics because it would reveal the fundamental rules governing how matter behaves under the most extreme conditions imaginable.

To hunt for this elusive point, scientists smash heavy gold atoms together at nearly the speed of light in massive particle accelerators. By adjusting the energy of these collisions, they can recreate different temperatures and densities, effectively scanning a map of the possible states of matter. In a recent study, researchers Agnieszka Sorensen and Paul Sorensen analyzed data from these high-energy collisions to see if the patterns of particles produced could reveal the location of this critical point. They focused on the fluctuations in the number of protons created during the collisions, looking for a specific mathematical signature that would indicate the system was approaching a critical threshold. Their approach was clever: instead of looking at the entire collision event at once, they examined smaller slices of the collision zone, treating these slices as if they were smaller versions of the whole system. By seeing how the fluctuations changed as they looked at larger or smaller slices, they hoped to spot the tell-tale signs of a critical point.

The researchers examined data from gold-gold collisions at a wide range of energies, from very low speeds up to the highest energies achieved at the Relativistic Heavy Ion Collider. They measured how the number of protons varied from one collision to the next, a quantity that acts like a thermometer for the system's internal pressure and density. When they looked at collisions with energies of 7.7 billion electron volts and higher, they found that the fluctuations behaved in a very specific way. As they changed the size of the slice they were observing, the data seemed to collapse onto a single, smooth curve, a pattern that theoretical models predict should happen if a critical point is nearby. Furthermore, when they plotted these fluctuations against the density of the matter, the data suggested a critical point located at a baryon chemical potential of roughly 700 MeV, a value that aligns with several recent theoretical predictions.

However, the story does not end with a simple discovery. The researchers were careful to test whether these patterns could be created by something other than a critical point. They ran computer simulations of collisions that deliberately excluded any critical behavior, using standard models of how particles interact. Surprisingly, these simulations, which contained no critical point at all, produced patterns that looked remarkably similar to the real experimental data. The simulations showed that the same mathematical scaling could emerge simply from the way particles are distributed and conserved during the collision, without any need for a phase transition. This finding casts a shadow of doubt over the initial excitement, suggesting that the patterns seen in the real data might be a coincidence of the experimental setup rather than a definitive signal of new physics.

The investigation also looked at a different mathematical tool called the Binder cumulant, which is designed to show a specific crossing point if a critical state is reached. The experimental data did show such a crossing, which initially appeared to be strong evidence. Yet, when the researchers applied the same analysis to their simulations, they found that the crossing could appear simply because the detectors used in the experiment view the collision from slightly different angles at different energies. This sensitivity to the experimental "window" means that the crossing might be an artifact of how the data was collected rather than a fundamental property of the matter itself.

Ultimately, the study concludes that while the experimental data is consistent with the existence of a critical point, it cannot yet be taken as proof. The patterns observed are compelling, but they are not unique; they can be mimicked by standard physics models that do not include a critical point. The researchers emphasize that the apparent signals are likely a mix of genuine critical behavior and background effects that are difficult to separate. Until these background effects can be fully understood and ruled out, the search for the QCD critical point remains an open and challenging chapter in our understanding of the universe. The work serves as a crucial reminder that in the quest to map the extremes of nature, seeing a pattern is only the first step; proving that the pattern is real requires ruling out every other possible explanation.

Problem Statement
The location of the Quantum Chromodynamics (QCD) critical point (CP), where the crossover transition between hadronic matter and quark-gluon plasma (QGP) becomes a first-order phase transition, remains an open question in nuclear physics. While lattice QCD calculations at zero baryon chemical potential ( $\mu_B = 0$ ) confirm a smooth crossover, they are hindered by the fermion sign problem at finite $\mu_B$ , preventing direct calculation of the phase diagram in this region. Experimental searches, such as the Beam Energy Scan (BES) at the Relativistic Heavy Ion Collider (RHIC), utilize heavy-ion collisions to vary temperature ( $T$ ) and $\mu_B$ . A primary observable for this search is the net-proton number cumulants, which are expected to exhibit non-monotonic behavior or specific scaling laws near the CP. However, distinguishing critical fluctuations from non-critical background effects (such as conservation laws, thermal motion, and acceptance effects) remains a significant challenge.

Methodology
The authors perform a finite-size scaling (FSS) analysis on net-proton number cumulants extracted from 0–5% central Au+Au collisions at center-of-mass energies ( $\sqrt{s_{NN}}$ ) ranging from 2.4 to 200 GeV.

System Size Proxy: Instead of varying collision centrality (which alters global system properties like temperature and volume), the authors vary the system size by analyzing subvolumes defined by the rapidity bin width ( $W$ ).
Observables: They extract the second-order susceptibility ( $\chi_2$ $χ_{2}$ ) and the Binder cumulant ( $U_4$ $U_{4}$ ).
- $\chi_2$ is derived from the second-order net-proton cumulant ( $C_2$ ), normalized by freeze-out temperature ( $T_{fo}$ ) and volume density.
- $U_4$ is calculated as $1 - \langle(\delta N)^4\rangle / (3\langle(\delta N)^2\rangle^2)$ .
Scaling Analysis: The study tests for data collapse onto universal curves by plotting scaled variables (e.g., $\chi_2 W^{-\gamma/\nu}$ vs. $t W^{1/\nu}$ , where $t$ is the reduced distance to the CP) using 3-D Ising critical exponents. The critical baryon chemical potential ( $\mu_{B,c}$ ) is tuned to achieve this collapse.
Control Simulations: To assess whether observed scaling is unique to critical physics, the authors compare experimental data with SMASH hadronic transport simulations run without a CP or mean fields. They also analytically test non-critical distributions (Poisson and Skellam) along the chemical freeze-out line.

Key Results

Susceptibility Scaling: For $\sqrt{s_{NN}} \geq 7.7$ $s_{N N} \geq 7.7$ GeV, the susceptibility $\chi_2$ $χ_{2}$ shows a weak dependence on the rapidity bin width $W$ $W$ but increases as a power law with respect to the distance from a critical point: $\chi_2 \sim (\mu_B - \mu_{B,c})^{-\gamma}$ $χ_{2} \sim (μ_{B} - μ_{B, c})^{- γ}$ .
- The fit yields $\mu_{B,c} = 700 \pm 130$ MeV and $\gamma = 1.6 \pm 0.5$ , consistent with 3-D Ising universality class estimates.
- When $\gamma$ is fixed to the 3-D Ising value ($1.237 $),$ \mu_{B,c}$ shifts to $605 \pm 11$ MeV.
- Crucial Counter-evidence: SMASH simulations without a CP reproduce a similar power-law divergence (at $\mu_B \approx 558$ MeV). The authors demonstrate that for a Skellam distribution (which describes BES net-proton multiplicities), the interplay between the freeze-out line and the ratio of cumulants ( $C_2/C_1$ ) can naturally generate an apparent power-law scaling that mimics critical behavior.
Binder Cumulant Crossing: The experimental data for $U_4$ $U_{4}$ at different $W$ $W$ values appear to cross at $\mu_{B,c} \approx 807$ $μ_{B, c} \approx 807$ MeV.
- However, SMASH simulations using asymmetric rapidity bins (mimicking the experimental setup at $\sqrt{s_{NN}} = 3.0$ GeV) also produce an apparent crossing, whereas symmetric bins do not.
- This suggests the crossing is sensitive to experimental acceptance and binning geometry rather than being a definitive signature of criticality.
Higher-Order Cumulants: Analysis of the third-order susceptibility ( $\chi_3$ ) yields a significantly higher $\mu_{B,c}$ ($811$ MeV) and exhibits poorer scaling behavior compared to $\chi_2$ , further complicating the interpretation.

Significance and Claims
The paper concludes that while the observed finite-size scaling of $\chi_2$ and the crossing of the Binder cumulant are consistent with the presence of a QCD critical point near $\mu_B \approx 700$ MeV, they cannot yet be taken as definitive evidence.

The authors emphasize that "non-critical baselines" containing standard QCD physics (specifically the chemical freeze-out line and Skellam-like multiplicity distributions) can reproduce the observed scaling and crossing signatures. Consequently, the apparent scaling of $\chi_2$ is not unique to critical fluctuations. The study highlights the necessity of understanding how non-critical ingredients of the scaling function evolve across the phase diagram before claiming the discovery of the critical point. The work serves as a cautionary analysis, demonstrating that signatures traditionally associated with criticality can emerge from non-critical dynamics under specific experimental conditions.

More like this