Quantifying fluctuation signatures of the QCD critical point using maximum entropy freeze-out

This paper utilizes a maximum entropy freeze-out approach to quantify how non-universal mapping parameters and the distance between the critical point and the freeze-out curve influence factorial cumulants of proton multiplicities, thereby linking QCD thermodynamics near a critical point to observable event-by-event fluctuations in heavy-ion collisions.

The Gist

The Big Picture: Hunting for a "Critical Point" in the Universe's Soup

Imagine the universe is a giant pot of soup. At the beginning of the Big Bang, this soup was a super-hot, chaotic liquid made of tiny, free-floating particles called quarks and gluons (this is called the Quark-Gluon Plasma or QGP). As the universe cooled down, this soup froze into solid chunks, forming protons and neutrons (the building blocks of atoms).

Physicists believe that if you heat this soup up again and squeeze it really hard (like in a particle collider), you can turn the solid chunks back into liquid. But there's a mystery: Is there a specific spot on the "temperature vs. pressure" map where this transition changes from a smooth melting (like ice to water) to a sudden, explosive boiling (like water to steam)?

That specific spot is called the Critical Point. Finding it is like finding the "Holy Grail" of nuclear physics. It would tell us exactly how matter behaves under extreme conditions.

The Problem: The Soup is Too Messy to Read

Scientists smash heavy atoms together at near-light speed to recreate this soup. They look for "signatures" of the Critical Point. One of the best signatures is fluctuation.

Think of it like this:

• Normal Soup: If you stir a pot of soup, the bubbles are random and small.
• Critical Soup: As you get close to the Critical Point, the soup starts to "stutter." Huge bubbles form and pop everywhere. The particles start acting like a giant, synchronized crowd rather than individuals.

Scientists measure these "bubbles" by counting how many protons come out of the collision. If the Critical Point exists, the number of protons will fluctuate wildly in a very specific pattern.

The Catch: The soup cools down and turns into particles so fast that the scientists can't see the soup itself; they only see the frozen chunks (the particles). It's like trying to figure out the weather inside a hurricane by only looking at the debris on the ground after the storm passes. You have to work backward to guess what the storm was doing.

The Solution: The "Maximum Entropy" Detective

This paper introduces a new, smarter way to work backward. The authors use a method called Maximum Entropy Freeze-Out.

The Analogy: The Perfectly Organized Moving Truck
Imagine you are packing a moving truck (the freeze-out). You have a chaotic pile of boxes (the fluid soup) and you need to pack them into a truck (the particles).

• Old Method: Scientists used to guess how the boxes were packed based on a few rules of thumb. They might say, "Well, heavy boxes go here," but they didn't know for sure how the wobbly boxes (fluctuations) were arranged.
• New Method (Maximum Entropy): This paper says, "Let's assume the most logical, least biased way to pack the truck." If you have a pile of energy and a pile of baryons (matter), and you want to turn them into particles without inventing any new rules, what is the most "random" (maximum entropy) way to do it?

This method acts like a mathematical translator. It takes the chaotic, wobbly fluctuations of the fluid just before it freezes and translates them into the specific patterns of particles just after it freezes. It ensures that the laws of physics (like conservation of energy) are respected in every single event, not just on average.

The Experiment: Testing Different Maps

The authors don't know exactly where the Critical Point is, so they created a "family" of possible maps (Equations of State). They imagined the Critical Point could be in different places and have different shapes.

They used four "knobs" (parameters) to tweak these maps:

1. Where is the point? (Location)
2. How strong is the wobble? (Strength)
3. What shape is the critical region? (Shape)
4. How is it tilted? (Orientation)

They ran their "Maximum Entropy Translator" on all these different maps to see what the proton fluctuations would look like for each scenario.

The Results: What the Patterns Look Like

The paper produces a series of graphs (Figs 6, 7, and 8) that show what scientists should look for in their data. Here is the "cheat sheet" they created:

1. The "Dip and Peak" Dance:
   If the Critical Point is nearby, the fluctuations won't just go up; they will do a specific dance.

   • The Dip: As you change the collision energy, the 4th-order fluctuation (a measure of how "spiky" the distribution is) might dip negative.
   • The Peak: Then, as you go further, it shoots up to a positive peak.
   • Analogy: It's like a rollercoaster that dips below the track before shooting up. If you see this specific "dip-then-peak" pattern in the proton counts, it's a very strong sign you've found the Critical Point.

2. The "Stretch" Effect:
   The paper shows that if the Critical Point is "wide" (a large critical region), the peak in the data will be wide and spread out. If it's "sharp," the peak will be narrow. This helps scientists figure out the shape of the Critical Point, not just its location.

3. The "Distance" Factor:
   The closer the collision happens to the Critical Point, the bigger the wobble. The paper calculates exactly how much the signal gets weaker if the collision happens a little further away.

Why This Matters

Before this paper, scientists had to make a lot of guesses about how the fluid turns into particles. They had to assume specific, often arbitrary, connections between the fluid and the particles.

This paper removes the guesswork. It says: "If the fluid is in equilibrium (calm) right before it freezes, here is the exact mathematical translation to the particles."

It provides a quantitative framework. It doesn't just say "Look for a wiggle." It says, "If the Critical Point is at location X with strength Y, you will see a wiggle of height Z."

The Caveat: The "Memory" Problem

The authors admit one big limitation: They assumed the fluid is calm (in equilibrium) right before it freezes. In reality, near a Critical Point, the fluid gets "lazy" (a phenomenon called critical slowing down). It takes longer to react to changes.

Think of it like a heavy truck trying to turn a corner. If you turn the wheel slowly, the truck follows. If you turn the wheel fast, the truck keeps going straight for a bit (it has "memory" of where it was).

• The paper assumes the truck turns perfectly with the wheel.
• In reality, the truck might lag.

However, the authors suggest a clever workaround: treat the "lag" as a variable. If the data shows a signal that looks like it's coming from a hotter temperature than the actual freeze-out, it means the fluid "remembered" the heat from earlier. This allows them to still use their method to find the Critical Point, even if the fluid isn't perfectly calm.

Summary

This paper is a user manual for finding the Critical Point.

1. It builds a translator (Maximum Entropy) that turns fluid wobbles into particle counts.
2. It tests this translator against many possible maps of the Critical Point.
3. It tells experimentalists exactly what pattern to look for (the dip-and-peak dance) and how the shape of that pattern tells them about the Critical Point's location and strength.

It's a crucial step in turning the "search" for the Critical Point into a precise "measurement" of the universe's most extreme states.

Technical Summary

1. Problem Statement

The existence and location of a critical point (CP) in the Quantum Chromodynamics (QCD) phase diagram, separating the hadronic phase from the quark-gluon plasma (QGP), remains one of the most significant open questions in nuclear physics. While lattice QCD calculations have confirmed a crossover at zero baryon chemical potential ($\mu_B = 0$), the nature of the transition at high $\mu_B$ is unknown.

• Experimental Challenge: Heavy-ion collision experiments (e.g., RHIC's Beam Energy Scan) search for the CP by measuring event-by-event fluctuations in particle multiplicities (specifically protons). Higher-order cumulants of these fluctuations are expected to diverge near the CP due to the divergence of the correlation length ($\xi$).
• Theoretical Gap: To interpret experimental data, a rigorous theoretical framework is needed to link the QCD Equation of State (EoS) near the CP to the observed particle spectra. Previous approaches often relied on phenomenological couplings between a critical scalar field ($\sigma$) and hadron masses, which lacked quantitative control over absolute magnitudes and ignored conservation laws. Furthermore, the non-universal mapping parameters that define the specific shape and strength of the critical region in QCD are unknown.

2. Methodology

The authors employ a two-pronged theoretical approach combining universality arguments with a novel statistical mechanics technique:

A. Equation of State (EoS) Construction via Ising Mapping

• Universality: The QCD critical point is assumed to belong to the 3D Ising universality class. The singular part of the QCD pressure ($P_{sing}$) is mapped to the Gibbs free energy of the 3D Ising model ($G_{Ising}$).
• Non-Universal Mapping: A parametric family of EoSs is constructed using four non-universal mapping parameters:
  • $\mu_c, T_c$: Location of the critical point.
  • $w, \rho$: Scaling parameters determining the size and shape of the critical region.
  • $\alpha_1, \alpha_2$: Angles defining the orientation of the Ising axes relative to the QCD phase diagram.
• Quadratic Mapping: The authors utilize an improved quadratic mapping (extending Ref. [72]) to ensure the Ising crossover line maps correctly to the curved chiral crossover in the QCD phase diagram.

B. Maximum Entropy Freeze-Out Procedure

• Core Concept: Instead of assuming phenomenological couplings, the authors use the Maximum Entropy (MaxEnt) method (developed in Ref. [36]) to translate hydrodynamic fluctuations (energy density $\epsilon$ and baryon density $n$) just before freeze-out into particle multiplicity fluctuations just after freeze-out.
• Principle: The method maximizes the relative entropy (minimizes information gain) subject to the constraints of local conservation laws (energy, momentum, baryon number) on an event-by-event basis.
• Result: This yields a direct, model-independent relation between the irreducible relative cumulants (IRCs) of hydrodynamic variables and the factorial cumulants of particle multiplicities.
• Assumption: The study assumes the system is in thermal equilibrium at freeze-out. This is a simplification, as critical slowing down (non-equilibrium effects) is expected to be significant near a CP. The authors treat the distance between the crossover and freeze-out curves ($\Delta T_f$) as a free parameter to account for potential "memory effects" where fluctuations lag behind equilibrium.

3. Key Contributions

1. Quantitative Framework: The paper provides the first quantitative application of the MaxEnt freeze-out procedure to QCD critical fluctuations, eliminating ad-hoc coupling parameters ($g_p$) used in previous studies.
2. Parameter Sensitivity Analysis: The authors systematically quantify how the non-universal mapping parameters ($w, \rho$) and the freeze-out distance ($\Delta T_f$) affect the magnitude, width, and position of peaks in proton factorial cumulants.
3. Scaling Relations: They derive explicit scaling laws connecting the cumulants of baryon density fluctuations to the mapping parameters, showing how the shape of the critical region in the phase diagram translates to observables.
4. Role of Energy Density: The analysis explicitly includes contributions from energy density fluctuations and mixed correlators, demonstrating that while baryon density dominates, energy density contributions are non-negligible for quantitative precision.

4. Key Results

The authors computed the second, third, and fourth factorial cumulants of proton multiplicity ($\Delta \omega_{2p}, \Delta \omega_{3p}, \Delta \omega_{4p}$) along freeze-out curves parallel to the chiral crossover.

• Dependence on Mapping Parameters ($w, \rho$):

  • $w$ (Strength/Scale): Increasing $w$ generally suppresses the peak heights of the cumulants (scaling roughly as $w^{-6/5}$) and squeezes the critical region toward the critical point.
  • $\rho$ (Shape/Anisotropy): Increasing $\rho$ stretches the critical region along the crossover curve. This results in wider peaks in the cumulants and shifts their location further away from the critical point in $\mu_B$.
  • Peak Height vs. Width: The height of the peaks is primarily controlled by $w$ and the distance to the critical point ($\Delta T_f$), while the width and position are strongly controlled by $\rho$ and $w$.

• Dependence on Freeze-out Distance ($\Delta T_f$):

  • Increasing the temperature difference between the crossover and freeze-out ($\Delta T_f$) significantly reduces the magnitude of the critical signal (scaling roughly as $\Delta T_f^{6/5 - k}$ for order $k$).
  • This highlights the sensitivity of the signal to the "memory" of the system; if non-equilibrium effects cause the system to freeze out further from the CP, the signal diminishes rapidly.

• Fourth Cumulant ($\Delta \omega_{4p}$) Signatures:

  • In many scenarios, $\Delta \omega_{4p}$ exhibits a characteristic negative dip followed by a positive peak as $\mu_B$ increases, a signature predicted for the CP.
  • However, for large values of $w$ and $\rho$, subleading contributions can dominate, resulting in a broad positive peak without a preceding negative dip. This implies that the absence of a negative dip does not rule out a CP, but rather constrains the non-universal parameters.

• Resonance Feed-down:

  • Calculations including "daughter protons" from resonance decays (Appendix B) show that while they introduce quantitative changes, the qualitative features (peak shapes, signs, and parameter dependencies) remain unchanged.

5. Significance and Outlook

• Framework for Data Analysis: This work establishes a rigorous, quantitative framework to interpret future experimental data from RHIC and the upcoming FAIR/CBM experiments. It allows for a Bayesian analysis where experimental measurements of cumulants can be used to constrain the non-universal parameters ($w, \rho, \mu_c$) of the QCD EoS.
• Beyond Phenomenology: By replacing phenomenological couplings with the MaxEnt procedure, the paper removes a major source of theoretical uncertainty, linking observables directly to the thermodynamic EoS.
• Future Directions: The authors identify the assumption of thermal equilibrium as the primary limitation. Future work must incorporate Hydro+ formalism to account for critical slowing down and non-equilibrium dynamics. Additionally, the inclusion of a hadronic afterburner (e.g., SMASH) is necessary to fully model kinetic freeze-out effects.

In summary, this paper bridges the gap between the universal critical behavior of the 3D Ising model and specific experimental observables in heavy-ion collisions, providing a tool to not only detect the QCD critical point but also to map the detailed structure of the QCD phase diagram near it.