[Paper Review] Fast simulation of jet quenching in ultrarelativistic heavy ion collisions

I. Lokhtin, A. Snigirev|ArXiv.org|Jun 3, 2004

High-Energy Particle Collisions Research3 citations

TL;DR

This paper presents PYQUEN, a fast Monte Carlo simulation tool integrated with PYTHIA6.2 to model jet quenching in ultrarelativistic heavy ion collisions by simulating medium-induced parton energy loss via collisional and radiative mechanisms. The method reproduces experimentally observed high-$p_T$ suppression and low-$p_T$ enhancement in hadron spectra, enabling efficient study of jet quenching at LHC energies with realistic QGP properties.

ABSTRACT

The method for simulation of medium-induced rescattering and energy loss of hard partons in ultrarelativistic heavy ion collisions is developed. The model is realized as fast Monte-Carlo tool implemented to modify standard PYTHIA jet event.

Motivation & Objective

• To develop a fast, event-by-event Monte Carlo simulation tool for jet quenching in ultrarelativistic heavy ion collisions at LHC energies.
• To model medium-induced energy loss of hard partons through collisional and radiative processes in a quark-gluon plasma (QGP) environment.
• To enable efficient testing of LHC observables' sensitivity to jet quenching and detector response by integrating with standard PYTHIA jet generation.
• To provide a user-accessible Fortran routine (PYQUEN.F) for researchers to simulate jet quenching with adjustable QGP parameters such as initial temperature and centrality.

Proposed method

• The model uses a modified PYTHIA6.2 framework to simulate jet production and subsequent parton rescattering in a dynamically evolving QGP medium.
• It incorporates collisional energy loss via the differential cross section (2) and momentum transfer regularization at the Debye scale, with $\sigma \propto \mu_D^2$.
• Radiative energy loss is modeled using the BDMS formalism, with the radiative spectrum (4) and $\hat{q}$-parameter dependence on $\mu_D^2$ and $\lambda_g$, the gluon mean free path.
• Parton trajectories are simulated by stochastically sampling scattering distances using an exponential distribution based on local $\lambda^{-1}(\tau)$, the inverse mean free path.
• Energy loss per scattering is computed as $\Delta E_{\text{tot},i} = \Delta E_{\text{col},i} + \Delta E_{\text{rad},i}$, with transverse momentum kicks calculated from elastic scattering kinematics.
• The simulation terminates when partons escape the medium, cool below $T_c = 200$ MeV, or lose too much energy ($p_T < 2T$), with final-state gluons added to the PYTHIA parton list for fragmentation.

Experimental results

Research questions

• RQ1 How can jet quenching be efficiently simulated in event-by-event fashion for LHC heavy ion collisions with realistic QGP evolution?
• RQ2 To what extent do collisional and radiative energy loss mechanisms reproduce the observed suppression of high-$p_T$ hadrons and enhancement at low-$p_T$?
• RQ3 Can a fast Monte Carlo tool be effectively integrated into standard PYTHIA6.2 to model jet quenching without sacrificing physical accuracy?
• RQ4 How do initial QGP temperature and centrality selection affect the resulting hadron $p_T$ spectra in simulated events?

Key findings

• The inclusion of jet quenching via PYQUEN leads to a clear suppression of high-$p_T$ hadron spectra compared to PYTHIA-only simulations, consistent with RHIC and LHC data.
• The model reproduces a low-$p_T$ enhancement due to medium-induced gluon radiation, particularly in semi-hard and soft channels.
• For central Pb–Pb collisions at $\sqrt{s_{\text{NN}}} = 5.5$ TeV, the simulation shows significant energy loss for partons with initial $E_T > 100$ GeV.
• The model's predictions are sensitive to initial QGP temperature and centrality, with results becoming unreliable for very peripheral collisions ($b \sim 2R_A$).
• The simulation framework successfully integrates with PYTHIA6.2, allowing full event generation including string formation and fragmentation after in-medium parton evolution.
• The tool is publicly available as the Fortran routine PYQUEN.F, enabling reproducible and accessible jet quenching studies across the LHC physics community.