Constraining the Phase-Transition EoS using the Energy Dependence of Directed Flow

By combining a hybrid VDF+MIT equation of state with the AMPT-HC transport model and recent experimental directed flow data, this study constrains the hadron-quark phase transition to occur between 5 and 6 times nuclear saturation density and proposes the energy derivative of the mid-rapidity $v_1$ slope as a robust observable for identifying the QCD critical point.

The Gist

Imagine the universe is made of a giant, invisible Lego set. Usually, these Legos (protons and neutrons) stick together to form the matter we see around us. But if you squeeze them hard enough—like using a cosmic hydraulic press—they might break apart into their tiny, fundamental building blocks: quarks and gluons. This transformation is called a phase transition, similar to how ice melts into water or water boils into steam.

The big mystery in physics right now is: At exactly what pressure does this "melting" happen? And what does the "melted" stuff (called Quark-Gluon Plasma) feel like? Is it squishy (soft) or stiff?

This paper is like a detective story where scientists try to solve this mystery by smashing atoms together in a lab and watching how they bounce off each other.

The Tools of the Trade

1. The Cosmic Squeeze (Heavy-Ion Collisions):
   Scientists use giant machines (like the RHIC collider) to smash gold atoms together at nearly the speed of light. This creates a tiny, super-hot, super-dense drop of matter that mimics the conditions of the early universe.

2. The "Traffic Flow" Analogy (Directed Flow):
   When these atoms smash, the debris doesn't just fly straight out; it gets pushed sideways, like cars in a traffic jam trying to squeeze through a narrow gap. This sideways movement is called Directed Flow ($v_1$).

   • The Analogy: Imagine a crowd of people running into a wall. If the wall is soft (like a mattress), they sink in and stop. If the wall is hard (like concrete), they bounce back hard. The "stiffness" of the matter inside the collision determines how hard the particles bounce sideways.

3. The New Recipe (The Hybrid EoS):
   The authors created a new "recipe" for how this matter behaves, called the VDF+MIT EoS.

   • VDF is like a rulebook for how normal nuclear matter acts.
   • MIT is a rulebook for what happens when matter turns into quark soup.
   • They combined them to create a smooth transition, like a recipe that tells you exactly how to turn a cake batter into a soufflé without it collapsing.

The Investigation

The scientists ran computer simulations using their new recipe and compared the results to real data from experiments. They looked at three different "recipes" that assumed the phase transition happens at different pressures:

• Recipe A: Transition happens early (at 3x normal density).
• Recipe B: Transition happens in the middle (at 4-5x normal density).
• Recipe C: Transition happens late (at 5-6x normal density).

The Verdict:
When they compared the "sideways bounce" (directed flow) of protons and Lambda particles in their simulations to the actual data from the RHIC lab, Recipe A was a bust. It predicted the particles would bounce differently than what was observed.

• Recipe C (5-6x density) matched the real-world data almost perfectly.
• Recipe B was okay, but not as good.

Conclusion: The "melting" of protons into quark soup likely doesn't happen until the matter is squeezed to 5 or 6 times its normal density. It rules out the idea that it happens at lower pressures (below 3x).

The "Magic Compass" (A New Discovery)

The most exciting part of the paper is a new tool they invented to find the exact "tipping point" of this phase transition.

• The Problem: It's hard to know exactly when the transition happens just by looking at one crash.
• The Solution: They looked at how the "sideways bounce" changes as you increase the energy of the crash.
• The Analogy: Imagine you are walking up a hill. As you get closer to the peak (the phase transition), your walking speed changes in a specific way. If you plot your speed against the height, there is a moment where the slope of your path hits zero (you stop going up and start going down, or vice versa).

The authors found that the rate of change of this sideways bounce hits zero exactly when the collision energy is high enough to cross the phase transition threshold.

• If you see this "zero crossing" in future experiments, it's a giant neon sign saying: "You have just crossed the line from normal matter to quark soup!"

Why Does This Matter?

1. Mapping the Universe: It helps us draw a better map of the "QCD Phase Diagram," which is like a weather map for the universe, showing where matter is solid, liquid, or plasma.
2. Neutron Stars: This helps us understand what's inside neutron stars. If we know how stiff the matter is, we can predict how big and heavy these stars can get before they collapse.
3. Future Experiments: This new "zero crossing" tool will be a guide for future experiments (like those at HIAF in China or FAIR in Germany) to hunt for the exact spot where the phase transition happens.

In a nutshell: The scientists built a better model of how matter behaves under extreme pressure, proved that the "melting" of atoms into quark soup happens at very high pressures (5-6 times normal), and invented a new "compass" to help future explorers find the exact spot where this cosmic magic happens.

Technical Summary

1. Problem Statement

The structure of the Quantum Chromodynamics (QCD) phase diagram, specifically the location of the critical point and the nature of the first-order phase boundary between hadronic matter and quark-gluon plasma (QGP), remains a central unsolved problem in nuclear physics and astrophysics.

• The Challenge: Deriving an accurate Equation of State (EoS) for nuclear matter at finite baryon chemical potential from first-principles QCD calculations is currently intractable.
• The Gap: While astrophysical observations (neutron stars) and heavy-ion collisions provide constraints, they are often insufficient to uniquely determine the density at which the hadron-to-quark phase transition occurs. Existing models struggle to describe the "softening" of the EoS (a reduction in pressure relative to density) expected during a phase transition without violating causality or failing to match experimental flow data.
• Objective: To constrain the density range of the hadron-quark phase transition and propose a robust, model-independent observable to detect the critical point in future experiments.

2. Methodology

The authors employed a multi-faceted approach combining theoretical modeling, transport simulations, and comparison with experimental data.

• Hybrid Equation of State (VDF+MIT):

  • Hadronic Phase: Utilized the Relativistic Vector Density Functional (VDF) model, based on Landau Fermi liquid theory, which incorporates vector-type interactions to describe normal nuclear matter and the liquid-gas transition.
  • Quark Phase: Incorporated the MIT Bag Model to describe high-density quark matter.
  • Hybridization: The authors constructed a hybrid EoS where the VDF model smoothly transitions into the MIT Bag Model pressure after crossing the spinodal region of the phase transition. This ensures causality is maintained at high densities.
  • Variants: Three specific EoS variants (VDF1, VDF2, VDF3) were created, differing primarily in the critical density ($\rho_c$) at which the phase transition begins (approx. $3\rho_0$, $4\rho_0$, and $6\rho_0$, respectively, where $\rho_0$ is nuclear saturation density).

• Transport Model Simulation (AMPT-HC):

  • The AMPT-HC (A Multi-Phase Transport - Hadronic Cascade) model was used. Unlike the standard AMPT, this mode operates purely in the hadronic sector (no partonic phase), making it suitable for intermediate-to-low energy collisions where hadronic dynamics dominate.
  • Simulations were performed for Au+Au collisions at center-of-mass energies ($\sqrt{s_{NN}}$) of 3.0, 4.5, and 7.7 GeV (the Beam Energy Scan range).
  • The hybrid VDF+MIT EoS was embedded into the transport framework to simulate the spatiotemporal evolution of the collision.

• Observables:

  • Directed Flow ($v_1$): The first Fourier coefficient of the azimuthal particle distribution, quantifying the sideward deflection of matter. The slope of $v_1$ at mid-rapidity ($dv_1/dy|_{y=0}$) is highly sensitive to the stiffness of the EoS.
  • Energy Derivative: A new observable was proposed: the energy derivative of the mid-rapidity $v_1$ slope, denoted as $d(dv_1/dy)/d(\sqrt{s_{NN}})$.

3. Key Results

• Constraints on Phase Transition Density:

  • Astrophysical Consistency: The VDF1+MIT EoS (transition at $\sim 3\rho_0$) was found to be inconsistent with both astrophysical constraints (neutron star mass-radius relations from NICER and GW190425) and heavy-ion data.
  • Flow Data Agreement: Comparisons with RHIC-STAR experimental data for proton and $\Lambda$ hyperon directed flow ($v_1$) showed that:
    • VDF1+MIT ($\sim 3\rho_0$) failed to reproduce the measured $v_1$ distributions.
    • VDF2+MIT ($\sim 4-5\rho_0$) provided a reasonable fit but was less optimal.
    • VDF3+MIT ($\sim 5-6\rho_0$) provided the best agreement with experimental data across the 3.0–7.7 GeV range.
  • Conclusion: The study rules out phase transitions occurring below $3\rho_0$ and strongly suggests the hadron-quark transition occurs in the range of $5\rho_0$ to $6\rho_0$.

• Mechanism of $v_1$ Suppression:

  • The presence of a phase transition introduces a "softening" of the EoS. When the maximum compression density of the collision approaches the transition density, the system spends significant time in the mixed phase.
  • This mixed phase generates a weak attractive potential during compression (enhancing compression) but an inhibitory effect during expansion (reducing transverse pressure).
  • This leads to a suppression of the directed flow slope ($v_1$). The effect is most pronounced when the collision energy is such that the maximum density closely matches the phase transition critical density.

• The "Zero Crossing" Observable:

  • The authors analyzed the rate of change of the $v_1$ slope with respect to collision energy: $d(dv_1/dy)/d(\sqrt{s_{NN}})$.
  • Behavior: This derivative exhibits a non-monotonic trend. It starts negative, crosses zero, becomes positive, and then drops back to negative at higher energies.
  • Significance: The zero crossing point corresponds to the collision energy where the maximum compression density aligns with the phase transition critical point.
  • Differentiation: The zero crossing occurs at significantly lower energies for VDF1 (early transition) and higher energies for VDF3 (late transition). This makes the zero crossing a direct, weakly model-dependent signature of the critical point.

4. Significance and Contributions

• Refined EoS Constraints: The paper provides a rigorous constraint on the QCD phase diagram, narrowing the likely location of the hadron-quark transition to $5\rho_0–6\rho_0$, effectively ruling out lower density scenarios.
• New Diagnostic Tool: The proposal of $d(dv_1/dy)/d(\sqrt{s_{NN}})$ as a "zero-crossing" observable offers a powerful, model-independent method to pinpoint the critical point in future Beam Energy Scan (BES) programs.
• Experimental Guidance: The findings provide clear targets for upcoming high-precision experiments at facilities like HIAF (China), FAIR (Germany), and the continued RHIC-STAR BES program. Identifying the zero crossing in experimental data will serve as a definitive signature of the phase transition.
• Theoretical Framework: The successful integration of the VDF and MIT Bag models into the AMPT-HC transport framework demonstrates a viable pathway for simulating complex phase transitions in heavy-ion collisions, bridging the gap between microscopic QCD properties and macroscopic observables.

In summary, this work utilizes directed flow data to exclude low-density phase transitions and proposes a specific derivative observable to precisely locate the QCD critical point, offering a roadmap for future experimental verification of the QCD phase structure.