Hydrodynamics and jets in dialogue

Energy and momentum loss of jets in heavy ion collisions can affect the fluid dynamic evolution of the medium. We determine realistic event-by-event averages and correlation functions of the local energy-momentum transfer from hard particles to the soft sector using the jet-quenching Monte-Carlo code Jewel combined with a hydrodynamic model for the background. The expectation values for source terms due to jets in a typical (minimum-bias) event affect the fluid dynamic evolution mainly by their momentum transfer. This leads to a small increase in flow. The presence of hard jets in the event constitutes only a minor correction.


Introduction
For a phenomenological, fluid dynamic description of heavy ion collisions it is usually assumed that the bulk of the medium produced after a heavy ion collision is in local thermal equilibrium. While this may be a reasonable approximation for the low-p ⊥ part, it is phenomenologically clear that high-p ⊥ particles do not originate from a locally equilibrated or thermal distribution. It is an interesting question how this non-equilibrated part influences the hydrodynamical evolution of the bulk. The most important effect in this regard might be due to the energy loss of high-p ⊥ particles propagating in the medium (jet quenching) which leads to a transfer of energy and momentum to the bulk part described by fluid dynamics. This could have important implications for the interpretation of soft observables (e.g. the anisotropic flow coefficients v n ) as well as jet measurements [1], which rely on background subtraction techniques assuming currently that the soft event and jets are uncorrelated. On the other side one expects that local fluid properties determine the strength of the energy loss. Generally speaking, energy loss is expected to be stronger for a denser medium. More a e-mail: stefan.floerchinger@cern.ch b e-mail: korinna.zapp@cern.ch specific, the jet-quenching parameterq is expected to depend on the temperature and other parameters of the medium.
The energy loss of hard partons due to induced gluon radiation in a hydrodynamic background has been studied in various approaches [2][3][4][5]. They are, however, not suited for quantifying the energy and momentum deposition into the bulk. Firstly, they operate in a high energy limit, where there is no collisional energy loss, and secondly, they do not keep track of radiated gluons. It thus has to be assumed that all radiated energy gets dissipated locally, which is clearly a bad approximation for energetic emissions. This is different in Monte Carlo codes aiming at a consistent description of the entire jet and its interactions in a background described by hydrodynamics [6][7][8]. They can trace all radiated partons, but also here the interactions between the jets and the bulk are accounted for in an effective way that cannot easily be translated into a local energy and momentum transfer between jets and background.
The influence of jets on the hydrodynamic evolution of the bulk was first discussed in the context of Mach cone formation [9][10][11][12][13][14] and has been extended recently to other observables [15,16]. Within AdS/CFT the interplay between the energy loss of a heavy quark and hydrodynamic excitations has been discussed in detail, see [17,18] for an overview. In that context the source term can be extracted unambiguously from the energy and momentum loss of a (single) heavy quark. Arguably, holographic models provide currently the only consistent field theoretic framework for calculating source functions from first principles. The weakness of this approach is, however, that it is unclear to what extent it captures dominant features of jets propagating in a QCD medium.
The studies in QCD, on the other hand, feature an elaborate treatment of the hydrodynamic side of the problem, but have very simplified models for the energy and momentum deposited by the jets. In general, a fully self-consistent description of the soft medium and high-p ⊥ part of the spectrum with its mutual interactions in QCD can be a rather difficult task (first steps in this direction are taken by transport codes [19][20][21]). On the other side, the phenomenological success of the current fluid dynamic model, which neglects the influence of non-equilibrated hard particles completely, suggests that the influence of the latter is not too large. In order to investigate this question more quantitatively, we employ here a nonself-consistent description where the bulk medium is first described in terms of conventional fluid dynamics neglecting non-thermal components. This leads in particular to a temperature and fluid velocity profile as a function of the space-time coordinates. In a second step we use these results to estimate the local transfer of energy and momentum from the hard particles to the medium. This results effectively in an additional source term in the fluid dynamic evolution equations. The influence of this source for fluid dynamics can then be estimated in a third step. The effect of a fourth step, namely re-calculating the jets in the modified background, is expected to be numerically small and can thus be neglected.
In this study the jets are simulated with Jewel [22], which employs a microscopic description for the interactions of the jets in the background. As a first approximation one can thus interpret the energy and momentum flow in the individual scattering processes as the energy-momentum exchange between the jets and the background. This provides a realistic and well constrained model for the local energy-momentum transfer to the bulk. Moreover, Jewel provides a realistic ensemble of jets including their spatial and kinematic distributions and event-to-event fluctuations, e.g. in the jet fragmentation. This allows us to study statistical properties, in contrast to earlier approaches that have mainly focused on the medium response to single (simplified) jets. This paper is organized as follows. In Sect. 2 we discuss the fluid dynamic evolution of the bulk including source terms for energy and momentum transferred from the hard sector. In Sect. 3 we introduce the description of jets and in the subsequent Sect. 4 we quantify the local energy and momentum transfer in terms of expectation values and correlation functions. Finally, we draw some conclusions in Sect. 5.

The hydrodynamic evolution
In this section we discuss our formalism on the fluid dynamic side in more detail. We start from the fluid dynamic expression for the energy-momentum tensor T μν bulk = ( + p + π bulk )u μ u ν + ( p + π bulk )g μν + π μν , (2.1) where is the energy density, p is the pressure, u μ is the fluid velocity, π μν is the shear stress tensor and π bulk is the bulk viscous pressure [we use signature (−, +, +, +)]. Equation (2.1) accounts for the bulk contribution to the total energy-momentum tensor and gets supplemented by a similar contribution from the non-equilibrated, hard part of the medium T μν hard . The total energy-momentum tensor is conserved, We define the effective source terms for the bulk evolution by such that the energy-momentum conservation equation becomes When projected in the direction of the fluid velocity, this leads to the evolution equations for the energy density (2.5) The third and fourth term on the left hand side of Eq. (2.5) give the increase in thermal energy due to the shear and bulk viscous dissipative effects, respectively. This can be understood as a transfer of energy from the mechanical motion of the fluid to its thermal energy. The second law of thermodynamics implies that this leads to an increase of entropy. In a first order fluid dynamic formalism, the constitutive relations is the projector to the transverse and traceless part and μν = u μ u ν + g μν is the projector orthogonal to the fluid velocity) ensure that this is indeed the case. In a very similar way, the term on the right hand side of Eq. (2.5) accounts for a change of the thermodynamic internal energy (and therefore enthalpy and entropy) due to energy loss of high-momentum particles. This is a dissipative process, as well, and a consistent fluid dynamic description with the second law of thermodynamics requires In the local fluid rest frame where u μ = (1, 0, 0, 0) this implies J 0 ≥ 0, i.e. energy must be transferred from the nonequilibrated particles to the fluid (and not the other way). We will check below to what extent Eq. (2.7) is indeed satisfied in currently used Monte-Carlo simulations of jet energy loss.
Let us now consider the remaining set of equations that follows from the conservation of energy-momentum, Eq. (2.4).
It is obtained by projecting to the direction orthogonal to the fluid velocity and yields an evolution equation for the latter, In this equation, the second term on the left hand side accounts for an acceleration of the fluid due to pressure gradients. The third term accounts for the change in the fluid velocity due to dissipation of macroscopic kinetic energy into thermal energy. This leads usually to a damping of the fluid motion. The term on the right hand side of Eq. (2.8) is a force term that accounts for the acceleration of the fluid due to high energetic particles that propagate through it. This is the force opposing drag.
The fluid dynamic evolution equations (2.5) and (2.8) have to be supplemented by constitutive equations for the shear stress tensor and bulk viscous pressure. In a first order (Navier-Stokes type) formalism these are of the form (2.6), in a second order formalism these equations get supplemented by relaxation time terms. To solve the evolution equations one also needs an equation of state that relates pressure and energy density as well as the transport coefficients η and ζ (and possible further coefficients such as relaxation times).
So far, we have not yet specified the source terms on the right hand side of Eqs. (2.5) and (2.8). If these correspond to high-momentum, non-equilibrated particles, they are in general different for each event. One might attempt at this point to implement an event-by-event description of fluid dynamics and a model for the high-momentum particles coupled to each other. On a technical level this becomes quickly rather involved. There is also a conceptual difficulty of drawing a line between the high-momentum part of the medium that is usually described in a microscopic way in terms of single particle excitations or partons and the low-momentum part that is described in a more macroscopic way in terms of fluid dynamics.
We follow here another approach that uses a statistical description also for the non-equilibrated high-momentum part. More specific, we describe the influence of the nonequilibrated part of the medium onto the fluid dynamic variables and evolution in terms of a statistical ensemble of sources J ν or, equivalently, of a source component parallel to the fluid velocity, and orthogonal to it, as they appear on the right hand side of Eqs. (2.5) and (2.8). A particular configuration for a single event corresponds to one element of this ensemble. One possibility to characterize such an ensemble is in terms of a functional probability density which is a functional of the source components J S (x) and J V (x) for a single event. Another possibility is in terms of the correlation functions or moments of this distribution, i.e.
and similar for J V and cross-terms. For fluctuations that are approximately Gaussian one can map the two characterizations to each other; the properties of a Gaussian distribution are fixed uniquely in terms of its expectation values and two-point correlation functions. In that case our description amounts to splitting and statistical Gaussian noise terms that are characterized in terms of the correlation functions 14) It is clear that these objects depend on the details of the ensemble of events considered, for example collision energy, centrality and so on.
Let us now specialize our considerations to a situation with Bjorken boost and azimuthal rotation symmetry. For the fluid dynamic fields (enthalpy density w = + p, fluid velocity u μ , shear stress π μν and bulk viscous pressure π bulk ) this implies that they can depend only on Bjorken time τ and radius r (but not on rapidity η and azimuthal angle φ). For the fluid velocity only the components u τ and u r can be non-zero and similar for the shear stress tensor π μν . For the ensemble of sources J S and J V we assume that Bjorken boost and azimuthal rotation invariance are realized in a statistical sense. For the expectation values in (2.13) this has the same implications as for the hydrodynamical fields. For the correlation functions in (2.14) the situation is more complicated since they can depend also on the differences in rapidity and azimuthal angle between the two space-time points.
We concentrate now on the fluid dynamic equations for an averaged situation where the source term on the right hand side of (2.5) is replaced by the expectation valueJ S = J S . It reads Similarly, the right hand side of Eq.
The components u φ and u η vanish due to symmetries and as a result of the constraint u μ u μ = −1 only one of the remaining equations (say for u r ) is independent. It reads Equations (2.15) and (2.16) determine the time evolution of energy density (τ, r ) and (the radial component of) the fluid velocity, respectively. They depend on the shear stress components π μν and the bulk viscous pressure. We neglect here the latter, while the former is determined from a time evolution equation that supplements the first line of Eq. (2.6) with a relaxation time term ∼ τ Shear . The shear viscosity η and the relaxation time τ Shear are chosen for concreteness according to their AdS/CFT values η = s/(4π), τ Shear = (2 − ln 2)/(2π T ). Somewhat larger values as they are presumably more realistic for QCD would not change our findings substantially. To solve Eqs. (2.15) and (2.16) one also needs a thermodynamic equation of state that relates pressure and energy density to the temperature. We take the parametrization s95p-PCE of Ref. [34]. For the initial conditions at time τ = τ 0 = 0.6 fm we follow Ref. [35] in assuming that the radial fluid velocity vanishes, u r = 0, that the shear stress assumes its Navier-Stokes value, and that the transverse energy density is determined by the nuclear overlap function for central collisions with maximal temperature T = 485 MeV.
Once the differential equations (2.15) and (2.16) are solved, one can use the equation of state to extract the temperature. The result is shown as a function of radius for different times τ in Fig. 1. The solid lines give the result without source terms, i.e. forJ S =J V = 0, while the dashed lines correspond to the full result with expectation values for sources for a 'typical' event calculated as described in Sect. 4.
In Fig. 2 we plot the radial fluid velocity as a function of radius in a similar way.
One observes that the source termsJ S andJ V have two effects: one is a slight increase in temperature at small radii and at early times, which is an expected effect from dissipation. The other is an increase of radial flow at intermediate times, the jets drag the fluid outwards. The slight decrease of the temperature in the center and the increase at large radii at later times are also a consequence of the larger flow. Both effects are relatively small. The change in the (averaged) temperature evolution seems to be negligible for practical purposes. The effect of the additional dissipated energy is hardly visible in Fig. 1, while the larger radial flow leads at larger radii to an increase in temperature of a few percent.
For the radial component of the fluid velocity this effect is more direct and leads to an increase up to about 10 %. In our setup the expectation valuesJ S andJ V are by construction symmetric under azimuthal rotations and can therefore not contribute to the harmonic flow coefficients v m . The effect of energy and momentum transfer from jets to the medium on these observables is encoded in correlation functions as in Eq. (2.14). We plan to investigate this more quantitatively in a separate publication.

Jet quenching in a hydrodynamic background
Let us now describe the formalism we use for the description of jets. Jets are simulated in Jewel [22] with the hydrodynamic calculation presented in Sect. 2 as background.
Geometrical aspects of the nucleus-nucleus collision are modeled using a Glauber model [23] with a Woods-Saxon nuclear potential. Once the impact parameter b chosen according to the geometrical cross section is fixed, the nuclear overlap can be computed Here, n A denotes the nuclear potential, z is chosen along the beam axis and r and the impact parameter b are orthogonal to this direction. The mean number of di-jets in the event is then given by where σ di-jet ( p ⊥,cut ) is the cross section per nucleon-nucleon collision for the production of a di-jet with p ⊥,jet > p ⊥,cut .
In Jewel the jet production matrix elements and initial state parton showers are simulated by Pythia 6.4 [24], which provides the leading order di-jet cross section with the EPS09 nuclear PDF set [25]. The number of di-jets per event is Poisson distributed. The jets are produced at z = 0 due to the strong Lorentz contraction of the colliding nuclei along the beam direction. In the transverse plane they are distributed according to the density of binary nucleon-nucleon collisions. In this setup jets from different nucleon-nucleon interactions are uncorrelated. The QCD evolution of the jets and re-scattering in the background take place on comparable time scales and are described in a common framework. It is assumed that the interactions of a jet resolve quasi-free partons in the medium and an infra-red continued version of the perturbative matrix elements can be used to describe them. Jewel uses leading order (LO) 2 → 2 matrix elements to describe the rescattering of hard partons in the medium and generates radiative corrections with the parton shower. The cross section for the re-scattering of a hard parton of type i with energy E in a background of temperature T and fluid velocity u μ is then given by where the partonic PDFs f i j (x, Q 2 ) encode possible initial state radiation off the energetic parton. 1 Keeping the leading terms only and introducing the infra-red regulator μ D ≈ 3T the partonic cross section reduces to where C j is the appropriate color factor. By adding the parton shower a systematic approximation to higher order 2 → n matrix elements is constructed. In this way both elastic and inelastic scattering processes are generated with the (leading log) correct relative rates. The parton shower thus generates all emissions-those associated to the QCD evolution of the jet (which would also take place in the absence of the background) and those initiated by re-scattering. In fact, it is generally impossible to assign an emission to a particular scattering process. The interplay of competing sources of radiation as well as the LPM interference are governed by the formation times of the emissions. When two emissions take place at the time the one with the shorter formation is formed as an independent particle and all scattering process within the formation time of an emission act coherently.
The local scattering rate is given by the product of the parton density and the scattering cross section (3.3) (taking care of the color factors for different parton species). When a scattering takes place a scattering center is generated from the local thermal distribution and the scattering process is simulated explicitly. The scattering centers are dynamical and recoil against the hard parton. This allows one to keep track of the energy and momentum exchange between the jet and the background.
When Jewel runs with the hydrodynamic background described in Sect. 2 it takes the temperature T (x) and transverse fluid rapidity β(x) (related to the radial component of the fluid velocity by u r = sinh β(x)) as input. The parton densities and momentum distributions are then computed assuming an ideal gas equation of state.
With the same parameter settings as used in [26] with the simple Bjorken-type background a very reasonable agreement with the jet-quenching data is found. Figure 3 shows as examples the nuclear modification factor of jets (R AA = (dσ jet /d p ⊥ )) and the di-jet asymmetry (A J = ( p ⊥,1 − p ⊥,2 )/( p ⊥,1 + p ⊥,2 )).

Characterizing the source term
The four-momentum exchange between the jet and the background can be regarded as source term in the hydrodynamic evolution. In Jewel an event contains a discrete set of scatterings with momentum transfers p μ i . The source term can thus be written as (4.1) We are now coming back to condition (2.7), namely that the energy transfer has to be positive (i.e. from the jets to the background) in order not to violate the second law of thermodynamics. The phenomenologically relevant cases concern the propagation of a jet, which is by construction harder than the thermal background and will thus on average lose energy and momentum so that the source term is positive. In principle, the framework also allows one to consider the propagation of very soft (compared to the thermal momentum scale) partons in the background. In this case the parton would on average gain energy through elastic scattering in the background (soft partons being unable to emit resolvable radiation) and the energy transfer becomes negative. This does not automatically violate the entropy condition (2.7), as the contribution of the jet or soft parton to the total entropy has to be taken into account as well. This is non-trivial since these partons are out of equilibrium, but it is a generic expectation that in both cases (hard or soft partons) their (out-of-equilibrium) entropy increases.
For the most central collisions (b = 0) the averaged source term J μ is azimuthally symmetric. It is, however, not boost invariant, since the jet production cross section is rapidity dependent 2 and the energy loss itself can in general also be rapidity dependent. For simplicity, we extract the source term only in the central unit of rapidity, where it varies only mildly, to preserve the symmetry of the background (the extension to a non-trivial rapidity dependence is straightforward). Consequently, J μ depends only on τ and r and not on φ and η. The projections of J μ parallel and orthogonal to the fluid velocity are computed using for u μ the solution to the hydrodynamic equations without the source term. This is a good approximation as long as the source term is small, when it is not small the procedure may be iterated using the new solution. Concerning the source term from jet energy loss, we assume the nucleon-nucleon collisions in a nucleus-nucleus event to be independent. Then the expectation valuesJ S and J μ V scale trivially with the number of di-jets. The results presented in this section are averaged over φ and the central unit in η.
We study Pb+Pb collisions at √ s NN = 2.76 TeV in two scenarios, namely (1) minimum-bias collisions representing a typical event and (2) events containing a O(100) GeV jet representing events in which a hard jet was triggered. In the former a major difficulty consists in defining the jet population. The perturbative jet cross section is infra-red diver- dN/dp T [GeV pert. mode with p t,cut = 3 GeV soft QCD mode thermal component Fig. 4 Distribution of jets entering the final state parton evolution in the perturbative and soft QCD mode compared to the thermal parton population with T = 485 MeV in the central unit of rapidity gent and has to be regularized, e.g. by a p ⊥ cut-off. Very low p ⊥ 'jets' should not be included anyway, because particles with momenta close to the thermal momentum at early times should be considered part of the background. The jet population should thus contain all jets that are harder than the thermal background. The p ⊥ -cut has to be placed on the jet production matrix element. In the distribution of jets entering the final state evolution it appears smeared out due to momentum re-distribution in the parton shower. With a p ⊥ -cut of p ⊥,cut = 3 GeV the initial jet distribution crosses the thermal distribution at its maximum, so that it populates predominantly the p ⊥ region where the jet population is larger than the thermal one (Fig. 4). With this value of the cut-off the leading order di-jet cross section is σ di-jet = 52.5 mb, which nearly saturates the inelastic proton-proton cross section of σ inel = (62.8 +2.4 −4.0 ± 1.2) mb [32]. It can thus be expected that in this regime some mechanism sets in that unitarizes the cross section (e.g. multi-parton interactions). This is not quantitatively under control and we do not attempt not model it here, but it adds to the already large uncertainty of the leading order di-jet cross section. In Pythia 6.4 a soft QCD mode is available, which regularizes the jet cross section in a different way [33] and also leads to a lower di-jet cross section of 41.3 mb. This scenario is studied as an alternative to our default setup to get an estimate of the uncertainties. The initial jet distribution (before final state parton showering) in the perturbative (default) and soft QCD scenario are compared to the thermal parton distribution in Fig. 4.
In the soft QCD mode the average jet p ⊥ is smaller than in the perturbative mode. The energy and momentum deposited in the medium per jet is thus lower in the former. In addition, the number of jets per event is also smaller due to the smaller cross section. Figure 5 shows the event average of the scalar sourceJ S as it enters the evolution equation of the fluid energy den- ). The physical effect is an increase of the thermal energy due to the dissipation of energy from hard particles or jets. Note thatJ S is typically maximal at small radii where the medium is densest. We compare again the perturbative and soft QCD modes. The results differ by a factor ∼1.8, but the shapes are very similar. The number of di-jets per event in the perturbative and soft QCD mode differs only by 30 % (1705 and 1342, respectively). The remaining difference between the source terms is due to the different p ⊥ . Initially,J S scales roughly as T 2 (τ, r ) · N coll (r ). This can be understood as follows: In the absence of strong flow (which is the case at mid-rapidity) the average energy loss per scattering is proportional to the temperature T . In Jewel, the scattering cross section in also temperature dependent since the infra-red regulator depends on temperature (μ D ≈ 3T ). Together with the density of scattering centers, n ∼ T 3 , this leads to a linear temperature dependence of the scattering rate per hard parton. The factor N coll (r ) comes from the initial distribution of jets. In practice, the source termJ S falls off faster with τ than T 2 (τ, r ): The average jet p ⊥ is only a few GeV so that the jet partons evolve through splitting and scattering quickly to nearly thermal momenta and do not contribute to the source terms any more. Figure 6 shows the non-vanishing components of the event-averaged vector source termJ μ V that enters the evolution equation for the fluid velocity (2.16). The two compo-nentsJ τ V andJ r V are related by u μJ μ V = 0 andJ φ V as well asJ η V vanish by symmetry reasons. The physical effect of J r V is an acceleration of the fluid in radial direction by the force opposing drag. By symmetry reasons this acceleration vanishes for r → 0, which explains the linear dependence for small r as seen in Fig. 6. While the difference between the perturbative and soft QCD modes is similar toJ S (shown in Fig. 5), the time dependence is more complicated due to the interplay of the time evolution in the fluid sector and jet In the hard di-jet scenario a cut of p ⊥,cut = 100 GeV is placed on the matrix element. The final jet population looks very different due to quenching of the jets. When comparing to experimental data one would have to place the cut on the final jet p ⊥ , which is straightforward but not necessary for this exploratory study. In Fig. 7 the momentum deposition of such a hard di-jet is compared to a minimum-bias di-jet. As expected, the source term of hard jets is much larger in magnitude and extends to significantly later times. The energy transferJ S from hard di-jets follows the approximate scaling with T 2 (τ, r ) · N coll (r ) for much longer since they do not reach thermal scales quickly. Deviations come from the dilution of the N coll (r ) distribution due to the propagation of the jets and an increase in the multiplicity of hard partons due to splitting. The same effects are also at work inJ V .
To obtain the source term for the entire event one has to add to the contribution of the hard di-jet N di-jet − 1 times that of a minimum-bias jet. A O(100) GeV di-jet deposits roughly a factor of 40 more energy and momentum than a minimum-bias jet. Since the number of minimum-bias di-jets per event is of the order 1500, the presence of a hard di-jet increases the energy transfer per event only by about 2-3 %.
The influence of the expectation valuesJ S andJ V in the perturbative minimum-bias scenario onto the fluid dynamic evolution are discussed in Sect. 2, see in particular Figs. 1 and 2. As the effect is already small in this case, the contribution of an additional hard di-jet is negligible for all practical purposes.
For the correlation functions of source terms as defined in Eq. (2.14) it is convenient to take the event average and subtract the disconnected contribution, In this form the correlators also scale linearly with N di-jet . They depend on τ 1 , τ 2 , r 1 , r 2 , φ and η = |η 1 − η 2 |. We again average over η in the rapidity window under consideration, but keep the dependence on φ. The scalar corre-latorC SS is shown in Fig. 8 for the perturbative minimumbias scenario. It is positive everywhere (except for statistical fluctuations) and strongly peaked at x = y. The fact that the correlation functions in Eq. (4.4) decay quickly with the separation between the arguments x − y implies that the corresponding fluctuations that influence the fluid dynamics are essentially local. Similarly to the averaged source terms, the correlation functions are largest at small times τ and radii r . They extend to radii of about 6 fm whereas the decay with time is very fast. The minimum-bias and hard di-jet samples look similar, for the latter the correlation functions are slightly broader. Figure  The correlators generate potentially sizeable contributions to correlation observables on the fluid dynamic side such as the anisotropic flow coefficients v n . The calculations are somewhat more involved than for the averaged source terms and will be discussed in a separate publication.

Conclusions
We have studied the influence of the energy deposition by jets onto the evolution of the medium by combing a real- istic microscopic jet quenching model with a fluid dynamic description of the bulk. The energy-momentum transfer from jets constitutes source terms in the hydrodynamic evolution equations, which we characterize in terms of event averages and correlation functions. The event-averaged source function that enters the time evolution of the energy density is largest at early times and for small radii. It leads to an increase of temperature due to the additional dissipated energy but the effect is numerically very small. The expectation value for the source function that enters the evolution of the radial fluid velocity peaks at intermediate times (a few fm/c) and for large radii (about 6 fm/c). This term gives the effect of the force opposing drag and leads to an increase of radial flow of up to about 10 %: The jets drag the fluid outwards.
The momentum transfer causing the increase in radial flow was shown to have a non-trivial functional form, which is not easily captured by simple parametrizations of the source term. This highlights the advantage of constructing a realistic source term using a model based on microscopic dynamics.
Our formalism allows one to study also event-by-event fluctuations in the source terms. Here we quantify connected two-point correlation functions and find that they are largest at early times and that they are local (the correlation functions peak strongly for equal space-time arguments).
A conceptual difficulty of a formalism that combines a microscopic description of jets with a macroscopic description of the medium is that the separation between the two components is to some extent arbitrary. This becomes appar- The η dependence is averaged over ent in the difficulties related to defining the jet population. We chose to regularize the perturbative jet cross section such that jets are produced predominantly in the phase space region where they dominate over the thermal distribution. This leads to a rather low p ⊥ cut-off of the order a few GeV. In this region the perturbative cross section has large uncertainties and multi-parton interactions may play a role. These difficulties are extenuated to some extent by the fact that-due to the fact that we employ a dynamical model of jet quenchingvery soft partons on average do not lose energy and thus do not contribute to the fluid dynamic source terms. Nevertheless, we estimate the resulting uncertainties for the latter to be of the order of a factor 2 or 3. We studied the source terms generated in a 'typical' event, i.e. without cuts on the jets, and the effect of a high-p ⊥ (O(100) GeV) di-jet. In events containing a hard di-jet the energy and momentum deposition is increased by only a few percent as compared to the minimum-bias scenario. The presence of a hard jet is thus negligible for global observables. This can be different for correlation observables such as harmonic flow coefficients, which receive potentially sizable contributions from jets. These will be studied in an upcoming publication.