Subleading power factorization with radiative functions

The study of amplitudes and cross sections in the soft and collinear limits allows for an understanding of their all orders behavior, and the identification of universal structures. At leading power soft emissions are eikonal, and described by Wilson lines. Beyond leading power the eikonal approximation breaks down, soft fermions must be added, and soft radiation resolves the nature of the energetic partons from which they were emitted. For both subleading power soft gluon and quark emissions, we use the soft collinear effective theory (SCET) to derive an all orders gauge invariant bare factorization, at both amplitude and cross section level. This yields universal multilocal matrix elements, which we refer to as radiative functions. These appear from subleading power Lagrangians inserted along the lightcone which dress the leading power Wilson lines. The use of SCET enables us to determine the complete set of radiative functions that appear to O\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$ \mathcal{O} $$\end{document}(λ2) in the power expansion, to all orders in αs. For the particular case of event shape observables in e+e−→ dijets we derive how the radiative functions contribute to the factorized cross section to O\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$ \mathcal{O} $$\end{document}(λ2).


Introduction
The simplicity of the soft and collinear limits of gauge theories allows for an all orders understanding of the behavior of amplitudes and cross sections, typically formulated in terms of factorization theorems [1][2][3]. Unlike for observables which are amenable to a local operator product expansion (OPE) [4], these general factorization theorems typically involve non-local matrix elements with Wilson lines. While the structure of these matrix elements is well understood at leading power, the structure of power corrections is much less well understood. In general, complicated non-local matrix elements, typically involving operators strung along the light cone dressing the leading power Wilson line structure, are required [5][6][7][8][9]. The soft collinear effective theory (SCET) [10][11][12][13], an effective field theory describing the soft and collinear limits of QCD, provides an operator and Lagrangian based formalism for deriving factorization theorems at subleading power. As an example, SCET has been used to systematically study power corrections to the leading power factorization for B → X s γ, B → X u lν [13,14] in the shape function region [15][16][17], and derive subleading factorization theorems in terms of universal non-local operators [8,9,[18][19][20][21][22][23]. In this case, the power corrections take the form of non-local operators describing both soft fluctuations at the scale Λ QCD , in terms of matrix elements of the B meson, as well as the coupling of soft and collinear modes.
Recently, there has been significant interest in understanding the subleading power soft and collinear limits of perturbative scattering amplitudes and event shape observables. This has been motivated at the amplitude level both by their relation to asymptotic symmetries (see e.g. [24][25][26][27][28][29][30][31]), as well as to better understand the structure of amplitudes by studying their limits (see e.g. [32][33][34][35][36]). At the cross section level an understanding of subleading power corrections will allow for the improved accuracy of perturbative predictions involving resummation, and improvements to next-to-next-to-leading order subtraction schemes [37][38][39] by analytically calculating subleading power corrections [40][41][42][43][44][45][46], amongst many other applications. From explicit calculations, there are hints for the simplicity of power corrections at higher loop order, for example in splitting functions [47], in the threshold limit [48][49][50][51][52][53], for event shape observables [40-42, 54, 55], for power corrections in quark masses [56,57], and in the Regge limit [58]. To obtain an all loop understanding, and identify universal structures which persist at subleading powers, it is desirable to formulate subleading power factorization theorems whose renormalization group structure allows the prediction of higher loop results from lower loop data, as has been successful at leading power. Recently this was used to derive the first resummation at subleading power for the thrust event shape observable in H → gg [59] and for threshold in [60].
In this paper we use SCET to derive an all orders gauge invariant factorization for subleading power soft emissions, focusing in particular on non-local corrections described by so called radiative functions. We use the SCET Lagrangian, formulated in terms of non-local gauge invariant quark and gluon fields to provide gauge invariant definitions of the radiative functions for the emission of both soft quarks and gluons. Gauge invariance is guaranteed by an intricate Wilson line structure, dictated by the symmetries of the effective theory. We show how these radiative functions appear in factorization formulas at subleading power, both at the level of the amplitude and the cross section, as multilocal matrix elements with convolutions of operators along the lightcone. These operator insertions correct the leading power Wilson line structure. This completes our derivation of all the required components for subleading power factorization initiated in [61,62], and we review in detail the complete factorization structure at subleading power, highlighting the role that radiative functions play.
If we consider the subleading power emission of a soft quark or gluon from a hard scattering vertex, there are two potential classes of contributions, as shown in figure 1. First, there are contributions where the soft emission localizes to the hard scattering vertex, as shown in figure 1. Here the power suppression arises due the lack of a nearly onshell propagator, and these contributions are described by local hard scattering operators, complete bases of which are known forqΓq [61,63] and gg currents [62], as well as recently for N -jet configurations [64,65]. Second, there are contributions from a non-local emission from the energetic parton, which arise from corrections beyond the eikonal limit to the dynamics of the interaction of the soft and collinear particles. Such contributions were studied in the abelian case in the work of Del Duca [66], extending the work of Low, Burnett and Kroll (LBK) [67,68], and were referred to as radiative jet functions. For the emission of a single soft gluon from an energetic quark, they have been extended to the non-abelian case in [69,70]. They were also studied in [30] using SCET, where a one-loop expression for soft emission was derived. Our work goes beyond this, by providing explicit all orders factorization in terms of gauge invariant soft and collinear matrix elements.
Here we will refer to the general class of such objects as radiative functions. We reserve "radiative jet function" for the analogous objects at cross section level.
To provide an all orders description, one must consider a subleading power soft emission in the presence of an arbitrary number of leading power soft, or collinear emissions. At leading power, the energetic partons emitted from the hard scattering eikonalize, and act as a source for the long wavelength soft radiation. In this limit, the dynamics of the energetic partons can be integrated out, and replaced with a Wilson line along their path. This is shown schematically in figure 2. This leads to the ubiquitous appearance of cusped light-like Wilson lines in the description of the soft and collinear limits of gauge theory amplitudes and cross sections, whose renormalization is controlled by the universal Γ cusp [71,72]. Beyond leading power we expect corrections to this picture associated with the breakdown of eikonalization, namely we expect the Wilson lines to be decorated with operators, which we will associate with radiative functions. To achieve subleading power factorization of amplitudes and cross sections, and to understand the universality of these factorizations, we would like to have a systematic approach to the construction of gauge invariant radiative functions in terms of well defined field theoretic objects. This is more difficult due to the nature of the operators, which possess intricate Wilson line structure to ensure gauge invariance in a non-abelian theory. In QED, this is not an issue since the gauge group is abelian.
To see how radiative functions naturally emerge from the effective theory, we consider the SCET Lagrangian (here we restrict ourselves to the case of SCET I ), which consists of both hard scattering operators, and a dynamical Lagrangian each of which is a power expansion in λ. The hard scattering operators, included in L hard , describe all the localized contributions in figure 1, while the non-local contributions are described by the dynamical Lagrangians, L (i) . The dynamical Lagrangian is universal, and known up to O(λ 2 ) [18,[73][74][75][76][77]. Finally, L G is the leading power Glauber Lagrangian [78]. The SCET Lagrangian is fixed by the symmetries of the theory, namely soft and collinear gauge symmetries and reparametrization invariance [74,75], and is known to not be renormalized to all orders in α s [18], which will allow us to prove the universality of the radiative functions.

JHEP11(2019)153
In the SCET framework, the leading power eikonalization of figure 2 is achieved at the level of the Lagrangian through the BPS field redefinition [79], B aµ n⊥ → Y ab n B bµ n⊥ , χ α n → Y αβ n χ β n , (1.2) which is performed in each collinear sector. Here Y n , Y n are fundamental and adjoint ultrasoft Wilson lines, where for a general representation, r, the ultrasoft Wilson line is defined by where P denotes path ordering. The BPS field redefinition decouples the ultrasoft degrees of freedom from the leading power collinear Lagrangian [79], so that they appear only in the hard scattering vertex. Beyond leading power, the BPS field redefinition does not decouple the ultrasoft and collinear interactions. However, since subleading power interactions can only occur a finite number of times at a given power, they can be viewed as decorating the Wilson line, and giving rise to radiative functions. Consider a subleading power emission, in the presence of an arbitrary number of additional soft emissions, as shown in figure 2. The insertion of the subleading power Lagrangian implies that the soft emissions cannot simply be pulled back into a Wilson line at the hard scattering vertex, since they become trapped at the subleading power Lagrangian insertion. The Wilson lines appearing in the BPS field redefinition can be used to sandwich the covariant derivative describing the soft gluon emission. For a derivative in an arbitrary representation, r, we have Y (r) † n i iD (r) µ us Y (r) n i = i∂ µ us + [Y (r) † n i iD (r) µ us Y (r) n i ] = i∂ µ us + T a (r) gB aµ us(i) , (1.4) which allows us to define the ultrasoft gauge invariant gluon building block field by gB aµ us(i) = 1 in i · ∂ us n iν iG bνµ us Y ba n i . (1.5) Furthermore, there will be a single remaining Wilson line at the hard scattering vertex. This is shown schematically in figure 2. After applying the BPS field redefinition, there are no further interactions between soft and collinear partons. The finite number of subleading power interactions between soft and collinear fields at a given power are represented by gauge invariant operator insertions along the lightcone, which dress the leading power Wilson lines. These give rise to universal non-local string operators appearing in subleading power factorization theorems. A similar picture also applies to soft quark emission. This provides a systematic way to provide gauge invariant operator definitions of radiative functions in the effective theory, which is the goal of this paper. The use of non-local gauge invariant fields is crucial to achieve factorization for these non-local operators, since it enables gauge invariant definitions of soft and collinear matrix elements tied together by convolution variables, that can be separately renormalized. This is non-trivial in a non-abelian gauge theory, where the soft emission carries a color charge.

JHEP11(2019)153
An important result of this paper is that we will show how to achieve an all orders gauge invariant factorization at the cross section level for the radiative contributions to event shape observables. This will allow us to express the cross section as a convolution of gauge invariant collinear and soft factors, each of which can in principle be separately renormalized. Unlike at leading power, these subleading power factorizations will involve an additional convolution over the gauge invariant momentum (or equivalently position along the light cone) of the insertion of the gauge invariant field.
As an example, we can consider the factorization involving a radiative contribution at cross section level involving a soft quark field. We will show that such a contribution can be writen as a convolution which can be shown schematically as This involves a factorization into a convolution over a collinear radiative function which emits the soft quark, and a soft function involving a gauge invariant soft quark field, in addition to the standard Wilson lines. This form of factorization allows us to describe systematically, and in a gauge invariant manner, all subleading power corrections to event shape observables. As noted earlier, non-local gauge invariant operators describing the coupling of soft and collinear particles have appeared in the SCET literature on subleading power corrections in B-physics [8,9,[13][14][15][16][17][18][19][20][21][22][23], many of which have a similar structure to those discussed in this paper. Here we provide a simplified approach by determining the subleading power Lagrangians in terms of gauge invariant building blocks, and then we focus on the application to perturbative scattering amplitudes and collider event shape observables. We will also emphasize the connection to the radiative functions approach which has been studied in the literature [69,70].
An outline of this paper is as follows. In section 2 we review the key features of SCET that will be important for our analysis, in particular, the definition of gauge invariant collinear and soft fields and the BPS field redefinition. In section 3 we discuss the different contributions to factorization at subleading power, and how they appear in SCET, placing in context the contributions from radiative functions. In section 4 we derive the form of the subleading power SCET Lagrangians describing the interactions of the non-local gauge invariant quark and gluon fields, using the BPS field redefinition and the equations of motion. In section 5 we use these Lagrangians to study subleading power factorization at the amplitude level, showing how radiative functions naturally emerge, and deriving in detail their structure. We also compare our radiative functions to those previously discussed in the literature. In section 6 we consider radiative functions at the level of the cross section for event shape observables, and derive the structure of convolutions between the radiative functions. In section 7 we classify all radiative functions contributing to SCET I type observables in e + e − → dijets. We conclude in section 8.

Review of SCET
In this section we briefly review the aspects of SCET that will be needed in this paper, in particular the use of gauge invariant quark and gluon fields. Reviews of SCET can be found in refs. [80,81], and [80] in particular contains a discussion of the formalism required for SCET at subleading power.

General formalism
SCET is an effective field theory of QCD describing the interactions of collinear and soft particles in the presence of a hard interaction [10][11][12][13]79]. Since the dynamics of collinear particles is concentrated along the light cone, it is convenient to use light cone coordinates. For each relevant light-like direction (jet direction) we define two reference vectors n µ i and n µ i such that n 2 i =n 2 i = 0 and n i ·n i = 2. Any momentum p can then be written as A particle is referred to as "n i -collinear" if it has momentum p close to the n i direction. More precisely, the components of its momentum must scale as (n i ·p,n i ·p, p n i ⊥ ) ∼n i ·p (λ 2 , 1, λ). Here λ 1 is a formal power counting parameter defining the expansion. Different reference vectors n i and n i , with n i · n i ∼ O(λ 2 ) provide equivalent descriptions. This enforces a symmetry on the effective theory known as reparametrization invariance (RPI) [74,75]. We will often use this symmetry to simplify our description, for example by choosing that the total ⊥ momentum of a particular collinear sector vanishes. Extensions of SCET to describe more complicated situations with collinear particles with soft momentum n i · p 1 are referred to as SCET + [82][83][84][85]. SCET is constructed as an expansion about the light cone in powers of λ. Explicitly, momenta are expanded into label and residual components where,n i ·p ∼ Q andp n i ⊥ ∼ λQ are the large label momentum components, with Q a characteristic scale of the hard interaction, while k ∼ λ 2 Q is a small residual momentum. A multipole expansion is then performed to obtain fields with momenta of definite scaling, namely collinear quark and gluon fields for each collinear direction, as well as soft quark and gluon fields. Independent gauge symmetries are enforced for each set of collinear or soft fields. As a consequence of the multipole expansion all fields and their derivatives acquire a definite power counting [12], shown in table 1. The detailed structure of the fields will be described shortly. Similarly, the SCET Lagrangian is expanded as 3)

JHEP11(2019)153
with each term having a definite power counting O(λ i ). As written, the SCET Lagrangian is divided into three different contributions. The L (i) hard contains local hard scattering operators, and is derived by a matching calculation. The L (i) describe the long wavelength dynamics of ultrasoft and collinear modes in the effective theory, and can be found in [80] to O(λ 2 ). The Glauber Lagrangian [78], L (0) G , describes interactions between soft and collinear modes in the form of potentials, which break factorization unless they can be shown to cancel. In this paper, we will considerably simplify the structure of the subleading power SCET Lagrangians not involving Glauber exchange, re-organizing them using the equations of motion of the effective theory, and writing them in terms of gauge invariant operators. Subleading power contributions to the Glauber Lagrangian that describe quark reggeization are also known [86].
We will write the SCET fields for n i -collinear quarks and gluons, ξ n i ,p (x) and A n i ,p (x), in position space with respect to the residual momentum and in momentum space with respect to the large momentum components. The large momentump, and the collinear direction then act as labels for the fields. Derivatives acting on the fields give the residual momentum dependence, i∂ µ ∼ k ∼ λ 2 Q, while label momentum operators P µ n i give the label momentum P µ n i ξ n i ,p =p µ ξ n i ,p . An important feature of the multipole expansion is that the propagator for collinear fields is independent ofn · p r and p µ r⊥ and hence is local in the residual x − , x ⊥ components. This will allow for factorized expressions involving convolutions to be reduced to a single variable convolution in the position along the lightcone, and will play an important role in our definitions of the radiative functions.
Soft degrees of freedom are described in the effective theory by separate quark and gluon fields, q us (x) and A us (x). We will assume that we are working in the SCET I theory where these soft degrees of freedom are referred to as ultrasoft. These fields do not carry label momentum, and have i∂ µ ∼ λ 2 Q.

Gauge invariant fields and the BPS field redefinition
In our study of radiative functions the use of gauge invariant soft and collinear fields will play a central role, allowing us to systematically construct gauge invariant non-local operators. Collinearly gauge invariant quark and gluon fields are defined as [11,12] is the collinear covariant derivative and is localized with respect to the residual position x, and we can therefore treat χ n i ,ω (x) and B µ n i ,ω (x) as local quark and gluon fields from the perspective of ultrasoft derivatives ∂ µ that act on x.
To formulate gauge invariant radiative functions, which describe the subleading power emission of a ultrasoft quark or gluon, it will be essential to define gauge invariant ultrasoft quark and gluon fields. This can be done using the BPS field redefinition of eq. (1.2), which we repeat here for convenience Gauge invariant ultrasoft operators will be necessarily non-local at the ultrasoft scale, involving ultrasoft Wilson lines. However, the form of this non-locality is completely determined by the BPS field redefinition. Matching calculations from QCD to the effective theory are performed prior to the BPS field redefinition, when the theory is local at the hard scale, and then the non-locality arises only from the BPS field redefinition. However, one can take a bottom-up approach and consider gauge invariant soft gluon fields as building blocks of the theory. This approach has been used in [61][62][63] to construct the operator bases at subleading power. Here we will show that the subleading power Lagrangians of SCET can be rewritten after BPS field redefinition purely in terms of gauge invariant ultrasoft quark and gluon building blocks and collinear fields. This will play an important role, allowing for the definition of gauge invariant radiative functions.
To define gauge invariant ultrasoft fields, we can group all Wilson lines arising from the BPS field redefinition with fields to form gauge invariant combinations. In particular, we can define an ultrasoft gauge invariant quark field as Similarly, we can group Wilson lines with gauge covariant derivatives in an arbitrary representation, r, to define gauge invariant ultrasoft derivatives, and gauge invariant ultrasoft gluon fields Here the square brackets indicate that the covariant derivative acts only on the Wilson lines within the brackets. In both eqs. (2.9) and (2.10), the fields are ultrasoft gauge invariant for an arbitrary lightlike direction, n i , and the soft fields themselves are not naturally associated with a given direction. Due to eq. (2.8), this direction is usually naturally taken to coincide with that of a collinear direction, for example the direction of the collinear

JHEP11(2019)153
particles emitting the soft quark or gluon. Note that the ultrasoft gauge invariant gluon field is the analogue of the gauge invariant collinear gluon field of eq. (2.5), which can also be written The gauge invariance of the ultrasoft fields in eqs. (2.9) and (2.10), is enforced by the presence of the soft Wilson line. This also implies that it has Feynman rules describing an arbitrary number of soft emissions. For example, expanded up to two emissions, we have As was described in the section 1, and as illustrated in figure 2, this particular combination sandwiching an ultrasoft emission will naturally appear when considering subleading power emission from a collinear sector, giving rise to the emission of a gauge invariant gluon field B us(i) . The ability to formulate the emission of a soft parton in a gauge invariant manner in a non-abelian theory relies crucially on the use of non-local gauge invariant fields of eq. (2.11), as it will allow the definition of separately gauge invariant soft and collinear matrix elements tied together by a convolution variable which represents the momentum flowing into the soft emission.

Components of subleading power factorization
In this section we discuss the different components contributing to the formulation of a subleading power factorization theorem in SCET, extending the brief discussion provided in section 1. Although some of the different components of the factorization, such as the construction of operator bases [61][62][63] and the renormalization of subleading power operators [59,60,64,65,87,88], have been studied extensively in the literature, we wish to provide here a self contained discussion, and we hope that it puts into context the role of radiative functions in subleading power factorization more generally. We will take as a concrete example an SCET I event shape in e + e − → dijets, however, the considerations are more general. Throughout this section, we will not consider possible contributions from leading power Glauber modes, which we assume decouple from the soft and collinear modes. For the case of e + e − this is reasonable, since all QCD particles are in the final state, where we expect leading Glauber effects can be absorbed into the direction of soft Wilson lines. Under this assumption, subleading power corrections to the Glauber Lagrangian of [78] could be constructed and analyzed in the same manner discussed here, although this is beyond the scope of this paper. 1 Consider a dimensionless SCET I event shape observable, τ , in e + e − → dijets, which is chosen to vanish in the dijet limit. As a concrete example, one can consider τ = 1 − T ,

JHEP11(2019)153
where T is the thrust observable [89]. We can write the cross section as a power expansion where dσ (n) /dτ ∼ τ n/2−1 due to the scaling relation λ ∼ √ τ . We wish to find factorized expressions for the non-zero contributions in this series, in terms of hard, jet and soft functions of the schematic form Here H are field theoretic matrix elements involving only collinear or soft fields, and Q is the center of mass energy. Here the superscipts denote the power suppression, and we have n = n Hj + n Jj + n Sj . (3.3) In general, there will be multiple distinct hard, jet and soft functions at each power.
In the effective theory approach, the first step towards this goal is to write the expression for the differential cross section in terms of full theory QCD matrix elements, where for e + e − → dijets through a virtual photon, O = J µ L µ , where L µ is the leptonic current which includes the photon propagator and couplings, and J µ =qγ µ q. Here we use the shorthand notationδ (4) q = (2π) 4 δ 4 (q−p X ) for the momentum conserving delta function. The summation over all final states, X, includes phase space integrations. Here |L denotes the e + e − leptonic initial state. The measurement of the observable is enforced by δ τ − τ (X) , where τ (X), returns the value of the observable τ as measured on the final state X.
For τ 1, we are in the dijet limit and can match onto SCET hard scattering operators with two collinear sectors Here the sum is over powers in λ indicated by the superscript (k), and at each power distinct operators are labeled by λ l and j which include helicity and color labels. Our labels are split such that λ l = ± indicates the helicity of the lepton current and the index j denotes all helicity and color labels of the QCD component of the current. The C (k) λ l j coefficients include the electromagnetic coupling and charges. We work to all orders in the strong coupling, α s , but to leading order in the electroweak couplings. We can therefore factorize out the leptonic component, J λl of the hard scattering operators (3.6)

JHEP11(2019)153
Evaluating the tree level matrix element involving the external electron states, the expression for the cross section can be written in terms of matrix elements in the effective theory as (3.7) After having calculated the leptonic matrix element we are left with a normalization factor N , whose explicit form is not relevant for the current discussion.
To achieve an expression with homogeneous power counting, as in eq. (3.1), we must systematically expand eq. (3.7) in λ, working to all orders in α s . At leading power, assuming that the action of the measurement function factorizes, this is simple. The BPS field redefinition decouples leading power soft and collinear interactions so that the Hilbert space factorizes, and the state |X can be written |X = |X n |Xn |X us . (3.8) Algebraic manipulations can then be used to organize eq. (3.7) into a form involving separate matrix elements of soft and collinear fields, and hence derive a bare factorization formula.
In the effective field theory organization it is then evident from eq. (3.7) that there are three sources of power corrections 1. Subleading power hard scattering operators.
2. Subleading power corrections to the measurement function.

Subleading power Lagrangian insertions.
We will briefly discuss the structure of each of these sources of power corrections in turn, but our primary focus will be on the factorization of subleading power Lagrangian insertions, as these give rise to the subleading power radiative functions which are the focus of this paper.

Hard scattering operators
The matching of QCD onto SCET gives rise to hard scattering operators, see eq. (3.5). These operators are local at the scale of the matching, as shown schematically in figure 1. Subleading power hard scattering operators with two collinear directions were recently discussed in detail in refs. [61][62][63] where complete bases were derived forqΓq and gg currents using the approach of helicity operators [90,91]. The leading order matching was also performed. In addition, operator bases for N -jet configurations were studied in [64].
Hard scattering operators at subleading power are similar to those at leading power in that they are formed from products of the SCET operator building blocks of table 1. These building blocks provide a complete basis of building blocks to all powers, as can be proven by the use of equations of motion and operator relations [92]. The difference between leading and subleading power comes from additional collinear or ultrasoft fields, or P ⊥ operators which are inserted into the hard scattering operator to give the power suppression. For JHEP11(2019)153 example, leading power hard scattering operators for more inclusive processes typically have a single collinear field in each collinear sector, but at subleading power can have multiple collinear fields in a single sector.
If the entire power suppression comes from the hard scattering operator then the factorization proceeds similar to at leading power. 2 In particular, in this case only the leading power Lagrangian is required since subleading power Lagrangian insertions would induce additional power suppression, and therefore factorization formulae can be derived through the BPS field redefinition. We include the case that the power suppression comes from a mixture of suppression in the hard scattering operator and suppression from the subleading power Lagrangian in section 3.3, where the treatment is more complicated due to the subleading power Lagrangian.

Measurement function factorization
The action of the measurement function τ (X), which is a function of the soft and collinear momenta, must also be expanded homogeneously in the soft and collinear limits. As shown in eq. (3.7), we expand the measurement function as Here we have assumed that any O(λ) = O( √ τ ) corrections to the measurement function vanish. 3 This was shown explicitly for the case of thrust in [61].
The measurement function enters eq. (3.7) as a delta function constraint on the final state. This constraint can be expanded as where the dots represent higher derivatives of delta functions. To achieve a factorization, one must show that the measurement operators at each power can be factorized into contributions from soft or collinear degrees of freedom. To be specific, we restrict ourselves to what we have referred to as "pseudo-additive observables" [61] which we defined as those observables with measurement functions that can be factorized into contributions from collinear and ultrasoft modes at each order in the power expansion in the form The factors G n , Gn, G us , which can enter the measurement function for any sector, are global properties of a sector, and must be defined independent of the order in perturbation 2 Although the final formulas typically have a richer convolution structure since subleading power hard scattering operators can have multiple fields per collinear sector. 3 This is true for most observables in the SCET formulation of refs. [10][11][12][13] in which label and residual momentum are exactly conserved. In the approach of ref. [93], where momentum is not strictly conserved, an O(λ) contribution to the measurement function does appear. As shown in ref. [94] this O(λ) contribution to the measurement function contributes as a product with an O(λ) operator arising from the expansion of momentum conserving delta functions, and contributes at λ 2 .

JHEP11(2019)153
theory. 4 In this case one can define field theoretic measurement functions M us from the energy momentum tensor of the theory [100][101][102][103][104]. The measurement functions act as The subleading power measurement function has been derived for the thrust observable in [61,94] to O(λ 2 ). Note that subleading power corrections can also arise for measurement functions of observables, such as kinematic factors, that are not small in the λ 1 limit. This is particularly important when multiple measurements are performed on the final states as occurs for Born measurements in fully differential cross sections at hadron colliders.
As an example, let us consider the beam thrust [105] event shape T 0 or the q T spectrum in color singlet production at the LHC. To obtain distributions that are fully differential in the momentum of the color singlet, one needs to include not only a measurement function for the observable T 0 or q T , but also a measurementδ Y ≡ δ(Y −Ŷ (X)) for the rapidity and one for the invariant mass of the color singletδ Q ≡ δ(Q 2 −Q 2 (X)). To be precise, if we call q µ the 4-momentum of the color singlet in the hadronic center of mass frame, the rapidity and the invariant mass measurements take the form At Born level the T 0 or q T measurement gives δ(T ) or δ(q T ), however the observables defined in eq. (3.13) are in general non trivial already at the Born level. Hence, they are referred to as Born measurements. As shown in the fixed order calculations of [40,42,44,45] and explained in detail in [44], the power corrections to the Born measurements contribute significantly to the power correction of the differential distribution. In particular they introduce new non-perturbative functions, namely derivatives of the Parton Distribution Functions (PDFs), which do not appear at leading power. Since we are interested in deriving factorization to O(λ 2 ), and the first subleading power correction to the measurement function appears at O(λ 2 ), contributions to the cross section whose power suppression arises from the measurement functions can be factorized just like at leading power by using the BPS field redefinition. Any insertion of subleading power Lagrangians or hard scattering operators would lead to further power suppression.

Factorization with Lagrangian insertions
The most non-trivial aspect of subleading power factorization is the factorization of the subleading power Lagrangians. At leading power this is achieved in SCET through the BPS JHEP11(2019)153 field redefinition, however, the BPS field redefinition does not decouple soft and collinear interactions beyond leading power. The Lagrangian governing the dynamics of the effective theory has the power expansion (3.14) When working to any fixed power in λ only a finite number of insertions of L (i) , i ≥ 1 need to be considered. Explicitly, if we consider a time-ordered product (T -product) in the effective theory and we are interested in its expansion to O(λ 2 ), we have where the dots represent higher power corrections. In the final expression all matrix elements are evaluated using the leading power SCET Lagrangian, and the subleading power Lagrangians appear only a finite number of times. From now on we will drop the L (0) subscript. This expression highlights that to achieve factorization of the dynamics at any finite power in the power expansion, it is sufficient to show a decoupling of the leading power interactions. The insertions of the subleading power Lagrangians in the matrix elements will lead to the radiative functions in which we are interested. The leading power interactions of ultrasoft and collinear degrees of freedom can be decoupled at the Lagrangian level using the BPS field redefinition of eq. (2.8). After the BPS field redefinition, the leading power SCET Lagrangian decomposes as where the sum is over distinct collinear sectors. Since the leading power Lagrangian defines the time evolution, states in the Hilbert space can then also be factorized as Here we work in the interaction picture defined by the leading power Lagrangian, and considering perturbations in the power expansion. Note that these perturbations are in λ, unlike the interaction picture defining the perturbative expansion in α s . Here corrections in α s are kept to all orders,. This allows hard-soft-collinear factorization to be achieved to any power in the effective field theory. Deriving the explicit structure of the factorization in the case of subleading power Lagrangian insertions will be the main focus of this paper, and will give rise to radiative functions.

Factorized cross section to O(λ 2 )
Having understood the different sources of power corrections in the effective theory, we can now achieve a homogenous power expansion for the e + e − → dijet cross section, and give expressions at each order in the power expansion in terms of matrix elements of hard scattering operators, Lagrangian insertions, and measurement functions. Here we consider only the terms which arise up to O(λ 2 ). At O(λ 0 ) we have the simple expression in terms of the different helicity configurations of the leading power operator dσ dτ (3.18) Since all matrix elements are now evaluated with the leading power Lagrangian, there are no interactions between soft and collinear degrees of freedom, and the factorization into collinear and soft matrix elements is simply an algebraic exercise, leading to the well known factorization for the thrust observable [106][107][108]. We will review this factorization in more detail in section 6.
At O(λ 1 ) we have potential contributions from O(λ 1 ) hard scattering operators, as well as subleading Lagrangian insertions, dσ dτ where here and in the following T(T) denotes (anti-)time ordering. The vanishing of dσ dτ (1) from hard scattering operators was explained for thrust in [61], and in section 7 we will discuss the analogous explanation for the vanishing of Lagrangian insertion contributions at this order At O(λ 2 ) all three sources of power corrections contribute dσ dτ 20) or more explicitly dσ dτ (2) = N X,iδ
The O(λ 2 ) power correction for thrust was computed at fixed order to O(α s ) and to O(α 2 s ) using SCET in [94] and [40,42], respectively. Since the interactions between soft and collinear degrees of freedom have been decoupled, by algebraic manipulation of eq. (3.21), the contributions to the cross section at each order in the power expansion can be expressed as a sum of vacuum matrix elements, involving a measurement function insertion, and each containing only collinear n, collinearn, or ultrasoft fields. To do this, we write the constraint on the final state as a sum of the measurement operators We can then perform the sum over the |X n X n |, |Xn Xn|, and |X us X us | states to simplify all the matrix elements to vacuum matrix elements. The Lorentz, Dirac, and color structure can be simplified using Fierz relations, and the symmetry properties of the vacuum matrix elements, such that each matrix element is a scalar, and there are no index contractions between the soft and collinear functions, namely a completely factorized form.
In previous papers we have given complete bases of hard scattering operators [61][62][63], as well as the expansion of the measurement function [61]. Here we focus on the radiative type contributions, namely those term involving additional integrals over the position of Lagrangian insertions. We will formulate the factorization of the radiative contributions to the cross section as products of gauge invariant soft and collinear matrix elements involving either one or two convolutions, corresponding to the one or two Lagrangian insertions which can exist when working to O(λ 2 ). Explicitly, we will be able to derive a representation of the form where we choose to make the arguments of the soft functions dimension 1 and the arguments of the jet functions dimension 2 analogously to leading power. Here ⊗ denotes the convolution in the thrust variable, τ ,

JHEP11(2019)153
and we have used the symmetry under n ↔n to combine several equivalent contributions. The derivation of this factorized form at the cross section level is the main goal of this paper. The factorization derived in this paper will be at the bare level, namely we do not consider the renormalization of the hard, jet and soft functions. To derive a renormalized factorization formula, one must show that the hard, jet and soft functions can be separately renormalized, and that the convolutions in the r + 1 and r + 2 variables are well defined. This is in general non-trivial, and even in simple cases the renormalization of the subleading power jet and soft functions involves mixing with additional operators that do not appear in the matching [59,60,109], with evanescent operators [110][111][112] and possibly with EOM operators [113] (though the particular EOM operators will differ from those found in [113] due to differences in the construction of the subleading power Lagrangians), and the convolutions do not naively converge [114]. However, the derivation of a bare factorization is the first step towards a complete, renormalized factorization.
After studying the structure of the subleading power Lagrangians in terms of gauge invariant quark and gluon fields in section 4, in section 6 we will work out explicitly the structure of the factorization of the matrix elements for those contributions involving Lagrangian insertions, which give rise to the radiative jet functions. This will provide the ingredients needed to construct eq. (3.23) explicitly. This in turn yields all the pieces needed to explore the full factorization for subleading power thrust, which we plan to pursue in future work.

Subleading Lagrangians for gauge invariant fields
Having identified the different sources of power corrections in SCET, we now focus on the structure of the radiative functions. As has been emphasized, to achieve factorization into separately gauge invariant soft and collinear factors, it is essential that the radiative functions be formulated in terms of non-local gauge invariant fields, namely B us and ψ us . We therefore will derive the subleading power Lagrangians describing the interactions of these non-local fields to all orders in α s .
The general form of the subleading power Lagrangians is quite complicated, since they describe the complete dynamics of the soft and collinear sectors to all orders in α s . Nevertheless, due to the power counting and locality of the effective theory, there are a finite number of terms in each Lagrangian. Operationally, at a fixed order in perturbation theory, the number of terms in the Lagrangian which actually contribute is relatively small since most terms involve higher numbers of fields. Before proceeding to the full derivation of the subleading power Lagrangians, we give the structure of the Lagrangian in terms of field content, ignoring the detailed Dirac, Lorentz, and color structures. This is useful for understanding the general structure of the Lagrangians, and the order in perturbation theory at which different terms can contribute.

JHEP11(2019)153
At O(λ) the field structure of the Lagrangian is given by where we have organized the structure according to the collinear field content. The number of fields appearing in the Lagrangian is fixed by power counting and locality, and at O(λ) the Lagrangian involves up to three collinear fields. The operators that involve multiple collinear fields will not contribute at tree level to the emission of a soft parton from a single collinear parton, but are necessary to correctly reproduce the complete subleading power expression at loop level, or for multiple collinear emissions.
At O(λ 2 ) the field structure of the Lagrangian is given by where we have again organized the terms based on their collinear field content, and we see that the O(λ 2 ) Lagrangian involves up to four collinear fields. In this section we derive the exact form of the Lagrangians given in eqs. (4.1) and (4.2). We begin in section 4.1 by summarizing the notation used in this section and the BPS transformations of different covariant derivative operators, which will allow us to write the subleading power Lagrangians in terms of gauge invariant quark and gluon fields. In section 4.2 we discuss our reorganization of the Lagrangians using the equations of motion in the effective theory. Then, in sections 4.3-4.5 we present our simplified results for the BPS redefined Lagrangians in terms of gauge invariant quark and gluon fields, as well as relevant Feynman rules. These will be used to derive the structure of the radiative functions in sections 5 and 6.

Field redefinitions for subleading Lagrangians
The subleading power Lagrangians in SCET are typically written in a local form, which still involve the interactions of soft and collinear partons [75][76][77]. To derive subleading power factorization formulas involving radiative functions, we would like to rewrite them

JHEP11(2019)153
in terms of the non-local gauge invariant quark and gluon fields. This can be achieved by performing the BPS field redefinition and manipulating the Wilson lines into gauge invariant combinations, which is the goal of this section.
Before BPS field redefinition the subleading power Lagrangians are written in terms of a variety of different covariant derivatives which we summarize here for convenience. The gauge covariant derivatives that we will use are defined by and their gauge invariant versions are given by (4.4) It is also useful to summarize the transformation of the different derivative operators under the BPS field redefinition. These are all derived using the defining relations of the Wilson line, In addition, the ultrasoft Wilson lines commute with the label momentum operators Denoting the BPS transformation of an operatorÔ as BPS[Ô], we then have the following transformations for the derivative operators (4.8) Additional useful relations are given in appendix A.
Given these identities, it is now a straightforward algebraic exercise to compute the BPS field redefinitions of the Lagrangians. By applying the unitarity condition on the ultrasoft Wilson lines, all ultrasoft Wilson lines can either be cancelled, or absorbed into gauge invariant soft quark or gluon fields, as defined in eqs. (2.9) and (2.11). To illustrate explicitly how this works, we consider two simple examples. First, consider a term from the leading power collinear gluon Lagrangian,

JHEP11(2019)153
In this case, all the soft Wilson lines explicitly cancel, decoupling the interactions of the ultrasoft and collinear gluons. As a second example we consider a term from L (1) which contains an explicit D us . Here we find that the ultrasoft gluons do not decouple In the last step we used the definition of the gauge invariant ultrasoft gluon field. The derivation of the BPS field redefinition for other terms in the Lagrangian proceeds similarly, so in the following sections we will simply state the final results for the BPS redefined Lagrangians.

Simplifications using the equations of motion
In addition to writing the subleading power Lagrangians in terms of the non-local gauge invariant quark and gluon fields, we can also simplify their structure using the equations of motion. Recall that when building bases of hard scattering operators, only the gauge invariant building blocks in table 1 are required. In particular, for the collinear gluon field, only the two degrees of freedom in B n⊥ appear explicitly, and not the other components of B n . In particular, the large components of the gauge fieldn · A n appear entirely in Wilson lines, and the small components have been eliminated using the equations of motion. We begin by reviewing how this is achieved, following the results of [92], and then apply the same simplifications to the subleading power Lagrangians.
In SCET the collinear gauge invariant covariant derivative is given by which can be broken into components as where we have defined the gauge invariant fields for the different components as Here theP n operators act only within the external square brackets. We can now eliminate the n · B n component of the gluon field using the equation of motion This allows the Lagrangian to be written entirely in terms of B n⊥ fields. From the form of eq. (4.14) we can see why this will lead to significant simplifications when studying soft emissions from a single collinear gluon, since all terms on the right hand side involve either a higher number of fields, or the P ⊥ operator. When studying soft emission at lowest order

JHEP11(2019)153
and lowest multiplicity, any term of the form B us B n n · B n can therefore be dropped, which will simplify our discussion of the radiative functions. Additionally, it is also possible to eliminate from the Lagrangian all instances of the ultrasoft derivative operator n · ∂ n acting on n-collinear fields. This is achieved for the collinear quark field using the equation of motion and for the collinear gluon field usinḡ These equations of motion, particularly for the gluon case are considerably more cumbersome. When writing the full Lagrangian, as well as for performing fixed order calculations, we therefore find it simpler to work with ultrasoft derivatives. However, we note that if we are interested in tree level soft emissions off of a single collinear line, an identical discussion as for n · B n applies, and we can ignore all appearances of n · ∂ us acting on n collinear fields in the Lagrangian. By using these equations of motion, we are therefore able to greatly simplify the structure of the radiative functions we consider.

Lagrangian at O(λ 0 )
For completeness, we begin by considering the leading power SCET Lagrangian. Those familiar with the leading power BPS field redefinition and SCET Lagrangian can skip to the next section. Before BPS field redefinition, the leading power Lagrangian involves interactions between collinear and ultrasoft particles. It can be written as [10][11][12][13] where and the ultrasoft Lagrangian, L us , is simply the QCD Lagrangian. Throughout this paper, we use a general covariant gauge with gauge fixing parameter ζ for the collinear gluons, and c n are the corresponding ghosts.
After performing the BPS field redefinition we have

JHEP11(2019)153
where the ultrasoft Lagrangian is unchanged. The collinear quark Lagrangian is given by 20) and the collinear gluon Lagrangian is given by 5 explicitly showing that ultrasoft and collinear interactions have been decoupled to leading power.

Lagrangian at O(λ)
Before BPS field redefinition, the O(λ) Lagrangian can be written where [74][75][76][77] describes the interactions between collinear quarks and gluons, and describes the dynamics of the pure gluon sector, including gauge fixing terms, 6 and L (1) χnqus =χ n g / B n⊥ q us + h.c., (4.25) describes the coupling of soft and collinear quarks. We now wish to express the subleading power Lagrangians in a simplified form in terms of the gauge invariant building blocks, which will be one of the main results of this paper. This organization of the Lagrangians after BPS field redefinition was also considered in [30], although there it was performed schematically. Here we will provide explicit expressions for all components, as well as use the equations of motion to simplify the result so that it can easily be used for subleading power factorization.
After performing the BPS field redefinition, we can perform the same division of the Lagrangian as above, 2 which is the form sometimes used in the literature to write down this term of the collinear leading power Lagrangian. 6 Note that the presence of power suppressed gauge fixing Lagrangians is necessary due to the fact that RPI symmetry connects Lagrangians at different orders in the power counting, and would be broken if they were not included. For example, these subleading power gauge fixing Lagrangians have been shown to give important contributions to the derivation of the LBK theorem for gluons in SCET, see appendix D of [30].

JHEP11(2019)153
where the collinear quark Lagrangian is given by the collinear gluon Lagrangian is divided into three pieces 28) which are given by Finally, the interaction of soft quarks is described by the Lagrangian The structure of the O(λ) Lagrangian is quite complicated, since it describes the complete dynamics of the subleading power corrections to the soft and collinear dynamics, including ghost and gauge fixing terms. In its current form, it also involves multiple polarizations of the collinear gluon field. To simplify its structure, we use the equations of motion, 7 as discussed in section 4.2. Simplifying the result to focus only on ultrasoft emissions out of two collinear fields we find the structure where K (1) Bus and K (1) ∂us contain ≥ 3 collinear fields, and are therefore not relevant for our current analysis. In this form, the Lagrangian is written entirely in terms gauge invariant fields, and due to the organization in terms of fields, it is clear at which order in perturbation theory each term can contribute. After performing the BPS field redefinition, and writing the result in terms of collinear and soft gauge invariant fields, the soft and collinear fields are only coupled through Lorentz and color indices, as well as through potential derivative operators. Since each of the building blocks appearing in the Lagrangian is separately gauge invariant, this will allow for a simple factorization into collinear and soft components, tied together through Lorentz and color indices, which will give rise to the radiative functions. The first three terms of eq. (4.31) describe the O(λ) emission of a soft gluon from a collinear gluon, a soft gluon from a collinear quark, and a soft quark from a collinear quark or gluon, respectively. Using the Lagrangian, we can derive the tree level Feynman rules, which are given in figure 3. Note that in accord with the LBK theorem, the single ultrasoft gluon Feynman rule of L (1) vanishes when the label ⊥ momentum of the collinear leg is set to zero. Unlike for the emission of a soft gluon, the Feynman rule for a soft quark emission does not vanish when the P ⊥ of the collinear line vanishes.
Since the Lagrangian is defined in terms of gauge invariant soft quark and gluon fields, which involve ultrasoft Wilson lines, they also give the Feynman rules for an arbitrary number of additional leading power soft gluon emissions. Similarly, the gauge invariant collinear fields also involve collinear Wilson lines, which describe collinear radiative corrections to the above Feynman rules. The K (1)

Bus and K
(1) ∂us in eq. (4.31) involve additional collinear fields. For a single soft emission from a collinear line, these can first appear at loop level. We will not work out the explicit form of these loop contributions in this initial paper, however, we will discuss their contributions in later sections.

JHEP11(2019)153
where for convenience, we further decompose the gluon Lagrangian as gf . (4.33) The different components of the Lagrangian are given by After performing the BPS field redefinition, and writing the result in terms of ultrasoft gauge invariant fields, we find that the Lagrangians involving quark fields can be written The Lagrangians describing the pure glue sector are more complicated, involving both ghost and gauge fixing terms. We find that they can be written To make the O(λ 2 ) Lagrangian more tractable, we can use the equations of motion to write it entirely in terms of our basis of gauge invariant building blocks. This is a straightforward, but tedious algebraic exercise, and therefore we simply present the final result. Using the equations of motion to rewrite the Lagrangian in terms of our operator basis, we find where, as in eq. (4.31), the K contain ≥ 3 collinear fields, which will not be relevant for the discussion in this paper. This gives the Lagrangian at O(λ 2 ) in terms of gauge invariant soft and collinear quark and gluon fields in such a way that it is clear at which order in perturbation theory each term can contribute. We have used the EOM to write it entirely in terms of the B n⊥ field, eliminating the other polarizations. For practical applications, we can also apply the EOM of eq. (4.16), however, this significantly complicates the structure of the Lagrangian, and therefore we have not written it out explicitly. This simplified form of the O(λ 2 ) Lagrangian is one of the key results of this paper. We again emphasize that its highly non-trivial non-local structure, involving a multitude of soft and collinear Wilson lines, is fully determined by the structure of the BPS field redefinition, and the local SCET I Lagrangians, allowing it to be constructed systematically. This form, in terms of gauge invariant building blocks linked only by Lorentz and color indices, will allow for a straightforward factorization into radiative functions. In eq. (4.37), terms appear involving 0, 1 and 2 B us fields, as well as the gauge invariant ultrasoft quark field. Since the O(λ 2 ) Lagrangians contain various terms we are going to give the Feynman rules under the common assumption of vanishing label perpendicular momentum of all collinear fields to zero P ⊥ = 0. Under this assumption the Feynman rules are given in figure 4.

JHEP11(2019)153
It is important to emphasize that since B bν us(n) has Feynman rules with an infinite number of soft emissions, the terms involving one and two B us fields will both contribute to the complete two gluon Feynman rule. In the Feynman rules in figure 4 we have given only the contribution from the Lagrangian insertion involving two B us fields. For simplicity we have not given the two soft gluon Feynman rule from L (2) B which can be straightforwardly derived using the two gluon Feynman rule of the ultrasoft gauge invariant gluon field in eq. (2.13). These contributions are separately gauge invariant.
The other terms in eq. (4.37) involve additional collinear fields. For soft emissions from a collinear line, they can first contribute when there are collinear loops. We will not explicitly compute their loop level contributions, but will further discuss their structure in later sections. Note that at O(λ 2 ), we also have contributions from two insertions of the O(λ) operator. These will be discussed when we consider the complete classification of radiative function for e + e − → dijets in section 7.

Amplitude level factorization with radiative functions
In this section we derive factorization formulas in terms of radiative functions for soft emissions at amplitude level. While our goal is to study the factorization of event shapes, and the structure of radiative functions at cross section level, initiating our studies at amplitude level is useful for several reasons. First, it is useful for connecting to the study of the subleading power soft behavior of amplitudes, which in itself is an interesting subject to which the factorization theorems that we derive can be applied. Second, it allows us to connect to other work in the literature on radiative functions, which have typically been formulated at amplitude level. Finally, it also provides a slightly simpler situation to illustrate the general features of factorization involving radiative functions, which will persist at the cross section level. In particular, we will illustrate how radiative functions are defined as integrals along the lightcone of Lagrangian insertions, which dress the leading power Wilson lines, giving rise to a breakdown of eikonalization.
In this section we will only consider those radiative functions that are relevant for describing tree level soft emissions. In particular, this eliminates all contributions involving more than two collinear fields. Furthermore, we will use RPI to take each collinear sector in the amplitude to have a total P ⊥ = 0, which eliminates O(λ) contributions to radiative soft gluon emission, as is guaranteed by the LBK theorem [67,68]. See [30] for a detailed discussion in the context of SCET. This leaves us with the following cases of interest • Single B us emission at O(λ 2 ), • Double B us emission at O(λ 2 ), each of which will be studied in this section. Single ψ us emission could also be studied at O(λ 2 ), however, due to fermion number conservation in the soft sector, this can first contribute at cross section level at O(λ 4 ), and is therefore not of interest to us here. Terms with JHEP11(2019)153

Radiative Function
Factorization Tree Level Feynman Rule Equation additional collinear fields, that do not contribute at tree level, can be treated in an identical manner, and we will briefly comment on loop level contributions at the end of this section. For convenience a summary of the radiative functions is given in table 2, which shows the schematic factorization of the amplitude, the tree level Feynman rule for the radiative function, as well as the equation number where its definition can be found. The derivation of the factorizations are given in the text.
After extending the formalism of this section to cross section level factorization in section 6, a complete classification of the field content of all radiative jet functions contributing to a physical observable, namely thrust in e + e − → dijets, will be given in section 7. This includes those which first contribute at loop level. By fixing a particular physical process, we will be able to exploit the symmetries of the problem to slightly simplify the structure of the radiative contributions.

Leading power amplitude factorization
Before considering radiative functions at the amplitude level, we begin by briefly reviewing the well known leading power amplitude level factorization. This will help to establish our notation, as well as to emphasize distinctions when we consider the subleading power case.

JHEP11(2019)153
To study factorization at the amplitude level, we can proceed as in section 3, however, we study only the matrix element instead of the squared matrix element. Here O is a full theory QCD operator, and X is an N -jet state in the full theory. In the soft and collinear limits, we can proceed identically to the factorization at the cross section level, namely we match to the leading power N -jet operator in the EFT, which we assume has a single collinear field, X k i n i , in each collinear sector Here the κ i labels the parton identity of the n i collinear field that can either be a quark jet field χ n i or a gluon jet field B ⊥ n i . Throughout this section, we will assume for simplicity that there is a single such operator, since the structure of the leading power operator will not play a significant role in our discussion. Furthermore, we will suppress explicit contractions of color indices, since they are standard. The BPS field redefinition factorizes the Hilbert space, and hence the state into collinear states X n i | and an ultrasoft state X us |. With the interactions in the Lagrangian decoupled, the leading power factorization of the matrix element is then simply an algebraic exercise, and we obtain the factorized expression (5.4) Here κ i labels both the parton identity of the n i collinear field, as well as the representation of the Wilson line, as determined by the BPS field redefinition, and T denotes time ordering. This gives rise to the familiar factorization into a hard matching coefficient, coefficient functions describing the collinear radiation along each lightlike direction, and a soft amplitude, The leading power collinear function and soft amplitude are defined as resolve the color and direction of the collinear legs. To simplify our notation, in eq. (5.4) we have left implicit the contraction of all color indices, and denote it simply by the "dot" symbol. No Lorentz or Dirac indices are passed between the jet and soft functions, and therefore the soft degrees of freedom have no sensitivity to the spin of the collinear particles. Furthermore, the factorization is multiplicative, with no convolution in the soft momentum. We will see that when we consider the factorization involving Lagrangian insertions at subleading power, these features no longer hold.

Definition of radiative functions
We now consider amplitude level factorization at subleading power. Here we will focus solely on contributions from Lagrangian insertions, which will give rise to radiative functions. We have studied the structure of subleading power hard scattering operators extensively in previous papers [61,62]. After performing the BPS field redefinition, the contributions from subleading power Lagrangian insertions to the amplitude take the form Here O (0)BPS N is the leading power BPS redefined N jet operator, and L (j) n i (x) is the O(λ j ), j ≥ 1, Lagrangian for the n i sector, after BPS field redefinition. The "rad." superscript on the matrix element emphasizes that this is only the radiative contribution to the amplitude. More generally, one must consider multiple Lagrangian insertions, or Lagrangian insertions with subleading power hard scattering operators, as detailed in eq. (3.21). These will factorize in a similar manner, and will be discussed in section 6.
Unlike the leading power case of eq. (5.4), where the amplitude factorized into a product of functions, at subleading power this factorization will include integration variables linking the jet and soft functions. These integration variables will parametrize the position along the light cone direction, and describes the position of the soft emission from the collinear sector. Furthermore, at subleading power the soft function no longer couples just to the color charge and direction of the jet functions, but can instead couple via Lorentz or Dirac indices in a manner which depends on the spin of the collinear particle.

Soft quark emission
The simplest case for which to consider the factorization is the L (1) emission of a soft quark. Unlike the L (1) emission of a soft gluon, this does not vanish when P ⊥ vanishes, and it has a simpler structure than the L (2) insertions. It therefore provides a simple demonstration of the convolution structure which will appear at subleading powers. At the amplitude level, fermionic soft theorems in supersymmetric field theories and supergravity, and their relation to asymptotic symmetries have been considered [115][116][117][118]. At the cross section level soft quarks were found to give a leading logarithmic contribution to B-physics process [9] and event shape observables [40][41][42]. The cross section level factorization will be discussed in section 6.

JHEP11(2019)153
We are therefore interested in the factorization of the matrix element The subscript ψ labeling the amplitude indicates that a soft fermion is radiated. Since the Lagrangian insertion appears only in the collinear sector n j , the factorization of the other collinear sectors proceeds exactly as at leading power, giving rise to the leading power jet functions discussed in section 5.1. For concreteness, we assume that the field in the n j collinear sector is a collinear quark. To simplify the notation in intermediate steps, we will drop the explicit time ordering, and reinstate it only in the final factorized formula. We then have 10) and it remains only to factorize the n j collinear sector from the soft sector. To facilitate a comparison with definitions of radiative functions given in the literature, it will be convenient to formulate the convolution between the jet and soft functions in momentum space.
In its current form, eq. (5.11) is written as a four dimensional convolution. To regulate this expression in dimensional regularization, one must extend this to a d = 4 − 2 dimensional convolution. This implies that one cannot separately consider the soft and collinear functions after expansion in dimensional regularization, and therefore that one cannot achieve a renormalized factorization. Furthermore, as written, the factorized expression has not yet made manifest the physical picture of decorating a Wilson line via an insertion of an operator along the light cone. To simplify the convolution structure, we make use of the multipole expansion that has been implemented in the effective theory. Due to the multipole expansion the collinear matrix elements in the effective theory are local in the n · x and x ⊥ components, where the ⊥ and + are in the light cone coordinates with respect to n j . This can also be seen from the multipole expanded propagator eq. (2.4). Using this property, we can simplify the eq. (5.11) to a single variable convolution. Focusing just on the n j and soft sectors, we have Note that here we use the lightcone definition for the x variable as This then gives the momentum space definition of the radiative jet function Jᾱ ψn j (k + ) i,s , and the corresponding soft function S α N ψn j (k + ) In position space, we have 16) where in this final form, we have explicitly reinstated the time ordering. The factorization for this contribution to the full amplitude is then given by This factorization gives the physical picture of dressing the Wilson line with an operator at a position x − along the light cone. Due to the multipole expansion, the soft sector still sees the collinear sector as lying exactly on the light cone. Unlike the leading power case, the soft function and jet function both carry a fermionic index. The soft degrees of freedom are therefore aware of the identity of the collinear partons. Other radiative functions will have an analogous structure, but can involve more complicated contractions between the soft and collinear sectors and additional integrals.

JHEP11(2019)153
The tree level result for the radiative function is given by In this initial investigation, we will not consider loop corrections. Although here we have treated an outgoing quark which emits a soft quark and becomes a collinear gluon, the opposite case of a gluon field converting to a quark field can be treated in an identical fashion.

Soft gluon emission from a collinear quark
We now consider the O(λ 2 ) emission of a soft gluon insertion. We will consider first the case of the emission from a collinear quark, where we will work through the convolution structure in detail. We will then state the result for the emission from a collinear gluon, which has a similar structure. For the emission of a soft gluon from a collinear quark, the relevant Lagrangian is We therefore must consider the factorization of the matrix element A (2),rad.
As in the previous section, we will drop the explicit time ordering until the final formula. Due to the presence of the ultrasoft derivative in the Lagrangian, the factorization is slightly more complicated than for the case of a soft quark emission. As motivation for the structure that the factorization should take, we can look at the LBK theorem for soft gluon emission at O(λ 2 ). For the emission of a soft gluon off of a collinear quark, the LBK theorem can be written as [30] S (2) (5.21) or in terms of the angular momentum generator J µν , as This expression holds only for as on-shell emission, which cannot be assumed of the group momentum flowing into the B us(n) field, and furthermore it is the complete result for the amplitude, not just the contribution from the radiative functions. Nevertheless, we would like our radiative functions to have a structure which mimics this as closely as possible. In particular, we would like that the ultrasoft derivatives appearing in the subleading power Lagrangian act only on the soft sector. Furthermore, it suggests that the radiative jet

JHEP11(2019)153
function should carry two Lorentz indices, which are contracted with a B us(n) field, and an ultrasoft derivative.
To perform the factorization, we assume for concreteness that there is a collinear quark field in the n j collinear sector. Since the Lagrangian insertion is in the n j sector we can immediately factorize the other collinear sectors, giving the leading power jet functions, and we obtain A (2),rad.
To achieve a factorization of the ultrasoft derivatives, we can choose the external states to have no ⊥ residual momentum. After BPS field redefinition, there are no soft collinear interactions, other than through the single Lagrangian insertions. Therefore, momentum is only passed between the soft and collinear sectors at this vertex. The ultrasoft derivative operator can either act on the group momentum that flows out through the B us(n) field, or on a residual momentum component coming from a collinear loop. In dimensional regularization, we have [80] q l where (q r ) j with j > 0 denotes positive powers of the q − r and q ⊥ r momenta, which are the only residual momenta which appear in the subleading power Lagrangians. Any residual momentum from a collinear loop momentum picked up by the derivative is therefore set to zero. Therefore the derivative picks up just the group momentum of the B us(n) field. Ultimately, this is due to the fact that the ultrasoft momentum of the B us(n) field is the only physical ultrasoft momentum flowing in the graphs. At loop level, this statement is slightly non-trivial since the subleading power Lagrangian can couple to closed fermion loops. However, given aχ n which produces a fermion in the hard scatter in the n collinear sector, it is possible to choose the momentum routing so that the soft momentum is routed only along the direction of fermion number flow of the collinear operator. Therefore, the ultrasoft derivative acts just on the soft gluon field, and only one of the tensor structures appears.
For the particular matrix elements of interest, we can therefore perform the following simplifications Here the square brackets indicate that the derivatives act only within those brackets. It is now straightforward to factorize this contribution to the amplitude, following the steps in section 5.2.1, so we will not repeat them explicitly. After simplifying the convolution JHEP11(2019)153 structure, we find A (2),rad.
Here, we have defined the radiative function J µνA q,n j (k + ) (5.27) We note that the subscript q denotes that this describes the radiative emission from a quark, with the fact that it is the radiative emission of a single gluon being specified by the single free adjoint color index.
This radiative function exhibits an explicit coupling of i∂ ν ⊥ B µ us to something that is reminiscent of the spin orbital angular momentum, inserted at a position along the light cone. In position space, we have where we have reinstated the explicit time ordering. The tree level result for the radiative function is In section 5.3 we will compare this with results in the literature. We have focused on the factorization for the particular insertion of the subleading power quark Lagrangian onto a collinear line, which contributes at tree level, and have shown how this gives rise to a radiative function. At loop level, one can also consider contributions from the other terms in the Lagrangian. For example, one can insert the JHEP11(2019)153 gluon component of the L (2) Lagrangian into an outgoing quark leg. This will contribute, for example through the diagram . (5.30) The calculation of contributions of this form were performed in the threshold limit in [70]. The factorization for this is identical to the case of the gluon, which will be discussed in the next section, section 5.2.3. As a matter of fact, since the coupling of the Lagrangian is the same, one must just exchange the field at the origin B n⊥ → χ n . Since the goal of this section is to show in detail for several examples how radiative functions at amplitude level arise in SCET, we will not perform a complete classification of the radiative functions at loop level. In section 7, we will perform this classification for the case of the thrust observable in e + e − → dijets, using symmetries specific to the problem to simplify the number of distinct contributions.

Soft gluon emission from a collinear gluon
We can also consider the emission of a gauge invariant soft gluon field from a collinear gluon. Since the derivation of the factorization is similar, here we skip most of the steps to quickly get to the final result. However, since the gluon fields carry both color and Lorentz indices and the resulting radiative functions definition will depend on them, in this paragraph we treat these indices explicitly. Even though the O(λ 2 ) Lagrangian in the gluon sector is somewhat complicated, it strongly simplifies under the assumption of no label perpendicular momentum flowing through the collinear fields, as explained in section 4.5. Therefore, for the emission of a soft gluon from a collinear gluon with P ⊥ = 0, the relevant Lagrangian is simply Since we want to study the factorization involving a gluon lagrangian insertion, our hard scattering operator O (0) N must contain at least one gluon field. For concreteness let's take O (0) N to have a collinear gluon field in the n j direction. We therefore must consider the factorization of the matrix element where α i are color indices of the external legs of the amplitude and they belong to the representation determined by the parton identity. 8 As usual, after BPS field redefinition we have

JHEP11(2019)153
It is convenient to isolate the gluon leg on which we want to insert the subleading lagrangian by singling out the gluon field B ρβ j n j ⊥ , where ρ is a Lorentz index and β j a color one, from the hard scattering operator, as follow where (Y g n j ) β j α j is the Wilson line 9 resulting from the BPS field redefinition of the gluon field B µβ j n j . Since all the gauge invariant gluon fields appearing in this section are perpendicular, we will drop the ⊥ label on B ⊥ n j and B ⊥ us to lighten the notation. Therefore, after BPS field redefinition, the amplitude reads (5.35) Having made the color indices explicit, we can start by factorizing all the n i collinear jets with i = j, where n j is the direction of the collinear jet on which we insert the subleading Lagrangian A (2),rad.
After simplifying the convolution structure, we find A (2),rad.
where we have defined As before, the subscript g indicates that this describes emission from a gluon, while the adjoint color index indicates that it is describing the emission of a single gluon field. The tree level Feynman rule for this radiative function is Note that this is an adjoint Wilson line, so that (Y g n j ) β j α j ≡ Yn j β j α j . In our notation this information is carried by the parton label g.

Double non-eikonal soft gluon emission
In addition to the radiative functions which couple a single B us(n) field to the collinear line, the L (2) Lagrangian also contains a term that couples to a product B us(n) B us(n) describing a double emission. It is important to emphasize that the radiative functions of section 5.2.2 with a single ultrasoft insertion also have two gluon Feynman rules arising from the Wilson lines present in the gauge invariant definition of the ultrasoft gluon field. However, such additional emissions are effectively eikonal. The double soft radiative function to be studied in this sections describes a genuinely double non-eikonal emission.
Double soft theorems have attracted some attention in the literature [119][120][121][122]. We note that as compared to some treatments in the literature, the SCET power counting is such that we always take a limit such that the soft gluons are becoming simultaneously soft. This is in contrast to consecutive limits, where the soft limit for two gluons is taken consecutively. See [123] for a discussion of the difference between these limits. As emphasized in [124], double soft limits are interesting as they allow one to probe the group structure of the theory, and are proportional to the structure constants of the theory. We will see that this is indeed true for the double radiative function.
Emissions from collinear quarks. We start with the emission from a fermionic leg. Therefore, we are interested in the factorization of the term in the L (2) χn Lagrangian involving two B us(n) fields L (2) χn,BusBus =χ n T a γ µ Since this matrix element contains no derivatives, the derivation of the convolution structure is identical to that for the soft quark emission in section 5.2.1. We therefore simply give the final result, skipping intermediate steps. We find A (2),rad.
We have defined the radiative function as

JHEP11(2019)153
The subscript q and the two free adjoint indices are meant to indicate that this radiative function describes the double emission of soft gluons from a quark field. In position space, we have A (2),rad.
The form of the J BB current is quite interesting. In particular, it involves both the symmetric and anti-symmetric structure constants, each coupling to different Lorentz structures. The antisymmetric color structure constants are associated with the perp components of the orbital momentum generator. Evaluating the radiative function at tree level, we find (5.44) It is important to emphasize that this radiative function does not completely describe double soft gluon emission at O(λ 2 ). One must also consider the two gluon Feynman rule of the radiative function of section 5.2.2 involving a single B us(n) . Interestingly, since the L (1) are directly proportional to P ⊥ , at tree level, there is no contribution to the double soft limit from a product L (1) · L (1) , or from the combination of a O(λ) hard scattering operator and a L (1) insertion.
Emissions from collinear gluons. The final tree level contribution is the emission of two gauge invariant gluon fields from a collinear gluon field. Since the derivation is similar, here we present only the final result. The relevant term in the O(λ 2 ) Lagrangian is From this we can derive the factorization of the amplitude using the same procedure as in section 5.2.3 to get A (2),rad.

JHEP11(2019)153
where we have defined Therefore we have just shown that the double emission radiative function is completely determined by the single emission one. Since we have proven this relation at the operator level, it is true at all orders in perturbation theory. Evaluating the radiative function at tree level, we have (5.49)

Comparison to the literature
In this section we make contact with other definitions of radiative functions given in the literature. Radiative functions have been studied in the context of threshold resummation [69,70,[125][126][127][128][129], where they have been defined for the case of the emission of a single soft gluon. If partonic initial states are used, then the radiative jet functions defined in the threshold limit are most similar to what we have referred to as amplitude level radiative functions in this section, but with incoming instead of outgoing conventions. This is due to the fact that the kinematics of the threshold limit imply that collinear emissions cannot cross the cut, and therefore the factorization is distinct from the cross section level factorizations we study in section 6. In refs. [66,69,70] a radiative function 10 was defined as J µ,a (p, n, k)u(p) = d d ye −i(p−k)·y 0|Φ n (∞, y)ψ(y)j µ,a (0)|p , (5.50) with the non-abelian current, j µ,a , given by [70] j µ Here Φ n (∞, y) is a Wilson line, so that the combination Φ n (∞, y)ψ(y) is equivalent to the gauge invariant field χ n in SCET. This structure is clearly recognized to be of the same general form as the amplitude level radiative functions defined in the previous sections. Following [69], the radiative function defined in eq. (5.50) can be expanded in powers of α s and at tree level its Feynman rule reads 10 Note that in [69] this is called radiative jet function. As explained in section 1 we have preferred to keep the term radiative jet function only for objects entering the factorization at the cross section level, in analogy with the leading power jet functions. The amplitude level functions are identified as radiative functions. Since the comparison in this section is done at the amplitude level we will use the term radiative function throughout the section.

JHEP11(2019)153
Using this radiative function, the next-to-leading power corrections in the threshold limit were computed to NNLO. We would now like to show that the Feynman rule of the radiative function we have derived in eq. (5.56) matches eq. (5.52) after expanding homogeneously in the SCET power counting and taking the ultrasoft emission to be on-shell and the external quark to be purely collinear. These two assumptions translate into the following kinematics We can then kinematically expand the current of eq. (5.52) as In the first line we neglected ω n µ 2 since we are focusing on the O(λ 2 ) expansion of this object. In the third and fourth line we used the properties of the Dirac projectors P n , Pn acting on outgoing collinear spinors P n = / n/ n 4 , Pn = / n/ n 4 , 1 = P n + Pn ,ū n Pn =ū n ,ū n P n = 0 . (5.55) We can now show how this result is reproduced in our framework. We first note, that in our definition of the radiative function for the emission of a soft gluon from a quark, as discussed in section 5.2.2, we have factorized the ultrasoft derivative and ultrasoft gluon field into a soft function. Therefore, we do not have a correspondence at the level of the radiative function itself with eq. (5.54). Instead, we must contract the radiative jet function with the corresponding soft function. Evaluating this at tree level, we have We see that this indeed matches eq. (5.56). This shows, that despite the slightly different organization, the radiative functions in the two approaches are describing the same physics. While the general form of the radiative functions in [70] and those defined in this paper are similar, and indeed we have shown they give the same tree level result, the radiative functions defined in this paper differ in many aspects from the ones of [70]. In addition to several trivial differences that are conventional, namely [70] uses an incoming, instead of outgoing collinear state, and that the Lagrangian insertion occurs at y = 0, there are some more major differences, which we now elaborate on.

JHEP11(2019)153
The definitions of refs. [66,69,70] involve a four dimensional convolution (which in dimensional regularization must be extended to a d = 4 − 2 dimensional convolution) instead of the one dimensional convolution along the light cone direction derived in the effective theory. The ability to simplify the definition to involve a single variable convolution along the light cone relied on the multipole expansion in SCET, which renders the collinear matrix elements local in certain directions. The multipole expansion also implies that no expansions need to be performed after performing any integrals, and that all results are automatically homogeneous in the power counting. We believe that this is essential for achieving a true factorization. In particular, to claim a factorization, it must be that no divergences are generated by the integral in the final convolution variable, and that the convolution variable must be gauge invariant. This seems difficult to achieve if the convolution variable is a d = 4 − 2 dimensional momenta in dimensional regularization. While we are not yet able to prove that the convolutions in the single scalar lightcone variables in our formulation converge, we believe that the reduction to single variable convolutions that are not dimensionally regularized is an essential first step.
A second major difference in the definitions relates to gauge invariance. The current of eq. (5.51) is not by itself a gauge invariant object, and hence neither is the jet function of eq. (5.50). However, as shown in refs. [66,69,70] it will give a gauge invariant result for the cross section with a single soft emission. This lack of gauge invariance in a radiative function defined in the full theory is perhaps not surprising, as one is attempting to factorize a gluon emission from a quark-antiquark current, but a full theory gluon is not gauge invariant. In the context of QED, where the radiative function was originally defined [66], this issue was not present, since the photon is not charged, and therefore the radiative function is itself gauge invariant. It is ultimately this difference which makes the extension to the QCD case more difficult. To achieve a true factorization into objects which can be separately renormalized, it seems desirable that each of the objects be separately gauge invariant.
In the effective theory approach to defining the radiative functions the radiative jet and soft functions are each separately gauge invariant, since they are constructed from the gauge invariant gluon, B us(n) and quark, ψ us(n) fields which couple to the radiative functions. This leads to a fairly intricate structure of Wilson lines in the definition of the currents, which ensures their gauge invariance. Without the presence of these soft Wilson lines, it is also not clear how the current of refs. [66,69,70] can describe multiple soft gluon emissions.
Finally, another difference between the two approaches is that due to the manifest power counting in the effective theory, there are a large number of distinct field structures present in the SCET O(λ) and O(λ 2 ) Lagrangians, as compared with only the two terms present in the current of eq. (5.51). These additional terms in the SCET Lagrangian have more than two collinear fields, or have an additional P ⊥ insertion, so that they first contribute when their is an additional collinear loop. For example, with a collinear loop, we have the distinct diagrams

JHEP11(2019)153
This is forced upon us in the effective theory by manifest power counting, which does not hold for the current of eq. (5.51), which must be expanded after performing loop integrals. While this admittedly requires more distinct terms to be considered, we believe that they will have a simpler structure, and will facilitate an all orders understanding. Furthermore, it is hoped that only a subset will contribute at a given logarithmic accuracy.
The additional terms involving multiple collinear fields will be discussed in more detail in section 7 where we will consider the field structure of the complete set of radiative functions which contribute to the thrust event shape.

Cross section level factorization with radiative jet functions
In this section we show in detail how to perform subleading power factorization at the cross section level for radiative contributions in a gauge invariant manner to all orders in α s in a non-abelian gauge theory. We derive the explicit structure of the factorization for the insertion of: • A single insertion of a B us field from the L (2) Lagrangian (section 6.2.1), • A double insertion of the ψ us field from two L (1) insertions (section 6.2.2).
Here we have given the sections where the factorization is given to aid the reader. Furthermore, we consider when this insertion is on both a collinear quark, or collinear gluon leg. These are the two contributions to a factorized cross section that are non-vanishing at lowest order in α s , and allow a factorized description for the contributions that have been computed to fixed order in [40,42]. Since they describe the lowest order contributions to the cross section, their renormalization should capture the LL series. The key result of this section is a factorized expression for these contributions to the cross section as a convolution between separately gauge invariant soft and collinear factors. The convolution is in terms of either one or two one-dimensional variables, which represent the position of operator insertions along the lightcone, similar to the convolution over the observable variable in the leading power factorization review in section 6.1. After having worked out these two examples in detail, and illustrating how a gauge invariant factorization can be achieved, in section 7, we then provide an enumeration of the different radiative jet functions which contribute to the event shape thrust in e + e − → dijets to O(λ 2 ), including those which first contribute at loop level. This will illustrate the complexity of subleading power factorization in an example of interest.

Convolution structure at leading power
We begin by reviewing the well known leading power factorization for an SCET I event shape, τ , in e + e − → dijets. Here the unique leading power operator is Our focus is on the convolution structure in the observable and momentum, as it is this aspect which will change at subleading power for the radiative jet functions. Therefore we will often suppress explicit Dirac or color indices.

JHEP11(2019)153
Following the discussion of section 3, after performing the BPS field redefinition, factorizing the observable, and Fierzing the Lorentz and Dirac structure (which we will suppress, as it is not important to our current discussion), the leading power cross section can be written as This is derived by considering the factorization of the squared matrix element of eq. (6.1), as was described in section 3.4. We would now like to simplify this expression, and remove the convolution over the variable x, which couples the different functions. To achieve this, we first define Fourier transforms of each of the functions, and writing the result in terms of light cone coordinates, we can express the cross section as · S(τ us , r)J n (τ n , l, Q)Jn(τn, k, Q) .

(6.4)
We now use that Jn depends only on k − , and J n depends only on k + , which follows from the multipole expansion in the effective theory. Performing the integrals in k + , k ⊥ , l − , l ⊥ , we then find which defines the standard leading power factorization Each of S, J n and Jn are infrared finite 11 and gauge invariant. The only coupling between the soft and collinear degrees of freedom is the convolution in the physical observable τ . 11 We assume that all functions are defined with their appropriate zero-bin subtractions [130].

JHEP11(2019)153
Diagramatically, we have where the blue lines indicated collinear fields, the double green lines indicate Wilson lines, and the ⊗ indicates the convolution in τ . We have not drawn explicitly the loop corrections, as it is understood that this holds to all orders in α s . Here, for simplicity, we have only illustrated the factorization of the soft and collinear pieces of the diagrams, since it is this aspect which is relevant for the radiative jet functions. This picture will be contrasted with that for the radiative functions derived in the next section.

Convolution structure at subleading power
We now consider the subleading factorization at cross section level for the contributions coming from Lagrangian insertions and show how they lead to radiative jet functions. We discuss in detail the convolution structure for how these radiative jet functions enter the factorization.

Soft quark emission
At the cross section level, two insertions of the L (1) soft quark Lagrangian are required to give a non-vanishing contribution. At tree level, such contributions will describe single soft quark emission. Subleading power shape functions arising from two insertions of the L (1) soft quark Lagrangian were considered in [8,9,18]. They have a similar structure to the matrix elements considered in this section, although the shape functions are defined as matrix elements of B meson states, instead of vacuum matrix elements. The Lagrangian of interest for soft quark emission is L (1)BPS χψ =χ n g / B n⊥ ψ us(n) +ψ us(n) g / B n⊥ χ n . (6.8) We first consider the insertion of this operator onto a collinear quark line, which at lowest perturbative order corresponds to the emission of a soft quark from a collinear quark. To illustrate this, we use as an example e + e − → dijets with the hard scattering operator of eq. (6.1). The derivation of the factorization requires both a factorization of the spin and color structures, as well as of the momentum convolutions. The factorization of the spin and color structures is in principle a straightforward excerise, which requires, in the present case, the repeated application of Fierz identities to the matrix element The details of this are presented in appendix B, and here we focus on the convolution structure. We note that in general one must consider the presence of evanescent operators,

JHEP11(2019)153
however, since in this paper we do not consider the renormalization of the matrix elements appearing in the factorization formula, we leave this to future work. After factorizing the spin and color structure, we arrive at the following expression for the radiative contribution to the cross section which we would like to factorize into hard, jet and soft functions. The superscripts and subscripts labeling the cross section indicate that this is the O(λ 2 ) contribution arising from the radiative emission of a soft quark from the n collinear sector. Fourier transforming each of the matrix elements, we have χnψ (τ us , r 1 , r 2 , r 3 , Q) χnψ (τ n , l 1 , l 2 , l 3 , Q)g µν Using the locality of the collinear matrix elements, as discussed in our derivation of the amplitude level radiative functions (see eq. (5.12) and the surrounding discussion), we find χnψ (Qτus, r + 2 , r + 3 )J (2) χnψ (τn, r + 2 , r + 3 , Q) .
Here J (2) nψ (τ n , r + 2 , r + 3 ) defines the radiative jet function at the cross section level. Note that it is defined as a collinear matrix element, so all loop corrections are collinear in nature.

JHEP11(2019)153
Diagramatically, we have At lowest order in α s , the soft function is proportional to δ(r + 2 − r + 3 ), which simplifies the structure of the convolution to a single variable. In the presence of radiative corrections, the full convolution structure is required. As expected at subleading power, this contribution is first non-vanishing crossing the cut, namely the soft quark. Purely virtual contributions are proportional to δ(τ ), and hence leading power. In the case that the soft quark is radiated from then collinear sector, an identical factorization applies.
The ability to formulate this factorization in a gauge invariant way relies on the use of non-local gauge invariant collinear quark and gluon and soft quark and gluon fields. These operators have a highly intricate Wilson line structure, involving both soft and collinear Wilson lines situated at a variety of positions along different light cones. However, this structure is completely dictated by the symmetries of the effective theory.
As compared with the leading power factorization of eq. (6.6), there is a convolution structure in the + component of the soft momentum, in addition to the convolution in the observable τ . This couples the collinear and soft sectors in a more non-trivial way, describing the "radiation" of a soft parton from the collinear jet at a position along the light cone, as was the case at the amplitude level. Note that while we can perform the factorization of the jet and soft functions into this convolution structure, it is not a priori guaranteed that such convolutions converge. This would indicate a naive breakdown of the factorized expression, or at least that a reorganization is required. This has been explicitly observed in subleading power factorization formulae for B-meson decays [114]. In defining the jet and soft functions, we have inserted factors of the convolution variables r + 2 and r + 3 . This is done so that the lowest order divergence appears entirely in the soft function, allowing it to be extracted using standard plus distributions. We can therefore make sense of the convolutions appearing in the factorization at lowest order in α s . An understanding of the convolutions appearing in the factorization formulas in this paper is left to future work. In this paper we will not consider radiative corrections to the soft and jet functions S (2) χnψ (τ us , r + 2 , r + 3 ), J χnψ (τ n , r + 2 , r + 3 ), leaving this to future work. The renormalization of these operators, and the corresponding renormalization group evolution will resum subleading power logarithms in the thrust variable.
We can also perform an identical factorization when the soft quark Lagrangian is inserted on a collinear gluon leg. Here we take as a concrete example the case of thrust in H → gg with the leading power operator ⊥n,ω 2 · B b ⊥n,ω 1 H . (6.14)
Schematically, this can be illustrated as This takes an identical form to the case of the emission from a collinear quark, with the exception of the detailed structure of the Wilson lines and fields in the jet function.

Soft gluon emission
At lowest order in α s , we also have contributions from soft gluon emission. Due to the LBK theorem, these first arise through a single insertion of the L (2) Lagrangian. Here we will

JHEP11(2019)153
show how these can be factorized to all orders at the cross section level, for the insertion of the L (2) Lagrangian on either a collinear quark or gluon leg.
We begin by considering the case of the insertion on a collinear quark leg, corresponding to soft gluon emission from a collinear quark. We again take as a concrete process thrust for e + e − → dijets, as mediated by the hard scattering operator of eq. (6.1). For convenience, we recall the form of the Lagrangian This contributes to the cross section through the interference with the leading power hard scattering operator. As was the case at amplitude level, the factorization for the ultrasoft gluon emission is more complicated due to the presence of the ultrasoft derivative. We would like to arrange the ultrasoft derivative such that it acts only on the ultrasoft fields in the soft function. This can be done using identical arguments as given in section 5.2.1.
Flipping the derivative to act in the soft sector, and performing the Fierzing of the matrix element to obtain a factorized form, we find

JHEP11(2019)153
we find nBus (τ us , r 1 , r 2 , Q) nBus (Qτ us , r + 2 )J (2) nBus (τ n , r + 2 , Q) . Here the green lines indicate Wilson lines, and the orange spring indicates the B us(n) field. As with the soft quark radiative function, this soft gluon radiative function is first nonvanishing with a single emission crossing the cut. We have done some rearrangement, multiplying and dividing by a factor of l + 2 . This choice is to make the tree level expression for the jet function independent of l + 2 . In particular, we find a radiative jet function J (2) nBus (τ n , r + 2 ) = dl + 1 2π J nBus (τ n , l + 1 , r + 2 ) · r + 2 , (6.24) which incorporates all the collinear dynamics in the n sector from which the B us(n) field is emitted. Since it is defined as a collinear matrix element, all loop corrections are collinear in nature, and are factorized from the soft loop corrections, which are described by the soft function.
We can also consider the insertion of an L (2) on a collinear gluon line, which corresponds at lowest order to the emission of a soft gluon from a collinear gluon leg. However, here we will find something interesting. If the collinear gluon field is taken to have no label perpendicular momentum, the O(λ 1 ) Lagrangian insertion vanishes. However, at O(λ 2 ), we have a contribution from the subsubleading gluon lagrangian L (2) ng , defined in eq. (4.36). By dropping the terms proportional to the label perpendicular momentum, the relevant insertion comes from the three field operator

JHEP11(2019)153
However, we note that (6.26) For the case of dijet production, where the two jets are back to back, all Wilson lines in the soft function are along the light cone directions n orn. Therefore, we see that the emission of a soft gluon from L ng gives a vanishing contribution to the cross section at lowest order in α s . This fact was used for the calculation of the leading log series in ref. [59]. We will therefore not consider its factorization explicitly, since in this section we have only considered the factorization for contributions at lowest order in α s . Beyond O(α s ), there are other contributions that will be discussed in more detail in the next section.

Discussion
In the section we have derived a factorized expression in terms of a convolution of gauge invariant soft and collinear matrix elements for the emission of a gauge invariant soft quark or gluon field at the cross section level. The main new feature at subleading power, are additional convolutions in the momentum passed between the soft and collinear sectors. Due to the multipole expansion in SCET this can be simplified to a single component k + i , which represents the position along the light cone. This also follows directly from the SCET power counting since only the +-momenta are the same size. To be able to separately renormalize the soft and collinear matrix elements, allowing for a resummation of subleading power logarithms, it is essential that both matrix elements are gauge invariant, and that no additional divergences appear in the single variable convolution. In a nonabelian gauge theory, the ability to formulate a soft emission in a gauge invariant manner is highly non-trivial since soft partons carry non-abelian charge, and it relies on our use of non-local gauge invariant soft quark and gluon fields, and an understanding of the all orders Lagrangian describing their interactions.
Some of the factorized expressions presented in this section are similar to those which have appeared in the B-physics literature [9]. In this case, instead of being vacuum matrix elements, the soft function is a matrix element of B meson states. These were studied for the case of heavy to light decays where the heavy quark is treated using HQET. Nevertheless, their structure as non-local operators with insertions along the light cone is similar. This can be viewed as another advantage of using the operator based approach, namely since the SCET Lagrangians are universal, the same structures will appear in a variety of distinct physical processes allowing for a universal framework. As an example, although we have used for concreteness here the case of dijet production, since the insertions appear in a single collinear sector, this could be extended to N-jet production in a straightforward manner.
In this section we have explicitly worked out the structure, including the Lorentz and Dirac structure for two examples, namely single soft quark and soft gluon emission, with the goal of showing in detail how the factorization with Lagrangian insertions can be performed, and how it gives rise to radiative functions at the cross section level. These

JHEP11(2019)153
two examples were chosen, since they contribute at lowest order. It should be clear that a similar factorization can be performed for an insertion of any term in the subleading power Lagrangians of eqs. (4.1) and (4.2), or for multiple Lagrangian insertions. In section 7 we will perform a systematic classification for the case of thrust in e + e − → dijets of the field structure of all possible radiative functions.

Classification of radiative contributions for dijet event shapes
Having shown how we can achieve a factorization of radiative contributions in terms of a convolution of gauge invariant soft and collinear matrix elements, in this section we extend this to a complete classification of the radiative contributions for an SCET I event shape observable. For concreteness, we will consider thrust in e + e − → dijets, however thrust for H → gg would have a very similar form. Unlike in the previous section, where we worked out in detail the factorization for those contributions that contribute at lowest order in α s , here we will not explicitly perform the Fierzing and factorization. Instead, we will content ourselves with an understanding of the field structure, leaving the details of the factorization and resummation to future work. Nevertheless, a study of the field content provides valuable insight into the structure and complexity of the factorization at subleading power. Despite the large number of different field structures in the Lagrangians of eqs. (4.1) and (4.2), we will see that due to the symmetries of the problem many terms can be shown not to contribute to all orders in α s . For example, we will show that there are no O(λ) ∼ O( (τ )) power corrections, as expected from perturbative calculations. A much larger reduction in the number of jet and soft functions is achieved if one works to a fixed order in α s , for example O(α 2 s ), and a study of the field content is sufficient to understand these aspects of the factorization.
In section 7.1, we show that O(λ) radiative contributions to the cross section vanish to all order in α s , and in section 7.2 we classify the O(λ 2 ) contributions. For those operators that can contribute, we will explicitly give the field structure and contributing diagrams to O(α 2 s ).

Vanishing at O(λ)
At O(λ), the only possible radiative contribution is an L (1) insertion into the leading power hard scattering operators. If this L (1) insertion involves the conversion of a collinear quark to an ultrasoft quark, it will vanish, since it will give rise to a vacuum matrix element involving an odd fermion number. To show that the remaining possible L (1) insertions also vanish, we note that by power counting they all involve either an ultrasoft gluon field, or an ultrasoft derivative operator, but not both. First consider the case of an ultrasoft derivative operator. With a single Lagrangian insertion into the leading power operators, after performing the BPS field redefinition, the ultrasoft derivative operator must act on collinear fields, and is therefore factorized into the jet function. We now use the dimensional regularization rule for residual momenta

JHEP11(2019)153
where (q r ) j denotes positive powers of the q − r and q ⊥ r momenta, which are the only residual momenta which appear in the subleading power Lagrangians, to show that this vanishes. After eliminating these terms, we can simplify the structure of the L (1) Lagrangian to L (1)BPS n ∼ B n⊥ B n⊥ B us(n)⊥ P ⊥ +χ n χ n B us(n)⊥ P ⊥ + B n⊥ B n⊥ B n⊥ B us(n)⊥ +χ n χ n B n⊥ B us(n)⊥ , (7.2) each term of which involves a B us(n) field. After performing the factorization, the B us(n) field can be factorized into the soft function. However, the soft function must be rotationally invariant about the n−n axes. Since the only other objects which appear in the soft functions at this power are Wilson lines, this implies that it is only possible to form a rotationally invariant soft function from then · B us(n) , or n · B us(n) components of the field. However, only the B us(n)⊥ component appears at this power. Therefore, there are no O(λ) power corrections from Lagrangian insertions. In previous papers [61,62], we have shown that O(λ) contributions from hard scattering operators vanish for both a qq and gg current, as do subleading power corrections to the measurement operator. Therefore, the present analysis provides a complete proof that there are no O(λ) power corrections to all orders in α s .

Contributions at O(λ 2 )
The O(λ 2 ) radiative contributions to the cross section do not vanish. As discussed in section 3, in addition to an insertion of the O(λ 2 ) Lagrangian, one must also consider two insertions of the O(λ) Lagrangian, as well as an insertion of the O(λ) Lagrangian into an O(λ) hard scattering operator. 12 The complete expression for the radiative contribution is therefore dσ dτ which is a subset of the complete set of O(λ 2 ) contributions given in eq. (3.21). We will consider the different contributions from each of these cases in turn. For simplicity, we will not separately discuss the different possible positions of the final state cut, since the operator structure is the same in all cases. 12 One must also consider subleading corrections to the measurement function. However, as discussed in section 3.2, the subleading power corrections to the measurement function first enter at O(λ 2 ), and therefore if we are working only to O(λ 2 ), we do not need to consider the interference of the subleading power measurement function with Lagrangian insertions. At higher powers, this contribution could be treated in a similar way to the contributions discussed in this paper.

Contributions from O(λ) hard scattering operators
We first consider term the first two lines of eq. (7.3) involving a single O(λ) hard scattering operators, O (1) , and a single O(λ) Lagrangian insertion, L (1) . A detailed derivation of the O(λ) hard scattering operators for e + e − → dijets was given in [61]. The operators are of two types, and either involve a collinear quark and a collinear gluon field in the same collinear sector, or two collinear quark fields in the same collinear sector. We can use arguments identical to those presented in section 7.1 to show that these terms do not contribute to all orders in α s . In particular, since the hard scattering operators at O(λ) do not involve additional ultrasoft fields, the rotational invariance of the soft function, combined with the structure of the O(λ) Lagrangian, which involves only B us(n)⊥ fields guarantees that all contributions involving ultrasoft field insertions vanish. Similarly, all contributions involving insertions of an ultrasoft derivative from the Lagrangian vanish for the same reasons as described in section 7.1.

Contributions from O(λ 2 ) Lagrangian insertions
Next we consider the third line of eq. (7.3), involving the T-products of leading power hard scattering operators and L (2) . When considering contributions to the cross section, the structure of the L (2) Lagrangian can be simplified, since one cannot have a single insertion involving an ultrasoft fermion, or involving just ultrasoft derivatives acting in the collinear sector. This reduces the structure down to 10 distinct field structures L (2)BPS n ∼ B n⊥ B n⊥ B us(n) B us(n) + B n⊥ B n⊥ ∂ us B us(n) + B n⊥ B n⊥ B n⊥ B n⊥ B us(n) (7.4) + B n⊥ B n⊥ B us(n) P 2 ⊥ + B n⊥ B n⊥ P ⊥ B n B us(n) +χ n χ n B us(n) B us(n) +χ n χ n ∂ us B us(n) +χ n χ n B n⊥ B n⊥ B us(n) +χχB us(n) P 2 ⊥ +χ n χ n P ⊥ B n B us(n) , the explicit form of which was given in eq. (4.37). Unlike at O(λ), these contributions do not vanish. For each of these different contributions, it is easy to understand at which order they can first contribute. Those with a P ⊥ insertion will first contribute at one-higher order, as they require an additional collinear emission to provide a non-zero ⊥ momentum. In particular, we see that while there are a large number of different radiative functions that are required for an all orders descriptions, many first contribute at high loops. Considering only those terms which contribute to O(α 2 s ), we have only 6 terms L (2)BPS n ∼χ n χ n ∂ us B us(n) +χ n χ n B us(n) B us(n) + B n⊥ B n⊥ ∂ us B us(n) + B n⊥ B n⊥ B us(n) P 2 ⊥ +χχB us(n) P 2 ⊥ +χ n χ n P ⊥ B n B us(n) . (7.5) The first term contributes at O(α s ). The second term contributes at O(α 2 s ), and has been analyzed in detail in section 6.2.2 where we worked out the complete structure of the factorization, and its tree level diagram is shown in table 3. The remaining terms first contribute at O(α 2 s ), and representative diagrams, along with the field structure of the factorization are given in table 4 for those with a single B us(n) field, and in table 5 for those with two B us(n) fields.

Contributions from [O(λ)] 2 Lagrangian insertions
At O(λ 2 ) one receives contributions from two L (1) insertions as given by the last two lines of eq. (7.3). The non-vanishing contributions can have either one, or two soft fields. Since the L (1) has the field structure L (1)BPS n ∼ B n⊥ B n⊥ ∂ us P ⊥ +χ n χ n ∂ us P ⊥ + B n⊥ B n⊥ B us(n) P ⊥ (7.6) +χ n χ n B us(n) P ⊥ + B n⊥ B n⊥ B n⊥ B us(n) +χ n χ n B n⊥ B us(n) + B n⊥ B n⊥ B n⊥ ∂ us +χ n χ n B n⊥ ∂ us +χ n B n⊥ ψ us(n) , this gives rise to a large number of distinct possibilities. While it is possible to write the radiative functions for all the different possible cross terms of eq. (7.6), here we focus only on those that contribute up to O(α 2 s ). Note that the two insertions do not need to be in the same collinear sector. However, those in distinct sectors first contribute at 3 loops. We therefore do not draw them explicitly.
The terms in the L (1) Lagrangian contain at most one ultrasoft field. Since radiative contributions where both L (1) insertions involve no collinear fields vanish, the double L (1) insertions can contain either one or two ultrasoft fields. Except for those terms involving the ultrasoft quarks, all terms in the L (1) Lagrangian involve either > 2 collinear fields, or a P ⊥ operator. In both cases, these first appear at O(α 2 s ). The only tree level contribution from L (1) insertions is the that of the soft quarks, which was discussed in detail in section 6.2.1 where the complete factorization structure was worked out, and whose leading order diagram is shown in table 3. In table 4 we summarize the field structures of the different radiative terms with a single emission arising from two L (1) insertions, as well as a representative diagram at O(α 2 s ). In all cases, these can be viewed as a propagator correction, marked with a cross, and a B us(n) emission from L (1) . The contributions with two B us(n) fields are shown in table 5, along with representative diagrams at O(α 2 s ). Note that here, and throughout this section, we will explicitly draw soft fields in the radiative jet functions to illustrate where they are attached. These fields will of course be factorized into the soft functions.
It is interesting to briefly consider the simplification of the L (1) · L (1) terms in the threshold limit. The threshold limit has received much attention in the power corrections literature [47,69,70,[125][126][127][128][129][131][132][133][134] due to the large amount of available perturbative data [48-50, 52, 135-138]. Up to crossing, the production of a color singlet Drell-Yan pair in a proton proton collision is identical to the case of e + e − → dijets considered here, and therefore power corrections in (1 − z) can be considered using the same formalism and operators. In the case of threshold production, collinear particles are kinematically forbidden from crossing the cut into the final state. This greatly reduces the number of possible contributions at O(α 2 s ). In particular, since for the L (1) · L (1) to contribute at O(α 2 s ), the collinear gluon must cross the cut, these terms will not contribute in the threshold limit to O(α 3 s ).

Discussion
The analysis of this section illustrates that the general structure of the Lagrangian contributions is quite complicated. In particular, while there are only several contributions at
lowest order in α s , there are a large number of contributions which first appear at higher orders, coming from Lagrangian insertions which involve multiple collinear fields. Indeed, there are terms whose first non-vanishing contribution is at 4-loops. These contributions do not seem in general to be related to the renormalization group evolution of those operators that contribute at lowest order, and thus appear to complicate the general structure. We believe that the leading logarithmic contributions are captured by the renormalization group evolution of the two radiative functions that contribute at lowest order, which are given in table 3, and that the Lagrangian insertions involving higher numbers of collinear fields are only required to reproduce subleading logarithms. This has been explicitly verified in the subleading power leading logarithmic resummation [59] for thrust in H → gg in pure Yang-Mills theory (i.e. without quarks), and here a proof of this statement was given under the assumption that the convolutions appearing in the factorization formula converge. However, we do not currently have a more general proof of this statement. This expectation comes from the fact that collinear loops from these operators involving multiple collinear fields are less singular than the identical collinear loops involving collinear Wilson line emissions. If this is true, then the leading logarithmic resummation may take a simple form, but beyond leading logarithmic order an increasing number of operators will be required. We plan to develop the leading logarithmic renormalization group evolution of these leading power operators in a future paper, and by comparing with explicit calculations it should be possible to understand the role of the operators involving additional collinear fields.
A key advantage of our approach is that the radiative functions are derived from factorized matrix elements of the SCET Lagrangians, which describe the interactions of the gauge invariant quark and gluon fields to all orders. This allows us to provide a complete classification of all radiative functions that will contribute at a given power in the expansion. The number of such operators at each given power is finite, and soft and collinear emissions to all orders in α s are described by the intricate Wilson line structure of the operators, dictated by the symmetries of the effective theory. As one particular example, we introduced a radiative function describing the emission of a soft quark. As was shown in the calculation of subleading power thrust [40][41][42] contributions from soft quarks give a leading logarithmic contribution at subleading power. Soft quarks are also required to compute subleading power corrections to the threshold limit for Drell-Yan production in the qg channel. We performed a classification of all the radiative functions which contribute for the event shape thrust in e + e − → dijets, and derived in detail the structure of the radiative functions which contribute at tree level. This provides the final missing ingredient in achieving factorization and resummation of event shapes at subleading power.
Combined with the subleading power operator bases forqΓq [61,63] and gg [62], and the analysis of the subleading power measurement function [61], all the sources of subleading power corrections in SCET (as were summarized in section 3.4) have now been characterized in detail. It will now be of significant interest to derive a complete subleading power factorization theorem for an event shape observable, and to study the renormalization group evolution of the different components, a direction which we intend to pursue in future work. The renormalization of higher twist matrix elements has been studied [5, 139-

JHEP11(2019)153
149], as has that of subleading power operators in B physics [150,151], and of several power suppressed operators for e + e − → dijets [87,88] and N -jet production [64,65]. We anticipate that combining all these ingredients will allow for the resummation of subleading power corrections for event shape observables.

JHEP11(2019)153
Reinstating the time ordering and measurement function, this corresponds to the result given in section 6.2.2.
Open Access. This article is distributed under the terms of the Creative Commons Attribution License (CC-BY 4.0), which permits any use, distribution and reproduction in any medium, provided the original author(s) and source are credited.