Event-by-Event Multiplicity Fluctuations in Heavy-Ion Collisions Using Modified HIJING Monte Carlo Generator

This paper utilizes a modified HIJING Monte Carlo generator to demonstrate that event-by-event multiplicity fluctuations serve as a sensitive diagnostic tool for distinguishing between hot and cold media, validating energy loss models, and identifying signatures of a first-order phase transition in relativistic heavy-ion collisions.

The Gist

Imagine you are a detective trying to figure out what happened inside a sealed, super-hot box just by looking at the debris flying out of it after an explosion. That is essentially what physicists do when they smash heavy atoms (like gold) together at nearly the speed of light. They are trying to recreate the conditions of the universe just microseconds after the Big Bang to see if they can create a "soup" of fundamental particles called Quark-Gluon Plasma (QGP).

This paper is about a new detective tool the authors built to solve this mystery. Here is the breakdown in simple terms:

1. The Big Mystery: Hot Soup vs. Cold Ice

When these atoms collide, they create a tiny, incredibly hot fireball.

• The Goal: Scientists want to know if this fireball turns into a super-hot, fluid-like "soup" (QGP) or if it stays more like a "cold" gas of particles.
• The Phase Transition: There is a specific moment where the matter might switch from one state to another, like water turning into steam. The scientists are looking for a "First-Order Phase Transition," which is a dramatic, sudden switch (like water boiling violently) rather than a slow, smooth change.

2. The Problem: The Box is Too Small to See Inside

You can't stick a thermometer inside this fireball because it lasts for a fraction of a second and is too small. So, how do you know what's inside?

• The Old Way: Scientists used computer simulations (like the HIJING generator) to guess what happens. But the old version of this computer program was a bit like a map that didn't show traffic jams or roadblocks. It didn't fully account for how particles lose energy when they fly through the hot soup.
• The New Tool: The authors in this paper took that computer program and "tuned" it. They added realistic rules about how particles lose energy when they crash into the hot soup (radiation) or bump into other particles (collisions). They also added a rule to simulate that dramatic "boiling" phase transition.

3. The Detective Work: Counting the Debris (Multiplicity Fluctuations)

Instead of just counting how many particles come out, the authors looked at how much the number of particles varies from one crash to the next.

• The Analogy: Imagine you are at a crowded concert.
  • Scenario A (Cold/Ice): If the crowd is calm and orderly, the number of people leaving through the exit every minute is very predictable. It's always about 50 people. There is very little fluctuation.
  • Scenario B (Hot Soup/Boiling): If the crowd is chaotic, panicked, or if the venue is "boiling over," the number of people leaving varies wildly. One minute 20 people leave, the next minute 80 leave. The fluctuation is huge.

The authors found that by measuring these "wild swings" in particle counts, they could tell if the collision created a hot, dense soup or a cold, calm gas.

4. The "Smoking Gun": The Phase Transition

The most exciting part of the paper is finding the "boiling point."

• The authors simulated collisions at different energy levels.
• They found that right at the specific energy where the "First-Order Phase Transition" (the dramatic switch) happens, the fluctuations go crazy.
• It's like watching a pot of water. As it gets close to boiling, you see big bubbles forming and popping. The surface becomes very unstable. Similarly, in the particle collisions, the number of particles coming out becomes extremely unpredictable right at the transition point.

5. Why This Matters

This study is like upgrading the detective's magnifying glass.

• Better Maps: By improving the computer simulation, they can now predict what real experiments (like those at the Large Hadron Collider or the upcoming FAIR facility) should see.
• Finding the Critical Point: The ultimate goal is to map the entire "phase diagram" of nuclear matter. This paper suggests that looking at how much the particle counts jump around is a very sensitive way to find the "Critical Point"—the specific spot in the universe's history where matter behaves in its most extreme way.

Summary

The authors took a computer model of particle collisions, made it more realistic by adding rules about how particles lose energy in hot and cold environments, and then used it to count particles. They discovered that chaos in the particle count is the best signal that a dramatic change in the state of matter (like a phase transition) is happening. It's a new, sensitive way to "feel" the temperature and state of the universe's smallest building blocks without ever touching them.

Technical Summary

Here is a detailed technical summary of the paper "Event-by-Event Multiplicity Fluctuations in Heavy-Ion Collisions Using Modified HIJING Monte Carlo Generator" by Y.A. Rusak and L.F. Babichev.

1. Problem Statement

The primary goal of high-energy nuclear physics is to map the Quantum Chromodynamics (QCD) phase diagram, specifically searching for the Quark-Gluon Plasma (QGP) and the Critical End-Point (CEP). While lattice QCD predicts a smooth crossover transition at zero baryon chemical potential ($\mu_B$), it suggests a first-order phase transition at high $\mu_B$, separated from the crossover region by a critical point.

The challenge lies in identifying experimental observables that are:

• Sensitive to the order of the phase transition (distinguishing between crossover and first-order).
• Capable of differentiating between "hot" (QGP) and "cold" hadronic media.
• Measurable in current and future experiments (e.g., PHENIX, STAR, MPD, FAIR).

Standard Monte Carlo generators often lack the specific mechanisms to simulate the complex interplay of parton energy loss and phase transition dynamics required to model these fluctuations accurately.

2. Methodology

The authors employed a modified version of the HIJING (Heavy Ion Jet INteraction Generator) v.1.41 Monte Carlo event generator to simulate central (0–5%) Au+Au collisions across a center-of-mass energy range of $\sqrt{s_{NN}} = 20–200$ GeV.

Key Modifications to HIJING:

• Parton Energy Loss Models: The generator was updated to incorporate energy loss mechanisms for partons traversing both hot and cold media.
  • Hot Medium (QGP): Implemented a total energy loss ($dE/dz$) comprising both collisional and radiative components.
    • Radiative Loss: Modeled using the BDMPS (Baier-Dokshitzer-Mueller-Peigné-Schiff) approach, distinguishing between finite ($L < L_{CR}$) and infinite ($L > L_{CR}$) medium lengths.
    • Collisional Loss: Modeled based on interactions with thermal partons, dependent on temperature ($T$) and Debye screening mass.
  • Cold Medium: Modeled using radiative energy loss only.
  • Temperature Profile: The medium temperature $T$ was defined as a function of collision energy and transverse radius, fitted to hydrodynamic data.
• First-Order Phase Transition Simulation: To simulate a first-order phase transition, the authors introduced a probabilistic mechanism based on the binodal point.
  • A probability distribution function, $\omega_i$, was applied to the center-of-mass energy of individual events.
  • This function determines the likelihood of a transition from cold matter to hot matter based on the proximity of the event's energy to the binodal energy $(\sqrt{s_{NN}})_b$, governed by a fluctuation width parameter $\sigma^2$.

Analysis Framework:

The study analyzed event-by-event multiplicity fluctuations of charged hadrons ($K^\pm, \pi^\pm, p, \bar{p}$). The analysis focused on cumulants of the multiplicity distribution:

• $C_1$: Mean multiplicity ($\langle N \rangle$).
• $C_2$: Variance ($\langle (\delta N)^2 \rangle$).
• $C_3$: Skewness-related moment ($\langle (\delta N)^3 \rangle$).
• Key Ratios: The study examined the ratios $C_2/C_1$ (related to scaled variance) and $C_3/C_2$ (related to skewness) as functions of collision energy and kinematic cuts ($|y|$ and $p_T$).

3. Key Contributions

1. Integration of Energy Loss Formalisms: Successfully integrated finite and infinite medium energy loss models (radiative and collisional) into the HIJING transport model, allowing for a realistic simulation of parton quenching in different matter states.
2. Probabilistic Phase Transition Modeling: Developed a method to simulate first-order phase transitions event-by-event by introducing a stochastic probability distribution for the phase state based on the binodal line.
3. Diagnostic Tool Validation: Demonstrated that multiplicity fluctuation ratios serve as a sensitive diagnostic tool to distinguish between:
   • Hot (QGP) vs. Cold media.
   • Scenarios with and without phase transitions.
   • Different kinematic regions.

4. Results

The simulations, involving 3.5 million events per data point, yielded the following findings:

• Sensitivity to Medium Temperature: Multiplicity fluctuations increase with rising temperature (and consequently, higher parton energy loss in hot, dense matter). This effect becomes more pronounced at higher collision energies ($\sqrt{s_{NN}}$).
• Differentiation of Models: The fluctuation ratios ($C_2/C_1$ and $C_3/C_2$) clearly distinguish between various models:
  • "No quenching" (no medium effects).
  • "Cold" matter models.
  • "BDMPS" (radiative only).
  • "Full" models (radiative + collisional).
• Detection of First-Order Phase Transition:
  • The presence of a binodal point and first-order phase transition creates a distinct signature in the fluctuation ratios.
  • The ratio of the third to the second cumulant ($C_3/C_2$) shows a more pronounced response to the phase transition than the second-to-first ratio.
  • Kinematic Dependence: Expanding the kinematic acceptance (e.g., increasing rapidity $|y|$ or lowering $p_T$ cuts) significantly enhances the difference in behavior near the phase transition region, making the signal clearer.

5. Significance

This work provides a robust theoretical framework for interpreting data from heavy-ion collision experiments aimed at mapping the QCD phase diagram.

• Critical Point Search: The identified sensitivity of higher-order cumulant ratios ($C_3/C_2$) to first-order phase transitions offers a concrete observable for experimentalists searching for the Critical End-Point (CEP).
• Model Validation: The modified HIJING generator allows researchers to test the adequacy of different energy loss models against experimental fluctuation data.
• Experimental Strategy: The results suggest that experiments should prioritize analyzing fluctuations in broader kinematic regions and focus on higher-order cumulants to maximize the detection probability of phase transitions and mixed-phase formation.

In conclusion, the authors establish that event-by-event multiplicity fluctuations, particularly when analyzed via higher-order cumulant ratios in extended kinematic regions, are a powerful and sensitive probe for diagnosing the state of matter (hot vs. cold) and identifying the nature of phase transitions in relativistic heavy-ion collisions.