Jet suppression in non-conformal plasma using AdS/CFT

In this paper, we study suppression of light quark in strongly coupled non-conformal plasmas using the AdS/CFT correspondence. The well-known falling string profile in the bulk is considered as light quark moving through the plasma. The maximum distance which string with energy E can travel before falling through the horizon is interpreted as thermalization distance of light quark in the hot-strongly coupled plasma. Our numerical results show that the thermalization distance of light quark increases by increasing deviation from conformal invariance. The relation between this distance and the energy of quark and the temperature of the plasma is analyzed numerically. The jet quenching parameter is also calculated in the non-conformal backgrounds and it is found that the jet quenching parameter is generally decreased by increasing the non-conformality. Our results are compared with the results of N = 4 SYM theory and also some available experimental data.

is not a fully conformal fluid and the bulk viscosity (a purely non-conformal effect) is needed for highly-precise extraction of the shear viscosity of the QGP [38]. Furthermore, hydrodynamics including non-conformal effects successfully described the smaller system such as p-Pb [39] and p-p [40,41] collisions [42].
There are many attempts to extend the original duality to the theories close to QCD or QGP either uses a top-down [43][44][45] or a bottom-up [46,47] construction. In the latter, the equations of motion of a five-dimensional supergravity action coupled to a matter content of fields will be solved. One may break the conformal invariance even at zero temperature by coupling a scalar field at pure gravity in AdS [48] which is dual to a CFT deformed by a source Λ for a dimension-three operator. This source breaks the scale invariance explicitly and triggers a non-trivial Renormalization Group (RG) flow from ultraviolet (UV) fixed point to an infrared (IR) fixed point. Usually, the scalar potential coupled to the gravity is chosen to mimic the detailed properties of QCD, e.g. [49]. However, the authors in [48] chose their potential by simplicity. The UV fixed point is needed to guarantee that we are in the regime that the holographic duality is best understood and the bulk metric is asymptotically AdS. The IR fixed point guarantees that the solutions are regular in the interior and the zero-temperature solution is smooth in the deep IR. Although turning on a source for a relevant operator is the simplest way to break conformal invariance, this simple model exhibits a rich phenomenology. For example, the relaxation of small-amplitude and homogeneous perturbations are studied by computing the spectrum of quasi-normal modes (QNM) and it is shown that the dominant channel for relaxation in this approximation depends on the value of the ratio T /Λ, with T the temperature of the system. At small T /Λ the system first EoSizes (the evolution of the energy density and the average pressure towards asymptotic values related to one another by the equation of state (EoS)) and then isotropises. While at large T /Λ, the system first isotropises and subsequently EoSizes [48].
One of the interesting properties of QGP is jet quenching processes in which a high energy parton propagates through the medium, interact with medium and consequently lose energy. Analyses the energy loss of these energetic partons as they travel throw QGP may reveal extremely valuable information about the dynamics of the plasma and exhibit distinctive properties such as jet-quenching which can clearly be observed at RHIC [1][2][3][4] and more recently LHC [50][51][52].
In the AdS/CFT correspondence, an external heavy quark with infinite mass can be introduced by adding a fundamental string attached to the flavor brane on the AdS boundary. The endpoint of the string indicates the heavy quark while the string itself can be considered as a gluonic cloud around the quark. The mass of quark is proportional to the length of the string, such that light quark or massless quark is mapped into a string attached to a flavor brane which is extended from the boundary to the horizon [53][54][55][56][57].
A pair of light quark-anti quark is modelled by an initially point-like open string created close to the boundary with endpoints that are free to fly apart [58,59]. The initial conditions of the string imply that the string extends in a direction parallel to the boundary as it falls toward the black hole horizon, so it is called the falling string. The maximum stopping distance of the falling string (the maximum distance which a quark with initial energy E can travel) can be used as a phenomenological guideline to estimate the stopping power of the strongly-coupled plasma. This quantity is calculated numerically for many different sets of string initial conditions in a thermal N = 4 SYM plasma at a temperature T and it is found that the maximum penetration depth scales as x max = C T ( E T √ λ ) n ef f with n ef f = 1/3 [58]. The analogous computation is done in a plasma with non-conformal effect [60] or with anisotropy effect [61] and it is shown that the power of n ef f deviates from its conformal isotropic cousin. The constant of proportionality is also important for phenomenological applications as it determines the overall strength of jet quenching. Here, we analyze the maximum stopping distance of light quark by numerical computation of stopping distance for many different sets of string initial conditions in the non-conformal background introduced in [48].
Another interesting experimental observables associated with quark energy lost in the dense hot QGP is the transport coefficientq which is called "jet quenching". This observable is defined as the average transverse momentum square transferred from the traversing parton, per unit mean free path [62]. As mentioned before, the quark-antiquark pair in the context of AdS/CFT is mapped to a fundamental string with both endpoints attach on the AdS boundary. The string hangs down to the bulk along the radial direction. The jet quenching parameterq is related to the thermal expectation value of the light-like Wilson loop operator [63] which correspond to the trajectory of two endpoints of the string. Several attempts were made to compute this parameter using the AdS/CFT correspondence [63][64][65][66][67][68][69][70][71][72][73][74][75][76][77][78]. In this paper, we also study the jet quenching parameter in the introduced nonconformal background. Our results are in good agreement with experimental expectation. This paper is organised as follows: In section 2 we briefly review the non-conformal background introduced in [48]. In section 3 we discuss the generic falling string solutions and calculate the maximum stopping distance of light quark in the strongly coupled nonconformal plasma. In section 4, we calculate the jet quenching parameter and section 5 is devoted to conclusion and discussion.

Non-conformal holographic model
In this section we review the non-conformal holographic model presented in [48]. The action for the five-dimensional Einstein gravity coupled to a scalar field is where κ 5 is is the five-dimensional Newton constant. The potential considered in this model has the following non-trivial rather simple form This potential possess a maximum at φ = 0 (UV fixed point) and a minimum at φ = φ M (IR fixed point) corresponding to two AdS solutions. L is the radius of the corresponding AdS solution at the UV whereas the radius of the IR AdS is 3) The number of degrees of freedom decreases along the flow (from a UV to an IR) because of the fact that L IR < L.
The vacuum solutions to the Einstein equations can be easily found Here, Λ is an arbitrary constant that breaks conformal invariance. The scalar field is dual to a scalar operator in the dual gauge theory Ø, with different dimension at the UV and IR fixed points. The dimension of this operator at the UV fixed point is ∆ U V = 3 while the operator Ø at the IR fixed point has dimension So, by increasing φ M the dimension of the scalar operator at the IR fixed point decreases and finally reaches ∆ IR = 4 at φ M → ∞. One can determine the vacuum expectation values (VEV) of the stress tensor and the scalar operator by studying the behavior of the metric and the scalar field near the boundary. Using the new variable u as u = L e −r/L , the metric can be rewritten as [79] The expectation values of the field theory operators can be readily determined by expanding the metric and the scalar field in powers of u in the u → 0 limit as follows [48] which imply the Ward identity for the trace of the stress tensor Since the scalar field is a monotonic function of r, we can use the φ as a coordinate and consider the following ansatz for black brane solutions of the action in the Eddington-Finkelstein form where h(φ) is the blackening factor which is zero at φ = φ H . The boundary of space is at φ = 0 and φ H is the location of the horizon. With this ansatz, the equations of motion can be written as where primes denote d/dφ. It is shown that these four equations can be combined into a master equation of the form [48] G (φ) where G(φ) = d dφ A(φ) is assumed to be a smooth generating function which can be integrated to obtain A(φ) as Integrating the first three equations of Eq. Eq. (2.13) leads to Using the fact that h has a simple zero at the horizon, the Eq. (2.13c) and Eq. (2.13d) give rise to which leads to the solution of G(φ) at the horizon as So, the near horizon solution of G can be read as On the other hand, near boundary (φ → 0) solution of G turns out to be where ∆ is the dual operator scaling dimension which is 3 for the choice of potential Eq. (2.2) as mentioned above. Substituting the Eq. (2.20) into the Eq. (2.15) and requiring that A must be finite at the boundary, the function A is obtained Comparing this equation with Eq. (2.15) at the limit of φ 0 → φ H lead to and finally A can be expressed as There is a residual gauge freedom in the ansatz of metric dr = ± e B dφ. As φ increases from zero to φ H , the radial coordinate r decreases from ∞ to a finite value. Therefore the minus sign is meaningful. The leading asymptotic behavior of metric functions at large r is Therefore, the asymptotic behavior of generating function at large r becomes (2.25) By considering the lower limit of the integral in the Eq. (2.16a) to be φ H and manipulating the equation, one can find the following relation and finally we get (2.28) In the last step, one can obtain the function h(φ) from Eq. (2.13c) and Eq. (2.13d) as follows The thermodynamical quantities of the background are expressed as Substituting the metric functions at the horizon leads to the following form for the temperature and entropy of the plasma The non-conformal behavior of the thermodynamics of the dual theory can be identified by studying the entropy, the speed of sound, the ratio of bulk viscosity to shear viscosity and the Violation of Buchels bound. In Fig. 1a, we plot the ratio of non-conformal to the conformal entropy in terms of temperature. At high and low temperature, the entropy of the theory is the same as the conformal theory and scale as T 3 .
The square of the speed of sound can be obtained from the inverse of the logarithmic derivative of the entropy as 1 The deviation of the speed of sound from its conformal value (c s = 1/ √ 3) is depicted in the Fig. 1b as a function of temperature. At low and high temperatures, this quantity reaches its conformal value while the deviation gets larger by increasing φ M at intermediate temperatures. The ratio of shear viscosity η to entropy has a universal value of η/s = 1/4π for all theories with a two-derivative gravity dual [6]. The bulk viscosity, ζ can be determined from the dependency of the entropy on the value of the scalar field at the horizon [80] ζ which is zero in any conformal field theory. This ratio has been shown in the Fig. 1c as a function of temperature for different values of φ M . The behavior of this quantity as a function of temperature is similar to that of the speed of sound. At low temperatures, the violation of this ratio from Buchels bound is shown in the Fig. 1d.

Light quark maximum stopping distance
We are interested to study the propagation of energetic excitations which resemble quark jets through the strongly coupled non-conformal plasma. These excitations are considered as open string configurations with high energy, moving through the bulk of AdS space. The maximum stopping distance of these strings, the quantity which is not sensitive to the precise initial conditions, can be used as a phenomenological guideline to estimate the stopping power of the strongly-coupled plasma.
In this section, we first briefly review the string configuration dual to a light quark jet in the non-conformal background and then numerically calculate the maximum stopping distance which a quark with initial energy E can travel.

String configuration
We use the black brane solution of Eq. (2.12) in Poincare coordinates to study the nonconformality effect on the dynamics of light quark jet moving in a strongly coupled plasma. According to the gauge/gravity duality, quarks moving through the medium are dual to open strings moving in the 10-d gravitational background. D-branes are added to the background in the probe limit to describe the fundamental matter of the field theory such that do not back-react on the background metric Eq. (2.12). These branes fill the 4D Minkowski space and extend along the radial coordinate from the boundary at φ = 0 down to the location of the horizon at φ = φ H . According to the AdS/CFT dictionary, quark-antiquark pairs in the field theory side are dual to open strings whose endpoints are attached to the D-brane. In this paper, we consider a back-to-back jet pair created in the medium which can be regarded as a string created at a point, extends in space-time and falls toward the horizon. The profile evolution of a typical string for various times has been shown in Fig. 2. In this setup, the two endpoints of the string move away from each other and the total spatial momentum of the string vanishes. The string embedding function is X µ (τ, σ) → (t(τ, σ), x(τ, σ), 0, 0, φ(τ, σ)) and for an open string created in a point φ = φ c at a time t c , the profile is given by where t c is the time for which the string becomes an extended object from a point and σ ∈ [0, π]. The Nambu-Goto action becomes singular for falling string at late times, therefore we use the Polyakov action in which additional degrees of freedom are being involved by a nontrivial worldsheet metric η ab and the equations of motion are well-behaved everywhere on the worldsheet [53,58,81]. The Polyakov action for the string has the following form Canonical momentum densities associated with the string can be obtained by varying the action with respect to the derivatives of the embedding functions, and the equations of motion are being obtained by variation of the Polyakov action with respect to the embedding functions X µ as We choose the following worldsheet metric where Σ is called the stretching function and generally can be an arbitrary function of worldsheet embedding functions. We found that choosing the stretching function of the form would cancel the singularities and equations of motion remain well behaved everywhere on the worldsheet. We choose values of a and b in the range of 1 to 5 depends on the string initial conditions. The constraint equation is obtained by variating the Polyakov action, Eq. (3.2) with respect to η ab which for a string with point-like initial condition reduces tȯ Combining the open string boundary condition and constraint equation, the initial profile must satisfy the following conditioṅ The below initial condition (IC) obeys all necessary conditions and assures the physical requirements for the non-conformal backgrounḋ where φ c and A are free parameters related to the energy and momentum of the dual quark in the field theory side. The IC are chosen such that the string is long-lived and yields stable numerical solutions. Also, most of the string's energy and momentum are concentrated near its endpoints. These IC yield a symmetric string profile about x = 0 at all times, becauseẋ(0, σ) is antisymmetric about σ = π/2 whileφ(0, σ) is symmetric. We performed the explicit check of any numerical solution, and our solutions respected the equation of constraint for all τ . Now, we proceed to finding the light quark energy and momentum. As the nonconformal background geometry G µν depends only on φ, hence for µ corresponding to (t, x) one can write ∇ a Π a µ = 0 . (3.9) Therefore, the momentum densities Π a µ are conserved Noether currents on the related worldsheet. The total energy of the falling string in the non-conformal geometry of Eq. (2.12) at any time is constant and equal to the initial energy of the string and the total energy of the light quark is obtained from The total momentum of the string is conserved in x-direction and can be find from the string IC as Here, we have considered that the total momentum of the string in x-direction is completely symmetric about σ = π/2 and equals to zero . Therefore, the quark momentum equals to the anti-quark momentum with a minus sign. Also, the string energy and momentum are determined by two parameters A and φ c in the string initial condition. These parameters are translated into the virtuality of the jet created in the non-conformal QGP which is obtained from The solution of Eq. (3.4) for a typical string with E = 100 GeV in the non-conformal background Eq. (2.12) with φ M = 3 is plotted in Fig. 2. The temperature of the plasma is about T = 0.265 GeV . The dynamics of the two half of the string are exactly the same, one moving in the opposite direction of the another one. As we expected, the two endpoints of the string move away from each other as the string extends along the (x, φ) direction and falls toward the horizon.

Maximum stopping distance
Numerical studies, Fig. 3 show that the distance which a quark with a specific energy traversed in the plasma before thermalizing highly depends on the chosen initial condition of the dual string, i.e. the initial distance from the boundary. Since there is no known map between the string IC and the actual field theory quantities, we aim to study a quantity which is insensitive to these IC in the string side. The stopping distance which is called the thermalization distance, x therm is defined as the length along the x-direction from the point of production of the original point-like string to the point at which the end of the string falls through the black hole horizon. On the field theory side of the duality, x therm corresponds to the length of the plasma traversed before the jet becomes completely thermalized (i.e. indistinguishable from the plasma). We are looking for the maximum distance which a  Figure 3: (Color online) The maximum stopping distance of the jet moving in the plasmas versus to its initial distance from the boundary located at φ = 0. The non-conformal plasma has the same temperature as the conformal plasma, T = 350 MeV. In all cases, the quark has 100 GeV energy but created at a different radial distance from the boundary. The dashed black line is the stopping distance in the AdS-Sch background, while the blue, red and green lines are stopping distance in the non-conformal plasmas with φ M = 1, 3 and 10, respectively. The stopping distance increases by increasing φ M in the medium ranges of φ c where it is maximum. quark with energy E can travel before thermalizing independent of the exact IC of the string. Although this quantity is not enough to calculate observables such as R AA or v 2 in general, it can be used as a phenomenological guideline to estimate the stopping power of the strongly-coupled plasma.
In order to compare the results of stopping distance of jet in non-conformal plasma with the results of conformal AdS-Sch, we can either fix the temperature or the entropy of the non-conformal plasma with the AdS-Sch ones. In Fig. 4 we plot the rescaled entropy of the non-conformal theory with different φ M compared to the conformal entropy in terms of temperature. The intersection of black dashed lines which shows the entropy of the AdS-Sch metric with T = 265 MeV is our reference point. This figure indicates that if we fix the temperature and increase φ M (moving along the vertical dashed black line), the entropy of the non-conformal theory decreases. On the other hand, increasing φ M on the line of fixed entropy (moving along the horizontal dashed black line) give rises to the non-conformal theories with higher temperature.
In Fig. 5, and Fig. 6 we numerically compute the stopping distance for many different sets of string initial conditions in backgrounds with different non-conformal parameter. In these figures, each point shows the stopping distance of a string with a specific initial condition in terms of its initial energy. Black dots show the strings in the AdS-Sch metric, while the blue triangles, red squares, and green stars correspond to strings in the nonconformal backgrounds with φ M = 1, φ M = 3, and φ M = 10, respectively. In the Fig. 5, all plasmas have the same temperature while in the Fig. 6 the entropy of the backgrounds is fixed. The interesting point is that the stopping distance increases with increasing the deviation from conformality when the temperature of the background is fixed while it decreases if the entropy of the background is fixed. It is expected, since as explained in Fig. 4, by fixing the entropy and increasing the φ M the plasma becomes hotter which leads to a larger jet suppression. As clearly shown in these figures, the dynamics of the string is a multi-variable quantity depends on both string IC and background thermodynamics. But, all data points for each set of data fall below the solid lines. These lines display the maximum stopping distance of the light quark in the corresponding medium, i. e. the maximum distance that a quark with energy E can travel in the medium before thermalization. It is shown that this distance in the conformal background scales like E 1/3 for enough large energies [58]. In fact, numerical results show that the maximum stopping distance for a given energy depends on the energy of the quark and temperature of the plasma as follows

AdS-Sch
which C and n ef f is estimated to 0.526 and 1/3 in the AdS-Sch case, respectively [58]. We explore this relation with our numerical results and estimate the values of C and n ef f in Figure 5: (Color online) The stopping distance of the jet in terms of its energy moving in the non-conformal backgrounds with different φ M . All of the non-conformal backgrounds have the same temperature as the AdS-Sch temperature sets to be 265 MeV. For each set of data, all data points fall below the solid line which determines the maximum stopping distance of the quark in terms of its energy.
the case of non-conformal plasma. We fit our data in the AdS-Sch with Eq. (3.14) by fixing C = 0.526. The estimated n ef f is less than 1/3 for the available energies at RHIC or LHC but it approaches this value at very high energies. In Fig. 7, the n ef f is calculated by fitting the numerical results to the Eq. (3.14) by considering the C = 0.45 for all non-conformal metrics and C = 0.526 for the AdS-Sch metric. In fact, these constants give rise to a more precise fit to the numerical results. Again, in Fig. 7a all backgrounds have the same T and in Fig. 7b they have the same entropy. In general, increasing the deviation from conformal invariance increases the stopping distance of jets of light quark.

Jet Quenching
To calculate the jet quenching parameter for the non-conformal theory, we use the thermal expectation value of a close light-like Wilson loop [82] W where L is the transverse distance (conjugate to the transverse momentum of the radiated gluons) and L − is the light-cone distance ( conjugate to partons with relativistic velocities). Figure 6: (Color online) The stopping distance of the jet in terms of its energy moving in the non-conformal backgrounds with different φ M . All of the non-conformal backgrounds have the same entropy as the AdS-Sch entropy. For each set of data, all data points fall below the solid line which determines the maximum stopping distance of the quark in terms of its energy. where S I is the normalized action after subtracting the self energy of the qq pair from the Nambu-Goto action of the string worldsheet. As in the large N c limit, Tr (Adj.) = Tr 2 (F und.) , one can writeq In the light-cone coordinates, x ± = (t±x 1 )/ √ 2, the non-conformal background of Eq. (2.12) becomes The string parametrization is x µ (τ, σ) where σ α = (τ, σ) is the worldsheet coordinates. The Nambu-Goto action of the string is where g αβ is the induced metric on the string worldsheet and (τ, σ) are set to be (x − , x 2 ). L is considered to be the contour length along x 2 -direction and L − to be its length along τ -direction. The boundary conditions are φ(± L 2 ) = 0 and x 3 (σ) and x + (σ) coordinates are constant. Then the action of Eq. (4.5) reads where prime denotes the derivative with respect to σ. Since the Lagrangian density is time independent, the Hamiltonian of the system is constant and one can obtain the following equation for φ Integrating equation Eq. (4.8) leads to where (4.10) Here, we have used the fact that for small length L, the constant Π φ is small and its higher order terms are negligible. Substituting Eq. (4.8) into Eq. (4.6) and expand for small Π φ yields This action diverges and we should subtract the self energy of two disconnected strings whose worldsheets are located at x 2 = ± L 2 and extended from φ = 0 to φ = φ H The normalized action is then written as Inserting Eq. (4.13) into Eq. (4.3) leads to the following expression for the jet quenching parameter of the non-conformal theorŷ where a 0 is the numerical integral of Eq. (4.10). Here, we have used the relation between the string tension and the tHooft coupling R 2 α = √ λ. For N = 4 supersymmetric Yang-Mills theory in the large N c and large λ limits, the jet quenching parameter is find to be [88] To understand how non-conformality affects the jet quenching parameter, we calculate Eq. (4.14) for different values of φ M numerically. Here, we have taken the 't Hooft coupling λ = 5.5 for numerical estimates. The ratio of jet quenching in non-conformal background to its conformal value in terms of temperature is shown in Fig. 8a for φ M = 1 (blue curve), φ M = 3 (red curve) and φ M = 10 (green curve). The plot shows that for each value of φ M , the ratio ofq N C /q SY M starts from a finite value and approaches to 1 at higher temperature where conformality is dominant. This behavior reflects the fact that deviations from conformality are magnificent at lower temperatures and suppress at higher temperatures. Also, increasing the values of φ M , decreases the values of jet quenching (and hence the ratio). In Fig. 8b, the temperature dependency of the jet quenching parameter is shown for φ M = 1 (blue curve), φ M = 3 (red curve), φ M = 10 (green curve) as well as N = 4 SYM theory (dashed curve). Two black circles with error bars are the corresponding absolute values forq for a 10 GeV quark jet in the most central Au-Au collisions at RHIC with highest temperature T = 0.37GeV and Pb-Pb collisions at LHC with highest temperature T = 0.47GeV [89]. Both temperatures are rescaled due to the fact that the number of degrees of freedom in the non-conformal theory is more than those of 3 favor  QCD and one can use T N C ≈ T SY M = 3 −1/3 T QCD [89].
In Table 1, the numerical values of the jet quenching parameter is shown for N = 4 SYM theory, non-conformal background (φ M = 1, φ M = 3 and φ M = 10) and experimental data from RHIC and LHC.q values are in the units of GeV 2 /f m. It could be found that the jet quenching results from non-conformal background are in good agreement with the experimental data. It could be found that at T = 0.37GeV varying φ M from 1 to 3, the decrease inq is ∼ 26% while varying φ M from 3 to 10, the decrease inq is ∼ 9%. On the other hand, at T = 0.47GeV varying φ M from 1 to 3, the decrease inq is ∼ 13% while varying φ M from 3 to 10, the decrease inq is ∼ 2%. Also, our results imply an increase of the jet-medium interaction at lower temperature where non-conformality effects are magnificent.
Holographic models such as [65,73] indicate that increasing nonconformality leads to decreasing the value of jet quenching parameter which is in consistency with our results from jet quenching parameter and also with results of light quark stopping distance obtained in the previous section. In general, increasing the deviation from conformal invariance decreases the capability of medium to quench the jets of quarks. However, in the case of holographic QCD models with phase transition, this parameter exhibits a peak around the QCD phase transition region due to the fact that the system degrees of freedom change rapidly from hadronic gas to QGP [90,91].

Summary
The quark-gluon plasma produced in heavy ion collisions is a strongly coupled, nonconformal plasma as indicated by lattice data [38] and hydrodynamics calculations [39][40][41][42]. The gauge/gravity duality gives us the opportunity of the study of hot plasma in a strongly coupled regime. Suppression of high energy partons produced at heavy ion collision with high transverse momentum is one of the interesting properties of this strongly-coupled plasma that can be investigated using the AdS/CFT correspondence. Although, previous studies show that suppression of jets with sensible energies in N = 4 SYM theory is much larger than its experimental expectation [59,61].
In this paper, we investigated the suppression of jet in the strongly-coupled nonconformal plasma in the framework of AdS/CFT correspondence. We considered a holographic five-dimensional model consisting of Einstein gravity coupled to a scalar field with a non-trivial potential which turns out a dual four-dimensional non-conformal gauge theory which exhibits a renormalization group flow between two different fixed points (located at UV and IR) at zero temperature [48]. The parameter φ M indicated the deviation from conformality as shown in Fig. 1.
We considered a point-like initial condition string created close to the boundary with endpoints that are free to fly apart. The initial conditions of the string are chosen such that the string extends in a direction parallel to the boundary as it falls toward the black hole horizon (falling string). The equations of motion for the string moving the non-conformal background have been solved numerically and plotted in Fig. 2 for a typical string moving in the background with φ M = 3. Since the dynamics of the string depends on the string initial conditions (see Fig. 3, Fig. 5 and Fig. 6 for clarifications), we studied the maximum stopping distance which is insensitive to these IC in the string side. Although, this quantity is not good enough to compute the observables like nuclear modification factor, still is a good quantity to estimate the stopping power of the plasma.
Our results show that the thermalization distance in non-conformal backgrounds compare to the N = 4 SYM theory depends on the fact that theories have the same temperature, Fig. 5 or the same entropy, Fig. 6. Actually, the Fig. 4 demonstrates that increasing non-conformality while keeping the entropy fixed leads to hotter plasma and then larger jet suppression. In order to avoid that, we fitted the maximum stopping distance to the well-known relation presented in Eq. (3.14) and used the n ef f as the stopping power of the plasma. Our numerical results, presented in Fig. 7 indicated that the thermalization distance increases by increasing the non-conformality of the plasma.
In the last section, we calculated the jet quenching parameter for the non-conformal geometry of section 2. For different values of φ M , our results are in a good agreement with experimental data. Deviations from conformality are magnificent at lower temperatures and suppress at higher temperatures where non-conformal results approach the corresponding value of the conformal background. As non-conformality increases, the value of jet quenching parameter decreases and this behavior is in consistency with experiments and also with results of light quark stopping distance obtained in section 3.2.