when an energetic quark or gluon propagates through a dense partonic system such as the quark - gluon plasma ( qgp ) , it interacts with the hot medium through multiple scattering and induced gluon bremsstrahlung , leading to transverse momentum broadening and energy loss of the propagating parton . as a consequence , spectra of large transverse momentum hadrons and jets in high - energy heavy - ion collisions should be suppressed if a dense partonic matter is formed in the early stage of the collisions . this phenomenon is called jet quenching and has been proposed as a powerful tool for the study of properties of the qgp in high - energy heavy - ion collisions @xcite .
it has been observed in experiments at both the relativistic heavy - ion collider ( rhic ) @xcite and the large hadron collider ( lhc ) @xcite as evidence for the formation of strongly interacting qgp in high - energy heavy - ion collisions .
investigators observed not only the suppression of single inclusive hadron spectra at large transverse momentum @xcite but also back - to - back dihadron @xcite and gamma - hadron correlations @xcite .
the jet quenching phenomenon becomes even more dramatic in heavy - ion collisions at lhc when one observes the onset of dijets and gamma jets with large asymmetries in transverse energy @xcite .
these patterns of jet quenching can all be understood quantitatively within the picture of multiple scattering and parton energy loss @xcite-@xcite .    within the picture of perturbative qcd ( pqcd ) , the interaction between an energetic parton and the medium is dominated by small - angle scattering and induced gluon radiation .
the transverse momentum broadening and energy loss of the propagating parton depend on the transverse momentum transfer from the medium during each scattering which in turn encodes the interaction among medium partons at the scale as probed by the energetic parton .
the jet transport parameter defined as the averaged transverse momentum broadening squared per unit length of propagation represents a fundamental property of the medium @xcite-@xcite .
recent phenomenological studies of single - hadron suppression in high - energy heavy - ion collisions within collinear factorized pqcd models give extracted values of the transport parameter @xmath0 gev@xmath1/fm and @xmath2 gev@xmath1/fm at an initial time @xmath3 fm/@xmath4 for a 10 gev quark at the center of most central au+au collisions at @xmath5 gev and pb+pb collisions at @xmath6 tev , respectively @xcite .
these values are about two orders of magnitude higher than in large cold nuclei @xcite .    given such large values of the jet transport parameter , the energy lost to the medium by a propagating parton is also very large .
dissipation of this large amount of lost energy in the medium can lead to jet - induced medium excitations such as supersonic waves or mach cones @xcite .
study of such jet - induced medium excitations can shed light on the transport properties of the bulk medium which complements that obtained from analyses of anisotropic flows due to collective expansion .
several types of models can be employed to study jet - induced medium excitations such as hydrodynamical models @xcite , string theory @xcite and parton transport models @xcite .
the linear boltzmann transport ( lbt ) model @xcite combines a kinetic description of parton propagation and hydrodynamic description of the underlying medium evolution .
since it also keeps track of thermal recoil partons from each scattering and their further propagation in the medium , the lbt model can be used to study both jet transport in medium and jet - induced medium excitations .
jet - induced medium excitations in the form of thermal recoil partons and their further propagation are also important for the study of jet quenching through reconstructed jets and the corresponding jet profiles . reconstructed jets defined through a jet - finding algorithm @xcite in experiments are collections of collimated showers of hadrons within a jet cone . neglecting nonperturbative effects through hadronization ,
jets reconstructed through jet shower partons and final hadrons are equivalent and have approximately the same energy .
since interaction in pqcd between an energetic parton and the thermal medium is dominated by small - angle scatterings , a large fraction of thermal recoil partons from jet - medium interaction fall within the jet cone and therefore should become a part of the reconstructed jet .
these thermal recoil partons within the jet cone can contribute to a non - negligible fraction of the total jet energy and modify both the jet transverse profile ( transverse energy distribution inside the jet cone ) and fragmentation functions ( particle distributions in longitudinal momentum inside the jet cone ) @xcite .
in addition , a diffusion wake generated behind a propagating parton effectively modifies the underlying fluid background and therefore also affects jet reconstruction , the final jet energy , and the jet profile .
these effects can all be studied within the lbt model and are important for understanding the medium modified jet cross section and jet profiles @xcite .    in this paper ,
we describe in detail the lbt model that has been developed over the last few years and improvements to the early versions that have been used for the study of jet - induced medium excitations @xcite , jet transport , and medium modifications @xcite .
we also employ the model to study elastic energy loss through multiple scattering , propagation of jet - induced medium excitations , and their contributions to parton distribution and jet profiles of reconstructed jets using the anti-@xmath7 algorithm in fastjet @xcite .    instead of using the small - angle approximation for parton scatterings in the early versions of the lbt model
, we have included the complete set of @xmath8  elastic " parton scattering processes with full scattering matrix elements including large - angle scatterings , annihilation processes , and flavor conversion .
these are necessary for more accurate description of parton transport in medium .
simulations of parton transport according to the boltzmann equation in the lbt model are based on local rates of scattering in a medium that evolves according to relativistic hydrodynamic equations .
we keep track of medium recoil partons from each parton - medium scattering and neglect interactions among jet shower partons and thermal recoil partons , wherein lies the meaning of  linear " in lbt .
the model is therefore lorentz covariant and can avoid difficulties of implementing the test - particle method in a full - fledged parton transport model @xcite .
furthermore , it also allows us to study jet transport and jet - induced medium excitations based on event - by - event simulations with full fluctuations and correlations .
this article is devoted to a general description and test of elastic scatterings in the lbt model and the corresponding jet transport and jet - induced medium excitations in a uniform qgp medium .
the final medium modification of jets and the effect of jet - induced medium excitations will depend strongly on the inelastic parton - medium scatterings and the dynamical expansion of the medium . these will be described in the future when we carry out realistic studies of jet transport and medium modification in high - energy heavy - ion collisions .
the remainder of this paper is organized as follows .
we first give a detailed description of the lbt model in sec .
we compare numerical and approximate analytic results for local parton scattering rates for different channels , parton energies , and local temperatures .
differential scattering rates as functions of the final - state momentum and scattering angle of the medium recoil parton from monte carlo simulations are also tested against results calculated with numerical integrations . in sec .
iii , we study the elastic energy loss from multiple scattering and in particular the energy and time or length dependence in a static and uniform qgp medium
. time evolution of the transverse momentum distribution of the leading parton is also investigated .
propagation of medium recoil partons in terms of jet - induced medium excitations are investigated in sec .
iv . a supersonic shock wave along with a diffusion wake is shown to arise from the propagation of the energy and momentum deposited by the propagating energetic parton . in sec .
v , we investigate contributions from medium recoil partons to the reconstructed jet energy using the anti-@xmath7 algorithm in fastjet jet finding package and study the influence of medium recoil partons to the jet energy profile and parton distributions . finally , a summary and discussions are given in sec .
considering only @xmath9 parton scattering , jet transport and propagation in a qgp medium can be approximately described by a set of linear boltzmann equations , @xmath10 \ , |m_{ab\rightarrow cd}|^2   \\ & \times s_2(s , t , u)(2\pi)^4\delta^4(p_1+p_2-p_3-p_4 ) , \end{split}\ ] ] where the sum is over the flavor of the initial thermal parton ( @xmath11 ) and all possible scattering channels @xmath12 for different flavors of the final partons ( @xmath4,@xmath13 ) .
@xmath14 is the corresponding matrix element @xcite that is averaged ( summed ) over the initial ( final ) spin and color as a function of standard mandelstam variables @xmath15 , @xmath16 , and @xmath17 .
we have neglected the effect of quantum statistics in this study .
the phase - space distributions for thermal partons in the qgp medium @xmath18 @xmath19 are the bose - einstein distribution for gluons and fermi - dirac distribution for quarks and antiquarks with a local temperature @xmath20 and fluid velocity @xmath21 . for @xmath22 number of colors in qcd , the spin and color degeneracy factors are @xmath23 and @xmath24 for gluons and each flavor of quark or antiquark , respectively . the propagating parton before and after each scattering is assumed to follow a classical trajectory with @xmath25 @xmath26 as the phase - space density .
most matrix elements @xmath14 for two - parton scattering diverge at small angle , @xmath27 , for massless partons .
such divergences disappear when quasi - particle modes due to hard - thermal - loop ( htl ) corrections are taken into account in the calculation of scattering rates and transport coefficients @xcite . while an exact formulation of the boltzmann transport within htl pqcd is difficult , we introduce a screening mass for light quarks and gluons in matrix elements for two - parton scattering .
this is the same approach as in refs .
this is equivalent to introducing a lorentz - invariant regularization condition , @xmath28 in the interaction kernel in the boltzmann transport equations [ eq .
( [ lbt ] ) ] , where @xmath29 is the debye screening mass , @xmath30 for @xmath31 number of active quark flavors in the qgp .
matrix elements for parton scattering we use in this paper are leading - order pqcd results .
the strong coupling constant @xmath32 will be kept at a constant value and should be determined through comparisons to the experimental data . since we will only consider parton propagation in a uniform medium in this study , we set @xmath33 throughout the paper .      for monte carlo simulations of the boltzmann transport equation , we need to evaluate the rate for a hard parton of type @xmath34 scattering with a thermal parton of type @xmath11 via a specific channel @xmath12 , @xmath35 summing over all types of initial thermal partons ( @xmath11 ) and all possible channels with different types of final - state partons ( @xmath4,@xmath13 ) , we can obtain the total scattering rate for an energetic parton @xmath34 with thermal partons in a thermal qgp medium : @xmath36 these scattering rates for different parton types and different channels will be used in monte carlo simulations of the lbt model . shown in fig .
[ fig : rates_vs_e ] are the total scattering rates for an energetic gluon ( upper panel ) and quark ( lower panel ) as a function of the incident energy @xmath37 in a thermal qgp medium with zero baryon chemical potential at a constant temperature @xmath38 mev . also shown
are contributions from different channels .
dominant contributions to the total scattering rate of a gluon come from gluon - gluon and gluon - quark scattering .
contributions from gluon - gluon to quark - antiquark conversion are about two orders of magnitude smaller .
the dominant channel to the quark scattering rate is from quark - gluon compton scattering followed by quark - quark scattering .
contributions from quark - antiquark annihilation are also many orders of magnitude smaller .    as shown in fig .
[ fig : rates_vs_e ] , scattering rates for @xmath15-channel annihilation processes decrease with the incident energy as @xmath39 .
total parton scattering rates are dominated by @xmath16-channel gluon - gluon , gluon - quark , and quark - quark scattering processes , whose rates are approximately independent of the incident energy @xmath37 .    in the high - energy limit @xmath40 ,
@xmath16-channel gluon and quark scattering cross sections can be approximated by their small - angle limits , @xmath41 scattering rates for energetic partons in a thermal qgp medium are then @xcite , @xmath42 where @xmath43 is the apry s constant , @xmath44 and @xmath45 . with the debye screening mass given in eq .
( [ cutoff ] ) , the above scattering rates are independent of the incident energy and proportional to the temperature @xmath20 . as shown in fig .
[ fig : rates_vs_e ] as solid lines , the above analytic estimates of total parton scattering rates agree with numerical results very well for @xmath40 .
we also show in fig .
[ fig : rates_vs_t ] both numerical results ( symbols ) and analytic estimates ( solid lines ) of total scattering rates and contributions from different channels as functions of the temperature in a thermal qgp medium for a fixed parton incident energy @xmath46 gev . as given in the above analytic estimates , total rates as well as rates for @xmath16 channels increase linearly with the temperature @xmath20 , while rates for @xmath15 channels increase quadratically with @xmath20 .
note that the parton scattering rates given here are only valid in a medium above the qcd phase transition temperature @xmath47 . for a rapid crossover
as indicated by recent lattice qcd simulations @xcite , the effective degrees of freedom around the phase transition temperature @xmath48 mev become ambiguous .
one has to rely on model assumptions to treat parton medium interaction @xcite .     and total scattering rates @xmath49 of ( a ) a gluon or ( b ) a quark ( antiquark ) as a function of the parton energy @xmath50 in a static and uniform qgp medium
. analytic estimates for the total rates discussed in the text are shown as solid lines.,title="fig:",width=321 ] +   and total scattering rates @xmath49 of ( a ) a gluon or ( b ) a quark ( antiquark ) as a function of the parton energy @xmath50 in a static and uniform qgp medium .
analytic estimates for the total rates discussed in the text are shown as solid lines.,title="fig:",width=321 ] +     and the total scattering rates @xmath49 of ( a ) a gluon or ( b ) a quark ( antiquark ) as a function of the temperature @xmath20 in a static and uniform qgp medium .
analytic estimates for the total rates discussed in the text are shown as solid lines.,title="fig:",width=321 ] +   and the total scattering rates @xmath49 of ( a ) a gluon or ( b ) a quark ( antiquark ) as a function of the temperature @xmath20 in a static and uniform qgp medium . analytic estimates for the total rates discussed in the text are shown as solid lines.,title="fig:",width=321 ] +      to sample the initial thermal and final parton four - momenta in a @xmath9 scattering , we choose the @xmath51 axis as the direction of the incident parton @xmath34 with its four - momentum @xmath52 .
one can integrate out the four - dimensional @xmath53 function in eq .
( [ singlerate ] ) and express the the scattering rate for a given channel in a differential form : @xmath54 where @xmath55 @xmath56 is the azimuth angle , @xmath57 the polar angle of a parton s momentum @xmath58 , and @xmath59 and @xmath60 are the azimuth and polar angles between two partons momenta @xmath58 and @xmath61 , respectively .     or ( b ) @xmath62 scattering processes with different initial parton energy @xmath50 in a thermal qgp medium at a temperature @xmath38 mev .
solid lines are results from eq .
( [ scattrate2 ] ) with direct integration while histograms are from the corresponding monte carlo simulations.,title="fig:",width=302 ]   or ( b ) @xmath62 scattering processes with different initial parton energy @xmath50 in a thermal qgp medium at a temperature @xmath38 mev .
solid lines are results from eq .
( [ scattrate2 ] ) with direct integration while histograms are from the corresponding monte carlo simulations.,title="fig:",width=302 ]     in ( a ) @xmath63 or ( b ) @xmath64 scattering processes.,title="fig:",width=302 ]   in ( a ) @xmath63 or ( b ) @xmath64 scattering processes.,title="fig:",width=302 ]    one can further carry out the integration over the direction of the initial thermal parton @xmath11 s momentum @xmath65 and obtain the final parton @xmath4 s angular ( @xmath66 ) and energy ( @xmath67 ) distributions . shown in figs .
[ fig - e - gluon ] and [ fig - e - quark ] as solid lines are the energy distributions of an outgoing parton in some selected channels of parton - medium interaction in a thermal qgp with a temperature @xmath38 mev .
the gluon - medium interaction is dominated by @xmath16 and @xmath17-channel gluon - gluon and gluon - quark(antiquark ) scattering .
the outgoing gluon can either be the deflected incoming gluon which carries most of the incident energy @xmath50 minus elastic energy loss @xmath68 or the thermal recoil gluon whose average energy is its initial thermal energy @xmath69 plus the elastic energy transfer from the scattering .
the energy spectrum of the outgoing gluon therefore has two peaks at @xmath70 and @xmath50 as shown in the upper panel of fig .
[ fig - e - gluon ] .
large - angle and @xmath15-channel scatterings contribute to the valley region between these two peaks which is suppressed by a factor of @xmath71 .
the energy spectrum of the outgoing gluon does not vanish at the incident energy @xmath50 .
it instead has an exponential fall - off with a slope as given by the temperature , indicating that the energetic incident parton can also gain energy due to its interaction with thermal partons .
the energy spectra of the outgoing quark in the quark - antiquark pair production @xmath72 ( lower panel of fig .
[ fig - e - gluon ] ) and the outgoing gluon in the quark - antiquark annihilation @xmath73 process ( upper panel of fig .  [ fig - e - quark ] ) have a similar structure with two peaks at @xmath70 and @xmath50 .
this structure , however , arises mainly from the @xmath17 and @xmath16 channels of the pair production or annihilation process in which the two outgoing partons have an equal probability to carry most of the incident energy .
contributions from the @xmath15 channels in these processes are suppressed by a factor of @xmath74 with a flat energy distribution between @xmath75 .    to illustrate the typical energy spectra of thermal recoil partons from parton - medium interaction we show in the lower panel of fig .
[ fig - e - quark ] the energy distribution of the recoil parton @xmath11 in @xmath76 ( @xmath77 ) processes . for an energetic quark @xmath78 ,
these processes are dominated by @xmath16 channels in which recoil partons carry an energy @xmath79 .
the recoil parton spectra should have a power law behavior @xmath80 with a peak at @xmath81 .
the recoil gluon spectra in the @xmath15 and @xmath17 channels in the quark - gluon compton scattering @xmath82 , however , still has a peak at @xmath83 , but is suppressed by a factor of @xmath71 relative to the @xmath16 channel .     or ( b ) @xmath62 scattering processes with initial parton energy @xmath84 gev in a uniform medium at a temperature @xmath38 mev .
solid lines are for results from the numerical integration and histograms for lbt monte carlo simulations . , title="fig:",width=302 ]   or ( b ) @xmath62 scattering processes with initial parton energy @xmath84 gev in a uniform medium at a temperature @xmath38 mev .
solid lines are for results from the numerical integration and histograms for lbt monte carlo simulations . , title="fig:",width=302 ]     for ( a ) @xmath85 or ( b ) @xmath86 scattering processes.,title="fig:",width=302 ]    for ( a ) @xmath85 or ( b ) @xmath86 scattering processes.,title="fig:",width=302 ]    because of the dominance of @xmath16 and @xmath17 channels of the small - angle scattering , only one of the outgoing partons in parton - medium scattering carries most of the incident energy for most of the time , whose flavor can be , however , different from the incident parton .
therefore , we will refer to the more energetic parton in the final state of each parton - medium scattering as the leading parton and the other less energetic one as the thermal recoil parton .    for monte carlo simulations in the lbt model
, we change the variable of the energy integration in eq .
( [ scattrate2 ] ) from the outgoing parton @xmath4 s energy @xmath67 to the initial parton @xmath11 s thermal energy @xmath87 .
the integrand then can be used to sample the thermal parton @xmath11 s energy @xmath87 , polar angle @xmath88 and the outgoing parton @xmath4 s direction ( @xmath89 ) with the rejection method .
the energy of the outgoing parton @xmath4 @xmath67 and parton @xmath13 s four - momentum @xmath90 can be determined by energy - momentum conservation . shown as histograms in figs .
[ fig - e - gluon ] and [ fig - e - quark ] are the energy distributions of the outgoing parton in some elected channels of parton - medium scattering .
the agreement between monte carlo results and those from direct integration in eq .
( [ scattrate2 ] ) provides numerical verifications of the monte carlo code for parton - medium scatterings in the lbt model .    for the purpose of discussing jet - induced medium excitation
later , we show in figs .
[ fig - theta - gluon ] and [ fig - theta - quark ] distributions in the polar angle of the thermal recoil parton from parton - medium scatterings via the same selected channels as in figs .
[ fig - e - gluon ] and [ fig - e - quark ] from both monte carlo simulations ( histograms ) and numerical integration of eq .
( [ scattrate2 ] ) ( solid lines ) .
the polar angle of the thermal recoil parton in two different energy ranges , @xmath910 - 1 and 1 - 2 gev , is calculated with respect to the direction of the incident parton with @xmath84 gev .
as shown in the figures , thermal recoil partons with energy significantly larger than the typical thermal energy gain their energy and momentum through elastic scattering and therefore are dragged along in the same direction of the energetic incident parton .
soft thermal recoil partons with @xmath37=0 - 1 gev , on the other hand , have large angles , almost perpendicular to the incident parton on the average .
these features of parton - medium scattering will determine the jet - induced medium excitation due to multiple parton scattering as we will describe later .      in the lbt model ,
we initiate the monte carlo program by calculating parton scattering rates according to eqs .
( [ singlerate ] ) and ( [ totalrate ] ) for a range of incident energies and medium temperatures which are stored in two - dimensional tables for later use .
we then approximate the propagation of partons in medium with a discretized time internal @xmath92 in the frame of hydrodynamic evolution of the qgp medium .
we sample the probability for @xmath93 number of parton - medium scatterings within a time interval @xmath92 according to a poisson distribution , @xmath94^n}{n ! }
e^{-\delta n_a(e,\delta t ) } , \\ \delta n_a(e,\delta t)&=&\delta t
\frac{p\cdot u}{e } \gamma_a ( p\cdot u , t),\end{aligned}\ ] ] where the scattering rate @xmath95 for a parton @xmath34 is evaluated at the local fluid comoving frame with incident energy @xmath96 and @xmath17 is the local fluid four - velocity . after determining the number of parton scatterings @xmath93 , one can assume these @xmath93 number of scatterings occur sequentially within the time interval @xmath92 and energy and momentum are strictly conserved in each scattering along the classical trajectory .    if we choose the value of time interval @xmath92 to be much smaller than the mean - free - path length @xmath97 , one can approximate the probability for at least one parton - medium scattering , @xmath98 as the probability for one parton - medium scattering during the time interval @xmath92 .    for each parton - medium scattering ,
we use fractional rates @xmath99 to decide the channel and flavors of the initial thermal parton @xmath11 and outgoing partons @xmath4 and @xmath13 . the kinematics of the parton - parton scattering @xmath100 is sampled according to the differential rate in eq .
( [ scattrate2 ] ) .
the more energetic parton of the two outgoing partons is chosen as the leading parton for the next parton - medium scattering . in the lbt model ,
all outgoing partons , both the leading and thermal recoil partons , in each parton - medium scattering are recorded and are allowed to go through further parton - medium scattering in the subsequent time intervals .
each parton - medium scattering could in principle accompanied by induced gluon bremsstrahlung .
we focus mainly on elastic scattering processes in the lbt model in this paper and leave the discussion on the implementation of induced gluon bremsstrahlung for the subsequent publication .    to take into account of the back - reaction in the boltzmann transport equation ( [ lbt ] )
, we also record the initial thermal parton @xmath11 and its four - momentum @xmath101 in each scattering process , which we denote as `` negative '' partons and are allowed to transport further according to the boltzmann equation .
these `` negative '' partons should be subtracted from the final parton spectra and energy - momentum density of the jet - induced medium excitation .
thermal recoil partons and thenegative " partons are collectively called _ jet - induced medium partons _ in this study .
in the lbt model , all channels of @xmath8 scattering are considered for parton - medium interaction in which the flavor of the most energetic parton in the final state can be different from the incident parton . in the lbt monte carlo simulations , we designate the more energetic one of the two final partons as the leading parton .
one can then follow the propagation of the leading parton and study its transverse momentum broadening and elastic energy loss .      during the propagation of an energetic parton in a qgp medium
, each parton - medium scattering contributes to the transverse momentum ( perpendicular to the initial direction of the propagating parton ) of the final leading parton according to matrix elements of the @xmath8 parton - parton scattering , leading to the increase of the averaged transverse momentum squared over time or the transverse momentum broadening .
the averaged transverse momentum squared per unit length is often defined as the jet transport parameter @xmath103 @xcite-@xcite which characterizes both the local gluon number density and the strength of jet - medium interaction .     or transverse momentum broadening squared per unit length of the leading parton for ( a ) an initial gluon or ( b ) quark with different initial energies going through a single ( open symbols ) or multiple scatterings ( solid symbols ) in a uniform and static qgp medium at different temperatures .
solid lines are analytic results from a single scattering with the small - angle approximation . , title="fig:",width=321 ] +   or transverse momentum broadening squared per unit length of the leading parton for ( a ) an initial gluon or ( b ) quark with different initial energies going through a single ( open symbols ) or multiple scatterings ( solid symbols ) in a uniform and static qgp medium at different temperatures .
solid lines are analytic results from a single scattering with the small - angle approximation . ,
title="fig:",width=321 ] +    employing the small - angle approximation of @xmath8 cross sections in eq .
( [ eq - small - el ] ) , one can estimate the transverse momentum broadening squared per mean - free - path length as , @xmath104 where @xmath105 is the transverse momentum transfer per scattering which is independent of the parton s flavor , @xmath106 for a quark and @xmath107 for a gluon .
the variable in the logarithm is defined as @xmath108 averaged over the kinematics of a single parton - medium scattering and @xmath15 is the center - of - mass energy squared for each parton - parton scattering .
one can assume , however , @xmath109 with the constant @xmath4 determined from the numerical calculations .    during multiple scatterings in the lbt model
, one can follow the propagation of the leading parton and study its final transverse momentum distribution and accumulated transverse momentum broadening . shown in fig .
[ qhat ] are the transverse momentum broadening per unit length from single and multiple scattering of an energetic parton with different initial energy @xmath50 in a static and uniform qgp at different temperatures over a length of @xmath110 fm .
gaussian fits are shown in dashed lines .
, title="fig:",width=321 ] +   mev .
gaussian fits are shown in dashed lines .
, title="fig:",width=321 ]    we also show in fig .
[ pt - broadening ] the transverse momentum distributions ( solid lines ) of the leading parton at different times for an initial gluon ( upper panel ) or quark ( lower panel ) with energy @xmath111 gev propagating in a uniform medium at a temperature @xmath112 mev .
the transverse momentum broadening manifests in these distributions as the increase of the widths of distributions with time .
such an increase of the width gives rise to an approximately time - independent jet transport parameter @xmath103 as shown in fig .
[ qhat ] . to illustrate the evolution of the transverse momentum distribution with time beyond the change of the width
, we also compare each distribution with a gaussian , @xmath113 where @xmath114 is the average transverse momentum squared from each corresponding transverse momentum distribution from lbt simulations .
the @xmath102 distribution from a single scattering follows a power law @xmath115 at small angles @xmath116 . at later times
, the distribution deviates from the power law behavior due to multiple scattering . according to eqs .
( [ analytic_gammag ] ) and ( [ analytic_gammaq ] ) , the total scattering rate for a gluon ( quark ) is @xmath117 ( 0.68 ) 1/fm .
the leading parton should experience several scatterings with medium partons after a few fm/@xmath4 of time .
its transverse momentum distribution will start to deviate from the power law and approach to a gaussian according to the central limit theorem because the transverse momentum transfers during each scattering are independent of each other .
since a gluon s scattering rate is more than twice that of a quark , its approach a gaussian form of the @xmath102 distribution is also twice sooner than a quark . in the high @xmath102 region , the power - law - like tail , however , will always remain due to a single hard or large - angle scattering .
this evolution of the transverse momentum distribution of the leading parton is also shown in a recent perturbative calculation within a htl pqcd approach @xcite .      following the propagation of the leading parton
, one can also calculate the elastic energy loss in the lbt model .
shown in fig .
[ fig : elosse ] is the calculated elastic energy loss per unit length from lbt simulations as a function of the initial parton s energy for a gluon ( upper panel ) or a quark ( lower panel ) propagating in a uniform qgp medium with a length @xmath110 fm at different constant temperatures . also shown in the figure are the simulated elastic energy loss per unit mean - free - path length from a single scattering , @xmath118 where @xmath119 is the energy transfer of the leading parton to a thermal medium parton which should depend on the energy @xmath120 of the thermal parton and the transverse momentum transfer squared @xmath121 in each scattering .
results from a single scattering are in good agreement with that from multiple scatterings over a long distance , especially for more energetic quarks at lower temperatures .
the difference between results from a single and multiple scatterings becomes bigger for lower energy gluons at higher temperature , which can be understood by the energy dependence of the energy loss in a single scattering .
using the small - angle approximation for elastic scattering cross sections in eq .
( [ eq - small - el ] ) and @xmath122 , the elastic energy loss per mean free path can be expressed as the product of the averaged transverse momentum transfer squared per scattering @xmath123 and the @xmath124-weighted scattering rate @xcite : @xmath125 one obtains the analytic form of the elastic energy loss per unit mean free path , @xmath126 under the small - angle scattering approximation .
this analytic result agrees with the full htl pqcd result @xcite in the leading logarithmic approximation and has a logarithmic dependence on the initial parton energy . as shown in fig .
[ fig : elosse ] ( solid lines ) , this analytic result under small - angle scattering approximation agrees well with lbt simulations of a single scattering .     of ( a ) a gluon or ( b ) a quark in a uniform and static medium with a temperature @xmath38 , 300 and 400 mev as a function of the initial energy @xmath50 from lbt simulations of a single scattering ( open symbols ) or multiple scatterings ( @xmath127=8fm ) ( solid symbols ) as compared to analytic results ( solid lines ) with a small - angle approximation.,title="fig:",width=321 ] +   of ( a ) a gluon or ( b ) a quark in a uniform and static medium with a temperature @xmath38 , 300 and 400 mev as a function of the initial energy @xmath50 from lbt simulations of a single scattering ( open symbols ) or multiple scatterings ( @xmath127=8fm ) ( solid symbols ) as compared to analytic results ( solid lines ) with a small - angle approximation.,title="fig:",width=321 ] +    the logarithmic energy dependence of the elastic energy loss can also explain the difference between single and multiple scatterings in fig .
[ fig : elosse ] . during multiple scatterings in the lbt model , energy and momentum
is conserved during each scattering .
the leading parton will lose its energy along its propagation path .
the reduced energy of the leading parton will then lead to smaller elastic energy loss for subsequent scatterings according to the logarithmic energy dependence .
this will lead to an overall reduction of the averaged energy loss per unit distance during multiple scattering over a long distance .
this reduction on the averaged energy loss per unit distance is more significant when the total energy loss @xmath128 is comparable to parton s initial energy .
since a gluon s energy loss is @xmath129 times larger than that of a quark , this degradation of energy and the energy loss is more significant for a gluon with lower energies in a qgp medium at higher temperatures as shown in fig .
[ fig : elosse ] . to illustrate this degradation of energy loss over time
, we show in figs .
[ fig : elosstimeg ] and [ fig : elosstimeq ] the averaged elastic energy loss per unit distance as a function of propagation time for different initial parton energies and medium temperatures .
it is clear that the degradation of parton energy loss with time is the strongest for a gluon with the lowest initial energy and in a medium at the highest temperature .
for an initial high energy quark , the energy loss remains almost a constant even at a later time .     of a gluon with different initial energies in a uniform medium with a constant temperature ( a ) @xmath38 mev and ( b ) @xmath130 mev as a function of the propagation time.,title="fig:",width=321 ] +   of a gluon with different initial energies in a uniform medium with a constant temperature ( a ) @xmath38 mev and ( b ) @xmath130 mev as a function of the propagation time.,title="fig:",width=321 ] +     for a quark.,title="fig:",width=321 ] +    for a quark.,title="fig:",width=321 ] +    in the lbt model , energy and momentum
are conserved in each parton - medium scattering ( with the energy and momentum of  negative " partons subtracted from the sum ) along the path of the parton propagation .
the energy lost by the leading parton as shown in the above calculations will be carried by thermal recoil partons , minus that of  negative " partons , in the form of jet - induced medium excitations . to illustrate the transfer of energy from the leading partons to jet - induced medium excitations , we plot in fig .
[ spectra - time ] the energy spectra of all partons including both leading and the thermal recoil partons ( with the  negative " partons subtracted ) at different times ( @xmath131 , 6 , and 10 fm/@xmath4 ) during the propagation of an initial gluon with energy @xmath111 gev in a uniform qgp medium at a constant temperature @xmath112 mev . at the beginning of the parton propagation , leading partons peak around the initial energy @xmath50 with power law tails towards small values due to the dominant @xmath16-channel scattering , while thermal recoil partons ( minus  negative " partons )
peak around their thermal energy @xmath132 with a power law tail . at later times ,
leading partons continue to lose energy due to multiple scattering and peak at smaller energies @xmath133 with increasingly broadened distributions . the number of jet - induced medium partons ( thermal recoil partons minus  negative " partons ) , on the other hand , will continue to increase with softened spectra which resemble a thermal distribution at later times . in principle , the spectra of leading and jet - induced medium partons will eventually merge and approach a thermal spectrum as a result of the equilibration . but
this limit is beyond the validity of the lbt model .
gev in a uniform qgp medium at a constant temperature @xmath112 mev.,width=321 ]
during the parton propagation in a qgp medium , the energy and momentum deposited in the medium should dissipate with time and propagate like a sound wave .
since the leading massless parton travels at the speed of light which is much faster than the maximum sound velocity @xmath134 in an ideal qgp , the sound wave induced by the propagating parton becomes supersonic .
it can form a mach - cone shock wave as shown in many recent studies @xcite .
simulations of jet propagation within parton transport models also show similar mach - cone - like features @xcite .
studies of such jet - induced medium excitations can help us to understand the experimental data on dissipation of energy lost by quenched jets @xcite and extract bulk medium properties such as sound velocity or the equation of state ( eos ) of the qgp medium .
since this article is devoted to a general description and test of the lbt model in a uniform medium , we will focus on the effect of dynamical scatterings of propagating partons in a uniform medium .
this is similar to the effect of viscosity on the propagation of shockwaves in a qgp medium @xcite .
we will discuss the propagation of jet - induced mach cones in an expanding qgp medium in the future when we carry out realistic studies of jet - medium interaction in high - energy heavy - ion collisions .      to study jet transport and jet - induced medium excitations at the same time in a qgp medium , both of the outgoing partons in each @xmath8 scattering
are tracked in lbt monte carlo simulations .
they are allowed to undergo further propagation and scattering in the medium .
of the two outgoing partons , the more energetic one is designated as the leading parton for further propagation and scattering .
the flavor of the outgoing leading parton can be different from incident parton though this is quite rare since the scattering is dominated by small - angle @xmath16-channel processes .
the less energetic or soft parton from each parton - medium scattering is considered as the recoil medium parton that is also allowed to further propagate and interact with the medium .
in addition , the initial parton from the thermal medium participating in the scattering is also tracked and allowed to further propagate in the lbt model whose four - momentum is subtracted from the final result to account for the back - reaction in the boltzmann transport equation .
these thermal partons from the back - reaction are denoted as  negative " partons in our study and are included in all final results on parton spectra and jet reconstruction . in general , we refer to the thermal recoil and  negative " partons as jet - induced thermal partons .
since jet - induced thermal partons are allowed to propagate in the qgp medium according to the boltzmann transport equation in the lbt model , we can study the linear response of the medium to jet - medium interaction along the path of a jet s propagation .
interactions among jet - induced medium partons are neglected .
such an approximation of linear response is valid when the jet - induced medium excitation is much smaller than the thermal background , @xmath135 .
gev along ( @xmath51 ) and transverse ( @xmath136 ) to the initial parton direction after ( a ) 4 and ( b ) 8 fm of propagation in a uniform qgp medium at a constant temperature @xmath112 mev .
the jet parton starts propagating along the @xmath51 direction at @xmath137 and @xmath138.,title="fig:",width=340 ] +   gev along ( @xmath51 ) and transverse ( @xmath136 ) to the initial parton direction after ( a ) 4 and ( b ) 8 fm of propagation in a uniform qgp medium at a constant temperature @xmath112 mev .
the jet parton starts propagating along the @xmath51 direction at @xmath137 and @xmath138.,title="fig:",width=340 ]    the energy and momentum lost to the medium by the propagating parton will act like a source for jet - induced medium excitation through the propagation and transport of the jet - induced medium partons . shown in fig .
[ fig : g - cone ] are two - dimensional and contour energy density distributions of the medium excitation induced by a propagating gluon with an initial energy @xmath111 gev along ( @xmath51 ) and transverse ( @xmath136 ) to the initial parton direction after 4 fm/@xmath4 ( upper panel ) and 8 fm/@xmath4 ( lower panel ) of propagation in a uniform qgp medium at a constant temperature @xmath112 mev .
the jet parton starts propagating along the @xmath51 direction at @xmath137 and @xmath138 .
note that in each event of parton propagation there are only a limited number of jet - induced medium partons .
figure [ fig : g - cone ] shows the energy density distributions of jet - induced medium partons or medium excitation averaged over many events .
the number of medium recoil and  negative " partons increases separately with time in lbt monte carlo simulations .
the net energy density , however , remains approximately constant because of energy - momentum conservation and the near constant energy - momentum transfer rate via parton - medium interaction .
one can see clearly the formation and propagation of a mach - cone - like shock wave trailing the leading parton .
the edge of the shock wave travels at a speed limited by the velocity of light while the peak of the shock wave propagates at the effective sound velocity in the medium . the shock wave is significantly diffused during its propagation because of the dissipation due to finite values of viscosities as a result of parton - parton collisions as implemented in the lbt model .
the energy density of jet - induced medium partons is negative behind the propagating parton .
this is the diffusion wake induced by the jet - medium interaction which essentially depletes the thermal medium parton density behind the path of the propagating parton .
the diffused shock wave in the shape of an elliptic paraboloid induced by jet - medium interaction within the lbt model as shown in fig .  [ fig : g - cone ] is sharply different from the mach - cone shape given by other linear response theories @xcite , hydrodynamic simulations @xcite or anti - de sitter / conformal field theory ( ads / cft ) correspondence studies @xcite with a source term that travels along a given direction .
the leading parton in the lbt monte carlo simulations , however , does not travel along a fixed direction .
it receives transverse momentum transfers in random directions and accumulates transverse momentum broadening along its propagation path .
such momentum broadening amounts to a smearing in the direction of the source term for jet - induced medium excitation and the width of the smearing increases with the propagation time .
this eventually leads to the elliptic paraboloid shape of the shock wave instead of a sharp mach - cone .
( solid ) , 1.5 - 1.75 ( dashed ) , 3.0 - 3.25 gev ( dotted ) induced by a gluon with an initial energy @xmath111 gev after ( a ) @xmath131 , ( b ) 6 and ( c ) 10 fm/@xmath4 of propagation in a uniform qgp medium at a constant temperature @xmath112 mev.,title="fig:",width=302 ] +   ( solid ) , 1.5 - 1.75 ( dashed ) , 3.0 - 3.25 gev ( dotted ) induced by a gluon with an initial energy @xmath111 gev after ( a ) @xmath131 , ( b ) 6 and ( c ) 10 fm/@xmath4 of propagation in a uniform qgp medium at a constant temperature @xmath112 mev.,title="fig:",width=302 ] +   ( solid ) , 1.5 - 1.75 ( dashed ) , 3.0 - 3.25 gev ( dotted ) induced by a gluon with an initial energy @xmath111 gev after ( a ) @xmath131 , ( b ) 6 and ( c ) 10 fm/@xmath4 of propagation in a uniform qgp medium at a constant temperature @xmath112 mev.,title="fig:",width=302 ] +      to illustrate consequences of the transverse momentum broadening of leading partons on the final spectra of jet - induced medium partons , we show in fig .
[ fig : g - angle100 ] angular distributions of jet - induced medium partons with different ranges of energy relative to the initial direction of a gluon with energy @xmath111 gev after @xmath131 ( upper panel ) , 6 ( middle panel ) , and 10 fm/@xmath4 ( lower panel ) of propagation in a uniform medium at a constant temperature @xmath112 mev . in a uniform medium , event - averaged parton distributions should be azimuthal symmetric with respect to the initial direction of the propagating parton .
for the convenience of future study of azimuthal distributions in a cylindrical frame along the beam direction in high - energy heavy - ion collisions , we project the angular distribution onto a plane aligned with the initial propagating parton , @xmath139\right),\ ] ] where @xmath140 and @xmath141 are the polar and azimuthal angles , respectively , relative to the initial jet parton direction .
we refer to the angle @xmath142 as the projected azimuthal angle and the distribution as the projected azimuthal distribution .
soft jet - induced medium partons ( @xmath143 gev ) are mostly  negative " partons from the the back - reaction whose contributions to the spectra are negative according to our subtraction scheme .
these  negative " partons have a broad angular distribution that does not change much over time .
contributions to more energetic jet - induced medium partons from  negative " partons decrease with time , especially along the initial direction of the propagating parton , and thermal recoil partons dominate .
intermediate energy thermal recoil partons during the early time have an angular distribution similar to that of a single scattering as shown in figs .
[ fig - theta - gluon ] and [ fig - theta - quark ] that has two peaks away from the initial direction . at later times
, however , the double - peak structure disappears due to transverse momentum broadening of the leading parton and diffusion of thermal recoil partons due to multiple scatterings .
the effect of the transverse momentum broadening on the angular distribution of thermal recoil partons is stronger for less energetic initial partons at later times , especially when the total energy loss becomes sizable relative to the initial parton energy .
in high - energy heavy - ion collision experiments , reconstructed jets have become a powerful tool for the study of jet quenching @xcite .
these jets are collimated clusters of hadrons and are reconstructed from the calorimetric energy of final hadrons within a jet cone , @xmath144 using a jet finding algorithm @xcite , where @xmath145 ( @xmath146 ) and @xmath142 ( @xmath147 ) are the pseudorapidity and azimuthal angle of the final hadron ( jet ) , respectively . in heavy - ion collisions , one should also subtract the background energy within the jet cone from hadrons in the underlying events .
such background subtraction is extremely nontrivial and may affect the final jet energy and jet spectra .
if we neglect nonperturbative effects of parton hadronization , we can approximate final hadronic jets with partonic jets reconstructed from final partons in theoretical simulations . in our lbt monte carlo simulations , we assume a perfect subtraction of the background from underlying events if the initial propagating parton is not present .
we achieve this by including both leading partons as well as jet - induced medium partons for the jet reconstruction .
we use a modified version of the fastjet code @xcite with the anti-@xmath7 algorithm for jet reconstruction in our study within the lbt model in which a  negative " parton is treated as a normal one when its pseudorapidity and azimuthal angle are calculated , but its energy and momentum are subtracted from the total energy and momentum of the reconstructed jet in each iteration of the jet - finding algorithm .
we only consider the leading jet per event through the fastjet jet finding package with the given algorithm .      during the propagation of a fast parton , it undergoes multiple scatterings with medium partons in the lbt model .
the leading parton loses energy and experiences transverse momentum broadening .
the lost energy and momentum transfer are then carried by jet - induced medium partons which will propagate and go through further scattering with the medium leading to jet - induced medium excitations .
some fraction of the energy and momentum carried by jet - induced medium partons will still be confined within the jet cone .
however , an increasing fraction will be transported outside the jet cone via large - angle scattering and diffusion of jet - induced medium partons over the course of the parton propagation .
this will lead to a reduction of the energy of the reconstructed jet . shown in fig .
[ jet - eloss ] is the energy loss of reconstructed jets with a cone size @xmath148 as a function of time for an initial gluon with energy @xmath149 and 100 gev in a uniform qgp medium with a constant temperature @xmath112 mev .
the solid symbols represent the energy loss of jets reconstructed with all of jet - induced medium partons while the open symbols represent jets without  negative " partons .
the jet energy loss increases linearly with time or distance of propagation similarly to the elastic energy loss of a propagating parton .
mev as a function of propagation time of an initial gluon with energy @xmath149 and 100 gev .
jets are reconstructed with all partons ( leading and jet - induced medium partons ) ( solid symbols ) or without  negative " partons ( open symbols ) .
, title="fig:",width=321 ] +    th effect of  negative " partons from the diffusion wake is negligible during the early stage of the parton propagation because there are very few of them . at the later stage , as the number of these  negative " partons becomes large , they deplete significantly the thermal medium behind the propagating parton and effectively modify the background underlying the jet .
when such modification of the underlying background is taken into account via the subtraction of  negative " partons from the jet cone , the effective jet energy loss becomes bigger as shown in fig .
[ jet - eloss ] . without subtraction of  negative " partons
, the jet energy loss tends to saturate at later times .
the linear time dependence of the total effective jet energy loss is restored only after the energy of  negative " partons is subtracted .    during the propagation of an energetic parton , the leading parton experiences transverse momentum broadening through interaction with medium partons .
the thermal recoil parton in the same scattering should also carry the same amount of transverse momentum as the leading parton but in the opposite direction . if a reconstructed jet contains both the leading parton and the thermal recoil partons ( minus the  negative " partons ) , its total transverse momentum should not be influenced by the scattering .
however , reconstructed jets do not contain all jet - induced medium partons because of the finite jet - cone size @xmath150 .
large - angle scatterings and diffusion of jet - induced medium partons should lead to a reduction of the jet energy , as shown in the above , as well as a net transverse momentum with respect to the initial parton direction . shown in fig .
[ fig - jet - pt ] are the transverse momentum distributions of reconstructed jets with respect to the initial direction of a propagating gluon with energy @xmath111 gev in a uniform medium at a constant temperature @xmath112 mev .
there is an apparent transverse momentum broadening at later times during jet transport , similar to the broadening of leading partons .
there is also a power law tail in the transverse momentum distribution .
this is caused mainly by a single large angle parton - medium scattering during the jet propagation .
gev and ( b ) @xmath111 gev in a uniform qgp medium at a temperature @xmath112 mev .
jets are reconstructed with all partons ( leading and jet - induced medium partons ) ( solid symbols ) or without  negative " partons ( open symbols).,title="fig:",width=321 ] +   gev and ( b ) @xmath111 gev in a uniform qgp medium at a temperature @xmath112 mev .
jets are reconstructed with all partons ( leading and jet - induced medium partons ) ( solid symbols ) or without  negative " partons ( open symbols).,title="fig:",width=321 ] +      reconstructed jets in lbt monte carlo simulations in this study contain both leading partons and jet - induced medium partons inside the jet cone .
to examine the composition of reconstructed jets and their evolution with time , we show parton distributions within reconstructed jets as functions of the longitudinal momentum fraction @xmath151 in fig .
[ jet - ff ] and @xmath152 in fig .
[ jet - ff0 ] , where @xmath153 is a parton s longitudinal momentum along the direction of a reconstructed jet with energy @xmath154 .
since we start with a single energetic parton , the initial parton distribution function should be a @xmath53 function at @xmath155 .
jet - induced medium partons from parton - medium interaction and leading partons with reduced energy populate the distribution at @xmath156 or @xmath157 after the jet transport begins .
further scatterings of both leading partons and jet - induced medium partons lead to the enhancement of soft partons at @xmath158 or @xmath159 as shown in figs .
[ jet - ff ] and [ jet - ff0 ] .
since the definition of momentum fraction @xmath160 is normalized by the reconstructed jet energy @xmath154 which is still dominated by the energetic leading parton , the parton distribution in @xmath160 continues to have a peak at @xmath161 .
the parton distribution in @xmath162 which is normalized by the initial parton energy @xmath50 is , however , suppressed at large @xmath163 and its initial peak at large @xmath163 disappears at later times as the leading parton continues to lose energy during its propagation through the medium . similar behavior was reported in the study of @xmath164-triggered jets in high - energy heavy - ion collisions @xcite .
the modified fragmentation functions in terms of momentum fraction of the original jet energy ( or the @xmath164 s energy ) are much more sensitive to parton energy loss than that in terms of the momentum fraction of the reconstructed jet energy .
since  negative " partons are subtracted from both the parton distribution and the reconstructed jet energy ( reducing @xmath154 ) , they will reduce the soft parton distributions at small @xmath165 or @xmath159 and enhance the parton distribution at large @xmath161 .     within the jet cone @xmath148 at different times for a propagating gluon with initial energy ( a ) @xmath149 and ( b ) 100 gev in a uniform qgp medium at a temperature @xmath112 mev .
solid symbols represent all partons including leading , recoil , and  negative " partons while the open symbols show results without  negative " partons.,title="fig:",width=321 ] +   within the jet cone @xmath148 at different times for a propagating gluon with initial energy ( a ) @xmath149 and ( b ) 100 gev in a uniform qgp medium at a temperature @xmath112 mev .
solid symbols represent all partons including leading , recoil , and  negative " partons while the open symbols show results without  negative " partons.,title="fig:",width=321 ] +     within the jet cone @xmath148 at different times for a propagating gluon with initial energy ( a ) @xmath149 and ( b ) 100 gev in a uniform qgp medium at a temperature @xmath112 mev .
solid symbols represent all partons including leading , recoil , and  negative " partons while the open symbols show results without  negative " partons.,title="fig:",width=321 ] +   within the jet cone @xmath148 at different times for a propagating gluon with initial energy ( a ) @xmath149 and ( b ) 100 gev in a uniform qgp medium at a temperature @xmath112 mev .
solid symbols represent all partons including leading , recoil , and  negative " partons while the open symbols show results without  negative " partons.,title="fig:",width=321 ] +     and ( b ) 100 gev at different times in a static and homogeneous medium at temperature @xmath112 mev .
results that include all ( leading+positive+negative ) partons are presented in solid symbols while results without  negative " partons are presented in open symbols .
, title="fig:",width=321 ] +   and ( b ) 100 gev at different times in a static and homogeneous medium at temperature @xmath112 mev .
results that include all ( leading+positive+negative ) partons are presented in solid symbols while results without  negative " partons are presented in open symbols .
, title="fig:",width=321 ] +    to study the evolution of the jet transverse structure , we show in fig .
[ jet - profile ] the transverse profile of reconstructed jets at different times for a propagating gluon with an initial energy @xmath149 gev ( upper panel ) and @xmath111 gev ( lower panel ) in a uniform qgp medium at a constant temperature @xmath112 mev .
the jet transverse profile is defined as the average fraction of jet energy inside an annulus in the @xmath166 plane , @xmath167 as a function of @xmath168 , where @xmath169 is the summed energy of all partons in the annulus between radius @xmath170 and @xmath171 inside the jet cone .
when a high energy parton propagates through the plasma , it loses energy to the medium .
the lost energy is carried away by thermal recoil partons whose further diffusion through multiple scatterings will transport some of the lost energy outside the jet cone leading to the jet energy loss and broadening of the transverse profile toward the edge of the jet cone , as seen in fig .
[ jet - profile ] .
the  negative " partons are mostly soft partons at large angles away from the initial parton s direction .
when their energy is subtracted from the reconstructed jet , the broadening of the transverse profile will be reduced , especially near the edge of the jet cone .
in this paper we have described in detail the basic elements of the lbt model for the study of jet transport in a qgp medium with the complete set of elastic @xmath8 scattering processes .
we verified the monte carlo code in the lbt model by direct comparisons to semi - analytic results of total and differential scattering rates via numerical integrations .
we calculated the transverse momentum broadening and elastic energy loss of the leading parton along the path of the parton propagation . within the pqcd approach to each elastic scattering between the leading and thermal medium parton , both the jet transport coefficient @xmath103 and the energy loss per mean - free - path
have a logarithmic dependence on the parton energy .
such an energy dependence leads to some nontrivial time or distance dependence of the average transverse momentum broadening and energy loss per unit distance due to the decrease of the leading partons energy over time .
this dependence is more significant when the total energy loss becomes comparable to the initial parton s energy .
we have also illustrated the dissipation of the energy and momentum lost by a propagating parton in the medium via transport of jet - induced medium partons ( thermal recoil partons and  negative " partons from the back - reaction ) in the lbt model .
the transport of these jet - induced medium partons within the lbt model effectively forms a supersonic shock wave and a diffusion wake behind the leading parton .
the supersonic wave has a diffused mach - cone shape because of the transverse momentum broadening of the leading parton along its path .
the energy spectra of jet - induced medium partons in the supersonic wave are shown to resemble that of a thermal distribution at later times of the jet propagation .    using a modified version of the fastjet jet - finding package with the anti-@xmath7 algorithm @xcite in which the energy of  negative " partons from the back - reaction is subtracted , we show the effect of jet - induced medium partons and their dissipation in medium on reconstructed jets .
transport of these jet - induced medium partons outside the jet cone constitutes a significant reduction of jet energy loss which should be taken into account for any realistic study of jet suppression in high - energy heavy - ion collisions .
their inclusion in the jet reconstruction also influences the transverse momentum broadening , parton distributions ( fragmentation functions ) and transverse profiles of reconstructed jets .
transverse momentum distributions of both leading partons and reconstructed jets are shown to have a power law tail underneath gaussian distributions from multiple scattering .
these power law tails come from a single large angle parton - medium scattering during the jet propagation .
study of these power law tails in heavy - ion collisions can shed light on the microscopic nature of the qgp medium at different scales .    for a realistic description of jet transport in the qgp medium in high - energy heavy - ion collisions
, we will have to implement inelastic processes such as gluon bremsstrahlung induced by multiple parton - medium scattering @xcite and couple lbt monte carlo simulations with bulk medium evolutions from 3 + 1d hydrodynamic model calculations , which should be constrained by bulk hadron spectra from existing experimental data .
these parts of the lbt model and the phenomenological study of jet suppression in heavy - ion collisions will be discussed in subsequent publications .
we thank wei chen for helpful discussions , in particular on the elliptic paraboloid shape of the shock wave induced by a propagating parton .
this work is supported by china most under grant no .
2014dfg02050 , the nsfc under grant no . 11221504 , the major state basic research development program in china ( grant no .
2014cb845404 ) , by the director , office of energy research , office of high energy and nuclear physics , division of nuclear physics of the u.s .
department of energy under contract no .
de - ac02 - 05ch11231 and within the framework of the jet collaboration .
y .- z . is also supported by european research council grant no . hotlhc erc-2011-stg-279579 .
m.  gyulassy and m.  plumer , phys .
b * 243 * , 432 ( 1990 ) .
x.  n.  wang and m.  gyulassy , phys .
* 68 * , 1480 ( 1992 ) .
j.  adams _ et al . _ [ star collaboration ] , nucl .  phys .  a * 757 * , 102 ( 2005 ) . k.  adcox _ et al . _ [ phenix collaboration ] , nucl .  phys .  a * 757 * , 184 ( 2005 ) .      c.  adler _ et al . _ , phys .  rev
.  lett .   * 89 * , 202301 ( 2002 ) .
k.  aamodt _ et al . _
[ alice collaboration ] , phys .
b * 696 * , 30 ( 2011 ) .
s.  chatrchyan _ et al . _
[ cms collaboration ] , eur .
j.  c * 72 * , 1945 ( 2012 ) . c.  adler _ et al .
_ , phys .
* 90 * , 082302 ( 2003 ) .
a.  adare _ et al .
_ [ phenix collaboration ] , phys .
c * 80 * , 024908 ( 2009 ) .
b.  i.  abelev _ et al . _ [ star collaboration ] , phys .
c * 82 * , 034909 ( 2010 ) .
x.  -n .
wang , nucl .  phys .   * a750 * , 98 - 120 ( 2005 ) . i.  vitev and m.  gyulassy , phys
.  rev .
* 89 * , 252301 ( 2002 ) .
x.  n.  wang , phys .
b * 595 * , 165 ( 2004 ) .
t.  renk , phys .
c * 74 * , 024903 ( 2006 ) .
h.  zhang , j.  f.  owens , e.  wang and x.  n.  wang , phys .
rev .  lett .
* 98 * , 212301 ( 2007 ) .
g.  y.  qin , j.  ruppert , c.  gale , s.  jeon , g.  d.  moore and m.  g.  mustafa , phys .
* 100 * , 072301 ( 2008 ) .
t.  renk , phys .
c * 74 * , 034906 ( 2006 ) .
h.  zhang , j.  f.  owens , e.  wang and x .-
wang , phys .
.  lett .   *
103 * , 032302 ( 2009 ) .
g.  y.  qin , j.  ruppert , c.  gale , s.  jeon and g.  d.  moore , phys .
c * 80 * , 054909 ( 2009 ) .          c.  young , b.  schenke , s.  jeon and c.  gale , phys .
c * 84 * , 024907 ( 2011 ) .
y.  he , i.  vitev and b.  -w .
zhang , phys .
b * 713 * , 224 ( 2012 ) .
t.  renk , phys .
c * 85 * , 064908 ( 2012 ) .
w.  dai , i.  vitev and b.  -w .
zhang , phys .
* 110 * , 142001 ( 2013 ) .
m.  gyulassy and x.  n.  wang , nucl .
b * 420 * , 583 ( 1994 ) .
r.  baier , y.  l.  dokshitzer , s.  peigne and d.  schiff , phys .
b * 345 * , 277 ( 1995 ) .
r.  baier , y.  l.  dokshitzer , a.  h.  mueller , s.  peigne and d.  schiff , nucl .
b * 484 * , 265 ( 1997 ) .
b.  g.  zakharov , jetp lett .
* 63 * , 952 ( 1996 ) .
m.  gyulassy , p.  levai and i.  vitev , phys .
* 85 * , 5535 ( 2000 ) .
m.  gyulassy , p.  levai and i.  vitev , nucl .
b * 594 * , 371 ( 2001 ) .
u.  a.  wiedemann , nucl .
b * 588 * , 303 ( 2000 ) .
u.  a.  wiedemann , nucl .
a * 690 * , 731 ( 2001 ) .
x.  f.  guo and x.  n.  wang , phys .
lett .   *
85 * , 3591 ( 2000 ) .
x.  n.  wang and x.  f.  guo , nucl .
a * 696 * , 788 ( 2001 ) .
p.  arnold , g.  d.  moore and l.  g.  yaffe , jhep * 0111 * , 057 ( 2001 ) .
p.  arnold , g.  d.  moore and l.  g.  yaffe , jhep * 0206 * , 030 ( 2002 ) .
g.  -y .
qin , a.  majumder , h.  song and u.  heinz , phys .
lett .   *
103 * , 152303 ( 2009 ) .
r.  b.  neufeld and b.  muller , phys .
lett .   * 103 * , 042301 ( 2009 ) .
p.  m.  chesler and l.  g.  yaffe , phys .
d * 78 * , 045013 ( 2008 ) .
s.  s.  gubser , s.  s.  pufu and a.  yarom , phys .
lett .   * 100 * , 012301 ( 2008 ) .
i.  bouras , a.  el , o.  fochler , h.  niemi , z.  xu and c.  greiner , phys .  lett .
b * 710 * , 641 ( 2012 ) .
h.  li , f.  liu , g.  -l .
ma , x.  -n .
wang and y.  zhu , phys .
* 106 * , 012301 ( 2011 ) .
x.  -n .
wang , h.  -l .
li , f.  -m .
liu , g.  -l .
ma and y.  zhu , nucl .
a * 855 * , 469 ( 2011 ) .
x.  -n .
wang and y.  zhu , phys .
lett .   * 111 * , 062301 ( 2013 ) .
x.  -n .
wang and y.  zhu , nucl . phys .
a * 932 * , 447 ( 2014 ) .
k.  zapp , g.  ingelman , j.  rathsman , j.  stachel and u.  a.  wiedemann , eur .
j.  c * 60 * , 617 ( 2009 ) .
j.  auvinen , k.  j.  eskola and t.  renk , phys .
c * 82 * , 024906 ( 2010 ) .
x.  n.  wang , phys .
rept .   *
280 * , 287 ( 1997 ) .
a.  bazavov _ et al .
_ [ hotqcd collaboration ] , phys .  rev .
d * 90 * , no . 9 , 094503 ( 2014 ) [ arxiv:1407.6387 [ hep - lat ] ] .
x.  f.  chen , c.  greiner , e.  wang , x.  n.  wang and z.  xu , phys .
c * 81 * , 064908 ( 2010 ) .
f.deramo , m.lekaveckas , h.liu , and k.rajagopal , jhep 1305 ( 2013 ) 031 .
i.  bouras , e.  molnar , h.  niemi , z.  xu , a.  el , o.  fochler , c.  greiner and d.  h.  rischke , phys .
lett .   * 103 * , 032301 ( 2009 ) .
i.  bouras , e.  molnar , h.  niemi , z.  xu , a.  el , o.  fochler , c.  greiner and d.  h.  rischke , phys .  rev .
c * 82 * , 024910 ( 2010 ) .