Aspects of quarkonium propagation in a thermal medium as seen by string models
Oleg Andreev



 L.D. Landau Institute for Theoretical Physics, Kosygina 2, 119334 Moscow, Russia
Arnold Sommerfeld Center for Theoretical Physics, LMU-München, Theresienstrasse 37, 80333 München, Germany

 
We use gauge/string duality to model a heavy quark-antiquark pair in a color singlet moving through a thermal plasma. In particular, we explore the effect of velocity on the string tension and Debye screening mass. Then we apply the results to the analysis of heavy quarkonium bound states. With some assumptions, we estimate the characteristic size of quarkonium and its dissociation temperature. 

11.25.Tq, 12.38.Lg, 12.38.Mh
LMU-ASC 05/19

Aspects of quarkonium propagation in a thermal medium as seen by string models
    Oleg Andreev
    December 30, 2023
==============================================================================







§ INTRODUCTION


Heavy-ion collision experiments being performed at BNL and CERN aim to recreate a state of matter known as the quark-gluon plasma (QGP). In this regard, the question arises: how to probe the properties of the plasma formed in such collisions? One of the ideas put forward by Matsui and Satz <cit.> is that quarkonium production rates, in particular that of the J/Ψ mesons, are a good indicator of whether the QGP is produced, and at what temperature.[The literature on quarkonium suppression in   heavy-ion collisions is vast. For a review, see, e.g., <cit.> and references therein.] In this case the understanding of the in-medium behavior of heavy quark bound states is the main ingredient of the analysis with the key feature that the QGP formation leads to a color screening and, as a consequence, no existence of bound states at high enough temperature. 

In reality, the physics is rather complicated, and as yet there is no complete understanding of quarkonium suppression in heavy-ion collisions. To be more accurate, the analysis has to include various confounding effects <cit.>. In this paper, we will focus on one which has a particularly simple origin: quarkonium mesons may be produced moving through a thermal medium. Then the question is what happens to the suppression of those mesons? 

Since the original paper in 1986, the lattice calculations of color screening <cit.> have played a pivotal role in the analysis of bound states. Although lattice QCD is a powerful computational tool to deal with strongly coupled gauge theories, its use is limited to the static case when heavy quark bound states are at rest. Meanwhile the gauge/string duality allows one to at least get a good intuition in situations where evolution in real time is an essential feature, as in the case of interest. This makes it an effective tool to tackle the problem. It is worth noting that such a duality was first proposed by Maldacena for conformal theories and then expanded to include QCD-like theories, and thus somehow relate string theory to heavy-ion collisions <cit.>. Because there still is no string theory dual to QCD, effective string models will only be used for our purposes.

The present paper continues a series of studies on five-dimensional string models started in <cit.>. For orientation, we begin by describing a thermal medium from the viewpoint of a black hole in five dimensions. We go on in section III to discuss general aspects of a quark-antiquark pair propagating in the medium. In section IV, we analyze two important quantities for characterizing the pair: the string tension and Debye screening mass. In doing so, we particularly focus on the velocity dependence of those. Then, in section V, we use the results to make two phenomenological estimates. These include the estimates of the characteristic size of quarkonium and its dissociation temperature in the presence of a thermal wind. Finally, we conclude in section VI. Some technical details are given in the appendices.




§ PRELIMINARIES




As a preliminary to discussing a heavy quark pair propagating in a thermal medium, we will first provide some background information on the medium itself. 

It is well known that the thermodynamics of N = 4 super Yang-Mills theory in four dimensions is related to the thermodynamics of the Schwarzschild black hole in a five-dimensional Anti-de Sitter space <cit.>. For QCD-like theories the situation is more involved, as the geometries become much complicated and often include additional background fields. In this paper we will discuss a pair of effective string models. Our reasons for doing so are: (1) Because it would be very good to gain some intuition on problems in QCD that can be done within the effective string models already at our disposal. (2) Because the string theory dual to QCD remains still unknown. (3) Because the results provided by the models are consistent with the lattice calculations and QCD phenomenology. In particular, this takes place for the equation of state <cit.>, expectation value of the Polyakov loop <cit.>, Debye screening mass <cit.>, spatial string tension and heavy quark diffusion coefficients <cit.>. (4) Because analytic formulas can be obtained by solving these models. (5) Because the aim of our work is to make predictions which may then be tested by means of other methods, e.g., effective field theories or numerically.

The five-dimensional metric (with Euclidean signature) is taken to be


    ds^2=w(r)R^2(f(r)dt^2+dx_i^2+f^-1(r)dr^2)
     ,
      with  
    w(r)=^ r^2/r^2 .

We assume that the blackening factor f is a decreasing function of r on the interval [0,] such that f(0)=1 at the boundary and f()=0 at the horizon. t is a periodic variable of period 1/T. The index i goes from 1 to 3. The boundary of this space, thus, is 𝐑^3×𝐒^1. Note that one can think of this geometry as a one-parameter deformation, parameterized by , of the Schwarzschild black hole on AdS_5 space of radius R. 

Two forms of f have been discussed in the present context. One, as originally postulated in <cit.>[Actually, this was the first attempt to extend the soft wall (metric) model <cit.> to finite temperature.], is given by 


    f(r)=1-(r/)^4
     ,

which is the blackening factor of the Schwarzschild black hole in AdS_5 space. With this choice, the Hawking temperature is simply


    T=1/π√(/h) ,

where h = ^2. This is the simplest case that enables one to write most of the resulting equations analytically. 

Another one was derived from a differential equation which relates the functions w(r) and f(r). It follows from the Einstein equations and can be solved explicitly <cit.>. For w given by (<ref>), one gets 


    f(r)=1-
    1-(1+32 r^2) ^-32 r^2/1-(1+32h) ^-32h .

Given this, the corresponding temperature is 


    T(h)=9/8π√() h^3/2/^32h-1-32h .

These blackening factors lead to results which coincide with each other at high temperatures, but display some differences at lower temperatures, especially near the critical point.

We will use both backgrounds in our analysis, suggesting that the worldsheet path integral for the expectation value of Polyakov loops can be evaluated on some classical solutions which describe the minimal surfaces. At that point there is no requirement on the background metric to be a solution of β_μν=0, with β_μν the worldsheet beta function <cit.>. Of course, such an approximation requires a caveat, because of quantum corrections (effects of order α').  



§ A MOVING QUARK-ANTIQUARK PAIR 




Now we turn to a quark-antiquark pair moving thought a thermal medium. In doing so, we assume that the quarks are heavy enough such that there is no relative motion between them. It is also worth keeping in mind one fact: in the dual formulation we are using the pair does not feel a drag force and, consequently, does not lose its energy <cit.>. In other words, a string configuration is moving uniformly through the medium. Of course, it is an assumption. In the reality there is energy lost that makes string configurations time dependent and, as a consequence, the analysis becomes extremely tedious and complicated. 

In terms of strings in five dimensions, the pair is described by a Nambu-Goto string attached to quark sources on the boundary. We place the quarks on the x-axis and then consider their uniform motion with velocity v in the z-direction, as sketched in Figure <ref>. The action is the sum of the string action and those of two point-like particles. For our purposes in this paper, what we need to know is how the energy of this system depends on the quark separation. This can be determined entirely from the Nambu-Goto action 


    S=-1/2πα'∫ d^2ξ √(-γ_αβ) ,

where α' is a string constant, (ξ_1,ξ_2) are worldsheet coordinates, and γ_αβ is an induced metric. 

In studying bound state properties, it is more convenient to choose the rest frame of the pair. In that case, the medium takes on the role of a hot wind blowing on the static pair. To describe such a wind geometrically, one first Wick rotates to Minkowski space and then boosts the metric in the z-direction <cit.>. With t=γ(t'+vz') and z=γ(z'+vt'), a simple algebra shows that the metric (<ref>) takes the form[To avoid clutter, we suppress the primes as we will always be interested in the pair rest (primed) frame.]   


    ds^2=w(r)R^2
    (- dt^2+2γ^2(1-f)dt dz
    +
    dx^2+dy^2+γ^2(1-f^2)dz^2
    +
     f^-1dr^2)
     ,

where γ=1/√(1-v^2). For convenience, we introduce a new blackening factor 


    (r)=γ^2(f-^2)
     .

This metric describes the hot wind felt by the static pair positioned at a right angle to the direction of the wind.




 §.§ A connected string configuration


We now want to consider a connected string configuration. It is similar to the configuration sketched in Figure <ref>, but with  replaced by . The analysis proceeds along the lines described in <cit.>.

We start by taking the static gauge t=ξ_1 and x=a ξ_2+b. Combining the latter with the boundary conditions at the string endpoints x(0)=-ℓ/2 and x(1)=ℓ/2, we obtain a=ℓ and b=-ℓ/2. Since the configuration is static, the lowest energy corresponds to z=const, and therefore r is a function of x only, as in the AdS case <cit.>. For such a configuration, the induced metric is given by 


    ds^2=w(r)R^2
    (- dt^2+(1+f^-1r'^2)dx^2)
     .

This metric has a horizon at r=, where the equation f(r)=v^2 has a solution. We call it an induced horizon as a shorthand for saying that it is related to the induced metric on the string worldsheet. Note that, since f is a decreasing function of r such that f()=0, it holds < for all v≠0. Also, it is useful for further analysis to introduce the notion of an effective string tension 


    (r)=^ r^2/r^2√() ,
 
and think of it as a function of r.

The Nambu-Goto action now becomes


    S=-∫ dt ∫_-ℓ/2^ℓ/2 dx √(1+f^-1r'^2) ,

where =R^2/2πα' and r' means ∂ r/∂ x. Since the Lagrangian does not depend on x explicitly, there is a conserved quantity. It is nothing else but the first integral of the equation of motion for r 


    I=/√(1+f^-1r'^2) .

On symmetry grounds, r'(0)=0 that allows one to express I in terms of  as I=(0̊), where r(0)=0̊. 

The string length along the x-axis can be found by further integrating (<ref>). Writing it as a relation between dx and dr and then performing the integral over [-ℓ/2,ℓ/2] over x, we obtain


    ℓ=2∫_0^0̊dr/√(f) [ ^2/I^2-1 ]^- ,

where the factor 2 comes from the reflectional symmetry of the configuration.

Given a string configuration (solution), we can compute its energy. Since the string is static, the energy is simply related to the Lagrangian. So, 


    E_=2∫_0^0̊dr/√(f) [ 1-I^2/^2 ]^- .

Here we used (<ref>) to reduce the integral over x to an integral over r. This expression is not well-defined, because the integral diverges at r=0. We regularize it by imposing a cutoff ϵ such that r≥ϵ. Then in the limit ϵ→ 0 the regularized expression behaves like 


    E_^ R=2/ϵ+E_+O(ϵ)
     .

Subtracting the 1/ϵ term and letting ϵ=0, we get a finite result


    E_=2∫_0^0̊ dr(/√(f)[ 1-I^2/^2 ]^--1/r^2)  -2/0̊+2c
     ,

where c is a normalization constant. Thus, the energy of the string is given by the parametric equations (<ref>) and (<ref>), with 0̊ a parameter.[When written in this form, the formulas are also applicable to other backgrounds, like those with w=^A(r)/r^2.] It can take values in the interval [0,]. The reason for this is that both E_ and E_ (as defined below) must be real-valued. As an immediate comment, we note that for v=0, the formulas are reduced to those of <cit.> obtained for the string free energy. 

It is now straightforward to investigate the properties of E_ at short distances. In doing so, an important fact will be that small ℓ's correspond to small 0̊'s. This means that a short string stays close to the boundary. The effective string tension and blackening factor behave for r→ 0 as 
∼ 1/r^2 and f∼ 1. Therefore, in this limit, the behavior of E_ is the same as in the AdS case at zero temperature and velocity <cit.>. We have


    E_=-α_/ℓ+2c+o(1)
    	 ,

with α_=(2π)^3/Γ^4(14), and this holds for a large class of backgrounds.

Let us now investigate the long distance behavior of E_ that is achievable with a little more effort. Following <cit.>, we start with the effective string tension. There is one very simple but important fact: the form of , as a function of r, is temperature and velocity dependent.[For the Reissner-Nordström blackhole whose blackening factor reduces to that of <cit.> at zero chemical potential, it was recently discussed in <cit.>.] At the qualitative level, it can be understood as follows. For small T and v, where f≈ 1,  has a local minimum near r=1/√() defined by the warp factor. For large enough T and v, where ≪ 1/√(), there is no local minimum as both factors are decreasing functions in the interval [0,]. This is illustrated in Figure <ref>.



We return to this issue in Section IV, after describing the behavior of E_ in the infrared region.

Such a behavior of  makes it more or less obvious that the moving string will have similar infrared properties as that at rest <cit.>. Indeed, if the effective tension has a local minimum at r=, then the string cannot get deeper than  in r direction because a repulsive force prevents it from doing so. One can think of this as a soft wall located at r =. In this case, the large distance behavior of the string is completely determined by the wall. Concretely, a short inspection shows that ℓ(0̊) and E_(0̊) are monotonically increasing in the interval [0,]. They become infinite at 0̊=, where


    ℓ(0̊)=-kln( -0̊) +O(1)
     ,    
    E_=-kσln( -0̊)+O(1)
     ,

with coefficients

    k=2√(/f”)() ,
        σ=()
     .

It follows from this that 


    E_=σℓ+O(1)
    	 .

Thus, this phase can indeed be interpreted as the phase of confinement. Note that σ is the physical string tension.

On the other hand, if the effective string tension has no local minimum, then the string can get deeper into the bulk and finally reach the induced horizon. Its large distance behavior now is determined by the near horizon geometry. The crucial point here is that the effective string tension vanishes on the induced horizon. So it seems reasonable to expect that there is no linear term in the expression for E_ and therefore interpret this phase as the phase of deconfinement. However, the reality does not meet the expectation. The problem is that the connected configuration exists only if ℓ does not exceed the critical value .[This fact motivates the use of an operational definition of the screening length ℓ_s= <cit.>. However, it remains to be seen to what extent such a length will be consistent with the Debye screening length of lattice QCD <cit.>.] One way to address this problem is to introduce one more diagram as suggested in <cit.>. But we will not do so. Instead, we will discuss the Debye screening mass along the lines of <cit.>.

We close this subsection with one example which illustrates the difference between the behavior of E_ in the confined and deconfined phases. In Figure <ref> we plot 



E_ as a function of ℓ. Although we restrict ourselves here to the blackening factor (<ref>), this picture holds for (<ref>) as well. 




 §.§ A disconnected string configuration


Since a disconnected string configuration is so important in the deconfined phase, we will now study it. This configuration is similar to that of Figure <ref> on the right, except that  is replaced by . It consists of two similar parts, each of which represents a single quark (antiquark) placed on the boundary and connected to the induced horizon by a Nambu-Goto string. 

We want, as in subsection A, to find the total energy of those strings which is twice larger than the energy of each string. The latter can be computed along the lines of <cit.>. Since the configuration in static, we choose the static gauge such that t=ξ_1 and r=aξ_2 +b. From the boundary conditions r(0)=0 and r(1)=, we have a= and b=0. The Nambu-Goto action is then


    S=-∫ dt ∫_0^ dr √(x'^2+f^-1) ,

where a prime stands for a derivative with respect to r. It is obvious that x(r)=const is a solution to the equation of motion that represents a straight string stretched between the boundary and the induced horizon. Moreover, this solution corresponds to the lowest string energy. 

In proceeding in this way, the important fact can be obscured by the gauge condition r=aξ_2+b. So it is instructive to consider another, namely x=aξ_2+b. Although this gauge condition is not convenient for studying the straight string, it is however helpful in understanding the force balance at the string endpoint located on the induced horizon. After a simple calculation, we find a boundary term in the variation of the Nambu-Goto action  


    δ S_b=-∫ dt  r'/√(1+f^-1r'^2) δ r |_0^ .

A prime now means a derivative with respect to x. The integrand provides the string tension in the r-direction.[By contrast, the string tension in the x-direction can easily be found by analyzing the action (<ref>).] For the straight string it is simply σ_r=-g√(f). This tension vanishes on the induced horizon that implies that there is no need for an external force to keep the string endpoint in place. In other words, the configuration is in mechanical equilibrium.

The string energy is simply obtained by evaluating the Lagrangian in Eq.(<ref>) on the solution. We have


    E_=∫_0^dr/√(f) .

The integral is is divergent at r=0. As before, we regularize it by imposing the cutoff ϵ. Therefore, the regularized expression behaves for ϵ→ 0 as


    E_^ R=/ϵ+E_+O(ϵ)
     .

After subtracting out the 1/ϵ term and letting ϵ=0, we get 


    E_=∫_0^dr(/√(f)-1/r^2)  -/+c
     ,

with the same normalization constant as in (<ref>). Clearly, E_ is also given by this expression. Note that for v=0, it reduces to the formula obtained in <cit.> for the heavy quark free energy. As such, E_ can be interpreted as a heavy quark energy in the presence of a hot wind.





§ MORE DETAIL ON THE STRING TENSION AND DEBYE SCREENING MASS




Here we will describe more precisely important subtleties that arise when a heavy quark-antiquark pair moves though the thermal medium. Equivalently, these are related to the effects of the hot wind described in section III. 

We start with the simplest model in which the blackening factor is given by (<ref>). In this case we can explicitly write down the expressions for the induced blackening factor 
and horizon 

    =1-r^4/^4 ,
        =√(1-^2) .

This allows a drastic simplification of the problem of determining local extrema of the function (r). A simple argument using the properties of the cubic equation shows that it depends on a parameter ^2. The extrema exist only if ^2 >√(27)/2 and occur at 


    =√(2/√(3)sinϕ) ,    =√(2/√(3)sin(π/3-ϕ)) ,

where ϕ=1/3arcsinγT^2/T_pc^2 and T_pc=1/π√(2/3√(3)). Accordingly, the condition for confinement can now be written as 


    √(γ) T/T_pc<1
     .

Thus v goes to zero as T approaches T_pc.

Now let us pause here for a moment to make a couple of remarks. First, T_pc is a temperature which was originally introduced in <cit.>. Although it can be interpreted as a critical temperature of the model, this requires a caveat: in contrast to lattice QCD <cit.>, the string tension does not vanish at T=T_pc (see the left panel of Figure <ref>). Second, it is unclear to us how to extend the above analysis to include the blackening factor (<ref>). Therefore we will examine this case numerically. 



 §.§ The string tension


Using the formula (<ref>) for the string tension and the expression (<ref>) for , we deduce a formula for the model with the blackening factor (<ref>)


    σ(,T)=1/3σγT^2/T_pc^2√(sin^-2ϕ-4/3) exp(3γ^-1T^2_pc/T^2sinϕ -1)
     ,

where σ= is the string tension at T==0 <cit.>. This result may be obtained in a simpler manner by rescaling T by a factor of √(γ) in the corresponding formula of <cit.>. It seems to be easy, but it is not without pitfalls, as we will see in the Appendix B. Notice that on the critical curve √(γ)T/T_pc=1 the tension is non-vanishing and approximately equal to 0.9σ. 

 It is instructive to examine the low-temperature behavior of σ(v,T). In view of Eq.(<ref>), we expand it in powers of √(γ)T/T_pc. Then its expansion reads[It is worth mentioning that the leading low-temperature correction found within the four-dimensional string model <cit.> at zero velocity is proportional to T^2.]


    σ(v,T)/σ≊ 1-2/27γ^2T^4/T_pc^4 .

Here we drop the higher order terms. It turns out that this is a very good approximation for √(γ)T/T_pc<0.6. For larger values the discrepancy becomes clear and reaches a maximal value on the critical curve, where it is of order 3%.

To complete our discussion of this model we plot the string tension as a function of temperature and velocity. From the left panel of Figure <ref> we see that σ(v,T) decreases with T and v. Moreover, it is a slowly varying function except 



near the critical curve (especially for high velocities). For T=0  the tension is independent of , as expected due to Lorentz invariance.

By contrast, algebra does not help much in understanding the model in which the blackening factor is given by (<ref>). In that case, numerics is the only way. There are two special features to note in this example. The first is that T_pc≈ 2· 10^-3√() which is two orders lower than that in (<ref>).[For comparison, note that the estimate T_c=√()/π made in <cit.> for the critical temperature seems more in line with the lattice value given in <cit.>.] The second is that f does not go to 1 as T goes to 0. In other words, this model does not reduce to that of <cit.> at T=0. Therefore, in practice, we can study this model only in a narrow temperature interval just above zero temperature. We take it from 4· 10^-7T_pc to T_pc. For v=0, the string tension is approximately σ_n=1.63 at the lower limit. In the right panel of Figure <ref>, we plot the result for the string tension. We see that σ(v,T) is a slowly varying function of temperature as before, but now it shows a steep decline as a function of velocity. In fact, the string tension vanishes near v≈ 0.27. We return to this in section V, after describing the Debye screening mass and gaining some information about its behavior in the temperature range close to T_c.

It is noteworthy that what we have discussed above is applicable only for confining backgrounds (confining gauge theories). For non-confining ones, including the AdS Schwarzschild black hole, it makes no sense as the string tension is vanishing. 




 §.§ The Debye screening mass


By now we understand that at zero velocity the Debye screening mass can be computed from the correlator of two Polyakov loops <cit.>. Moreover, it can be done for both channels (singlet and octet) with the same result. Here there is an important issue that arises when one tries to split the correlator in two parts non-perturbatively.[This issue is still under debate <cit.>.] We won't discuss it here, in part because in string models the octet channel is not understood even on the very basic level, and in part because for our purposes it is enough, as shown before for v=0 <cit.>, to use the singlet channel.  

From here, we consider only the singlet channel. This means that the quark-antiquark pair is in the color singlet state. For this channel the string configuration is given by the connected configuration of section III. We define a binding energy of the pair as a difference between the energies of the connected and disconnected string configurations


    Δ E_=E_-E_-E_ ,

where E_ and E_(E_) are given by Eqs.(<ref>) and (<ref>), respectively. The above definition implies that the quark contributions, if any, cancel each other. In particular, this happens for the normalization coefficient c. Thus, the binding energy is expressed in terms of the string energies only. In the limit v→ 0, it reduces to Δ F_=F_-F_-F_, with F's free energies. Δ F_ plays a special role in determining the Debye screening mass via the singlet channel <cit.>. For more explanations, see the Appendix A.

Now we want to find the behavior of the binding energy at short distances. From Eqs.(<ref>) and (<ref>), we immediately obtain that 


    Δ E_=-α_/ℓ+
    2[1/+
    ∫_0^dr(1/r^2-/√(f))
    ]  +o(1)
     .

Motivated by the analysis in the Appendix A, we parameterize the binding energy as


    Δ E_=-4α/3ℓ^-mℓ+E_0
     ,

where E_0 is an ℓ-independent constant which reduces to that of Eq.(<ref>) as v→ 0. 

Expanding both sides of this equation in powers of ℓ and keeping only two leading terms in the expansion, at leading order we obtain α_=4/3α which is the same relation as in the static case. But at the next order we have that


    m=Γ^4(14)/4π^3 [1/+
    ∫_0^dr(1/r^2-/√(f))
    +𝐜]
     ,

with 𝐜=-E_0/2. Once we set the velocity to zero, the integral can be performed with the final result given by (<ref>). 

One can think of E_0 as encoding information about quantum corrections to the classical energy of the string as well as about possible ambiguities in lattice QCD (at v=0). To this we have nothing to add, and will proceed with a couple of operational definitions for the screening mass. Consider the simplest possible forms of 𝐜. On dimensional grounds, 𝐜 may be proportional either to T or to 1/. In the first case the simplest option is to take 


    𝐜=𝔴T
     ,

while in the second


    𝐜=
    𝔴/·π^-1     for    (<ref>) ,
    
     T
        for    (<ref>) .

In both cases 𝔴 is a proportionality constant which coincides with that of Eq.(<ref>). So, we treat it as the model parameter. 

In contrast to the static case, we are unfortunately unable to perform the integral in (<ref>) analytically. But it is still possible to do so in certain limiting cases. For f defined by (<ref>), we take the limit v→ 1, getting 


    √(1-v^2) m/T=C
     ,

where the constant C is given by 

    C=Γ^2(14)/√(2π)·
    1     for    (<ref>) ,
    
    
    1+Γ^2(14)/(2π)^3/2𝔴/π    for    (<ref>) .

For the blacking factor (<ref>), the story is less simple. To get the behavior (<ref>), one has to take a double limit in which v→ 1 and T→∞. The latter condition implies that both f's coincide. Notice that a similar formula arises in hot N = 4 supersymmetric QCD <cit.> without any restrictions on v and T. 

We go on further with a SU(3) pure gauge theory, keeping in mind our goal of getting some intuition for phenomenology. Like in <cit.>, we write the critical temperature T_c=0.629√(σ) <cit.> as T_c=0.435√(), where we used σ= and in the last step set =0.176 <cit.>. Then, the condition T>T_c automatically implies that the models are in the deconfined phase as T_c>T_pc in both cases. In Figure <ref> we plot the screening mass versus temperature and velocity. 



For none-zero velocity, it is convenient to normalize it also by a factor of √(γ). [Thus, √(1-v^2)m/T∼m/m_AdS, as follows from (<ref>).] We set 𝔴=-2.05 to fit the lattice data at v=0 <cit.>. A common feature of these string models is that they fail to reproduce the lattice results in the near vicinity of the critical point as well as for high temperatures. Because of this, we take a small cutoff at T=T_c and restrict T to be less than 5T_c. We see that the Debye screening mass is an increasing function of both variables, except for v≳ 0.72 in the model with (<ref>). It increases steeply near the edges T=T_c and v=1. For non-relativistic velocities, the difference due to the two forms of 𝐜 is negligible. As we will see in the next section, charmonium bound states may propagate through plasma only with low (non-relativistic) velocities, while bottomonium with v≲ 0.82. Therefore, the behavior of m near v=1 could matter only for the analysis of topomonium bound states that seems hypothetical.
 



§ SIMPLE ESTIMATES




Having understood how a wind velocity affects the string tension and Debye screening mass, we can formally redo the estimates of <cit.> for quarkonia moving through a thermal medium. In this section, we will do so in the context of finite temperature pure SU(3) gauge theory. We assume that a relative motion of quarks is non-relativistic in the boosted (primed) frame.

We begin with an estimate of the characteristic size of quarkonium (ground state) in the confined phase at low temperature. Taking the quark-antiquark potential of the Cornell form V(r)=-α_/r+σ r +V_0 and minimizing the ground state energy with respect to r, we see that it has a local minimum at 


    ρ(v,T)=
    (1+√(1-4μ^2/27σα^3_)/2μσ )^1/3
    +
    α_/3σ(1+√(1-4μ^2/27σα^3_)/2μσ )^-1/3 ,

with μ the reduced mass. The above expression is well defined for σ≥4μ^2/27α_^3 that yields σ≥ 9.67× 10^-4  GeV^2 for charmonium states and σ≥ 1.05× 10^-2  GeV^2 for bottonium states. These values are obtained by setting α_=0.253, as follows from (<ref>) at =0.176, along with m_c=1.27 GeV and m_b=4.18 GeV <cit.>. This means in practice that in the vicinity of the critical line the above estimate is not reliable.

For illustrative purposes, in Figure <ref> we show the plots of ρ(v,T) for both forms of σ(v,T) discussed in section IV. Here we  



assume that the string tension mildly depends on the wind direction. In the first case, when f is given by (<ref>), ρ(v,T) can be found analytically by using the expression (<ref>). The resulting formula is rather cumbersome to work with, so it seems natural to look for a simple approximation. We found that it can be well approximated by 1.25(μσ(v,T))^-1/3. So, just like the string tension, the size of the ground state slowly varies with velocity and temperature, as also seen from the left panel of the figure. In the second case, with f coming from (<ref>), the restrictions of section IV remain and, in addition, it is now prohibited to approach close to the critical line where the approximation (<ref>) breaks down. We are unable to find a simple approximate formula for all velocities, but for low ones it is as described above. It is seen from the right panel that ρ(v,T) increases rather slowly with temperature and steeply with velocity.

What we have learned above is that the characteristic size of quarkonium states increases with increasing wind velocity. So, it is natural to expect that quarkonium states will now dissociate at lower temperature. Our remaining goal will be to give another explanation of this and to make a simple estimate of the dissociation temperature.

First, let us recall one fact about the Debye screening mass. An effective parameterization for the lattice results is given by 


    m/T=a-b T_c^2/T^2 .

It is motivated by the results obtained within the models of section II <cit.>. For SU(3) pure gauge theory a=0.71 and b=0.24. Such a parameterization is quite good in the temperature range 1.06  T_c≲ T≲ 3 T_c. 

In the deconfined phase, the analysis proceeds as before, but this time the heavy quark potential is taken to be of the form V(r)=-α_/r^-mr+V_0. It can be shown that a bound state exists if m/μα_≤ 0.84 <cit.>. Using the above expression for m, the dissociation temperature at zero velocity is 


    =T_c/a
    (0.42 μα_/T_c+
    √(ab+(0.42 μα_/T_c)^2) )
     .

Let us make a simple estimate of . With the same parameter values as quoted above, the dissociation temperature is then =1.05  T_c (cc̅), =1.37 T_c (cb̅), and =2.54  T_c (bb̅). As expected, bottomonium states dissociate at higher temperature. 

The formula (<ref>) for m can be interpreted similarly to that for the pressure <cit.>. In a temperature range between the critical temperature and some characteristic temperature of perturbative QCD, the ratio m/T is given by a series in powers of 1/T^2. In the absence of the 1/T^2 corrections this formula reduces to m/T=a, which is also known in the context of AdS/CFT <cit.>. In this case the estimate gives =0.73T_c (cc̅), =1.12T_c (cb̅), and =2.41T_c (bb̅). This indicates that a reasonable formula has to include not just a constant term but also corrections. Notice, however, that the bottonium is less sensitive (at least) to the leading correction because of its large mass. Indeed, it follows from (<ref>) that =0.84μα_/a as μ→∞. 
 
The generalization to non-zero velocities does not change much of what we have said above as it follows from the analysis in the Appendix B. We start with a formula similar to Eq.(<ref>) 

    m/T= a(v)-b(v) T_c^2/T^2 ,

but this time with velocity-dependent coefficients which reduce to those of (<ref>) at v=0. 

Since the Debye mass increases with increasing velocity, it is natural to expect that the dissociation temperature will decrease with it. As an example of this, consider the simplest model, i.e. that with (<ref>). In this case, the coefficients take the form 


    a(v)=m_0(v)
    	 ,    
    	b(v)=b m_1(v)/m_1(0) ,

where the m's are given by (<ref>). The value of 𝔴 is determined from the condition a(0)=a. With a=0.71, it is 𝔴=-2.63. A simple but crude approximation to (<ref>) is to assume that 


    a(v)=aγ^ ,    
    	b(v)=bγ^- .

This follows from a rescaling argument (see the Appendix B). We have also assumed that the Debye screening mass mildly depends on the wind direction.

Now, using Eqs.(<ref>)-(<ref>), we can make a simple estimate of the dissociation temperature at finite velocity. In Figure <ref>, we plot the results  



for charmonium and bottomonium states. We see that in the presence of a thermal wind  quarkonium states dissociate at lower temperature, as expected. In particular, charmonium states turn out to be more sensitive to the wind than bottomonium ones. This suggests that charmonium states may propagate through the plasma only with non-relativistic velocities. From the left panel, we obtain an upper bound v≲ 0.13. For these low velocities, the difference due to the two forms of 𝐜 turns out to be unimportant.

Finally, let us note that the approximation (<ref>) results in a very simple formula for the dissociation temperature: 


    (v)=√(1-v^2) (0)
     .

It is excellent in the context of AdS/CFT <cit.>, but rather crude in comparison with ours (more phenomenologically motivated). 




§ CONCLUDING COMMENTS




The purpose of this paper has been to further develop the five-dimensional effective string models as a theoretical tool for studying strongly coupled gauge theories. In particular, we feel that it could be a correct way for dealing with situations in QCD where evolution in real time is an essential feature. This is an optimistic view. Of course, there are many open questions left which deserve further study, and the main question: what is the string dual to QCD?

 Of particular interest for heavy-ion phenomenology is those backgrounds which describe the medium not only at finite temperature but also at finite baryon density. Furthermore, there is an important aspect related to the anisotropic expansion of the quark-gluon plasma. We believe that these and similar issues are worthy of future study.



We would like to thank I. Aref'eva, P. de Forcrand, R.R. Metsaev, P. Weisz, and U.A. Wiedemann for useful and encouraging discussions. We also thank the Arnold Sommerfeld Center for Theoretical Physics and CERN Theory Division for the warm hospitality. This work was supported in part by RFBR Grant 18-02-40069.



 



§ THE DEBYE SCREENING MASS




One of the features of the quark-gluon plasma is the phenomenon of Debye screening. Its numerical characteristic is a screening mass which provides the scale over which interactions are effective. Here we briefly recall how the five-dimensional string models can be used to estimate the Debye mass along the lines of lattice QCD <cit.>.

In lattice QCD, the Debye mass can be determined from the exponential fall-off of correlation functions at asymptotically large separation. The simplest possible correlator is that of two Polyakov loops <cit.>. Since a quark-antiquark pair could be either in a singlet or in an octet state, it seems natural to write the correlator as a sum of the singlet and octet contributions 


    ⟨  L^†(x_1) L(x_2) ⟩/|⟨  L ⟩|^2=1/9^-F^1_T+8/9^-F^8_T .
 
Here F^1_ and F^8_ are the singlet and octet free energies, respectively. 

One parameterization of the F's, motivated by high temperature perturbation theory, is simply


    F^1_(ℓ)=-4α/3ℓ^-mℓ ,    
    F^8_(ℓ)=α/6ℓ^-mℓ ,

where ℓ=| x_1-x_2|. Obviously, for small ℓ the singlet channel dominates and therefore determines the correlator 

 
    ⟨  L^†(x_1) L(x_2) ⟩/|⟨  L ⟩|^2=1/9^-F^1_T ,    with    
    F^1_(ℓ)=-4α/3ℓ+4/3α m+o(1)
     .


In discussing the correlator within five-dimensional effective models, one has to evaluate a worldsheet path integral obeying the boundary condition that the string worldsheet has the loops for its boundary. In practice, the integral can only be evaluated semiclassically in terms of minimal surfaces. The result is written as 


    ⟨  L^†(x_1) L(x_2) ⟩=∑_n w_n ^-S^(n) .

Here S^(n) is a Nambu-Goto action evaluated on a classical solution (string configuration) whose relative weight is w_n. For small ℓ, the configuration depicted in the left panel of Figure <ref> dominates. 



In this case S^(1) behaves as S^(1)=1T(-α_/ℓ+2c+o(1)), with a universal Coulomb coefficient.[This statement can be expressed as the space becomes asymptotically AdS near its boundary.] Such a behavior allows one, after fitting the Coulomb coefficients, to associate configuration (1) with the singlet channel. The normalization is provided by configuration (2) shown in the right panel. So, we have
 
 
    ⟨  L^†(x_1) L(x_2) ⟩/|⟨  L ⟩|^2=
    ω ^-Δ F_T .

Here Δ F_=F_-F_-F_, with F_=TS^(1) and F_=F_=TS^(2). The factor ω is a ratio of the weights.

Equating the expressions (<ref>) and (<ref>), we get 


    F^1_=Δ F_-Tln(9ω)
      .

If we expand both sides of this equation in powers of ℓ and equate coefficients, we find α_=4/3α at leading order, while at the next order we have α_m=-2F_-Tln(9ω). Taking w=^ r^2r^2 for the warp factor and using the explicit expression for F_, we arrive at


    m=Γ^4(14)/4π^3(
    ^h/-√(π) (√(h) )+𝔴T
    )
     ,

with 𝔴=-ln(9ω)/2. Note that the screening mass does not depend on the actions for the heavy quark sources, and in particular on the normalization constant c. For our purposes, we think of m as a function of T, and treat 𝔴 as a model parameter. The reasoning for this is that at the moment a detailed calculation of the weight factors is beyond our abilities.




§ TAYLOR SERIES OF M




Our goal here is to compute the first few coefficients of the Taylor series of m/T(h) about h=0. This provides the basis for making the estimate of the dissociation temperature in Section V. The higher order coefficients can be computed in a similar way.

For both forms of f and 𝐜 we see explicitly that m/T is regular at h=0. So we expand it in powers of h 


    m/T(h)=m_0-m_1 h+O(h^2)
    	 .

This formula is completely equivalent to our earlier formula (<ref>), with 1/T^2 expansion.

First, we consider the model whose blackening factor is of the form (<ref>). After expending it in powers of h and performing the corresponding integrals over r, we get


    m_0=
    Γ^4(14)/(2π)^2γ^(
    C_ A
    + v^2 B(34,12)
    F(32;34,54;γ^-2)
    )
     ,
        
    m_1=Γ^6(14)/6(2π)^5/2γ^-F(12;14,74;γ^-2)
     ,

with

    C_ A(v)=
    𝔴/π·γ^-   for  (<ref>) ,
    
    1    for  (<ref>) .

Here F(a; b, c; x) and B(a,b) are the hypergeometric and beta functions <cit.>.[As usual, the subscripts in _2F_1 are omitted.] In the limit v→ 1 the coefficients become


    m_0=Cγ^ ,
        
    m_1=0
     ,

as expected from (<ref>). Here we used F(a; b, c; 0)=1 and B(34,12)=(2π)^3/2Γ^2(1/4).

The v-dependence of these coefficients merits some comment. At first sight, it may seem that it is completely determined by Eqs.(<ref>) and (<ref>) so that the result may be obtained from that of <cit.> by rescaling T→√(γ)T. This is partially true, but substantially false as seen from the explicit expressions of the coefficients. The reasons for this are as follows. First, the form of the metric (<ref>) is not simple in sense that it includes not only , but also f and v. While the former does lead to the rescaling, the others do not. Second, the uncertainty related to the form of 𝐜 may cause, as seen from (<ref>).  

For the model defined by (<ref>), the leading coefficient is given by (<ref>), while the next to leading is 


    m_1=
    Γ^4(14)/16π^2γ^-(
    γ C_ B
    +
    B(14,32)
     F(12;14,74;γ^-2)
    +v^2/2B(54,12)
     F(32;54,74;γ^-2)
    -3/2v^2γ
    B(34,12)
     F(32;34,54;γ^-2)
    )
     ,

with

    C_ B(v)=𝔴/π·
    0    for  (<ref>) ,
    γ^-1-1    for  (<ref>) .

In this form we can easily take the limit v→ 1 so that we find  


    m_1=
    Γ^4(14)/16π^2γ^(C_ B(1)-3/2B(34,12))
     .

Thus, if 𝔴=-2.05, then m_1 is negative. This implies that for relativistic velocities the ratio m/T decreases with the increase of temperature, as also seen from Figure <ref>. 

We conclude this appendix with some further remarks on the v-dependence of the coefficients. To understand deviation from the scaling T→√(γ)T, it is convenient to introduce two test functions 


    m̅_0(v)=m_0γ^- ,    m̅_1(v)=m_1γ^ .

In Figure <ref> we plot these functions normalized to their values at v=0.  



For both C_ A, m̅_0 is regular in the interval. It increases more rapidly for C_ A defined by (<ref>). By contrast, m̅_1 is regular for (<ref>) but singular for (<ref>), where m_1∼γ as follows from the asymptotic formula (<ref>). It decreases with v, although very slowly for (<ref>). Thus, we see that treating both functions as if they were constants may be a rather crude approximation. 



99

satz
T. Matsui and H. Satz, Phys.Lett.B 178, 416 (1986).
bram-rev
N. Brambilla et al., Eur.Phys.J.C 74, 2981 (2014).
uaw-book
J. Casalderrey-Solana, H. Liu, D. Mateos, K. Rajagopal, and U.A. Wiedemann, Gauge/String Duality, Hot QCD and Heavy Ion Collisions, Cambridge University Press, 2014.
satz-book
H. Satz, Lect. Notes Phys. 945 (2018) 1.
 lattice-rev
P. Petreczky, J.Phys. G39, 093002 (2012).
ssg
S.S. Gubser, Found.Phys. 43, 140 (2013).


 az3
O. Andreev and V.I. Zakharov, J. High Energy Phys. 0704 (2007) 100.
a-screen
O. Andreev, Phys.Rev.D 94, 126003 (2016).

wittenT
E. Witten, Adv.Theor.Math.Phys. 2, 505 (1998).
a-T2
O. Andreev, Phys.Rev.D 76, 087702 (2007).
a-pol
O. Andreev, Phys.Rev.Lett. 102, 212001 (2009).
a-D
O. Andreev, Mod.Phys.Lett.A 33, 1850041 (2018).
az2
O. Andreev and V.I. Zakharov, Phys.Lett.B 645, 437 (2007).
soft
R.R. Metsaev, IIB supergravity and various aspects of light cone formalism in AdS space-time, hep-th/0002008; A. Karch, E. Katz, D.T. Son, and M.A. Stephanov, Phys.Rev.D 74 (2006) 015005; O. Andreev,  Phys.Rev.D 73 (2006) 107901.
noro
S.I. Finazzo and J. Noronha, Phys.Rev.D 90, 115028 (2014).
gsw
M.B. Green, J.H. Schwarz, and E. Witten, Superstring Theory, Volume I, Cambridge University Press, 1987.

arg
P.C. Argyres, M. Edalati, and J.F. Vazquez-Poritz, J. High Energy Phys. 0701 (2007) 105.
uaw-wind
H. Liu, K. Rajagopal, and U.A. Wiedemann, Phys.Rev.Lett. 98, 182301 (2007).
chernic
M. Chernicoff, J.A. Garcia, and A. Guijosa, J. High Energy Phys. 0609 (2006) 068.
malda
J.M. Maldacena, Phys.Rev.Lett. 80, 4859 (1998).
ch
X. Chen, S.-Q. Feng, Ya-F. Shi, and Y. Zhong, Phys.Rev.D 97, 066015 (2018). 
bak
	D. Bak, A. Karch, and L.G. Yaffe, J. High Energy Phys. 0708 (2007) 049.
az1
O. Andreev and V.I. Zakharov, Phys.Rev.D 74, 025023 (2006).
pis
R.D. Pisarski and O. Alvarez, Phys.Rev.D 26, 3735 (1982).
boyd
G. Boyd, J. Engels, F. Karsch, E. Laermann, C. Legeland, M. Lutgemeier, and B. Petersson, Nucl.Phys.B 469, 419 (1996).
octet
O. Jahn and O. Philipsen, Phys.Rev.D 70, 074504 (2004); G.C. Rossi and M. Testa, PoS LATTICE 2013, 369 (2014).
3q
O. Andreev, Phys.Rev.D 93, 105014 (2016).
pdg
C. Patrignani et al. (Particle Data Group), Chin.Phys.C, 40, 100001 (2016). 
pis-fb
R.D. Pisarski, Phys.Rev.D 74, 121703 (2006).









 gr
I.S. Gradshteyn and I.M. Ryzhik, Table of Integrals, Series, and Products, Academic Press, 1994. 



