Particle spectra in the integrated hydrokinetic model at RHIC Beam-Energy-Scan energies

This study utilizes an extended Integrated HydroKinetic Model to analyze light-hadron production in Au+Au collisions at RHIC Beam-Energy-Scan energies, demonstrating that both crossover and first-order phase transition equations of state can similarly describe soft particle spectra when parameters are adjusted, with the most significant distinctions appearing in proton and kaon yields at the lowest energy of 7.7 GeV.

The Gist

Imagine you are trying to understand what happens when two massive, high-speed trains crash into each other. But instead of steel and glass, these trains are made of atomic nuclei (gold atoms), and the crash happens at nearly the speed of light.

When they smash together, they don't just shatter; they melt. For a split second, the matter inside them turns into a super-hot, super-dense soup called Quark-Gluon Plasma (QGP). This is the state of matter that existed just microseconds after the Big Bang.

This paper is like a detective story where scientists are trying to figure out the "rules of the road" for this cosmic soup. They are using a super-computer model called iHKMe (Integrated HydroKinetic Model) to simulate these crashes and see if their theory matches what real experiments at the RHIC (Relativistic Heavy Ion Collider) actually see.

Here is the breakdown of their investigation, explained simply:

1. The Big Question: Is the Soup Smooth or Chunky?

Scientists have two main theories about how this plasma turns back into normal matter (like protons and neutrons) as it cools down:

• The "Smooth Slide" (Crossover): Imagine sliding down a gentle, grassy hill. The transition from hot soup to cold solid is gradual and seamless.
• The "Cliff Edge" (First-Order Phase Transition): Imagine driving off a cliff. The transition is sudden and dramatic. This theory suggests there might be a "critical point" in the universe where matter behaves very strangely before settling down.

The researchers wanted to see which of these two "maps" (Equations of State) fits the data better.

2. The Challenge: Timing is Everything

In high-energy crashes (like at the LHC), the trains overlap for a tiny fraction of a second, and the soup forms almost instantly. It's like a flashbang.

But in this study, they are looking at lower energy crashes (the "Beam Energy Scan"). Here, the trains move slower. They take longer to fully crash into each other.

• The Analogy: Think of it like two people trying to hug. At high speed, they collide and lock arms instantly. At low speed, they might be brushing shoulders for a long time before they actually hug.
• The Problem: Because they take longer to overlap, the "soup" doesn't form instantly. It takes time to "thermalize" (get hot and mixed). The scientists had to figure out exactly when the chaos turns into order. They treated this "mixing time" as a dial they could turn to make their model fit reality.

3. The Experiment: Tuning the Radio

The team ran thousands of simulations, tweaking the "dials" (parameters) of their model. They were looking for the perfect setting that would make their simulated particle speeds match the real data collected by the STAR experiment.

They focused on three main "dials":

• When does the mixing start? (Thermalization start time).
• How fast does it mix? (Relaxation time).
• When does the soup stop acting like a fluid and start acting like individual particles? (Switching time).

4. The Findings: What Did They Discover?

A. The "Mixing Time" is Consistent
Surprisingly, no matter how slow the crash was (within the range they studied), the time it took for the chaos to turn into a smooth fluid was always about 1 femtosecond (that's a quadrillionth of a second). It's like a universal rule: once the crash really starts, the "mixing" takes the same amount of time.

B. Both Maps Work (Mostly)
Whether they used the "Smooth Slide" theory or the "Cliff Edge" theory, they could get the model to match the data pretty well. It's like driving from New York to Boston; you can take the highway (smooth) or the scenic route (cliff), and you'll get there.

• However: At the lowest energy (7.7 GeV), the two maps started to look different. The "Cliff Edge" map predicted slightly different numbers of protons and kaons (heavy particles) than the "Smooth Slide" map. This suggests that at these lower energies, the "shape" of the transition matters more.

C. The "Antiproton" Mystery
The model struggled a bit with antiprotons (the "anti-matter" twins of protons). The computer model kept predicting too few of them compared to reality. The scientists suspect this is because their model for the "aftermath" (the stage where particles fly apart) isn't perfect at handling how matter and anti-matter destroy each other.

5. The Bottom Line

This paper is a success story of "tuning." The scientists proved that their complex computer model can successfully describe these low-energy crashes.

• Key Takeaway: The universe behaves like a fluid even at these lower energies, but the transition from "fluid soup" to "solid particles" is sensitive to the energy level.
• Future: The model works great for the energies they tested, but if they go even lower (to the energies of future experiments), the "mixing" might never finish, and the soup might never fully form. That's the next mystery they plan to solve.

In a nutshell: They built a high-tech video game engine to simulate atomic crashes. They found that the "physics engine" works well for both smooth and sudden transitions, but the details of the crash (how many particles are made) change depending on how hard the nuclei hit each other. They are now ready to use this engine to hunt for the "Critical Point"—the holy grail of nuclear physics.

Technical Summary

Here is a detailed technical summary of the paper "Particle spectra in the integrated hydrokinetic model at RHIC Beam-Energy-Scan energies."

1. Problem Statement

The Relativistic Heavy Ion Collider (RHIC) Beam Energy Scan (BES) program aims to explore the Quantum Chromodynamics (QCD) phase diagram at high baryon chemical potentials ($\mu_B$) and moderate temperatures ($T$). A primary objective is to locate the critical endpoint and identify signals of a first-order phase transition between the Quark-Gluon Plasma (QGP) and the Hadron Resonance Gas (HRG).

However, analyzing data from these lower-energy collisions ($\sqrt{s_{NN}} = 7.7 - 39$ GeV) presents significant challenges:

• Breakdown of Standard Assumptions: At lower energies, the nuclear overlap time is significantly longer ($>1.5$ fm/c) compared to LHC energies, leading to strong longitudinal inhomogeneities and a violation of boost invariance (Bjorken picture).
• Non-Equilibrium Dynamics: The system may not reach full thermal equilibrium before hadronization. The pre-equilibrium stage can dominate the evolution, making standard hydrodynamic models insufficient.
• Model Sensitivity: Distinguishing between different Equations of State (EoS)—specifically a smooth crossover versus a first-order phase transition—requires precise modeling of the system's evolution from initial impact to final particle detection.

2. Methodology

The authors utilize an extended Integrated HydroKinetic Model (iHKMe), a hybrid framework designed to handle the full space-time evolution of heavy-ion collisions from non-equilibrium initial states to final hadronic decoupling.

Key Components of the Model:

• Initial State & Pre-Equilibrium: Instead of assuming immediate thermalization, the model uses UrQMD (Ultra-relativistic Quantum Molecular Dynamics) cascade simulations to describe the initial non-equilibrium dynamics.
• Thermalization Stage: The model employs a continuous transition mechanism where the system evolves from a non-equilibrium UrQMD phase to a hydrodynamic phase. This is governed by a weight function $P(\tau)$ based on a relaxation time approximation, blending the stress-energy tensors of the non-equilibrium and hydrodynamic components.
  • Key parameters include the onset of thermalization ($\tau_0$), the end of thermalization ($\tau_{th}$), and the relaxation time ($\tau_{rel}$).
• Hydrodynamic Expansion: Once thermalized, the system evolves via viscous relativistic hydrodynamics (solved using the vHLLE code). The model includes shear viscosity ($\eta/s$) but neglects bulk viscosity and heat flow for this study.
• Particlization & Afterburner: When the energy density drops below a switching threshold ($\epsilon_{sw}$), the fluid converts to hadrons via a Cooper-Frye-like algorithm. The system then undergoes a hadronic cascade (UrQMD) to simulate the final "afterburner" stage.
• Equations of State (EoS): Two distinct EoS are tested:
  1. Crossover (CO): Features a smooth transition (chiral EoS).
  2. First-Order Phase Transition (PT): Features a soft EoS with a mixed phase.

Calibration Strategy:

• The authors performed a sensitivity analysis using Latin Hypercube Sampling on 750 parameter sets for Au+Au collisions at $\sqrt{s_{NN}} = 14.5$ GeV (20–30% centrality).
• Parameters varied included $\tau_0$, $\tau_{rel}$, $\eta/s$, $\epsilon_{sw}$, and the initial smearing parameter $R$.
• The model was subsequently fine-tuned against experimental data (STAR Collaboration) across the full BES energy range ($7.7$to$39$ GeV).

3. Key Contributions

• Unified Framework for Low Energies: The paper demonstrates the application of iHKMe to the intermediate/low-energy regime where the separation between pre-equilibrium, hydrodynamic, and hadronic stages is blurred.
• Parameter Sensitivity Analysis: The study quantifies how specific model parameters influence observables. Notably, it identifies that antiproton yields are highly sensitive to the switching energy density ($\epsilon_{sw}$) and the initial smoothing parameter ($R$).
• Thermalization Timescale Determination: The authors establish that thermalization begins slightly before full nuclear overlap ($\tau_0 \approx 0.75 \tau_{overlap}$) and lasts approximately 1 fm/c across all BES energies.
• EoS Discrimination: The work clarifies the conditions under which the choice of EoS (crossover vs. first-order) significantly impacts observables.

4. Results

• Thermalization Duration: The thermalization stage ($\tau_{th} - \tau_0$) is consistently found to be $\approx 1$ fm/c, regardless of collision energy. However, the start time ($\tau_0$) shifts later at lower energies, scaling with the nuclear overlap time.
• Equation of State Comparison:
  • High Energies ($>19.6$ GeV): Both the Crossover and First-Order EoS provide similarly good descriptions of transverse momentum ($p_T$) spectra. The rapid expansion at these energies masks the differences in the phase transition region.
  • Low Energy ($7.7$ GeV): Significant differences emerge. The First-Order EoS (softer) and Crossover EoS yield distinct results for proton and kaon yields. The softer EoS generally produces higher mid-rapidity yields but weaker transverse flow.
• Model vs. Data:
  • The model successfully reproduces pion and kaon spectra across the energy range.
  • Discrepancies: The model systematically overestimates proton yields and underestimates antiproton yields at low $p_T$ (deviations up to 40%). The authors attribute this to an imbalance in baryon creation/annihilation during the hadronic afterburner stage or inaccuracies in the baryon density distribution at the particlization hypersurface.
• Future Implications: Extrapolating the trends suggests that at energies below $\sim 4.5$ GeV (accessible at STAR FXT and FAIR CBM), the system will spend so much time in the transition region that full hydrodynamization will not occur for most of the fireball, leading to spectra dominated by non-thermal components.

5. Significance

• Validation of Hybrid Models: The study confirms that hybrid models integrating microscopic transport and hydrodynamics are essential for interpreting RHIC BES data, as pure hydrodynamics fails to capture the extended pre-equilibrium dynamics at lower energies.
• Constraints on QCD Phase Diagram: By showing that both EoS types can fit the data with slight parameter adjustments, the paper suggests that transverse momentum spectra alone are insufficient to definitively distinguish between a crossover and a first-order phase transition at these energies.
• Future Directions: The authors conclude that to detect the critical point or phase transition, future analyses must focus on bulk observables such as elliptic flow ($v_2$) and femtoscopy (HBT radii), which are likely more sensitive to the softening of the EoS and the specific trajectory of the system in the $T-\mu_B$ plane.
• Methodological Benchmark: The identification of the thermalization timescale ($\sim 1$ fm/c) and the specific sensitivity of antiprotons to the switching condition provides a crucial benchmark for future theoretical developments in heavy-ion collision modeling.