Open heavy-flavor transport and hadronization in heavy-ion collisions
Yu Fu, Tharun Krishna, Weiyao Ke, Steffen A. Bass, Ralf Rapp
TL;DR
This work addresses the need for a unified, end-to-end description of heavy-quark transport through the QGP, from production in ultrarelativistic collisions to hadronization and hadronic rescattering. It combines a $T_{\text{RENTo}}$-based bulk medium with (2+1)D viscous hydrodynamics, nonperturbative $T$-matrix diffusion including medium-induced radiation, and two hadronization schemes—iSCM with fragmentation and RRM with HQET fragmentation—implemented within the LIDO framework. The study demonstrates that heavy-quark observables for $D$ mesons in Pb–Pb collisions at $\sqrt{s_{NN}} = 5.02$ TeV are sensitive to the choice of in-medium potential, with Vector confinement potential (VCP) generally yielding better agreement with $R_{AA}$ and $v_2$ data than Wilson-line correlator (WLC) constraints. The results support the viability of extracting HQ transport coefficients and QGP properties from open heavy-flavor data and set the stage for future event-by-event simulations and bottom-quark extensions. Overall, the framework advances the integration of nonperturbative HQ dynamics with realistic bulk evolution and hadronization, enabling more precise comparisons with LHC measurements.
Abstract
We develop a comprehensive model for heavy-quark evolution in a realistic QGP, from their production in the initial collision to hadronic freeze-out. Heavy-quark transport is described by a Langevin approach including medium-induced radiation, coupled to a 2+1D viscous hydrodynamic bulk evolution. Transport coefficients are obtained from non-perturbative $T$-matrix calculations with resonant correlations near the transition temperature. Hadronization is implemented via two fragmentation+recombination schemes: an improved sudden coalescence model and a resonance recombination model. We present results for key open heavy-flavor observables, i.e., the nuclear modification factor and elliptic flow, and compare to LHC Pb-Pb data at $\sqrt{s_{\NN}}$=5.02\,TeV.}
