Charm quark evolution in the early stages of heavy-ion collisions

This study investigates the impact of early-time charm quark dynamics on D-meson observables in Pb+Pb collisions at $\sqrt{s_{NN}}=5.02$ TeV using the IP-Glasma+MUSIC+UrQMD framework and finds that despite significant momentum broadening in the pre-equilibrium stage, the nuclear modification factor ($R_{AA}$) and elliptic flow ($v_2$) of D-mesons remain only weakly sensitive to these early interactions.

The Gist

The Big Picture: Smashing Atoms to See the Universe's "Soup"

Imagine you are trying to understand how a giant, hot bowl of soup behaves. To do this, you take two massive, frozen blocks of ice (representing atomic nuclei) and smash them together at nearly the speed of light.

When they collide, they don't just shatter; they melt instantly into a super-hot, super-dense liquid called Quark-Gluon Plasma (QGP). This is the state of matter that existed microseconds after the Big Bang. Scientists want to know: How does this soup flow? How sticky is it? How hot does it get?

To answer this, they drop tiny, heavy "marbles" into the soup. In the world of physics, these marbles are Charm Quarks. Because they are heavy, they don't get swept away easily; they plow through the soup, and by watching how they slow down or change direction, scientists can learn about the soup's properties.

The Mystery: Did the Soup Exist Before It Was a Soup?

The main question this paper asks is about the very first split-second of the collision.

1. The "Before" State: Before the soup forms, there is a chaotic, fluctuating mess of energy fields called the Glasma. It's like the moment the ice hits the heat but before it has fully melted into a smooth liquid. It's turbulent and wild.
2. The Question: Do the heavy charm quarks interact with this wild "Glasma" phase? Or do they just wait until the smooth "Soup" (QGP) forms before they start moving?

Some scientists thought the Glasma phase was so short-lived that the charm quarks wouldn't notice it. Others thought it might be a crucial moment where the quarks get a "head start" on their journey.

The Experiment: A Digital Time Machine

The authors built a massive computer simulation (a "digital time machine") to test this. They used a multi-stage framework:

• IP-Glasma: Simulates the chaotic initial crash (the Glasma).
• MUSIC: Simulates the smooth, flowing soup (the QGP).
• UrQMD: Simulates the cooling down phase where the soup freezes back into particles.

They tracked millions of charm quarks, watching how they moved through the Glasma and then the QGP, and compared their final positions to real data from the Large Hadron Collider (LHC).

The Findings: The "Head Start" Didn't Matter Much

Here is the surprising result, explained with an analogy:

The Analogy of the Sprinter:
Imagine a sprinter (the charm quark) running a race through a stadium.

• Scenario A: The sprinter starts running immediately when the gun goes off, even though the track is muddy and chaotic (the Glasma).
• Scenario B: The sprinter waits until the track is perfectly smooth and paved (the QGP) before starting to run.

You might think, "If the sprinter runs through the mud first, they will be exhausted and slower by the time they hit the smooth track!"

The Paper's Conclusion:
The simulation showed that while the sprinter did get muddy and jostled around in the Glasma (they gained some "momentum broadening," meaning their path got a bit wobbly), it didn't really change the final result of the race.

Whether the quark started in the chaotic Glasma or waited for the smooth QGP, the final outcome looked almost identical. The "D-mesons" (the final particles the charm quarks turn into) ended up with the same speed and direction.

Why Was This Surprising?

The authors found that the Glasma phase actually contributed almost as much to the "jostling" of the quarks as the main soup phase did. It was a very energetic, high-temperature environment.

However, because the Glasma phase is so short and the flow of the medium hasn't developed yet, the quarks get "scrambled" randomly. They don't pick up a specific direction from the Glasma. By the time the smooth soup forms, the quarks are just as ready to interact with the soup as if they had waited.

The "So What?" for Everyday People

1. Robustness of Models: This is good news for physicists. It means their current models are robust. They don't need to know the exact complex physics of the very first trillionth of a second to predict the final results accurately.
2. The "Noise" Factor: The paper shows that the "noise" of the early chaotic phase gets washed out by the time the particles are detected. It's like trying to hear a whisper in a hurricane; the hurricane (the main QGP soup) is so loud that the whisper (the early Glasma interaction) doesn't change the final message you hear.
3. Future Work: While the final results didn't change much, the paper notes that we still need to understand the early phase better to be 100% sure. It's like knowing the car engine works, but still wanting to know exactly how the spark plugs fire in the first millisecond.

Summary in One Sentence

Even though heavy quarks get a rough ride in the chaotic first split-second of a particle collision, it turns out that this early turbulence doesn't significantly change how they behave by the time the "soup" cools down, suggesting our current models of the universe's earliest moments are on the right track.

Technical Summary

1. Problem Statement

Heavy quarks (charm and bottom) are produced predominantly in the initial hard scatterings of relativistic heavy-ion collisions. Consequently, they are expected to traverse the entire evolution of the Quark-Gluon Plasma (QGP), potentially retaining imprints of the earliest, pre-equilibrium stages of the collision. However, standard phenomenological models often neglect the pre-equilibrium phase (the "Glasma" stage) due to its short duration ($\tau < 0.4$ fm/c), assuming that heavy quark interactions only begin once the medium thermalizes into a hydrodynamic QGP.

The central question addressed in this study is: Do the interactions of charm quarks during the pre-equilibrium Glasma phase significantly alter final-state observables, specifically the nuclear suppression factor ($R_{AA}$) and the elliptic flow coefficient ($v_2$) of D-mesons? While recent theoretical work suggests the pre-equilibrium phase contributes significantly to momentum broadening, its phenomenological impact on final observables remains unclear.

2. Methodology

The authors employ a sophisticated multi-stage hybrid framework to simulate the evolution of Pb+Pb collisions at $\sqrt{s_{NN}} = 5.02$ TeV.

• Bulk Medium Evolution:

  • Initial State: The IP-Glasma model is used to describe the fluctuating initial state based on the Color Glass Condensate (CGC) framework. It solves classical Yang-Mills equations to generate the energy-momentum tensor of the Glasma.
  • Hydrodynamic Phase: The system transitions to viscous hydrodynamics at $\tau = 0.4$ fm/c using the MUSIC solver. This utilizes a second-order DNMR formulation including shear and bulk viscosity.
  • Hadronic Phase: The system freezes out and evolves through the hadronic phase using the UrQMD transport model.

• Heavy Quark Production and Evolution:

  • Production: Charm quark pairs are generated using PYTHIA at binary collision locations derived from the IP-Glasma events. Nuclear shadowing is accounted for via EPS09 nuclear parton distribution functions (nPDFs).
  • Pre-Equilibrium (Glasma) Evolution: The authors approximate the Glasma as a locally thermalized medium composed of gluons. Charm quarks interact with this medium via Langevin dynamics within the MARTINI event generator. The interaction is governed by drag ($\eta$) and diffusion ($\kappa$) coefficients derived from lattice QCD spatial diffusion coefficients ($D_s$).
  • QGP Evolution: Once the system enters the hydrodynamic phase, charm quarks continue to evolve via Langevin dynamics using the same transport coefficients, now coupled to the fluid velocity field.
  • Formation Time: A formation time $\tau_F = 1/m_c$ (in the quark's rest frame) is implemented. Quarks are bound by a Cornell potential before this time and interact with the medium only after.

• Hadronization:

  • At $T_{th} = 165$ MeV, charm quarks hadronize via a combination of coalescence (recombination with light quarks) and fragmentation (Peterson function).
  • Coalescence probabilities are determined by Wigner functions for mesons and baryons, with a normalization factor $N_{coal}$ tuned to data.

3. Key Contributions

• First Phenomenological Integration: This study provides the first comprehensive phenomenological analysis of charm quark evolution starting from the fluctuating IP-Glasma initial state through to the final hadronic stage, explicitly including the pre-equilibrium phase.
• Quantification of Pre-Equilibrium Impact: The authors rigorously test the sensitivity of $R_{AA}$ and $v_2$ to the inclusion of the Glasma phase, comparing scenarios with and without pre-equilibrium interactions.
• Uncertainty Analysis: The paper systematically evaluates uncertainties arising from:
  • Event-by-event fluctuations vs. smooth averaged backgrounds.
  • Charm quark formation time ($\tau_F$).
  • Coalescence normalization ($N_{coal}$).
  • Momentum dependence of transport coefficients.
  • Lattice QCD uncertainties in the spatial diffusion coefficient ($D_s$).

4. Key Results

• Momentum Broadening vs. Observables:

  • The study confirms that the pre-equilibrium phase contributes significantly to the total transverse momentum broadening of charm quarks. The ratio of non-equilibrium to equilibrium momentum broadening is found to be $\Delta p_T^2|_{\text{non-eq}} / \Delta p_T^2|_{\text{eq}} \approx 0.8 - 1$.
  • Crucially, despite this large contribution to momentum broadening, the final-state observables ($R_{AA}$ and $v_2$) show only weak sensitivity to the pre-equilibrium interactions.

• $R_{AA}$ and $v_2$ Sensitivity:

  • $R_{AA}$: Including pre-equilibrium evolution slightly enhances $R_{AA}$ in the intermediate momentum range ($2 < p_T < 8$ GeV). This counter-intuitive enhancement is attributed to high-momentum kicks (diffusion) in the high-temperature Glasma phase, which compensate for energy loss. However, the magnitude of this effect is small and comparable to systematic model uncertainties.
  • $v_2$: The elliptic flow is virtually unaffected by pre-equilibrium evolution. The authors explain that while the Glasma scrambles the isotropic charm quark distribution via diffusion, the background medium flow has not yet developed anisotropies at this early stage. Thus, charm quarks remain uncorrelated with the spatial anisotropy of the medium, resulting in negligible impact on $v_2$.

• Sensitivity to Model Parameters:

  • Fluctuations: Event-by-event fluctuating backgrounds yield higher $R_{AA}$ and $v_2$ compared to smooth averaged backgrounds, highlighting the necessity of realistic initial conditions.
  • Formation Time: Shorter formation times lead to greater energy loss (lower $R_{AA}$) and higher $v_2$ in the high $p_T$ region, as quarks interact with the medium for a longer duration.
  • Transport Coefficients: The momentum dependence of drag and diffusion coefficients is essential for describing the data. Lattice QCD uncertainties in $D_s$ significantly affect the error bands of $v_2$, particularly in central collisions.

5. Significance and Conclusion

The paper concludes that while the pre-equilibrium Glasma phase is a dynamic environment that significantly broadens the momentum of heavy quarks, it does not serve as a strong discriminator for final-state D-meson observables ($R_{AA}$ and $v_2$) within current experimental precision and theoretical uncertainties.

• Implication: The "heavy quark puzzle" (the simultaneous description of $R_{AA}$ and $v_2$) is likely dominated by the properties of the thermalized QGP and the hadronization mechanism rather than the pre-equilibrium dynamics.
• Future Outlook: The authors suggest that to truly probe the pre-equilibrium stage, one must look beyond standard $R_{AA}$ and $v_2$. They propose that two-particle azimuthal correlations or more sophisticated probes of the non-Abelian gauge fields (beyond the thermalized approximation used here) are necessary. Future work will explore radiation-improved Langevin models and the "dead-cone effect" to refine these predictions.

In summary, this study establishes a robust baseline for heavy-flavor phenomenology, demonstrating that current heavy-ion observables are largely insensitive to the specific details of the pre-equilibrium charm quark evolution, provided the bulk medium evolution is modeled realistically.