Variations of the crossover and first-order phase transition curve in modeling the QCD equation of state

This paper presents a method to construct a smooth QCD equation of state that incorporates a critical point and a first-order phase transition line terminating at it, while ensuring consistency with lattice QCD results, experimental heavy-ion collision data, and the expected 3D Ising universality class critical behavior.

The Gist

Imagine the universe's most fundamental building blocks—quarks and gluons—as a giant, cosmic pot of soup.

At extremely high temperatures (like the moments right after the Big Bang or inside a particle collider), this soup is a free-flowing, chaotic liquid called Quark-Gluon Plasma. The ingredients (quarks) are swimming freely. But as the soup cools down, it "freezes" into solid chunks called Hadrons (like protons and neutrons), where the ingredients are locked together in tight little cages.

The big question physicists are asking is: How does this soup change from liquid to solid?

The Two Scenarios

There are two ways this change can happen, depending on how much "stuff" (baryon density) is in the pot:

1. The Smooth Transition (Crossover): At low density, the soup doesn't suddenly snap into ice. It just gets thicker and thicker, like honey cooling down, until it's eventually solid. There's no sharp line where it changes; it's a gradual blur.
2. The Sharp Snap (First-Order Transition): At high density, the change is violent. It's like water suddenly boiling into steam or freezing into ice. One moment it's liquid, the next it's gas/solid. There is a sharp boundary line between the two states.

The Missing Piece: The Critical Point

Physicists suspect that somewhere in the middle of the temperature and density map, these two scenarios meet at a special spot called the Critical Point.

Think of it like the top of a mountain pass. On one side, the path is a gentle slope (the smooth transition). On the other side, the path is a steep cliff (the sharp snap). The very top of the pass is the Critical Point. If you stand there, the rules of physics get weird and "critical."

The Problem: We can't easily simulate this Critical Point on computers because the math gets too messy (a problem called the "sign problem"). So, we have to build a model—a "map"—that guesses where this point is and what the path looks like around it.

What This Paper Does

The authors, Joseph Kapusta and Shensong Wan, are cartographers. They are trying to draw the most accurate map possible of this QCD soup.

1. The "Smooth Background" (The Terrain)
First, they build a smooth, boring background map using known data from supercomputers (Lattice QCD) and theoretical physics. This represents the soup when it's behaving normally, far away from the Critical Point.

2. The "Critical Patch" (The Storm)
Next, they take a mathematical recipe based on the 3D Ising Model (a famous model used to describe how magnets flip or how water boils). They use this recipe to create a "storm zone" around the Critical Point. This storm zone has special properties (called critical exponents) that tell us how the soup behaves right at the edge of the change.

3. The "Window" (The Blend)
The tricky part is stitching the "Storm Zone" onto the "Smooth Background" without creating a tear in the fabric of reality. They use a mathematical "window function." Imagine a spotlight that shines brightly on the Critical Point but fades out smoothly as you move away, blending the weird critical physics into the normal physics seamlessly.

The New Twist: Redrawing the Boundary

In previous maps, the line separating the liquid and solid phases (the phase boundary) looked like an upside-down "U" that stopped abruptly. It didn't quite make sense intuitively because it didn't connect back to the temperature axis in a way that matched what we see in experiments.

In this paper, the authors tried two new ways to draw that boundary line:

• Condition A: They used a mix of entropy (disorder) and density to draw the line.
• Condition B: They used energy density to draw the line.

The Result: These new lines look more like a gentle curve that swoops down and connects smoothly to the temperature axis. It's like realizing the mountain pass doesn't just end; it flows naturally into the valley.

Why Does This Matter?

The ultimate goal is to help scientists at the Relativistic Heavy Ion Collider (RHIC) and the Large Hadron Collider (LHC).

These labs smash heavy atoms together to recreate the Big Bang soup. They are looking for the "Critical Point" by scanning different energy levels. If the soup behaves in a specific, wobbly way (fluctuations) as they change the energy, it's a sign they've hit the Critical Point.

The equation of state (the map) created in this paper is the instruction manual for the computer simulations that model these collisions. By feeding this new, more flexible map into the simulations, physicists can:

• Predict exactly what signals to look for.
• Compare their computer models with real experimental data (the "freeze-out" points where the soup stops changing).
• Finally pinpoint if the Critical Point exists and exactly where it is hiding on the map.

In a Nutshell

This paper is about refining the GPS navigation system for the subatomic world. The authors took an existing map, fixed the way the "Critical Point" connects to the rest of the world, and made sure the route looks realistic. This helps experimentalists know exactly where to drive their particle colliders to find the most mysterious feature of the universe: the point where matter changes from a smooth flow to a sharp snap.

Technical Summary

1. Problem Statement

The Quantum Chromodynamics (QCD) phase diagram describes the transition between hadronic matter (confined quarks and gluons) and the quark-gluon plasma (deconfined state).

• Knowns: Lattice QCD calculations confirm that at vanishing baryon chemical potential ($\mu_B = 0$), the transition is a rapid crossover at $T \approx 155\text{--}160$ MeV.
• Unknowns: At high baryon chemical potential, many effective models predict a first-order phase transition. Theoretical consensus suggests these two regimes meet at a Critical Point (CP). However, the location of this CP and the precise shape of the phase boundary (the line separating the crossover and first-order regions) remain unknown.
• Challenge: Constructing a consistent Equation of State (EoS) that:
  1. Matches lattice QCD results at low $\mu_B$.
  2. Matches perturbative QCD at high $T$ and/or high $\mu_B$.
  3. Incorporates a critical point with correct universality class properties (3D Ising/liquid-gas).
  4. Defines a phase boundary ($\mu_x(T)$) that behaves physically (e.g., intersecting the temperature axis) and aligns with experimental freeze-out data from heavy-ion collisions (RHIC, LHC).

Previous work (Ref. [10]) successfully embedded a critical point but resulted in a phase coexistence curve with an "inverted U-shape" in the $T-n$ plane, where the analytic extension of the phase boundary did not intersect the temperature axis as intuitively expected.

2. Methodology

The authors build upon the framework established in Ref. [10] (based on Schofield's parametrization) to embed a critical singularity into a smooth background EoS.

A. Embedding the Critical Point

The singular part of the EoS is parametrized using Ising model variables ($R, \theta$) mapped to QCD thermodynamic variables:

• Mapping: Baryon density ($n$) and chemical potential ($\mu$) are expressed in terms of Ising magnetization ($M$) and field ($H$) using scaling variables $R$ (distance from CP) and $\theta$ (angular coordinate).
• Critical Exponents: The model utilizes exponents from the 3D Ising universality class ($\beta \approx 0.326$, $\gamma \approx 1.237$).
• Window Function: A window function $W(\mu, T)$ is introduced to smoothly merge the singular critical part with the background EoS. This function ensures the critical behavior is confined to a specific region near the CP and vanishes in the vacuum. An additional factor is added to suppress residual critical effects above $T_c$.

B. Background Equation of State

The background pressure $P_{BG}$ is constructed by switching between two regimes:

1. Hadronic Phase: Modeled using an excluded volume model (exI).
2. Quark-Gluon Phase: Modeled using perturbative QCD (pQCD) with a 3-loop running coupling.
   A switching function $S(T, \mu)$, dependent on a renormalization scale $M$, determines the contribution of each phase. The parameters are fitted to lattice QCD data for pressure, trace anomaly, and higher-order cumulants ($\chi_2, \chi_4$) at $\mu_B=0$.

C. Defining the Phase Boundary ($\mu_x(T)$)

The core innovation of this paper is redefining the condition that determines the phase coexistence line $\mu_x(T)$.

• Previous Approach (Ref. [10]): Defined by constant background density ($n_{BG} = n_c - n_0$). This yielded a curve that did not intersect the $T$-axis.
• New Approaches (This Paper): The authors propose two new criteria to ensure the crossover line intersects the temperature axis and the resulting EoS is an even function of $\mu$:
  • Condition A: Based on the sum of squared normalized entropy and density deviations:
    $$\left(\frac{\tilde{s}}{\tilde{s}_c}\right)^2 + \left(\frac{\tilde{n}}{\tilde{n}_c}\right)^2 = 2$$
  • Condition B: Based on normalized energy density:
    $$\frac{\tilde{\epsilon}}{\tilde{\epsilon}_c} = 1$$

3. Key Contributions

1. Revised Phase Boundary Criteria: The paper introduces Conditions A and B, which resolve the geometric inconsistency of previous models. These conditions produce phase boundaries that intersect both the temperature and chemical potential axes, aligning with intuitive expectations of a liquid-gas type transition.
2. Flexible EoS Construction: The authors demonstrate a flexible framework where the critical point location ($T_c, \mu_c$), critical exponents, and background EoS can be varied independently.
3. Experimental Alignment: By tuning the critical point to $T_c = 120$ MeV and $\mu_c = 670$ MeV, the resulting crossover curves (dashed lines) pass through experimental chemical freeze-out data from central heavy-ion collisions at RHIC and LHC (Andronic et al.).
4. Sensitivity Analysis: The study investigates how variations in $T_c$ and $\mu_c$ affect the shape of the phase boundary, finding that the overall qualitative shape is robust despite shifts in the critical point's exact location.

4. Results

• Phase Diagrams:
  • The new phase boundaries ($\mu_x(T)$) are concave functions with globally negative derivatives.
  • Unlike the "inverted U" shape of the previous model, the coexistence curves in the $T-n$ plane are now right-skewed. This is because the background density $n_{BG}$ increases as $\mu$ increases along the critical curve.
• Thermodynamic Quantities:
  • For $T > T_c$, the EoS and its derivatives are smooth (crossover).
  • For $T < T_c$, discontinuities in pressure, entropy, and density appear at $\mu = \mu_x(T)$, signaling a first-order phase transition.
  • Isotherms of pressure vs. density show the characteristic van der Waals loops below $T_c$ and smooth curves above $T_c$.
• Comparison with Data:
  • The crossover curves derived from Conditions A and B agree well with experimental freeze-out points at low $\mu$.
  • At high $\mu$, the critical curves lie above the freeze-out points, a behavior consistent with previous studies and potentially indicative of changing collision dynamics at lower beam energies.
• Window Function: The inclusion of the extra factor (Eq. 9) successfully eliminates residual critical effects above $T_c$, ensuring the EoS remains smooth in the crossover region.

5. Significance

• Hydrodynamic Simulations: The constructed EoS is ready for implementation in hydrodynamic simulations of heavy-ion collisions. This is crucial for the Beam Energy Scan II (BES-II) at RHIC, which aims to locate the QCD critical point.
• Astrophysical Applications: The EoS can also be applied to numerical simulations of neutron star mergers, where high baryon densities are relevant.
• Theoretical Robustness: By providing multiple valid criteria (A and B) for defining the phase boundary, the paper highlights that while the exact location of the critical point is uncertain, the methodology for embedding it into a consistent thermodynamic framework is robust. The results suggest that the qualitative features of the phase diagram (skewed coexistence curves) are stable across different definitions of the boundary.
• Universality: The work reinforces the utility of the 3D Ising universality class in modeling QCD critical behavior, bridging the gap between lattice QCD (low $\mu$) and effective models (high $\mu$).

In summary, this paper refines the theoretical modeling of the QCD phase diagram by correcting geometric inconsistencies in the phase boundary definition, thereby providing a more reliable tool for interpreting experimental data from heavy-ion collisions and searching for the elusive QCD critical point.