No Arabic abstract
We review the statistical hadronization picture for charmonium production in ultra-relativistic nuclear collisions. Our starting point is a brief reminder of the status of the thermal model description of hadron production at high energy. Within this framework an excellent account is achieved of all data for hadrons built of (u,d,s) valence quarks using temperature, baryo-chemical potential and volume as thermal parameters. The large charm quark mass brings in a new (non-thermal) scale which is explicitely taken into account by fixing the total number of charm quarks produced in the collision. Emphasis is placed on the description of the physical basis for the resulting statistical hadronization model. We discuss the evidence for statistical hadronization of charmonia by analysis of recent data from the SPS and RHIC accelerators. Furthermore we discuss an extension of this model towards lower beam energies and develop arguments about the prospects to observe medium modifications of open and hidden charm hadrons. With the imminent start of the LHC accelerator at CERN, exciting prospects for charmonium production studies at the very high energy frontier come into reach. We present arguments that, at such energies, charmonium production becomes a fingerprint of deconfinement: even if no charmonia survive in the quark-gluon plasma, statistical hadronization at the QCD phase boundary of the many tens of charm quarks expected in a single central Pb-Pb collision could lead to an enhanced, rather than suppressed production probability when compared to results for nucleon-nucleon reactions scaled by the number of hard collisions in the Pb-Pb system.
The fast simultaneous hadronization and chemical freeze out of supercooled quark-gluon plasma, created in relativistic heavy ion collisions, leads to the re-heating of the expanding matter and to the change in a collective flow profile. We use the assumption of statistical nature of the hadronization process, and study quantitatively the freeze out in the framework of hydrodynamical Bjorken model with different quark-gluon plasma equations of state.
Charmonium directed flows are studied based on transport model coupled with the realistic three dimensional expansions of the bulk medium. The non-central symmetric nucleus-nucleus collisions can generate the rotating quark-gluon plasma (QGP) with symmetry-breaking longitudinal distributions. In $sqrt{s_{NN}}=200$ GeV Au-Au semi-central collisions, charmonium are primordially produced in the initial hard process, they are mainly dissociated by the initial tilted source with high temperatures and then move out of the bulk medium to keep the early information of the medium. Charmonia are less affected by the hydrodynamic expansions of QGP where its tilted shape is being diluted. This biased dissociation can generate directed flows of $J/psi$ and $psi(2S)$ which are much larger than the values of light charged hadrons and open heavy flavor. They are sensitive and clean to the effect of QGP rapidity-odd initial energy densities.
The interaction of heavy flavor with the quark-gluon plasma (QGP) in relativistic heavy-ion collisions is studied using JETSCAPE, a publicly available software package containing a framework for Monte Carlo event generators. Multi-stage (and multi-model) evolution of heavy quarks within JETSCAPE provides a cohesive description of heavy flavor quenching inside the QGP. As the parton shower develops, a model becomes active as soon as its kinematic region of validity is reached. Two combinations of heavy-flavor energy-loss models are explored within a realistic QGP medium, using parameters which were tuned to describe {it light-flavor} partonic energy-loss.
The suppression and modification of high-energy objects, like jets, in heavy-ion collisions provide an important window to access the degrees of freedom of the quark-gluon plasma on different length scales. Despite increasingly precise and differential measurements of the properties of jets in heavy-ion collisions, however, it has remained challenging to use jets to make unambiguous and model-independent statements about the quark-gluon plasma. Here I will give a personal take on some origins of these challenges, including the difficulty of modelling and biases from jet selection that obfuscate the direct interpretation of jet modification measurements. I will discuss a few model studies that have helped to disentangle the source of non-intuitive effects in measurements, and finally highlight data-driven approaches as an interesting opportunity toward studying the quark-gluon plasma in a model-independent way using jets.
We study charm production in ultra-relativistic heavy-ion collisions by using the Parton-Hadron-String Dynamics (PHSD) transport approach. The initial charm quarks are produced by the Pythia event generator tuned to fit the transverse momentum spectrum and rapidity distribution of charm quarks from Fixed-Order Next-to-Leading Logarithm (FONLL) calculations. The produced charm quarks scatter in the quark-gluon plasma (QGP) with the off-shell partons whose masses and widths are given by the Dynamical Quasi-Particle Model (DQPM) which reproduces the lattice QCD equation-of-state in thermal equilibrium. The relevant cross section are calculated in a consistent way by employing the effective propagators and couplings from the DQPM. Close to the critical energy density of the phase transition, the charm quarks are hadronized into $D$ mesons through coalescence and/or fragmentation depending on transverse momentum. The hadronized $D$ mesons then interact with the various hadrons in the hadronic phase with cross sections calculated in an effective lagrangian approach with heavy-quark spin symmetry. Finally, the nuclear modification factor $rm R_{AA}$ and the elliptic flow $v_2$ of $D^0$ mesons from PHSD are compared with the experimental data from the STAR Collaboration for Au+Au collisions at $sqrt{s_{rm NN}}$ =200 GeV. We find that in the PHSD the energy loss of $D$ mesons at high $p_T$ can be dominantly attributed to partonic scattering while the actual shape of $rm R_{AA}$ versus $p_T$ reflects the heavy quark hadronization scenario, i.e. coalescence versus fragmentation. Also the hadronic rescattering is important for the $rm R_{AA}$ at low $p_T$ and enhances the $D$-meson elliptic flow $v_2$.