Quarkonium production in pp and heavy-ion collisions

Ultra-relativistic heavy-ion collisions produce an extremely hot and dense matter which is supposed to cross over the phase boundary between hadron gas and quark-gluon plasma (QGP). Quarkonium is a flavorless bound state of heavy quark and heavy antiqurk, whose suppression in heavy-ion collisions was suggested as a signature for QGP formation, because the binding of heavy quark pair is sensitive to the environment through the color screening and scattering with thermal partons Matsui:1986dk . However, it was discovered that $J/\psi$ is less suppressed in mid-rapidity as well as in higher energy collisions, which indicates the regeneration of quarkonium from open heavy flavors Manceau:2012ka . So the complete description of quarkonium production in heavy-ion collisions requires suppression as well as regeneration.

Quarkonium production takes two steps. First, heavy quark pair is produced and then the produced heavy quark pair forms a bound state by emitting soft gluons to be a color singlet. The first step is a hard process which can be described by pQCD, while the second step is a soft process which requires a model. In heavy-ion collisions, the first step is modified by so-called cold nuclear matter effects such as shadowing effects and the second step by hot nuclear matter effects. In this study we try to describe quarkonium production not only in pp collisions but also in heavy-ion collisions by using the same Remler’s formalism.

The probability that a quarkonium eigenstate $i$ with momentum ${\bf P}$ is produced at ${\bf R}$ is given by Remler:1975fm ; Song:2023zma

where ${\bf R}=({\bf r_{1}}+{\bf r_{2}})/2$, ${\bf r}={\bf r_{1}}-{\bf r_{2}}$, ${\bf P}={\bf p_{1}}-{\bf p_{2}}$, ${\bf p}=({\bf p_{1}}-{\bf p_{2}})/2$ with $\bf r_{i}$ and $\bf p_{k}$ being respectively the position and momentum of $k$’th (anti)heavy quark; . $\Phi^{W}_{i}({\bf r},{\bf p})$ is the Wigner density of quarkonium wavefunction which is nonvanishing only for flavorless and colorless combination; $W^{(N)}({\bf r_{1},p_{1},r_{2},p_{2},...,r_{N},p_{N}})$ is density matrix in Wigner representation of the $N$ (anti)heavy quarks produced in a proton-proton or a heavy-ion collision, which is simply taken for the classical phase space density distribution Song:2023zma . The summation in the right hand side implies that all possible combinations of heavy quark pairs are considered to form quarkonium state $i$.

In pp collisions the momentum of heavy quark and that of heavy antiquark are generated by the PYTHIA event generator followed by rescaling the transverse momentum and rapidity to imitate the FONLL calculations Song:2015sfa . The spatial separation between heavy quark and heavy antiquark is given by the uncertainty principle Song:2023zma .

Fig. 1 shows the $p_{T}$ spectra and rapidity distributions of $\Upsilon(1S)$, $\Upsilon(2S)$ and $\Upsilon(3S)$ including the feed-out from excited states in pp collisions at $\sqrt{s}=5.02$ TeV. One can see that the results are consistent with the experimental data from CMS and ALICE Collaborations CMS:2018zza ; ALICE:2021qlw

A hot dense matter produced in relativistic heavy-ion collisions delays the formation of quarkonium. We assume that the free energy of heavy quark pair is the heavy quark potential Kaczmarek:2003ph . In this case the dissociation temperature of $\Upsilon(1S)$ is roughly 3 $T_{c}$, while those of excited states are near $T_{c}$ Song:2023zma . The first Wigner projection of Eq. (1) is carried out when the local temperature of heavy quark pair is below the dissociation temperature of quarkonium state $i$. After that, whenever scattering happens to heavy quark or heavy antiquark, the Wigner production is updated Song:2023ywt , considering temperature-dependent quarkonium radius, which is obtained by solviing Schrödinger equation with the heavy quark potential.

Th left panel of Fig. 2 displays $R_{AA}$ of $\Upsilon(1S)$ and $\Upsilon(1S)$ as a function of transverse momentum in Pb+Pb collisions at $\sqrt{s}=5.02$ TeV. The upper limit of the magenta band in the left panel indicates $R_{AA}$ of $\Upsilon(1S)$ at the dissociation temperature of $\Upsilon(1S)$ and the lower limit of the band $R_{AA}$ of $\Upsilon(1S)$ at $T_{c}$. One can see that Wigner density of $\Upsilon(1S)$ decreases with time, because bottom and antibottom quarks are more and more separated from each other. The red solid line indicates $R_{AA}$ of $\Upsilon(1S)$, assuming only 10 % of (anti)bottom quark scatterings update the Wigner density of $\Upsilon(1S)$, which is consistent with the experimental data from the CMS Collaboration CMS:2018zza . In this case the interaction rate of $\Upsilon(1S)$ is roughly estimated to be 40-100 MeV at T=0.2-0.4 GeV Song:2023zma ; Andronic:2024oxz . As for $\Upsilon(2S)$ its dissociation temperature is close to $T_{c}$. That is why the blue band is so narrow that it looks like a line. The right panel shows $R_{AA}$ as a function of the number of participants. The solid lines include all contributions, while the dashed lines only diagonal contribution which originates from initial bottom quark pairs. One can see that the off-diagonal contribution which originate from two different initial bottom quark pairs is little even in central Pb+Pb collisions at $\sqrt{s}=5.02$ TeV, because the number of produced bottom quark pairs is not so large.

Normally quarkonium production in heavy-ion collisions is separately described from the production in pp collisions, though the nuclear modification factor ($R_{AA}$) is defined as their ratio. In this study we have used an unified method, the Remler’s formalism, for quarkonium production in pp as well as heavy-ion collisions. Assuming that the interaction rate of bottom (anti)quark bound in quarkonium is 10 % that of unbound bottom (anti)quark, our results well describe the experimental data in Pb+Pb collisions at $\sqrt{s}=5.02$ TeV. Furthemore, we have found that the off-diagonal contribution is little even in LHC.

We acknowledge support by the Deutsche Forschungsgemeinschaft (DFG, German Research Foundation) through the grant CRC-TR 211 ’Strong-interaction matter under extreme conditions’ - Project number 315477589 - TRR 211. This work is also supported by the European Unionâs Horizon 2020 research and innovation program under grant agreement No 824093 (STRONG-2020). The computational resources have been provided by the LOEWE-Center for Scientific Computing and the "Green Cube" at GSI, Darmstadt, and by the Center for Scientific Computing (CSC) of the Goethe University, Frankfurt.