Species-dependent viscous corrections at particlization: A novel relaxation time approximation approach

This paper demonstrates that a novel generalized relaxation time approximation introducing species-dependent viscous corrections significantly alters identified hadron yields and transverse momentum spectra in p-Pb and Pb-Pb collisions while preserving bulk flow descriptions, thereby offering new sensitivity for Bayesian inference without compromising existing constraints.

The Gist

Imagine a massive, chaotic dance party where millions of tiny particles are swirling around at nearly the speed of light. This is what happens when scientists smash heavy atoms (like lead) together in giant machines called particle accelerators. For a split second, they create a "soup" of energy and matter called the Quark-Gluon Plasma (QGP). This soup behaves like a perfect, super-fluid liquid.

As this hot soup cools down, it freezes into individual particles (like protons, pions, and kaons) that fly out to be detected. Scientists use complex computer simulations to predict what this "freezing" process looks like.

This paper is about fixing a specific glitch in the computer code used to simulate that freezing process. Here is the story in simple terms:

1. The Old Problem: The "One-Size-Fits-All" Mistake

Imagine you are a bouncer at a club trying to let people out. In the old simulation method, the bouncer used a single rule for everyone: "Everyone takes exactly 10 seconds to leave the room."

In physics, this rule is called the Relaxation Time Approximation (RTA). It assumes that every particle, regardless of its weight or speed, relaxes (calms down) at the same rate.

The Flaw: In reality, a heavy particle (like a proton) is like a tired elephant trying to leave a room—it moves slowly. A light particle (like a pion) is like a hyperactive mouse—it zips out quickly. The old "one-size-fits-all" rule ignored these differences. It also broke the laws of physics (specifically, the conservation of energy and momentum) when scientists tried to make the exit time depend on how fast the particle was moving.

2. The New Solution: The "Personalized Exit Plan"

The authors of this paper introduced a new, smarter rule (a "Generalized RTA"). Instead of a single rule for everyone, the bouncer now has a personalized plan for every guest based on their weight and speed.

• The Innovation: They added "counter-terms" (think of these as balancing weights) to the math. This allows the simulation to say, "Heavy particles take longer to settle, light particles take less time," without breaking the laws of physics.
• The Result: The simulation now knows that a heavy proton and a light pion should behave differently when the soup freezes.

3. What Happened When They Tested It?

The team ran their new simulation on two types of "parties":

1. Pb-Pb: A huge collision (like a massive stadium crowd).
2. p-Pb: A smaller collision (like a small gathering).

They looked at the "guest list" (the particles produced) and found some fascinating changes:

• The "Flavor" Shift: Because the new rule treats heavy and light particles differently, the ratio of particles changed.

  • Old Simulation: Predicted a certain number of pions (light) vs. protons (heavy).
  • New Simulation: Predicted fewer protons and more pions (or vice versa, depending on the specific settings).
  • Analogy: It's like if the old bouncer let out 100 mice and 10 elephants, but the new bouncer, realizing elephants are slower, lets out 90 mice and 15 elephants. The total number of animals leaving might look similar, but the mix is completely different.

• The "Invisible" Total: Interestingly, if you just count the total number of charged particles (ignoring what type they are), the difference is tiny. The heavy particles and light particles cancel each other out. It's only when you look at the specific types (the "flavors") that you see the new rule making a big difference.

4. Why Does This Matter?

This might sound like a small technical fix, but it's actually a game-changer for two reasons:

1. Better Detective Work: Scientists use these simulations to figure out the properties of the Quark-Gluon Plasma (like how "sticky" or viscous it is). If the simulation has a glitch in how it freezes the particles, the scientists might get the wrong answer about the soup's properties. This new method removes that glitch.
2. Bayesian Inference (The "Tuning" Knob): Scientists use a method called Bayesian inference to tune their models to match real-world data. The old method was like trying to tune a radio with only one knob. This new method adds a new knob (the parameter $\gamma$) that specifically controls how particle mass affects the outcome. This gives scientists a new way to fine-tune their models to match the real universe more accurately.

The Bottom Line

The authors didn't just find a bug; they found a way to make the simulation "speak" the language of different particle species. They showed that how heavy a particle is matters when the universe cools down from a hot soup into solid matter.

By fixing this, they ensure that when we look at the data from particle colliders, we aren't being misled by a clumsy computer model. It's a crucial step toward understanding the fundamental building blocks of our universe.

Technical Summary

1. Problem Statement

In the hybrid modeling of relativistic heavy-ion collisions, the transition from the continuous fluid description of the Quark-Gluon Plasma (QGP) to the discrete hadronic gas (particlization) is a critical step. This process typically employs the Cooper-Frye formula with viscous corrections ($\delta f$) to account for non-equilibrium effects.

The standard approach uses the Anderson-Witting Relaxation Time Approximation (RTA). However, this method suffers from a fundamental inconsistency: if the relaxation time ($\tau_R$) depends on particle momentum (as expected from microscopic interactions) or varies between particle species, the traditional RTA violates the macroscopic conservation laws of energy and momentum. To avoid this violation, previous simulations have been forced to assume constant, momentum-independent relaxation times, ignoring realistic energy dependencies.

2. Methodology

The authors address this limitation by implementing a generalized, multi-species RTA (nRTA) proposed in recent theoretical works [24, 25] into a state-of-the-art hybrid simulation framework.

• Theoretical Framework: The nRTA introduces specific counter-terms to the collision kernel. These terms are constructed from an orthogonal basis of conserved quantities (energy and momentum). This mathematical structure allows the relaxation time $\tau_{R,i}(p)$ to depend arbitrarily on both the particle species ($i$) and its momentum ($p$) while rigorously preserving local energy-momentum conservation.
• Viscous Corrections: The resulting first-order viscous corrections ($\delta f_i$) to the phase-space distribution functions become explicitly dependent on the particle mass ($m_i$) and the momentum dependence parameter $\gamma$. The relaxation time is parametrized as $\tau_{R,i} = t_R (E_i/T)^\gamma$.
• Numerical Implementation:
  • Hybrid Chain: The study utilizes a standard hybrid chain: TRENTo (initial state) $\to$ Free-streaming (pre-equilibrium) $\to$ MUSIC (viscous hydrodynamics) $\to$ iSS (particlization with nRTA corrections) $\to$ UrQMD (hadronic cascade).
  • Systems Studied: Simulations were performed for Pb-Pb ($\sqrt{s_{NN}} = 2.76$ TeV) and p-Pb ($\sqrt{s_{NN}} = 5.02$ TeV) collisions.
  • Parameters: The study varies the phenomenological parameter $\gamma$ (values: 0.0, 0.5, 1.0, 1.5), which controls the power-law dependence of the relaxation time on energy.
  • Baseline: The authors use previously validated parameter sets (MAP values from JETSCAPE and Duke groups) without refitting, isolating the effect of the new particlization prescription.

3. Key Contributions

• First Full Numerical Implementation: This is the first study to integrate the generalized nRTA prescription into a realistic, full-scale hybrid simulation of heavy-ion collisions.
• Species-Dependent Corrections: The work demonstrates that the new formalism naturally generates viscous corrections that differ significantly between particle species (pions, kaons, protons) due to mass dependence, a feature absent in standard RTA.
• Conservation Compliance: It proves that momentum-dependent relaxation times can be used without violating conservation laws, provided the specific counter-terms of the nRTA are employed.

4. Key Results

The study analyzes observables both immediately after particlization (before hadronic decays) and after the hadronic cascade (UrQMD).

• Identified Hadron Yields:
  • The multiplicity of light hadrons is significantly altered by $\gamma$.
  • Pions ($\pi$): Generally show an enhancement in yield as $\gamma$ increases.
  • Kaons ($K$) and Protons ($p$): Generally show a depletion in yield as $\gamma$ increases (though proton trends reverse slightly after the hadronic cascade due to resonance decays).
  • Relative Ratios: Consequently, ratios like $K/\pi$ and $p/\pi$ are strongly modified. For instance, the $p/\pi$ ratio decreases with increasing $\gamma$ before the cascade, offering a new mechanism to tune baryon-to-meson enhancement.
• Transverse Momentum ($p_T$) Spectra:
  • The corrections are $p_T$-dependent. For pions, increasing $\gamma$ enhances low-$p_T$ production ($< 1.5$ GeV/c) and suppresses high-$p_T$ production.
  • Protons exhibit different spectral shapes, with suppression at low $p_T$ in Pb-Pb but a recovery or enhancement at higher $p_T$ in p-Pb collisions.
• Inclusive Charged Particles:
  • The effects on the total charged-particle multiplicity are significantly reduced compared to identified species. This is due to compensating effects: the enhancement of pions is partially canceled by the suppression of kaons and protons.
  • The total $p_T$ spectrum remains relatively flat, indicating that the new prescription redistributes momentum among species without drastically altering the bulk collective flow profile.
• Persistence after Cascade:
  • While the hadronic cascade (UrQMD) dilutes the magnitude of these effects, the species-dependent signatures (particularly in relative yields and spectral shapes) remain visible in the final state.
  • Effects are more pronounced in small systems (p-Pb) due to steeper gradients and larger relative viscous corrections.

5. Significance and Implications

• Bayesian Inference: The controlled deformation of identified hadron observables (without spoiling bulk flow constraints) makes this new prescription ideal for future Bayesian analyses. It introduces new sensitivity directions (species-dependent viscous effects) to constrain QCD transport coefficients more precisely.
• Disentangling Mechanisms: The mass-dependent $\delta f$ provides a kinetic mechanism that directly impacts the $p/\pi$ ratio. This allows researchers to distinguish between effects arising from kinetic freeze-out sampling and those arising from hadronization mechanisms (like coalescence), which is crucial for understanding the "baryon anomaly."
• Future Directions: The authors suggest that this framework should be extended to study electromagnetic radiation, heavy-flavor hadronization, and conserved-charge fluctuations, as these observables are sensitive to the microscopic phase-space distribution at the hydrodynamic interface.

In conclusion, the paper establishes that species-dependent viscous corrections arising from a generalized RTA are not merely theoretical artifacts but leave significant, observable imprints on identified hadron production, necessitating their inclusion in precision analyses of heavy-ion collision data.