⚛️ nuclear experiments

Temperature fluctuations in a realistic Polyakov-loop extended Nambu--Jona-Lasinio Model along the freeze-out line

By reparameterizing the Polyakov-loop extended Nambu--Jona-Lasinio model to align with lattice data and the BES-II energy range, this study demonstrates that a distinct dip in the second-order temperature cumulant along the freeze-out line near 7.7 GeV supports the hypothesis that the non-monotonic transverse momentum correlations observed by the STAR Collaboration are linked to the QCD critical endpoint.

Original authors: He Liu, Peng Wu, Hong-Ming Liu, Peng-Cheng Chu

Published 2026-07-22✓ Author reviewed ⓘ

📖 5 min read🧠 Deep dive

CC BY 4.0

Original authors: He Liu, Peng Wu, Hong-Ming Liu, Peng-Cheng Chu

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 by the authors. For technical accuracy, refer to the original paper. Read full disclaimer

Imagine the universe as a giant, cosmic soup. In the very first moments after the Big Bang, this soup was so hot and dense that the tiny building blocks of matter—protons and neutrons—couldn't stick together. Instead, they melted into a free-flowing, chaotic sea of even smaller particles called quarks and gluons, known as the Quark-Gluon Plasma (QGP). As the universe cooled, this soup "froze" into the solid matter we see today, like the atoms in your body. But what if this freezing process wasn't just a simple, smooth transition? What if, under the right conditions of heat and pressure, the universe hit a "critical point"—a specific spot on the map where the rules of matter change dramatically, like water turning into ice but with a sudden, wild burst of activity?

Scientists are obsessed with finding this "Critical Endpoint" (CEP). To do this, they smash heavy atoms together at nearly the speed of light in massive machines called colliders. These collisions recreate the hot, dense conditions of the early universe for a split second. By changing the energy of the smash, scientists can explore different parts of the "phase diagram," which is basically a map showing how matter behaves under different temperatures and pressures. The big question is: Does this map have a special spot where the transition from a hot soup to solid matter becomes chaotic and unpredictable? If we find it, it would be a massive discovery about how the universe works.

The Cosmic Treasure Hunt

In this study, a team of physicists led by He Liu and Peng Wu decided to play detective with a theoretical map to see if they could spot the clues to this Critical Endpoint. They used a sophisticated mathematical tool called the "Polyakov-loop extended Nambu–Jona-Lasinio" (PNJL) model. Think of this model as a highly detailed video game engine that simulates how the cosmic soup behaves. However, previous versions of this engine had some glitches; they didn't quite match the data from real-world experiments or the supercomputer simulations (called Lattice QCD) that act as the "gold standard" for these calculations.

So, the authors upgraded their engine. They tweaked the rules of the game—adding new interactions between the particles and adjusting the settings—until their simulation perfectly matched the known data at zero pressure and placed their predicted "Critical Endpoint" right in the zone where the STAR experiment at the Relativistic Heavy Ion Collider (RHIC) is currently looking. This new, "realistic" model is their trusty compass for the journey.

The Dip in the Data

With their upgraded model ready, the team set out to investigate a specific mystery. The STAR experiment had recently noticed something weird in their data: as they changed the energy of their collisions, the way particles moved sideways (their transverse momentum) didn't just go up or down smoothly. Instead, it showed a "dip"—a sudden drop and then a rise again—around a collision energy of 7.7 GeV. It was like driving a car and feeling the road suddenly get bumpy and then smooth out again at a specific mile marker.

The researchers wanted to know: Is this dip a sign that they are passing near the Critical Endpoint? To find out, they calculated the "temperature fluctuations" along the path where the cosmic soup freezes (the "freeze-out line"). In their simulations, they tracked a number called $C_2$ , which measures how much the temperature jitters.

The results were exciting. Their model showed a distinct "dip" in the temperature fluctuations right near the Critical Endpoint and the boundary where the matter changes phase. When they traced this path using the experimental data from the STAR experiment, the dip in their simulation appeared right around the same energy: 7.7 GeV. This suggests that the strange dip STAR saw might indeed be a fingerprint of the Critical Endpoint.

Ruling Out the Noise

However, the authors were careful not to jump to conclusions. They acknowledged that other scientists have suggested this dip might just be a trick of the math, caused by mixing different types of particles (like protons and pions) rather than a sign of a critical point. Their study doesn't completely disprove that idea, but it strongly suggests that the Critical Endpoint is a very plausible explanation. They also noted that their model is a "mean-field" simulation, which means it ignores some tiny, chaotic quantum jitters that happen in real life. Because of this, while their model gets the shape of the dip and the location (7.7 GeV) right, the exact size of the dip might be different in reality.

Looking Deeper with Ratios

To get a clearer picture and filter out the "noise" of the experiment (like the fact that every collision is slightly different in size), the team looked at more complex numbers called "cumulant ratios," such as $C_3/C_2^2$ and $C_4/C_3^2$ . They found that these ratios showed even more dramatic, non-monotonic patterns—sharp peaks and dips—that are very sensitive to the Critical Endpoint. They argue that measuring these specific ratios in future experiments would be a much better way to confirm if the Critical Endpoint is real, as these numbers are less likely to be fooled by initial volume fluctuations.

The Road Ahead

In short, this paper suggests that the mysterious dip in particle movement observed by the STAR experiment could be the long-sought signal of the Critical Endpoint. The authors' upgraded model predicts this dip at the exact energy where it was seen. While this is a strong hint, it's not a final proof. The authors suggest that future experiments need to measure even higher-order fluctuations and that scientists should use advanced simulations that include critical dynamics to confirm these findings. The treasure map looks promising, but the final treasure hunt is just beginning.

Technical Summary: Temperature Fluctuations in a Realistic Polyakov-loop Extended Nambu–Jona-Lasinio Model along the Freeze-out Line

Problem Statement
The search for the QCD Critical Endpoint (CEP) and the first-order phase transition boundary is a primary objective of relativistic heavy-ion collision experiments, particularly the Beam Energy Scan (BES) program at RHIC. Recent high-statistics data from the STAR Collaboration have revealed a significant non-monotonic dip structure in the energy dependence of the two-particle transverse momentum correlation, $C_{pT}$ , in central collisions at $\sqrt{s_{NN}} \approx 7.7$ GeV. While this behavior is hypothesized to be a signal of critical fluctuations near the CEP, alternative interpretations suggest it may arise from the mixing of particle distributions (e.g., protons and pions) during the transition from a baryon-dominated to a meson-dominated system. Furthermore, previous theoretical studies using the 2+1 flavor LEFT-fRG framework failed to reproduce a non-monotonic dip in the second-order temperature cumulant ( $C_2$ ) along the freeze-out line. There is a need for a theoretical framework that simultaneously reproduces lattice QCD data at zero baryon chemical potential ( $\mu_B$ ) and positions the CEP within the experimental energy range to clarify the origin of these fluctuations.

Methodology
The authors extend and reparameterize the three-flavor Polyakov-loop extended Nambu–Jona-Lasinio (PNJL) model to create a "realistic" model capable of addressing these issues. The methodology involves:

Model Extension: The standard PNJL Lagrangian is augmented with two eight-quark interaction terms ( $\mathcal{L}_{8q}^1$ and $\mathcal{L}_{8q}^2$ ) to resolve vacuum instabilities and lower the crossover temperature ( $T_c$ ) at $\mu_B=0$ , ensuring agreement with Lattice QCD results (HotQCD and Wuppertal-Budapest collaborations).
CEP Tuning: A scalar-vector coupling interaction ( $\mathcal{L}_{SV}$ ) is introduced. The coupling constant $G_{SV}$ is treated as a free parameter (specifically set to negative values) to shift the position of the CEP to higher temperatures and lower baryon chemical potentials, aligning it with the BES-II experimental energy range (3–7.7 GeV).
Thermodynamic Framework: Temperature fluctuations are analyzed using a systematic framework where the $n$ -th order temperature cumulant $C_n$ is derived from the $n$ -th derivative of the thermodynamic potential. The second-order cumulant $C_2$ is identified as the inverse of the dimensionless heat capacity ( $C_V/T^3$ ).
Freeze-out Trajectory: Calculations are performed along experimentally extracted chemical freeze-out lines (based on fits from Lu et al. and Andronic et al.) covering the $\mu_B$ range of approximately 400–760 MeV.
Analysis: The study computes the behavior of $C_2$ and higher-order cumulants ( $C_3$ through $C_6$ ) and their dimensionless ratios (e.g., $C_3/C_2^2$ , $C_4/C_3^2$ ) along these freeze-out lines to investigate their connection to the phase diagram structure.

Key Contributions and Results

Realistic Phase Diagram: The reparameterized model successfully reproduces the lattice QCD crossover temperature ( $T_c \approx 156.6$ MeV) at $\mu_B=0$ while placing the CEP in the region probed by BES-II. The chiral crossover transition is found to be largely insensitive to the $G_{SV}$ parameter, whereas the CEP position shifts significantly with variations in $G_{SV}$ .
Dip Structure in $C_2$ : The study reveals a distinct dip structure in the second-order temperature cumulant $C_2$ within the chiral phase transition region, particularly near the CEP and the first-order phase boundary. This dip arises because the entropy discontinuity at the first-order transition induces a peak in heat capacity, suppressing $C_2$ .
Correlation with STAR Data: Along the experimentally extracted freeze-out lines, the model predicts a minimum in $C_2$ occurring around $\sqrt{s_{NN}} \approx 7.7$ GeV. This trend is consistent with the non-monotonic dip observed by the STAR Collaboration in $C_{pT}$ measurements. Since $C_{pT}$ is proportional to $\sqrt{C_2}$ , the authors suggest the observed dip may indeed be related to the CEP.
Higher-Order Cumulants and Ratios: While $C_3$ shows no significant non-monotonic behavior, higher-order cumulants ( $C_4, C_5, C_6$ ) exhibit complex non-monotonic structures. Crucially, the dimensionless cumulant ratios (e.g., $C_3/C_2^2$ , $C_4/C_3^2$ ) display pronounced peaks or dips. The authors argue that these ratios are superior indicators of critical fluctuations because they eliminate the influence of initial volume fluctuations, which are significant in low-energy collisions but not accounted for in the mean-field calculation of raw cumulants.

Significance and Claims
The paper claims to provide a theoretical basis linking the non-monotonic behavior of transverse momentum correlations to the QCD CEP. By employing a realistic PNJL model that satisfies both lattice constraints and experimental freeze-out conditions, the authors demonstrate that a dip in temperature fluctuations naturally occurs near the CEP.

The authors maintain a modest stance regarding their findings:

They acknowledge that their model relies on equilibrium thermodynamics and the mean-field approximation, neglecting quantum fluctuations and initial volume effects. Consequently, while the trend and minimum position of the dip agree with experimental data, the numerical magnitudes deviate.
They do not claim definitive proof of the CEP's existence but suggest that the observed dip in $C_{pT}$ is consistent with CEP-induced critical fluctuations.
They propose that future experimental measurements of higher-order event-mean transverse momentum fluctuations, particularly the dimensionless cumulant ratios, could provide further verification. They suggest that hydrodynamic or transport simulations incorporating critical dynamics would be necessary to fully validate these predictions.

In summary, the work offers a priori theoretical predictions for future experimental measurements, supporting the hypothesis that the non-monotonic energy dependence of $C_{pT}$ is a signature of the QCD phase structure near the critical endpoint.

More like this