Full jet in quark-gluon plasma with hydrodynamic medium response
Yasuki Tachibana Ning-Bo Chang Guang-You Qin



yasuki.tachibana@mail.ccnu.edu.cn
Institute of Particle Physics and Key Laboratory of Quark and Lepton Physics (MOE), Central China Normal University, Wuhan, Hubei, 430079, China

changnb@mail.ccnu.edu.cn
Institute of Theoretical Physics, Xinyang Normal University, Xinyang, Henan 464000, China
Institute of Particle Physics and Key Laboratory of Quark and Lepton Physics (MOE), Central China Normal University, Wuhan, Hubei, 430079, China

guangyou.qin@mail.ccnu.edu.cn
Institute of Particle Physics and Key Laboratory of Quark and Lepton Physics (MOE), Central China Normal University, Wuhan, Hubei, 430079, China






We study the nuclear modifications of full jets and their structures in relativistic heavy-ion collisions including the effect of hydrodynamic medium response to jet quenching.
To study the evolutions of the full jet shower and the traversed medium with energy and momentum exchanges between them,
we formulate a coupled jet-fluid model consisting of a set of jet transport equations and relativistic hydrodynamics equations with source terms.
In our model, the full jet shower interacts with the medium and gets modified via collisional and radiative processes during the propagation.
Meanwhile, the energy and momentum are deposited from the jet shower to the medium and then evolve with the medium hydrodynamically.
The full jet defined by a cone size in the final state includes the jet shower and the particles produced from jet-induced flow.
We apply our model to calculate the full jet energy loss and the nuclear modifications of jet rate and shape in Pb+Pb collisions at 2.76 A TeV.
It is found that the inclusion of jet-induced flow contribution leads to stronger jet-cone size dependence for jet energy loss and jet suppression.
Jet-induced flow also has a significant contribution to jet shape function and dominates at large angles away from the jet axis.



Full jet in quark-gluon plasma with hydrodynamic medium response
    Guang-You Qin
    December 30, 2023
================================================================





§ INTRODUCTION

In relativistic heavy-ion collisions at the Relativistic Heavy Ion Collider (RHIC) and the Large Hadron Collider (LHC), the deconfined state of quarks and gluons, namely quark-gluon plasma (QGP), has been created.
One of the remarkable properties of the produced QGP is the strong collective motion which has been well described by relativistic hydrodynamics with the extremely small viscosity to entropy ratio <cit.>.
The hydrodynamic behavior of the QGP has been confirmed by the large anisotropic flows as measured by experiments, which implies strong interactions among the QGP constituents.

In addition to the collective phenomena, one can study the novel features of the QGP as strongly interacting matter via jet quenching <cit.>.
In high-energy nucleus-nucleus collisions, the energetic jet showering partons may be produced in hard processes at very early time.
During their propagation through the QGP, the jet partons interact with the medium constituents via collisional and radiative processes, which change the momenta of the shower partons in the hard jet. The phenomena related to the modification of the jets and their inner structure by the QGP medium effect are commonly referred to as jet quenching in a broad sense.

In recent experiments of Pb+Pb collisions at the LHC, detailed measurements of fully reconstructed jets with very high transverse momentum have become available, thanks to the high center of mass energy.
One of the most general consequence of jet quenching seen in the experiments is the suppression of jet rates due to the energy loss of the full jets with a given cone of the size R = √(Δη_p^2 + Δϕ_p^2) <cit.>.
Furthermore, the measurements of substructures inside the full jets have provided us more detailed information on jet quenching <cit.>.
Motivated by the detailed measurements, a lot of theoretical effort has been devoted to the study of jet-medium interaction and the medium effect on the full jets with showering structures in relativistic heavy-ion collisions
<cit.>


The interactions between the energetic partons in the jet shower and the QGP medium constituents exchange energies and momenta between them.
For the QGP fluid, the jet shower is a bunch of fast-moving energy-momentum deposition sources, which are supposed to excite the medium fluid, and induce flows propagating with the jet as a collective medium response <cit.>.
The jet-induced flow carries the energy and momentum deposited by the jet and enhances the hadron emission from the medium around the direction of the jet axis.
Some of these enhanced hadrons are detected as part of the jet together with the fragments from the jet shower and affect the final state full jet energy and structures <cit.>.
Since the structure of the jet-induced flow is characterized by the bulk properties of the QGP, e.g., sound velocity, and viscosities, the detailed investigations of the contribution of the collective medium response to the jet structure may provide us not only the precise interpretation of the experimental data, but also unique opportunities to study the bulk properties of the QGP through the jet events in relativistic heavy-ion collisions.

In this work, we study the nuclear modification of full jet structures in the QGP medium including the contributions from the hydrodynamic medium response to jet-deposited energy and momentum.
We employ a coupled full jet shower and QGP fluid model which is composed of a set of transport equations for the jet shower evolution and the hydrodynamic equations with source terms for the QGP medium evolution.
The transport equations describe the evolution of the three-dimensional momentum distributions of partons in the jet shower <cit.>, including the collisional energy loss, the transverse momentum broadening, and medium-induced partonic splittings for all partons within the showering jet.
The space-time evolution of the QGP medium is described by (3+1)-dimensional relativistic ideal hydrodynamic equations with source terms <cit.>.
The source terms account for the transfer of the deposited energy and momentum from the jet shower to the QGP fluid, and are constructed with the evolving distributions of the partons in the jet shower
obtained as solutions of the jet transport equations.
Based on our coupled jet-fluid model, we perform simulations of jet events in Pb+Pb collisions at 2.76A TeV, investigate the flow induced in the medium as a hydrodynamic response to jet quenching, and study how the jet-induced flow affects the full jet structures in the final state.
Our study shows that the contribution of the particles originating from the jet-induced flows increases the jet-cone size dependence of the full jet energy loss.
It is also found that jet-induced flow has a significant contribution to the final state full jet shape, and dominates jet shape function at very large angles away from the jet direction.


The paper is organized as follows. In Sec. <ref>, we present the formulation of our coupled full jet shower and QGP fluid model used in this work.
In Sec. <ref>, we present and discuss the results from the simulations of jet events in Pb+Pb collisions at 2.76A TeV.
We will focus on the effect of the hydrodynamic medium response to the jet quenching on the final state full jet observables. Section <ref> is devoted to the summary and concluding remarks of this work.







§ THE COUPLED JET-FLUID MODEL




 §.§ Jet Shower Evolution in Medium


Jets are collimated clusters of the particles originating from high-p_T partons produced in early stage partonic hard scatterings.
The produced high-p_T parton successively radiates partons and develops a shower of partons due to its high virtuality.
In relativistic heavy-ion collisions, the interactions between the propagating jet and the constituents of the QGP medium, including collsional and radiative processes, change
the momenta of the jet shower partons, and thus modify the energy
as well as the structure of the full jet defined by a cone size.

In this work, to describe the time-evolution of the jet shower structure driven by the interaction with the medium, we employ a set of coupled transport equations for the energy and transverse momentum distributions of the partons contained in the jet shower,

    f_i(ω_i, k_i⊥^2,t)=dN_i(ω_i, k_i⊥^2,t)/dω_i dk_i⊥^2,

where the index i denotes the parton species (quark or gluon), ω_i is its energy, and k_i⊥ is its transverse momentum with respect to the jet axis.
The transport equations have the following generic form <cit.>:

    d/dtf_j(ω_j, k_j⊥^2, t) = (ê_j ∂/∂ω_j
      + 1/4q̂_j ∇_k_⊥^2)f_j(ω_j, k_j⊥^2, t)    
          +∑_i∫ dω_idk_i⊥^2 dΓ̃_i→ j(ω_j, k_j⊥^2|ω_i, k_i⊥^2)/dω_j dk^2_j⊥dt f_i(ω_i, k_i⊥^2, t)
          -∑_i ∫ dω_idk_i⊥^2 dΓ̃_j→ i(ω_i, k_i⊥^2|ω_j, k_j⊥^2)/dω_i dk^2_i⊥dtf_j(ω_j, k_j⊥^2, t).

The first two terms on the right hand side in Eq. (<ref>) describe the energy-momentum changes of the jet partons via
scatterings with the medium constituents.
The first term accounts for the energy-momentum changes in the longitudinal direction with respect to the jet axis; such contribution is called collisional energy loss and its strength is determined by the longitudinal momentum loss rate ê=d⟨ E⟩/dt.
The momentum changes in the transverse direction, whose effect is called transverse momentum broadening, is taken into account through the second term, with the exchange rate of transverse momentum squared q̂=d⟨Δ p_⊥^2⟩/dt <cit.>.
The last two terms are the contributions from the medium-induced partonic splitting processes, including the gain term coming from the radiation of parton j with ω_j and k_j⊥ from the parton i with ω_i and k_i⊥ and the loss term for the radiation of parton i from parton j.
It should be noted that the transport equations (<ref>) are coupled to each other through the terms for the medium-induced parton splittings, i.e., the gain term for the process i→ j in the transport equation for parton j also appears as a loss term in the transport equation for parton i.
The medium-induced parton splittings change the number of partons and the internal momentum distribution of the jet, but do not change the total energy and momentum inside the full jet (for sufficiently large cone size) because they are conserved in the splitting processes.

For the rates of the medium-induced partonic splittings, we employ the results from higher-twist jet energy loss formalism <cit.>:

    dΓ_i→ j(ω_j, k_j⊥^2|E_i)/dω_j dk_j⊥^2 dt = 2 α_s/πq̂_g  x P_i→ j(x)  /ω_j k_j⊥^4sin^2 (t - t_i/2τ_f).

Here, P_i→ j(x=ω_j/E_i) is the vacuum splitting function for the process i→ j with ω_j the energy of the radiated parton and E_i the energy of the parent parton, τ_f = 2E_ix(1-x)/k_j⊥^2 is the formation time of the radiated parton with k_j⊥ the transverse momentum with respect to the propagation direction of the parent parton, and t_i is the production time of the parent parton.
To obtain the splitting rates of the form dΓ̃_i→ j(ω_j, k_j⊥^2|ω_i, k_i⊥^2)/dω_j d^2k_j⊥dt
in Eq. (<ref>), the multiplication with the Jacobian 𝒥 = |∂ k_ij⊥^2/∂ k_j⊥^2| is necessary (see Ref. <cit.> for more details).
In the calculation, we put a constraint that the medium-induced radiations are allowed only for the partons whose formation times are achieved (t>τ_f).

In our model, the strength of the each contribution of the medium modification in Eq. (<ref>) is controlled by two transport coefficients: ê for collisional energy loss,
and q̂ for transverse momentum broadening and medium-induced partonic splittings.
Assuming that the medium is in almost local thermal equilibrium, the fluctuation-dissipation theorem can relate these two transport coefficients to each other: q̂ = 4T ê <cit.>.
Also since q̂ for quarks and gluons are related by a color factor, q̂_g/q̂_q = C_ A/C_ F, only one transport coefficient (we choose q̂_q for quark) determines the sizes of all the medium effects in Eq. (<ref>).
In this work, we use q̂_q of the form,

    q̂_q (x) = q̂_q,0[T(x)/T_0]^3  p· u(x)/p_0.

where T_0=T(τ=τ_0,x=0,y=0,η_ s=0) is the initial temperature at the center of the QGP medium produced in central Pb+Pb collisions at 2.76A TeV at the LHC, p^μ is the four-momentum of the propagating parton, u_μ is the flow four-velocity of the medium, and q̂_q,0 is the exchange rate of transverse momentum squared for quark in a static QGP medium at T_0. The factor p· u/p_0 is introduced to account for the flow effect in a non-static medium <cit.>.

We solve the transport equations (<ref>) numerically to describe the jet evolution in terms of the energy and transverse momentum distributions of the partons contained in the jet shower.
The initial condition for full jet shower is constructed from  simulation <cit.>. 
The jet shower evolves according to the transport equations 
during the propagation through the QGP medium with the local temperature 
above T_c = 160  MeV.
We introduce a minimum energy ω_ cut for the partons in the jet shower.
When the energy of the partons in the jet shower becomes below ω_ cut during their evolution, they are considered to be absorbed in the medium.
We put the energy and momentum of the absorbed partons into the medium together with the collisional energy loss and transverse momentum broadening terms in Eq. (<ref>).
The same cut-off energy is also used for the medium-induced partonic splittings, i.e, only partons with energy above ω_ cut can be radiated; this is to take into account some effect for the balance between radiation and absorption.
In this study, we set q̂_q,0=1.7  GeV^2/fm and ω_ cut=1  GeV. These values provide a good description to the nuclear modification factor R_ AA for inclusive jet spectrum in most central Pb+Pb collisions at 2.76 A TeV measured by the ATLAS, ALICE, and CMS Collaborations <cit.>.
Our value of q̂_q,0 is also consistent with the one obtained by JET Collaboration <cit.>.




 §.§ Hydrodynamic Equations with Source Terms


The conventional way to describe the space-time evolution of the QGP is to use the hydrodynamic equations ∂_μ T_ QGP^μν=0,
where T_ QGP^μν is the energy-momentum tensor of the QGP fluid.
This equation represents the energy-momentum conservation only in the fluid.
However, in the case that jets propagate through the QGP, the QGP and the jets exchange their energies and momenta through the scatterings between their constituents.
Therefore, the energy-momentum conservation is satisfied not for the QGP only, but for the combined system of the QGP and the jets:

    ∂_μ[T_ QGP^μν(x)+T_ jet^μν(x)]=0.

Here
T_ jet^μν is the energy-momentum tensor of the jet shower.
Assuming that the energy and momentum deposited by the jets are quickly thermalized, we model the QGP as an ideal fluid in local equilibrium whose energy-momentum tensor can be decomposed as:

    T^μν_ QGP=(ϵ+p)u^μu^ν-pη^μν,

where ϵ is the energy density, p is the pressure, u^μ is the flow four-velocity, and η^μν= diag(1,-1,-1,-1) is the Minkowski metric.
If we define the source terms by

    J^ν(x)=-∂_μT_ jet^μν(x),

Equation (<ref>) becomes the hydrodynamic equations with source terms:

    ∂_μT_ QGP^μν(x)   =   J^ν(x).

In the numerical hydrodynamic calculations for the medium evolution in the relativistic heavy-ion collisions, we use the relativistic τ-η_ s coordinates which is convenient to describe the longitudinal dynamics.
To do so, we employ the equations in which the partial derivative in the Eq. (<ref>) are the replaced with the covariant derivative in the relativistic τ-η_ s coordinates:

    D_μ̅T_ QGP^μ̅ν̅   =   J^ν̅,

where μ̅ and ν̅ are the suffixes for the components in the τ-η_ s coordinate system: x^μ̅=(τ,x,y,η_ s).
Equation (<ref>) is numerically solved to describe the space-time evolution of the QGP medium.
Using this framework, the collective response to the jet quenching in the expanding QGP fluid can be properly described.

To close the system of hydrodynamic equations (<ref>), an equation of state is necessary.
Here we employ the equation of state from the lattice QCD calculation <cit.>.
As the fluid expands, the QGP fluid cools down and finally turns into a hadronic matter according to the equation of state.
We keep using Eq. (<ref>) to describe the space-time evolution of the hadronic matter until the temperature drops to the freeze-out temperature T_ FO.



 §.§ The Source Terms


The source terms in Eq. (<ref>) describe the transfer of the energy-momentum between the jet and the QGP fluid.
Here we construct the source terms from the evolving parton distributions in the jet shower obtained as the solutions of the transport equation (<ref>).
We recall the kinetic definition of the energy-momentum tensor of the jet shower <cit.>:

    T^μν_ jet(x)
         =   ∑_j∫d^3k_j/ω_j k_j^ν k_j^μ f_j([k]_j,[x],t),

where f_j([k]_j,[x],t) is the phase-space distribution of the parton j in the jet shower.
Referring to Eq. (<ref>), the source term can be written as:

    J^ν(x)
     =
     -dP_ jet^ν/dt d^3x
     =
      -∑_j∫d^3k_j/ω_j k_j^ν k_j^μ∂_μ f_j([k]_j,[x],t),

where P_ jet^ν is the total four-momentum of the jet shower
and dP_ jet^ν/d^3x is its 3-dimensional space density.
As mentioned above, there are three types of processes, i.e., the collisional energy loss, the transverse momentum broadening, and the medium-induced partonic splitting, contributing to the derivative of the phase-space parton distribution in Eq. (<ref>).
However, since the total energy and momentum in the jet shower are conserved in the splitting processes, their contribution turns out to vanish by the integration and summation in Eq. (<ref>).
Therefore, the energy and momentum are exchanged between the jet and the QGP fluid through the processes of first two terms on the right hand side of Eq. (<ref>), i.e., the collisional energy loss and transverse momentum broadening.
We also have to estimate the position of each energetic parton in the full jet. For a parton with energy ω_j and momentum [k]_j at time t, its position is estimated as: [x]=[x]_0^ jet+([k]_j/ω_j)t, where [x]_0^ jet is the production point of the jet.
Then the source term is obtained as:

    J^ν(x)
         =   
      -∑_j∫d^3k_j
           k^ν_j.df_j ([k]_j,t)/d t|_ col.δ^(3)([x]-[x]_0^ jet-[k]_j/ω_jt).

Here .df_i ([k],t)/d t |_ col. is the part of the time derivative of the momentum distribution corresponding to the first and second terms on the right hand side in Eq. (<ref>):

    .df_j (
        ω,k^2_⊥,
        t)/dt|_ col.   =   (ê_j ∂/∂ω_j+1/4q̂_j ∇_k_⊥^2)f_j(ω_j, k_j⊥^2, t).

Considering the rotational symmetry for the distributions of the shower partons along the jet axis and using the collimated shower approximation (r = θ≈sinθ = k_⊥/ω), the source term can be written as:

    J^ν(x)   ≈   -
      1/2π r t^3(x^ν- x^ν_ jet,0).d E^ jet/d t dr|_ col.δ(|[x]-[x]_0^ jet|-t),

where

    .d E^ jet/d t dr|_ col.   =   ∑_j∫dω dk^2_j⊥ω_j.df_j(
               ω_j,k^2_j⊥,
               t)/dt|_ col.δ(r -k_j⊥/ω_j).

Here r is defined as r=√((η_ s-η_ jet)^2+(ϕ'-ϕ_ jet)^2),
with ϕ'=arctan[(y-y_0^ jet)/(x-x_0^ jet)] denoting the azimuthal angle with respect to the initial position of the jet center [x]_0^ jet, η_ jet and ϕ_ jet the pseudorapidity and azimuthal angle of the jet, and [x]^ jet(t) the position of the jet center at the time t.
Using Eqs. (<ref>), (<ref>) and (<ref>), we can construct the source terms from the energy and transverse momentum distributions of the partons inside the jet shower.
To obtain the source terms
in the relativistic τ-η_ s coordinates for Eq. (<ref>),
we perform Lorentz transformation:

    J^ν̅(τ, x, y, η_s)
       =   
    -dP_ jet^ν̅/τ dτ dx dy dη_s
       =   Λ^ν̅_μ
    J^μ(x)
     =
    - Λ^ν̅_μdP_ jet^μ/dt d^3x,

where

    Λ^ν̅_μ   =   [     coshη_ s            0            0    -sinhη_ s;            0            1            0            0;            0            0            1            0; -1/τsinhη_ s            0            0  1/τcoshη_ s ].







 §.§ Initial Condition of the Medium


We assume that the QGP medium is locally thermalized at the proper time τ=τ_0 (set as τ_0 = 0.6 fm/c in our calculation), and then apply relativistic hydrodynamics to describe its space-time evolution.
We set the initial entropy density distribution in the transverse plane at midrapidity η_ s=0 as:

    s_T([x]_⊥)
    	=
    	C/τ_0 [
    		(1-α)/2
    n^0.8[b]_ part([x]_⊥)
    +α n^0.8[b]_ coll([x]_⊥)
    	],

where n_ coll^0.8[b] and n_ part^0.8[b] are the number densities of nucleon-nucleon binary collisions and participating nucleons generated from the optical Glauber model with impact parameter [b], respectively.
The parameters C=41.4 and α=0.08 are chosen for Pb+Pb collision at 2.76A TeV by fitting the centrality dependence of the multiplicity at midrapidity from the ALICE Collaboration <cit.>.
In η_ s direction, we initialize the medium profile using the following function form:

    H(η_ s)   =   exp[-(|η_ s|-η_ flat/2)^2/2σ_η^2θ(|η_ s|-η_ flat/2) ].

The function H(η_s) has a shape consisting of a flat region around the mid-rapidity,
like the Bjorken scaling solution <cit.>,
and two halves of a Gaussian connected smoothly to the vacuum at both ends of the flat region.
The parameters η_ flat=3.8 and σ_η=3.2 are fitted to reproduce the pseudorapidity distribution for the multiplicity for central Pb+Pb collisions at 2.76A TeV.
Finally, the full 3-dimensional profile of the initial entropy density is given as:

    s(τ_0,[x]_⊥,η_ s)
    	   =   
    s_T([x]_⊥)
    H(η_ s).

We further assume that there is no initial flow at τ=τ_0 in the transverse direction, which means that the radial expansion of the medium is caused solely by the initial pressure gradient incorporated in the initial profile condition of Eq. (<ref>).
For the flow velocity in the longitudinal direction, the space-time rapidity component is initially set to zero: u^η_ s(τ=τ_0)=0 <cit.>.

In this work, we employ the smooth averaged initial profile of the medium obtained from the optical Glauber model, which does not include the initial geometrical fluctuation of the nucleons and their internal structures in the colliding heavy ions.
This initial state fluctuations in principle can affect both the jet evolution and the medium response to jet quenching; we would like to leave the event-by-event studies including the initial state fluctuations as a future work.



 §.§ Freeze-out


The hydrodynamic medium response to the energy and momentum deposited from the jets will induce additional flows in the QGP medium.
These jet-induced flows propagate and enhance hadron emissions in directions around the jet axes.
Some of the hadrons produced from jet-induced flows remain in the jet cone after the background subtraction and are counted as part of the full jets.
To obtain the momentum distributions of particles produced from the medium, we use the Cooper-Frye formula <cit.>:

    E_idN_i/d^3p_i
    =
    g_i/(2π)^3∫_Σp_i^μ dσ_μ(x)/exp[p_i^μu_μ(x)/T(x)]∓ 1,

where g_i is the degeneracy, ∓ corresponds to Bose or Fermi distribution for particle species i,
and Σ is the freeze-out hypersurface.
The freeze-out is chosen to occur at a fixed temperature T_ FO.
Here a typical value T_ FO=140  MeV (e.g., Ref. <cit.>) is used. 
For the sake of simplicity, when calculating the particle spectra via the Cooper-Frye formula, we consider only one species of bosons with the same mass as the charged pions.
The contributions from other hadron specicies are taken into account by replacing g_i by the effective degrees of freedom d_ eff defined as follows:

    e_ lat(T)   =   d_ eff(T)/(2π)^3∫ d^3p
    √(m_π^±^2+[p]^2)/exp[√(m_π^±^2+[p]^2)/T]-1,

where e_ lat is the energy density obtained from lattice QCD calculations <cit.>.

The contribution from the jet-induced medium flow to the full jet is obtained by removing the contribution from the background medium without jet propagation:

    ΔdN/d^3 p   =   .dN/d^3 p|_
    -
    . dN/d^3 p|_.

The contribution Δ dN/d^3 p is added to that of the jet shower calculated from the transport equations (<ref>) to obtain the final state full jet.
For the particles included in Eq. (<ref>), we impose the transverse momenta cut p^ trk, hyd_T>1  GeV/c for inclusive jet analyses, and p^ trk, hyd_T>0.5  GeV/c for dijet analyses following the measurements by CMS Collaboration <cit.>.



§ SIMULATIONS AND RESULTS


In this work, we initialize the jet production points in the transverse plane η_ s=0 according to the distributions of the binary nucleon-nucleon collisions which are calculated by using the Glauber model <cit.>.
Jet spectra and the momentum distributions of the shower partons inside the full jets are obtained via  simulation <cit.> with   package <cit.> used for the full jet reconstruction.
Hard jets are assumed to be created at τ = 0, and travel freely until the thermalization proper time of the QGP, τ=0.6  fm/c.
Then the jet shower starts to interact with the QGP and evolves according to the transport equations (<ref>).
The jet-medium interaction is turned off when the local temperature of the medium is below T_c = 160  MeV.
The interaction in the hadronic matter is usually small compared to the QGP phase, and is neglected in this work.

The medium profile at initial time τ=τ_0 is calculated by employing Eq. (<ref>) with the impact parameter [b]=0 (central collisions).
The evolution of the medium is governed by the ideal hydrodynamic equations with source terms (<ref>).
As the system expands and cools down, it transits from QGP phase to hadronic phase, and finally the freeze-out occurs.
The momentum distributions of hadrons produced from the medium are calculated via the Cooper-Frye formula (<ref>).
After the subtraction of the background (without jet), the remaining part  (<ref>) contributes to final state full jets.

Hereafter we call the part of jet described by the transport equations (<ref>) as shower part, and that coming from the jet-induced flow via freeze-out as hydro part.
In this study, we mainly focus on the contribution of the hydro part of the full jet.
Detailed studies of the contribution of each medium modification effect on the jet shower part can be found in Ref. <cit.>.
Also since we perform the simulation independently for each single jet shower, the effect of possible interference between the flows induced by multiple jet showers in one event is not included.

  




 §.§ Flow Induced by Full Jets


Figure <ref> shows the snapshots of the energy density distribution of the medium in the transverse plane at η_ s=0 at different proper times for an example event.
For this event, the single jet 
with the initial p_T^ jet = 150  GeV/c 
is produced at (x_0^ jet, y_0^ jet)=(0  fm, 6.54  fm) and travels in the direction of ϕ_p=5π/8.
The upper panels show the whole energy density of the medium, and the lower panels show the energy density after the subtraction of the energy density in the events without jet propagation.
From these figures, we can see that the V-shaped wave fronts (shown by higher energy density region) are induced by the jet propagation, and develop with time in the medium.
This V-shaped wave front is the Mach cone <cit.>, a conical shock wave that appears as an interference of sound waves caused by an object moving faster than the medium sound velocity.
Here the highly collimated jet shower deposits its energy and momentum and induces a Mach cone whose vertex is the center of the jet <cit.>.
This wave front of the Mach cone carries the energy and momentum, propagates outward and also causes the lower energy density region behind the wave front.
During the propagation, the Mach cone and the radial flow of the medium are pushed and distorted by each other.
One can see that the Mach cone is asymmetrically deformed in this example because the jet travels through the off-central path in the medium.

In this work, we neglect the effect of the finite small shear viscosity of the QGP and model the medium created in relativistic heavy-ion collisions as an ideal (non-viscous) fluid.
The finite viscosities are important for more precise description of the medium evolution and the collective anisotropic flows observed in the final states <cit.>.
It can also affect the shape of the medium response to the jet-deposited energy and momentum, e.g., the Mach cone can be smeared by the finite shear viscosity <cit.>.
In our study, we assume the instantaneous thermalization of the energy and momentum deposited by the jet; the finite relaxation time effects may be included in the source terms <cit.> (note that the smearing due to the finite grid size in the hydrodynamic simulation mimics some relaxation effect).
Since the relaxation times for the deposited energy and momentum are closely related to the transport coefficients of the QGP, and the inclusion of such effects would provide further information on the QGP's properties, which we would like to leave as a future work.


  

  



 §.§ Full Jet Energy Loss and Suppression


In our framework, the final full jets are contributed from two parts: jet shower part and hydrodynamics response part.
The shower part of the jet loses energy due to three mechanisms: the collisional energy loss and the absorption of the soft partons by the medium, the transverse momentum broadening which kicks the partons out of the jet cone, and the medium-induced radiation outside the jet cone.
The hydro part of the jet comes from the lost energy and momentum from the jet shower which thermalize into the medium and induce connical flow; some of the energy is still inside the jet cone.
Thus the hydro part will partially compensate the energy loss experienced by the jet shower part.
Here we study the effect of jet-induced medium flow on full jet energy loss and full jet suppression.


Figure <ref> shows the mean value of the total energy (transverse momentum) loss for inclusive jets with and without the hydro part contribution as a function of initial full jet transverse momentum.
The left panel shows the results for the jet-cone sizes R=0.2, 0.3, and 0.4, and the right shows for R=0.3, 0.6, and 0.9.
One can see the general feature that the amount of the energy loss increases with increasing initial jet transverse momentum while the fractional energy loss decreases.
The total p^ jet_T loss for the full jets with the inclusion of the hydro part contribution is smaller than that without the contribution from the medium response.
For jets with the cone size R=0.3, about 10 % of the lost p^ jet_T from the jet shower part is recovered by the hydro part.

We can also see the jet cone size dependence of jet energy loss from Figure <ref>.
For the shower part without the hydro part contribution, the jet cone size dependence is rather weak.
This is due to the reason that the shower part of the jet is quite collimated, i.e., most of the energy in the shower part is covered by a narrow jet cone, therefore, jet energy does not change much with increasing jet cone sizes.
On the contrary, jet-induced flow evolves with medium, diffuses, and can spread quite widely around jet axis. As a result, the jet cone size dependence becomes much stronger when adding the hydro part contribution.

The effect of full jet energy loss in the relativistic heavy-ion collisions can be quantified by the measurements of nuclear modification factor R_ AA for single inclusive jet spectrum, defined as:

    R_ AA= 1/⟨ N_ coll⟩d^2N_ jet^ AA/dη_p dp^ jet_T/d^2N_ jet^ pp/dη_p dp^ jet_T,

where ⟨ N_ coll⟩ is the number of binary nucleon-nucleon collisions averaged over events in a given centrality class, N_ jet^ AA is the number of jets in nucleus-nucleus collisions, and N_ jet^ pp is that in p+p collisions.
One important result of jet energy loss is that jet p_T spectrum in nucleus-nucleus collisions is shifted to lower p_T^ jet compared to that in p+p collisions.
Since the jet spectrum is a steeply decreasing function of p_T^ jet, jet R_ AA will become smaller than unity in high-p_T^ jet region.

Figure <ref> shows the nuclear modification factor R_ AA for single inclusive jets as a function of p_T^ jet for different jet cone sizes:
the left panel for the jet-cone sizes R=0.2, 0.3, and 0.4, and the right for R=0.3, 0.6, and 0.9.
We also compare the results with and without the inclusion of the contribution from the jet-induced flow.
We find that without the hydro part contribution, the jet cone size dependence for jet R_ AA is very week, which is consistent with the weak dependence for jet energy loss as seen in Figure <ref>.
The inclusion of the contribution from jet-induced flow decreases the total energy loss and thus increase the value of R_ AA; it also increases the jet-cone size dependence of R_ AA.
Our results are comparable with CMS measurements with the jet-cone sizes R=0.2, 0.3, and 0.4, which show relatively small jet cone size dependence (but with large error bars).


  



  



 §.§ Full Jet Shape Function


One of the advantages of studying fully reconstructed jets in relativistic heavy-ion collisions is that one may investigate not only the full jet energy loss and suppression, but also their internal structures which provide us the detailed information on how the energy is distributed inside the full jets and how the energy distribution is modified by the interaction with the QCD medium.
Jet shape function describes how the energy inside (and outside) the full jets is distributed in the radial direction (transverse to the jet axis) and is defined as follows:

    ρ_ jet(r)    =   1/N_ jet∑_ jet[
        1/p_T^ jet∑_ trk∈(r-δ r/2,r+δ r/2) p_T^ trk/δ r],

where r = √((η_p - η_ jet)^2 + (ϕ_p - ϕ_ jet)^2) is the radial distance of the jet constituents from the jet axis, δ r is the bin size,
and the sum is taken over all constituents (tracks) of the full jets in the bin at r.

The left panel of Figure <ref> shows our result for the jet shape function inside the jet cone for inclusive jets with p^ jet_T>100  GeV/c and R=0.3 in central Pb+Pb collisions and in p+p collisions, compared to the experimental data from CMS Collaboration <cit.>.
To see the medium effect on the jet shape function more clearly, the nuclear modification factor for the jet shape function R^ρ_ AA(r)= ρ_ AA(r)/ρ_ pp(r) is shown in the right panel of Figure <ref>.
We can see that our results (both with and without the contribution from the hydrodynamic response part) show similar nuclear modification pattern for the jet shape function to the experimental data from CMS Collaboration, i.e., little change for small r, a dip at r∼0.1 and an enhancement at large r.
In other words, the inner hard core of the jet is more collimated while the tail (the outer soft part) of the jet is broadened, in central Pb+Pb collisions compared to pp collisions.

The medium modification feature for the shower part of the full jet has been extensively studied in Ref. <cit.> which shows that the collisional energy loss and the thermalization of the soft shower partons (into the medium) make the jet narrower with more collimated hard core, while the transverse momentum broadening and medium-induced radiation transport the energy from the inner to the outer sides of the jet and broaden the tail of the jet shape function.
After the inclusion of the contribution from jet-induced medium flow, the jet shape function at small r is not modified much, but for large r region (r>0.2-0.25), there is a significant enhancement of the jet broadening effect.
This seems to be quite natural considering the jet cone size dependence of full jet energy loss as seen in Figure <ref>, i.e., the energy loss from the shower part of the jet induces conical flow and medium excitation which evolve with the medium and diffuse to larger angles with respect to the jet axis.


To see more clearly the contribution from the hydrodynamic response part (jet-induced medium flow) to jet broadening effect, we show in Figure <ref> the jet shape function ρ(r) for inclusive jets with an extended radial distance 0 < r<1.
The trigger p_T threshold for the inclusive full jets is set to be p_T^ jet> 100  GeV/c.
Here we still use p_T^ jet defined by the jet-cone size R=0.3 as the normalization factor for the jet shape function at r>R=0.3, to be consistent with the experimental results from CMS Collaboration <cit.>. 
The red solid line shows the result for jets with both shower and hydro parts, 
the orange dashed dotted solid line shows the contribution from the hydro part,
and the blue dashed line shows the result for jet without hydro part. 
The green dotted line shows the result from   simulation.
As we can see, the shower part of jet shape function is a deep falling function of r, while the energy (momentum) from the hydrodynamic response part is a quite flat distribution in a wide range of r.
This is because the energy loss from the shower part is carried away by the jet-induced flow which evolves with the medium and diffuses to large distances <cit.>.
Compared to   simulation, the broadening of the shower part of the jet continues to large r by the transverse momentum kicks and medium-induced radiation, but the contribution from the hydro part to the jet shape function is quite flat and finally dominates over the shower part in the region with r>0.5.



CMS Collaboration has recently measured the jet shape functions with a wide range of r (up to r=1) for both leading and subleading jets in asymmetric dijet events in Pb+Pb collisions at 2.76A TeV <cit.>.
We also perform the calculation for the jet shape functions in dijet events, and the comparison with CMS data is shown in Figure <ref>.
In the calculation, we chose dijet events with 
the leading jet with p^ jet_T,1 > 120  GeV/c, 
the subleading jet with p^ jet_T,2 > 50  GeV/c, 
and the azimuthal angle between the leading and subleading jets Δφ_1,2>5π/6.
In the figure, we do not show the jet shape function for p+p collisions from CMS <cit.> since the data contain the contamination from the underlying event (and therefore are quite different at large r region as compared to the jet shape function obtained from  simulation).
In Pb+Pb collisions, such background effect is supposed to be small at large r, since the jet shape function is dominated by the shower part at small r and by the jet-induced medium flow at large r (see below).

From Figure <ref>(a), we can see that the jet shape function for leading jets in central Pb+Pb collisions is quite similar to that for inclusive jets shown in Figure <ref>.
The shower part dominates the jet shape function at relatively small r and is broadened by the medium effect starting from r=0.2-0.3, while the hydro part starts to dominate the jet shape function at large r region (r>0.5).
Our full result on jet shape function with the contributions from both shower and hydro parts reproduces the experimental result quite well throughout the entire r range (up to r=1).
For subleading jets as shown in Figure <ref>(b), we can see that the jet shape function is much broader than that of leading jets due to larger jet-medium interaction for subleading jets.
As a result, the shower part of the jet shape function is wider and the hydro part contribution is also larger and more widely distributed than that in leading jets.
Our full result also provides a good description of the jet shape function for subleading jets except the middle r region.
The above results clearly show that the hydrodynamic medium response to jet-medium interaction plays an important role in the study of fully reconstructed jets, especially at large r region.
This means that the jet shape function at large angles with respect to the jet direction provides a good opportunity to study the hydrodynamic medium response to jet quenching.





§ SUMMARY


In this work, we have studied the nuclear modifications of full jet structures in relativistic heavy-ion collisions with the inclusion of the contribution from the medium exitations induced by the propagating hard jets.
We have formulated a model which consists of a set of transport equations to describe the full jet shower evolution in medium and relativistic ideal hydrodynamic equations with source terms to describe the dynamical evolution of the QGP medium.
The transport equations control the evolutions of energy and transverse momentum distributions of the shower partons within the full jet.
The contribution from the momentum exchange with the medium via
scatterings with medium constituents is taken into account by collisional energy loss and transverse momentum broadening terms. 
The partonic splitting terms account for the contribution of the medium-induced radiations; the rates for the induced splittings were taken from higher-twist jet energy loss formalism.
The local temperature and flow velocity of the medium are embedded in the jet quenching parameter q̂ which controls the amplitudes of all the medium modification processes.
The relativistic ideal hydrodynamic equations with source terms determine the space-time evolution of the medium which exchanges the energy and momentum with the propagating jet shower.
The energy and momentum deposited by the jet shower into the medium fluid are included via the source terms, which may be constructed from the solutions of jet shower transport equations based on the energy-momentum conservation for the combined system of the QGP medium and the jet shower.

Based on our coupled jet shower and QGP fluid model, we have performed the simulations of jet events in central Pb+Pb collisions at 2.76A TeV.
The transport equations for jet shower were numerically solved with the initial conditions generated by .
We kept track of both the energy and transverse momentum distribution of all shower partons in the full jet until the interaction with the QGP ceases.
The relativistic ideal hydrodynamic equations with source terms were numerically solved in the (3+1)-dimensional τ-η_ s coordinates, with the initial profile of the medium obtained from the optical Glauber model.
We found that the additional energy density flow can be induced by the jet shower propagation in the medium, and the jet-induced conical flow is pushed and distorted by the medium radial flow (and vice versa) during the propagation.

To study how the hydrodynamic medium response contributes to the full jet observables, we calculated the particles produced from the hydrodynamic medium response by using the Cooper-Frye formula, which are combined with the jet shower part to obtain the final state full jets.
We generated (di)jet events in central Pb+Pb collisions by using  for the momentum distribution and the Glauber model for the spatial distribution.
We calculated the total energy (p_T) loss of the full jets for different jet-cone sizes and found that the contribution of the hydro part partially compensates the energy loss of the jet shower.
Such compensation effect increases with increasing jet cone size.
As a result, one obtains stronger jet cone-size dependence for the single inclusive jet R_ AA when taking into account hydrodynamic response contribution.
The effect of jet-induced medium flow on jet shape functions were studied for inclusive jets as well as for dijets in central Pb+Pb collisions at 2.76A TeV.
Jet shape functions for leading jets in dijet events and single inclusive jets are quite similar, while the nuclear modification for subleading jets is larger due to more jet-medium interaction.
Our results showed that the particles produced from jet-induced medium flow do not affect very much the jet shape function at small r, but significantly enhance the broadening of the jet shape function, and finally dominate at large r region.
Our full results for jet shape functions can reproduce quite well the experimental data from CMS Collaboration <cit.>, after taking into account both jet shower and hydrodynamic response contributions.
In summary, we have found that jet-induced flow plays significant roles in the study of jet structure in relativistic heavy-ion collisions, especially at large angles with respect to the jet axis; detailed studies of full jet structure at large r should provide us much information about the medium response effect in the process of jet-medium interaction.



The authors would like to thank X.-N. Wang for discussions and  helpful comments.
Also, Y.T. is grateful to Y. Hirono for discussions regarding numerical implementations.
This work is supported in part by the Natural Science Foundation of China (NSFC) under Grant Nos. 11375072 and 11405066, Chinese Ministery of Science and Technology under Grant No. 2014DFG02050, and the Major State Basic Research Development Program in China under Grant No. 2014CB845404.






