# Energy and momentum deposited into a QCD medium by a jet shower

## Abstract

Hard partons moving through a dense QCD medium lose energy by radiative emissions and elastic scatterings. Deposition of the radiative contribution into the medium requires rescattering of the radiated gluons. We compute the total energy loss and its deposition into the medium self-consistently within the same formalism, assuming perturbative interaction between probe and medium. The same transport coefficients that control energy loss of the hard parton determine how the energy is deposited into the medium; this allows a parameter free calculation of the latter once the former have been computed or extracted from experimental energy loss data. We compute them for a perturbative medium in hard thermal loop (HTL) approximation. Assuming that the deposited energy-momentum is equilibrated after a short relaxation time, we compute the medium’s hydrodynamical response and obtain a conical pattern that is strongly enhanced by showering.

###### pacs:

12.38.Mh, 11.10.Wx, 25.75.DwJet quenching (the modification of hard jets in dense media) is one of the most studied discoveries at the Relativistic Heavy-Ion Collider (RHIC) RHIC_Whitepapers (). It is expected to play a key role in the study of the quark-gluon plasma (QGP) produced in heavy-ion collisions at the Large Hadron Collider (LHC). Numerous experiments highpt () have established the suppression of hadrons with high transverse momenta; others indicate that the lost energy manifests itself as conical flow in the soft sector back_to_back ().

Calculations of jet modification tend to focus on one of two separate questions: the modification of the final hadron distribution from the hard parton due to its energy loss, or the response of the medium to the energy deposited. Numerous studies of the former, based on perturbative QCD (pQCD), have yielded near-rigorous measures of the two non-perturbative transport coefficients and which codify the transverse (to the jet axis) momentum diffusion and longitudinal drag experienced by a fast parton Majumder:2007iu (). Computations of the medium response consist of two parts: an ansatz for the space-time profile of the energy-momentum deposition, and a calculation of the dynamical response to this “source” of excess energy and momentum. Based on its success at RHIC, ideal fluid dynamics has been used to compute this medium response hydro_mach (), assuming that the energy lost by the jet is entirely deposited into the medium at a constant rate and thermalizes instantaneously.

So far there exists no first principles calculation of the magnitude and space-time profile of the energy-momentum deposited in a medium by a hard parton that can be considered on par with the pQCD energy loss calculations string_mach (). A noteworthy attempt to calculate the deposition profile in pQCD is the semi-phenomenological approach of Neufeld and Müller Neufeld:2009ep () who use the differential single gluon emission spectrum of Ref. radiative () and interpret this as the rate of gluon emission in the medium. A non-diffusive Fokker-Planck equation is then motivated to compute how this distribution changes due to elastic energy loss of the emitted gluons. As anticipated in Majumder:2004pt (), they find that not all of the energy lost to gluon radiation is deposited in the medium. However, since the underlying formalism radiative () lacks information about virtuality evolution, this calculation does not include gluon multiplication by showering, i.e. the splitting of a radiated gluon into two lower virtuality gluons. The transverse momentum deposition thus cannot be computed and, due to the strict eikonal limit used in radiative (), the parent parton does not lose energy after radiation. We present a new formalism in which the radiative and elastic energy loss of the fast parton, its virtuality evolution by radiation, the showering and multiplication of the radiated gluons, and the energy deposited by them in the medium are all calculated consistently in the same approach.

Hard jets in vacuum or in heavy ion collisions are produced with considerable virtuality. As the jets proceed through vacuum or medium, this virtuality is lost by sequential radiative emissions. The effect of this perturbative shower on the non-perturbative hadronization process is computed using DGLAP evolution equations AP () for the fragmentation function. These equations express the radiation of multiple partons, which hadronize independently, via an evolution in virtuality of the parent parton. In a medium, one can derive analogous equations where the gluon radiation probability is modified by the scattering of hard parton and emitted gluons off medium constituents. These are referred to as “medium modified evolution equations”. In addition to stimulating gluon emission, the scattering of the hard parton causes it to lose forward light-cone momentum by elastic exchanges with the medium elastic (); Majumder:2008zg (). At the same time the parton gains transverse momentum from the medium Majumder:2007hx () and imparts to it an equal amount in return. In an arbitrary medium these effects are encoded in two non-perturbative transport coefficients, and , defined in terms of in-medium gluon field correlation functions Majumder:2007hx (); Majumder:2008zg (). The medium modification of the standard vacuum evolution depends on these transport coefficients.

In-medium evolution equations where the medium modified fragmentation function (MMFF) is affected only by were derived in Majumder:2009zu (). We point out that the same processes can be used to compute the amplification of the energy deposited through multiple radiations stimulated by transverse broadening. Formally, this can be computed by replacing the operator expression for the fragmentation function with that for the energy deposited; this is identical to . Using this calculation of the modified as the energy deposited through all elastic scatterings of the shower places it on the same footing as energy loss. The diagrams involved and the resulting expressions for the in-medium splitting functions (IMSF) are identical. The solution of the in-medium evolution equation for no longer represents the elastic energy loss by one parton, but rather the energy deposited by the jet shower.

Imagine a hard quark or gluon with large light cone momentum (and thus energy ) and virtuality entering a medium of fixed length held at a constant temperature . Let us assume that the rate of energy deposition by this jet in the medium as a function of length , denoted as , is known (i.e., can be calculated or measured). Note that both the deposited energy and are actually the light-cone quantities and . For brevity we refer to these simply as deposited energy and distance travelled. Given the above function, the total energy deposited by a jet originating at location and propagating to is given as

(1) |

where the last approximate equality is solely for the case of a single parton propagating without radiation.

If the scale is much larger than , the change with virtuality in the partonic shower pattern may be calculated perturbatively: a leading quark at the higher virtuality may split into a quark and a gluon with lower virtuality, and similarly for a gluon. As a result, there is change in the energy deposited in the medium due to the increase of the number of partons depositing energy. Using the IMSF from Majumder:2009zu (), the change in the energy deposition by a quark with energy from to due to the increase in virtuality can be expressed as tobepublished ()

(2) | |||||

Here the first term in square brackets represents the energy deposited by a quark with energy and virtuality , from the initial location to the intermediate location ; the second and third terms represent the energy deposited by the quark and the emitted gluon with reduced energies and , respectively, from the intermediate location to the final location . In Eq. (2) the quark IMSF is given as Majumder:2009zu (); HT ()

(3) |

The increase in the energy deposited due to the splitting of the parton is reduced by the virtual correction which restores unitarity to the evolution equations. The effect of such corrections on Eq. (2) is incorporated by subtracting from it the virtual term

(4) |

Along with the energy deposition from a quark jet one has to evolve the one from a gluon jet of virtuality , using a similar evolution equation that includes the splitting of a gluon into two gluons or a pair. Similar to the MMFFs, one solves a coupled set of evolution equations for and both of which are functions of three variables at the scale .

The evolution equations (2,3,4) for a quark jet and the coupled equations for gluon jets are motivated by existing rigorous derivations of the medium modification of the fragmentation functions due to gluon emission HT (), the accumulation of transverse momentum Majumder:2007hx () and longitudinal drag Majumder:2008zg () by propagating hard partons in a QCD medium, and the effect of such accumulated momentum on radiative processes Majumder:2007ne (). The IMSF (3) accounts for interference between diagrams where the gluon is emitted at the origin or at the location . In propagating up to the quark loses a fraction of its energy; while this is included in the total energy deposited, its effect on the interference pattern in Eq. (3) is ignored; this is justified in the eikonal limit for the propagating parton. Yet another approximation is the neglect of the energy lost by the radiated (reabsorbed) gluon in the virtual correction. Since the radiated gluon in the virtual correction exists in only one amplitude, with a single parton in the complex conjugate, its energy loss is balanced by the quark propagating in the loop.

In the eikonal approximation, the hard jet loses light cone momentum and remains close to on-shell, thus the -component of the deposited light-cone momentum is approximately equal to the energy deposited (). Note that the negative light cone momentum () is not conjugate to and thus it is not inconsistent to compute the dependence of the deposited. The remaining two components that may be computed are the transverse momentum deposited by the jet as a function of . This can again be directly estimated from a pQCD calculation: A parton traversing a medium gains transverse momentum squared with length as

(5) |

By momentum conservation this equals the deposited in the medium by the same parton.

For a hard virtual quark the total transverse momentum deposited increases due to parton splitting. This can be calculated using an equation similar to that for light-cone momentum deposition. For a quark with energy and virtuality , traversing a medium from to , the change of the transverse momentum deposited with virtuality is obtained as

(6) | |||||

The splitting function here is identical to that in Eq. (3), and the meaning of the three terms in the bracket is analogous to Eq. (2). Further, one must include a virtual correction and couple Eq. (6) to a similar equation for the deposited by a virtual gluon.

Using Eqs. (2,6) (along with the coupled ones for gluon jets), we can compute the 3-momentum deposited by a hard virtual parton, disintegrating into a shower of partons, in a dense medium as a function of the length traversed. Similar to the case of in-medium evolution equations for the MMFF Majumder:2009zu (), these equations require an initial condition. For the case of the MMFF, the only possible choice was to insist that the part of the jet with virtuality below a minimum exited the medium and use the known vacuum FF at that scale as an input. For the deposited part of the energy-momentum we here assume that the medium is weakly coupled, thus when the virtuality of the parton is the deposited energy and can be obtained from the expressions for and in an HTL plasma elastic ():

(7) |

Here is the Debye screening length and is the representation specific Casimir. The integrated energy and deposited from Eq. (7), as a function of the length traveled, is plotted for gluons (circled) and quarks as solid lines in Figs. 1 and 2. For a consistent description, we impose that partons with an energy become part of the thermal medium. This condition is maintained through the evolution equations.

Using Eq. (7) as input, we may calculate the increase in the energy and deposition in the medium as a function of for initially highly virtual hard partons that evolve into a radiative shower. Starting from the scale of (in all calculations we pick MeV and a partonic plasma with 3 quark flavors) we evolve up to an initial scale . These are plotted as dashed lines in Figs. 1 and 2. One notes immediately that both quantities increase as we evolve up in virtuality. For comparison, we also estimate the total energy lost by the hard parton due to elastic, radiative inelastic and flavor changing interactions (dash-dotted lines in Fig. 1). The last type of energy loss refers to the case where a quark splits with the gluon carrying a larger fraction of the momentum, or a gluon splits into a quark-antiquark; in this case we assume that the entire energy of that parent parton has been lost. This leads to a somewhat artificial enhancement of the total energy loss.

As an illustration of the effect of this energy-momentum deposition in the medium, we compute its hydrodynamic response to the following source term:

(8) |

In this first attempt we ignore the transverse momentum contribution to the source current. Following Refs. linear_hydro (), we assume that the energy deposited is a small perturbation and solve for the linear response of the medium:

(9) |

is the unperturbed energy-momentum tensor of a homogeneous and static partonic medium in equilibrium. The small excess is decomposed as

(10) | |||||

is the excess energy density, is the momentum current density and is the sound attenuation length. For the specific shear viscosity we took . We delay the response to the source by a time to account for thermalization of the deposited energy.

In Fig. 3 we show the azimuthal projection of the energy density at fm/c after the parton is created, for a single non-radiating parton (left) and a parton-initiated jet shower (right). A gluon (bottom row) deposits more energy than a quark (top row), due to its larger color factor that enters both in the elastic energy loss and shower production rate. One immediately notes that, while the basic Mach cone structure is not changed, showering leads to an enhancement by a factor of 3 in the overall magnitude of the response. For quark jets our results are qualitatively similar to Ref. Neufeld:2009ep ().

In this Letter, we have presented a consistent pQCD based calculation of the light-cone and transverse momentum () deposited by a jet in a medium, as a function of distance traversed. Assuming a short thermalization time for the deposited energy we also computed the hydrodynamic response. The pQCD shower has the effect of a large part of the energy being deposited later in the history of the jet Neufeld:2009ep () which tends to enhance the Mach cone like structure formed.

We thank B. Müller and R. B. Neufeld for helpful discussions. This work was supported by the U.S. Department of Energy under grant DE-FG02-01ER41190.

### References

- I. Arsene et al., Nucl. Phys. A 757, 1 (2005); B. B. Back et al., ibid. 757, 28 (2005); J. Adams et al., ibid. 757, 102 (2005); K. Adcox et al., ibid. 757, 184 (2005).
- K. Adcox et al., Phys. Rev. Lett. 88, 022301 (2002); C. Adler et al., Phys. Rev. Lett. 89, 202301 (2002).
- A. Adare et al., Phys. Rev. Lett. 98, 232302 (2007); B. I. Abelev et al., Phys. Rev. Lett. 102, 052302 (2009).
- A. Majumder, J. Phys. G 34, S377 (2007).
- A. K. Chaudhuri and U. Heinz, Phys. Rev. Lett. 97, 062301 (2006); T. Renk and J. Ruppert, Phys. Lett. B 646, 19 (2007); B. Betz, et al. arXiv:0812.4401 [nucl-th].
- For an alternate approach that assumes that the hard parton couples strongly to the medium see P. M. Chesler and L. G. Yaffe, Phys. Rev. Lett. 99, 152001 (2007); S. S. Gubser, S. S. Pufu and A. Yarom, ibid. 100, 012301 (2008).
- R. B. Neufeld and B. Müller, Phys. Rev. Lett. 103, 042301 (2009).
- U. A. Wiedemann, Nucl. Phys. B 588, 303 (2000); C. A. Salgado and U. A. Wiedemann, Phys. Rev. D 68, 014008 (2003).
- A. Majumder, E. Wang and X. N. Wang, Phys. Rev. Lett. 99, 152301 (2007).
- V. N. Gribov and L. N. Lipatov, Sov. J. Nucl. Phys. 15, 438 (1972); Yu. L. Dokshitzer, Sov. Phys. JETP 46, 641 (1977); G. Altarelli and G. Parisi, Nucl. Phys. B 126, 298 (1977).
- E. Braaten and M. H. Thoma, Phys. Rev. D 44, 2625 (1991). S. Wicks, et al., Nucl. Phys. A 784, 426 (2007); G. Y. Qin, et al., Phys. Rev. Lett. 100, 072301 (2008).
- A. Majumder, arXiv:0810.4967 [nucl-th].
- A. Majumder and B. Müller, Phys. Rev. C 77, 054903 (2008).
- A. Majumder, arXiv:0901.4516 [nucl-th].
- A. Majumder, to be published.
- X. F. Guo and X. N. Wang, Phys. Rev. Lett. 85, 3591 (2000); Nucl. Phys. A 696, 788 (2001).
- A. Majumder, R. J. Fries and B. Müller, Phys. Rev. C 77, 065209 (2008).
- J. Casalderrey-Solana, E. V. Shuryak and D. Teaney, J. Phys. Conf. Ser. 27, 22 (2005); R. B. Neufeld, Phys. Rev. C 79, 054909 (2009).