arxiv: v1 [nucl-th] 20 Apr 2018

Transcription

1 Calculatin hard probe radiative enery loss beyond soft-luon approimation: how valid is the approimation? arxiv: v [nucl-th] Apr 8 Bojana Blaojevic, Madalena Djordjevic, and Marko Djordjevic Institute of Physics Belrade, University of Belrade, Belrade, Serbia Faculty of Bioloy, Institute of Physioloy and Biochemistry, University of Belrade, Belrade, Serbia Dated: April 3, 8 The soft-luon approimation, which implies that radiated luon carries away a small fraction of initial parton s enery, is a commonly used assumption in calculatin radiative enery loss of hih momentum partons traversin QGP created at RHIC and HC. While soft-luon approimation is convenient, different theoretical approaches reported sinificant radiative enery loss of hih p partons, thereby questionin its validity. To address this issue, we relaed the soft-luon approimation within DGV formalism. The obtained analytical epressions are quite distinct compared to the soft-luon case. However, numerical results for the first order in opacity fractional enery loss lead to small differences in predictions for the two cases. The difference in the predicted number of radiated luons is also small. Moreover, the effect on these two variables has an opposite sin, which when combined results in almost overlappin suppression predictions. Therefore, our results imply that, contrary to the commonly held doubts, the soft-luon approimation in practice works surprisinly well in DGV formalism. Finally, we also discuss eneralizin this relaation in the dynamical QCD medium, which suests a more eneral applicability of the conclusions obtained here. PACS numbers:.38.mh; 4.85.p; q I. INTRODUCTION One of the main assumptions in the radiative enery loss calculations of eneretic parton in the further tet referred to as jet traversin the QGP medium, is the soft-luon approimation which assumes that radiated luon carries away a small portion of initial jet enery, i.e. = ω/e, where E is the enery of initial jet and ω is the radiated luon enery. Such assumption was widely used in various enery loss models: i in multiple soft scatterin based ASW model [ 3]; ii BDMPS [4, 5] and BDMPS-Z [6, 7]; iii in opacity epansion based GV model [8, 9] and iv in multi-luon evolution based HT approach [, ], etc. These various enery loss models predict a sinificant medium induced radiative enery loss, questionin the validity of the soft-luon approimation. To address this issue, a finite or lare limit was introduced in some of these models [, 3] or their etensions [4]. However, introduction of finite lead to different conclusions on the importance of relain the soft-luon approimation, which was assessed from relatively small [4], but noticeable, to moderately lare [3]. The soft-luon approimation was also used in the development of our dynamical enery loss formalism [5 7], specifically in its radiative enery loss component. This formalism was comprehensively tested aainst anular averaed nuclear modification factor R AA [8, 9] data, where we obtained robust areement for wide rane of probes [7, ], centralities [] and beam eneries [, ], includin clear predictions for future eperiments [3, 4]. This miht stronly suest that our enery loss formalism can well eplain the jet-medium interactions in QGP, makin this formalism suitable for the tomoraphy of QCD medium. However, the soft-luon approimation obviously breaks-down for: i intermediate momentum ranes 5 < p < GeV where the eperimental data are most abundant and with the smallest error-bars, and ii luon enery loss, since due to the color factor of 9/4 luons lose sinificantly more enery compared to quark jets, therefore questionin the reliability of our formalism in such cases. Due to this, and for precise predictions, it became necessary to rela the soft-luon approimation, and consequently test its validity in dynamical enery loss formalism. This paper presents our first step toward this oal. Since the dynamical enery loss is computationally very demandin, we will, in this study, start with relain this approimation on its simpler predecessor, i.e. DGV [5] formalism. Within this, we will concentrate on luon jets, since, due to their color factor, the soft-luon approimation has the larest impact for this type of partons. For the luon jets, we perform the radiative enery loss calculation, to the first order in the number of scatterin centers opacity, where we consider that the radiation of one luon is induced by one collisional interaction with the medium. Our calculation is done within the pqcd approach for a finite size, optically thin QCD medium and since it is technically demandin it will be divided in several steps: i First, the calculation will be done in the simplest case of massless luons in the system of static scatterin centers [6] within GV, ii Then it will be etended towards the luons with the effective mass [7], which presents epansion of DGV [5] toward larer loss of jet enery via radiated luon, and iii Finally, we will discuss the impact of finite on the radiative enery loss, when

2 dynamical medium [7] i.e. a recoil with the medium constituents is accounted. In that manner we will assess the validity of the softluon assumption for luon jets, and this will also provide an insiht into whether or not a finite has to be implemented in quark-jet radiative enery loss calculations within our formalism. Namely, if the relaation of the soft-luon approimation only slihtly modifies luonjet radiative enery loss, then even smaller modification would be epected in quark-jet case, thus makin this relaation redundant. Otherwise, if the effect of a finite appears to be a sinificant in luon-jet case, then the relaation in quark-jet case may also be required, and would represent an important future task. Secondly, as stated above, the relaation of the softluon approimation is needed in order to etend the applicability of our model [7] towards intermediate momentum reion. Thus, the another benefit of this relaation would be to etend the p rane in which our predictions are valid. The sections are oranized as follows: In section II, we provide the theoretical framework. In section III, we outline the computation of the zeroth order in opacity luonjet radiative enery loss in static QCD medium, beyond soft-luon approimation, in the cases of both massless and massive luons. For the results from [9, 5] are reproduced. Section IV contains concise description of relain the soft-luon approimation in calculatin the first order in opacity radiative enery loss for massless luon jet in static QCD medium. In a limit of very small result from [9] is recovered. In section V we eplain the computation of the first order in opacity luon-jet enery loss in static QCD medium, with effective luon mass [7] included, and beyond soft-luon approimation. This presents an etension of the calculations from [5] toward finite, so that results from [5] can be recovered in limit. The detailed calculations correspondin to sections III - V are presented in the Appendices C-J. In section VI we outline the numerical estimates based on our beyond soft-luon calculations for luon jet and the comparison with our previous results from [5], i.e. the results with soft-luon approimation. Particularly, we investiate the effect of finite on luon-jet fractional radiative enery loss, number of radiated luons, differential radiative enery loss, sinle luon radiation spectrum and luon suppression [8]. Conclusions and outlook are presented in section VII. II. THEORETICA FRAMEWORK In this work, we concentrate on relain soft-luon approimation in calculatin the first order in opacity radiative enery loss of hih p eikonal luon jets within GV DGV [5] formalism. That is, we assume that hih p luon jet is produced inside a thin finite QGP medium at some initial point t, z,, and that the medium is composed of static scatterin centers [6]. Therefore, we model the interactions in QGP assumin a static Debye colored-screened Yukawa potential, whose Fourier and color structure acquires the followin form [9, 5, 6]: V n = V q n e iqnn = πδq nv q n e i qn n T an R T an n, v q n = 4πα s q n µ, where n denotes time-space coordinate of the n th scatterin center, µ is Debye screenin mass, α s = s/4π is stron couplin constant, while T an R and T an n denote the enerators in SUN c = 3 color representation of luon jet and taret scatterin center, respectively. For consistency with [9, 5], we use the same notation for 4D vectors e.. momenta, which is described in detail in Appendi A and proceed throuhout usin iht-cone coordinates. The same Appendi contains alebra manipulation and identities for SUN c enerators, as well as the Feynman rules, used in these calculations. The approimations that we assume throuhout the paper are stated in Appendi B. The small transverse momentum transfer elastic cross section for interaction between luon jet and taret parton in GW approach [8, 6] is iven by: dσ el d q = C GC T d G v, q π, 3 where q corresponds to transverse momentum of echaned luon, C G represents Casimir operator in adjoint representation G of luons SUN c = 3 with dimension d G = 8, whereas C T denotes Casimir operator in taret T representation. Since this formalism assumes optically thin plasma, the final results are epanded in powers of opacity, which is defined as the mean number of collisions in the medium: /λ = Nσ el /A [9], where is the thickness of the QCD medium, λ is a mean free path, while N denotes the number of scatterers tarets in transverse area A. Note that, we restrict our calculations to the first order in opacity, which is shown to be the dominant term [9, 3]. III. ZEROTH ORDER RADIATIVE ENERGY OSS To radually introduce technically involvin beyond soft-luon calculations, we first concentrate on massless luons traversin static QCD medium. We start with M Feynman diaram, which corresponds to the source J that produces off-shell luon with momentum p k, that further, without interactions with

3 3 QCD medium, radiates on-shell luon with momentum k and emeres with momentum p. We will further refer to these two outoin luons as the radiated k and the final p luon. Note that, both in this and consecutive sections that involve interactions with one and two scatterin centers, we consistently assume that initial jet propaates alon the lonitudinal z ais. The detailed calculation of M for finite in massless case is presented in Appendi C, with all assumptions listed in Appendi B. We also assume that luons are transversely polarized particles and althouh we work in covariant aue, we can choose any polarization vector for the eternal onshell luons [4], so in accordance with [9, 4, 5] we choose n µ = [,, ] i.e. ɛk k =, ɛk n = and ɛp p =, ɛp n =. ikewise, we assume that the source has also the same polarization as real luons [4] i.e. ɛp k p k =, ɛp k n =. Thus, for massless luon s momenta we have: p k = [E, E, ], k = [E k, E, k], p = [ E p,, p], 4 E where E = p k p z k z, E = p k p z k z and due to 4-momentum conservation: The polarization vectors read: p k =. 5 ɛ i k = [, ɛ i k E, ɛ ɛ i p i], ɛ i p = [, E, ɛ i], ɛ i p k = [,, ɛ i ], 6 where i =,, and we also make use of Eq. 5. So, the amplitude that luon jet, produced at inside QCD medium, radiates a luon of color c without final state interactions reads: M = J a p ke ipk i s ɛ k k T c da. 7 The radiation spectrum is obtained when Eq. 7 is substituted in: d 3 N d 3 N J Tr M d 3 p d 3 k π 3 p π 3 ω, 8 where ω = k, and where d 3 N J reads: d 3 N J = d G Jp k d 3 p J π 3. 9 E J Here E J = E = p k and p J denotes enery and 3D momentum of the initial luon jet, respectively. The jet part can be decoupled by usin the equality: d 3 p d 3 k π 3 p π 3 ω = d3 p J d k π 3 E J π 3, which is obtained by substitutin p z, k z p J z, E. Finally, enery spectrum acquires the form: d 3 N dk = α s π C G k, which recovers the well-known Altarrelli-Parisi [3] result. We now briefly concentrate on eneratin result in finite temperature QCD medium, since in [7], it was shown that luons in finite temperature QGP can be approimated as a massive transverse plasmons with mass m = µ/, where µ is the Debye mass. In this case, M amplitude becomes: M = J a p ke ipk i s ɛ k k m T c da, leadin to: IV. d 3 N dk = α s π C G k k m. 3 FIRST ORDER RADIATIVE ENERGY OSS IN MASSESS CASE In accordance with [5], we compute the first order in opacity radiative enery loss of luon jet for finite startin from the epression: d 3 N d 3 N J = Tr M Re Tr M M d T d T d 3 p π 3 p d 3 k π 3 ω, 4 where M corresponds to the diaram without final state interactions with QCD medium, introduced in previous

4 4 section, M is the sum of all diarams with one scatterin center, M is the sum of all diarams with two scatterin centers in the contact-limit case, while d T denotes the dimension of the taret color representation for pure luon medium d T = 8. In obtainin the epression for differential enery loss, we aain incorporate Eqs. 9, in Eq. 4. The assumption that initial jet propaates alon z- ais, takes the followin form in the two cases stated below:. One interaction with QCD medium M : p k q = [E q z, E q z, ], 5 where pk q corresponds to the initial jet, while k and p retain the same epressions as in Eq. 4, with the distinction that now p k, since due to 4-momentum conservation, the followin relation holds: q = p k; 6 The rest of the notation is the same as in Eq. 4.. Two interactions with QCD medium M : p k q q = [E q z q z, E q z q z, ], 7 where p k q q corresponds to the initial jet and q i = [q iz, q iz, q i ] to echaned luons, i =, with qi =, while p, k retain the same epressions as in Eq. 4. Also, due to 4-momentum conservation, the followin relation between luon transverse momenta holds: p k = q q, 8 which in the contact-limit case when q q = reduces to p k =. Note that Eq. 6 has to be satisfied for M diarams in order to claim that initial jet propaates alon z-ais, i.e. for M diarams p k is different from. This is an important distinction between the calculations presented in our study, and the calculations done within SCET formalism see e.. [4], where p k = was used in calculation of both M and M diarams, thouh the assumption that initial jet propaates alon z-ais was used in that study as well. The transverse polarization vectors ɛ i k and ɛ i p for both: M and M amplitudes are iven by the same epression as in the previous section with an addition that in M case: p k, as discussed above, while ɛ for initial jets consistently has the same form as in Eq. 6, i.e. ɛ i p k q = [,, ɛ i ] for M amplitudes, and ɛ i p k q q = [,, ɛ i ] for M amplitudes. The detailed calculation of the remainin Feynman diarams, under the approimations stated in Appendi B, contributin to the first order in opacity radiative enery loss, is iven in Appendices D-H, whereas thorouh derivation of the sinle luon radiation spectrum beyond soft-luon approimation in massless case is iven in Appendi I and reads enery loss epression can be straihtforwardly etracted by usin de / ωdn / EdN /: dn = C Gα s π { q λ π 4 E k q k q 4 µ q µ k k k q 4 E k 4 k q dk k k q k k q k q k q k q } k, 9 where we assumed a simple eponential distribution e z of lonitudinal distance between the luon-jet production site and taret rescatterin site, emerin as 4 E in the denominators of the interand. Note that, Eq. 9 reduces to massless case of Eq. from [5] in the limit, as epected. It is straihtforward to show that our result is symmetric under the echane of radiated k and final p luon, as epected beyond soft-luon approimation, and due to inability to distinuish between these two identical luons.

5 5 V. GUON RADIATIVE ENERGY OSS IN FINITE TEMPERATURE QCD MEDIUM Net, we note that in ultra-relativistic heavy ion collisions, finite temperature QCD medium is created, that modifies the luon self eneries, and can consequently sinificantly influence the radiative enery loss results. It is therefore essential to include finite temperature effects in luon radiative enery loss calculations beyond soft-luon approimation, which is the main oal of this section. To address this issue, we note that in [7], it was shown that luons can be approimated as massive transverse plasmons with effective mass m for luons with the hard momenta k T equal to its asymptotic value. The assumption of initial jet propaatin alon z-ais, for massive case, leads to the followin form of momenta, in the three cases stated below:. No interaction with QCD medium M : p k = [E, E, ], k = [E, k m E, k], p = [ E p m,, p], E where Eq. 5 holds;. One interaction with QCD medium M : k and p retain the same epressions as in Eq., with addition that as in the previous section Eq. 6 holds due to conservation of 4-momentum, while initial jet has the momentum of a same form as in Eq Two interactions with QCD medium M : p, k have the same epressions as in Eq.. Also, due to 4-momentum conservation Eq. 8 holds and in the contact-limit case reduces to p k =, while initial jet momentum retains the same form as in Eq. 7. The transverse polarization vectors remain the same as in the massless case. We retain all approimation from the previous section, which are reviewed in Appendi B, and recalculate the same diarams from Appendices C-H, also beyond soft-luon approimation. The overview of all intermediate results is contained in Appendi J. Thus, Eq. 9 in the massive case acquires more comple form iven by: dn q λ π = C Gα s µ π q µ dk { k q χ k q 4 E k q χ k q χ k k q k χ k χ k 4 E k χ k χ k k q k q χ k q k q k q χ k q k q χ k } k χ, dn where χ = m. It can easily be verified that, in the soft-luon limit, we recover Eq. from [5] note that for luon jet M m, so that the term M from [5] should be nelected, and that in the massless limit Eq. reduces to Eq. 9. To our knowlede, this result presents the first introduction of effective luon mass beyond-soft-luonapproimation radiative enery loss. Additionally, we aain verified that sinle luon radiation spectrum is symmetric to substitution of p and k luons, as necessary see the previous section and Appendi J. Furthermore, note that the analytical form of Eq. is quite different from the correspondin epression with the soft-luon approimation Eq. from [5]. In the net section, we will evaluate the etent of numerical differences to which these two different analytical epressions lead. In particular, we are interested in what is the effect of finite on luon fractional radiative enery loss E E, number of radiated luons N and on the suppression R AA. We accordinly note that de ω E from which we can further straihtforwardly nu- dn merically evaluate E E ated luons N. VI., as well as the number of radi- NUMERICA RESUTS We net assess how the relaation of soft-luon approimation modifies luon-jet enery loss to the st order in opacity. We consequently compare the predictions based on the results derived in this paper, with the one obtained in the soft-luon limit from [5] applied to luons - the comparison is done for luons with effective mass m = µ/, where µ = 4πα s n f /6T and n f = 3 is the number of the effective liht-quark flavor. For all fiures, we assume constant stron couplin [3]

6 6.8 3 E E N p GeV FIG. : The effect of relain the soft-luon approimation on interated variables to the st order in opacity of DGV formalism, as a function of p. The top left panel compares luon s fractional radiative enery loss without the solid curve and with the dashed curve soft-luon approimation. The top riht panel represents relative chane of the radiative enery loss when the soft-luon approimation is relaed with respect to the soft-luon limit. The bottom left panel compares number of radiated luons without the solid curve and with the dashed curve soft-luon approimation, whereas the bottom riht panel provides a percentae of radiated luon number chane when soft-luon approimation is relaed. α s = s 4π =.3, and use = 5 fm, λ = fm, T = 3 MeV, to mimic standard HC conditions. The top left panel of Fi. shows comparison of the fractional radiative enery loss E E, for calculations beyond the soft-luon approimation, and with the softluon approimation, as a function of initial jet transverse momentum p. More specifically, the curve correspondin to beyond soft-luon approimation bs case is obtained from Eq. multiplied by E and interated over, while the curve correspondin to soft-luon approimation s case is obtained by numerically interatin Eq. from [5]. These two curves almost overlap, even convere towards one another at hiher p. Note that, the upper limit of interation is equal to / instead of, in order to avoid double countin. The upper limits of interation for k and q, determined kinematically, are E and 4ET, respectively [5]. The bottom left panel of Fi. presents comparison of number of radiated luons in bs and s cases. These two curves also nearly overlap, with a sliht disareement at hiher p. Quantitative assessment of relain the soft-luon approimation on these two variables can be observed from the two riht panels of Fi.. We see that finite values of slihtly increase fractional radiative enery loss by maimum of 3% up to p GeV compared to s case. Afterwards, the difference between bs and s E E steeply decreases towards %. Additionally, finite also decreases number of radiated luons for a small amount up to 5% compared to s case for very low transverse momenta. Further the relative difference reaches a peak of % also at p GeV, and for hiher transverse momenta remains nearly constant somewhat below %. Consequently, the overall conclusion from Fi. is that the effect on both variables is small and with opposite sins. The effect of finite value is further assessed on the differential radiative enery loss de, and on sinle dn luon radiation spectrum and it s relative chane. These effects are shown as a function of on Fi., for

7 7 deh 4 dnh 3.5 H dn,bs dnh,s - H% 5 p =5GeV dnh 5 H.5 p =GeV dnh H.5.3. dn,bs dnh,s - H% de.5.5 H.5 dn,bs dnh,s - H%.5 deh.4 4 p =5GeV FIG. : The effect of relain the soft-luon approimation on differential variables to the st order in opacity of DGV formalism, as a function of. The comparison of: i differential luon radiative enery loss de /; ii sinle luon radiation spectrum distribution in momentum fraction dn / between bs the solid curve and s the dashed curve case, for different values of initial jet transverse momenta 5 GeV, GeV, 5 GeV, as indicated in panels is shown in the first and second column, respectively. The relative chane of the sinle luon radiation spectrum with respect to soft-luon limit is shown in the third column. different values of initial jet transverse momentum p ; bs curves for de are obtained by numerically interatin Eq. multiplied by E over k and q, whereas s curves correspond to Eq. in [5]. From Fi., we observe a small difference between bs and s results for..3 rouhly up to.4, i.e. for smaller, as epected. We also reconize.3 as a cross-over value, below which differential radiative enery loss and sinle luon radiation spectrum are somewhat lower in bs compared to s case, and above which the opposite is true. At hih value of, i.e..4 <.5, the differences between our bs differential radiative enery loss and previously obtained s [5] ascend to notable values 5% and increases with increasin p. To investiate the effect of relain the soft-luon approimation on the sinle luon radiation spectrum in more detail, the third column is added in Fi. see also dn Fi. 3, showin relative chane of. This quantitative estimation difference smaller than % for..4 is in areement with the previous discussion. In particular, at hiher values, there is a notably larer spectra in bs compared to s case, and this difference enhances up to 6% at p = 5 GeV with increasin p. Never-

8 8 / dn,s dn,bs - % 6 4 p =5 GeV p = GeV p = GeV p =5 GeV R AA is defined as the ratio of the quenched A A spectrum to the p p spectrum, scaled by the number of binary collisions N bin : R AA p = dn AA/dp N bin dn pp /dp. In order to obtain luon quenched spectra, we use eneric pqcd convolution [33]: E f d 3 σ dp 3 f = E id 3 σ dp 3 i P E i E f, FIG. 3: The effect of relain the soft-luon approimation on dn / for different p values. The relative chane of the sinle luon radiation spectrum with respect to soft-luon case, calculated to the st order in opacity of DGV formalism, for different values of initial p as indicated in the leend is depicted as a function of. The curves fade as transverse momentum increases. dn theless, for both variables: de and bs and s cases lead to similar results for.4. The effect of relain the soft-luon approimation on sinle luon radiation spectrum for different transverse momentum values of initial luon jet is further addressed in Fi. 3. We observe that a notable, that is, tenfold increase of p leads to a modest increase less than 5% of dn in bs compared to s case. Note that the same dependence is obtained for de bs dn / de s de since, so that cancels when takin the relative ratio. Therefore, we conclude that the relaation of the softluon approimation has nearly the same effect on and de dn across the whole reion independently on p of the initial jet. Althouh we showed that relain the soft-luon approimation has small numerical impact on both interated E E, N, across the whole reion and differential de dn,, up to.4 variables, the difference between bs and s cases can o up to % and with different sins, and moreover can be quite lare for >.4. This, therefore, leads to a question, how the relaation of the soft-luon approimation affects predictions for measured observables, such as the anular averaed nuclear modification factor R AA [8, 9]. Comparin R AA with and without soft-luon approimation allows assessin how adequate is this approimation in obtainin reliable numerical predictions. To that end, we net concentrate on eneratin the predictions for bare luon R AA, based only on radiative enery loss, with and without soft-luon approimation. where Eid3 σ denotes initial luon spectrum, which is dp 3 i computed accordin to [34, 35], while P E i E f denotes radiative enery loss probability, which includes multi-luon [9] and path-lenth [33] fluctuations. The path-lenth distributions for 5% most central collisions are implemented as described in [9]. Note that we omitted framentation and decay functions, because we are considerin the parton s quenchin, as we are primarily interested in how the relaation of the soft-luon approimation in enery loss affects R AA. Thereupon, we will also investiate how the initial luon distribution influences R AA. Therefore, the left panel of Fi. 4 compares R AA predictions with and without soft-luon approimation accounted, while the percentae chane arisin from relain the approimation is iven by the riht panel of Fi. 4 as a function of the final p. We observe that this relaation barely modifies R AA, in particular the relative chane drops to somewhat less than % at p GeV and further rises to the constant value of %, with increasin p. This very ood areement with even smaller differences compared to previously studied variables between bs and s R AA raises questions of: i why relain the soft-luon approimation has neliible effect on R AA and ii why the lare discrepancy observed in Fis., 3 for hih values does not lead to larer difference in R AA? Reardin i above, we arue that this pattern is epected, as it is well-known that in suppression calculations both E E and N non-trivially affect the R AA. Namely, by comparin the two riht panels of Fi. with the riht panel of Fi. 4 we observe that relain the softluon approimation has opposite effects on E E and, while their interplay is responsible for the nelii- N ble effect on R AA - i.e. the effect on R AA is qualitatively a superposition of the effects on E E and N. To answer ii above, it is convenient to recall that suppression of luon jet see Eq. 3 depends not only on the enery loss probability, but also on the initial luon distribution. In order to intuitively interpret the role of the initial luon distribution, we refer to a descriptive Fi. 5, which represents its dependence on initial transverse momentum. The concept considered is the followin: Some parent luon with unknown initial momentum traverses QGP, looses its enery by luon bremsstrahlun, and emeres with final

9 9 % RAA RAA -.8 % bs RAA s % p HGeV p HGeV FIG. 4: The effect of relain the soft-luon approimation on luon nuclear modification factor RAA versus p. The suppression of luon jet beyond soft-luon approimation the solid curve is compared to soft-luon RAA the dashed curve as a function of transverse momentum in the left panel. The riht panel quantifies the effect and epresses it in percentae. dσ/dp mbgev =. =. 5 =.3 dn RAA. =.4 5 bility correlated with? From Fi. 5 we infer that, due to the eponentially decreasin initial luon momentum distribution, the initial luon correspondin to =. has the hihest probability to be the parent one, and as increases the probability sharply decreases i.e. for &.4 it diminishes for orders of manitude compared to the =. case. This makes only the reion below the crossover noted above..3 relevant for eneratin the predictions, and also eplains why the lare inconsistency between bs and s de or equivalently dn = curves from Fis., 3 at hiher does not affect Additionally, the effect of relain the soft-luon 7 p GeV FIG. 5: The role of initial luon distribution in constrainin relevant reion. The solid black curve represents initial luon distribution as a function of p at the HC [34, 35]. The dot-dashed ray line marks the final luon transverse momentum, while dotted arrows link parent luons, that lost momentum fraction equal to, with their correspondin initial transverse momenta. The arrows fade as increases, as indicated in the leend. VII. momentum p = 3 GeV. This final luon can descend from the parent luon with any p hiher than its own, but we restrict ourselves to 5 different initial momenta, correspondin to different fractional momentum loss {.,.,.3,.4,.5}. For instance, =.5 corresponds to initial luon momentum of 3/.5 GeV= 6 GeV, i.e. to the parent luon that lost half of its momentum etc. The question is which of these 5 luons is the most likely to be the parent one, and how is this proba- approimation on and de is practically insensitive to initial transverse momentum see Fi. 3, which is the reason why finite affects equivalently luon RAA reardless of it s transverse momentum, as observed in Fi. 4. Finally, we also recalculated our finite results, when runnin couplin αs Q, as defined in [36], instead of constant value αs =.3, is introduced in radiative enery loss formula. The obtained predictions lead to the same conclusions as obtained above and are consequently omitted, which supports the enerality of the obtained results. CONCUSIONS AND OUTOOK The main theoretical oal of this paper was to investiate what effect relain of the soft-luon approimation has on radiative enery loss, and consequently on suppression, which depends only on initial distribution and enery loss of hih-momentum parton in QGP. Particularly we chose hih p luon, as due to the color factor of 9/4 compared to the quarks, this assumption affects luons the most. To this end, we analytically calculated all Feynman diarams contributin to the first order in opac-

10 ity radiative enery loss beyond soft-luon approimation, first within GV [9] massless case, and later within DGV [5] massive case, formalism, and numerically predicted: fractional and differential enery loss, number of radiated luons, sinle luon radiation spectrum and luon s suppression. Unepectedly we obtained that, althouh the analytic results sinificantly differ from the correspondin soft-luon results, the numerical predictions are nearly indistinuishable, i.e. within few percents. We then eplained that, due to eponentially decreasin initial luon distribution, only.3 reion effectively contributes to the interated variable predictions. We also showed that neliible suppression chane is due to an interplay between the finite effects on i fractional enery loss and ii number of radiated luons, that have opposite sin. The presented comparisons are done under the assumption of fied stron couplin constant, but also tested with runnin couplin leadin to the same conclusions. Since we showed that luon quenchin in QCD medium composed of static scatterin centers is not affected by the soft-luon assumption, quark radiative enery loss is even less likely to be notably altered, thouh this remains to be further tested. This, to our knowlede, presents the first opportunity to assess the effect of relain the soft-luon approimation on radiative enery loss within DGV formalism. Some other radiative enery loss formalisms, which also imply static scatterers, enerated their results on a finite. However, contrary to the conclusions derived for these formalisms, where sinificant difference in the radiative enery loss was obtained, we found that relain soft-luon approimation brins neliible chane to the results. Consequently, this implies that, within DGV formalism, there is no need to o beyond the soft-luon approimation. Based on the results of this paper, we also epect that the soft-luon approimation can be reliably applied to the dynamical enery loss formalism, as implicitly suested by the previous robust areement [7, ] of our theoretical predictions with a comprehensive set of eperimental data. In particular, the effective cross section vq which corresponds to interaction between the jet and echaned luon [37] does not depend on, so introduction of finite will not affect this term. We also epect that the rest of the enery loss epression i.e. fk, q,, which corresponds to interaction between the jet and radiated luon [37] will be modified in the similar manner as in the static case, since when, these two epression coincide. However, relain the soft-luon approimation in dynamical enery loss model is out of the scope of this paper, and this claim still remains to be riorously tested in the future. Acknowledments: This work is supported by the European Research Council, rant ERC-6-COG: 7574, and by the Ministry of Science and Technoloical Development of the Republic of Serbia, under project numbers ON74 and ON735. Appendi A: NOTATIONS AND USEFU FORMUAS In this paper we used the followin notation for vectors, in consistency with both [9, 5]: p denotes momentum 3D vector p denotes transverse momentum D vector p z denotes component of momentum vector alon the initial jet p = p, p z, p = [p, p, p] denotes momentum 4D vector in Minkowski and iht Cone coordinates, respectively, where p = p p z and p = p p z. For simplicity, we here consider QCD medium consistin of static partons and model the interactions of the luon jet with the medium via static color-screened Yukawa potential, whose Fourier and color structure acquires the followin form [9, 6]: V n = V q n e iqnn = πδq nv q n e i qn n T an R T an n, A v q n = 4πα s q n µ, A where n denotes space-time coordinate of the n th scatterer taret, T an R and T an n denote enerators in SUN c = 3 color representation of jet and taret, respectively, while µ is Debye screenin mass and α s = s/4π is stron couplin constant. In the followin lines we will briefly display the identities and alebra that SUN c = 3 enerators meet: TrT a n =, A3 TrT a it b j = δ ij δ ab C id i d G, A4 where d G = 8 is the dimension of the adjoint representation G. We assume that all taret partons are in the same d T dimensional representation T with Casimir operator C T, while the luon jet is in the adjoint representation G, with Casimir operator C G. In SUN c = 3 color alebra, the followin identities hold as well: [T a, T b ] = if abc T c, A5 while in the adjoint representation we have: T b ab = if abc, A6 T a GT a G = C GI, A7

11 where I denotes identity matri of dimension d G and the SUN c = 3 structure constants f abc are completely antisymmetric to indices permutations, which we frequently apply. In the adjoint representation the followin equalities also stand: CG = C G = N c = 3, A8 Since our etensive calculations are done in pqcd at finite temperature and include only luon interactions, below we list the necessary Feynman rules in covariant aue that we employ: massless luon propaator in Feynman aue: TrT a GT a G = d G C G. A9 And finally, in our computations we frequently make use of the fact that trace is invariant under cyclic permutations and that enerators are Hermitian matrices. a,µ p 3-luon verte: b,ν = iδ ab µν p iɛ, A b,ν a,µ p p = s f abc µρ p p 3 ν µν p p ρ νρ p 3 p µ. A p 3 c,ρ Note that, all diarams in this paper are obtained by usin [38]. Since only physical transverse luon states must be accounted, summin is done accordin to Eq. 57 from [3]: ɛ i kɛ j k = δ ij ki k j, A k pol where i, j =,, 3. Appendi B: ASSUMPTIONS Throuhout the paper we assume that initial luon jet propaates alon the z-ais, i.e. has transverse momentum equal to zero, while radiated luon carries away a finite rate of initial luon lonitudinal momentum and enery, and final luon emeres with momentum p. Therefore, instead of assumin soft-luon approimation, as it was done in [9, 5], we allow to acquire finite non-zero values, thus relain the soft-luon approimation. Since we are calculatin radiative enery loss within the GV DGV formalism apart from abandonin the soft-luon approimation, the followin assumptions remain: The soft-rescatterin approimation. Consistently with [9, 5] we assume that partons eneries and lonitudinal momenta are hih compare to their transverse momenta, which disables the radiated and the final luon to diress much from the initial lonitudinal direction the eikonal approimation. E E E p, k, q i, B The first order approimation. The luon-jet radiative enery loss is calculated up to the first order in opacity epansion, as arued in [9, 9, 3]. Scatterin centers distribution and ensemble averae. We consider that all scatterin centers i are distributed with the same transversely homoeneous density: ρ = N A ρz, B where z ρz = and also that impact parameter i.e. relative transverse coordinate b = i alters within a lare transverse area A compared to the interaction area µ. Therefore, the ensemble averae over the scatterin center locations reduces to an impact parameter averae:... = b A..., B3 which in our case is mainly used in the followin form: e iqiqjb = π A δ q i q j. B4

12 We also assume that the enery of initial hard probe is lare compared to the potential screenin scale: E, E, E µ. B5 Net, we assume that the distance between the source J and the scatterin centers is lare relative to the interaction lenth: z i z µ, B6 then, that source current varies slowly with momentum: Jp k q Jp k, B7 and that the source current can be written eplicitly in terms of polarization vector: J µ a p k q J a p k qɛ µ p k q J a p kɛ µ p k q. B8 In the followin sections first we assume that luons are massless GV in order to make the comprehensive derivations more straihtforward and easier to follow, but later we recalculate all the results with luon mass [7] included DGV Appendi J. Appendi C: Gluon jet M First we calculate luon-jet radiation amplitude to emit a luon, carryin a finite fraction of initial jet enery, with momentum, polarization and color k, ɛk, c and without interactions with the medium M. We assume that initial luon p k propaates alon z-ais. By usin M amplitude as an eample, we will implement the aforementioned assumptions in order to acquire momentum and polarization epressions. Thus, the initial luon 4-momentum reads: p k = p k, p z k z,, p k = [E, E, ], C where E = p k p z k z and E = p k p z k z. Assumin massless real luons for simplicity, the momentum vectors of the radiated k and the final p luons acquire the followin form: k = k = [E, p = p = [ E, k E, k], C p, p]. C3 E We also assume that luons are transversely polarized particles. Althouh we work in covariant aue, we can choose any polarization vector for the eternal on-shell M c,ρ k J µ a z p k z a,µ p d,σ FIG. 6: Zeroth order diaram that includes no interaction with the QCD medium, and contributes to luon radiation amplitude to the first order in opacity /λ. The dashed circle represents the source J, which at lonitudinal coordinate z produces an off-shell luon jet, propaatin alon z-ais. z denotes lonitudinal coordinate at which the luon is radiated. atin indices denote color chares, while Greek indices denote components of 4-vectors. k denotes 4-momentum of the radiated luon carryin the color c, and p denotes 4-momentum of the final luon jet carryin the color d. luons, so in accordance with [9, 4, 5] we choose n µ = [,, ], as stated above: ɛk k =, ɛk n =, ɛk =, ɛp p =, ɛp n =, ɛp =, C4 while we assume that the source has also the physical polarizations as the real luons [4]: ɛp k p k =, ɛp k n =, ɛp k =. C5 Usin Eqs. C- C5 we can now obtain the followin epressions for the luon polarizations: ɛ i k = [, ɛ i k E, ɛ ɛ i p i], ɛ i p = [, E, ɛ i], ɛ i p k = [,, ɛ i ], C6 where i =, counts for polarization vectors. Note that the 4-momentum is conserved, which leads to the relation: p k =, C7 that we implement in Eqs. C3, C6 in order to ensure that everythin is epressed in terms of k. Also, E E, E k = E, where E = p k is the enery of initial jet. Usin the notation from Fi. 6 we obtain:

13 3 M = ɛ σpɛ ρk s f acd µσ p k ρ µρ p k σ ρσ p k µ iδ aa µµ p k iɛ ij a p keipk ɛ µ p k J a p ke ipk s ɛ k k f acd = J a p ke ipk i s ɛ k k T c da. C8 Eq. C8 after summation by usin Eq. A ives: M = Jp k 4 s C Gd G k. Net we substitute the Eq. C9 in: d 3 N d 3 N J Tr M d 3 p d 3 k π 3 p π 3 ω, C9 C Note that, contrary to the soft-luon approimation [5], where: d 3 N J d G Jp k d 3 p π 3 p, C now p, denotin the momentum of the final luon jet, is not approimately equal to the momentum of initial luon jet i.e. the radiated luon can carry away a substantial amount of the initial jet enery and lonitudinal momentum. Thus instead of Eq. C throuhout this paper we use the eneral one: d 3 N J = d G Jp k d 3 p J π 3, E J C where E J = E and p J denotes enery and 3D momentum of the initial luon jet, respectively. Knowin that the substitution of variables p z, k z p J z, E ives: d 3 p d 3 k π 3 p π 3 ω = d3 p J π 3 E J d k π 3, C3 and by substitutin Eqs. C9, C, C3 in Eq. C, for radiation spectrum we now obtain: d 3 N dk = α s π C G k, C4 which recovers well-known Altarelli-Parisi result [3] and for reduces to the massless soft-luon limit of Eq. 9 from [5]. We have also checked that, the same result can be obtained by directly implementin polarization vectors Eq. C6 in Eq. C8, instead of usin Eq. A when averain. Appendi D: Diarams M,,, M,,, M,, In this section we provide a detailed calculations of Feynman amplitudes, correspondin to luon-jet interaction with one scatterin center, which are depicted in Fi. 7. Aain for consistency, we assume that initial jet p k q propaates alon z-ais. Throuhout this section momentum and polarization vector for initial luon read: p k q = [E q z, E q z, ], ɛ i p k q = [,, ɛ i ], D D where q = [q z, q z, q ], with q =, denotes the momentum of echaned luon, while p, k and correspondin polarization vectors retain the same epression as in Eqs. C, C3, C6, with the distinction that the followin relation between luon transverse momenta, due to 4-momentum conservation, holds: q = p k.. Computation of M,, diaram D3 We chose to start with thorouh derivation of the epression for M,, amplitude, simply because it has no counterpart reardin the symmetry under p k,, c d substitutions, and it provides all necessary steps for calculatin the remainin two amplitudes from this chapter, apart from havin one less sinularity compared to the amplitudes M,, and M,,. Thus, usin the notation from the left diaram of Fi. 7, we write:

14 4 M,, k c,ρ M,, k,c M,, J µ a z z z q a,µ b,ν b,ν p d,σ J z z z p,d J z z z p,d k,c a, q,a q,a FIG. 7: Three diarams, correspondin to interaction with one static scatterin center, that contribute to luon-jet radiation amplitude to the first order in opacity /λ. z denotes lonitudinal coordinate of the interactions with one scatterin center. Crossed circle represents scatterer that echanes 3D momentum q with the jet. Note that, all three diarams assume equivalently ordered atin and Greek indices as indicated by the first diaram. Remainin labelin is the same as in Fi. 6. M,, = 4 q π 4 ɛ σpɛ ρk s f bcd νσ p k ρ νρ p k σ ρσ p k ν iδ bb νν p k iɛ f ab a µ p k q ν µν p k q ν p k q µ T a V q e iq iδ aa µµ p k q iɛ ij a p k q ɛ µ p k q e ipk q J a p ke ipk f bcd f aab T a i E dqz π vq z, q e iqzz z p k q iɛp k iɛ, q e iq ɛ k ɛ p π s D4 where we used the equation: p k = k p, D5 and assumed that J varies slowly with momentum q, i.e. Eq. B7. The lonitudinal momentum transfer interal: I p, k, q, z z dqz π vq z, q e iqzz z p k q iɛ D6 has to be performed in the lower half-plane of the comple plain, since z > z. In order to determine the pole arisin from potential, we rewrite Eq. A in a more appropriate form: v q n = 4πα s q nz iµ n q nz iµ n, D7 where µ n = µ q n, with n denotin the correspondin scatterin center. Aside from the pole oriinatin from Eq. D7 q z = iµ, there is also a sinularity emerin from the luon propaator: q = k E p E iɛ = k ω k q iɛ, ω D8 where ω = k E. The residue around the pole at q is computed as the neative sin is due to the clock-wise orientation of the closed contour in the comple plain: Res q v k E p E, q i = v k ω k q ω k E ei, q i E e i p z E E z ω k k q z z. D9

15 The pole oriinatin from the potential q z = iµ does not contribute to the lonitudinal interal, since residue around that pole is eponentially suppressed due to Eq. B6, i.e. µz z = µλ and µ µ : where we assumed that E approimation. 5 µ and soft-rescatterin 4πα s Res iµ i iµ E iµ e µ z z, D This makes only q sinularity relevant for calculatin lonitudinal interal. Therefore I coincides with Eq. D9, i.e.: I p, k, q, z z v k E p E, q i = v, q i E e i k E ei ω k k q z z, p z E E z v, q i k E ei p z E E z D where we used eikonal approimation i.e. for a finite amplitude reads: k : p E and E. Finally, M,, M,, = J a p ke ipk i f bcd f aab q T a π v, q ɛ k p e iq b i s k p k i e p z E E z =J a p ke ipk i T c T a da T a q π v, q e iq b i s ɛ k q k q e i ω k k q z z, D where we denoted b. In this subsection, we constantly make use of Eq. D3 in the followin form:. Computation of M,, and M,, diarams p = k q, D3 and also manipulate with SUN c = 3 structure constants by usin Eqs. A5, A6. Note from Fi. 7 that, as epected, M,, is symmetric under the substitutions: p k,, c d, where the symmetry can be straihtforwardly verified by implementin these substitutions in the first two lines of Eq. D. Applyin the same procedure as in the previous subsection, we proceed with calculatin M,,. Note that the order of the color and Dirac indices denotin vertices is the same for all three diarams in Fi. 7, and are therefore omitted in the last two diarams. M,, = 4 q π 4 ɛ σpf bda ν p q σ νσ p q σ p q ν T a V q e iq iδ bb νν p q iɛ s f acb µν p k q ρ µρ p k q ν ρν p k q µ ɛ iδ aa µµ ρk p k q iɛ ij a p k q ɛ µ p k q e ipk q J a p ke ipk f bda f acb T a i E q π e iq b s ɛ k I, D4

16 where: I p, k, q, z z dqz vq z, q e iqzz z π p k q iɛ p q iɛ. D5 6 we aain have to close the contour below the real ais. Similarly as in M,, amplitude, aain only the poles oriinatin from the propaators contribute to the interal: k p E E iɛ and k p E iɛ, while iµ is eponentially suppressed due to µz z. Therefore we obtain: In order to calculate the previous interal, due to z > z I p, k, q, z z i E k v, q i E k v, q k i e p z E E z e i p k e i ω k k q z z e i E z z ω k q k z z, D6 leadin to: M,, = J a p ke ipk i f bda f acb q T a π v, q e iq b i s ɛ k k k i e p z E E z e i p k z E z = J a p ke ipk i T a T c q da T a π v, q e iq b i s ɛ k k e i ω k k q z z e i ω k k q z z. D7 By applyin similar procedure for M,, we obtain: M,, = J a p ke ipk i f bca f adb q T a π v, q e iq b i s ɛ p p k i e p z E E z e i k p z E z = J a p ke ipk i [T c q, T a ] da T a π v, q e iq b i s ɛ k q k q e i ω k k q z z e i ω k k q z z. D8 Notice from Fi. 7 that M,, and M,, are symmetric under the followin substitutions: p k,, c d; it can be straihtforwardly verified that Eqs. D7, D8 are symmetric under these substitutions. Appendi E: Diaram M,, Net we concentrate on the diarams containin two interactions with the static scatterin centers, since they also contribute to the luon radiative enery loss to the first order in opacity, when multiplied by M. There are seven such diarams, that we ather into four roups,

17 each of which contains two or one diarams symmetric under p k,, c d substitutions. For consistency the initial luon jet with momentum p k q q propaates alon z-ais, i.e.: luon transverse momenta holds: p k = q q. 7 E3 p k q q = [E q z q z, E q z q z, ], E ɛ i p k q q = [,, ɛ i ], E where q i = [q iz, q iz, q i ], i =, with qi = denote momenta of echaned luons, while p, k and correspondin polarizations retain the same epressions as in Eqs. C, C3, C6, with distinction that, due to 4- momentum conservation, the followin relation between Aain, from seven diarams we chose one model diaram M,,, based on the same reason as in Appendi D, for thorouh derivation of the final amplitude epression. From Fi. 8, where luon jet after two consecutive interactions with scatterin centers radiates a luon with momentum k, we observe that there are two limitin cases that we consider. Usin the notation from Fi. 8 we write: M,, = 4 q π 4 d 4 q π 4 ɛ σpɛ ρk s f ecd ξσ p k ρ ξρ p k σ ρσ p k ξ iδ ee ξξ p k iɛ f be a ν p k q ξ νξ p k q ξ p k q ν T a V q e iq iδ bb νν p k q iɛ f ab a µ p k q q ν µν p k q q ν p k q q µ T a V q e iq iδ aa µµ p k q q iɛ ij a p k q q ɛ µ p k q q e ipk q q ij a p ke ipk f ecd f bea f aba T a T a i e iq b e iq b E dqz π dq z π q d π i q π i ɛ k p s k p vq z, q vq z, q e iqzz z e iqzz z p k q q iɛp k q iɛ, E4 where b i i, i =, denote transverse impact parameters. We used Eq. D5 and assumed that J varies slowly with momentum q i, i.e. Jp k q q Jp k. Reardin the lonitudinal q z interal, we introduce a new variable: q z = q z q z throuhout this, and the followin sections involvin Feynman amplitudes which include interactions with two scatterin centers. Therefore, we rewrite the eponent in the followin manner: e iqzz z e iqzz z = e iqzz z e iqzz z. Rewritin q z lonitudinal interal in terms of q z, i.e. chanin the variables, we obtain: I p, k, q, q, z z = qz π vq z q z, q e iqzz z p k q q. iɛ E5 Aain, due to z > z, the contour must be closed in the lower half-plane of comple q z plain, so additional minus sin arises from the neative orientation of the contour and also we nelect the pole at q z = iµ q z, since it is eponentially suppressed due to Eq. B6. Thus, only one pole, oriinatin from the luon propaator, contributes to the first lonitudinal interal: q = k E p E iɛ = k ω k q q iɛ, ω E6 where we used, as well as throuhout the Appendices F-H the relation between transverse momenta Eq. E3. The residue at Eq. E6 then ives: