Transverse momentum broadening of medium-induced cascades in expanding media

In this work, we explore the features of gluonic cascades in static and Bjorken expanding media by numerically solving the full BDIM evolution equations in longitudinal momentum fraction x and transverse momentum k\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${{\varvec{k}}}$$\end{document} using the Monte Carlo event generator MINCAS. Confirming the scaling of the energy spectra at low-x, discovered in earlier works, we use this insight to compare the amount of broadening in static and expanding media. We compare angular distributions for the in-cone radiation for different medium profiles with the effective scaling laws and conclude that the out-of-cone energy loss proceeds via the radiative break-up of hard fragments, followed by an angular broadening of soft fragments. While the dilution of the medium due to expansion significantly affects the broadening of the leading fragments, we provide evidence that in the low-x regime, which is responsible for most of the gluon multiplicity in the cascade, the angular distributions are very similar when comparing different medium profiles at an equivalent, effective in-medium path length. This is mainly due to the fact that in this regime, the broadening is dominated by multiple splittings. Finally, we discuss the impact of our results on the phenomenological description of the out-of-cone radiation and jet quenching.

In recent years a lot of studies were focused on the out-of-cone radiation or equivalently transverse-momentum dependence of produced mini-jets [39].On a theoretical level, this problem leads to the generalization of the BDMPS-Z framework [11][12][13][14][15][16][23][24][25][26], and to the formulation of the BDIM equations for gluons [40] and combined system of quarks and gluons [35].One of the actively investigated problems is to account for the expansion of the medium together with turbulent dynamics of the jet losing energy.In particular, in Ref. [43], the medium-modified gluon splitting rates for different profiles of the expanding partonic medium, namely the profiles for the static, exponential, and Bjorken expanding medium were calculated.This generalized the rate equation for the energy distribution that accounts for multiple soft scatterings with a time-dependent transport coefficient characterizing the expanding medium.This allowed quantifying the sensitivity of the inclusive jet suppression on the way how the medium expands.In the subsequent study [41], the framework was generalized to account for quarks.Apart from that, the initialstate nuclear effects, vacuum-like emissions as well as coherence effects were taken into account.This allowed to achieve a reasonable description of the nuclear modification of jets.
In this paper, we would like to generalize the framework of Ref. [41,43] by accounting for the transverse-momentum dependence of mini-jets produced during the jet quenching.This will amount to generalizing the BDIM equations to the expanding medium-dependent case.The result will allow us to determine the effect of expansion on transverse-momentum spectra of fragmenting jets and to investigate the interplay of dynamics of expansion and re-scatterings.One of the objectives of this work is to study the impact of different effects of the medium expansion on the radiation and broadening revealed by the scaling analysis in the effective evolution parameters.Secondly, we study the radial distributions of gluons in bins of x (or energy) to show where the hard and soft emissions dominate.This reveals the medium kinematical scales where the splitting and the broadening dominate, respectively, as well as differences between various medium profiles.
In a very recent paper [42], the authors studied similar physical quantities as in this work through a kinetic theory approach with thermalization effects included for an infinite static medium.In this paper, we present results for a more realistic description using a timedependent splitting rate that captures the interplay of the formation time of the radiation and the medium length for the static media.Next, we highlight qualitative differences between the media profiles by including the Bjorken expanding media and explore the sensitivity of the physical quantities on the time for the onset of the jet quenching.
The paper is organized as follows.In Sec. 2, we revisit the single-gluon distributions and collinear splitting kernels for the expanding as well as static medium.We also introduce the in-medium gluon evolution equation.Next, in Sec. 3, we discuss the scaling features of a purely radiative cascade for different medium profiles, present formal solutions, and compare them to numerical evaluations in a dedicated Monte Carlo event generator MIN-CAS.The established scaling allows corroborating the concept of an equivalent effective in-medium path length L eff , that leads to the same low-x spectrum in different medium scenarios.In Sec. 4, we study fully differential spectra of mini jets both in the static medium as well as in the medium undergoing the Bjorken expansion, obtained via numerical evaluations in MINCAS.In particular, we analyze the dependence of the fragmentation function on a polar angle.From the full distributions in the longitudinal momentum fraction x and the transverse momentum k T , alternatively x and the polar angle θ, we discuss three key quantities: a) the average k T distribution; b) the x-distribution within a certain cone limited by θ < θ max ; and c) the angular distribution in specific bins of x.Finally, in Sec. 5, we perform an in-depth study of the fraction of gluons in a given x-range that remain within a cone.In Sec. 6, we summarize our work and explore the possible qualitative impact on jet quenching phenomenology.

Medium-induced cascades in expanding media
In this work, we study two classes of medium profiles: a static, non-expanding medium and a medium following the Bjorken expansion.

Static medium:
In this case, we consider the temperature of the medium to be independent of time, T = T 0 .
Bjorken expanding medium: Assuming one-dimensional longitudinal expansion, we use the following form for the temperature evolution as a function of the proper time t, where the reference time t 0 is a free parameter.
In addition, we allow for two different initial time values t 0 in the Bjorken scenario, resulting in a total of three studied cases.Further discussion about the expanding media can be found in Refs.[28,41,[43][44][45].
We consider the gluon distribution where x is the longitudinal momentum fraction and k = (k x , k y ) is the transverse momentum.The evolution equation for this distribution in a dense medium, neglecting quark contributions, reads [46] (2.3) The above evolution equation describes the interplay between the collinear splittings (first two terms on the r.h.s. of the equation) and diffusion in momentum space (the last term).In this equation, we consider the case when the transverse-momentum transfer and broadening during the branching is neglected.The collinear branching kernel K(z, t) will be discussed in detail below.The transverse-momentum dependence comes, nevertheless, from both the elastic scattering, as described by the elastic collision kernel C(l, t), and multiple splittings, see the terms in the first line of Eq. (2.3).Finally, the characteristic time with ᾱ = α s N c /π and p + 0 being the (light-cone) energy of the initial jet particle, determines the stopping time of the jet at the initial, or reference, value of the jet quenching parameter q0 , which will be specified in more detail below.In this work, we have not included the thermalization effects for x ≤ 0.01 as in Ref. [42].
The elastic collision kernel C(l, t) is given by where w(l, t) ∝ n(t)d 2 σ el /d 2 l is proportional to the density of the medium n(t) and the in-medium elastic cross section d 2 σ el /d 2 l.When the medium expands, the density drops reducing the impact of elastic momentum transfer.Although it is less important, screening effects, contained in the elastic cross section, will also typically be hampered by the expansion 1 .
Considering purely gluon systems, we use the so-called Hard Thermal Loop (HTL) approximation for w(l, t) at the leading order: where the medium is characterized by the temperature T and the Debye mass m D , and we have not written explicitly their time dependence.We have also introduced the jet transport coefficient, which is given by q = α s N c m 2 D T with g 2 = 4πα s .The Debye mass in the QCD medium at the leading order is given by (2.7) In addition to the above perturbative collision kernels, one can also study the effect of using scattering potentials from non-perturbative extraction from the lattice QCD [48] as well as the NLO HTL contributions [36].In a static medium, m 2 D = m 2 D0 ≡ 3/2g 2 T 2 0 and q = q0 ≡ α s N c m 2 D0 T 0 are simply constant, whereas in the Bjorken model these parameters, together with the temperature, are time-dependent and given by q(t) = q0 where the time evolution of the temperature is given in Eq. (2.1).
Next, we turn to the splitting kernel.It is related to the in-medium splitting rate as K(z, t) = t * dI/(dz dt).This rate has recently been evaluated numerically in a static medium in Ref. [38].It has also been investigated in analytical resummation schemes that cover the whole kinematical phase space [20-22, 38, 49].We shall use it together with the so-called "harmonic oscillator" approximation which is valid for soft emissions in a large medium [49].The main features of soft, multiple emissions are captured by this approximation.It is however worth pointing out that the scattering kernel, discussed in the preceding paragraphs, includes hard scatterings as well.A systematic study of improving the splitting kernel to include such rare interactions, which affects the early-time behavior of the rate, is left for the future.
For a time-dependent splitting kernel, a pertinent question is whether the relevant time should refer to the time between subsequent splittings or whether it should be counted from the beginning of the cascade.The former scenario is significantly more complicated, since it would involve a complicated interplay between subsequent emissions, and remains an open challenge (see Ref. [50] for an exploratory study).This problem applies both to the static and time-dependent media, see below.In Eq. (2.3), subsequent splittings are assumed to be independent and the global time is always counted from the beginning of the cascade up to its end.
In the static medium, it is easily identified from the medium-induced spectrum [13,14,41,43,51,52], and reads where ] and the time-independent part is where the (un-regularized) Altarelli-Parisi splitting function reads In the Bjorken expanding medium, the spectrum is given by Refs.[28,53,54], resulting in the rate [41,43] where τ 0 ≡ t 0 /t * , J α (•) and Y α (•) are the Bessel functions of the first and second kind, respectively, and ) At late times, the factor Re . . .behaves similarly to the static case, saturating at 1.In this case, the main effect of the expansion comes from the factor τ 0 /τ which arises directly from considering the static rate with the time-dependent q.Note that, in this case, the final time τ corresponds to τ ≡ (t 0 + L)/t * , where L refers to the path length in the medium.It is also interesting to note the additional dependence of the rate on t 0 in the Bjorken case.A study of the role of the t 0 has recently been done in Ref. [41] on the rapidity ratio and v 2 of the jets.In the subsequent sections, we shall further explore the possible physical implications of t 0 regarding the scaling laws.
We consider a plasma with the initial temperature T 0 = 0.4 GeV (N c = N f = 3) and the initial momentum to be p + 0 = 100 GeV.The coupling constant is fixed by ᾱ = 0.3, which corresponds to g ≈ 1.9869.For the above input parameters we get m D0 = 0.97 GeV, q0 = 1.81 GeV 2 /fm and, finally, t * = 11.0 fm.In addition, to facilitate a numerical solution, we impose the minimal value on the x variable at x min = 10 −4 .

Since we present the results in terms of
such that corresponds to the final longitudinal momentum, or in short energy, spectrum.
We solve the full evolution equation (2.3) for a collinear kernel using the initial condition of the evolution of the single-particle D(x, k, t = 0) = δ(1 − x)δ(k).For the results presented in the rest of the paper, we use the numerical solutions of the evolution equation that takes into account the transverse-momentum broadening as well as the splittings in an effective way as already demonstrated in Refs.[34,35,55] for the static, infinite media.To this end, we apply the Markov Chain Monte Carlo (MCMC) algorithm implemented in the event generator MINCAS, as described in Ref. [34], with the necessary extensions.These extensions account for the time dependence of the splitting kernels given in Eqs.(2.9) and (2.11) as well as of the collision kernel according to Eq. (2.6).Such a dedicated MCMC algorithm proved to be efficient in solving the evolution equation (2.3), even for the complicated Bjorken-expanding medium case.

Scaling in the energy spectrum
Now, integrating Eq. ( 2.3) over the transverse momentum k, we obtain the evolution equation for the energy distribution D(x, t) = d 2 k D(x, k, t) = x dN/dx of the gluon cascade [39,56]: in which the collision term has been integrated out to zero.In this work, we attempt to use the static in-medium kernel rates defined in Eqs.(2.9) and (2.11) for which the above equation can only be solved numerically.This was previously also analyzed in Ref. [57], albeit neglecting the finite-size effects in the rate, as contained in Eq. (2.11).Inclusion of the finite-size effects in the rate was done in Ref. [43] for the pure-gluon cascade, followed by the multipartonic cascades in Ref. [41].Here, we use the dedicated MCMC algorithm implemented in the MINCAS program [34] to solve the equations numerically.
Before we turn to the numerical evaluations, let us first re-visit the scaling laws in the evolution variable among different medium profiles.The quenching parameter q(t) for expanding medium is time-dependent.The single-gluon emission spectrum for different medium-expansion scenarios possesses scaling features for an average transport coefficient [28,44].However, as demonstrated in Ref. [43], the average scaling in q is valid only in the hard part of the single-gluon spectrum and not relevant for the soft part of the spectrum which contributes the most to the multiplicity, and therefore quenching [58], of emitted gluons.Instead, to establish the optimal scaling in the soft sector of the spectra, an "effective" quenching parameter has been identified [43].In terms of the effective inmedium evolution time L eff , this translates to where the temperature profile is given by Eq. (2.1).This establishes a relation between the length traversed in a static medium (on the left-hand side of the equation) and in an expanding medium (on the right-hand side), which corresponds to the equivalent amount of quenching.For the Bjorken model considered here, the relation reads If we assume L t 0 , we can write Since we are interested in the regime of multiple soft scatterings, we consider only the effective scaling.
The different profiles have been implemented in the Monte Carlo code MINCAS, and the results are shown in Fig. 1 for the evolution parameters given in Table 1.In the upper row, we plot the resulting spectra for a fixed length of the medium L = 4 fm and 6 fm, with two choices for the Bjorken initial time t 0 = 0.6 fm and 1 fm.Naturally, one obtains a much more pronounced evolution for the static medium as compared to the expanding-medium cases.However, we see the appearance of the turbulent-like cascade ∼ 1/ √ x for all the medium profiles.In the lower row of Fig. 1, the spectra have instead been generated at the equivalent effective medium length L eff ; see the third column in Table 1.One can observe the scaling occurring in the low-x part of the spectra for different medium profiles, thus re-confirming the findings in Ref. [43].These scaling features can also be found analytically within a simplified evolution equation, see e.g.Ref. [57].

Fully differential spectra
We now move to the full solutions of the evolution equation in the energy fraction x and the transverse momentum k T ≡ |k|.Actually, the full solution of Eq. (2.3) is 3-dimensional in the variables (x, k), but its dependence on the azimuthal angle in the transverse-momentum plane φ k is trivial, so it is integrated out.In order to facilitate comparisons between the different medium profiles, we shall henceforth consistently plot the distributions at their equivalent effective medium lengths L eff , see Table 1.The reasoning behind this particular choice for presenting our numerical results will become clear in a moment.Figure 1: The comparison of the gluon energy spectra for different medium profiles.The upper row shows the spectra evaluated at the medium lengths L = 4 fm (left) and 6 fm (right), respectively.The lower panels show the scaling of the gluon spectra for the scaling variable L eff corresponding to the two above values of L (see Table 1).
expected to be dominated by multiple soft-gluon scatterings at small k T , leading to the Gaussian profile.At high-k T , this becomes a power-law suppression due to rare hard medium interactions.
• In the low-x (x 1) regime, on the other hand, the k T distribution is narrower and approximately Gaussian.No distinct transition to a power-law behavior is observed for the range of plotted k T values.
These features are qualitatively in agreement with the discussion in Ref. [59] for the static media.There, the fully differential spectrum was postulated to factorize as D(x, k, t) D(x, t)P (k, t).Let us, therefore, parallel the qualitative observations made directly from the obtained numerical results in Fig. 2 with a discussion based on these analytical estimates.
radiative corrections.This distribution is defined where v(x, t) = N c n(t) d 2 q (1 − e −iq•x )σ(q)/(2π) 2 is the so-called dipole cross section, σ(q) ≡ dσ el /d 2 q being the medium elastic scattering potential.This distribution is characterized by two distinct regimes, see e.g.Ref. [19].At small k, we can approximate v(x, t) q(t)x 2 /4, leading to where Q 2 s (t) = t 0 dt q(t ) is the characteristic scale, sometimes called the saturation scale.In a static medium, Q 2 s = q0 L. This regime is therefore dominated by multiple soft scat-terings, leading to the Gaussian broadening.At large k, on the other hand, follows from the generic behavior of t-channel exchange.As already indicated in the formulas, this regime is easily generalized to expanding media via the time-dependent quenching parameter q(t) (or the density n(t)).Also, it becomes clear that the broadening of the leading fragments is significantly less effective when the medium is becoming dilute compared to when it remains at constant density.
At low-x, the situation becomes more complicated.The energy distribution can approximately be described by a turbulent cascade [39].The broadening is, however, not easily disentangled in this case due to the structure of the evolution equation (2.3).For a static medium and in the low-k regime, it was found [59] that P (k, t) again becomes Gaussian, however with a width given by k 2 ∼ √ xE q.This particular behavior comes about since the k T distribution in this regime is driven mainly by multiple parton splittings.Note also that the k T distribution is much narrower than in the large-x regime, in line with the results in Fig. 2. In angular variables, where in the soft limit we can approximate θ ∼ k T /(xp + 0 ), the width is enhanced in the small-x regime.In an expanding medium, we have not been able to obtain an analytically closed formula in these approximations.However, assuming that the broadening continues to be dominated by multiple splittings, we expect to see a similar amount of broadening in the soft sector for the static and expanding medium profiles when evaluated at the equivalent effective evolution time L eff .In the next steps, see also Sec. 5, we will study the behavior in the numerical data from MINCAS and reveal how wide the jet becomes in expanding media compared to a static one.We start with presenting the transverse-momentum spectra of the gluons as obtained from the evolution equation (2.3) in the case of the static medium and the Bjorken expan-sion.First, we consider the k T distribution defined as2 D(k T , t) = 1 0 dx D(x, k T , t).This corresponds roughly to the k T distribution of the average x fragments in the cascade, and is plotted in Fig. 3 for different media.Due to the Jacobian, the distribution is suppressed at small k T .At intermediate k T we recognize the characteristic Gaussian peak, followed by a strong, power-law suppression.We have plotted the distributions for the equivalent effective in-medium path lengths, see Table 1.Since L eff was calculated with the scaling of the soft sector in mind, these distributions, which are mostly sensitive to the large-x fragments, do not scale.Rather, we observe that the interactions occurring in the expanding medium are significantly less efficient in generating large-k T fragments than in the equivalent static medium.
As mentioned above, the simultaneous process of multiple splittings and scatterings, as encoded in Eq. ( 2.3), sets up a dynamical picture where one has to look for leading and sub-leading contributions of the radiative and collisional processes.At this point, it is useful to switch variables from the transverse momentum k T to the polar angle θ.This is achieved by the following transformation: where we have used the small-angle approximation k T = xp + 0 θ with θ being the polar angle 3 .It then naturally follows that the fully inclusive angular distribution is given by D(θ, t) = 1 0 dx D(x, θ, t).In order to explore the scales where different effects are at play, we separate the problem by looking at the x distribution in specific bins of the angle θ and vice versa.
To understand the behavior of soft gluons in the low-x regime, where the scaling estimates should work, we analyze the distribution of gluons in x for a particular angular range.This corresponds to the angular integrated energy distribution D(x, t; {θ max }) = θmax 0 dθ D(x, θ, t) . (4.5) In Fig. 4, we show the comparisons of the medium evolved spectra for different angular choices of θ max = 0.2, 0.4 and 0.6.It is worth pointing out that the fully integrated over θ spectrum (dotted lines) was analyzed in the literature for purely radiative cascades in the cases of infinite media [24,39,66,67], finite media [41,43,49] as well as including inmedium collisional elastic scattering [34,68,69].It was realized that a turbulent cascade, transferring energy from hard to soft modes, was responsible for the characteristic behavior in the low-x regime.
In Fig. 4, one can see the turbulent behavior in the region 0.3 ≤ x ≤ 1 for gluons within an opening angle θ max .Since the hard part of the spectra (x ∼ 1) is mainly driven  by collinear splittings, we observe that the spectra remain the same with increasing the opening angle to a much broader cone of θ max = 0.4 and 0.6.Thus, the hard and the intermediate energetic gluons are largely confined to a narrow cone size.On increasing the medium size to L eff = 6 fm (right panels), we observe the depletion of the peak at high x, resulting in the transfer of more softer particles to the intermediate and small x for all the angles.This is true for both the static and expanding finite media.The finite-size effect in the medium-induced emissions shows up around x ∼ 0.3 (a dip) which would otherwise be much flatter for the infinite media.Next, in the softer region (0.001 ≤ x ≤ 0.3), the energy distribution caused by the accumulation of soft gluons towards the medium scale leads to a significant broadening effect on increasing the angular region from the narrow one of θ max = 0.2 to the wider regions of θ max = 0.4 and 0.6.This region is driven mainly by radiative processes.Let us note that opening up the polar angle as well as increasing the medium size recovers more softer gluons for all the medium profiles.Interestingly, we observe the universality of the scaling between different medium profiles for all the jet opening angles as compared to the total angular integrated ones.In the following, we further probe the region of the soft gluons in angular space.
Next, we analyze the angular dependence of the spectra for all the media profiles in four different x ranges, restricting us again to the same equivalent in-medium path L eff for different medium scenarios.These are: 0.5 ≤ x ≤ 0.9 (hard) and 0.3 ≤ x ≤ 0.5 (semi-hard), which recover the gluons undergoing less in-medium quenching, 0.1 ≤ x ≤ 0.3 (semi-soft) and 0.01 ≤ x ≤ 0.1 (soft).For our purposes, we have chosen the angular region from 0.1 to 0.6 in θ to recover as much of the phase space of interest as possible.More importantly, ■ Bjorken, t 0 = 0.6 fm, L = 10.7 fm ■ Bjorken, t 0 = 1.0 fm, L = 8.0 fm 0.1 0.2 0.3 0.4 0.5 0.6 10 -6 0.001 ■ Static, L = 6.0 fm ■ Bjorken, t 0 = 0.6 fm, L = 21.0 fm ■ Bjorken, t 0 = 1.0 fm, L = 15.0 fm 0.1 0.2 0.3 0.4 0.5 0.6 10 -6 0.001 ■ Bjorken, t 0 = 0.6 fm, L = 10.In Fig. 5, for the hard and semi-hard regimes (upper panels) we observe that the angular distribution drops sharply for hard gluons as the energy-momentum conservation imposes a restriction on the available phase space for large angle scatterings and tends towards the expected Coulomb tail ∼ 1/θ 4 [42].In all the panels, we observe the transverse momentum broadening effect where the energy is re-distributed to larger angles due to elastic scatterings with the medium quasi-particles as well as subsequent gluon splittings.In a larger length of the medium (right column), we observe a greater magnitude of the broadening as the jet spent more time within the media.More strikingly, for the hard, semi-hard, and semi-soft regimes included, the broadening effect is significantly larger in the static medium than in the expanding scenarios.Only in the genuinely soft regime do we observe a similar amount of broadening.
The qualitative features of the results presented in Fig. 5 are as follows: • Collinear radiation causes depletion of the particles around the hard momentum sector x ∼ 0.5 -0.9.Subsequent splittings from hard particles to softer fragments undergo a successive broadening by which the energy is deposited in the soft sector of x ∼ 0.1 -0.01.
• The broadening effect is a combination of a decrease in momentum due to subsequent splittings as well as elastic collisions with the medium.The broadening effect is visualized as the energy deposited at small angles, θ = 0.2, being transported to larger angles, θ = 0.4 and 0.6, out of the jet cone.
• We observe insignificant broadening of the hard partons (a sub-leading effect compared to splittings), such that only the broadening near the medium scales x ∼ 0.01 contributes to the energy at larges angles.
• The broadening effect in the soft sector is a nearly universal effect with respect to the static or expanding medium and the finite medium-size corrections.It is also universal with respect to the starting time for quenching of the Bjorken profile.
In order to estimate the magnitude of the transverse-momentum broadening effect, in the next section, we estimate more precisely the amount of the in-cone gluon fraction in the soft as well as the hard limit of the gluon momentum fraction.

Soft versus hard gluons in a jet cone
In this section, we study the features of the jet in-cone energy loss for the opening angle θ and as a function of the momentum fraction x.As before, we restrict ourselves to  considering the same L eff in all three medium scenarios.In order to more precisely pin down the effect of medium expansion, let us define the fraction of the gluon energy contained inside a cone of the opening angle θ within the x-range x min ≤ x ≤ x max .This can be expressed as In the above equation, the denominator counts the total amount of energy contained in the chosen x-range for any angle.The numerical results are shown in Fig. 6, where we analyze the angular structure of the energy distribution of the gluons by dividing the medium evolved spectra into two regions corresponding to the large energy fraction: 0.5 ≤ x ≤ 1.0 (square markers) and the soft sector: 0.01 ≤ x ≤ 0.1 (round markers).
For the hard sector (square markers), the fraction of energy in expanding media is already quite close to unity for very small cone angles.Nevertheless, the static medium recovers most of the energy already at θ = 0.2.Hence, for phenomenological studies, one does not expect to be very sensitive to the details of medium expansion.
In the soft sector (round markers), the gluon cascade has developed a significant width and one needs to go to large angles θ 0.6 to recover the full energy content.However, for L eff = 6 fm we clearly observe that the cascade is narrower in the expanding medium than in the static one (the situation for L eff = 4 fm is less clear).The potential to be sensitive to the details of the medium expansion through the width of the angular distribution at small-x is interesting from a phenomenological point of view.

Conclusions and outlook
In this work, we have obtained spectra of jet particles that undergo scatterings and induced coherent splittings in a time-dependent medium, which is parameterized by the timedependent transport coefficient q and the corresponding splitting kernels [43] of Eqs.(2.9) and (2.11), corresponding to the static and Bjorken expanding medium, respectively.Numerical results have been obtained using the MINCAS Markov Chain Monte Carlo algorithm which calculates for the first time the corresponding fragmentation functions due to evolution via scatterings as well as medium-induced coherent splittings in the expanding medium.
Our findings can be summarized as follows: • The scaling approximations extracted from the singular limits of the purely radiative cascade work reasonably well for the medium-evolved spectra.
• We find that during the evolution of the energetic parton inside the media, hard partons remain effectively collinear and the momentum broadening is mainly caused by consecutive splittings rather than medium collisions and transverse-momentum exchanges.This type of behavior for hard partons is visible in the distributions of x and k T in Fig. 2, which follow at large x and small k T largely the Gaussian distribution (due to the predominance of multiple soft scatterings), while at large x and k T , rare hard medium interactions have a significant influence, leading to a power-law behavior.
• Subsequent decrease of the momentum of these energetic partons into the soft momentum sector causes momentum broadening by transverse-momentum exchanges with the medium through elastic scattering as well as subsequent splittings into softer fragments.However, the broadening due to subsequent gluon splittings in the softer sector contributes to the out-of-cone energy loss at larger angles.This kind of behavior becomes apparent from Fig. 4, where the differences between the full fragmentation functions D(x, t) and its respective contributions within jet cones of different sizes are shown in comparison.It appears that with increasing the jetcone angle, contributions of more and more soft particles are considered, while the behavior at large x has almost no contributions due to large angles.
• The scaling approximation minimizes the difference for the soft jet fragments (small x) and for the same value of L eff for the medium profiles as evident in Fig. 4. The harder (large x) jet fragments differ among each other for different medium profiles for the same scaling approximation.
• The fraction of energy within a x-bin is also analyzed to distinguish between different medium profiles.In Fig. 6, we observe persistent differences in the soft sector, implying that the cascades in the expanding media are still relatively more collimated than their counterparts in the static scenarios.
• In the quantities we have chosen to analyze, we have not observed any sensitivity to the hydro initialization time t 0 , in contrast to v 2 of jets [41], see also Ref. [70].
Naturally, for a similar in-medium path length in the static and expanding media, a leading parton will evolve much less due to the rapid diluting of medium density.This will lead to much fewer low-x gluons and a narrower profile in the polar angle θ.Comparing the distributions at the equivalent effective path length L eff should reduce this "trivial" effect.Nevertheless, our results indicate subtle differences between the developing profiles of the cascade in the static and expanding media in both the large-and small-x regimes.
One of the natural extensions of the results presented in this paper is introducing the transverse-momentum-dependent splitting kernels for the individual medium profiles to solve the evolution equations for both the quark and gluon-initiated jets.Secondly, one should include the effects of thermal re-scattering whenever the energy of the fragments becomes comparable to the local medium temperature, i.e. xp + 0 ∼ T , see e.g.Ref. [69].Furthermore, one can study the role of these modifications on the radial distributions, or the so-called jet shape functions, of the jet-quenching by comparison with the recent data from the CMS [71] and ATLAS [72] experiments at the LHC.However, these studies are beyond the scope of the present paper and will be reported in a separate upcoming work.

Figure 2 :
Figure 2: The D(x, k T ) distributions for the effective length L eff = 4 fm (left column) and 6 fm (right column), see Table 1.

Figure 3 :
Figure 3: The transverse-momentum distribution D(k T , t) for the effective medium length L eff = 4 fm (left panel) and 6 fm (right panel).

Figure 4 :
Figure 4: The distributions D(x, t; {θ max }) for the θ max choices of 0.2, 0.4 and 0.6 for effective medium length L eff = 4 fm (left column) and 6 fm (right column).The faded dotted lines show the corresponding distribution for the full angular region and correspond to the curves shown in Fig. 1.

Figure 5 :
Figure 5: The angular distributions D(θ, t) in different ranges of x for the effective medium length L eff = 4 fm (left column) and 6 fm (right column).The upper panels show the distributions for the semi-hard and hard gluons, while the lower panels show the ones for the soft and semi-soft gluons.

0 fm•■Figure 6 :
Figure 6: Comparisons of the fraction of the gluon energy in a cone P (θ, t; {x min , x max }) as a function of its opening angle θ within different x-ranges, x min ≤ x ≤ x max , for the effective medium length L eff = 4 fm (left panel) and 6 fm (right panel).

Table 1 :
Values of the initial time t 0 , the effective length L eff , and the actual length L of the evolution for different types of the medium used in Monte Carlo simulations.