Next Article in Journal
Effective Online Knowledge Distillation via Attention-Based Model Ensembling
Previous Article in Journal
Prediction Model of Wastewater Pollutant Indicators Based on Combined Normalized Codec
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Review

Anisotropic Flows of Charmonium in the Relativistic Heavy-Ion Collisions

1
Culver Academies, Academy Rd., Culver, IN 46511, USA
2
Department of Physics, Tianjin University, Tianjin 300354, China
*
Author to whom correspondence should be addressed.
Mathematics 2022, 10(22), 4284; https://doi.org/10.3390/math10224284
Submission received: 4 October 2022 / Revised: 24 October 2022 / Accepted: 28 October 2022 / Published: 16 November 2022
(This article belongs to the Section Mathematical Physics)

Abstract

:
We review recent studies about anisotropic flows ( v 1 , v 2 , v 3 ) of charmonium in the quark-gluon plasma produced in relativistic heavy-ion collisions. Collective flows of the bulk medium are developed due to the anisotropic pressure gradient of the medium. Strongly coupled with the bulk medium, charm quarks carry collective flows from the expanding medium, which will be inherited by the regenerated charmonium via the coalescence process. In event-by-event collisions where nucleon positions fluctuate from the smooth distribution, there is triangularity in the medium initial energy density. Triangular flows of the bulk medium and heavy flavor particles can be developed due to the initial fluctuations. In the longitudinal direction, the rapidity-odd distribution of the initial energy density is induced by the rotation of the medium in non-central heavy-ion collisions. Charmonium suffers biased dissociation along positive and negative x-directions in forward (backward) rapidity. The directed flow of charmonium becomes non-zero. The directed, elliptic and triangular flows ( v 1 , v 2 , v 3 ) of charmonium come from the anisotropic initial distributions of the medium energy density in the transverse and longitudinal directions.

1. Introduction

In the relativistic heavy-ion collisions, an extremely hot deconfined medium consisting of quarks and gluons, also called “quark-gluon plasma” (QGP), is believed to be generated in the collision area [1]. Studying the signals and the properties of the deconfined matter, which is similar to the state of our early universe, is one of the important tasks in relativistic heavy-ion collisions. Deconfined matter is created when the energy density is above a critical value ∼ 1 GeV / fm 3 [2]. Then, the hot medium expands outside violently due to a large gradient of the energy density and turns into a hadronic medium with a first-order or crossover phase transition when the temperature drops below a critical temperature T c [3]. The anisotropy of the energy density profile is transformed into the anisotropy of the azimuthal distribution of light hadrons [4]. Based on the theoretical and experimental studies about the anisotropic flows of light hadrons [5,6,7,8], QGP turns out to be nearly a perfect fluid with a small ratio of shear viscosity over entropy density η / s . Heavy quarks and quarkonium, which are produced at the beginning of nuclear collisions, have been suggested to be a clean probe of an early stage of the medium [9]. Compared with light hadrons, which are usually produced at the hadronization hypersurface, heavy quarkonium is more sensitive to the initial profiles of the hot medium. The color-screening effect and random scatterings with thermal partons dissociate heavy quarkonium states in the hot medium. It suppresses the nuclear modification factor R A A , which is defined as the ratio of the quarkonium production in AA collisions and the production in proton-proton (pp) collisions scaled by the number of binary collisions N c o l l [10,11,12,13,14,15,16,17,18,19,20,21]. The suppression of R A A is proportional to the cold and hot nuclear matter effects, which have been observed in the experiments at the Super Proton Synchrotron (SPS) [22], the Relativistic Heavy-Ion Collider (RHIC) [23] and the Large Hadron Collider (LHC) [24,25].
The charmonium nuclear modification factor is enhanced in AA collisions at the LHC energies, where the medium temperatures are higher than the situation at RHIC energies. The increase in charmonium R A A at LHC energies is suggested to come from the coalescence of charm and anti-charm quarks in the QGP [26,27,28,29,30,31,32]. This coalescence contribution is proportional to the square of charm pair numbers N c c ¯ in the medium and is significantly enhanced when N c c ¯ becomes large. At LHC collision energies, N c c ¯ is much larger than the value at RHIC and SPS energies. The p T spectrum of regenerated charmonium depends on the distribution of charm quarks in phase space. In the expanding QGP, charm quarks are strongly coupled with the medium to dump energy and carry collective flows in the medium [33,34,35,36,37,38,39]. Heavy quark energy loss is determined by the coupling strength between heavy quarks and the medium. With stronger coupling strength, heavy quarks tend to reach momentum thermalization and diffuse into a larger volume with the expanding medium. Therefore, the regenerated charmonium from the kinetically thermalized charm quarks are mainly distributed at the low p T region [40]. Charm spatial diffusions with the expanding medium reduce the probability of their coalescence into a charmonium-bound state. At the LHC energies, where the regeneration mechanism dominates charmonium production, charmonium’s final distribution is closely connected with the dynamical evolutions of charm quarks in the medium.
In this manuscript, we review the recent studies about the effects of charm quark diffusions on the charmonium collective flows. To describe the dynamical evolutions of charm quarks and charmonium, the Boltzmann-type transport model [12] and the Langevin-plus-coalescence model [41] will be introduced. In the expanding medium, the diffusion of charm quarks reduces the spatial densities of charm quarks and suppresses the coalescence probability per one c and c ¯ . In semi-central collisions, the elliptic flows of the medium are developed. Additional triangular flows are also observed, which is induced by the fluctuations in the distribution of initial energy density. The rotation of the hot medium also results in the rapidity-odd distribution of the medium along the longitudinal direction. Charmonium suffers from biased dissociation along the positive and negative x−directions in the tilted medium. The directed flow of charmonium is non-zero in forward (backward) rapidities. All of charmonium anisotropic flows in AA collisions are due to the anisotropic initial distribution of the medium energy density in the transverse and longitudinal directions, as well as the charm quark diffusions in the medium.

2. Charm Thermalization and Charmonium Regeneration

2.1. Charm Quark Evolution

As heavy quarks are stongly coupled with the quark-gluon plasma, they suffer significant energy loss and carry collective flows in the medium. Large collective flows of D-mesons have been observed in experiments where D-mesons are close to kinetic equilibrium when they move out of the medium. The dynamical evolutions of charm quarks in QGP and D-mesons in hadronic gas are regarded as Brownian motion. Heavy quarks lose energy via random scatterings with thermal partons and the gluon radiation induced by the medium. Including both hot medium effects, the time evolution of heavy quark momentum is described with the Langevin equation [42],
d p d t = η p + ξ + f g ,
where p is the three-dimensional momentum of heavy quarks. The drag coefficient depends on the medium temperature T and the momentum η = κ / ( 2 T E ) with the definition of the energy E = m c 2 + | p | 2 . Charm quark mass is taken as m c = 1.5 GeV [41]. The dynamical evolutions of the medium local temperatures T ( x , t ) will be given by the hydrodynamic model. The momentum diffusion coefficient κ is connected with the spatial diffusion coefficient D s via the relation κ = 2 T 2 / D s . From theoretical and experimental studies about the elliptic flows and nuclear modification factors of D-mesons in nuclear collisions, the value of D s ( 2 π T ) is estimated to be 5 D s ( 2 π T ) 7 with weak dependence on temperature [43,44]. It becomes larger in the hadronic medium due to a weaker coupling strength.
In the noise term ξ , the correlation between random scatterings at different time steps have been neglected, and ξ is approximated to be a white noise satisfying the relation
ξ i ( t ) ξ j ( t ) = κ δ i j δ ( t t ) ,
where the momentum dependence has been neglected in ξ . Three directions are represented by the index, i , j = ( 1 , 2 , 3 ) . t and t are the time points of the evolution.
In the third term, the gluon radiation induced by the medium is represented by the term f g = d p g / d t , which contributes a recoil force on heavy quarks. p g is the three-dimensional momentum of an emitted gluon. Neglect the contribution from subsequent gluon emissions; in each time step Δ t , the number of emitted gluons N g in the time period t t + Δ t is calculated with the formula [42]
P rad ( t , Δ t ) = N g ( t , Δ t ) = Δ t d x d k T 2 d N g d x d k T 2 d t ,
where the spectrum of emitted gluons are taken from Refs. [45,46]. When the time period Δ t is small enough, the number of emitted gluons N g is always smaller than the unit, and the value of P rad can be interpreted as the probability of emitting one gluon. x = E g / E is the ratio of emitted gluons and the heavy quark energies. k T is the transverse momentum of the emitted gluons. The radiation term dominates the energy loss of heavy quarks at high p T region, while at low p T regions, heavy quark energy loss is dominated by the first two terms on the right-hand side of Equation (1). In this work, we focus on the regeneration of charmonium via the coalescence of charm and anti-charm quarks, where most regenerated charmonium is located at low p T bins.
The initial momentum distribution of heavy quarks can be calculated with the perturbative QCD model, such as the PYTHIA and FONLL model [47]. The initial spatial distribution of heavy quarks is proportional to the product of two thickness functions, d N c c ¯ / d x T = σ p p c c ¯ T A ( x T b / 2 ) T B ( x T + b / 2 ) . σ p p c c ¯ is the production cross-section of charm pairs in proton-proton collisions. T A ( B ) = d z ρ A ( B ) ( x T , z ) is the thickness function defined as the integration of the nucleon density over the longitudinal direction. When heavy quarks move in the medium according to the Langevin equation, their positions are updated at each time step, x ( t + Δ t ) = x ( t ) + p / E · Δ t . The charm quark mass is larger than the medium temperature. The thermal production of charm pairs can be neglected, and the total number of charm pairs are conserved when they undergo Brownian motion in the medium. Heavy quarks carry the information of the initial stage of the hot medium.

2.2. Hot Medium Evolution

To explain the large collective flows of light hadrons in experiments, the QGP medium is suggested to be a nearly perfect fluid. The dynamical evolution of the medium temperatures and the velocities are well described in the hydrodynamic models. The equation of state (EoS) of the QGP and hadronic medium are taken from lattice QCD calculations and the Hadron resonance gas model, respectively [48]. In Figure 1, time evolution of the temperatures at the center of the medium are given by a MUSIC calculation [49]. The collision energy is taken as s N N = 5.02 TeV. The initial energy density profile is given by the Glauber model, where the maximum initial temperature is determined with the hadron multiplicity observed in heavy-ion collisions. The evolution of the temperature stops when it decreases to a kinetic freeze-out temperature ∼ 0.12 GeV.

3. Charmonium Thermal Production

3.1. Diffusion of Charm Quarks

With multiple charm quarks undergoing Brownian motion in the medium, uncorrelated c and c ¯ , which are usually from different c c ¯ pairs, may collide with each other to form a new charmonium state. At low temperatures close to T c , heavy quark potential is partially restored, so the newly formed charmonium state may survive from the hot medium. At LHC energies, where the number of charm pairs is large, the coalescence of c and c ¯ even dominates the final production of J / ψ and ψ ( 2 S ) in Pb-Pb collisions. The coalescence probability between one random c and c ¯ quarks is written as
P c c ¯ ψ ( x M , p M ) events = g M d x 1 d x 2 d p 1 ( 2 π ) 3 d p 2 ( 2 π ) 3 d 2 N 1 norm d x 1 d p 1 d 2 N 2 norm d x 2 d p 2 W ( x r , q r ) × δ ( 3 ) ( p M p 1 p 2 ) δ ( 3 ) ( x M x 1 + x 2 2 ) ,
where x 1 , 2 and p 1 , 2 are the positions and momentum of c and c ¯ quarks. d 2 N 1 norm d x 1 d p 1 is the normalized distribution of one charm in phase space. The delta function ensures the momentum conservation in the reaction c + c ¯ ψ + g , where the gluon momentum has been neglected for simplicity, x M = ( x 1 + x 2 ) / 2 and p M = p 1 + p 2 are the position and momentum of the formed meson. The degeneracy factor g M = 1 / 12 is from the color and spin degrees of freedom. The Wigner function W is the Weyl transform of the charmonium wave function,
W ( r , q ) = d 3 y e i q · y ψ r + y 2 ψ * r y 2 ,
where ψ ( r ) is a charmonium wave function calculated from the Schrodinger equation with the in-medium heavy quark potential. When the wave function is approximated to be a Gaussian function, the corresponding Wigner function becomes W ( x r , q r ) = 8 exp [ x r 2 σ 2 σ 2 q r 2 ] [27]. x r and q r are the relative coordinates and momentum between c and c ¯ in the center of mass frame of charmonium. The information of the charmonium wave function is encoded in the Gaussian width, which is connected with the root-mean-square radius of charmonium, σ 2 = 4 3 ( m 1 + m 2 ) 2 m 1 2 + m 2 2 r 2 M [27]. The production N M A A of charmonium in AA collisions is proportional to the square of the charm pair number N c c ¯ A A [41],
N M A A = d x M d p M ( 2 π ) 3 P c c ¯ ψ ( x M , p M ) events ( N c c ¯ A A ) 2 ,
N c c ¯ A A = d x T T A ( x T b 2 ) T B ( x T + b 2 ) R S d σ p p c c ¯ d y Δ y c c ¯ ,
where the factor R S represents the nuclear shadowing effect. It changes the initial distribution of charm pairs in AA collisions, especially at the LHC energies.
From Equations (4)–(6), charmonium final production depends on the dynamical evolutions of charm quarks in the medium. Strongly coupled with the medium, charm quark density is reduced by the expansion of the medium. This indicates that, in the hotter medium, it takes a longer time to cool down and restore the heavy quark potential. As a result, the charm quarks diffuse into a larger volume before they combine to form a charmonium state. This diffusion effect will suppress the coalescence probability of one c and one c ¯ .
The thermal production of J / ψ at s N N = 2.76 TeV and 5.02 TeV Pb-Pb collisions are studied based on the coalescence model in Refs. [11,12,14,15,27,32]. At 5.02 TeV, J / ψ R A A is enhanced due to more significant number of charm pairs in the medium compared with the situation at 2.76 TeV. However, the coalescence probability per one c and c ¯ is reduced when the time of the medium expansion is longer. This effect is characterized by the ratio N J / ψ / ( N c c ¯ ) 2 of J / ψ production over the square of the charm pair number, as shown in Figure 2. The ratio N J / ψ / ( N c c ¯ ) 2 is calculated as a function of N p at the collision energies of s N N = ( 2.76 , 5.02 , 39 ) TeV. Different time scales of charm momentum thermalization are considered in the figure. In upper and lower panels of the figure, charm quarks are assumed to reach momentum thermalization at τ τ t . Charm quarks perform free streaming before the time scale τ t and diffuse according to the diffusion equation in τ > τ t . In more central collisions, the initial energy density of the medium becomes larger. It takes a longer time for the medium to decrease to T c , where the coalescence of c and c ¯ happens. The spatial density of charm quarks is diluted to reduce the coalescence probability of one c and c ¯ . Similarly, at higher collision energies, with the production of higher energy density, the diffusion of charm quarks is more vital, and charm spatial density is more reduced.
The charm diffusion in the expanding medium tends to suppress the nuclear modification factor R A A . In Figure 3, when employing the same hydrodynamic evolution at both 2.76 TeV and 5.02 TeV as a test, R A A ( 5.02 T e V ) is enhanced because of a larger production cross-section of charm pairs d σ p p c c ¯ / d y at 5.02 TeV compared with the situation at 2.76 TeV. Their ratio is plotted with blue dotted line in Figure 3. When taking the realistic evolutions for the medium produced at each collision energy, the coalescence probability per one c and c ¯ is reduced at 5.02 TeV due to a stronger diffusion of charm quarks before the coalescence process. The difference between the blue dotted line (with the same hydrodynamic evolution) and the black solid line (with two realistic hydrodynamic evolutions) are due to the different diffusions of charm quarks in two collision energies.

3.2. Elliptic Flow v 2

In non-central nuclear collisions, the collision area of two nuclei is elliptic. The gradient of the initial energy density becomes anisotropic in the transverse plane. With different accelerations in the x- and y-directions, the azimuthal distribution d N / d ϕ of the final hadrons is anisotropic as well, which can be expanded as
d N / d ϕ 1 + 2 n = 1 v n cos n ( ϕ Φ n ) ,
where v n and Φ n are the magnitude and the phase of the nth-order harmonic flows in the azimuthal angle distribution of final particles. As charm quarks are strongly coupled with the medium, charm quarks diffuse in the expanding medium to carry elliptic flows. At temperatures close to T c where heavy quark potential is partially restored, charmonium is regenerated via the coalescence of charm quarks and inherits the collective flows of charm quarks.
In AA collisions, final charmonium comes from two mechanisms: the primordial production at the beginning of nuclear collisions and the regeneration from the coalescence of charm quarks in the QGP expansion. Charmonium bound states suffer dissociation in the hot medium. On the other hand, the regeneration process increases the final production of charmonium. Both contributions can be included in the transport equation. The charmonium distribution in phase space f Ψ evolves with time. It can be simplified as [11,12,50,51],
cosh ( y η ) τ + sinh ( y η ) τ η + v T · T f Ψ = α Ψ f Ψ + β Ψ
where τ and η are the proper time and the pseudo-rapidity, respectively. Here, y is the rapidity in momentum space. v T = p T / E T = p T / m Ψ 2 + p T 2 is the transverse velocity of charmonium moving in the medium. The terms on the left-hand side represent the charmonium distribution f Ψ changing with time τ and coordinates ( η , x T ). On the right-hand side, α and β represent the decay rate and the regeneration rate of charmonium. The decay rate α Ψ is proportional to the density of thermal gluons and the inelastic cross-section. The regeneration rate β Ψ in the reaction c + c ¯ Ψ + g is proportional to the densities of charm and anti-charm quarks and their coalescence probability. The formulas of α Ψ and β Ψ are given in Equations (10) and (11),
α Ψ ( p , x , τ | b ) = 1 2 E T d 3 p g ( 2 π ) 3 2 E g W g ψ c c ¯ ( s ) f g ( p g , x , τ ) Θ ( T ( x , τ | b ) T c ) ,
β Ψ ( p , x , τ | b ) = 1 2 E T d 3 p g ( 2 π ) 3 2 E g d 3 p c ( 2 π ) 3 2 E c d 3 p c ¯ ( 2 π ) 3 2 E c ¯ × W c c ¯ g ψ ( s ) f c ( p c , x , τ | b ) f c ¯ ( p c ¯ , x , τ | b ) × ( 2 π ) 4 δ ( 4 ) ( p + p g p c p c ¯ ) Θ ( T ( x , τ | b ) T c ) .
In the above formulas, f g is the gluon density, and W g ψ c c ¯ is the dissociation rate calculated from pQCD [52,53]. The step function Θ ( T T c ) ensures that the parton dissociation only exists in the QGP. In β Ψ , W c c ¯ g ψ is connected with W g ψ c c ¯ via the detailed balance. δ functions indicate the energy-momentum conservation in the coalescence process. f c and f c ¯ are the density of c and c ¯ quarks. Their distribution can be obtained from the Langevin model. As D meson elliptic flows are large and close to the flows of light hadrons [34,54], in the thermalization limit, charm spatial density satisfies the equation μ ( ρ c u μ ) = 0 , where u μ is the four-velocity of the fluid. Charm momentum distribution is taken as a normalized Fermi distribution as an approximation.
The initial distribution of primordially produced charmonium in nuclear collisions is obtained from the production in proton-proton collisions scaled with the number of binary collisions N c o l l . In proton-proton collisions, the p T -spectrum of J / ψ can be parametrized with the form
f p p J / ψ d 2 σ pp J / ψ d y 2 π p T d p T = d σ pp J / ψ d y × ( n 1 ) π ( n 2 ) p T 2 p p [ 1 + p T 2 ( n 2 ) p T 2 p p ] n ,
where the rapidity-differential cross-section d σ p p J / ψ / d y is determined by experiments [55,56]. The parameters n and p T 2 characterize the shape of the J / ψ p T -distribution in pp collisions [57]. In AA collisions, the cold nuclear matter (CNM) effect also changes the production of charmonium. The nuclear shadowing effect modifies the parton densities in the nucleus, which affects the charmonium production. The partons from nucleons scatter with other nucleons before fusing into a charmonium and charm pair. This process increases the p T of produced charmonium compared with the situation in pp collisions. It is also called the Cronin effect, which can be included by modifying f p p J / ψ with a Gaussian smearing method. All of the cold nuclear matter effects happen before the hot medium effects. Therefore, those CNM effects are included in the initial distribution of Equation (9) [58],
f Ψ ( x , p T , y , τ 0 ) = ( 2 π ) 3 δ ( z ) T A ( x T b / 2 ) T B ( x T + b / 2 ) R S A R S B f p p J / ψ ¯ ,
where the shadowing factor R S A ( B ) is calculated with the EPS09 NLO calculation [59]. f p p J / ψ ¯ is the momentum distribution of primordially produced J / ψ with the modification of Cronin effect.
Solving the transport model, the J / ψ nuclear modification factor including both primordial production and regeneration is plotted in the left panel of Figure 4. Dotted, dashed and solid lines represent the primordial production, regeneration, and inclusive production, respectively. In semi-central and central collisions, the regeneration dominates the final production of J / ψ . The regenerated charmonium is mainly located in the low p T region because of the charm quark energy loss. The primordial production dominates the production at high p T . Therefore, in the right panel of Figure 4, v 2 of the final inclusive J / ψ is significant in low p T regions because of the collective flows of charm quarks. At high p T , v 2 is small because J / ψ as a color-singlet state is not strongly coupled with the medium, as with single charm quarks. Primordially produced charmonium performs free streaming in the medium, and the small non-zero v 2 of J / ψ is induced by the path-length difference of trajectories when they move along different directions in the transverse plane. At the high p T region in the right panel of Figure 4, theoretical calculations underestimate the experimental data. This may be due to the non-thermal distribution of charm quarks in the medium [60].

3.3. Triangular Flow v 3

In the event-by-event heavy-ion collisions, the nucleon positions fluctuate and deviate from the smooth distribution of the nucleus [61]. The energy density produced by the collisions also fluctuate event-by-event. The realistic fluctuating profiles of the medium energy density is simulated with the Monte Carlo Glauber model (MC-Glauber) [62] and Monte Carlo fKLN model (MC-KLN) [63]. In the dynamical expansion, where the anisotropy in the initial distribution of medium energy density is transformed into the anisotropy in the momentum distribution of hadrons, there is a non-zero triangular flow in the azimuthal distribution of hadrons. In experiments, the v 3 of light hadrons have been observed and explained well with the hydrodynamic model plus the fluctuating initial conditions. The value of v 3 is proportional to the triangularity in the initial fluctuating distribution. The regenerated charmonium from the coalescence of charm quarks also carry collective flows. The magnitude of the charmonium v 3 is determined by the initial fluctuations and the degree of charm kinetic thermalization.
Firstly, one can generate a large set of fluctuating initial conditions with the Monte Carlo method. Then, we can take these fluctuating medium evolutions into the transport model to calculate the v 3 of J / ψ in each event. This calculation needs a huge computation time. Instead, one can take single-shot hydrodynamics where triangular distributions of the initial energy density from fluctuations are encoded in the initial conditions. To introduce the triangularity in the initial energy density of the medium, a deformation factor is included in the formula of the energy denstiy [64],
ϵ ( x , τ 0 | b ) ϵ ( x ˜ , τ 0 | b ) ,
where the deformed coordinate is x ˜ = ( x T 1 + ε 3 cos [ 3 ( ϕ s Ψ 3 ) ] , ϕ s , η ) . x T = x 2 + y 2 and ϕ s = arctan ( | y | / | x | ) are the transverse radius and azimuth angle. Ψ 3 = 0 is regarded as the reaction plane angle in the smooth hydrodynamic evolution. The value of the parameter ϵ 3 characterizing the magnitude of the fluctuation is determined to be ϵ 3 = 0.08 , 0.1 , 0.2 in the centrality bins 0–10%, 10–30% and 30–50%, respectively, in 5.02 TeV Pb-Pb collisions according to the simulation from the MC-Glauber or MC-KLN model [64,65]. The spatial distributions of the energy density with and without the fluctuations are plotted in Figure 5.
In the left panel of Figure 5, the shape of the initial energy density is elliptic, which comes from the geometry size of the collision area. In the right panel of Figure 5, the triangular shape in the profile of the energy density is encoded, which characterizes the magnitude of the fluctuations. The corresponding elliptic and triangular flows ( v 2 , v 3 ) of J / ψ are calculated based on the transport equation by taking the corresponding hydrodynamic evolutions; please see Figure 6.
In Figure 6, elliptic flows depend sensitively on the collision centrality. In a central collision, the pressure gradient of the medium energy density is nearly isotropic in the transverse plane. Elliptic flows of the medium and heavy flavor particles are small in the centrality range of 0–10%. In a centrality range of 10–30%, the elliptic flow becomes larger due to the elliptic profile of the collision area, where medium acceleration is different along the x- and y-directions in the transverse plane. In the high p T region, theoretical calculations underestimate the experimental data. This may be due to the non-thermal distribution of charm quarks [60]. For triangular flows, they mainly come from the fluctuation of the medium energy density, which shows relatively weak dependence on the impact parameter. The charmoniums v 3 in three centralities are close to each other without significant differences. The observation of J / ψ   v 3 indicates clearly that charm quarks are strongly coupled with the medium to be affected by the fluctuations.

3.4. Directed Flow v 1

In the previous discussions, we have discussed the anisotropy of the medium initial energy density in the transverse plane. In this section, we discuss the rapidity-odd distribution of the energy density induced by the rotation of the medium in semi-central nuclear collisions. In the collision area, the nucleus moving with the positive (negative) rapidity tends to tilt the produced bulk medium to positive (negative) x with respect to the beam axis for positive (negative) z. The magnitude of the tilt in the bulk medium is connected with the directed flows of light-charged hadrons, which have been observed in experiments [67,68]. In the Au-Au collisions at RHIC energy, the rapidity distribution of the medium energy density is plotted in the centrality of 0–80%; please see the left panel of Figure 7. In the forward (backward) rapidity, the medium is tilted towards the positive (negative) x-direction. The initial entropy density of the medium is parametrized with the formula
s ( τ 0 , x T , η ) = s 0 × exp [ θ ( | η | η 0 ) ( | η | η 0 ) 2 2 σ 2 ] × [ c h a r d N c o l l + ( 1 c h a r d ) ( N p a r t + ζ + ( η ) + N p a r t ζ ( η ) ) ] ,
where s 0 and the parameters in the exponential factor characterize the initial entropy density and the rapidity-odd distribution of the medium, respectively [69]. The initial entropy density depends on both soft and hard processes, where the fractions are ( 1 c hard ) and c hard = 0.05 , respectively. N p a r t + and N p a r t are the number of participants in the forward and backward rapidities. The feature of the rapidity-odd distribution is introduced via the function ζ + , ( η ) . The hydrodynamic evolutions of the medium in both transverse and longitudinal directions are given by the MUSIC calculation with the input of the rapidity-odd initial entropy density.
The directed flows of D mesons have been observed [70]. Compared with light hadrons, charmonium is produced in the initial parton hard scatterings where the tilted shape of the entropy distribution is apparent. Different from charm quarks, primordially produced charmonium with zero net color charge is weakly coupled with the medium. Therefore, charmonium suffers biased dissociation along the positive and negative x-directions in the tilted medium, and this information is not contaminated by the following QGP evolutions where the tilted shape of the medium is less evident. In the right panel of Figure 7, directed flows of J / ψ and ψ ( 2 S ) as a function of pseudorapidity η are calculated in different p T bins in s N N = 200 GeV Au-Au collisions [69]. As one can see, in the forward rapidity, charmonium moving along positive x-direction suffer stronger medium dissociation due to the longer path-length in the medium. This results in anisotropy in the azimuthal distribution of the final charmonium, where charmonium v 1 is negative in forward rapidity and negative in backward rapidity. For an excited state ψ ( 2 S ) , the binding energy is smaller than the value of the ground state. They are more easily dissociated and sensitive to the profile of the medium entropy density. In Figure 7, charmonium directed flows are bigger than the values of the light hadrons. In the above calculations regarding v 1 , the regeneration process has been neglected for simplicity, which becomes less important, especially in the high p T region where the directed flow is strong, as shown by the lines with p T > 4 GeV/c in Figure 7. Confirming and quantifying the charmonium v 1 help to extract the rotation and the rapidity-odd distribution of the hot medium.

4. Summary

In the present work, we review recent studies about charmonium anisotropic flows in relativistic heavy-ion collisions. The pressure gradients of the medium energy density become different along different directions, where the acceleration of the medium expansion is also anisotropic. With strong coupling between heavy quarks and the bulk medium, charm quarks approach kinetic thermalization to carry collective flows from the expanding medium. The anisotropic flows of charm quarks will be inherited by the regenerated charmonium via the coalescence process. The elliptic flow of charmonium originates from the shape of the QGP produced in the collision area of two nuclei in semi-central collisions. It depends on the impact parameter and becomes almost zero in most central collisions. The triangular flow of charmonium comes from the fluctuations of the medium initial energy density. As the fluctuations of nucleon positions in the nucleus show weak dependence on the impact parameter, charmonium v 3 are similar to each other without evident differences in central and semi-central collisions. In the longitudinal direction, the rapidity-odd distribution of the energy density is introduced due to the rotation of the hot medium produced in non-central collisions. This results in biased charmonium dissociation along the positive and negative x-directions due to the different path-length in QGP. The directed flows of J / ψ and ψ ( 2 S ) are larger than the case of light hadrons. As primordially produced charmonium is mostly dissociated by the medium in the initial stage, where the tilted shape of the medium is evident, charmonium directed flow is proposed as a sensitive probe of the profile of the medium initial energy densitiy in relativistic heavy-ion collisions.

Funding

This work is supported by the National Natural Science Foundation of China (NSFC) under Grant No. 12175165.

Data Availability Statement

Not applicable.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Bazavov, A.; Bhattacharya, T.; Cheng, M.; DeTar, C.; Ding, H.-T.; Gottlieb, S.; Gupta, R.; Hegde, P.; Heller, U.M.; Karsch, F.; et al. The chiral and deconfinement aspects of the QCD transition. Phys. Rev. D 2012, 85, 054503. [Google Scholar] [CrossRef] [Green Version]
  2. Satz, H. Colour deconfinement and quarkonium binding. J. Phys. G 2006, 32, R25. [Google Scholar] [CrossRef] [Green Version]
  3. Aoki, Y.; Endrődi, G.; Fodor, Z.; Katz, S.D.; Szabó, K.K. The order of the quantum chromodynamics transition predicted by the standard model of particle physics. Nature 2006, 443, 675–678. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  4. Heinz, U.; Snellings, R. Collective flow and viscosity in relativistic heavy-ion collisions. Ann. Rev. Nucl. Part. Sci. 2013, 63, 123–151. [Google Scholar] [CrossRef] [Green Version]
  5. Song, H.; Bass, S.A.; Heinz, U.; Hirano, T.; Shen, C. 200 A GeV Au+Au collisions serve a nearly perfect quark-gluon liquid. Phys. Rev. Lett. 2011, 106, 192301, Erratum in Phys. Rev. Lett. 2012, 109, 139904. [Google Scholar] [CrossRef] [Green Version]
  6. Gale, C.; Jeon, S.; Schenke, B.; Tribedy, P.; Venugopalan, R. Event-by-event anisotropic flow in heavy-ion collisions from combined Yang-Mills and viscous fluid dynamics. Phys. Rev. Lett. 2013, 110, 012302. [Google Scholar] [CrossRef] [Green Version]
  7. Adams, J.; Aggarwal, M.M.; Ahammed, Z.; Amonett, J.; Anderson, B.D.; Arkhipkin, D.; Averichev, G.S.; Badyal, S.K.; Bai, Y.; Balewski, J.; et al. Azimuthal anisotropy in Au+Au collisions at s(NN)**(1/2) = 200-GeV. Phys. Rev. C 2005, 72, 014904. [Google Scholar] [CrossRef] [Green Version]
  8. Aamodt, K.; Abelev, B.; Abrahantes Quintana, A.; Adamová, D.; Adare, A.M.; Aggarwal, M.M.; Aglieri Rinella, G.; Agocs, A.G.; Agostinelli, A.; Aguilar Salazar, S.; et al. Higher harmonic anisotropic flow measurements of charged particles in Pb-Pb collisions at s N N = 2.76 TeV. Phys. Rev. Lett. 2011, 107, 032301. [Google Scholar] [CrossRef] [Green Version]
  9. Matsui, T.; Satz, H. J/ψ Suppression by Quark-Gluon Plasma Formation. Phys. Lett. B 1986, 178, 416–422. [Google Scholar] [CrossRef]
  10. Braun-Munzinger, P.; Stachel, J. (Non)thermal aspects of charmonium production and a new look at J/ψ suppression. Phys. Lett. B 2000, 490, 196–202. [Google Scholar] [CrossRef]
  11. Zhu, X.L.; Zhuang, P.F.; Xu, N. J/ψ transport in QGP and p(t) distribution at SPS and RHIC. Phys. Lett. B 2005, 607, 107–114. [Google Scholar] [CrossRef] [Green Version]
  12. Yan, L.; Zhuang, P.; Xu, N. Competition between J/ψ suppression and regeneration in quark-gluon plasma. Phys. Rev. Lett. 2006, 97, 232301. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  13. Liu, Y.; Chen, B.; Xu, N.; Zhuang, P. Υ Production as a Probe for Early State Dynamics in High Energy Nuclear Collisions at RHIC. Phys. Lett. B 2011, 697, 32–36. [Google Scholar] [CrossRef] [Green Version]
  14. Grandchamp, L.; Rapp, R.; Brown, G.E. In medium effects on charmonium production in heavy ion collisions. Phys. Rev. Lett. 2004, 92, 212301. [Google Scholar] [CrossRef] [Green Version]
  15. Zhao, X.; Rapp, R. Medium Modifications and Production of Charmonia at LHC. Nucl. Phys. A 2011, 859, 114–125. [Google Scholar] [CrossRef] [Green Version]
  16. Blaizot, J.P.; Boni, D.D.; Faccioli, P.; Garberoglio, G. Heavy quark bound states in a quark–gluon plasma: Dissociation and recombination. Nucl. Phys. A 2016, 946, 49–88. [Google Scholar] [CrossRef] [Green Version]
  17. Krouppa, B.; Ryblewski, R.; Strickland, M. Bottomonia suppression in 2.76 TeV Pb-Pb collisions. Phys. Rev. C 2015, 92, 061901. [Google Scholar] [CrossRef] [Green Version]
  18. Yao, X.; Ke, W.; Xu, Y.; Bass, S.A.; Müller, B. Coupled Boltzmann Transport Equations of Heavy Quarks and Quarkonia in Quark-Gluon Plasma. J. High Energy Phys. 2021, 2021, 46. [Google Scholar] [CrossRef]
  19. Yao, X.; Mehen, T. Quarkonium Semiclassical Transport in Quark-Gluon Plasma: Factorization and Quantum Correction. J. High Energy Phys. 2021, 2021, 62. [Google Scholar] [CrossRef]
  20. Wen, L.; Du, X.; Shi, S.; Chen, B. Probe the color screening in proton-nucleus collisions with complex potentials. Chin. Phys. C 2022, 46, 114102. [Google Scholar] [CrossRef]
  21. Shi, W.; Zha, W.; Chen, B. Charmonium Coherent Photoproduction and Hadroproduction with Effects of Quark Gluon Plasma. Phys. Lett. B 2018, 777, 399–405. [Google Scholar] [CrossRef]
  22. Abreu, M.C.; Alessandro, B.; Alexa, C.; Arnaldi, R.; Atayan, M.; Baglin, C.; Baldit, A.; Bedjidian, M.; Bellaiche, F.; Beolè, S.; et al. Transverse momentum distributions of J/ψ, ψ′, Drell–Yan and continuum dimuons produced in Pb–Pb interactions at the SPS. Phys. Lett. B 2001, 499, 85–96. [Google Scholar] [CrossRef] [Green Version]
  23. Adare, A.; Afanasiev, S.; Aidala, C.; Ajitanand, N.N.; Akiba, Y.; Al-Bataineh, H.; Alexander, J.; Al-Jamel, A.; Aoki, K.; Aphecetche, L.; et al. J/ψ Production versus Centrality, Transverse Momentum, and Rapidity in Au+Au Collisions at s N N = 200 GeV. Phys. Rev. Lett. 2007, 98, 232301. [Google Scholar] [CrossRef] [Green Version]
  24. Adam, J.; Adamová, D.; Aggarwal, M.M.; Aglieri Rinella, G.; Agnello, M.; Agrawal, N.; Ahammed, Z.; Ahn, S.U.; Aimo, I.; Aiola, S.; et al. Differential studies of inclusive J/ψ and ψ(2S) production at forward rapidity in Pb-Pb collisions at s N N = 2.76 TeV. J. High Energy Phys. 2016, 2016, 179. [Google Scholar] [CrossRef] [Green Version]
  25. Adam, J.; Adamová, D.; Aggarwal, M.M.; Aglieri Rinella, G.; Agnello, M.; Agrawal, N.; Ahammed, Z.; Ahmad, S.; Ahn, S.U.; Aiola, S.; et al. J/ψ suppression at forward rapidity in Pb-Pb collisions at s N N = 5.02 TeV. Phys. Lett. B 2017, 766, 212–224. [Google Scholar] [CrossRef]
  26. Thews, R.L.; Schroedter, M.; Rafelski, J. Enhanced J/ψ production in deconfined quark matter. Phys. Rev. C 2001, 63, 054905. [Google Scholar] [CrossRef] [Green Version]
  27. Greco, V.; Ko, C.M.; Rapp, R. Quark coalescence for charmed mesons in ultrarelativistic heavy ion collisions. Phys. Lett. B 2004, 595, 202–208. [Google Scholar] [CrossRef] [Green Version]
  28. Andronic, A.; Braun-Munzinger, P.; Redlich, K.; Stachel, J. Statistical hadronization of charm in heavy ion collisions at SPS, RHIC and LHC. Phys. Lett. B 2003, 571, 36–44. [Google Scholar] [CrossRef] [Green Version]
  29. Fries, R.J.; Greco, V.; Sorensen, P. Coalescence Models For Hadron Formation From Quark Gluon Plasma. Ann. Rev. Nucl. Part. Sci. 2008, 58, 177–205. [Google Scholar] [CrossRef] [Green Version]
  30. Du, X.; Rapp, R. Sequential Regeneration of Charmonia in Heavy-Ion Collisions. Nucl. Phys. A 2015, 943, 147–158. [Google Scholar] [CrossRef]
  31. Chen, B.; Zhao, J. Bottomonium Continuous Production from Unequilibrium Bottom Quarks in Ultrarelativistic Heavy Ion Collisions. Phys. Lett. B 2017, 772, 819–824. [Google Scholar] [CrossRef]
  32. Zhao, J.; Chen, B. Strong diffusion effect of charm quarks on J /ψ production in Pb–Pb collisions at the LHC. Phys. Lett. B 2018, 776, 17–21. [Google Scholar] [CrossRef]
  33. Chatrchyan, S.; Khachatryan, V.; Sirunyan, A.M.; Tumasyan, A.; Adam, W.; Bergauer, T.; Dragicevic, M.; Erö, J.; Fabjan, C.; Friedl, M.; et al. Suppression of non-prompt J/ψ, prompt J/ψ, and Y(1S) in PbPb collisions at s N N = 2.76 TeV. J. High Energy Phys. 2012, 2012, 63. [Google Scholar]
  34. Abelev, B.; Adam, J.; Adamová, D.; Adare, A.M.; Aggarwal, M.M.; Rinella, G.A.; Agnello, M.; Agocs, A.G.; Agostinelli, A.; Ahammed, Z.; et al. D meson elliptic flow in non-central Pb-Pb collisions at s N N = 2.76TeV. Phys. Rev. Lett. 2013, 111, 102301. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  35. Cao, S.; Qin, G.Y.; Bass, S.A. Heavy-quark dynamics and hadronization in ultrarelativistic heavy-ion collisions: Collisional versus radiative energy loss. Phys. Rev. C 2013, 88, 044907. [Google Scholar] [CrossRef]
  36. Cao, S.; Coci, G.; Das, S.K.; Ke, W.; Liu, S.Y.F.; Plumari, S.; Song, T.; Xu, Y.; Aichelin, J.; Bass, S.; et al. Toward the determination of heavy-quark transport coefficients in quark-gluon plasma. Phys. Rev. C 2019, 99, 054907. [Google Scholar] [CrossRef] [Green Version]
  37. He, M.; Fries, R.J.; Rapp, R. Ds-Meson as Quantitative Probe of Diffusion and Hadronization in Nuclear Collisions. Phys. Rev. Lett. 2013, 110, 112301. [Google Scholar] [CrossRef] [Green Version]
  38. He, M.; Rapp, R. Hadronization and Charm-Hadron Ratios in Heavy-Ion Collisions. Phys. Rev. Lett. 2020, 124, 042301. [Google Scholar] [CrossRef] [Green Version]
  39. Chen, B. Elliptic flow as a probe for ψ(2S) production mechanism in relativistic heavy ion collisions. Phys. Rev. C 2017, 95, 034908. [Google Scholar] [CrossRef] [Green Version]
  40. Liu, Y.P.; Qu, Z.; Xu, N.; Zhuang, P.F. J/ψ Transverse Momentum Distribution in High Energy Nuclear Collisions at RHIC. Phys. Lett. B 2009, 678, 72–76. [Google Scholar] [CrossRef] [Green Version]
  41. Chen, B.; Jiang, L.; Liu, X.H.; Liu, Y.; Zhao, J. X(3872) Production in Relativistic Heavy-Ion Collisions. Phys. Rev. C 2022, 105, 054901. [Google Scholar] [CrossRef]
  42. Cao, S.; Qin, G.Y.; Bass, S.A. Energy loss, hadronization and hadronic interactions of heavy flavors in relativistic heavy-ion collisions. Phys. Rev. C 2015, 92, 024907. [Google Scholar] [CrossRef] [Green Version]
  43. Rapp, R.; Gossiaux, P.B.; Andronic, A.; Averbeck, R.; Masciocchi, S.; Beraudo, A.; Bratkovskaya, E.; Braun-Munzinger, P.; Cao, S.; Dainese, A.; et al. Extraction of Heavy-Flavor Transport Coefficients in QCD Matter. Nucl. Phys. A 2018, 979, 21–86. [Google Scholar] [CrossRef] [Green Version]
  44. Zhao, J.; Zhou, K.; Chen, S.; Zhuang, P. Heavy flavors under extreme conditions in high energy nuclear collisions. Prog. Part. Nucl. Phys. 2020, 114, 103801. [Google Scholar] [CrossRef]
  45. Guo, X.f.; Wang, X.N. Multiple scattering, parton energy loss and modified fragmentation functions in deeply inelastic e A scattering. Phys. Rev. Lett. 2000, 85, 3591–3594. [Google Scholar] [CrossRef] [Green Version]
  46. Zhang, B.W.; Wang, E.; Wang, X.N. Heavy quark energy loss in nuclear medium. Phys. Rev. Lett. 2004, 93, 072301. [Google Scholar] [CrossRef] [Green Version]
  47. Cacciari, M.; Frixione, S.; Nason, P. The p(T) spectrum in heavy flavor photoproduction. J. High Energy Phys. 2001, 2001, 6. [Google Scholar] [CrossRef] [Green Version]
  48. Huovinen, P.; Petreczky, P. QCD Equation of State and Hadron Resonance Gas. Nucl. Phys. A 2010, 837, 26–53. [Google Scholar] [CrossRef] [Green Version]
  49. Schenke, B.; Jeon, S.; Gale, C. (3+1)D hydrodynamic simulation of relativistic heavy-ion collisions. Phys. Rev. C 2010, 82, 014903. [Google Scholar] [CrossRef] [Green Version]
  50. Chen, B. Detailed rapidity dependence of J/ψ production at energies available at the Large Hadron Collider. Phys. Rev. C 2016, 93, 054905. [Google Scholar] [CrossRef] [Green Version]
  51. Chen, B.; Liu, Y.; Zhou, K.; Zhuang, P. ψ′ Production and B Decay in Heavy Ion Collisions at LHC. Phys. Lett. B 2013, 726, 725–728. [Google Scholar] [CrossRef]
  52. Peskin, M.E. Short Distance Analysis for Heavy Quark Systems. 1. Diagrammatics. Nucl. Phys. B 1979, 156, 365–390. [Google Scholar] [CrossRef]
  53. Bhanot, G.; Peskin, M.E. Short Distance Analysis for Heavy Quark Systems. 2. Applications. Nucl. Phys. B 1979, 156, 391–416. [Google Scholar] [CrossRef]
  54. Abelev, B.; Adam, J.; Adamová, D.; Adare, A.M.; Aggarwal, M.M.; Rinella, G.A.; Agocs, A.G.; Agostinelli, A.; Salazar, S.A.; Ahammed, Z.; et al. [ALICE], Anisotropic flow of charged hadrons, pions and (anti-)protons measured at high transverse momentum in Pb-Pb collisions at s N N = 2.76 TeV. Phys. Lett. B 2013, 719, 18–28. [Google Scholar] [CrossRef]
  55. Abelev, B.; Adam, J.; Adamová, D.; Adare, A.M.; Aggarwal, M.M.; Rinella, G.A.; Agocs, A.G.; Agostinelli, A.; Salazar, S.A.; Ahammed, Z.; et al. [ALICE], Inclusive J/ψ production in pp collisions at s = 2.76 TeV. Phys. Lett. B 2012, 718, 295–306, Erratum in Phys. Lett. B 2015, 748, 472–473. [Google Scholar] [CrossRef]
  56. Abelev, B.; Adam, J.; Adamova, D.; Adare, A.M.; Aggarwal, M.; Rinella, G.A.; Agocs, A.G.; Agostinelli, A.; Salazar, S.A.; Ahammed, Z.; et al. Measurement of prompt J/ψ and beauty hadron production cross sections at mid-rapidity in pp collisions at s=7 TeV. J. High Energy Phys. 2012, 2012, 65. [Google Scholar]
  57. Chen, B. Thermal production of charmonia in Pb-Pb collisions at s N N = 5.02 TeV. Chin. Phys. C 2019, 43, 124101. [Google Scholar] [CrossRef] [Green Version]
  58. Chen, B.; Guo, T.; Liu, Y.; Zhuang, P. Cold and Hot Nuclear Matter Effects on Charmonium Production in p+Pb Collisions at LHC Energy. Phys. Lett. B 2017, 765, 323–327. [Google Scholar] [CrossRef]
  59. Eskola, K.J.; Paukkunen, H.; Salgado, C.A. EPS09—A New Generation of NLO and LO Nuclear Parton Distribution Functions. J. High Energy Phys. 2009, 2009, 65. [Google Scholar] [CrossRef] [Green Version]
  60. He, M.; Wu, B.; Rapp, R. Collectivity of J/ψ Mesons in Heavy-Ion Collisions. Phys. Rev. Lett. 2022, 128, 162301. [Google Scholar] [CrossRef]
  61. Alver, B.; Roland, G. Collision geometry fluctuations and triangular flow in heavy-ion collisions. Phys. Rev. C 2010, 81, 054905, Erratum in Phys. Rev. C 2010, 82, 039903.. [Google Scholar] [CrossRef]
  62. Miller, M.L.; Reygers, K.; Sanders, S.J.; Steinberg, P. Glauber modeling in high energy nuclear collisions. Ann. Rev. Nucl. Part. Sci. 2007, 57, 205–243. [Google Scholar] [CrossRef]
  63. Kharzeev, D.; Levin, E. Manifestations of high density QCD in the first RHIC data. Phys. Lett. B 2001, 523, 79–87. [Google Scholar] [CrossRef] [Green Version]
  64. Alver, B.H.; Gombeaud, C.; Luzum, M.; Ollitrault, J.Y. Triangular flow in hydrodynamics and transport theory. Phys. Rev. C 2010, 82, 034913. [Google Scholar] [CrossRef]
  65. Qiu, Z.; Heinz, U.W. Event-by-event shape and flow fluctuations of relativistic heavy-ion collision fireballs. Phys. Rev. C 2011, 84, 024911. [Google Scholar] [CrossRef] [Green Version]
  66. Zhao, J.; Chen, B.; Zhuang, P. Charmonium triangular flow in high energy nuclear collisions. Phys. Rev. C 2022, 105, 034902. [Google Scholar] [CrossRef]
  67. Abelev, B.I.; Aggarwal, M.M.; Ahammed, Z.; Anderson, B.D.; Arkhipkin, D.; Averichev, G.S.; Bai, Y.; Balewski, J.; Barannikova, O.; Barnby, L.S.; et al. System-size independence of directed flow at the Relativistic Heavy-Ion Collider. Phys. Rev. Lett. 2008, 101, 252301. [Google Scholar] [CrossRef]
  68. Adamczyk, L.; Adkins, J.K.; Agakishiev, G.; Aggarwal, M.M.; Ahammed, Z.; Alekseev, I.; Alford, J.; Anson, C.D.; Aparin, A.; Arkhipkin, D.; et al. Beam-Energy Dependence of the Directed Flow of Protons, Antiprotons, and Pions in Au+Au Collisions. Phys. Rev. Lett. 2014, 112, 162301. [Google Scholar] [CrossRef]
  69. Chen, B.; Hu, M.; Zhang, H.; Zhao, J. Probe the tilted Quark-Gluon Plasma with charmonium directed flow. Phys. Lett. B 2020, 802, 135271. [Google Scholar] [CrossRef]
  70. Adam, J.; Adamczyk, L.; Adams, J.R.; Adkins, J.K.; Agakishiev, G.; Aggarwal, M.M.; Ahammed, Z.; Alekseev, I.; Anderson, D.M.; Aoyama, R.; et al. First Observation of the Directed Flow of D0 and D 0 ¯ in Au+Au Collisions at s N N = 200 GeV. Phys. Rev. Lett. 2019, 123, 162301. [Google Scholar] [CrossRef]
Figure 1. The time evolution of the local temperature at the center of the medium in the central rapidity in s N N = 5.02 TeV Pb-Pb collisions. The collision centrality is chosen as 0–20%, 20–40%, 40–60% and 60–80%. This temperature evolution is calculated with a MUSIC calculation. This figure is cited from Ref. [41].
Figure 1. The time evolution of the local temperature at the center of the medium in the central rapidity in s N N = 5.02 TeV Pb-Pb collisions. The collision centrality is chosen as 0–20%, 20–40%, 40–60% and 60–80%. This temperature evolution is calculated with a MUSIC calculation. This figure is cited from Ref. [41].
Mathematics 10 04284 g001
Figure 2. The ratio of J / ψ production to the square of the total charm pair number (with N c = N D + N Λ c + including all charm flavor hadrons) as a function of the number of participants N p a r t in Pb-Pb collisions. The collision energies are chosen as s N N = 2.76 , 5.02 , 39 TeV, respectively. Charm quark momentum distribution is assumed to reach equilibrium at the beginning of medium expansion τ t = 0.6 fm/c in the upper panel and τ t = 4 fm/c in the lower panel. This figure is cited from Ref. [32].
Figure 2. The ratio of J / ψ production to the square of the total charm pair number (with N c = N D + N Λ c + including all charm flavor hadrons) as a function of the number of participants N p a r t in Pb-Pb collisions. The collision energies are chosen as s N N = 2.76 , 5.02 , 39 TeV, respectively. Charm quark momentum distribution is assumed to reach equilibrium at the beginning of medium expansion τ t = 0.6 fm/c in the upper panel and τ t = 4 fm/c in the lower panel. This figure is cited from Ref. [32].
Mathematics 10 04284 g002
Figure 3. The ratio of inclusive nuclear modification factors R A A at 2.76 TeV and 5.02 TeV Pb-Pb collisions, plotted with a black solid line. The blue dotted line is a situation where the same hydrodynamic evolution is used as a test in both collision energies. The reduction from dotted to black solid lines is due to the stronger diffusion of charm quarks in the hotter medium produced at 5.02 TeV. This figure is cited from Ref. [32].
Figure 3. The ratio of inclusive nuclear modification factors R A A at 2.76 TeV and 5.02 TeV Pb-Pb collisions, plotted with a black solid line. The blue dotted line is a situation where the same hydrodynamic evolution is used as a test in both collision energies. The reduction from dotted to black solid lines is due to the stronger diffusion of charm quarks in the hotter medium produced at 5.02 TeV. This figure is cited from Ref. [32].
Mathematics 10 04284 g003
Figure 4. (Left) J / ψ nuclear modification factor R A A as a function of number of participant N p at s N N = 5.02 TeV Pb-Pb collisions. Dotted, dashed and solid lines represent the primordial production, regeneration and inclusive production, respectively. The band corresponds to the uncertainty in the production cross-section of charm pair d σ p p c c ¯ / d y . (Right) Elliptic flows of inclusive J / ψ as a function of p T in centrality 0–20%. This figure is cited from Ref. [57].
Figure 4. (Left) J / ψ nuclear modification factor R A A as a function of number of participant N p at s N N = 5.02 TeV Pb-Pb collisions. Dotted, dashed and solid lines represent the primordial production, regeneration and inclusive production, respectively. The band corresponds to the uncertainty in the production cross-section of charm pair d σ p p c c ¯ / d y . (Right) Elliptic flows of inclusive J / ψ as a function of p T in centrality 0–20%. This figure is cited from Ref. [57].
Mathematics 10 04284 g004
Figure 5. The initial energy density in the transverse plane in s N N = 5.02 TeV Pb-Pb collisions. In the left panel, no fluctuations are considered. In the right panel, a deformation characterizing the triangularity in the energy density is included. This figure is cited from Ref. [66].
Figure 5. The initial energy density in the transverse plane in s N N = 5.02 TeV Pb-Pb collisions. In the left panel, no fluctuations are considered. In the right panel, a deformation characterizing the triangularity in the energy density is included. This figure is cited from Ref. [66].
Mathematics 10 04284 g005
Figure 6. (Left) Elliptic flows of J / ψ as a function of p T in s N N = 5.02 TeV Pb-Pb collisions. The dotted and dashed lines and shadow bands represent the elliptic flows of initially produced, regenerated and inclusive J / ψ , respectively. The collision centrality is 0–10%, 10–30%, and 30–50%, respectively. (Right) Triangular flows of charmonium as a function of p T . This figure is cited from Ref. [66].
Figure 6. (Left) Elliptic flows of J / ψ as a function of p T in s N N = 5.02 TeV Pb-Pb collisions. The dotted and dashed lines and shadow bands represent the elliptic flows of initially produced, regenerated and inclusive J / ψ , respectively. The collision centrality is 0–10%, 10–30%, and 30–50%, respectively. (Right) Triangular flows of charmonium as a function of p T . This figure is cited from Ref. [66].
Mathematics 10 04284 g006
Figure 7. (Left) Rapidity-odd distribution of initial entropy density in s N N = 200 GeV Au-Au collisions. The collision centrality is 0–80%, corresponding to the averaged impact parameter b = 8.3 fm. (Right) The directed flows of J / ψ and ψ ( 2 S ) as a function of rapidity η . Different p T bins ( p T > 0 and p T > 4 GeV/c) are calculated. This figure is cited from Ref. [69].
Figure 7. (Left) Rapidity-odd distribution of initial entropy density in s N N = 200 GeV Au-Au collisions. The collision centrality is 0–80%, corresponding to the averaged impact parameter b = 8.3 fm. (Right) The directed flows of J / ψ and ψ ( 2 S ) as a function of rapidity η . Different p T bins ( p T > 0 and p T > 4 GeV/c) are calculated. This figure is cited from Ref. [69].
Mathematics 10 04284 g007
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Li, C.; Chen, B. Anisotropic Flows of Charmonium in the Relativistic Heavy-Ion Collisions. Mathematics 2022, 10, 4284. https://doi.org/10.3390/math10224284

AMA Style

Li C, Chen B. Anisotropic Flows of Charmonium in the Relativistic Heavy-Ion Collisions. Mathematics. 2022; 10(22):4284. https://doi.org/10.3390/math10224284

Chicago/Turabian Style

Li, Chenyu, and Baoyi Chen. 2022. "Anisotropic Flows of Charmonium in the Relativistic Heavy-Ion Collisions" Mathematics 10, no. 22: 4284. https://doi.org/10.3390/math10224284

APA Style

Li, C., & Chen, B. (2022). Anisotropic Flows of Charmonium in the Relativistic Heavy-Ion Collisions. Mathematics, 10(22), 4284. https://doi.org/10.3390/math10224284

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop