Drag force on heavy quarks from holographic QCD

We study the drag force of a relativistic heavy quark using a holographic QCD model with conformal invariance broken by a background dilaton. The effects of the chemical potential and the confining scale on this quantity are analyzed. The drag force in this model is shown to be larger than that of supersymmetric Yang-Mills (SYM) plasma. In particular, the inclusion of the chemical potential and confining scale both enhance the drag force, in agreement with earlier findings. Moreover, we discuss how the chemical potential and confining scale influence the diffusion coefficient.


Introduction
The ultra-relativistic heavy-ion experimental programs at the Relativistic Heavy Ion Collider (RHIC) and the Large Hadron Collider (LHC) have created a new type of matter, namely the quark gluon plasma (QGP) [1][2][3]. One of the interesting features of QGP is jet quenching (parton energy loss): high energy partons propagate through the hot and dense medium, and lose their energy by interacting with the medium. Usually, the jet quenching phenomenon can be characterized by the jet quenching parameter, defined by the average transverse momentum square transferred from the traversing parton per unit mean free path [4]. Alternately, this phenomenon could be studied from the drag force: heavy quarks move through the QGP and experience a drag force, loosing energy as a consequence. In weakly coupling theories, this jet quenching has already been investigated in many studies [5][6][7][8][9][10][11]. However, a lot of experiments find strong evidence that QGP does not behave as a weakly coupled gas, but rather as a strongly coupled fluid. Consequently, computational techniques for strongly coupled, real-time QCD dynamics are required. Such techniques are now available via the gauge/gravity duality or AdS/CFT correspondence [12][13][14].
AdS/CFT, a conjectured duality between a string theory in the AdS space and a conformal field theory on the N = 4 dp/dt = −up + f dp/dt = 0 AdS space boundary, has provided many important insights for exploring various aspects of QGP [15]. In the framework of AdS/CFT, the drag force of a heavy quark moving in SYM plasma was proposed in [16,17]. Therein, this force could be related to the damping rate u (or referred to as the friction coefficient), defined by the Langevin equation, , subject to a driving force f. For , i.e., constant speed trajectory, f equals to the drag force. It was found that [16,17] where v, T, are the quark velocity, plasma temperature and 't Hooft coupling, respectively. Subsequently, this has been generalized to various cases, such as chemical potential [18][19][20], finite coupling [21], non-commutativity plasma [22] and AdS/QCD models [23][24][25]. Other studies in this direction can be found in [26][27][28][29][30][31][32][33][34].

T,µ T,µ
In this study, we investigate the drag force in a soft wall-like model, i.e., the SW model, which is motivated by the soft wall model of [35]. The SW model was proposed in [36] to investigate the heavy quark free energy and the QCD phase diagram. It was shown that such a model provides a good phenomenological description of the quark-antiquark interaction as well as some hadronic properties. Recently, the authors of Ref.
[37] investigated the imaginary potential in this model and found that the presence of the confining scale tends to reduce the quarkonia dissociation, and hence reverse the effect of the chemical potential. Moreover, there are other investigations of models of this type [38][39][40][41][42][43]. Motivated by this fact, we consider the drag force in this model to understand the possible physical implications of our results in this dual plasma. Another motivation for this study is that previous works either discuss the drag force in some strongly coupled thermal gauge theories (not holographic QCD) with chemical potential [18][19][20] or holographic QCD without chemical potential [23][24][25]. Here, we provide such analysis in holographic QCD with chemical potential.

T,µ T,µ
The structure of this paper is as follows. In the next section, we introduce the geometry of the SW model given in Ref. [36]. In Section 3 and Section 4, we analyze the behavior of the drag force and the diffusion coefficient for the SW model. Summary and discussions are given in the final section.

Setup
To begin with, we introduce the holographic models in terms of the action [44] where is the five-dimensional (5D) Newton constant. g denotes the determinant of the metric . represents the Ricci scalar. denotes the scalar that induces the deformation away from conformality.
refers to the potential, which contains the cosmological constant and some other terms. and are the field strength tensor and gauge kinetic function, respectively.
For Eq. (2), the equations of motion are where represents the Levi-Civita covariant derivative related to .

With
, the AdS black hole arises as a solution of Eq. (3) - (6). Subsequently, by considering a non-vanishing gauge-field component with the conditions and solving Eq. (3) -(6), one can obtain the AdS-RN spacetime, whose metric is given by with where R is the AdS radius. Q denotes the black hole charge, constrained in . z stands for the fifth coordinate with as the horizon and as the boundary. represents the chemical potential. The chemical potential implemented here is not the quark (or baryon) chemical potential of QCD, but a chemical potential corresponding to the R-charge of SYM. Nevertheless, it could serve as a simple way of introducing finite density effects into the system.
Moreover, the temperature of the black hole is Explicitly, for given in the boundary theory, of the bulk theory can be expressed as T,µ To emulate confinement in the boundary theory at vanishing temperature, one can introduce a warp factor to the background, similar to [38,45]. Then the metric of the SW model is given by [36] where c is the deformation parameter (related to the confining scale) with the dimension of energy. Here, we will not focus on a specific model, but rather study the behavior of the drag force in a class of models parametrized by c. Thus, we make c dimensionless by normalizing it to fixed temperatures and express other quantities in units of c. In Ref. [46], it was shown that the range of may be most relevant for a comparison with QCD. We employ this range in the current study.
To proceed, if one uses is used as the radial coordinate, then metric Eq. (13) turns into with 3 Drag force T,µ In this section, we will follow the argument in Refs. [16,17] to analyze the drag force in the SW model. The string dynamic is described by the Nambu-Goto action, given by with g the determinant of the induced metric and g µν X µ with the brane metric and the target space coordinates.
x 1 One considers a heavy quark moving with constant speed v in one direction, (e.g., the direction) and takes the gauge (18) then, the induced metric can be obtained by plugging Eq. (18) into Eq. (14), . (19) Given that, the Lagrangian density reads . As the action does not explicity depend on , the energy-momentum current is a constant, ξ ′ where can be regarded as the configuration of the string tail.
resulting in  The denominator also changes sign at , yielding In contrast, the current density reads (27) and the drag force reads where the minus indicates that the drag force is against the direction of motion, as expected. Next, applying the relations T,µ one arrives at the drag force in the SW model To proceed, we compare the drag force in the SW model with its counterpart of the SYM as the following (31) Fig. 1, we plot against for two fixed values of . The left panel is for and the right one is for . From these figures, one sees that the drag force in the SW model is larger than that of SYM plasma. Furthermore, increasing leads to an enhanced drag force. Namely, the chemical potential enhances the drag force, consistently with [18][19][20]. Comparing the two panels, one finds that increasing enhances the drag force as well. Thus, we can conclude that the confining scale also enhances the drag force, in accordance with Ref. [24]. Fig. 2, we plot versus v for various cases. From the left panel (small case), one finds that increasing v enhances f, and this effect is especially more pronounced for large . However, the right panel (large case) is different: for small , decreases with v, whereas for large the situation is reversed. Interestingly, a similar non-monotone behavior appears in the studies of the drag force with curvature-squared corrections [21].
The effect of c and on the viscosity can be analyzed. A stronger force implies a more strongly coupled medium, closer to an ideal liquid. As the drag force in the SW model is larger than that of SYM, one could infer that the plasma is less viscous in the SW model than in SYM plasma.
µ Recently, Ref. [26] has addressed the drag force in the soft wall model and estimated the spatial string tension at finite T and .

Diffusion coefficient
In this section, we investigate the behavior of the diffusion coefficient in the SW model. First, we remember the results of SYM. The drag force is given by Assuming , Eq. (37) becomes Integrating Eq. (33), one obtains where is the relaxation time of SYM. In contrast, the diffusion coefficient is related to the temperature, the heavy quark mass and the relaxation time as [16,17] Using the above approach, one readily gets T,µ Likewise, from Eq. (30) one can rewrite the drag force in the SW as Chinese Physics C Vol. 43, No. 11 (2019) 113103 In a similar manner, one arrives at the diffusion coefficient in the SW model Fig. 3, we plot as a function of . One can see that and c both reduce the diffusion coefficient. The inclusion of and c may influence the heavy t D = m T D quark mass. As discussed in Ref. [47], the relaxation time should be larger than the inverse temperature µ one can see that and c may enhance the heavy quark mass.

Conclusion
We studied the drag force and diffusion coefficient in a soft wall model with finite temperature and chemical potential. Our aim is to understand the interplay between the presence of confining scale (such that the plasma is not conformal) and the chemical potential (the QGP is assumed to carry a finite, albeit small, baryon number density) when estimating the drag force of a heavy quark in strongly coupled plasma. It is shown that with fixed c, the presence of increases the drag force and decreases the diffusion coefficient, in agreement with the findings of [18][19][20]. Moreover, with fixed the inclusion of c increases the drag force as well, in accordance with Ref. [24]. Therefore, our results ( and c exist at the same time) confirm the results of , or c stands alone. In contrast, the results show that the plasma is less viscous in the SW model than in SYM.

T,µ
Nevertheless, the SW model considered here has several drawbacks. First, it is not a consistent model, T,µ cz 2 since it does not solve the Einstein equations. It would be interesting to provide such analysis in some consistent models, e.g., [48][49][50][51][52] (generally, the metrics of those models are only known numerically, such that the calculations are quite challenging. Further, in these models the dilaton field would be non-trivial, such that the coupling of the dilaton to the word-sheet should be taken into account [53][54][55]). Moreover, the SW model may omit one part in the phase transition, since Refs. [56][57][58][59] argued that there may be a first order phase transition if one sets the warp factor to be and solves the equation of motion self-consistently. Consider these effects would be significant.

T,µ
Finally, it is also of interest to study the jet quenching parameter [60] and stopping distance [61][62][63] in the SW model and compare it with the results of this work. We leave this as a further study.
We would like to extend our gratitude to Dr. Zi-qiang Zhang for helpful and encouraging discussions.