Higher-order corrections to heavy-quark jet quenching

We calculate higher-order corrections to the quenching factor of heavy-quark jets due to hard, in-medium splittings in the framework of the BDMPS-Z formalism. These corrections turn out to be sensitive to a single mass-scale $m_\ast = (\hat q L)^{1/2}$, where $\hat q$ is the medium transport coefficient and $L$ the path length, and allow to draw a distinction between the way light, with $m<m_\ast$ (in contrast to massless $m=0$), and genuinely heavy, with $m>m_\ast$, quark jets are quenched in the medium. We show that the corrections to the quenching factor at high energies are double-logarithmic and qualitatively of the same order as for the massless quark jet.


I. INTRODUCTION
Jets are formed in the process of soft and collinear QCD radiation that results in a spray of collimated hadrons and energy deposition in the detector [1]. In heavy-ion collisions, partons traverse a hot and dense nuclear medium that leaves an imprint on the subsequent jet formation, for reviews see [2][3][4]. Currently it is widely accepted that the BDMPS-Z formalism of radiative energy loss [5][6][7] describes the propagation and multiple scattering of quark and gluon jets in the nuclear QCD medium that is produced in heavy ion collisions at LHC, for a review see Ref. [8].
Jet quenching is a multi-scale problem. Even for massless partons, there is an convoluted interplay between the intrinsic jet scales, such as the mass of the jet, and the scales of the medium, including typically the medium transport coefficient and the medium size.
For quarks, the non-zero mass introduces another scale. It is well known that the collinear divergence is regulated by the characteristic dead-cone angle where m is the mass of the heavy quark and E ≡ p T is the jet energy. As a consequence of the strong suppression of gluon radiation inside the dead-cone angle, heavy-quark jets fragment differently from jets originating from their massless counterparts or from gluons [1,[9][10][11].
Radiative energy loss was calculated by BDMPS-Z for massless partons [5][6][7]. It was first pointed out in Ref. [12] that the quenching of massive quarks would be different from massless ones because of the dead-cone effect. The resulting restriction of the phase space for radiation, and hence energy loss, leads to a systematically smaller suppression of single-inclusive hadron spectra the larger the mass of the constituent quarks. This was followed up by a more thorough analysis in [13,14], where it was shown that the heavy quark quenching factors get further corrections when the correct phase space constraints are taken into account. For results within the limit of dilute media, see also [15,16].
In summary, based on radiative processes alone one expects a smaller rate of emissions off massive quarks compared to massless ones, that brings about a mass-hierarchy of the suppression. In contrast, low-p T heavy mesons have a similar modification as the pions [17][18][19]. This has prompted many investigations of additional elastic energy loss processes, for a review see [20]. It is however worth keeping in mind that the final suppression of heavy mesons and heavy flavor jets depends also on the details of the partonic cross sections and the problem is still an open one.
While most of the previous contributions has focussed on small p T , where the cross sections are the largest, we will mainly focus on the genuinely high-p T regime where perturbative corrections play a crucial role. This regime is within the reach of the experiments at LHC, see e.g. [19,21] . Recently higher-order corrections to the quenching of independent, massless quark/gluon jets were calculated [22,23]. The results demonstrate how these contribution lead to the enhanced quenching of massless quark/gluon jets as compared to single partons. The role of in-medium jet splittings and their color coherence properties has also been emphasized in other contexts, see e.g. [24,25]. Consequently, it will be of interest to extend the previous efforts to include mass effects. In this work, we consider higher-order corrections to the quenching of a heavy-quark jet, i.e. a jet formed as a result of the fragmentation of a leading massive quark.
Our main result is that higher-order corrections lead to an enhanced suppression for heavy-quark jets relative to the leading BDMPS-Z result, corresponding to the quenching of a single parton. The magnitude is determined by the phase space available for the radiation of hard gluons within the jet, and is of similar magnitude as for quark/gluon jets in general. However, due to the restricted phase space determined by the dead-cone (1), the mass sets the scale where significant deviations between massive and massless jets can be observed. We identify a critical mass scale that permits to observe such discriminating features in the high-p T regime.
The paper is organized in the following way. In Sec. II, we introduce the generalized quenching weight and discuss its expansion in terms of the strong-coupling constant.
We calculate the radiative energy loss due to multiple, soft BDMPS-Z radiation off a single heavy quark and a heavy quark-gluon dipole in Sec. III, and obtain the evolution equations and expressions for related quenching factors. The details of the calculations of the associated spectra and rates are given in Appendices A and B, respectively, where the basic formulae for the interference contributions to antenna radiation are derived in detail.
In Sec. IV, we finally map out the logarithmic phase space for higher-order corrections and present explicit expressions for the collimator function of heavy-quark jet together with numerical results. We summarize our results and give a brief outlook in Sec. V.

II. GENERALIZED QUENCHING WEIGHT
Assuming small energy losses in the medium, p T , and accounting for a steeply falling hard spectrum, the spectrum of heavy-quark jets in heavy-ion collisions can be written as dσ dp 2 where dσ 0 /dp 2 T is the Born-level jet production cross section, P( , L|m) is an energy-loss probability distribution associated with a massive particle and L is the medium length (below we shall suppress the arguments L and m, unless it is unclear from the context).
The jet suppression factor Q(p T ), introduced in the second step, is the Laplace transform of the energy loss distribution P( ), i.e. Q(p T ) ≡P(n/p T ) = ∞ 0 d e −n /p T P( ) [26], where the effective power of the steeply falling spectrum is n = d d ln p T ln dσ 0 dp 2 T [12,27]. The jet suppression factor permits an expansion in the strong-coupling constant that accounts for the energy loss of in-medium jet splittings, The first term in the expansion is the quenching of the jet total charge which, for a heavyquark initiator, is given by Q (0) (p T ) =P q (n/p T ). This distribution is dominated by soft gluon radiation that transfers energy from the jet axis to large angles.
The resummation of higher-order terms leads to an additional suppression factor which was referred to as the "collimator" function in Ref. [23]. These corrections correspond to the energy-loss of composite, partonic systems created inside the medium during the jet formation. Hard splittings in the jet cone can be described by vacuum splitting functions.
Hence the next-to-leading correction to the jet quenching factor, that involves the (real and virtual) emission and subsequent quenching of an additional gluon [23], takes the form where P gq (z) is the Altarelli-Parisi splitting function and Q gq (θ, p T ) is the quenching factor of a composite quark-gluon system propagating in the medium [22]. This equation holds whenever the splitting takes place early in the medium. This enforces the formation time, , to be short compared to that of any process in the medium, in particular the medium length L. As will be discussed in more detail later, an important time-scale is the so called decoherence time t d which corresponds to the time when a dipole of size x ⊥ ∼ θt, characterized by its opening angle θ, is resolved by medium fluctuations.
The characteristic wave-length of the latter can be estimated via diffusive broadening as λ ⊥ ∼ (qt) −1/2 . The two length-scales become of the same order at t d ∼ (qθ 2 ) −1/3 . In this limit, t f t d L, the splitting process completely factorizes out on the level of the cross section [28] and effectively forms a color-charged antenna. This composite system undergoes further medium-induced radiation in the medium that turns out to be sensitive to its opening angle [22]. However, there can also be strong cancellations between the two quenching factors in the squared brackets in (4) for small-angle emissions, when t d > L or θ < θ c , where θ c ∼ (qL 3 ) −1/2 , due to interference effects.
Higher-order corrections naturally follow a similar logic, becoming sensitive to more complicated radiation patterns. In the large-N c limit the picture is simplified further, since a jet in this case can be decomposed into a set of mutually independent color-singlet dipoles, whose radiation is added to that of the total charge radiation [29]. The quenching of the total charge can therefore be factorized out, and the total quenching factor becomes, where C(p T , R|m) is the collimator function that accounts for the quenching of higher-order (real and virtual) jet emissions. The resummation of such emissions takes, in the general case, the form of a non-linear evolution equation for the collimator function but, in the limit of strong quenching, one can neglect all real emissions and resum the virtual terms, i.e. the second term in the squared brackets of Eq. (4). We will discuss the collimator in more detail in Sec. IV. In the remaining part of the paper, we will describe the radiative quenching of a heavy-quark system and identify the relevant time-scales that play a role in this problem in order to compute and resum these corrections. We focus on the high-p T regime and (relatively) large quark masses, where elastic energy losses, see e.g. [20], can be neglected. Our results at low-p T are therefore not completely realistic. However, we emphasize that for genuinely heavy quarks the high-p T regime (meaning p T 20 − 50 GeV, see below for more details) is relevant for the phenomenology of heavy-quark jets.

III. COMPUTING THE QUENCHING FACTORS
In this section we compute the quenching weights, that is energy loss probability distributions that resum multiple soft, gluon radiation responsible for transporting energy from the leading particle to large angles. As mentioned in the Introduction, it will be convenient to work directly in Laplace space, defined as where the contour C runs parallel to the imaginary axis in the complex-ν plane, Reν = const., to the right of any singularity ofP(ν, L|m). To recap, P( , L|m) acts as a probability distribution for radiating gluons that in total carry an energy off a particle with mass m after propagating through a medium of length L, and the quenching factor Q(p T ) =P(n/p T , L|m). It will be convenient to define a "regularized" splitting rate, where γ i (ν, t) ≡ γ ii (ν, t) and we have already anticipated the possibility of interference contributions between two different particles labeled "i" and "j" that refer to quarks, antiquarks or gluons. Here Γ ij (ω, t) is the rate of (interference) emissions in the medium, where the soft gluon is emitted by a parton i and absorbed by the parton j in the complex conjugate amplitude. We derive the generic interference spectrum off a color-charged antenna in Appendix A, and derive concrete expressions for the direct and interference rates in Appendix B within the multiple-soft scattering approximation.
A. Quenching of a single parton Let us start by considering a single propagating particle. Medium interactions can enhance the probability of gluon emissions. In Laplace space, the resummation of soft, medium-induced gluons takes the form of a rate equation, with initial conditionP(ν, 0) = 1, whose solution is simply given byP i (ν, L) = e L 0 dt γ i (ν,t) . For a massive quark, the emission rate of soft gluons was derived in Eq. (B11), and reads where ψ 0 (x) is the digamma function,ᾱ ≡ α s C F /π and the expression in the squared brackets is a function of the scaling variable ζ ≡ ω/ω DC , where The expression (9) is valid only for ζ < 1. It turns out the spectrum is strongly suppressed at ω > ω DC , where we also observe negative contributions owing to the treatment of the high-energy behavior which goes beyond the leading-logarithmic accuracy of our calculation [30]. In order to avoid these unphysical contributions, and retain the information about the physical scales, we approximate the rate by In Laplace space, this becomes where erf(x) is the error function. This is qualitatively similar to what is obtained in Ref. [12], although the precise form of the cut-off at ω DC determines a numerical constant in front of the second term in the brackets (according to the authors of Ref. [12], this factor is ∼ 1.5). The corresponding rates for a massless quark is found by taking m → 0 and for gluon by futher replacing the color factor C F → N c , e.g.
While these rates are time-independent in the limit of soft gluon emissions, this is violated at large energies. This approximation breaks down for emissions with formation times of the order of the medium length, corresponding to a critical energy ω c ∼qL 2 that brings about a power-like cut-off of the spectrum and, therefore, the rate as well.
The constraint from the dead-cone angle is stronger than this absolute limit whenever ω DC < ω c which, in turn, implies that Θ 0 > θ c . This marks the regime where the mass of the quark should start affecting the general properties of radiative energy loss that is dominated by LPM interference effects.
B. Quenching of a two-parton system (color-charged dipole) Let us now turn to the higher-order corrections to this picture, that arise from a quarkgluon antenna propagating in dense QCD media. Considering for the moment the energy loss of a quark-gluon dipole that is formed quasi-instantaneously after the hard vertex, in Laplace space the joint energy loss distribution factorizes in the large-N c limit into the product of energy loss off a total charge (triplet) and a color-singlet dipole, whereP q 0 (ν, t) ≡P q 0 (ν, t|m = 0) and we have decomposed the gluon into a quarkantiquark pair g = (q 0 , q 1 ) (where q 1 is an antiquark). Recall that only the quark that forms part of the dipole q = q 2 is massive. Note that the quenching of the total (quark) charge is not sensitive to the mass of the initial particle, since it is inherited from the radiated gluon. Instead, the mass controls the energy loss of the additional irreducible singlet P q 1 q 2 (ν, t|m) ≡P sing (ν, t|m). We will confirm below that, in the completely decoherent limit, the mass will be associated with the total color charge, as expected.
The two factors in (15) satisfy two separate evolution equations. First, the singleparticle quenching is given by Eq. (8), where the splitting rate γ q is explicitly given by (14), keeping in mind that this fictitious quark is massless. The singlet, dipole quenching weight is determined by solving the differential equation [22] ∂ tPsing (ν, The initial condition at t = 0 (corresponding to the time when the antenna was formed in the medium) is again trivial,P sing (ν, 0|m) = 1. The direct and interference rates were derived in Appendix B, and are given by The two direct terms correspond to emissions off the two legs, and similarly the interference terms correspond to emitting a gluon from one leg and "absorbing" it (in the complex conjugate amplitude) on the other (see Fig. 4 for details). For the singlet dipole Note that Eq. (16) contains a inhomogeneous term arising from the possibility of interferences between the dipole consituents. The interference spectra have a more complex structure since they involve both color and quantum decoherence processes [31][32][33]. The interference spectrum associated with a massive dipole is given explicitly in Eq. (A9), see also [34], and evaluated in the multiple-soft scattering approximation in Eqs. (B20) and (B21). As discussed further in Appendix B, color decoherence is related to the survival probability of a color-singlet dipole and is explicitly factorized out in the so-called decoherence parameter S 2 (t) in (16). This factor is responsible for the previously introduced time scale t d ∼ (qθ 2 ) −1/3 . At late times, t > t d , the dipole decoheres and the particles can radiate independently. Conversely, a coherent splitting corresponds to the situation when t d > L and the pair remains coherent during the passage through the medium. This applies to the regime of small angles, θ < θ c where θ c ∼ (qL 3 ) −1/2 .
In our formulation, the medium-induced interference rates (in energy-space) are themselves suppressed at a time-scale t quant ∼ (θ 2 ω) −1 . This scale can nevertheless be neglected is the branching angle of medium-induced gluons. Since energy loss is governed by soft gluons, ω ∼ᾱ 2 ω c , that parametrically go to large angles, in particular out of the jet cone θ br (ω) ∼ᾱ −3/2 θ c > R > θ, this implies that t d < t quant [22,32]. Since the multiplicity of hard gluons is small, and correspondingly their contribution to energy-loss is small [27] we will neglect such quantum effects in the following. Hence, for our purposes, i.e. at Eqs. (B22) and (B23). This property is independent of the mass.
The solution to the rate equation can be written symbolically as P sing (ν, L|m) =P 1 (ν, L)P 2 (ν, L|m) The extension of the time-integral of the second term is limited by the shortest timescale where interferences are suppressed. In the leading-logarithmic approximation it is sufficient to consider only large-angle radiation where, parametrically, the energy radiated via medium-induced gluons leave the jet cone. In this case the integral is limited by t < t d , as discussed above, and the singlet distribution can then be approximated bỹ where we have reinstated the mass dependence. Hence, the decoherence time acts as a "delay" for when energy loss processes start affecting the irreducible dipole and, in the limit t d L, the dipole constituents decohere early in the medium and lose energy independently along the whole medium length.
To summarize, the delay effect is strictly associated with the color dynamics of the dipole and, since this involves the shortest relevant time-scale, does not depend on the mass of the constituents. It might, at first look, seem strange that the mass-effect on quenching is delayed although it is intimately linked with the quark-initiator and, hence, the total charge. For instance, considering long decoherence times, t d ≥ L, applying to small-angle emissions θ ≤ θ c , the color-charged antenna is quenched as a massless quark, rather than a massive one. This effect gives rise to a mismatch between real and virtual emissions at small angles. This turns nevertheless out to be a sub-leading effect, see Eq. (29).

IV. HEAVY-QUARK COLLIMATOR FUNCTION
A. Quenching of total charge The first term in the expansion in (3) corresponds to the quenching of a single, massive quark. After implementing the result in (13), we find that where the first term corresponds to the quenching of a massless color parton, while the second is a mass-dependent enhancement factor. The mass-independent term implies that the regime of strong quenching of massless quarks, − ln Q (0) ∼ O(1), is given by up to numerical factors, where ω s ∼ᾱ 2q L 2 is a soft scale corresponding to large multiplicity of medium-induced emissions, N (ω < ω s ) > O(1). We can rewrite the quenching factor as Q (0) (p T |m = 0) ≈ exp − √ πN (p T /n) [27], which is interpreted as a Sudakov suppression factor for medium-induced gluons with energies ω > p T /n.
The regime with an additional strong enhancement of heavy compared to massless quarks arises for where θ s ∼ θ br (ω s ) ∼ᾱ −3/2 θ c . This condition is equivalent to demanding that θ s < Θ 0 , which implies that the regime of multiple, soft gluon emissions is cut off by the deadcone angle. So there exists a regime of strong massless-quark quenching with additional strong effects from heavy-quark quenching whenever for masses m < nᾱ 1/2 Q s where Q s ∼ (qL) 1/2 . To limit the scope of our qualitative analysis, we will assume that the index of the steeply falling spectrum n, combined with the medium parametersq and L, is large enough to work in the regime of strong quenching effects.

B. Scale analysis
Let us now turn to the next-to-leading correction (4), coming from the emission of an additional hard gluon off the initiating heavy quark early in the medium. We will currently focus on the leading-logarithmic contributions, leaving an analysis of sub-leading logarithmic contributions for the future. In this context we only consider strong ordering of scales and we will therefore not typically keep track of numerical factors that are anyway beyond the precision of this analysis. We will also work in the large-N c limit, where we can exploit the factorization of the color-charged dipoles, as in Eq. (15).

Massless quarks
Before turning to effects related to the mass of the jet particles, let us summarize the scale analysis for massless partons. In terms of angles, we have several characteristic scales: the jet radius R, the minimal medium resolution angle θ c ≡ θ br (ω c ) and the typical angle for soft gluon emissions θ s ≡ θ br (ω s ). Note that the two medium scales are parametrically separated by the smallness of the coupling constant, θ s ∼ᾱ −3/2 θ c .
If R > θ s the energy loss of jet is small since all radiated BDMPS gluons remain inside the jet. On the other hand, if θ c > R, i.e. the jet angle is less than decoherence angle, the propagation of the jet is not influenced by subjet structure, and it is equivalent to propagation of a total color charge, i.e. one parton through the medium. The typical ordering that interests us, where the quenching could be substantial and where higherorder effects are non-trivial, is therefore In what remains, we will assume that this hierarchy holds and, besides, that it is also the phenomenologically most relevant one.
However, note that the minimal angle θ c only is relevant for high-p T jets, p T > ω c .
Conversely, for p T < ω c the decoherence time is necessarily always shorter than the medium length, t d < L. In this case there is still the possibility for a regime of short formation times, t f < t d , but in this case this condition implies that θ > θ d , where For massless quarks, this regime is double-logarithmic in the jet scale [23], see below, but the p T -range is automatically limited by ω c . The window for a regime of short formation times closes whenever θ d = R, or p T ∼ (q/R 4 ) 1/3 . In the following, we will therefore distinguish between high-p T (with p T > ω c ) and low-p T (with p T < ω c ) jets.

Massive quarks
For massive quarks, the dead-cone angle (energy, etc.) introduces another physical scale to the problem. For a finite dead-cone, QCD radiation is no longer genuinely collinearly enhanced which necessitates a scale-dependent scheme to properly include mass-effects for resummed observables, see e.g. [1]. We will only stick to the leading-logarithmic approximation and only consider emissions (4). Hence, for R ∼ Θ 0 the heavy-quark jet only contains a single quark.
Comparing the mass scale to other relevant medium scales, in particular comparing Θ and θ c , can become involved because of the p T -dependence of the former. In order to organize the discussion, it will be useful to introduce a critical value of the mass, namely From now on we will call quarks with m > m * genuinely heavy, and quarks with m < m * for light (in contrast to massless, m = 0). We have also summarized the discussion about the relevant scales in Table I.
For heavy quarks, the dead-cone angle becomes comparable to the coherence angle, Rewriting the same condition, this happens at a critical energy p T = mθ −1 c ∼ m(qL 3 ) 1/2 . Hence, we expect the heavy-quark jet quenching to deviate from the light-quark jet quenching at a scale that is parametrically larger, by a factorᾱ −3/2 , than the soft scale identified for the quenching of the total charge, cf.
Eq. (23). In other words, while the quenching of a single heavy quark starts deviating from the massless one at relatively low p T due to the enhancement factor in (21), a jet initiated by a heavy-quark should start deviating from the behavior of a massless quark or gluon jet already in the high-p T regime, since by definition m(qL 3 ) 1/2 > ω c . Considering high-p T single-inclusive mesons to be proxies of single-parton dynamics, see e.g. [35], this analysis therefore predicts a different behavior of heavy-quark jets and heavy-quark mesons over a large range in p T . We will come back to a possible experimental signature for this effect in Sec. IV D.
For light quarks, Θ 0 = θ c for p T ≥ ω c and we have instead to consider the low-p T regime, i.e. p T < ω c . In this regime, the condition t f < t d implies that θ > θ d which is estimated in (25). Therefore Θ 0 = θ d when p T ∼ m 4 /q, which can be considered a relatively soft energy-scale (comparing it to the soft scale in (23) gives m ∼ᾱ 1/2 (qL) 1/2 which is compatible with the prior assumption about the smallness of the mass). We conclude therefore that the light-quark jets behave similarly to massless jets, as far as the higher-order corrections go, and start deviating from this behavior only when Θ 0 ∼ θ s , where the quenching of the total charge gets suppressed. This follows very closely the trend of single-parton, or single-inclusive meson, quenching.
The corresponding kinematical Lund planes for high-and low-p T heavy-quark jets are illustrated in Fig. 1, where we have spanned the plane in the logarithmic variables 1/z and 1/θ. At fixed coupling, the plane is equally filled with splittings with probability 2ᾱ, up to a color factor. The two diagonal lines, with slopes −4/3 and −2, delineate the conditions t f = t d and t f = L, respectively. The area between the two lines corresponds to in-medium radiation with t d t f < L, or k ⊥ √q ω, which is strongly influenced by medium interactions and broadening.
The high-p T regime, p T > ω c , is plotted on the left side of Fig. 1, where we have marked the location of the critical angle θ c with a (red) dotted line. Similarly, the dead-cone angle is marked, and corresponds to an energy scale ω DC at t f = t d . The low-p T regime, p T < ω c is conversely plotted on the left in Fig. 1. One observes immediately that the critical angle θ c is replaced by θ d . In both figures we have assumed that the dead-cone is appreciable, i.e. Θ > θ c for p T > ω c and Θ > θ d for p T < ω c , and marked out the phase space available for hard, in-medium splittings of the heavy-quark.

C. Higher-order contributions to quenching
The analysis in the preceding section allows us to calculate the higher-order contributions to jet quenching. Isolating the quenching of the initiating parton, that corresponds to the total color charge of the jet, into an overall pre-factor, see Eq. we see that the first-order correction the collimator function is where we have treated the splitting vertex in the leading-logarithmic approximation and adopted the notations of the previous section.
It is worth pointing out two limits of this equation. For t d L, we can neglect the decoherence times in the real term, i.e. the first term in the squared brackets, to obtain which, when taking m → 0, is equal to the contribution of massless quark quenching.
When Θ 0 > max(θ c , θ d ) the angular phase space is more restricted for heavy-quark jets, and therefore we expect a relatively smaller impact of the collimator function than in the massless case. Here it is worth pointing out that the quenching factor on the righthand side of (28) arises due to the quenching of the additional (massless) gluon since, at large-N c , Q 2 q (p T ) = Q g (p T ), which is a generic property of Sudakov suppression factors. Before continuing, we point out a new contribution in the small-angle limit in the regime that is unique to massive-quark jets. It appears for t d > L, or θ < θ c , relevant for p T > ω c , where Q(p T , L − t d |m) = 1. We are left with However, Θ 0 < θ c for p T > mθ −1 c , which leaves the factor Q(p T )/Q(p T |m) − 1 ᾱO (1), and therefore the contribution in this regime is sub-leading ∼ O(ᾱ 2 ). We will therefore altogether neglect this regime when working in the leading-logarithmic approximation.
Let us now evaluate the next-to-leading contributions for massless, light and heavy quarks. For completeness, we repeat here the resulting collimator function for massless quarks at first order in α s , that reads [23] Turning now to the new results, for light quarks we obtain and for heavy quarks we get instead, Equations (30) Let us briefly comment on further contributions to the collimator at higher-order (nextto-next-to-leading, and higher). Examining the structure of Eq. (15), one realizes that the dipole that "contains" the heavy-quark is distinct from further dipoles in the sense that it is massive while further dipoles, originating from other gluon emissions, are massless. However, as discussed in detail above, this distinction gives rise only to sub-leading corrections and for our purposes, having separated out the specific quenching factor of the originating parton (that also carried the total color charge), it is adequate to treat all dipoles on equal footing.
The problem then reduces to the massless case with a modified phase space, as detailed above. The resummation of higher-order contributions to the collimator involves solving a non-linear evolution equation and was derived in Ref. [23]. It goes beyond the scope of our investigation to solve this equation here for the massive case. Furthermore, since we are interested in a relatively modest p T range in order to be sensitive to the dead-cone, the phase space is limited and the first, non-trivial term should provide a good estimate of the effects. Using the same arguments as in [23], we therefore expect that the full collimator function is well approximated by In particular, the strong quenching limit, Q(p T ) 1, returns the correct exponentiation of the virtual terms and also the fixed point at p T → ∞, where Q(p T ) → 1 and therefore C(p T ) → 1, is reproduced.

D. Numerics
To emphasize the effects of higher-order contributions we propose the following phenomenological quantity, which is a ratio of nuclear modification factors of heavy-quark jets to heavy-quarks.
Within our approximations, this ratio is simply the collimator function for massless and massive quarks J AA (p T , R|m) ≈ C(p T , R|m), where we have utilized that R jet AA (p T , R|m) Q(p T , R|m) and R meson AA Q q (p T |m). The latter can be justified in the sense that, although the quenching factor of a heavy-quark is not a direct measurable, it is closely related to the quenching factor of the corresponding heavy-meson [35]. The underlying idea builds on the assumption of a similar path-length dependence for the two observables.
Before we then present our numerical calculations, it is worth emphasizing our approximations in computing the single-quark quenching factors based on the soft-gluon approximation to the full BDMPS-Z spectrum. It was already pointed out in Ref. [27], that the sub-leading logarithmic and numerical factors play an important role for computing the right order of quenching effects and this was also adopted in [12]. For the moment we focus only on the effect of massless quark quenching, which enters the dynamics of the collimator functions, cf. (4). Keeping these corrections, the massless quenching factor in The sub-leading terms result in a faster approach of the quenching factor to unity. We have plotted the quenching factors for massless quarks that enter the calculation of the collimator function in Fig. 2, where the parameters were chose as described below and the uncertainty arising from modeling radiative energy loss is marked with the shaded region.
As we have done throughout, we will assume the medium to be static and described by averaged parametersq and L. We have chosenq = 1 GeV 2 /fm and L = 2.5 fm, and chosenᾱ = 0.15 [36]. These are are qualitatively in the same range as the values obtained in more sophisticated extractions from comparisons to experimental data. Furthermore, we compute the collimator function R = 0.4 jets and the power of the steeply falling heavy-quark spectrum n = 5, that was extracted from a fit of the jet data [37]. With these choices m * 1.5 GeV, and therefore the charm quark (m = 1.3 GeV) can be considered light while the bottom quark (m = 5 GeV) is heavy.
We have adopted the approximations and the two ways of estimating the associated, massless quenching factors as discussed above in computing J AA (p T ) in Fig. 3. As a result, the ratio is simply given by the collimator functions that have been explicitly computed to On the left we have used only the leading term of soft-gluon approximation of the BDMPS-Z spectrum; on the right side, we included numerical corrections to the quenching, as in Refs. [12,27].
next-to-leading order in Eqs. (30)- (32), with the associated quenching factors computed in Fig. 2, and whose full resummation is given in Eq. (33). The solid, black (dashed, organe) curves represent the collimator function for charm (bottom) quark jets. On the left side, we used the massless quenching factor associated with the first leading term of the BDMPS-Z soft-gluon spectrum while, on the right side, we included the sub-leading logarithmic and numerical factors, as was done in Refs. [12,27]. It is clear from Fig. 3 that the order of magnitude of the effect does depend on the inclusion of these sub-leading BDMPS-Z corrections since the full expression (35) approaches unity much faster than the leading term alone. However, the characteristic p T where the collimator function of charm and bottom starts deviating on the level of ∼ 10% is roughly the same in the two cases (for our choice of parameters, this corresponds to roughly p T ∼ 50 GeV).

V. CONCLUSIONS
We have calculated a subset of higher-order corrections for massive-quark jet propagating in the quark-gluon plasma that are enhanced by logarithms of the jet energy. We considered the corrections due to rapid split of the leading particle into hard dipoles well within the medium. This contribution also plays an important role in the context of jet substructure [24]. We have shown that these corrections lead to the enhancement of jet energy loss, and consequently to the decrease of the jet quenching factor, like it was found for the case of massless quark/gluon jets previously.
The additional suppression can be factorized into a collimator function, that enhances the quenching factor associated with the leading particle and is here evaluated in the leading-logarithmic asymptotics. We have demonstrated that these corrections are essentially determined by phase space restrictions available for dipole creation -on the one hand related to the criterium of early splitting, in particular with formation times t f < t d , and, on the other hand, at angles larger than the associated dead-cone angle. These semi-quantitive estimates show that these corrections are of the same order of magnitude as for massless quark at high-p T , contrary to leading order results [12], where the substantial difference between quenching factors of massless and heavy quark jets was found at relatively low p T . The reason is that in significant part of the parameter space the heavy and massless quark corrections are just the same, and in the remaining region they only differ by the argument of the logarithm.
Our main results are given by Eqs. (31) and (32), where we demonstrate that, from the point of view of higher corrections, we can divide the quarks into light, with m < (qL) 1 /2 and heavy, m > (qL) 1 /2 . For light quarks the corrections are exactly the same as for the massless quark, see Eq. (30), and start to be weakly mass dependent only for rather small energies of order p T ≤ m 4 /q. For heavy quarks, given in Eq. (32), already at large p T > ω c the dead-cone effect start to play role.
Our work shows that antenna corrections are essentially the same for massless and heavy quarks, leading to correct from the experimental point of view decrease of quenching factors of heavy quark jets relative to the Dokshitzer-Kharzeev result [12]. Consequently they had little influence on the problem first raised in [12], namely that, contrary to calculation focussing on radiative energy loss, there is not much difference between heavy and massless quark quenching factors experimentally [17,18] (although the situation can improve in the future [19]). Nevertheless, since the BDMPS-Z calculation is the leading mechanism available for jet energy loss at high-p T it is of great interest to test it for different parts of the parameter space.
We have also derived the spectra and rates related to a color-charged, massive dipole We did not include the restrictions on the phase space that were found to be important for calculation of massless and massive quark quenching weights, see [38] and [13,14]. We do not expect that these restrictions will make qualitative influence on our results. The reason is that the corrections that we calculated were dominated by small frequencies regime outside the dead-cone while the corrections discussed in [13,14] are mostly important for frequencies in, or close to, the dead-cone region. Besides, the analysis of these constraints needs detailed numerical investigation.
Our results are valid in the leading-logarithmic approximation at high energies. Going to higher logarithmic precision and for applications at very high energies we expect the corrections to satisfy nonlinear integral equations similar to the one that was suggested in [23]. In this context, the effects of the dead-cone suppression in secondary heavy-quark production, i.e. from the splitting of a gluon into a massive qq pair, are still largely unexplored. For a clearer interpretation, these contributions could perhaps be suppressed using state-of-the-art grooming techniques, employed in [39].
The results worked out here demonstrate the factorization between the leading term in the quenching weight, calculated in [12] and higher order corrections, and introduces a new way to address jet observables involving massive quark, e.g. the "leading particle effect" [40][41][42], that we plan to address in forthcoming works. From a phenomenological point of view, it will be interesting to study the ratio of quenching factors of heavy quark jet and heavy mesons as a check of the current approach, cf. Eq. (34).

ACKNOWLEDGMENTS
KT is supported by a Starting Grant from Bergen Research Foundation ("Thermalizing jets: novel aspects of non-equilibrium processes at colliders") and the University of Bergen.
BB was supported by the Israel Science Foundation (ISF) grant no. 2025311. The authors thank the CERN Theory Department for hospitality, where most of the work was done.
Appendix A: Derivation of the interference spectrum Let us derive the rates of gluon emissions of massless and massive particles. Consider a parton splitting at very short time-scales that quasi-instantaneously forms a dipole inside the medium. Using conventional arguments in perturbation theory, one can thereafter define a subsequent emission spectrum off this system. For further applications, it will be sufficient to calculate the interference spectrum between the two constituents of the dipole since direct emissions can be recovered by setting the dipole opening angle to zero.
For concreteness, and in order to clarify the color structure of the process, consider a color-charged dipole originating from a q( p 0 ) → q( p 1 ) + g( p 2 ) splitting, where we have identified the momentum flow of the splitting [43]. We assume, for the time being, that the gluon is massive; this allows us to easily generalize the formula for arbitrary dipoles.
The amplitude describing the emission of a soft gluon inside a medium of size L is the sum of two terms that correspond to the possibility of being radiated by dipole constituents.
They are given by up to factors that cancel in the cross section and where the fundamental color matrix t c jk accounts for the color conservation of the system (for completeness, we recall that trajectories. In contrast, a soft gluon emission in the medium can experience momentum broadening due to transverse momentum exchanges with medium constituents. These interactions are encapsulated in the dressed propagator G ab , that is given by where the trajectory of the gluon r is explicitly time-dependent. We proceed now with the calculation of the interference spectrum. It was first computed for a color-charged dipole involving a massive quark in [34] but an explicit expression is not available in the literature. It is therefore meaningful to rederive the spectrum in more generality here. The spectrum involves in total four terms, corresponding to the different time orderings of emission (absorption) in the amplitude and the complex-conjugate amplitude. We depict the two main possibilities in Fig. 4. The process in Fig. 4a depicts an emission from the heavy-quark (in the amplitude) and later absorption by a gluon (in the complex conjugate), while the second diagram, Fig. 4b, describes an emission from gluon and subsequent absorption by the heavy-quark. The remaining contributions can be found by adding the complex conjugate of these terms to the final answer, and will be automatically included in the expressions below.
We now have for the diagram Fig. 4a where gluon is emitted at time t and absorbed at a later timet, i.e.t > t, by a heavy quark. Then, the double-differential spectrum reads where dΩ k = dωd 2 k/[(2π) 3 2ω] is the Lorentz-invariant phase-space element, TrΣ ≡ if abc tr(t a t b t c ) = −N c (N 2 c − 1)/4 and 2Re accounts for the two possible time orderings of this process, and we have used the short-hand notation u ≡ d 2 u. The normalization factor (N c C F ) −1 takes care of the averaging over the colors of the initial quark and divides out the color structure related to the dipole splitting. The in-medium two-point correlators appearing in (A4) are known from the literature and, after some simplifications, we rewrite (A4) as and where we have assumed that the medium is described by a static density n. The overall color factor in (A5), −N c /2, corresponds to the correct (negative) interference charge.
After performing the derivatives and shifting the coordinates, we finally obtain The calculation of the second diagram in Fig. 4b is completely analogous. From symmetry, it turns out that dI/dω| b = dI/dω| a (Θ 1 ↔ Θ 2 ). Hence, summing up all four contributions, our final result for the emission spectrum therefore reads, We have checked that this formula reproduces the double-differential interference spectrum in vacuum. It was first discussed in Ref. [34]. Note, however, that our formula differ from the related Eqs. (4.1) and (4.2) in Ref. [34] by the second term in the square brackets.
The in-medium energy spectrum is found by integrating out the transverse momentum, in which case we obtain which is the main formula for analyzing the direct and interference rates in this work.
Appendix B: Direct and interference radiation rates in the multiple-soft scattering approximation The formulas in Eqs. (A9) and (A10) can easily be generalized to any splitting process involving a gluon emission (for the time being let us disregard photon splitting into quarkantiquark). The general formula for the energy spectrum, cf. Eq. (A10), reads, where Q ij ≡ Q i · Q j = (Q 2 0 − Q 2 1 − Q 2 2 )/2 and Q i corresponds to the color charge-vector of the emitter (e.g. for a quark Q 2 q = C F and for a gluon Q 2 g = N c ). The particle i is associated to the dead-cone angle Θ i and direction n i = p i /E i (and similarly for particle j). The two terms describe two possible emission processes, namely one that is initiated by particle i and one initiated by particle j. The case under study, the Q → Q+g splitting that contributes to the fragmentation of a heavy-quark jet, corresponds to Θ j = 0 and Q ij = −N c /2, as derived explicitly in (A9).
In (B1), S 2 (t) is a two-point function describing the decoherence of the antenna before the splitting occurs and is often referred to as the decoherence parameter [31,32,44,45].
It depends explicitly on the opening angle of the pair, in particular S 2 (t) = 1 for n 2 12 = 0. The formula in (B1) generalizes the result previously obtained in [34]. It is worth pointing out that, due to the Galilean symmetry of the problem, the spectrum does not depend on the directions of the emitters apart from the dipole opening angle n 12 .
Working in the multiple-soft scattering limit of medium interactions, the n-point functions can be found exactly. It corresponds to the following approximation on the interaction term, N c nσ(r) ≈qr 2 /2, whereq is the jet quenching parameter. For the decoherence parameter, we immediately obtain S 2 (t) = exp[−qn 2 12 t 3 /12], and for the three-point function describing interactions during the formation of the soft gluon, we find where τ ≡t − t and Ω ≡ (1 + i) q/ω/2. The propagator in (B2) constraints the extent of the time difference between gluon emission and absorption τ q/ω ≡ t br . Hence, in the limit t br L which holds for soft gluon emissions ω qL 2 that are most important for energy loss at not too high energies [27], we can approximate the time integration over the time-difference of the emissions in the amplitude and in the complex-conjugate as see, e.g., [46]. Formally, this allows to treat multiple radiation as independent with a constant rate.
Focussing again on a q → q + g splitting, we therefore obtain the formulas for the direct and interference rates. In order to separate the dynamics of the dipole before and during the emission-time of medium-induced gluon, we will explicitly define the emission rate Γ ij (ω, t) as where the decoherence parameter is responsible for the long-distance color decoherence processes that are happening prior to the emission. Setting the dipole opening angle to zero, n 2 12 → 0, we find the independent heavy-quark and gluon spectra.
where we have restored the correct mass and color factors, and introduced the definition of the dead-cone angle Θ 0 ≡ θ q used throughout. The interference spectrum contains two terms, given bỹ Furthermore, we can altogether neglect the phases involving the terms Θ 2 0 t since, in the double-logarithmic approximation the hard gluon splitting takes place at angles much larger than the dead-cone angle.
Note that this expression, at finite n 12 = 0, diverges both in the limit of small t and τ .
However, this contribution is completely independent of medium parameters, ∼ ω 2 /[πtτ ], and can therefore be regularized by subtracting the vacuum contribution. This is in contrast to the regularization of the direct terms, which only exhibit a divergence at small τ , ∼ ω 2 /[πτ ], but is cured in the same way. It follows that all medium-induced contributions vanish in theq → 0 limit. Hence the rates that are used to compute energy loss processes are therefore regularized as Γ ij (ω, t) ≡Γ ij (ω, t) −Γ ij (ω, t) q→0 . (B9) After applying this regularization, we immediately find the rate of direct emissions by a massless gluon, recovering the well-known LPM rate valid for soft-gluon emission. For the corresponding rate off a massive quark, we find where ψ 0 (x) is the digamma function.
Although the formulas derived above can be evaluated numerically, we aim at understanding the problem through the fundamental scales appearing in these expressions. As a final step, let us therefore spend some time on discussing the different time-scales appearing, starting with the interferences. Neglecting the dead-cone angle compared to the dipole opening angle, we find that one of the exponentials in Eqs. (B12) and (B16) start oscillating at times This scale is related to quantum coherence [32,45], imposing that the wavelength of emitted quanta resolve the dipole x ⊥ (t) > λ ⊥ (t), where λ 2 ⊥ ∼ ω/t. In the vacuum t ∼ t f , which leads immediately to the angular ordering condition. The second time-scale appearing in the exponentials is t 2 = (qn 4 12 ω) −1/4 .
This scale is not an independent scale since we can rewrite it as t 2 = (t 3 d t 1 ) 1/4 , where t d is the color decoherence scale introduced above. This only allows for two possible orderings, t d < t 2 < t 1 and t 1 < t 2 < t d , and in neither of the cases does t 2 constitute the shortest, and therefore most relevant, time-scale.