FERMILAB-PUB-24-0005-T-V,IQuS@UW-21-072,IFJPAN-IV-2024-4
Comparative Study of Quarkonium Transport in Hot QCD Matter
Abstract
This document summarizes the efforts of the EMMI Rapid Reaction Task Force on “Suppression and (re)generation of quarkonium in heavy-ion collisions at the LHC”, centered around their 2019 and 2022 meetings. It provides a review of existing experimental results and theoretical approaches, including lattice QCD calculations and semiclassical and quantum approaches for the dynamical evolution of quarkonia in the quark-gluon plasma as probed in high-energy heavy-ion collisions. The key ingredients of the transport models are itemized to facilitate comparisons of calculated quantities such as reaction rates, binding energies, and nuclear modification factors. A diagnostic assessment of the various results is attempted and coupled with an outlook for the future.
1 Introduction
High-energy nucleus-nucleus collisions, such as those at the Large Hadron Collider (LHC) and Relativistic Heavy-Ion Collider (RHIC), provide the unique opportunity to create in the laboratory the conditions of the microsecond-old Universe, in particular the color-deconfined state called the Quark-Gluon Plasma (QGP). The microscopic understanding of the remarkable properties of the QGP and its hadronization remains a key challenge in nuclear physics.
The quarkonium families, charmonia and bottomonia, play a critical role in this objective and have long been prominent observables used to probe the fundamental color force in the hot QCD medium. Theoretically, the vacuum heavy-quark (HQ) potential provides a well-calibrated starting point for the study of quarkonium interactions in medium, see Refs. [1, 2, 3, 4, 5] for recent reviews. In particular, the string term in the HQ potential characterizes the long-range nonperturbative part of the force and is associated with the confining property of QCD. It is expected to play a central role in the transition from hadronic to partonic degrees of freedom, and may well be responsible for the strongly coupled properties of the QGP as evidenced by its transport properties, up to temperatures of 2-3 times the pseudo-critical temperature, , of the QCD crossover transition [6].
Similar to the way quarkonium spectroscopy provides information about the QCD force in vacuum, a systematic investigation of the in-medium force must involve the study of different quarkonium states in matter, e.g., as they subsequently dissolve with increasing temperature. The complexity in describing the in-medium properties of quarkonia and their implementation into transport calculations in ultra-relativistic heavy-ion collisions (URIHCs) prevents their use as a straightforward thermometer of the medium produced in these reactions. On the contrary, using available information on the space-time and temperature evolution in URHICs from other sources (e.g., hydrodynamic simulations or electromagnetic radiation), one can utilize quarkonium observables to deduce their in-medium properties and infer the fundamental interactions in QCD matter. In the vacuum, the 1S ground-state bottomonia ( and ) are small enough in size to be sensitive to the color-Coulomb part of the HQ potential. However, excited bottomonia and all charmonia are predominantly bound by the non-perturbative string term, which is characterized by a linearly increasing potential with a “string tension” 1 GeV/fm (and/or residual mesonic forces). Thus, charmonia and excited bottomonia are excellent probes of the in-medium confining force, as originally envisioned for the meson [7].
In the cooling of the expanding fireball, quarkonia can also be (re-)generated through recombination of individual heavy quarks and antiquarks diffusing through the medium. In particular, quarkonium formation can also occur through quarks and antiquarks from different initial pairs. This mechanism [8, 9, 10] has been shown to be critical for understanding the rise of production from RHIC to the LHC where (re)generation seems to constitute the main part of the yield observed in central Pb–Pb collisions [11]. The data are also compatible with production of exclusively through statistical hadronization at the QCD crossover phase boundary [12]. Precise measurements of the production cross section and the extraction of the charm-quark diffusion coefficient [13] are key objectives for the experimental program with heavy ions at the LHC in Runs 3 and 4 (2021-2029) [14] and will be important for drawing more definite conclusions. The measured elliptic flow, characterized by the flow coefficient , is significant for the mesons; it is consistent with transport model predictions at relatively low transverse momentum () and requires additional ingredients (such as space-momentum correlations of diffusing charm quarks) at higher [15]. The recent measurement of a rather large of mesons in high-multiplicity p-Pb events [16] came as a surprise and cannot be reproduced by transport model calculations [17], suggesting that initial-state effects could be the origin of the azimuthal correlations.
Regarding bottomonia, the current understanding suggests that (re)generation is less important for , but possibly figures as a major component in the strongly suppressed yield of excited states [18, 19]. However, the interplay of the reaction rates and in-medium binding energies, and thus the strength of the underlying HQ potential (i.e., its in-medium screening) remains a matter of debate [20, 21, 22, 23] and is also amenable to an interpretation based on comover interactions modeled with cross sections extracted from p-Pb data. It is therefore of great importance to obtain additional information about the typical time at which quarkonia are produced, in particular through spectra and elliptic flow, which contain information about the fireball collectivity imprinted on the quarkonia at the time of their decoupling. No significant elliptic flow has been observed experimentally for bottomonia [24, 25].
On the theoretical side, the basic objects are the quarkonium spectral functions which encode the information on the quarkonium binding energies and the (inelastic) reaction rates as well as melting temperatures. Steadily-improving constraints on the determination of the quarkonium spectral functions are available from thermal lattice QCD (lQCD), see Refs. [4, 26] for reviews. The information from the spectral functions can then be utilized in heavy-ion phenomenology via transport models. The latter provide the connection between first-principles information from lQCD and experiment that greatly benefits the extraction of robust information on the in-medium QCD force and its emergent transport properties, most notably the inelastic reaction rates of quarkonia. Thus far most transport models are based on rate equations and/or semiclassical Boltzmann equations [11]. In recent years transport approaches utilizing open-quantum system (OQS) frameworks have been developed in a Schrödinger-Langevin [27, 28, 29] or density-matrix [30, 31, 32] formulation, see also the review in Ref. [33]. These approaches can test the classical approximation underlying the Boltzmann and rate equation treatments and possibly quantify the corrections. Quantum effects may be particularly relevant at high in connection with finite formation times of quarkonia which are enhanced by the Lorentz time dilation in the moving frame; schematic treatments of this effect in semiclassical approaches suggest that varying formation times can lead to observable differences for high-momentum charmonia and bottomonia yields [34, 22, 18, 35, 21].
Finally, the implementation of phase-space distributions of heavy quarks diffusing through the QGP into quarkonium transport is being investigated by a few groups [36, 37, 38], which, as mentioned above, has direct impact on the magnitude and dependence of (re)generation processes [39, 40]. In particular, the role of non-perturbative effects in the HQ interactions in the QGP (which are believed to be essential to explain the large elliptic flow observed for D mesons) [13] needs to be accounted for; the associated large scattering widths are likely to require quarkonium transport implementations beyond semiclassical (or perturbative) approximations, which reiterates the motivation for quantum treatments of recombination processes.
The larger experimental data samples in Runs 3 and 4 at the LHC, combined with improved detector performance and measurement techniques, will lead to significant improvements over the current measurements, with extended kinematic coverage (in ) and the possibility to reach currently unobserved quarkonium states [14]. But even at present, the uncertainties in the quarkonium data are in some cases significantly smaller than model uncertainties. Improvements are needed in the conceptual aspects discussed above, as well as on the input parameters of the models. Clearly, the modeling of quarkonium production has reached a stage where the complexity of the problem, in connection with the precision reached by experiment, can no longer be handled by individual group efforts. A broadly-vetted consensus on the future path can only be achieved in a setting such as the one offered by a Rapid Reaction Task Force (RRTF) within the ExtreMe Matter Institute (EMMI) at GSI Darmstadt, which allows for much more in-depth discussions and working sessions, with the main protagonists from around the world in one room. We considered it critical to initiate a broad initiative now, to enable the theoretical progress that is needed for the interpretation of Runs 3 and 4 data at the LHC.
Based on the physics case outlined above, the most important issues where consensus and progress is most urgent were identified to be: i) the identification and model comparisons of transport parameters; ii) the controlled implementation of constraints from lattice QCD, and iii) the significance of quantum transport treatments. More concretely, these were addressed guided by the following questions:
1) To what extent are the currently employed transport approaches (mostly carried out in semiclassical approximations) consistent in their treatment of quarkonium dissociation and regeneration?
2) What are the equilibrium limits of the transport approaches and how do the former compare to the results of the statistical hadronization model?
3) What is the significance of the effects on quantum transport of the quarkonium wave packets, and what is needed to develop quantum transport into a realistic phenomenology?
4) How can the abundant information from lattice QCD (quarkonium correlation functions, heavy-quark free energies and susceptibilities, and the open heavy-flavor sector) be systematically implemented into transport approaches?
5) What are the ultimate model uncertainties, and will those allow conclusions regarding the fundamental question of the existence of hadronic correlations in a deconfined medium?
It is clear that the major experimental effort placed on measuring quarkonium with increasing precision needs to be matched by the effort on the theory front. Our RRTF was intended to stimulate a coordinated pooling of resources and lead the way towards a “unified” extraction of the fundamental QCD heavy-quark potential. We also reiterate that the modeling of heavy-quark diffusion in quarkonium production in a deconfined medium has important implications for the extraction of the diffusion coefficients of heavy quarks in QGP [13, 26].
Our article is structured as follows: we start out with a brief overview of the experimental status of quarkonium measurements at RHIC and the LHC in Sec. 2 and then collect some basic theoretical concepts for both the description of heavy quarkonia in medium as well as their implementation into transport models for heavy-ion collisions in Sec. 3. In Sec. 4 a rather detailed yet compact account of all transport approaches employed in the present effort is given; this is organized into a total of 19 categories (addressing model components and ingredients) that each research group was asked to elaborate on within their approach. In Sec. 5 we compare and discuss the numerical results of targeted calculations that all transport models were asked to carry out as far as their approach allows. We summarize and conclude our effort in Sec. 6.
2 Experimental overview
Quarkonium yields and the nuclear modification factor have been measured in heavy-ion collisions (including p/d-A collisions, which we will not discuss here) at the LHC and RHIC as a function of centrality and transverse momentum for both the charm and the bottom sector, both for ground states and for radially-excited states (the meson has not been measured yet in Au-Au collisions at RHIC). We briefly outline the status of quarkonium measurements in nucleus-nucleus collisions, with emphasis on the RHIC and LHC data.
2.1 Charmonium
After the studies of charmonium production at the SPS ( GeV) [41, 42] remained somewhat inconclusive, there was a question as to whether or not the observed suppression pattern of charmonium states is to be related to QGP effects. Since then, measurements at RHIC ( GeV) revealed a strong suppression of production in Au-Au collisions compared to pp [43, 44, 45, 46], of a similar magnitude as measured at SPS and stronger at forward rapidity than at midrapidity.
In U-U collisions at RHIC, the production yield of was measured to be higher than in Au-Au collisions for the same number pf participating nucleons, values, in central collisions [48]. Another observation at RHIC was that the suppression is similar at and 62.4 GeV compared to 200 GeV [49, 50]. A strong suppression was observed in Cu-Au collisions [51], while no suppression was measured in Cu-Cu [45], albeit with large uncertainties.
A clearly different behavior was measured in Pb-Pb collisions at the LHC, where a significantly reduced suppression was observed (larger ) [52, 53]. The current LHC data [54] are compared to the RHIC data in Fig. 1 as a function of centrality, expressed as the charged-particle pseudorapidity density (at ), a proxy for the energy density of the fireball. As at RHIC, the values are larger at midrapidity than at forward rapidity.
Another marked difference between the LHC and RHIC data can be observed in Fig. 2. While the RHIC data exhibit no significant dependence, the LHC data show significantly larger values at low , even with a hint at exceeding unity, at midrapidity. The data are compared to those of the D mesons and pions for central Pb-Pb collisions in the the left panel of Fig. 3. The observed ordering is a fingerprint of the different production mechanisms of the charm and light quarks and suggests as well a (re)generation mechanism for the mesons in QGP and/or at hadronization (QCD crossover phase transition).
While the elliptic flow of was found to be compatible with zero at RHIC (albeit with large uncertainties) [55], it is large at the LHC [56, 57]. This is shown in the right panel of Fig. 3, where the data are compared to those of D mesons and pions for mid-central collisions.
The relative production of and mesons as a function of collision energy is shown in Fig. 4. The data points for pp collisions are from experiments at the SPS [58], HERA (pA) [59], RHIC [60, 61], and the LHC [62, 63, 64]. The average value of the pp measurements is represented by the black horizontal line with the 1-uncertainty represented by the dashed lines. The point for central Pb–Pb collisions at SPS energy is from the NA50 experiment [65] and the point at the LHC is from ALICE [66].
2.2 Bottomonium
The early observation by the CMS collaboration at the LHC of the suppression of mesons [67], which was clearly stronger for the radially-excited states, was followed by the current precise measurements of for the states by CMS [68, 69, 70] and ATLAS [71] at midrapidity and ALICE at forward rapidity [72]. A similar and suppression pattern (and magnitude, for ) was measured at RHIC by STAR [73]. A summary of the data is presented in Fig. 5.
The suppression shows a similar magnitude at forward rapidity and midrapidity. It is gradually stronger going from the 1S to the 2S and 3S states; this significant pattern is denoted as “sequential suppression”, with the picture of the melting of the excited states in QGP and their ”missing” feed-down to the 1S state.
3 Theoretical background
One of the advantages of quarkonium physics in QCD matter is the availability of first-principle lQCD computations for a variety of finite-temperature quantities, such as quarkonium correlation functions, heavy-quark free energies and susceptibilities. This is particularly relevant for a strongly coupled medium where non-perturbative effects are expected to be prevalent and perturbative calculations must be interpreted with care. However, since lQCD computations are carried out in Euclidean space-time, the information is usually not readily applicable for use in transport simulations of heavy-ion collisions.
In Sec. 3.1, we discuss some of the methods that have been used to assess in-medium quarkonium properties and their implementation into transport models, and in Sec. 3.2 we give a more detailed account of pertinent lQCD results.
3.1 Theoretical Methods
A key concept in facilitating the connection between lQCD and effective-model calculations are in-medium spectral functions. Quite generally, the latter encode information on the degrees of freedom in the medium (through the presence or absence of well-defined peaks), their in-medium masses and (if applicable) binding energies and the elastic and inelastic reaction rates represented by their widths. Both binding energies and reaction rates are related to the in-medium interaction between heavy quarks and between the heavy quarks and the thermal partons of the medium, respectively. This connection also highlights the importance of utilizing a reasonably realistic model of the thermal medium in which the heavy quarks and quarkonia are embedded.
In practice one, therefore, needs effective descriptions of the in-medium quarkonium dynamics. In the present context, this brings in another asset of the HQ sector, namely the use of 1/ expansion schemes, such as HQ effective theory (which amounts to neglecting the antiparticle components of the Dirac spinors), non-relativistic QCD (NRQCD), and potential NRQCD (pNRQCD). The latter, in particular, is based on neglecting the energy transfer in HQ scatterings, which allows one to convert the 4D Bethe-Salpeter equation into 3D Lippmann-Schwinger type equations (or, in coordinate space, a Schrödinger equation), rendering the problem much more tractable. One further needs appropriate many-body approaches to implement the in-medium physics. Among the ones being used are: (i) perturbative Lagrangian-based hard-thermal loop (HTL) perturbative approaches, (ii) the dynamical quasi-particle model (DQPM) which utilizes a running coupling and accounts for parton spectral functions which can describe the lQCD equation of state (EoS) down to temperatures close to [74], and (iii) the thermodynamic -matrix, which is a non-perturbative Hamiltonian-based approach that solves the 1-and 2-body correlation functions self-consistently based on a unified 2-body potential for both heavy and light-parton interactions [75, 76].
To make contact with experiment, the equilibrium physics of quarkonia has to be implemented in transport descriptions. This has traditionally been done in semiclassical models, including the comover interaction model (CIM) [77], simulations of the Boltzmann equation [78, 79, 80] and kinetic rate equations [9, 81, 82]) based on dissociation rates in QGP and/or hadronic matter as well as suitably computed equilibrium limits to account for regeneration processes. Constraints from lQCD have been considered by, e.g., using the results for the HQ free or internal energies as underlying potentials to compute in-medium binding energies [34, 83, 21], and to test the masses and widths of charmonium spectral functions via pertinent Euclidean correlator ratios [34].
Due to the small energy scale (or long time scale) provided by the in-medium quarkonium binding energies, especially at temperatures in the vicinity of the quarkonium “dissociation temperature”, several efforts are underway to develop quantum transport equations based on the framework of open-quantum systems. Most of them utilize a classification of regimes where the binding energy, relative to the temperature, is either small (where quantum effects are expected to be relevant) or large (where semiclassical approximations are reliable); these are further combined with scale hierarchies from HQ effective field theories. Information from lQCD can be implemented through transport coefficients, such as the HQ diffusion coefficient, albeit without explicit 3-momentum dependencies (which will be challenging to obtain from lQCD). The question of regeneration, especially in the case of multiple quark pairs, has not been addressed yet at a level suitable for phenomenological applications. The equilibration of a single pair has been studied in Ref. [84] and in semiclassical limits for multiple pairs in Refs. [85, 86]. All of these considerations will figure in the model descriptions given in Sec. 4.
3.2 Lattice QCD Results
Lattice-QCD calculations can provide first-principles input into theoretical modeling of quarkonium production in heavy-ion collisions. Many quantities of interest, like in-medium quarkonium masses and widths, or transport coefficients, are encoded in the spectral functions, defined as the imaginary part of the retarded meson correlation functions [87, 4]. For example, the in-medium widths of quarkonia are closely related to the reaction rates used in transport models. If the widths are reasonably small, quarkonium states can be identified by peaks in the spectral functions. As temperature increases, the peaks become broader and ultimately can no longer be used to identify quarkonium states. For example, if the width of the peak is much larger than the energy splitting between different quarkonium states it is no longer possible to extract well defined quarkonium states. Obtaining the spectral functions from lQCD is challenging because the latter is formulated in Euclidean time, and the correlation functions are given in terms of integrals over spectral functions. Temporal correlation functions are related to spectral functions via a Laplace transformation, while spatial correlation functions are related to spectral functions via a double integral transformation [4]. Lattice QCD calculations can also be combined with EFT approaches. For example, information about quarkonium spectral function can be obtained using a lattice formulation of NRQCD [88, 89, 90, 91, 92, 93]. In this way one avoids large discretization effects due to HQ masses An additional benefit arises from the fact that meson correlators in NRQCD do not obey periodic boundary conditions, which effectively implies that information on meson correlators can be obtained from doubling the temporal extent in the Euclidean time direction. Heavy-Quark Effective Theory (HQET) can be used for lQCD calculations of the HQ diffusion coefficient [94, 95, 96].
Most lattice studies of quarkonium spectral functions use point meson operators, i.e., meson operators with the quark and antiquark field located at the same spatial point. It turns out that temporal meson correlation functions with point meson operators have limited sensitivity to the in-medium properties of quarkonia [1, 4]. This is due to large contributions from the continuum part of the spectral function to the correlators of point meson operators, as well as the rather small temporal lattice extent at high temperatures [97, 98]. Therefore, no conclusive results on the in-medium properties of the quarkonium states could be obtained from the temporal correlation functions of point meson operators. There only seems to be a consensus that the 1S bottomonium state can survive in the QGP for MeV, with a small mass shift [89, 90, 99]. The study of spatial meson correlation functions is not limited to small separations, rendering them more sensitive to the in-medium properties of quarkonia [100, 101, 102]. In particular, indications were found that 1S charmonia states may dissolve for temperatures of 200-300 MeV [100, 101], while 1S bottomonium states will dissolve for temperatures above 500 MeV [102]. The latter finding is consistent with the analysis of bottomonium spectral function in lattice NRQCD.
Correlators of extended meson operators, i.e., meson operators with quark and antiquark fields separated by some spatial distance, are more sensitive to the in-medium modification of quarkonia, since the contribution of the continuum part to the spectral function is reduced. Using NRQCD with extended operators it was possible to analyze in-medium masses and widths of different bottomonium states [91, 92]. Interestingly, it was found that the in-medium mass shift of all bottomonium states is small and compatible with zero within estimated errors. The in-medium width of different bottomonium states was found to increase with temperature, and that the magnitude of the width follows a hierarchy in the sizes of the different states [91, 92].
The in-medium modification of interactions in QGP has been traditionally studied in terms of the free energy and singlet free energy of a static pair. The latter quantity can be defined in Coulomb gauge. State-of-the-art calculations in 2+1 flavor QCD with physical quark masses suggest that color screening in the free energy sets in at distances [103]. Previous studies of the free energy for two [104] and three [105] degenerate quark flavors with unphysical masses have been used as input potentials in some phenomenological models.
The free energy characterizes the interactions at time scales much larger than the inverse temperature. For quarkonia physics, it is more relevant to consider a complex potential defined in terms of Wilson loops [106]. The first calculation of the complex potential along these lines with 2+1 flavor QCD with unphysical quark masses found [107] that the real part of the potential is screened. A parametrization of these results using a generalized Gauss law model [108] has also been used in some phenomenological models. The corresponding lattice calculations are performed on lattices and with limited statistics. Another lattice study that also uses lattices, but with much larger statistics, extracted a different result [109]. Here, simple but physically motivated parametrizations of the spectral functions were used to obtain the real and imaginary parts of the potential: the real part of the potential turns out not to be screened in general [109]. The only way to obtain a screened potential from the lattice results of Ref. [109] is to use a perturbative HTL inspired representation of the spectral function [109]. However, it turns out that, although HTL results for the equation of state of in good agreement with lattice calculations at high temperatures [110, 111], HTL calculations are in disagreement with the lattice results on the Wilson line correlators even for very high temperatures [109]. More recently, the analysis of Wilson line correlators from Ref. [109] has been repeated using lattices with much larger temporal extent, , and it was again found that the real part of the potential is not screened [112]. Very recently, a microscopic calculation of the Wilson line correlators has been carried out in the thermodynamic -matrix approach [113]; it was found that, with a refined input potential which also shows little screening at small and intermediate distances, the lattice results can be reasonably well reproduced, but with very large quarkonium widths.
The HQ momentum diffusion coefficient, , is also used as an input to a few phenomenological models of quarkonia production. It can be calculated on the lattice using HQET, as mentioned above, in terms of a chromo-electric correlation function with fundamental Wilson lines [94]. Most lQCD calculations of have been performed in quenched QCD [114, 115, 116, 117, 118, 119]. Very recently the first lattice calculation of the HQ diffusion constant in 2+1 flavor QCD appeared [120], which includes the HQ mass dependence of and the related HQ spatial diffusion coefficient, [121]. Within the pNRQCD approach one can consider a chromo-electric correlator with adjoint Wilson lines. This correlation function defines another transport coefficient , which enters the pNRQCD-based open quantum system approach of quarkonium production [122]. To leading order in perturbation theory, , but this is not true in general. Unfortunately, no lattice QCD calculations of are available so far.
As will be discussed below, phenomenological models of quarkonium production in heavy-ion collisions require knowledge of the in-medium HQ masses. Generalized charm and bottom susceptibilities can be useful for constraining these masses. Lattice calculations of these quantities exist both in quenched [98] and full QCD [123, 124, 125, 126, 127]. It was inferred that the in-medium charm- and bottom-quark masses decrease with increasing temperature [98, 123, 121].
4 Model Descriptions
In this section we provide a synopsis of each of the phenomenological approaches to describe heavy-quarkonium transport in URHICs that has been included in the present effort. By summarizing each approach, guided by the same list of key ingredients, assumptions and inputs, we attempt to lay out the strengths and areas of future improvement in the various model calculations.
4.1 The Duke-MIT approach
The Duke-MIT approach is based on a set of coupled Boltzmann equations that describe the in-medium evolution of both heavy quarks and quarkonia [128]. The Boltzmann equation for quarkonia describes their dissociation and recombination. It is derived by using the pNRQCD effective field theory in the hierarchy of energy scales
(1) |
and the open quantum system approach in the quantum optical limit [129, 130, 131]. The symbol means that the temperature is on the order of or smaller than the binding energy; e.g., is still in the validity region of the framework. The dissociation and recombination processes are formulated in a factorized way such that their collision terms can be written as a convolution of a dipole transition term of the heavy quark-antiquark wave function and a generalized (chromoelectric) gluon distribution function of the QGP [130]. The generalized distribution function (GGD) is a chromoelectric field correlator defined in a gauge-invariant and non-perturbative way, and it naturally generalizes the gluo-induced and inelastic scattering processes to account for the situation where the medium partons inducing these processes can be off-shell and non-perturbative. When , the GGD has been calculated to next-to-leading order accuracy [132]. When , the GGD has been calculated by holography [133, 134] and can be studied by Euclidean lattice QCD calculations [135]. The consistency between the dissociation and recombination implementations has been tested in a static QGP box [37]. The bulk dynamics of the QGP is described by a 2+1D viscous hydrodynamic equation that has been implemented in the “VISHNU” package and calibrated against the production of light particles [136]. The Boltzmann equation for open heavy quarks describes their diffusion and energy loss in the QGP, which has been implemented in the “Lido” package and has been calibrated against open heavy-flavor production observables [137]. The effect of nPDFs is taken into account by using the EPPS16 parameterization [138] and the feed-down network is considered in the hadronic stage. The framework has been used in phenomenological studies of bottomonium production in heavy ion collisions [128].
Model ingredients
-
-
In-medium potential: An unscreened Coulomb potential is used as the in-medium potential.
-
-
Vacuum limit of potential/spectroscopy: The Duke-MIT approach reproduces the (1S) mass in vacuum by construction but not the 2S and 3S excited states, which is a consequence of using a Coulomb potential.
-
-
Reaction rates: For dissociation, both gluo-dissociation and inelastic scattering contributions are included in a pNRQCD based calculation. The pNRQCD calculation takes into account the interference between the heavy quark scattering and the heavy antiquark scattering. Therefore, it goes beyond the quasi-free approximation. For recombination, similar processes (gluo-induced and inelastic) are included [139]. For unbound heavy quark-antiquark pairs, the diffusion and energy loss are included independently for each heavy quark.
Three-momentum dependence comes in since the reaction rates are calculated in the rest frame of the heavy quark-antiquark pair for dissociation and recombination and boosted back to the laboratory frame [140]. The rest frame is a natural frame for the non-relativistic treatment in the effective field theory. However, for sufficiently large transverse momenta, the non-relativistic treatment breaks down in the rest frame since the partons from the medium get significantly boosted.
-
-
Assumptions about the medium (degrees of freedom, etc.): The calculations of the gluo-induced and inelastic processes mentioned above assume the medium consists of free quarks and gluons. However, the Boltzmann equation for quarkonium formulated in our framework can easily go beyond this. The nonperturbative nature of the medium partons relevant for quarkonium is encoded in terms of the chromoelectric field correlator mentioned above. It has been shown that, beyond leading-order in the strong coupling, this chromoelectric field correlator for quarkonium is different from the one used to define the heavy quark diffusion coefficient [132, 141].
-
-
Temperature dependence of heavy quark masses: The heavy quark mass is treated as independent of temperature here.
-
-
Equilibrium limits in transport: The heavy quark Boltzmann equation can drive the system of unbound heavy quarks to kinematic equilibrium. The quarkonium Boltzmann equation can drive the system of bound and unbound pairs into chemical equilibrium [37].
-
-
Constraints from lattice QCD: At the moment, there is no constraint from lattice QCD calculations. However, the next systematic step is to apply lattice QCD techniques to calculate the chromoelectric field correlator mentioned above [135].
-
-
Range of applicability and how is this range established: The range of applicability is established by the power counting of the effective field theory pNRQCD. When the temperature is much bigger than the inverse of a typical quarkonium size, quarkonium cannot exist as a well-defined bound state inside the medium, so the dynamics are effectively described in terms of an unbound but possibly correlated pair. When the temperature is much smaller than the inverse of the typical quarkonium size, the effective dynamics is described as quarkonium dissociation and recombination and these processes can be systematically calculated in pNRQCD. When the medium becomes nonperturbative, a nonperturbative determination is needed. Due to the use of the optical limit there is also a upper limit on the temperature since the optical limit assumes that the binding energy of states is small compared to .
-
-
Quantum features: The Boltzmann equation for quarkonium is derived from the open quantum system in the quantum optical limit, with a systematic semi-classical (gradient) expansion. Quantum corrections to the recombination term have been worked out but not yet implemented. There is no quantum correction to the dissociation term [130, 131].
-
-
Regeneration: The recombination is formulated in the open quantum system approach (quantum optical limit) and thus it is consistent with the description of dissociation [129, 130, 131]. The recombination term in the Boltzmann equation is calculated in pNRQCD, which is valid when the quarkonium size is small, compared to the medium temperature. When the temperature is much bigger, the effective dynamics are described in terms of an unbound but possibly correlated heavy quark-antiquark pair, which is consistent with the quantum Brownian motion limit of the open quantum system framework.
-
-
Coupling to open heavy-flavor sector: The open heavy-flavor transport is an important part of the coupled Boltzmann equations. It provides the heavy quark-antiquark distribution that is needed as an input in the recombination term, which allows the study of both correlated and uncorrelated recombination processes. The open heavy-flavor transport is necessary for the system to reach kinetic equilibrium.
-
-
Hadronic-Phase Transport: No hadronic-phase transport except for feed down.
- -
-
-
Are cold nuclear matter effects, nPDF effects, etc. taken into account?: nPDF parameterized by EPPS16.
-
-
Constraints from pA and dA collisions: EPPS16 uses some pA data to fit the nPDF. We also use p-Au data from STAR to fix the cold nuclear matter effects in 200 GeV Au-Au collisions.
-
-
Medium evolution model: 2+1D viscous hydro (“VISHNU”) that has been calibrated against light particle production observables.
-
-
Feed down implementation: Feed down contributions from 2S, 3S, 1P, and 2P states for states.
-
-
Comparisons to experimental data: Can describe (nS) data from both LHC and RHIC collisions, except for the recent (3S) data from CMS since recombination of (3S) has not been included.
-
-
Phenomenological breadth: We included (nS) and in the reaction network. We studied both RHIC and LHC collision energies.
4.2 The Munich-KSU approach
We utilize a set of coupled quantum evolution equations for the reduced density matrix derived in Refs. [31, 122, 32] using the formalism of open quantum systems (OQS) and the effective field theory (EFT) potential nonrelativistic QCD (pNRQCD) [143, 144]. The resulting set of coupled evolution equations describes the time evolution of singlet and octet Coulombic states propagating in a thermal medium at temperature realizing the hierarchy of scales
(2) |
where is the heavy quark mass, is the Bohr radius of the quarkonium state, , is the temperature of the strongly-coupled QGP medium and is the Coulombic binding energy of the bound state. In this regime, the evolution equations take the form of a Lindblad equation
(3) |
where and are the quarkonium density matrix and in-medium Hamiltonian and the collapse operators encode interactions with the medium. The anticommutator term is responsible for the in-medium width. The first term in parentheses ensures that the evolution is trace preserving. For a detailed discussion, see Eqs. (2)-(5) and accompanying text of Ref. [19].
We work in the bottom sector, which realizes the dilute limit described by the above Lindblad equation linear in the quarkonium density . We, furthermore, work up to order in the finite temperature pNRQCD multipole expansion and implement an additional expansion to linear order in . The heavy quark mass is taken in the scheme: GeV where is taken from the Particle Data Group [145]. The Bohr radius is computed by solving its self-defining relation
(4) |
where is the quadratic Casimir of the fundamental representation of , is the strong coupling evaluated at the energy scale using the 1-loop, 3-flavor running coupling and MeV [145]. This gives . In the OQS+pNRQCD formalism, interactions with the medium are calculated systematically and occur as in-medium modifications to the heavy quarkonium potential (the real and imaginary parts of which cause a thermal mass shift and thermal width , respectively) and color singlet to octet, octet to octet and octet to singlet transitions (the latter implementing recombination of the heavy quark-antiquark pairs). For the hierarchy of scales realized in Eq. (2), the imaginary part of the in-medium potential is controlled by the heavy quarkonium momentum diffusion coefficient and the real part by , which is the dispersive counterpart of ; these transport coefficients are fixed from independent, unquenched lattice measurements of and [146].
Model ingredients
-
-
In-medium potential: We use a vacuum Coulomb potential plus a quadratic in-medium modification through the transport coefficient . The Coulombic binding energy MeV sets the scale of the spacing of the energy levels.
-
-
Vacuum limit of potential/spectroscopy: In the kinematical situation described by Eq. (2), the real part of the in-medium quarkonium potential is the sum of a Coulomb part, , and a term proportional to [147, 31]. Taking the vacuum limit , the in-medium potential reduces to the Coulomb potential plus a correction proportional to [143]. We emphasize that to accurately reproduce the vacuum spectrum (in particular fine and hyperfine splittings), higher-order loop corrections in to the Coulomb potential may be necessary to include (cf. Ref. [144]). In the adopted hierarchy of energy scales, they are parametrically subleading with respect to the considered in-medium quarkonium potential. The imaginary part of the in-medium quarkonium potential is proportional to [147, 31].
-
-
Reaction rates: Using the OQS+pNRQCD formalism, we systematically calculate interactions with the medium at strong coupling: our self-energy contributions induce a thermal shift in the quarkonium energy, a thermal quarkonium width and singlet-octet, octet-octet and octet-singlet transitions that describe recombination and account for gluodissociation [148] parton dissociation [149] and screening effects (in the short-distance limit). We note that the thermal decay width can be computed as either (1) an expansion in with the state of the art expression extending up to order or (2) using the eigenvalue of the non-Hermitian Lagrangian which emerges at a given order in the expansion. For a detailed discussion of method (1) and explicit expressions, see Sec. 4.1 of Ref. [32]. For this report we use method (2) since, in practice, it provides a better description of the dynamical evolution than using widths based on the vacuum states. We consider a quarkonium state comoving with the medium; the width is, therefore, independent of .
-
-
Assumptions about the medium (degrees of freedom, etc.): We assume a strongly-coupled QCD plasma in local thermal equilibrium. The degrees of freedom are singlet and octet quarkonium states. Interactions with in-medium gluons are integrated out and described using the pNRQCD Lagrangian. In a medium realizing the hierarchy of scales of Eq. (2), the relevant transport properties of the medium are encoded in the heavy quarkonium momentum diffusion coefficient and its dispersive counterpart .
-
-
Temperature dependence of heavy quark masses: The heavy quark mass is treated as independent of temperature. Similar effects are already included in the chromoelectric field correlator mentioned above [130]. We work in the scheme for the heavy-quark mass in which . It is important to note that the mass of the bottom quark is not an observable quantity and is, furthermore, a scheme-dependent quantity. In our framework, the heavy-quark mass is a parameter entering the evolution equations; it receives no medium corrections. The bound-state mass, on the other hand, is an observable quantity, and its thermal correction up to order in the multipole expansion of pNRQCD at finite temperature is given by .
-
-
Equilibrium limits in transport: We do not enforce and have not observed (on time scales relevant to HICs) an equilibrium limit; this is a subject of ongoing investigation.
-
-
Constraints from lattice QCD: The transport coefficients and , which fully encode the interaction of the quarkonium state with the medium, are calculated from lattice measurements of the thermal width [92, 109] and the thermal mass shift [109] of the . For a presentation of the method (with older lattice data), see Ref. [146].
-
-
Range of applicability and how is this range established: We do not extend the medium evolution to temperatures below MeV. As discussed in our executive summary, our evolution equations are derived assuming the hierarchy of scales ; furthermore, our in-medium width operator is computed order-by-order as an expansion in . At zeroth-order in this expansion [31, 122, 150, 151], we do not extend our analysis below MeV. Our state of the art results include terms up to and including order ; a phenomenological comparison against the zeroth-order results (see discussion at the end of sec. 4.2 of Ref. [32]) leads us to select MeV as a lower bound when including these terms. Finally, to ensure that is fulfilled at all times, we do not extend our analysis above MeV.
-
-
Quantum features: Our in-medium EFT description of the system is fully quantum. In practice, we solve the Lindblad equation by computing the real-time evolution of the Schrödinger equation for a large ensemble using the quantum trajectories algorithm. This provides a numerical solution of the evolution of the reduced density matrix that describes both singlet and octet states.
-
-
Regeneration: pNRQCD octet to singlet chromoelectric dipole transitions implement quantum regeneration. The quantum mechanical wave function includes both bound and free pairs and transitions between these two types of states.
-
-
Coupling to open heavy-flavor sector: We work in the dilute limit and thus only in the bottom sector; open bottom states are implemented via the octet (unbound) contribution.
-
-
Hadronic-Phase Transport: When the local temperature of the medium falls below , we terminate the coupling with the medium and evolve the state in the vacuum; in practice, this corresponds to termination of the evolution. We apply our feed-down procedure to the resulting survival probabilities and do not implement hadronic-phase transport.
-
-
Initial quark/quarkonium distributions: As the heavy quark mass , which sets the scale of heavy quarkonium formation, is the largest scale of the combined system, our initial state is point-like with respect to the medium. We thus utilize a peaked Gaussian as a numerically tractable approximation of a position-space delta function. We sample the initial positions of these quantum mechanical wave-packets from the AA binary collision profile for a given impact parameter and sample their initial transverse momentum from a distribution where is the average mass of a bottomonium state.
-
-
Are cold nuclear matter effects, nPDF effects, etc. taken into account?: We do not include cold nuclear matter or nPDF effects.
-
-
Constraints from pA and dA collisions: We do not include constraints from pA and dA collisions.
-
-
Medium evolution model: We make use of a 3+1D dissipative hydrodynamics code based on the relativistic quasiparticle anisotropic hydrodynamics (aHydroQP) framework [152, 153, 154]. The aHydroQP framework has been tuned to reproduce a large set of experimental soft-hadronic observables, such as the total charged-hadron multiplicity, identified hadron spectra, integrated and identified hadron elliptic flow, and HBT radii at both RHIC and LHC nucleus-nucleus collision energies [155, 156, 157, 158, 159, 160]. The evolution uses a realistic equation of state determined from lattice QCD measurements [161] and self-consistently computed second- and higher-order dissipative transport coefficients [154].
-
-
Feed down implementation: We implement a feed-down procedure on our results for the survival probability before comparing against experimental measurements. The feed-down procedure consists of multiplying a vector containing the survival probabilities by a feed-down matrix, the entries of which are determined from the branching fractions of the considered states obtained from the Particle Data Group listings. For a more detailed description, including the states considered and the precise values of the entries of the feed down matrix, see sec. 6.4 of Ref. [150].
-
-
Comparisons to experimental data: We compare with experimental measurements of the nuclear modification factor of the , and states as functions of the number of participants and of transverse momentum; double ratios of the aforementioned nuclear modification factors; and the elliptic flow of the aforementioned states from the ALICE [72], ATLAS [162] and CMS [69, 163] collaborations.
-
-
Phenomenological breadth: We work exclusively in the bottom sector and compare with experimental data for the , and states obtained in Pb-Pb collisions at TeV.
4.3 The Nantes model
The Nantes model is based on earlier work by Blaizot and Escobedo [85]. Using NRQCD and the open quantum systems framework, coupled quantum master equations are derived in the quantum Brownian motion regime, describing the dynamics of a single quark-antiquark pair in the QGP. Our work extended their equations in order to ensure that the positivity of the density operator is preserved [164, 165]. The static screening in the medium and the dynamical processes occurring between the pair and the medium constituents are described by a complex potential. The quantum master equations (QMEs) derived take into account both diffusive and dissipative effects in the medium. The dissociation and diagonal recombination are described by the singlet octet color transitions of the pair. All together, this results in a quantum master equation of the Lindblad type expressed as
(5) |
where is the density operator containing a singlet and an octet component. and represent the kinetic energy and the (screened) real potential, while , , and respectively account for the fluctuations, the dissipation and the preservation of positivity; they depend on the assumed imaginary part of the complex potential.
The off-diagonal recombination is not treated as only one pair is considered in our modeling. The medium is assumed to have a uniform temperature, which is either fixed during the whole evolution or is dynamical. The dynamical medium follows an average temperature profile obtained from EPOS4, corresponding to the average temperature along the path of multiple pairs. The resolution of the quantum master equations is performed in one dimension. We have developed a new complex potential, specifically tailored to model the three-dimensional phenomenology [166]. In addition, to obtain tractable equations for the evolution of a pair, the center of mass degrees of freedom were integrated out.
Model ingredients
- -
-
-
Vacuum limit of potential/spectroscopy: In the vacuum, the real part of the potential is linear with a saturation corresponding to twice the mass of the lightest D or B meson. The mass spectra of charmonium and bottomonium states are reproduced with a precision of 20 MeV.
-
-
Reaction rates: The reaction rates are obtained from an imaginary potential, which is one-dimensional and fitted to a three-dimensional potential [108]. Each term , , and brings it own contribution, with increasing power of the relative velocity. The dominant one stemming from the term simply leads to .
-
-
Assumptions about the medium (degrees of freedom, etc.): The medium constituents are light quarks and gluons. The medium is assumed to be in equilibrium.
-
-
Temperature dependence of heavy quark masses: The temperature dependence of quark mass is taken as , where is the bare mass – taken as 1.47 GeV for c quarks and 4.88 GeV for b quarks and is the asymptotic value of the 1D real potential.
-
-
Equilibrium limits in transport: The QME (5) admits a unique asymptotic state for each temperature and such a state is reached as the final outcome of any initial state. For large relative distance, the complex potential saturates and the QME for the pair splits into two independent QMEs for each quark. In the semi-classical limit, these QMEs are equivalent to Fokker-Planck equations (with a momentum diffusion coefficient expressed as ) and naturally lead to Boltzmann distributions for individual quarks. In the quantum regime, the HQ effective temperature is 10% to 20% higher than the medium temperature. This is due to the UV regularization procedure adopted for the imaginary potential. When the evolution is performed in a box, the asymptotic spatial densities of both and are finite and some spatial correlation emerges at small relative distance. The corresponding asymptotic distribution of quarkonia is presently under investigation.
-
-
Constraints from lattice QCD: Lattice QCD is used to constrain the complex potential.
-
-
Range of applicability and how is this range established: The model is developed in the quantum Brownian regime due to the assumed scale hierarchy.
-
-
Quantum features: The treatment within the open quantum system framework is fully quantum (under the usual assumptions made in the Quantum Brownian regime).
-
-
Regeneration: “Diagonal” regeneration is included but not “off-diagonal” regeneration as only one pair is considered.
-
-
Coupling to open heavy-flavor sector: As explained previously, the dynamics naturally factorize when the distance between and is large. This has however little influence on the quarkonia probability measured at finite time since only pairs which do not evolve too far apart during the evolution contribute.
-
-
Hadronic-Phase Transport: No transport in the hadronic phase.
-
-
Initial quark/quarkonium distributions: The initial distribution of heavy quarks is taken from the EPOS4 framework. The ensuing transverse distribution of charmed hadrons is found in good agreement with the experimental distributions for the cases investigated up to now [167]. In AA collisions, the spatial distribution of heavy quarks and heavy quarks pairs follows a scaling. In our approach, there is no initial quarkonia per se but pairs. The initial internal degrees of freedom for these pairs are described as a pure-state density operator either in the singlet or in the octet representation. As for the relative distance dependence, compact Gaussian states are privileged, although it is also possible to initiate the evolution with any eigenstate from the vacuum or thermal potential.
-
-
Are cold nuclear matter effects, nPDF effects, etc. taken into account?: No Cold Nuclear Matter effects are included.
-
-
Constraints from pA and dA collisions: No constraints from pA/dA collisions.
-
-
Medium evolution model: Two cases are considered: a static medium with fixed temperature and an evolving medium following an average temperature profile from EPOS4.
-
-
Feed down implementation: No realistic feeddown is included.
-
-
Comparisons to experimental data: No comparison to experimental data for the moment.
-
-
Phenomenological breadth: Charmonium and bottomonium systems can be studied within our model. As our model is one-dimensional, we only consider even (S-like) and odd (P-like) eigenstates of the real part of the complex potential.
4.4 Parton-hadron string dynamics
Quarkonium production in the Parton-Hadron-String Dynamics (PHSD) is composed of the production of heavy quark pairs, the interaction and dynamics of heavy (anti)quarks in the medium and the Wigner projection of heavy quark pairs onto the quarkonium wave function based on the Remler formalism [168, 169, 170, 171, 86, 172].
First, the Wigner production is carried out in pp collisions without any nuclear matter effects. Heavy quarks and heavy antiquarks are separated from each other in space by Monte Carlo such that the average distance is on the order of the inverse mass and each momentum is provided by the PYTHIA event generator [173] after rescaling the rapidity and transverse momentum of heavy quark to mimic the FONLL calculations [174]. Then the Wigner projection finds out the most suitable radius of each quarkonium state from the comparison to the experimental data in pp collisions [175].
In heavy-ion collisions the (anti)shadowing effects are implemented using the EPS09 [176] which modifies heavy quark production and affects quarkonium production as well. Produced heavy quarks and antiquarks interact in the QGP within the Dynamical Quasi-Particle Model (DQPM) [74]. Quarkonium properties in the QGP, such as the dissociation temperature and radius are obtained by solving the Schrödinger equation with the heavy quark potential from lattice QCD. Once the local temperature is lower than the dissociation temperature, quarkonium begins to form through the Wigner projection with a temperature-dependent radius. Since each state of quarkonium has a different dissociation temperature, this projection takes place at different times and positions in heavy-ion collisions. We take into account 1S, 1P, and 2S states for charmonia and 1S, 1P, 2S, 2P, and 3S states for bottomonia and also feed down from excited states to a lower state as well as -hadron decay to charmonium.
Whenever a heavy (anti)quark scatters in QGP, a Wigner projection is carried out and the density of each quarkonium state is updated. If a heavy quark pair forms a bound state, its cross section with thermal partons will be smaller than the sum of the cross sections of heavy quarks and of heavy antiquarks due to the interference terms. Therefore, we introduce an ad-hoc suppression factor to the heavy quark scattering such that only part of the scatterings affects quarkonium production and dissociation. The suppression factor is roughly 10 % for bottomonia in Pb+Pb collisions at 5.02 TeV, which produces results consistent with experimental data [177].
Model ingredients
-
-
In-medium potential: The present form of PHSD uses the lattice free energy as heavy quark potential. But we can switch the potential to the internal energy or a combination of both energies. By solving the Schrödinger equation with the heavy quark potential, one obtains the dissociation temperatures and the binding energies of charmonia (, , ) and of bottomonia (), ). If a dissociation temperature is lower than , it is taken to be .
-
-
Vacuum limit of potential/spectroscopy: The heavy quark mass in vacuum is adjusted such that the ground state of quarkonium has the physical mass in the Particle Data Book. As a result, the charm quark mass is taken to be 1.26 GeV and the bottom quark mass 4.62 GeV. These are only used to solve the Schrödinger equation and the masses of open charm and bottom quarks in the PHSD are, respectively, 1.5 GeV and 4.8 GeV.
-
-
Reaction rates: The reaction rate of quarkonium is directly related to the reaction rate of heavy quarks, based on the Remler formalism [168, 169, 170, 171, 172]. It is similar to a quasi-free approach. However, if the heavy quark and heavy anti-quark form a bound state, the reaction rate will be smaller than twice the reaction rate of heavy quarks due to interference terms of heavy quark scattering and heavy antiquark scattering. So an adhoc suppression factor is introduced to the reaction rate of heavy quark, which is about 0.1 for bottomonia.
-
-
Assumptions about the medium (degrees of freedom, etc.): A QGP in PHSD is composed of off-shell massive partons whose pole mass and spectral width depend on temperature and baryon chemical potential and are fitted to the equation-of-state (EoS) from the lattice QCD calculations [74]. The scattering cross sections of heavy quarks and light partons are calculated at leading order in the strong coupling constant, but the propagator of massive offshell partons implements a resummation of the parton propagator and removes all singularities, making it unnecessary to introduce an explicit Debye screening mass [178].
-
-
Temperature dependence of heavy quark masses: In principle heavy quark will also be off-shell in medium. But that is not taken into account in the present form of the PHSD, since its effects on the observables are expected to be small. So the heavy quark mass is constant and does not depend on temperature in the PHSD simulations [179]. The charm quark mass is taken to be 1.5 GeV and bottom quark mass 4.8 GeV [180].
-
-
Equilibrium limits in transport: We have recently tested the Remler formalism which the PHSD adopts in thermalized and thermalizing boxes and found that the results are in good agreement with the statistical model [172].
-
-
Constraints from lattice QCD: In the PHSD approach, the pole mass and spectral width of massive off-shell partons which interact with heavy quark are fitted to lattice EoS through the strong coupling that depends on temperature and baryon chemical potential [74]. And the heavy quark spatial diffusion coefficient from lattice calculations is also described in the DQPM with an adhoc multiplication factor of 2 [180, 181].
-
-
Range of applicability and how is this range established: In principle there is no limitation for the application of the Remler formalism, because it is calculating the quarkonium density in quantum statistics. In practice, the dynamics of heavy quarks are calculated resorting to semi-classical trajectories which are a good approximation to the exact quantum mechanical evolution when the local temperature is larger than the binding energies. At lower temperatures, quantum corrections should however be applied.
-
-
Quantum features: Once the local temperature is lower than the dissociation temperature of quarkonium in heavy-ion collisions, there is a possibility for and which are close in phase space to form a bound state, based on the projection of the pair on the quarkonium wave function. Whenever a heavy (anti)quark scatters in the QGP, the projection is carried out and updated.
-
-
Regeneration: In the Remler formalism each scattering of a (anti)heavy quark brings about quarkonium dissociation and regeneration simultaneously, based on the change of 2-body Wigner density. Since the formalism provides only the change of quarkonium density with time, it is not simple to separate dissociation and regeneration.
-
-
Coupling to open heavy-flavor sector: Quarkonium production and dissociation are closely related to the dynamics of heavy (anti)quarks, since the quarkonium exists in QGP as a density probability, which is calculated from the distribution of heavy (anti)quarks in coordinate and momentum space.
-
-
Hadronic-Phase Transport: The dissociation cross section of quarkonium into a pair by light pseudoscalar or vector mesons is obtained by assuming a constant transition amplitude, which is fitted to experimental data at SPS energy [182]. The transition amplitude squared is multiplied by 3, if there is a vector meson such as or in the final state due to spin degeneracy, and divided by 3, if the light meson has strangeness. The reverse reaction is realized by the detailed balance.
-
-
Initial quark/quarkonium distributions: The initial spatial distribution of heavy quark pairs is given by the distribution of primary nucleon+nucleon scatterings in the Glauber model, and the initial momentum distribution is generated by the PYTHIA event generator [173] and then tuned by rescaling rapidity and/or transverse momentum of heavy quark pairs such that the distribution becomes similar to that from the FONLL calculations [183]. In pp collisions the produced heavy quark pair is projected on the quarkonium wave function after separating the heavy quark and antiquark by a distance which is inversely proportional to heavy quark mass on average. In heavy-ion collisions the projection is carried out when the local temperature is lower than the dissociation temperature of each quarkonium state.
-
-
Are cold nuclear matter effects, nPDF effects, etc. taken into account?: In the PHSD approach (anti)shadowing effects on heavy quark production are included using EPS09 [176, 184]. This suppresses the production of heavy flavor at low and at mid-rapidity, depending on collision energy. It also affects the production of quarkonium.
-
-
Constraints from pA and dA collisions: There is no additional constraint from pA and dA collisions except introducing EPS09 for the (anti)shadowing effects.
-
-
Medium evolution model: The production and time evolution of a nuclear medium in heavy-ion collisions are described using PHSD [185, 186, 187, 188, 189] where initial partons are produced through the string fragmentation or string melting, and then evolve and interact in off-shell dynamics. Around the critical temperature the off-shell partons hadronize to off-shell hadrons which eventually freeze out as on-shell hadrons.
-
-
Feed down implementation: We first realize the feed down in pp collisions, because experimental data are available there. The branching ratio of each excited state to a lower energy state is given by the Particle Data Book, and its contribution to the ground state is controlled by adjusting the radius of the excited states, because the yield of each state depends on its radius in the Wigner projection [175]. The same method is applied to heavy-ion collisions. We also take into account B-hadron decay to charmonium which is not negligible at LHC energies.
-
-
Comparisons to experimental data: We have tested our model first in pp collisions and found good agreement with the experimental data on the rapidity distributions and -spectra of , and their excited states [175, 177]. Based on these successful results in pp collisions, the Remler formalism is applied to bottomonia production in Pb-Pb collisions at 5.02 TeV and it produces the results consistent with experimental data [177].
-
-
Phenomenological breadth: Our model is applicable to any quarkonium state, if its wave function is given. This method does not depend on collision energy or system size.
4.5 The Saclay model
The aim of the “Saclay model” is to highlight the importance of the finite energy gap between singlets and octets (which essentially corresponds to the binding energy). It is an elaboration of a model that was developed in [190] for illustration purposes. We studied there the conditions that a master equation should fulfill in order to lead to thermal equilibrium. We found that taking into account the finite energy gap between energy levels was crucial, and we derived a general set of equations (valid in the limit) that describe an evolution in which the free energy monotonically decreases.
As a simple model we consider equations that consist of a rate equation for the singlet combined with a Langevin equation for the octet component. The rate equation that governs the populations of singlets ( and octets () is given by
(6) |
where , with the relative coordinate, and is a correlator of fields in the Coulomb gauge, evaluated in the hard-thermal-loop (HTL) approximation. The octet can also decay into a singlet. However, we checked in [190] that the probability for that to happen during the lifetime of the fireball is negligible, so we ignore it, which is achieved in practice by taking in eq. (6). We note that in this model, the decay width of the bound state depends both on its binding energy and its wave function.
In [191] we made a phenomenological application of this model. We did that considering two different choices for the real part of the potential. In the first case, we considered a Yukawa potential with the running of the coupling constants that would correspond to a computation in which . In the second scenario, we considered the real part of the potential that was fitted to the lattice data of [108]. The results presented here extend some of the results reported in [191].
Model ingredients
-
-
In-medium potential: As just mentioned, two different scenarios are considered, which, for simplicity, we name the “perturbative” and the “lattice-inspired” scenarios [191]. In both cases, we solve the Schrödinger equation using the real part of the potential. In the perturbative case, we use a Yukawa potential in which is evaluated at the Bohr radius of () and the for the calculation of the Debye mass is evaluated at . In the lattice inspired case, we use the real part of the potential fitted to the lattice data of [108].
-
-
Vacuum limit of potential/spectroscopy: In the perturbative case, we would get a Coulomb potential at zero temperature. In the lattice inspired scenario, we recover a Cornell potential. Note however that with the parameters used this Cornell potential (composed of a Coulomb plus a linear confining part) reproduces the spectroscopy only approximately.
-
-
Reaction rates: We assume that the leading mechanism is inelastic scattering with medium particles. This is computed in the HTL one-gluon exchange approximation using the wave functions and binding energies obtained by solving the Schrödinger equation with the real part of the potential. The octet wave function is modeled in the large limit. Our reaction rate can be encoded in a frequency-dependent imaginary potential (which reduces to the commonly used static imaginary potential when the frequency dependence is ignored). It was studied and computed in [190].
-
-
Assumptions about the medium (degrees of freedom, etc.): We assume the accuracy of the one-gluon exchange approximation and HTLs. The only exception is the real part of the potential in the lattice-inspired scenario, which is taken from a fit to lattice data.
-
-
Temperature dependence of heavy quark masses: Our heavy quark masses do not depend on temperature.
-
-
Equilibrium limits in transport: Strictly speaking, the equations we use do not bring the system to equilibrium because we did not include regeneration. However, we checked in [190] that this was an accurate approximation for the typical time scales that we are interested in at the LHC.
-
-
Constraints from lattice QCD: In the lattice inspired scenario we use as input lattice data on the real part of the potential.
-
-
Range of applicability and how is this range established: One condition for the model to be valid is that the decay width be much smaller than the energy gap between singlets and octets. We start with a quantum non-Markovian master equation and only in this limit we arrive to a Markovian rate equation that we can solve. The validity of other approximations used (one-gluon exchange and validity of the HTL propagators) is difficult to establish.
-
-
Quantum features: Quantum features are included in the rate (wave function and binding energy, as well as the gap between energy levels).
-
-
Regeneration: We do not consider regeneration since we work in the dilute limit (valid for bottomonium) and we have checked that in this limit regeneration effects are not important for the relevant LHC time-scales.
-
-
Coupling to open heavy-flavor sector: There is no direct relation to open heavy-flavor. In our model the octet evolves following a Langevin equation that uses the heavy quark diffusion coefficient as input. However, since we do not include regeneration, this has little impact on the results that we present.
-
-
Hadronic-Phase Transport: We do not consider hadronic effects.
-
-
Initial quark/quarkonium distributions: In the present simple approach, the survival probability coincides with . The spatial distribution of quarkonium at the initial time is proportional to the density of collisions (unintegrated number of collisions in the Glauber model). We do not consider the momentum distribution.
-
-
Are cold nuclear matter effects, nPDF effects, etc. taken into account?: No
-
-
Constraints from pA and dA collisions: No
-
-
Medium evolution model: We assume a Bjorken expansion using the equation-of-state of a free gas. The initial energy density scales as the density of participants.
-
-
Feed down implementation: We do not include feed down.
-
-
Comparisons to experimental data: Our motivation was to underline the importance of the finite energy gap between singlets and octets. No attempt was made to quantitatively reproduce experimental data.
-
-
Phenomenological breadth: Our main focus is on bottomonium since the dilute approximation and the condition are more likely to be fulfilled in that case.
4.6 The Santiago comover interaction model
The comover interaction model (CIM) [192, 193, 194, 195, 77] was originally proposed to explain the suppression of quarkonium states by final-state interactions with a hadronic or partonic medium. In this framework, the suppression arises from the scattering of the nascent quarkonia with comoving particles, i.e. particles with similar rapidities whose density is directly connected to the particle multiplicity measured at that rapidity for the corresponding colliding system. The comover interaction is governed by the Boltzmann gain and loss differential equations in a transport theory for a quarkonium state. The main parameter of the model is the comover-quarkonium cross section that results from the convolution of the comover momentum distribution in the transverse plane and the momentum-dependent comover-quarkonium cross section, proportional to the quarkonium geometrical cross section [196]. This model includes the initial-state nuclear effects through nuclear shadowing, i.e., the nuclear modification of the parton distribution functions. This model takes into account not only the dissociation of quarkonium, but also the possibility of recombination of into secondary quarkonium states [197, 198]. We recall two basic features of the comover approach. First, larger quarkonia are more affected by dissociation, due to larger interaction cross sections. As a consequence, excited states are more suppressed than the ground states. Second, the suppression increases with comover densities, which are proportional to particle multiplicities. Thus, the suppression increases with centrality in nucleus-nucleus collisions and is stronger in the nucleus direction for proton-nucleus collisions [199].
Model ingredients
-
-
Reaction rates: The rates are based on gluon and pion dissociation; they depend on the comover momentum distribution in the transverse plane (assumed to be characterized by an effective temperature, ) and the momentum dependent comover-quarkonium cross section as . The average over the comover energies is carried out using
(7) where corresponds to the threshold energy to break the quarkonium bound state (and as such is sensitive to the binding energies of the different states), is the energy of the comover in the quarkonium rest frame (with =0 and 140 MeV for gluons and pions, respectively), and is the transverse density of comovers, proportional to their multiplicities. The geometric cross section, , depends on the bound-state size; the phase space parameters 0.5-2 and the effective temperature, 200-300 MeV control the energy dependence and are fitted to the CMS and ATLAS p-Pb data.
-
-
Assumptions about the medium (degrees of freedom, etc.): The medium is formed by the comoving particles, considered to be partons or hadrons, i.e., gluons or pions. The relevant degrees of freedom are mostly hadrons in proton-nucleus collisions, whereas the gluons become relevant in the hotter nucleus-nucleus environment. As for the momentum distribution of the comovers in the transverse plane, we take a Bose-Einstein distribution .
-
-
Temperature dependence of heavy quark masses: The heavy quark masses are independent of temperature. They do not appear explicitly in the model.
-
-
Equilibrium limits in transport: We use the Boltzmann equation to describe the dissociation of quarkonium,
(8) where is the the energy-averaged quarkonium-comover-interaction cross section of quarkonium dissociation due to interactions with the comoving medium characterized by transverse density at . By integrating this equation from to , one obtains the survival probability of a quarkonium interacting with comovers, where the argument of the logarithm comes from converted to ratios of densities where we assumed that the interaction stops at when the densities have diluted, reaching the value of the proton-proton density at the same energy and rapidity, .
-
-
Constraints from lattice QCD: No lattice constraints are used.
-
-
Range of applicability and how is this range established: Our approach can be used to describe quarkonium production in both proton-nucleus and nucleus-nucleus collisions at SPS, RHIC, and LHC energies without the need to invoke any other phenomena. It may also be at play in high-density proton-proton collisions.
-
-
Quantum features: N/A
-
-
Regeneration: We consider the regeneration contribution for charmonium states. The regeneration through the uncorrelated charm and anticharm quark is represented by the gain term in the Boltzmann equation
(9) No regeneration is considered for bottomonium states.
-
-
Coupling to open heavy-flavor sector: N/A
-
-
Hadronic-Phase Transport: The Boltzmann transport equations apply for both the partonic and hadronic phases.
-
-
Initial quark/quarkonium distributions: The initial quarkonium and heavy quark distributions are proportional to the number of binary collisions.
-
-
Are cold nuclear matter effects, nPDF effects, etc. taken into account?: We take into account nCTEQ15 or EPPS16 shadowing. Nuclear absorption is included at low (SPS) energies.
-
-
Constraints from pA and dA collisions: The parameter that characterizes how quickly the cross section approaches the geometrical cross section, is fixed through a fit to proton-nucleus data. Both ground- and excited-quarkonium states are taken into account.
-
-
Medium evolution model: We assume a dilution in time of the comover densities due to longitudinal motion, which leads to a dependence on proper time, since Bjorken expansion is included. Transverse expansion is neglected.
-
-
Feed down implementation: Realistic feed-down contributions are included.
-
-
Comparisons to experimental data and phenomenological breadth: The model has been used to reproduce the nuclear modification factor of ground and excited charmonium and bottomonium states in proton-nucleus and nucleus-nucleus collisions at SPS, RHIC, and LHC energies.
4.7 The statistical hadronization model
The Statistical Hadronization Model (SHM) [200, 201, 202] assumes full dissociation of all quarkonium states in the QGP and exclusive generation at the QCD (crossover) phase boundary, in a concurrent hadronization process of all quark flavors. While the production of the heavy quarks occurs in hard parton scattering processes at initial stages of the collision, their thermalization, an essential condition for the applicability of SHM, is realized in the expansion of the deconfined medium. In the SHM, the absence of chemical equilibrium for heavy quarks is accounted for by introducing a fugacity factor ( in case of charm quarks). The fugacity is not a free parameter but is obtained from the balance equation [202] that accounts for the distribution of all initially produced heavy quarks into hadrons at the phase boundary, with a thermal weight constrained by exact charm conservation. With this approach, the knowledge of the heavy quark production cross section along with the thermal parameters obtained from the analysis of the yields of hadrons composed of light quarks [12], is sufficient to determine the total (-integrated) yield of hadrons containing heavy quarks in ultra-relativistic nuclear collisions.
With the assumption of the kinetic freeze-out taking place also at the phase boundary and employing hydrodynamics to determine the collective expansion velocity (at the freeze-out hypersurface), the transverse momentum distribution can be calculated as well. A corona contribution is added (both for the total and -differential yields), based on measurements in pp collisions.
Model ingredients
-
-
In-medium potential: Full dissociation (screening) is assumed for all quarkonium states.
-
-
Vacuum limit of potential/spectroscopy: Not directly relevant, but quarkonium (vacuum) masses are essential inputs for the model.
-
-
Reaction rates: N/A
-
-
Assumptions about the medium (degrees of freedom, etc.): Not explicit, but full (partial) thermalization of charm (bottom) quarks in the quark-gluon phase is assumed.
-
-
Temperature dependence of heavy quark masses: N/A
-
-
Equilibrium limits in transport: N/A
-
-
Constraints from lattice QCD: Not considered in an explicit way. Full dissociation of all quarkonium states is assumed.
-
-
Range of applicability and how is this range established: Applicable for AA collisions. Given the assumption of (full) heavy-quark thermalization, the SHM is most justified at LHC energies, but it was applied for lower energies too, down to SPS energies [203].
-
-
Quantum features: Not explicit (quantum numbers are of course considered and very relevant for SHM).
-
-
Regeneration: Exclusive generation at the QCD crossover boundary.
-
-
Coupling to open heavy-flavor sector: Full open heavy-flavor chemistry is predicted and plays a major role in the model. In fact, the knowledge of the full spectrum of the heavy-quark resonances is crucial for the quarkonium production in SHM. While the PDG hadron spectrum is the default implementation in the model, a version with enhanced charm baryon resonances was also explored [202]. This leaves the charmonium predictions unchanged, under the assumption that the possible existence of such additional states (not yet measured) leads to a commensurately larger (by 1%) total charm production cross section.
-
-
Hadronic-Phase Transport: None.
-
-
Initial quark/quarkonium distributions: No explicit treatment of initial heavy-quark distributions.
-
-
Are cold nuclear matter effects, nPDF effects, etc. taken into account?: Yes, currently based on the ALICE measurements of D mesons in (central) Pb-Pb collisions [204]. The charm production cross section in Pb-Pb collisions is the value describing in SHM the D data in central Pb-Pb collisions. The equivalent shadowing value, 30%, is used for the bottom sector as well.
-
-
Constraints from pA and dA collisions: Such constraints apply in SHM only if interpreted exclusively as shadowing on the total initial heavy-quark production. This was employed earlier, but is since recently superseded by the direct knowledge via the measurement of D-meson production in (central) Pb-Pb collisions.
- -
-
-
Feed down implementation: Feed down from excited quarkonia is considered, based on statistical production, more recently also from bottom hadrons.
-
-
Comparisons to experimental data: Full comparisons to the centrality dependence of integral yields () and to the spectra and vs. in central and semi-central collisions.
-
-
Phenomenological breadth: Yes, with the remark that treatment is based on an ad-hoc thermalization fraction in the bottom sector [208].
4.8 The Texas A&M University (TAMU) model
The starting point of the TAMU transport approach is a set of kinetic rate equations which describe the time evolution of quarkonium yields, , according to [81, 34]
(10) |
This provides a versatile tool that has been uniformly applied to a wide variety of quarkonia, , including both ground- and excited-state charmonia [209, 210] and bottomonia [211, 18], as well as mesons [212] and the exotic [213], over a large range of collision energies, from SPS via RHIC to the LHC [11]. The key transport parameter is the inelastic reaction rate, . In the QGP it is calculated based on in-medium HQ masses and binding energies obtained from a thermodynamic -matrix approach [75, 214] constrained by lattice QCD. The rate is dominated by “quasifree” dissociation processes [10, 34, 18], which are computed using perturbative diagrams with an effective but universal coupling constant as the main parameter. In the hadronic phase, the reaction rate is obtained from effective interactions with a large set of hadronic states (currently restricted to charmonia). The long-time equilibrium limit of each state, , is manifest in the regeneration term and computed from relative chemical equilibrium between the HF states in the system at given temperature. Transverse-momentum spectra have been computed from the Boltzmann equation for the suppressed primordial yields supplemented with a blast-wave approximation for the regenerated yields [40]. The latter has recently been improved by accounting for transported (off-equilibrium) HQ spectra obtained from relativistic Langevin simulations [15, 38].
Model ingredients
- -
-
-
Vacuum limit of potential/spectroscopy: The vacuum potential, as the zero-temperature limit of a screened Cornell potential, reproduces quarkonium spectroscopy (currently without spin-induced interactions). This provides an important benchmark for a smooth transition to realistic binding energies at moderate temperatures.
-
-
Reaction rates: The inelastic reaction rates for quarkonia are evaluated from both leading-order gluo-dissociation, , and next-to-leading order (NLO) inelastic parton scattering, where . The NLO processes are evaluated in a quasi-free approximation, as a half-off-shell inelastic scattering off the heavy quark (plus antiquark) in the bound state, thereby accounting for the in-medium binding and conserving four-momentum. The rate in the QGP has been multiplied by an additional phenomenological -factor of 3 extracted from d-Au collision data at RHIC (see the pertinent item below). For bottomonia, interference effects between the scattering off the and quarks are accounted for using an interference factor, , that depends on the size, , of the bound state via the 3-momentum transfer, . This leads to an -dependent reduction of the width, also referred to as an imaginary part of the potential; in practice, it is mostly relevant for the [18] and therefore not implemented for charmonia. Reaction rates in hadronic matter are included for charmonia as described below under Hadronic-Phase Transport.
-
-
Assumptions about the medium (degrees of freedom, etc.): The thermal QGP medium is modeled with massive quasiparticles that are used to compute the inelastic dissociation reactions. The quasiparticle masses are taken as , to approximately describe the energy density of the QGP down to temperatures of approximately 190 MeV. A transition to a hadron resonance gas is performed using a mixed-phase construction at a temperature of MeV (with charmonium reaction rates estimated from effective hadronic models). For bottomonia, hadronic dissociation is currently neglected. The QGP EoS has been updated with a lQCD parameterization [215], matched to a hadron resonance gas below a transition temperature of =170 MeV. The impact of this modification on bottomonium transport has been found to be small [18].
-
-
Temperature dependence of heavy-quark masses: The HQ mass is composed of a bare mass and a temperature-dependent in-medium contribution determined from the infinite-distance limit of the potential. The bare mass is fitted to the quarkonium masses in vacuum where the potential saturates at a string-breaking distance of about 1 fm. The temperature corrections to the HQ mass are constrained together with the in-medium potential by results from lattice QCD for the heavy-quark internal energy [75].
-
-
Equilibrium limits in transport: The quarkonium equilibrium limits are calculated from the statistical model. In the QGP, the total abundance of heavy quarks (with their in-medium masses) is assumed to be conserved (given by their production in primordial collisions), using fugacity factors, , throughout the fireball evolution. The quarkonium equilibrium number then follows from the standard thermal density multiplied with a factor of . In addition, a thermal-relaxation time correction is accounted for in the equilibrium limits to simulate the presence of non-equilibrated HQ distributions [39], which imply a reduced phase space for quarkonium production and thus a smaller equilibrium limit [36, 38].
-
-
Constraints from lattice QCD: The quarkonium spectral functions in the thermodynamic -matrix approach, from which the combination of in-medium quarkonium binding energies and HQ masses are taken, have been constrained by Euclidean correlation functions computed in lQCD [75], and cross-checked using the widths as employed in the rate equation [34]. Furthermore, the in-medium charm-quark masses have been checked against charm-quark susceptibilities in the QGP [214], cf. also Ref. [216] for more recent work.
-
-
Range of applicability and how is this range established: Both (LO) gluo-dissociation and (NLO) quasifree dissociation have been computed, with the former only relevant (albeit numerically small) in a small temperature window close to for and (1S,2S,1P); this follows the expected applicability for temperature ranges, and , respectively, although with a rather large coefficient for the temperature, , due to nonperturbative effects. A conceptual drawback is the current use of tree-level amplitudes to compute the quasifree reaction rates. A more realistic and consistent implementation directly using the nonperturbative in-medium -matrix amplitudes is currently underway.
-
-
Quantum features: The TAMU transport model is generally a semi-classical transport approach. A quantum feature is implemented in the early evolution after primordial production via bound-state formation times, ; they are different for the various quarkonia, scaled by their binding energies as (0.2-2 fm), thus increasing for higher excited states. The build-up of wave packets for heavy quarkonia is assumed to reduce their dissociation rates by a factor for , with the rates growing from zero to the equilibrium value linear in time, with an additional Lorentz-time dilation at finite .
-
-
Regeneration: Regeneration is included for all quarkonia, manifestly enforcing the statistical-equilibrium value in the long-time limit. For the spectra of the regeneration component thermal blast-wave spectra are employed following the fireball’s flow profile at an average regeneration temperature which is smaller for more loosely bound states which emerge at lower temperatures but at higher flow velocities.
-
-
Coupling to open heavy-flavor sector: In the baseline implementation, the effect of the thermalization process of HQ spectra is implemented using a thermal-relaxation time approximation for the quarkonium equilibrium limits. Explicitly transported charm-quark spectra have recently been employed using the resonance recombination model (RRM) with space-momentum correlations (SMCs) [15], which extends the relevance of regeneration contributions by about a factor of two in ; this much improved the description of the at intermediate in Pb-Pb collisions at LHC energies. Furthermore, the Boltzmann equation for charmonia has recently been solved by implementing time-dependent (transported) charm-quark distributions from Langevin simulations into the regeneration term [38]; these calculations have demonstrated the sensitivity of the (-dependent) regeneration yield to the degree of charm-quark equilibration (corresponding to the thermal-relaxation time factor in the equilibrium limit, cf. the item Equilibrium limits in transport).
-
-
Hadronic-Phase Transport: Transport in the hadronic phase is currently only implemented for charmonia. Starting from effective SU(4)-symmetric Lagrangians, reaction rates are calculated for charmonium dissociation induced by - and -mesons [39]. In Ref. [210] these calculations have been extended to include a large set of hadron resonance states through suitable changes in the available phase space. The rates in hadronic matter are generally small (up to MeV), while those for the are larger (up to a few tens of MeV) and phenomenologically relevant.
-
-
Initial quark/quarkonium distributions: The initial momentum distributions of primordial quarkonia in pp collisions are from fits to experimental spectra, while the initial spatial distributions follow a collision profile taken from the Glauber model.
-
-
Are cold nuclear matter effects, nPDF effects, etc. taken into account? Various CNM effects are included in the TAMU model: Nuclear shadowing with both and dependence; Cronin effect via a Gaussian smearing to simulate nuclear broadening [40], and nuclear absorption through effective (high-energy) - cross sections extracted from pA data at SPS and RHIC energies [40].
-
-
Constraints from pA and dA collisions: Data from p/dA collisions at SPS and RHIC are primarily used to constrain the CNM effects. However, the medium assumed to be formed in d-Au reactions at RHIC has also been used to determine a nonperturbative correction to the QGP suppression rate of the in terms of a phenomenological factor of 2-3 multiplying its quasifree reaction rate [210].
-
-
Medium evolution model: The fireball model we used in the calculation is an isentropically and cylindrically expanding isotropic fireball [81, 34, 40, 10]. This expansion model reproduces the hadron spectra at thermal freezeout that are consistent with the empirically extracted light-hadron spectra (, , ) similar to hydro calculations. The entropy is estimated from the multiplicities of observed charged particles and assumed to be conserved during the adiabatic expansion.
- -
-
-
Comparisons to experimental data: TAMU calculations have been used to compare to charmonium data at the SPS (Pb-Pb, S-U, In-In at 17 GeV), RHIC (Au-Au at 39, 62 and 200 GeV and Cu-Cu at 200 GeV), and the LHC (Pb-Pb at 2.76 and 5.5 TeV), bottomonium data at RHIC (Au-Au at 200 GeV) and the LHC (Pb-Pb at 2.76 and 5.5 TeV), and data at the LHC (Pb-Pb at 5.5 TeV) (as well as in p/dA systems, see above).
-
-
Phenomenological breadth: The TAMU model has been broadly applied to charmonium, bottomonium, , and phenomenology in AA and p/dA collisions at SPS, RHIC, and LHC energies, using the same input quantities within the same formalism. Specifically: , (spin averaged), [81, 34, 209, 210] and [213] for charmonia, , , , and for bottomonia [18], as well as and for charm-bottom mesons [212].
4.9 The Tsinghua model
In the Tsinghua transport model, the evolution of quarkonium in hot and dense QCD matter is described by a relativistic Boltzmann transport equation [83, 218, 219, 220, 221, 222]. The phase-space distribution of quarkonium states, , is controlled by the relativistic Boltzmann transport equation (with Bjorken coordinates),
(11) |
where represents different states (e.g., for charmonia); and are the rapidities in coordinate and momentum space, respectively, and is the quarkonium transverse velocity with the transverse energy . The anomalous suppression and regeneration mechanisms are reflected in the loss term and the gain term . Hot medium effects such as color screening effect and gluo-dissociation are included in . The regeneration process, which is the inverse of the gluo-dissociation, is represented byn ,
(12) | |||||
where is the gluon distribution taken as a Bose distribution; is the energy of gluon (heavy quark) and is the momentum of gluon (heavy quark). The dissociation rate, , contains the in-medium binding energy which is reduced by the color screening effect, and also the gluo-dissociation cross section in the reaction . In the term, the regeneration rate for the inverse reaction of gluo-dissociation is connected with the dissociation rate via detailed balance. Quarkonium regeneration also depends on the densities of heavy quarks, . The HQ density in the expanding medium is controlled by the diffusion equation, since heavy quarks are strongly coupled with the QGP. With this collisional term, the transport equation can be solved analytically [222].
In nucleus-nucleus collisions, the quarkonium initial distribution is treated as a superposition of quarkonium distribution in pp collisions. Cold-nuclear-matter effects, such as nuclear absorption, Cronin effect, and shadowing effect, are included by modifying the initial distribution extracted from pp collisions. For the non-prompt from decays, the Langevin equation is employed to simulate the energy loss of bottom quarks in the medium [223]. The hadronization of bottom quarks into mesons is described with the instantaneous coalescence model [223].
Model ingredients
-
-
In-medium potential: Charmonia and bottomonia experience a color screening effect in the hot deconfined medium, which reduces the heavy-quark potential and the in-medium binding energies. The color screening effect increases with the temperature and distance. As quarkonium dissociation mainly happens in the early stage of the medium evolution with high temperatures, we employ an effective constant in-medium binding energy in the calculation of quarkonium dissociation [219, 220].
-
-
Vacuum limit of potential/spectroscopy: The vacuum limit of potential is Cornell potential, with and . The charm and bottom quark masses are taken as 1.5 GeV and 4.5 GeV, respectively.
-
-
Reaction rates: Gluo-dissociation, , is considered as the dominant dissociation process in the QGP for tightly bound quarkonium. The cross-section, , in vacuum can be derived through the operator product expansion method and was calculated firstly by Peskin and Bhanot [224, 225]. The gluon density is taken as the Bose distribution , where and are the temperature and the four-velocity of the medium given by hydrodynamic models. At different temperatures, the channging density of thermal gluons gives the temperature dependence in the quarkonium dissociation rates. For excited quarkonium states, their dissociation rates are obtained via the geometric scaling of their radii over the ground-state one.
-
-
Assumptions about the medium (degrees of freedom, etc.): The hot deconfined medium generated in high-energy nuclear collisions is treated as an ideal gas consisting of massless / quarks and gluons, and strange quarks with a mass of MeV. There is a first-order phase transition between QGP and the hadronic gas, where the critical temperature is MeV at the zero baryon chemical potential.
-
-
Temperature dependence of heavy quark masses: The heavy quark mass is independent of temperature.
-
-
Equilibrium limits in transport: We use the Boltzmann equation to describe the dissociation and regeneration of quarkonium in relativistic heavy-ion collisions. The dissociation and regeneration are related to each other via the detailed balance. Thus, the equilibrium limit is naturally satisfied.
-
-
Constraints from lattice QCD: The in-medium heavy-quark potential is close to the internal energy extracted from the lattice QCD calculations. The quarkonium in-medium binding energy can be calculated with the two-body Schrödinger equation where the in-medium potential is taken as the internal energy extracted from the free energy evaluated in Ref. [226]. An effectively constant binding energy is then extracted and used in the calculation of the quarkonium decay rate.
-
-
Range of applicability and how is this range established: As the formation time/decoherence of quarkonia was not considered, the Tsinghua model mostly applies in the quantum optical limit. The charm quark phase-space distribution is assumed to be a kinetically thermalized distribution as mentioned in “initial quark/quarkonium distribution” section. This approximation is based on the observation of the large of mesons in experiment. It is assumed to be a good approximation for charm with low transverse momentum and/or in central collisions. In addition, regeneration is not considered for bottomonium due to the scarcity of bottom quark in the hot medium.
-
-
Quantum features: There are no quantum features in the Tsinghua model.
-
-
Regeneration: Regeneration is considered for charmonium states, via uncorrelated charm and anti-charm quark as represented by the gain term in the Boltzmann equation. The regeneration process is related to the dissociation process via detailed balance.
-
-
Coupling to open heavy-flavor sector: For prompt charmonium, the charm-quark is assumed to be kinetically thermalized as used in charmonium regeneration. For bottom quarks, energy loss in the medium is simulated with the Langevin equation. The final momentum distribution of mesons is used to calculate the production of non-prompt .
-
-
Hadronic-Phase Transport: Charmonium experience additional suppression through scattering with and mesons in the hadron gas via , and , with inelastic cross-sections taken from Ref. [227]; and mesons are assumed to be thermalized in the hadronic phase.
-
-
Initial quark/quarkonium distributions: The quarkonium initial distribution in nucleus-nucleus collisions is treated as a superposition of quarkonium distributions in pp collisions. The quarkonium momentum distribution is obtained by fitting the experimental data of distribution in pp collisions, with a further modification from CNM effects. In the regeneration part, where the charm-quark momentum distribution is needed, it is taken as a kinetically thermalized distribution. In this case, the initial momentum distribution of charm quarks does not affect charmonium regeneration. The initial spatial distributions of charm quarks and quarkonium are proportional to the density of binary collisions, .
-
-
Are cold nuclear matter effects, nPDF effects, etc. taken into account?: Yes, CNM such as nuclear shadowing, the Cronin effect, and nuclear absorption are included. The survival probability is related to the nuclear absorption cross-section [227], . The absorption cross sections for the excited states are obtained from through geometric scaling, , where the mean-square-radius can be obtained by solving the two-body Schrödinger equation [222]. The Cronin effect is included in the initial distribution of quarkonium via the Gaussian smearing method [228]. The nuclear shadowing factor is calculated with the EPS09 package [176].
-
-
Constraints from pA and dA collisions: The nuclear absorption cross-section, , and Cronin momentum-broadening parameter, , are determined by fitting experimental data from pA and dA collisions.
-
-
Medium evolution model: 2+1D ideal hydrodynamics is mostly used to describe the medium evolution in Tsinghua model for many years [229]), while a 3+1D viscous hydro (MUSIC package) has been used in past two years. In the present context paper, the results are obtained with the 2+1D ideal hydro. The initial condition (entropy density) of the hydrodynamics in the transverse plane is given by the two-component model [230]. The maximum entropy density is determined by the charged-hadron multiplicity observed in the experiment. The start time of hydrodynamics is and the maximum temperature is for central Pb-Pb collisions with .
-
-
Feed down implementation: For , the feeddown contributions from and are considered. For , the feeddown contributions from higher states such as , , , and are considered. The feeddown branching ratios are taken fromthe PDG [231].
-
-
Comparisons to experimental data: The Tsinghua model has been used to explain the nuclear modification factor, , anisotropic flow, , and mean transverse-momentum squared, , of charmonium and bottomonium states.
-
-
Phenomenological breadth: The Tsinghua model has been applied in the studies of charmonium, bottomonium, and in both small (p-Pb) and large (Pb-Pb, Au-Au) collision systems, at the collision energies of SPS, RHIC, and LHC.
4.10 Cold nuclear matter effects
Cold nuclear matter effects generally refer to all modifications present in pA collisions, when there is a nuclear target, but not in pp collisions. These effects are referred to as being “cold” in the sense that it is assumed that no quark-gluon plasma is created. Note that this assumption may be challenged in high-multiplicity pp and pA collisions. These effects are also presented in AA collisions when quark-gluon plasma is also created. The nuclear modifications of the parton distributions (nPDF effects) are typically included by a parameterization such as the EPPS16 [138] parameterization. The centrality dependence of nPDF effects may also be taken into account, see Ref. [232]. Enhanced broadening in the nucleus relative to a proton, due to multiple scattering or a Cronin-type effect, is also sometimes included [233]. The broadening may or may not be related to energy loss in cold matter. A number of energy loss models have been proposed, see for example Ref. [234]. Such models can account for quarkonium suppression at high . Other models that include intrinsic charm without energy loss can also provide a good description of this suppression, see Ref. [233].
Quarkonium absorption by nucleons has been studied by many different groups but has been suggested to be negligible at the LHC. See Ref. [235] for a discussion of the energy dependence of absorption. Dissociation of quarkonium by comovers has also long been suggested, see Sec. 4.6. The comovers can be considered to be partons or hadrons, see Refs. [236, 237, 238] for some early discussions of hadronic comover dissociation. It was found that, in this description, the nuclear dependence of comover dissociation is similar to nuclear absorption [239].
In the traditional color evaporation model (CEM), the quarkonium production cross section is some fraction, , of all pairs below the threshold where is the lowest mass heavy-flavor hadron. The color of the octet state is ‘evaporated’ through an unspecified process which does not significantly change the momentum. The quarkonium yield may be only a small fraction of the total cross section below . Schematically, the production cross section of quarkonium state in a pp collision is
(13) |
where or and is the cross section and and are the factorization and renormalization scales, respectively. In pA collisions, the cold nuclear matter effects on quarkonium production are
(14) |
where
(15) | |||||
(16) |
where is the survival probability for nucleon absorption, is the nuclear modification of the parton distributions (nPDF), and and account for transverse momentum broadening in the proton and nucleus, respectively [233].
The CEM has been improved (ICEM) to better account for feed down and the quarkonium mass [240]. The unpolarized direct quarkonium production cross section in pp collisions in the ICEM is
(17) |
Note that the change in integration range changes to while the change in the momentum range modifies the distribution of the quarkonium states relative to each other.
Model ingredients
-
-
In-medium potential: N/A
-
-
Vacuum limit of potential/spectroscopy: N/A
-
-
Reaction rates: N/A
-
-
Assumptions about the medium (degrees of freedom, etc.): N/A
-
-
Temperature dependence of heavy quark masses: N/A
-
-
Equilibrium limits in transport: N/A
-
-
Constraints from lattice QCD: N/A at the moment. If lattice could provide constraints on feed down, then it would be useful. Also, for NRQCD-type formulations, calculating the LDMEs could be useful.
-
-
Range of applicability and how is this range established: The calculation of initial production is applicable over all center of mass energies.
-
-
Quantum features: N/A
-
-
Regeneration: N/A at the moment, however, in the past the HVQMNR code was used by Thews and Mangano to calculate regeneration by “non-diagonal quarkonium production” [241].
-
-
Coupling to open heavy-flavor sector: N/A
-
-
Hadronic-Phase Transport: N/A
-
-
Initial quark/quarkonium distributions: The charm quark mass and scale parameters used to calculate the distributions are [242] . In the case of production, . The values of are fixed for the central parameter set in each case and all calculations employing other masses and scales use the same value of to obtain the extent of the and mass and scale uncertainty bands. The normalization factors for the CEM are for the central result with and with . The calculations use the CT10 proton parton distributions [243]. The quark mass makes a larger contribution to the cross section uncertainty than does the scale choice [243].
-
-
Are cold nuclear matter effects, nPDF effects, etc. taken into account?: The EPPS16 [138] parameterization is employed for the nPDF calculation. The intrinsic broadening employed in production is augmented in nuclei according to multiple scattering in the nucleus [233]. While absorption [235] is included at lower energies, it is considered negligible at LHC energies.
-
-
Constraints from pA and dA collisions: There are very few constraints from pA and dA collisions in the basic model. There are constraints from these collisions in the global analyses of the nPDFs used in the calculations. For example, EPS09 [176] used data from RHIC to constrain the gluon distribution and EPPS16 [138] used LHC p-Pb data on and production to separate the antiquark distributions in the sea at high and on dijet data to further constrain the gluon distribution. The centrality dependence of was used to study the centrality dependence of shadowing [232]. Fixed-target pA data were used to constrain the absorption cross section [235]. However, absorption is considered negligible at the LHC collider energies and is not included in the calculation [233]. The intrinsic employed in pp collisions was obtained from comparison to data but the broadening in pA collisions is based on a model and not a fit within these calculations [233].
-
-
Medium evolution model: N/A
-
-
Feed down implementation: The traditional CEM does not distinguish between the states, all of the distributions are assumed to be the same. Thus, the feed down distributions are all identical to that of the ground state modulo emission of soft particles, photons or pions, in the decays. In the improved CEM, the feed down distributions depend on the specific mass and momentum of the particular quarkonium state and are thus realistic in terms of the model [240, 244, 245, 246]. For example, the model reproduces the dependence of the ratio [240].
-
-
Comparisons to experimental data: The model agrees well with experimental data, both for the traditional CEM and the improved CEM [246].
-
-
Phenomenological breadth: Both the traditional CEM and the improved CEM can be applied to charmonium and bottomonium production, including both S and P states, , (2S), and (S) as well as and states respectively [243, 244, 245, 246]. Charmonium production in particular has been studied from GeV to TeV [233].
The improved CEM has been developed to also calculate quarkonium polarization [247, 248, 244, 245, 246]. The calculation was also extended to Pb-Pb collisions including only cold nuclear matter effects [249] and was shown to agree with the data. Note that because the polarization calculation involves a ratio of cross sections with a certain spin, such effects generally cancel in the ratios [249]. Indeed, only the charm quark mass generally affects the polarization but is not a strong effect [244].
The ability of the CEM to consistently cover the entire energy range, from near production threshold to the highest available energies can be contrasted with the NRQCD approach which typically requires a cut to fit the LDMEs [250] to the data and has only been matched to the low part of the spectrum by including a color-glass condensate contribution, matching to the high part [251]. Such an approach cannot work at lower energies, such as fixed-target energies. Additionally, the NRQCD LDMEs that have been fit to total cross section quarkonium data [252] also cannot describe distributions [253].
5 Comparisons of Model Ingredients and Interpretation of Results
In this section we will confront the theoretical inputs and results of the transport models to quantitatively analyze key components of in-medium quarkonium kinetics and how they manifest themselves in phenomenologically-relevant outcomes. We start by collecting the temperature evolution of the medium expansion models in Sec. 5.1 for the case of central Pb-Pb ( = 5 TeV) collisions and then turn to the arguably most fundamental transport parameter, i.e., the reaction rate, for various charmonia and bottomonia as a function of temperature in Sec. 5.2. To aid in the interpretation of these results, we inspect the inputs for binding energies and HQ masses in Sec. 5.3, the resulting 3-momentum dependence of the rates in Sec. 5.4 and the spatial dependence of the imaginary part of the potentials (as applicable) in Sec. 5.5. We then turn to more phenomenologically-oriented studies, by testing the medium evolution models with a prescribed temperature dependent reaction rate to compute an in Sec. 5.6, and by studying the impact of formation time effects in early quarkonium evolution on their suppression factor in Sec. 5.7.
5.1 Medium Evolution
In Fig. 7 we collect calculations of the temperature evolution for the central cell in 0-10% central Pb-Pb collision at TeV for the expanding QCD medium using the approaches included in the present effort. For the hydrodynamical models (Duke, Munich-KSU, Nantes, and Tsinghua), the results shown are for the central cell (), which is usually the hottest region of the plasma; the 4 results agree quite well, in particular for the former three, while the latter one cools faster due to the EoS based on a massless parton gas (and a first-order transition) which is “harder” (generating a faster transverse expansion) and has a smaller temperature at given entropy density, compared to the lQCD-based EoS used in the other 3 models. A marked difference in the hydro models is the assumption of the initial thermalization time, , which controls the initial temperature, , varying between close to 700 MeV at = 0.2 fm for Munich-KSU and 400 MeV at 1.5 fm for Nantes; this may also be affected by the initial temperature profiles used in the calculations (e.g., two- vs. three-dimensional). For the Saclay (ideal Bjorken hydrodynamics), Santiago (comover), Tsinghua (2+1D ideal hydrodynamics), and the TAMU fireball (isentropically expanding blastwave with lQCD EoS) models, spatial averages lead to lower temperatures, and the primarily longitudinal expansion in the former two leads to a significantly slower cooling at later times (which affects the more weakly bound quarkonia). In the SHM, the MUSIC 3+1D hydrodynamic calculations [256] with IP-Glasma initial conditions are used in a blast-wave parametrization, although the pertinent quarkonium production is only evaluated at the pseudo-critical temperature, 155 MeV. For the Nantes approach additional information is provided in terms of the average temperature encountered at the position of a HQ pair propagating through the fireball, which is significantly smaller than the values in the central cell but agrees fairly well with the spatially uniform profile of the TAMU fireball evolution. The temperature evolution extracted from the PHSD transport model has a rather steep decrease, starting at about 660 MeV and dropping down to 150 MeV after about 9 fm.
5.2 Reaction Rates
5.2.1 Charmonia
The available results for charmonium reaction rates, i.e., the 1S and 2S vector states, are compiled in Fig. 8 as a function of temperature for 2 different momenta. We note that most OQS approaches work with an EFT scale hierarchy which is not applicable to charmonia. In the Nantes approach, results are only shown for temperatures below the dissociation temperature, while for TAMU the rates above the dissociation temperature simply correspond to two times the HQ collision rate (as a means to characterize the decorrelation of a primordially produced HQ pair that would end up in a quarkonium state in a pp collision). The Tsinghua model assumes a constant binding energy and geometric scaling for the rates of excited states, while the comover model does not evaluate charmonium rates microscopically but via a constant dissociation cross section fit to data for each state.
At vanishing momentum, the rates from TAMU, Tsinghua and Nantes agree well at temperatures below 350 MeV. This bodes well for the comparison between TAMU and Nantes, as both rates are computed with a similar mechanism (inelastic parton scattering closely related to the imaginary part of the potential) and a potential that leads to similar in-medium binding energies, cf. Fig. 10. For Tsinghua and Nantes the further increase beyond MeV is markedly stronger, but likely due to different mechanisms. For Tsinghua, this is rooted in the gluo-dissociation cross section with a constant in-medium binding energy which is peaked near the binding energy; at higher momenta this reduces the rate at higher . For TAMU, once the binding energy vanishes for MeV, and for Santiago with constant comover cross section, the dependence becomes rather weak, close to linear or weaker. This can also be seen for the . For the latter, the Nantes and TAMU rates are again comparable, albeit in this case with rather different input for the bound-state properties. The Santiago and Tsinghua rates are obtained from the rate by geometric scaling of the radii. Especially the Santiago rates are much larger than all other rates at low (we recall that they can explain the suppression in pA collisions), while the TAMU rates include an extra -factor of 3 to mimic nonperturbative effects and enable a better description of pA results.
At finite charmonium momentum (taken as =5 GeV in the bottom panels of Fig. 8), the Nantes and comover results remain unchanged for both and . On the other hand, the TAMU results increase (as consequence of a perturbative matrix element for the heavy-light coupling and a larger phase space available for dissociation) while the Tsinghua results decrease (as a consequence of the gluo-dissociation cross section being probed beyond its peak structure).
5.2.2 Bottomonia
Next, we turn to the dependence of the inelastic bottomonium rates, compiled in Fig. 9, again for 2 different 3-momenta. The Duke, Munich-KSU, Nantes, TAMU, Tsinghua, PHSD and Santiago (comover) rates for at =0 approximately agree in the phenomenologically most relevant range of 300-400 MeV. While the underlying mechanism (inelastic parton scattering) and in-medium binding for Nantes and TAMU are similar (see Fig. 11), the Tsinghua (gluo-dissociation) and Duke results utilize constant quark masses and binding energies, which are also quite different from each other being based on either a Cornell potential (large binding and large HQ mass) or color-Coulomb potential (small binding and small HQ mass), respectively. The widths for Saclay (with a binding similar to Nantes and TAMU) are rather small. Finally, the PHSD and comover results have the weakest -dependence, and they are quite large at relatively low . Moving on to =10 GeV, the Santiago, Munich-KSU, Nantes, PHSD and Saclay widths do not change significantly (or at all), while the Duke and TAMU widths increase and the Tsinghua width decreases. The overall spread in the results increases. For the , the calculated widths generally increase substantially relative to the , especially at low temperatures, with the exception of PHSD (which utilize the same reaction rates for and ). The large spread (the agreement between Nantes and Tsinghua for =0 must be considered a coincidence given the very different mechanisms), suggests that phenomenological constraints currently suffer from large uncertainty; e.g., even small contributions from regeneration can make a big difference in a phenomenological extraction based on a small , whereas regeneration is expected to be relatively less important for the (except for the SHM).
5.3 Quarkonium binding energies and heavy-quark masses
As already mentioned above, the in-medium binding energies of quarkonia are an important ingredient to compute their dissociation rates, as this determines the available phase space for inelastic reactions, and more compact states (at larger binding) are also subject to interference effects (that depend on the wave function) which further suppress the rate. We define the binding energy in the standard way as the difference between the bound-state mass and the threshold (defined as twice the in-medium HQ mass).
For charmonia, shown in Fig. 10, we can reiterate what we discussed before, i.e., a rather good agreement between Nantes and TAMU, while Tsinghua uses a constant , in part motivated by the restriction to gluo-dissociation which becomes increasingly inefficient for small [10] but is extended to excited states by geometric-size scaling. For the excited states the Nantes binding energies are somewhat larger, even exceeding the commonly quoted vacuum values (relative to the open-charm threshold), see, however, Sec. 4.3 for caveats on the model applicability at low .
For the bottomonium sector, more results are available. One may broadly classify them as constant vs. -dependent and being based on either a Cornell or a color-Coulomb only potential. Specifically, the Duke, Munich-KSU and Tsinghua groups employ constant binding energies that are at about half of the vacuum value (similarly for and for Duke and Munich-KSU), while Nantes, TAMU, Saclay and PHSD have a strong -dependence which, for the former three, starts out near the vacuum value around . Furthermore, only Nantes and TAMU have an in-medium -quark mass, which drops with and results in a near-constant mass close to the vacuum value. Most other approaches have a significantly smaller mass even close to . Similar observations also hold for the excited states, where in particular the color-Coulomb based approaches have much smaller bound-state masses compared to the vacuum values, mostly due to the relatively small -quark masses, cf. Fig. 12.
5.4 Momentum Dependence of Reaction Rates
Transverse-momentum spectra are a pivotal observable in heavy-ion collisions. In the context of quarkonium transport models one needs to calculate the 3-momentum dependence of the dissociation rate, which is directly reflected in the nuclear modification factor of the primordial-production component. However, because of detailed balance, the 3-momentum dependence also affects the regeneration of quarkonia, which is additionally influenced by the interplay with the phase space distribution of heavy quarks as they diffuse through the expanding fireball.
We start again by inspecting available model results for charmonia, for 3 different temperatures, see Fig. 13. The Santiago (comover) and Nantes models are currently restricted to =0 and therefore have no -dependence. The Tsinghua results for the exhibit the typical momentum dependence generated by the peaked structure of the gluo-dissociation cross section, namely initially large values that quickly fall off once the center-of-mass energy exceeds the peak position. This effect is further accelerated with increasing temperature due to the increase in thermal-parton motion. As mentioned before, for the excited states, the rates are geometrically scaled by the bound-state size. In the TAMU model, the rate at =200 MeV exhibits a strong increase with momentum which is driven by steep increase in phase space due to a relatively large binding energy of about 550 MeV, requiring a rather large threshold momentum of the thermal partons. For the , where 10 MeV at = 200 MeV (and vanishing at higher ), the phase-space suppression is absent and the remaining increase with of about 25% is due to the perturbative-scattering matrix element, as can also be seen for the at = 400 MeV.
The 3-momentum dependence of bottomonium reaction rates is compiled in Figs. 14, 15 and 16 for temperatures of 200, 300 and 400 MeV, respectively. At = 200 MeV, the spread in the rates is substantial, spanning more than an order of magnitude for all states. For the , this spread is mitigated by the overall smallness of the rates (except for PHSD and the Santiago models, recall Fig. 9) which renders them phenomenologically of little relevance. Nevertheless, one finds a qualitative ordering according to the binding energy (recall Fig. 11), with decreasing rates following the increasing binding energy from PHSD to Duke to Munich-KSU to Tsinghua to Nantes and TAMU. While the Santiago, Munich-KSU and Nantes models have currently no -dependence, PHSD and Duke show an increasing trend which is even more pronounced in TAMU due to the large binding (suppressing the phase space at small ), while the gluo-dissociation of Tsinghua is decreasing (roughly comparable to the TAMU gluo-dissociation for the which has approximately the same binding at this temperature). For the the spread is also large (with little -dependence in most models), which, as noted before, is likely due to current uncertainties in the phenomenological implementations. Similar remarks apply to and .
At higher temperatures, the model agreement for the is better (as noted in the context of the -dependence). The TAMU and Tsinghua rates for the excited states at MeV are probably underestimated, since once the binding vanishes, one is basically dealing with heavy-light couplings, which are underestimated when using perturbative diagrams (recall, e.g., the -factor for the ) [23, 26]. This is probably more realistic in models based on the dipole expansion where the rate is taken as (with the HQ momentum diffusion coupling, , usually taken from lattice-QCD results), which incorporates non-perturbative interaction strength (which is large [13]). However, this expansion breaks down for small where becomes large, and similarly for the quantum-optical limit in the Duke model (where, e.g., GeV at = 400 MeV).
5.5 Imaginary Potential in Coordinate Space
Next we turn to models that include an explicit (-dependent) imaginary part in their heavy-quark potentials. This enters into transport models based on open quantum systems where the quarkonium wave function is coupled to a QGP medium consisting of light degrees of freedom 111An -dependent dissociation rate can also be obtained by incorporating an interference factor to diagrammatic computations with subsequent use in semiclassical transport [257, 18]. Upon integrating out the light degrees of freedom, one can obtain an evolution for the heavy-quarkonium reduced density matrix. In the Markovian limit that emerges when the time scale for medium relaxation, , is much faster than the time for internal transitions (parametrically given by ), the resulting evolution equation is of Lindblad form [166, 19]. The real part of the potential can be combined with the transition rates, defined in terms of jump operators, into a complex effective Hamiltonian. The Lindblad equation is then solved by using either the method of quantum trajectories or quantum diffusion with a stochastic potential where the correlators for the stochastic terms in the potential are set by the imaginary part of the potential. This leads to a non-unitary dynamics of the quantum mechanical wave function via the complex effective Hamiltonian, which results in suppression of quarkonium states.
In Fig. 17 we compare the imaginary part, , of the singlet effective potential, , used in the Nantes and Munich-KSU OQS approaches. At low temperatures (upper panels), the Nantes and Munich-KSU results are rather different, while at higher temperatures (lower panels) they are closer, although the functional form remains qualitatively different. In the Munich-KSU approach, the imaginary part of the singlet effective potential is manifestly independent of the HQ mass and, at leading order in a expansion, given by , where is the scaled heavy-quarkonium transport coefficient which can be defined via chromoelectric correlators; here, a temperature-dependent parameterization is used that was fit to lQCD data in Ref. [116]. We note that beyond leading order in the imaginary part of the effective Hamiltonian cannot be expressed as a local potential as there appear additional contributions proportional to the anti-commutator between the HQ relative momentum and their distance. For purposes of comparison, both the Munich-KSU and Nantes groups did not include such momentum-dependent contributions to the width.
The dipole expansion in the Munich-KSU approach, which causes the behavior, limits the reliability of the imaginary part of the potential to relatively small distances. In practice, this problem is mitigated since the imaginary part of the effective potential causes the wave function to be rather localized near the origin when using the effective Hamiltonian evolution. From Fig. 17 one sees that the imaginary part in the Nantes approach has a weak dependence on the HQ input mass. With increasing temperature its structure is compressed to smaller distances, a feature that is also seen in the Munich-KSU framework.
5.6 Nuclear Modification Factor with a Common Reaction Rate
Another comparison of the models was conducted for the observable for direct production by imposing a simple but uniform parametrization for (assumed to be constant vs. momentum), consisting of a linear increase from 0 at =200 MeV to 0.2 GeV at =600 MeV, which every group implemented in their respective evolution model222For models not directly relying on as an input for the dynamical evolution, a proper rescaling of the corresponding quantity – like the imaginary potential – has been applied.. Initial formation time effects for the bound states, as well as CNM effects and feeddown, have been neglected. The results are shown in Fig. 18 as a function of centrality and . For the TAMU model the regeneration component is shown separately as well as its sum with the (suppressed) initial production. The results of the SHM, included for the case of 50% thermalized bottom quarks, constitute pure generation at .
At first sight, models vary significantly in the way they describe centrality dependence. However, closer inspection reveals that the Nantes results are anomalously high as the underlying EPOS4 background is assumed to thermalize only after which leads to a rather low initial temperature, recall Fig. 7. The TAMU result is reasonably close to the 3 curves by Saclay, Munich-KSU and Tsinghua, with a shape close to Munich-KSU but with a higher yield in peripheral collisions which is presumably caused by larger medium thermalization times leading to smaller initial temperatures, cf. Fig. 7, while Tsinghua has a smaller initial temperature in central collisions but a larger one at ”intermediate” centralities which could explain the stronger suppression in more peripheral collisions. The Saclay and the Munich-KSU results are closest to each other, which may be due to due to longer lifetime in the Bjorken model (Saclay) and escape effects in the Munich-KSU which render the latter’s higher in peripheral but lower in central collisions where the initial temperature is higher. Both models generate a stronger suppression than TAMU, which is largely consistent with the temperature evolution shown in Fig. 7.
5.7 Quarkonium Formation Time Effects
To study the impact of quarkonium formation times calculations were performed starting from a “realistic” initial state (the one used in the respective dynamical model, usually reported as a ”point-like initial state” in the OQS and the ground state in semiclassical approaches). This state was evolved in a QGP at fixed temperature =300 MeV, neglecting regeneration. The models provided the ”survival” probability as a function of time to find this pair at =0 in an eigenstate of the in-medium potential.
Figure 19 illustrates how suppression mechanisms underlying the calculations of the decay rate are realized in a basic time evolution scenario at constant . We focus the discussion on the bottomonium case which was addressed by most of the groups. In the TAMU approach, where the initial state is an in-medium state and regeneration mechanisms were discarded for the purpose of this study, one finds as expected a survival probability , where agrees with the reaction rate displayed in Fig. 9 (including the gluo-dissociation mechanism); the inclusion of a formation time typically delays the evolution, with an offset of . The same exact agreement with the exponential decay law is obtained in the Saclay calculation as the regeneration was not considered in this implementation of the model. In the Duke approach the regeneration component was not removed, leading to a slight deviation with respect to the exponential decay initiated with a vacuum state of the 333Note that ”in-medium states” are not defined in the Duke approach owing to the hierarchy., of the order of 5% after 8 fm/. In the Munich-KSU calculation, the evolution starts from a compact state close to a Dirac -function peak. While the evolution of the survival probability of the in-medium state decreases nearly exponentially, the associated decay rate is found to be twice the imaginary part of the eigenvalue ( GeV) corresponding to the fundamental eigenstate of the non-hermitian effective Hamiltonian, which is an indication that this mode governs the evolution starting from very early time. In the Nantes approach, the initial condition is taken as a compact state as well. At early time (during typically the first 1 fm/), the dynamics is impacted by quantum interference which leads to a non-exponential decay of the survival probability and could be interpreted as an effective formation time. After 1 fm/, one observes an exponential decrease, with a decay rate close to the one evaluated in Fig. 9 for an in-medium state.
6 Conclusions
In this report we have summarized the common effort of a task force, composed of various theoretical groups, to scrutinize their models for the description of quarkonium production in high-energy heavy-ion collisions. To begin with, a synopsis of each model was given guided by 19 items to specify the inputs and spell out underlying assumptions. This revealed a large variety of the approaches, e.g., input potentials ranging from perturbative color-Coulomb to Cornell potentials with different degrees of screening, reaction rates based on gluo-dissociation and/or inelastic scattering, implemented into semiclassical vs. quantum transport models, deviations in the treatment of regeneration (from none to diagonal to multiple independent recombinations), and different accounting for cold-nuclear-matter effects, to name a few. A set of calculations of derived quantities has been defined to study how these differences manifest themselves in key quantities of the transport approach, specifically in-medium binding energies and HQ masses, inelastic reaction rates in momentum space (and/or imaginary parts of the potential in coordinate space), and nuclear modification factors with controlled rates and medium conditions.
Let us summarize the results in light of the 5 basic questions posed in the introduction. First, concerning the model consistency, the semiclassical approaches employ rather different inputs for the in-medium binding energy and the reaction rate (which is strongly affected by the binding energy), while all of them include regeneration in a way that accounts for multiple heavy-quark pairs in the system; in current quantum-transport approaches, which mostly focus on bottomonia, regeneration, if any, is only based on a single pair, sometimes further restricted to the configurations that in a pp collision would form an state (which is a small subset). Second, concerning the equilibrium limits relevant for regeneration, they have not been explicitly compared, but from the model descriptions it appears that currently only the semiclassical treatments have control over them, although they are expected to be quantitatively different due to different HQ masses and binding energies. Third, concerning the significance of quantum effects, the comparisons of suppression factors appear to confirm that these are mostly relevant in the early stages; the long-time behavior of suppression can be characterized by exponential decays that correspond to the pertinent reaction rate in semiclassical approaches and the lowest eigenvalue in the quantum approaches. Fourth, concerning lQCD constraints, they have been implemented in a number of models, either in terms of directly computed quantities, i.e., transport coefficients (which, however, are restricted to vanishing 3-momentum), or more indirectly by computing lQCD quantities (e.g., free energies or Euclidean correlators) within a model approach to constrain its input quantities (like the potential or HQ masses); the latter variant usually offers broader phenomenological flexibility as well as microscopic insights. Finally, concerning the ultimate model uncertainties, it will be necessary to go beyond a diagnostic level of comparing the model calculations of specific quantities as conducted in this work. Clearly, a more systematic implementation of lQCD constraints on the input quantities (such as the in-medium potential) on an equal footing across model approaches is desirable. Then one could envision that a comparison of semiclassical to quantum transport approaches with the same microscopic input (such as in-medium potential and HQ masses, evaluated within the same process for the reaction rate, e.g., inelastic parton scattering) could reveal systematic uncertainties in the transport part. This might also identify in how far simplifications in certain components of the model calculations are justifiable. In doing all this, it will be important to account for the strongly coupled nature of the QGP, which manifests itself in both binding energies and reaction rates. While the former are a direct consequence of the in-medium QCD force (input potential), the latter are closely related to the single HQ transport coefficient, which has been experimentally established to be in the strong-coupling regime. Since the HQ interactions with the QGP which govern HQ transport can also be expected to be operative in the bound-state properties of quarkonia, a self-consistent treatment will provide the strongest constraints.
Acknowledgements
We thank the Extreme Matter Institute for support for this Rapid Reaction Task Force. The work of XD was supported by the Deutsche Forschungsgemeinschaft (DFG, German Research Foundation) through the CRC-TR 211 ’Strong-interaction matter under extreme conditions’– project number 315477589 – TRR 211; XD, MAE and EGF acknowledge support from the Xunta de Galicia (Centro singular de investigacion de Galicia accreditation 2019-2022), the European Union ERDF, the “Maria de Maeztu” Units of Excellence program under projects CEX2020-001035-M and CEX2019-000918-M, the Spanish Research State Agency under projects PID2020-119632GB-I00 and PID2019-105614GB-C21, and the European Research Council under project ERC-2018-ADG-835105 YoctoLHC; the work of MAE has also been supported by the Generalitat de Catalunya under grant 2021-SGR-249. MS was supported in part by U.S. Department of Energy award No. DE-SC0013470; PP was supported by the U.S. Department of Energy, Office of Science, Office of Nuclear Physics through Contract No. DE-SC0012704; XY was supported in part by the U.S. Department of Energy, Office of Science, Office of Nuclear Physics, InQubator for Quantum Simulation (IQuS) under Award Number DE-SC0020970 via the program on Quantum Horizons: QIS Research and Innovation for Nuclear Science; SD is supported by the Centre national de la recherche scientifique (CNRS) and Région Pays de la Loire and acknowledges the support of Narodowe Centrum Nauki under grant no. 2019/34/E/ST2/00186. PBG was supported by the European Union’s Horizon 2020 research and innovation program under grant agreement No 824093 (STRONG-2020). PVG was supported by the U.S. Department of Energy Award No. DE-SC0019095 and is grateful for the support and hospitality of the Fermilab theory group; RR and BW have been supported by the U.S. National Science Foundation under grant nos. PHY-1913286 and PHY-2209335; PP, RR, MS, RV and BW have been supported by the U.S. Department of Energy, Office of Science, Office of Nuclear Physics through the Topical Collaboration in Nuclear Theory on Heavy-Flavor Theory (HEFTY) for QCD Matter under award no. DE-SC0023547
References
- [1] R. Rapp, D. Blaschke, and P. Crochet Prog. Part. Nucl. Phys. 65 (2010) 209–266, arXiv:0807.2470 [hep-ph].
- [2] P. Braun-Munzinger and J. Stachel Landolt-Bornstein 23 (2010) 424, arXiv:0901.2500 [nucl-th].
- [3] L. Kluberg and H. Satz Landolt-Börnstein I/23 (2010) , arXiv:0901.3831 [hep-ph].
- [4] A. Mocsy, P. Petreczky, and M. Strickland Int. J. Mod. Phys. A 28 (2013) 1340012, arXiv:1302.2180 [hep-ph].
- [5] Y. Liu, K. Zhou, and P. Zhuang Int. J. Mod. Phys. E 24 no. 11, (2015) 1530015.
- [6] S. Y. F. Liu and R. Rapp Eur. Phys. J. A 56 no. 2, (2020) 44, arXiv:1612.09138 [nucl-th].
- [7] T. Matsui and H. Satz Phys. Lett. B 178 (1986) 416–422.
- [8] P. Braun-Munzinger and J. Stachel Phys. Lett. B 490 (2000) 196–202, arXiv:nucl-th/0007059.
- [9] R. L. Thews, M. Schroedter, and J. Rafelski Phys. Rev. C 63 (2001) 054905, arXiv:hep-ph/0007323.
- [10] L. Grandchamp and R. Rapp Phys. Lett. B 523 (2001) 60–66, arXiv:hep-ph/0103124.
- [11] R. Rapp and X. Du Nucl. Phys. A 967 (2017) 216–224, arXiv:1704.07923 [hep-ph].
- [12] A. Andronic, P. Braun-Munzinger, K. Redlich, and J. Stachel Nature 561 no. 7723, (2018) 321–330, arXiv:1710.09425 [nucl-th].
- [13] A. Beraudo et al. Nucl. Phys. A 979 (2018) 21–86, arXiv:1803.03824 [nucl-th].
- [14] Z. Citron et al. CERN Yellow Rep. Monogr. 7 (2019) 1159–1410, arXiv:1812.06772 [hep-ph].
- [15] M. He, B. Wu, and R. Rapp Phys. Rev. Lett. 128 no. 16, (2022) 162301, arXiv:2111.13528 [nucl-th].
- [16] ALICE Collaboration, S. Acharya et al. Phys. Lett. B 780 (2018) 7–20, arXiv:1709.06807 [nucl-ex].
- [17] X. Du and R. Rapp JHEP 03 (2019) 015, arXiv:1808.10014 [nucl-th].
- [18] X. Du, R. Rapp, and M. He Phys. Rev. C 96 no. 5, (2017) 054901, arXiv:1706.08670 [hep-ph].
- [19] N. Brambilla, M. A. Escobedo, A. Islam, M. Strickland, A. Tiwari, A. Vairo, and P. Vander Griend arXiv:2302.11826 [hep-ph].
- [20] K. Zhou, N. Xu, and P. Zhuang Nucl. Phys. A 931 (2014) 654–658, arXiv:1408.3900 [hep-ph].
- [21] B. Krouppa, A. Rothkopf, and M. Strickland Phys. Rev. D 97 no. 1, (2018) 016017, arXiv:1710.02319 [hep-ph].
- [22] J. Hoelck, F. Nendzig, and G. Wolschin Phys. Rev. C 95 no. 2, (2017) 024905, arXiv:1602.00019 [hep-ph].
- [23] X. Du, S. Y. F. Liu, and R. Rapp Phys. Lett. B 796 (2019) 20–25, arXiv:1904.00113 [nucl-th].
- [24] ALICE Collaboration, S. Acharya et al. Phys. Rev. Lett. 123 no. 19, (2019) 192301, arXiv:1907.03169 [nucl-ex].
- [25] CMS Collaboration, A. M. Sirunyan et al. Phys. Lett. B 819 (2021) 136385, arXiv:2006.07707 [hep-ex].
- [26] M. He, H. van Hees, and R. Rapp Prog. Part. Nucl. Phys. 130 (2023) 104020, arXiv:2204.09299 [hep-ph].
- [27] J.-P. Blaizot, D. De Boni, P. Faccioli, and G. Garberoglio Nucl. Phys. A 946 (2016) 49–88, arXiv:1503.03857 [nucl-th].
- [28] R. Katz and P. B. Gossiaux Annals Phys. 368 (2016) 267–295, arXiv:1504.08087 [quant-ph].
- [29] S. Kajimoto, Y. Akamatsu, M. Asakawa, and A. Rothkopf Phys. Rev. D 97 no. 1, (2018) 014003, arXiv:1705.03365 [nucl-th].
- [30] Y. Akamatsu Phys. Rev. D 91 no. 5, (2015) 056002, arXiv:1403.5783 [hep-ph].
- [31] N. Brambilla, M. A. Escobedo, J. Soto, and A. Vairo Phys. Rev. D 96 no. 3, (2017) 034021, arXiv:1612.07248 [hep-ph].
- [32] N. Brambilla, M. A. Escobedo, A. Islam, M. Strickland, A. Tiwari, A. Vairo, and P. Vander Griend JHEP 08 (2022) 303, arXiv:2205.10289 [hep-ph].
- [33] Y. Akamatsu Prog. Part. Nucl. Phys. 123 (2022) 103932, arXiv:2009.10559 [nucl-th].
- [34] X. Zhao and R. Rapp Phys. Rev. C 82 (2010) 064905, arXiv:1008.5328 [hep-ph].
- [35] S. Aronson, E. Borras, B. Odegard, R. Sharma, and I. Vitev Phys. Lett. B 778 (2018) 384–391, arXiv:1709.02372 [hep-ph].
- [36] T. Song, K. C. Han, and C. M. Ko Phys. Rev. C 85 (2012) 054905, arXiv:1203.2964 [nucl-th].
- [37] X. Yao and B. Müller Phys. Rev. C 97 no. 1, (2018) 014908, arXiv:1709.03529 [hep-ph]. [Erratum: Phys.Rev.C 97, 049903 (2018)].
- [38] X. Du and R. Rapp Phys. Lett. B 834 (2022) 137414, arXiv:2207.00065 [nucl-th].
- [39] L. Grandchamp and R. Rapp Nucl. Phys. A 709 (2002) 415–439, arXiv:hep-ph/0205305.
- [40] X. Zhao and R. Rapp Phys. Lett. B 664 (2008) 253–257, arXiv:0712.2407 [hep-ph].
- [41] NA50 Collaboration, B. Alessandro et al. Eur. Phys. J. C 39 (2005) 335–345, arXiv:hep-ex/0412036.
- [42] NA60 Collaboration, R. Arnaldi J. Phys. G 35 (2008) 104133.
- [43] PHENIX Collaboration, A. Adare et al. Phys. Rev. Lett. 98 (2007) 232301, arXiv:nucl-ex/0611020.
- [44] PHENIX Collaboration, A. Adare et al. Phys. Rev. C 84 (2011) 054912, arXiv:1103.6269 [nucl-ex].
- [45] STAR Collaboration, L. Adamczyk et al. Phys. Rev. C 90 no. 2, (2014) 024906, arXiv:1310.3563 [nucl-ex].
- [46] STAR Collaboration, J. Adam et al. Phys. Lett. B 797 (2019) 134917, arXiv:1905.13669 [nucl-ex].
- [47] ALICE Collaboration arXiv:2211.04384 [nucl-ex].
- [48] PHENIX Collaboration, A. Adare et al. Phys. Rev. C 93 no. 3, (2016) 034903, arXiv:1509.05380 [nucl-ex].
- [49] PHENIX Collaboration, A. Adare et al. Phys. Rev. C 86 (2012) 064901, arXiv:1208.2251 [nucl-ex].
- [50] STAR Collaboration, L. Adamczyk et al. Phys. Lett. B 771 (2017) 13–20, arXiv:1607.07517 [hep-ex].
- [51] PHENIX Collaboration, C. Aidala et al. Phys. Rev. C 90 no. 6, (2014) 064908, arXiv:1404.1873 [nucl-ex].
- [52] ALICE Collaboration, B. Abelev et al. Phys. Rev. Lett. 109 (2012) 072301, arXiv:1202.1383 [hep-ex].
- [53] ALICE Collaboration, B. B. Abelev et al. Phys. Lett. B 734 (2014) 314–327, arXiv:1311.0214 [nucl-ex].
- [54] ALICE Collaboration, S. Acharya et al. Phys. Lett. B 849 (2024) 138451, arXiv:2303.13361 [nucl-ex].
- [55] STAR Collaboration, L. Adamczyk et al. Phys. Rev. Lett. 111 no. 5, (2013) 052301, arXiv:1212.3304 [nucl-ex].
- [56] ALICE Collaboration, S. Acharya et al. Phys. Rev. Lett. 119 no. 24, (2017) 242301, arXiv:1709.05260 [nucl-ex].
- [57] CMS Collaboration, A. Tumasyan et al. JHEP 10 (2023) 115, arXiv:2305.16928 [hep-ex].
- [58] NA51 Collaboration, M. C. Abreu et al. Phys. Lett. B 438 (1998) 35–40.
- [59] HERA-B Collaboration, I. Abt et al. Eur. Phys. J. C 49 (2007) 545–558, arXiv:hep-ex/0607046 [hep-ex].
- [60] PHENIX Collaboration, A. Adare et al. Phys. Rev. C 95 no. 3, (2017) 034904, arXiv:1609.06550 [nucl-ex].
- [61] PHENIX Collaboration, U. A. Acharya et al. Phys. Rev. D 101 no. 5, (2020) 052006, arXiv:1912.13424 [hep-ex].
- [62] LHCb Collaboration, R. Aaij et al. Eur. Phys. J. C 71 (2011) 1645, arXiv:1103.0423 [hep-ex].
- [63] LHCb Collaboration, R. Aaij et al. Eur. Phys. J. C 72 (2012) 2100, arXiv:1204.1258 [hep-ex].
- [64] ALICE Collaboration, S. Acharya et al. Eur. Phys. J. C 77 no. 6, (2017) 392, arXiv:1702.00557 [hep-ex].
- [65] NA50 Collaboration Collaboration, B. Alessandro et al. Eur. Phys. J. C 49 (2007) 559–567, arXiv:nucl-ex/0612013 [nucl-ex].
- [66] ALICE Collaboration, S. Acharya et al. Phys. Rev. Lett. 132 no. 4, (2024) 042301, arXiv:2210.08893 [nucl-ex].
- [67] CMS Collaboration, S. Chatrchyan et al. Phys. Rev. Lett. 109 (2012) 222301, arXiv:1208.2826 [nucl-ex]. [Erratum: Phys.Rev.Lett. 120, 199903 (2018)].
- [68] CMS Collaboration, A. M. Sirunyan et al. Phys. Rev. Lett. 120 no. 14, (2018) 142301, arXiv:1706.05984 [hep-ex].
- [69] CMS Collaboration, A. M. Sirunyan et al. Phys. Lett. B 790 (2019) 270–293, arXiv:1805.09215 [hep-ex].
- [70] CMS Collaboration, A. Tumasyan et al. arXiv:2303.17026 [hep-ex].
- [71] ATLAS Collaboration, G. Aad et al. Phys. Rev. C 107 no. 5, (2023) 054912, arXiv:2205.03042 [nucl-ex].
- [72] ALICE Collaboration, S. Acharya et al. Phys. Lett. B 822 (2021) 136579, arXiv:2011.05758 [nucl-ex].
- [73] STAR Collaboration, B. Aboona et al. Phys. Rev. Lett. 130 no. 11, (2023) 112301, arXiv:2207.06568 [nucl-ex].
- [74] P. Moreau, O. Soloveva, L. Oliva, T. Song, W. Cassing, and E. Bratkovskaya Phys. Rev. C 100 no. 1, (2019) 014911, arXiv:1903.10257 [nucl-th].
- [75] F. Riek and R. Rapp Phys. Rev. C 82 (2010) 035201, arXiv:1005.0769 [hep-ph].
- [76] S. Y. F. Liu and R. Rapp Phys. Rev. C 97 no. 3, (2018) 034918, arXiv:1711.03282 [nucl-th].
- [77] A. Capella and E. G. Ferreiro Eur. Phys. J. C 42 (2005) 419–424, arXiv:hep-ph/0505032.
- [78] C. Spieles, R. Vogt, L. Gerland, S. A. Bass, M. Bleicher, H. Stoecker, and W. Greiner Phys. Rev. C 60 (1999) 054901, arXiv:hep-ph/9902337.
- [79] B. Zhang, C. M. Ko, B.-A. Li, Z.-W. Lin, and S. Pal Phys. Rev. C 65 (2002) 054909, arXiv:nucl-th/0201038.
- [80] L. Yan, P. Zhuang, and N. Xu Phys. Rev. Lett. 97 (2006) 232301, arXiv:nucl-th/0608010.
- [81] L. Grandchamp, R. Rapp, and G. E. Brown Phys. Rev. Lett. 92 (2004) 212301.
- [82] T. Song, K. C. Han, and C. M. Ko Phys. Rev. C 84 (2011) 034907, arXiv:1103.6197 [nucl-th].
- [83] Y. Liu, B. Chen, N. Xu, and P. Zhuang Phys. Lett. B 697 (2011) 32–36, arXiv:1009.2585 [nucl-th].
- [84] T. Miura, Y. Akamatsu, M. Asakawa, and A. Rothkopf Phys. Rev. D 101 no. 3, (2020) 034011, arXiv:1908.06293 [nucl-th].
- [85] J.-P. Blaizot and M. A. Escobedo JHEP 06 (2018) 034, arXiv:1711.10812 [hep-ph].
- [86] D. Y. A. Villar, J. Zhao, J. Aichelin, and P. B. Gossiaux Phys. Rev. C 107 no. 5, (2023) 054913, arXiv:2206.01308 [nucl-th].
- [87] A. Bazavov, P. Petreczky, and A. Velytsky, Quarkonium at Finite Temperature, pp. 61–110. World Scientific, 2010. arXiv:0904.1748 [hep-ph].
- [88] G. Aarts, S. Kim, M. P. Lombardo, M. B. Oktay, S. M. Ryan, D. K. Sinclair, and J. I. Skullerud Phys. Rev. Lett. 106 (2011) 061602, arXiv:1010.3725 [hep-lat].
- [89] G. Aarts, C. Allton, S. Kim, M. P. Lombardo, M. B. Oktay, S. M. Ryan, D. K. Sinclair, and J. I. Skullerud JHEP 11 (2011) 103, arXiv:1109.4496 [hep-lat].
- [90] S. Kim, P. Petreczky, and A. Rothkopf Phys. Rev. D 91 (2015) 054511, arXiv:1409.3630 [hep-lat].
- [91] R. Larsen, S. Meinel, S. Mukherjee, and P. Petreczky Phys. Lett. B 800 (2020) 135119, arXiv:1910.07374 [hep-lat].
- [92] R. Larsen, S. Meinel, S. Mukherjee, and P. Petreczky Phys. Rev. D 100 no. 7, (2019) 074506, arXiv:1908.08437 [hep-lat].
- [93] R. Larsen, S. Meinel, S. Mukherjee, and P. Petreczky Phys. Rev. D 102 (2020) 114508, arXiv:2008.00100 [hep-lat].
- [94] S. Caron-Huot, M. Laine, and G. D. Moore JHEP 04 (2009) 053, arXiv:0901.1195 [hep-lat].
- [95] A. Bouttefeux and M. Laine JHEP 12 (2020) 150, arXiv:2010.07316 [hep-ph].
- [96] M. Laine JHEP 06 (2021) 139, arXiv:2103.14270 [hep-ph].
- [97] A. Mocsy and P. Petreczky Phys. Rev. D 77 (2008) 014501, arXiv:0705.2559 [hep-ph].
- [98] P. Petreczky Eur. Phys. J. C 62 (2009) 85–93, arXiv:0810.0258 [hep-lat].
- [99] S. Kim, P. Petreczky, and A. Rothkopf JHEP 11 (2018) 088, arXiv:1808.08781 [hep-lat].
- [100] F. Karsch, E. Laermann, S. Mukherjee, and P. Petreczky Phys. Rev. D 85 (2012) 114501, arXiv:1203.3770 [hep-lat].
- [101] A. Bazavov, F. Karsch, Y. Maezawa, S. Mukherjee, and P. Petreczky Phys. Rev. D 91 no. 5, (2015) 054503, arXiv:1411.3018 [hep-lat].
- [102] P. Petreczky, S. Sharma, and J. H. Weber Phys. Rev. D 104 no. 5, (2021) 054511, arXiv:2107.11368 [hep-lat].
- [103] TUMQCD Collaboration, A. Bazavov, N. Brambilla, P. Petreczky, A. Vairo, and J. H. Weber Phys. Rev. D 98 no. 5, (2018) 054511, arXiv:1804.10600 [hep-lat].
- [104] O. Kaczmarek and F. Zantow Phys. Rev. D 71 (2005) 114510, arXiv:hep-lat/0503017.
- [105] P. Petreczky and K. Petrov Phys. Rev. D 70 (2004) 054503, arXiv:hep-lat/0405009.
- [106] A. Rothkopf, T. Hatsuda, and S. Sasaki Phys. Rev. Lett. 108 (2012) 162001, arXiv:1108.1579 [hep-lat].
- [107] Y. Burnier, O. Kaczmarek, and A. Rothkopf Phys. Rev. Lett. 114 no. 8, (2015) 082001, arXiv:1410.2546 [hep-lat].
- [108] D. Lafferty and A. Rothkopf Phys. Rev. D 101 no. 5, (2020) 056010, arXiv:1906.00035 [hep-ph].
- [109] HotQCD Collaboration, D. Bala, O. Kaczmarek, R. Larsen, S. Mukherjee, G. Parkar, P. Petreczky, A. Rothkopf, and J. H. Weber Phys. Rev. D 105 no. 5, (2022) 054513, arXiv:2110.11659 [hep-lat].
- [110] N. Haque, A. Bandyopadhyay, J. O. Andersen, M. G. Mustafa, M. Strickland, and N. Su JHEP 05 (2014) 027, arXiv:1402.6907 [hep-ph].
- [111] N. Haque and M. Strickland Phys. Rev. C 103 no. 3, (2021) 031901, arXiv:2011.06938 [hep-ph].
- [112] A. Bazavov, D. Hoying, O. Kaczmarek, R. N. Larsen, S. Mukherjee, P. Petreczky, A. Rothkopf, and J. H. Weber arXiv:2308.16587 [hep-lat].
- [113] Z. Tang, S. Mukherjee, P. Petreczky, and R. Rapp arXiv:2310.18864 [hep-lat].
- [114] D. Banerjee, S. Datta, R. Gavai, and P. Majumdar Phys. Rev. D 85 (2012) 014510, arXiv:1109.5738 [hep-lat].
- [115] A. Francis, O. Kaczmarek, M. Laine, T. Neuhaus, and H. Ohno Phys. Rev. D 92 no. 11, (2015) 116003, arXiv:1508.04543 [hep-lat].
- [116] N. Brambilla, V. Leino, P. Petreczky, and A. Vairo Phys. Rev. D 102 no. 7, (2020) 074503, arXiv:2007.10078 [hep-lat].
- [117] L. Altenkort, A. M. Eller, O. Kaczmarek, L. Mazur, G. D. Moore, and H.-T. Shu Phys. Rev. D 103 no. 1, (2021) 014511, arXiv:2009.13553 [hep-lat].
- [118] D. Banerjee, R. Gavai, S. Datta, and P. Majumdar arXiv:2206.15471 [hep-ph].
- [119] TUMQCD Collaboration, N. Brambilla, V. Leino, J. Mayer-Steudte, and P. Petreczky Phys. Rev. D 107 no. 5, (2023) 054508, arXiv:2206.02861 [hep-lat].
- [120] HotQCD Collaboration, L. Altenkort, O. Kaczmarek, R. Larsen, S. Mukherjee, P. Petreczky, H.-T. Shu, and S. Stendebach Phys. Rev. Lett. 130 no. 23, (2023) 231902, arXiv:2302.08501 [hep-lat].
- [121] L. Altenkort, D. de la Cruz, O. Kaczmarek, R. Larsen, G. D. Moore, S. Mukherjee, P. Petreczky, H.-T. Shu, and S. Stendebach arXiv:2311.01525 [hep-lat].
- [122] N. Brambilla, M. A. Escobedo, J. Soto, and A. Vairo Phys. Rev. D 97 no. 7, (2018) 074009, arXiv:1711.04515 [hep-ph].
- [123] RBC-Bielefeld Collaboration, P. Petreczky, P. Hegde, and A. Velytsky PoS LAT2009 (2009) 159, arXiv:0911.0196 [hep-lat].
- [124] A. Bazavov et al. Phys. Lett. B 737 (2014) 210–215, arXiv:1404.4043 [hep-lat].
- [125] R. Bellwied, S. Borsanyi, Z. Fodor, S. D. Katz, A. Pasztor, C. Ratti, and K. K. Szabo Phys. Rev. D 92 no. 11, (2015) 114505, arXiv:1507.04627 [hep-lat].
- [126] S. Mukherjee, P. Petreczky, and S. Sharma Phys. Rev. D 93 no. 1, (2016) 014502, arXiv:1509.08887 [hep-lat].
- [127] A. Bazavov, D. Bollweg, O. Kaczmarek, F. Karsch, S. Mukherjee, P. Petreczky, C. Schmidt, and S. Sharma arXiv:2312.12857 [hep-lat].
- [128] X. Yao, W. Ke, Y. Xu, S. A. Bass, and B. Müller JHEP 01 (2021) 046, arXiv:2004.06746 [hep-ph].
- [129] X. Yao and T. Mehen Phys. Rev. D 99 no. 9, (2019) 096028, arXiv:1811.07027 [hep-ph].
- [130] X. Yao and T. Mehen JHEP 02 (2021) 062, arXiv:2009.02408 [hep-ph].
- [131] X. Yao Int. J. Mod. Phys. A 36 no. 20, (2021) 2130010, arXiv:2102.01736 [hep-ph].
- [132] T. Binder, K. Mukaida, B. Scheihing-Hitschfeld, and X. Yao JHEP 01 (2022) 137, arXiv:2107.03945 [hep-ph].
- [133] G. Nijs, B. Scheihing-Hitschfeld, and X. Yao JHEP 06 (2023) 007, arXiv:2304.03298 [hep-ph].
- [134] G. Nijs, B. Scheihing-Hitschfeld, and X. Yao arXiv:2310.09325 [hep-ph].
- [135] B. Scheihing-Hitschfeld and X. Yao Phys. Rev. D 108 no. 5, (2023) 054024, arXiv:2306.13127 [hep-ph].
- [136] J. S. Moreland, J. E. Bernhard, and S. A. Bass Phys. Rev. C 92 no. 1, (2015) 011901, arXiv:1412.4708 [nucl-th].
- [137] W. Ke, Y. Xu, and S. A. Bass Phys. Rev. C 98 no. 6, (2018) 064901, arXiv:1806.08848 [nucl-th].
- [138] K. J. Eskola, P. Paakkinen, H. Paukkunen, and C. A. Salgado Eur. Phys. J. C 77 no. 3, (2017) 163, arXiv:1612.05741 [hep-ph].
- [139] X. Yao and B. Müller Phys. Rev. D 100 no. 1, (2019) 014008, arXiv:1811.09644 [hep-ph].
- [140] X. Yao, Application of Effective Field Theory in Nuclear Physics. PhD thesis, Duke U., 2019. arXiv:1911.08500 [nucl-th].
- [141] B. Scheihing-Hitschfeld and X. Yao Phys. Rev. Lett. 130 no. 5, (2023) 052302, arXiv:2205.04477 [hep-ph].
- [142] J. E. Bernhard, J. S. Moreland, S. A. Bass, J. Liu, and U. Heinz Phys. Rev. C 94 no. 2, (2016) 024907, arXiv:1605.03954 [nucl-th].
- [143] N. Brambilla, A. Pineda, J. Soto, and A. Vairo Nucl. Phys. B 566 (2000) 275, arXiv:hep-ph/9907240.
- [144] N. Brambilla, A. Pineda, J. Soto, and A. Vairo Rev. Mod. Phys. 77 (2005) 1423, arXiv:hep-ph/0410047.
- [145] Particle Data Group Collaboration, P. A. Zyla et al. PTEP 2020 no. 8, (2020) 083C01.
- [146] N. Brambilla, M. A. Escobedo, A. Vairo, and P. Vander Griend Phys. Rev. D 100 no. 5, (2019) 054025, arXiv:1903.08063 [hep-ph].
- [147] N. Brambilla, J. Ghiglieri, A. Vairo, and P. Petreczky Phys. Rev. D 78 (2008) 014017, arXiv:0804.0993 [hep-ph].
- [148] N. Brambilla, M. A. Escobedo, J. Ghiglieri, and A. Vairo JHEP 12 (2011) 116, arXiv:1109.5826 [hep-ph].
- [149] N. Brambilla, M. A. Escobedo, J. Ghiglieri, and A. Vairo JHEP 05 (2013) 130, arXiv:1303.6097 [hep-ph].
- [150] N. Brambilla, M. A. Escobedo, M. Strickland, A. Vairo, P. Vander Griend, and J. H. Weber JHEP 05 (2021) 136, arXiv:2012.01240 [hep-ph].
- [151] N. Brambilla, M. A. Escobedo, M. Strickland, A. Vairo, P. Vander Griend, and J. H. Weber Phys. Rev. D 104 no. 9, (2021) 094049, arXiv:2107.06222 [hep-ph].
- [152] M. Alqahtani, M. Nopoush, and M. Strickland Phys. Rev. C 92 no. 5, (2015) 054910, arXiv:1509.02913 [hep-ph].
- [153] M. Alqahtani, M. Nopoush, and M. Strickland Phys. Rev. C 95 no. 3, (2017) 034906, arXiv:1605.02101 [nucl-th].
- [154] M. Alqahtani, M. Nopoush, and M. Strickland Prog. Part. Nucl. Phys. 101 (2018) 204–248, arXiv:1712.03282 [nucl-th].
- [155] M. Alqahtani, M. Nopoush, R. Ryblewski, and M. Strickland Phys. Rev. Lett. 119 no. 4, (2017) 042301, arXiv:1703.05808 [nucl-th].
- [156] M. Alqahtani, M. Nopoush, R. Ryblewski, and M. Strickland Phys. Rev. C 96 no. 4, (2017) 044910, arXiv:1705.10191 [nucl-th].
- [157] M. Nopoush, R. Ryblewski, and M. Strickland Phys. Rev. C 90 no. 1, (2014) 014908, arXiv:1405.1355 [hep-ph].
- [158] D. Almaalol, M. Alqahtani, and M. Strickland Phys. Rev. C 99 no. 4, (2019) 044902, arXiv:1807.04337 [nucl-th].
- [159] M. Alqahtani and M. Strickland Phys. Rev. C 102 no. 6, (2020) 064902, arXiv:2007.04209 [nucl-th].
- [160] M. Alqahtani and M. Strickland Eur. Phys. J. C 81 no. 11, (2021) 1022, arXiv:2008.07657 [nucl-th].
- [161] A. Bazavov J. Phys. Conf. Ser. 446 (2013) 012011, arXiv:1303.6294 [hep-lat].
- [162] ATLAS Collaboration, G. Aad et al. Phys. Rev. C 107 no. 5, (2023) 054912, arXiv:2205.03042 [nucl-ex].
- [163] CMS Collaboration, “Observation of the meson and sequential suppression of states in PbPb collisions at ,” tech. rep., CERN, Geneva, 2022. https://cds.cern.ch/record/2805926.
- [164] S. Delorme, Theoretical description of quarkonia dynamics in the Quark-Gluon Plasma with a quantum master equation approach. PhD thesis, École Nationale Supérieure Mines-Télécom Atlantique Bretagne Pays de la Loire - IMT Atlantique, 2021.
- [165] S. Delorme, T. Gousset, R. Katz, and P. B. Gossiaux Acta Phys. Pol. B Proc. Suppl. 16 (2023) 1–A112.
- [166] R. Katz, S. Delorme, and P.-B. Gossiaux Eur. Phys. J. A 58 no. 10, (2022) 198, arXiv:2205.05154 [hep-ph].
- [167] J. Zhao, J. Aichelin, P. B. Gossiaux, and K. Werner arXiv:2310.08684 [hep-ph].
- [168] E. A. Remler and A. P. Sathe Annals Phys. 91 (1975) 295–324.
- [169] E. A. Remler Annals Phys. 95 (1975) 455–495.
- [170] E. A. Remler Annals Phys. 136 (1981) 293.
- [171] M. Gyulassy, K. Frankel, and E. a. Remler Nucl. Phys. A 402 (1983) 596–611.
- [172] T. Song, J. Aichelin, and E. Bratkovskaya Phys. Rev. C 107 no. 5, (2023) 054906, arXiv:2302.14001 [hep-ph].
- [173] T. Sjostrand, S. Mrenna, and P. Z. Skands JHEP 05 (2006) 026, arXiv:hep-ph/0603175.
- [174] M. Cacciari, S. Frixione, N. Houdeau, M. L. Mangano, P. Nason, and G. Ridolfi JHEP 10 (2012) 137, arXiv:1205.6344 [hep-ph].
- [175] T. Song, J. Aichelin, and E. Bratkovskaya Phys. Rev. C 96 no. 1, (2017) 014907, arXiv:1705.00046 [nucl-th].
- [176] K. J. Eskola, H. Paukkunen, and C. A. Salgado JHEP 04 (2009) 065, arXiv:0902.4154 [hep-ph].
- [177] T. Song, J. Aichelin, J. Zhao, P. B. Gossiaux, and E. Bratkovskaya Phys. Rev. C 108 no. 5, (2023) 054908, arXiv:2305.10750 [nucl-th].
- [178] H. Berrehrah, E. Bratkovskaya, W. Cassing, P. B. Gossiaux, J. Aichelin, and M. Bleicher Phys. Rev. C 89 no. 5, (2014) 054901, arXiv:1308.5148 [hep-ph].
- [179] T. Song, H. Berrehrah, D. Cabrera, J. M. Torres-Rincon, L. Tolos, W. Cassing, and E. Bratkovskaya Phys. Rev. C 92 no. 1, (2015) 014910, arXiv:1503.03039 [nucl-th].
- [180] T. Song, H. Berrehrah, J. M. Torres-Rincon, L. Tolos, D. Cabrera, W. Cassing, and E. Bratkovskaya Phys. Rev. C 96 no. 1, (2017) 014905, arXiv:1605.07887 [nucl-th].
- [181] T. Song, P. Moreau, J. Aichelin, and E. Bratkovskaya Phys. Rev. C 101 no. 4, (2020) 044901, arXiv:1910.09889 [nucl-th].
- [182] E. L. Bratkovskaya, W. Cassing, and H. Stoecker Phys. Rev. C 67 (2003) 054905, arXiv:nucl-th/0301083.
- [183] M. Cacciari, P. Nason, and R. Vogt Phys. Rev. Lett. 95 (2005) 122001, arXiv:hep-ph/0502203.
- [184] T. Song, H. Berrehrah, D. Cabrera, W. Cassing, and E. Bratkovskaya Phys. Rev. C 93 no. 3, (2016) 034906, arXiv:1512.00891 [nucl-th].
- [185] W. Cassing and E. L. Bratkovskaya Phys. Rev. C 78 (2008) 034919, arXiv:0808.0022 [hep-ph].
- [186] W. Cassing Eur. Phys. J. ST 168 (2009) 3–87, arXiv:0808.0715 [nucl-th].
- [187] W. Cassing and E. L. Bratkovskaya Nucl. Phys. A 831 (2009) 215–242, arXiv:0907.5331 [nucl-th].
- [188] E. L. Bratkovskaya, W. Cassing, V. P. Konchakovski, and O. Linnyk Nucl. Phys. A 856 (2011) 162–182, arXiv:1101.5793 [nucl-th].
- [189] O. Linnyk, E. L. Bratkovskaya, and W. Cassing Prog. Part. Nucl. Phys. 87 (2016) 50–115, arXiv:1512.08126 [nucl-th].
- [190] J.-P. Blaizot and M. A. Escobedo Phys. Rev. D 98 no. 7, (2018) 074007, arXiv:1803.07996 [hep-ph].
- [191] J.-P. Blaizot and M. A. Escobedo Phys. Rev. D 104 no. 5, (2021) 054034, arXiv:2106.15371 [hep-ph].
- [192] A. Capella, A. Kaidalov, A. Kouider Akil, and C. Gerschel Phys. Lett. B 393 (1997) 431–436, arXiv:hep-ph/9607265.
- [193] N. Armesto and A. Capella Phys. Lett. B 430 (1998) 23–31, arXiv:hep-ph/9705275.
- [194] N. Armesto, A. Capella, and E. G. Ferreiro Phys. Rev. C 59 (1999) 395–404, arXiv:hep-ph/9807258.
- [195] A. Capella, E. G. Ferreiro, and A. B. Kaidalov Phys. Rev. Lett. 85 (2000) 2080–2083, arXiv:hep-ph/0002300.
- [196] E. G. Ferreiro and J.-P. Lansberg JHEP 10 (2018) 094, arXiv:1804.04474 [hep-ph]. [Erratum: JHEP 03, 063 (2019)].
- [197] A. Capella, L. Bravina, E. G. Ferreiro, A. B. Kaidalov, K. Tywoniuk, and E. Zabrodin Eur. Phys. J. C 58 (2008) 437–444, arXiv:0712.4331 [hep-ph].
- [198] E. G. Ferreiro Phys. Lett. B 731 (2014) 57–63, arXiv:1210.3209 [hep-ph].
- [199] E. G. Ferreiro Phys. Lett. B 749 (2015) 98–103, arXiv:1411.0549 [hep-ph].
- [200] P. Braun-Munzinger and J. Stachel Phys. Lett. B 490 (2000) 196–202, arXiv:nucl-th/0007059.
- [201] A. Andronic, P. Braun-Munzinger, K. Redlich, and J. Stachel Nucl. Phys. A 789 (2007) 334–356, arXiv:nucl-th/0611023.
- [202] A. Andronic, P. Braun-Munzinger, M. K. Köhler, A. Mazeliauskas, K. Redlich, J. Stachel, and V. Vislavicius JHEP 07 (2021) 035, arXiv:2104.12754 [hep-ph].
- [203] A. Andronic, P. Braun-Munzinger, K. Redlich, and J. Stachel Phys. Lett. B 659 (2008) 149–155, arXiv:0708.1488 [nucl-th].
- [204] ALICE Collaboration, S. Acharya et al. JHEP 01 (2022) 174, arXiv:2110.09420 [nucl-ex].
- [205] B. Schenke, S. Jeon, and C. Gale Phys. Rev. Lett. 106 (2011) 042301, arXiv:1009.3244 [hep-ph].
- [206] S. Floerchinger, E. Grossi, and J. Lion Phys. Rev. C 100 no. 1, (2019) 014905, arXiv:1811.01870 [nucl-th].
- [207] A. Andronic, P. Braun-Munzinger, H. Brunßen, J. Crkovská, J. Stachel, V. Vislavicius, and M. Völkl arXiv:2308.14821 [hep-ph].
- [208] A. Andronic, P. Braun-Munzinger, K. Redlich, and J. Stachel Acta Phys. Polon. Supp. 16 no. 1, (2023) 107, arXiv:2209.14562 [hep-ph].
- [209] X. Zhao and R. Rapp Nucl. Phys. A 859 (2011) 114–125, arXiv:1102.2194 [hep-ph].
- [210] X. Du and R. Rapp Nucl. Phys. A 943 (2015) 147–158, arXiv:1504.00670 [hep-ph].
- [211] L. Grandchamp, S. Lumpkins, D. Sun, H. van Hees, and R. Rapp Phys. Rev. C 73 (2006) 064906, arXiv:hep-ph/0507314.
- [212] B. Wu, Z. Tang, M. He, and R. Rapp Phys. Rev. C 109 no. 1, (2024) 014906, arXiv:2302.11511 [nucl-th].
- [213] B. Wu, X. Du, M. Sibila, and R. Rapp Eur. Phys. J. A 57 no. 4, (2021) 122, arXiv:2006.09945 [nucl-th]. [Erratum: Eur.Phys.J.A 57, 314 (2021)].
- [214] F. Riek and R. Rapp New J. Phys. 13 (2011) 045007, arXiv:1012.0019 [nucl-th].
- [215] M. He, R. J. Fries, and R. Rapp Phys. Rev. C 85 (2012) 044911, arXiv:1112.5894 [nucl-th].
- [216] S. Y. F. Liu and R. Rapp Phys. Rev. C 106 no. 5, (2022) 055201, arXiv:2111.13620 [hep-ph].
- [217] A. Andronic et al. Eur. Phys. J. C 76 no. 3, (2016) 107, arXiv:1506.03981 [nucl-ex].
- [218] K. Zhou, N. Xu, Z. Xu, and P. Zhuang Phys. Rev. C 89 no. 5, (2014) 054911, arXiv:1401.5845 [nucl-th].
- [219] B. Chen Chin. Phys. C 43 no. 12, (2019) 124101, arXiv:1811.11393 [nucl-th].
- [220] B. Chen, T. Guo, Y. Liu, and P. Zhuang Phys. Lett. B 765 (2017) 323–327, arXiv:1607.07927 [nucl-th].
- [221] J. Zhao and P. Zhuang Phys. Rev. C 105 no. 6, (2022) 064907.
- [222] J. Zhao, K. Zhou, S. Chen, and P. Zhuang Prog. Part. Nucl. Phys. 114 (2020) 103801, arXiv:2005.08277 [nucl-th].
- [223] B. Chen, L. Jiang, X.-H. Liu, Y. Liu, and J. Zhao Phys. Rev. C 105 no. 5, (2022) 054901, arXiv:2107.00969 [hep-ph].
- [224] M. E. Peskin Nucl. Phys. B 156 (1979) 365–390.
- [225] G. Bhanot and M. E. Peskin Nucl. Phys. B 156 (1979) 391–416.
- [226] P. Petreczky J. Phys. G 37 (2010) 094009, arXiv:1001.5284 [hep-ph].
- [227] Z.-w. Lin and C. M. Ko Phys. Rev. C 62 (2000) 034903, arXiv:nucl-th/9912046.
- [228] J. Huefner and P.-f. Zhuang Phys. Lett. B 559 (2003) 193–200, arXiv:nucl-th/0208004.
- [229] P.-f. Zhuang and X. L. Zhu Phys. Rev. C 67 (2003) 067901, arXiv:nucl-th/0304019.
- [230] D. Kharzeev and M. Nardi Phys. Lett. B 507 (2001) 121–128, arXiv:nucl-th/0012025.
- [231] Particle Data Group Collaboration, M. Tanabashi et al. Phys. Rev. D 98 no. 3, (2018) 030001.
- [232] D. C. McGlinchey, A. D. Frawley, and R. Vogt Phys. Rev. C 87 no. 5, (2013) 054910, arXiv:1208.2667 [nucl-th].
- [233] R. Vogt Phys. Rev. C 106 no. 2, (2022) 025201, arXiv:2207.04347 [hep-ph].
- [234] F. Arleo, R. Kolevatov, S. Peigné, and M. Rustamova JHEP 05 (2013) 155, arXiv:1304.0901 [hep-ph].
- [235] C. Lourenco, R. Vogt, and H. K. Woehri JHEP 02 (2009) 014, arXiv:0901.3054 [hep-ph].
- [236] J. Ftacnik, P. Lichard, and J. Pisut Phys. Lett. B 207 (1988) 194–198.
- [237] S. Gavin, M. Gyulassy, and A. Jackson Phys. Lett. B 207 (1988) 257–262.
- [238] R. Vogt, M. Prakash, P. Koch, and T. H. Hansson Phys. Lett. B 207 (1988) 263–268.
- [239] S. Gavin and R. Vogt Nucl. Phys. B 345 (1990) 104–124.
- [240] Y.-Q. Ma and R. Vogt Phys. Rev. D 94 no. 11, (2016) 114029, arXiv:1609.06042 [hep-ph].
- [241] R. L. Thews and M. L. Mangano Phys. Rev. C 73 (2006) 014904, arXiv:nucl-th/0505055.
- [242] R. E. Nelson, R. Vogt, and A. D. Frawley Phys. Rev. C 87 no. 1, (2013) 014908, arXiv:1210.4610 [hep-ph].
- [243] R. Vogt Phys. Rev. C 92 no. 3, (2015) 034909, arXiv:1507.04418 [hep-ph].
- [244] V. Cheung and R. Vogt Phys. Rev. D 98 no. 11, (2018) 114029, arXiv:1808.02909 [hep-ph].
- [245] V. Cheung and R. Vogt Phys. Rev. D 99 no. 3, (2019) 034007, arXiv:1811.11570 [hep-ph].
- [246] V. Cheung and R. Vogt Phys. Rev. D 104 no. 9, (2021) 094026, arXiv:2102.09118 [hep-ph].
- [247] V. Cheung and R. Vogt Phys. Rev. D 95 no. 7, (2017) 074021, arXiv:1702.07809 [hep-ph].
- [248] V. Cheung and R. Vogt Phys. Rev. D 96 no. 5, (2017) 054014, arXiv:1706.07686 [hep-ph].
- [249] V. Cheung and R. Vogt Phys. Rev. C 105 no. 5, (2022) 055202, arXiv:2203.10154 [hep-ph].
- [250] N. Brambilla et al. Eur. Phys. J. C 74 no. 10, (2014) 2981, arXiv:1404.3723 [hep-ph].
- [251] Y.-Q. Ma and R. Venugopalan Phys. Rev. Lett. 113 no. 19, (2014) 192301, arXiv:1408.4075 [hep-ph].
- [252] F. Maltoni et al. Phys. Lett. B 638 (2006) 202–208, arXiv:hep-ph/0601203.
- [253] Y. Feng, J.-P. Lansberg, and J.-X. Wang Eur. Phys. J. C 75 no. 7, (2015) 313, arXiv:1504.00317 [hep-ph].
- [254] HotQCD Collaboration, A. Bazavov et al. Phys. Lett. B 795 (2019) 15–21, arXiv:1812.08235 [hep-lat].
- [255] S. Borsanyi, Z. Fodor, J. N. Guenther, R. Kara, S. D. Katz, P. Parotto, A. Pasztor, C. Ratti, and K. K. Szabo Phys. Rev. Lett. 125 no. 5, (2020) 052001, arXiv:2002.02821 [hep-lat].
- [256] B. Schenke, S. Jeon, and C. Gale Phys. Rev. C 82 (2010) 014903, arXiv:1004.1408 [hep-ph].
- [257] T. Song, Y. Park, S. H. Lee, and C.-Y. Wong Phys. Lett. B 659 (2008) 621–627, arXiv:0709.0794 [hep-ph].