{'fraction_non_alphanumeric': 0.03536419492519113, 'fraction_numerical': 0.009403874114722575, 'mean_word_length': 4.243206688798721, 'pattern_counts': {'":': 0, '<': 0, '<?xml version=': 0, '>': 2, 'https://': 0, 'lorem ipsum': 0, 'www.': 0, 'xml': 0}, 'pii_count': 0, 'substrings_counts': 0, 'word_list_counts': {'cursed_substrings.json': 2, 'profanity_word_list.json': 0, 'sexual_word_list.json': 0, 'zh_pornsignals.json': 0}}
The Large-x Factorization of the Longitudinal Structure Function

A leading-twist factorization formula is derived for the longitudinal structure function in the x -->1 limit of deeply inelastic scattering. This is achieved by defining a new jet function which is gauge independent and probes the transverse momentum of the struck parton in the target. In moment space, terms of order (\ln^k N)/N, which are the leading ones for F_L, are shown to be resummable through the cusp anomalous dimension \gamma_K and the anomalous dimension \gamma_{J^\prime} of the new jet function. This anomalous dimension is computed to O(\alpha_s). The suggested factorization for F_L reproduces the fixed order results known to O(\alpha_s^2). The general ideas for resumming the terms of order (\ln^k N)/N in moment space may be extended to the other structure functions and to other inclusive processes near the elastic limit.


Introduction
Typical applications of QCD perturbation theory involve processes with two mass scales, the momentum transfer Q and the scale of strong interactions Λ QCD with Q 2 ≫ Λ 2 QCD . Within this large family of hard processes a particular class is of great theoretical and phenomenological importance. These are the quasi-elastic processes, characterized by a third scale M 2 , the invariant mass of the final states, with Q 2 ≫ M 2 . Here perturbation theory can be used to study the transition from short-distance dynamics (M 2 ≫ Λ 2 QCD ) to long-distance dynamics as the elastic limit is approached (M 2 → Λ 2 QCD ). Deeply inelastic scattering (DIS) structure functions at large Bjorken-x, Drell-Yan production near the partonic threshold and thrust in e +e − annihilation near the two-jet limit are all characteristic examples of quasi-elastic processes. In this paper we focus on the large-x DIS structure functions for definiteness, but we expect our discussion to be generalizable to all quasi-elastic processes.
Near the elastic phase-space region x → 1 of DIS the coefficient functions receive corrections in 1−x. For F 2 the leading power in x corrections are terms of the form α n s [ln m (1−x)/(1−x)] + , with m ≤ 2n − 1. The formalism of resumming this class of corrections has been developed in Refs. [1,2,3] and recently reviewed in Ref. [4]. A key step is the factorization of the perturbation theory diagrams at large x in terms of a hard amplitude, a universal soft function, and a jet function. We call this factorization "light-cone expansion", distinct from the Wilson operator product expansion. The soft function and the jet functions contain singular pieces like ln k (1 − x)/(1 − x) + and 1/(1 − x) + respectively, which are resummed. All these corrections correspond to O(α n s ln k N) terms in moment space. The question we consider here is whether a formalism beyond the light-cone expansion can be constructed to resum terms that do not correspond to power singularities in x and are still leading power in Q. A class of such terms are the simple logarithmic corrections of order α n s ln m (1 − x), which in moment space generate O(α n s ln k N/N) power suppressed contributions. The cleanest context where such corrections are encountered is in the longitudinal structure function F L . F L , as distinct from F 2 , does not contain any singular distributions in x. Its leading power in x corrections are all logarithmic and they start at order α 2 s . A phenomenological motivation for studying the resummation of F L near the elastic limit comes from considering the ratio R of the longitudinal and the transverse cross sections : where F i , i = 2, L denote the structure functions, M is the target mass and x = −q 2 /(2p · q). The longitudinal structure function F L vanishes due to helicity conservation in the naive parton model and thus R provides us with a clean test of perturbative QCD. There already exist some experimental data [6] on this quantity and more is expected in the near future. Leading O(α s ) and next to leading O(α 2 s ) perturbative corrections to F L have been calculated [7]- [9] including the target mass corrections. Comparison with the data indicates that the theoretical predictions lie systematically below the experimental points, especially in the large x region. This may be indicative of rather large higher twist corrections [10]. In order to shed more light on this issue it is worth studying the effects of resummation at the leading twist level.
The resummability of quasi-elastic corrections that do not correspond to singular distributions was considered in Ref. [5]. There it was argued that final state power suppressed corrections of order α n s (1 − x) p ln k (1 − x), with p ≥ 0, can be identified and resummed through an expansion of the final state interactions in terms of higher dimensional bilocal operators.
The expectation values of these operators define classes of new jet functions that generalize the jet function of the light-cone expansion. In this new expansion operators of different dimensionality all contribute to the leading power in Q (leading twist) after integration over the final state phase-space. However, higher dimensional operators produce terms suppressed by higher powers of 1 − x near the elastic limit x → 1. In Ref. [5] rules were given for identifying independent sets of these higher dimensional operators. Once they are constructed, their anomalous dimensions can be used to resum the power suppressed logarithmic corrections, just like in the case of the light cone expansion.
Within the context of the above formal construction, in this work we consider the longitudinal structure function in the elastic limit as a specific example of a quasi-elastic process. We shall show that by analyzing the perturbative diagrams in the appropriate infrared limits we, indeed, obtain the operator expansion suggested in Ref. [5]. We accomplish this not by dimensional arguments but through a diagrammatic analysis by explicitly considering the subleading terms in the eikonal approximation for the final state. The result is a factorization formula for F L that is used to resum in moment space all the O(ln k N/N) corrections. We find it remarkable that such power suppressed corrections can be treated systematically using the standard methods for analyzing the infared region of perturbation theory. Simple logarithms of 1 − x come from the pinch singular surfaces of Feynman diagrams. This observation, we believe has not been considered in the literature so far.
The general ideas for resumming terms of order ln k N/N in moment space can be extended to other structure functions and to other inclusive processes near the elastic limit. For example, it has been noticed [11] in the context of Higgs production that such corrections can be phenomenologically important. They may also be of relevance in renormalon phenomenology [12]. Thus, the analysis of the resummation of the O(ln k N/N) corrections is easiest done for the case of the longitudinal structure function, but as discussed, the applications go well beyond that to phenomenologically more interesting situations which we hope to address in the future.
The structure of this paper is as follows. In section 2 we set notations and review known fixed order results. In section 3 we briefly summarize the light-cone expansion for F 2 , which is used as a basis for what follows. Then we go on to show how the new jet function, characteristic of F L , arises from the infrared singularities of Feynman diagrams, thus explicitly constructing the relevant jet operator of Ref. [5]. Here we also study the renormalization group properties of the functions entering the F L factorization formula. In section 4 we present an explicit one-loop calculation of the anomalous dimension of the new jet function, which is subsequently used in the F L resummation formula. The resummed coefficient functions are found in agreement with the fixed order results. We summarize our conclusions in the final section. Certain useful formulae are listed in the Appendix.

Notations and known results
We consider the DIS structure functions F 2 , F L of a massless quark of electric charge Q f = 1 and with collinear divergences dimensionally regularized in D = 4 − 2ǫ dimensions. Their perturbative expansion is The DIS quark tensor for an electromagnetic probe is defined as where P S n is the Lorenz invariant phase space for n particles in the final state and M L,R are the left and right production amplitudes. Structure functions are obtained via the projections F r = P µν r W µν , for r = 2, L, given at D = 4 by The normalization in Eq. (2) is chosen so that F . The light cone is defined by the two dimensionless vectors v µ andv µ normalized as v ·v = 1 and we fix the v µ (+) direction to be parallel to the incoming quark momentum p. Then thev µ (−) direction is set by the vectorq = xp + q. We denote by p f the total final state momentum, p f = p + q. The light cone (Sudakov) decomposition of the external momenta is and for every momentum vector k µ The two massless particle phase space in this basis and for azimuthically symmetric integrands becomes dP S 2 (k 1 , k 2 ) = 1 8π with the Sudakov variable β 2 defined as β 2 = (v · k 2 )/(v ·q). The three massless particle phase space in terms of Sudakov variables is given in the appendix. At high momentum transfer, Q 2 ≫ Λ 2 QCD , and fixed x all structure functions satisfy the Wilson operator product expansion. At leading twist this can be conveniently written as a factorization theorem in terms of the convolution symbol Specifically for F L we have with µ the factorization scale and f the quark distribution function of the quark target. The coefficient functions are determined given a factorization scheme, which will be the MS scheme for this paper. The coefficient functionC L has been calculated to O(α 2 s ) and the results are [8,9,13] (1) L (x, 1), with the well known functions andC (2) L (x, 1) = In the above expression only terms that are singular in the x → 1 limit have been retained. It is these logarithms of 1 − x that are the focus of the present study.
3 The factorization formula for F L We shall show that F L satisfies a novel factorization formula in the x → 1 limit which can be used to sum the singular ln k (1 − x) contributions in the coefficient function, (see Eq. (12)). This factorization formula involves a new jet function that contains these leading corrections to F L in the large x limit, which in moment space are of O(ln k N/N), N being the moment variable. Such contributions are also present in F 2 , although they are subleading to the ones that are of O(ln k N) in moment space.

The new jet function
Before presenting the factorization formula for F L in the x → 1 limit it is worth reviewing the light-cone expansion for F 2 , since the latter will serve as a prototype for the former. The leading power in x factorization of F 2 can be cast in the following form [1] is the hard scattering function on either side of the final state cut. |H 2 (Q 2 )| 2 is by construction UV dominated, i.e. all its internal lines are far from their mass shell. V is the universal soft radiation function defined as FT (1) −wp stands for the Fourier transformation along the light-cone direction −v. The Wilson line operator starting in the far light-cone past and defined along the direction v µ is given by Its perturbative expansion generates the Feynman rules for the eikonal approximation. φ is the light-cone parton distribution function. Its operator definition, for the scattering of quarks, is This distribution function is not identical to f used in Eq. (8), but they both share the feature that the only singularities they contain come from collinear propagation along the initial (v) light cone direction. The role of V −1 in Eq. (4) is to remove soft contributions, which would otherwise be double-counted in both φ and V . Finally, J is the current jet function that describes the stream of almost collinear particles with invariant mass (x ′ − x − w)Q 2 moving along thev light-cone direction. J is defined as a gauge invariant and soft subtracted fermion propagator where FT (4) q+zp denotes the four-dimensional Fourier transformation with respect to the total jet momentum q + zp.
We note that the factorization theorem of Eq. (1) is written in terms of UV renormalized effective non-local operators. When one analyzes the contributions to the characteristic functions |H 2 | 2 , V , J and φ coming from classes of graphs of the full theory, it is often the case that the assignment to soft and jet functions of these contributions is gauge dependent with gauge independence occuring only for their convolution, i.e. the final result. Here, however, the characteristic functions are defined in terms of gauge invariant operators, thus we can classify singularities and enhancements as soft or collinear in an unambiguous way. The cost for this is the need to perform infrared subtractions, as denoted for instance by the denominator V in Eq. (5) for the jet function.
In Ref. [5] it was argued that F L satisfies the following factorization formula: The factors V and φ are the same as in the F 2 case, see Eqs.
(2), (4). The new jet function J ′ , is defined as the expectation value of a bi-local operator that probes the transverse momentum of the outgoing struck quark. Explicitly, The conventional normalization factors of Eqs. (2), (3) have been included in this definition. The operator defining J ′ is the only gauge invariant bi-local operator that is non-vanishing by the massless quark equations of motion or the F L projection of Eq. (3), and although it is of higher dimension than the corresponding one that defines J, Eq. (5), it contributes to leading twist. All the terms that it generates are of order (1 − x) 0 , i.e. only logarithms. It is therefore the appropriate operator for resumming the leading logarithmic quasi-elastic corrections in F L . In Ref. [5] the new jet operator was identified by requiring it satisfied all the above properties.
In the next subsection we show that it actually comes out naturally from the infrared properties of the Feynman diagrams that contribute to F L at leading twist. The jet function J ′ is depicted in fig. 1(a), and the factorization expression for F L is depicted in fig. 1(b). We emphasize that the factor V ⊗ φ is common to every DIS observable, and that the only difference between F 2 and F L at their respective leading power in x level is that they are different Lorentz projections of the same tensor W µν (p, q), hence they differ in the jet functions entering their respective factorization formulae.

Operators beyond the eikonal approximation
In order to prove the factorization formula in Eq. (6) in this subsection we follow the constructive approach. We analyze the perturbation theory diagrams and show how the operator in Eq. (7) emerges once one goes beyond the eikonal approximation for F L . Infrared singularities and ensuing corrections arise from characterisitc regions in momentum space called pinch surfaces. These are the same for all observables, therefore the classification of internal lines in an F L diagram are the same as in any other observable, i.e. we have collinear, soft and ultraviolet lines. The factorization of the soft lines from the incoming v-jet is as in F 2 . We, therefore, concentrate on the factorization ofv-collinear lines from the initial state v-jet since this is where the difference between F L and F 2 lies.
Consider the configuration shown in fig. 2, in which the depicted gluon belongs to thev-jet. Clearly the fermion line carrying momentum (p − l) is far off-shell and belongs to the hard scattering subdiagram in the notation of reduced diagrams. If, however, we were to factorize the collinear gluon as usual by eikonalizing the fermion propagator then this, in general leading, IR contribution would be projected out in F L because of the contraction with p µ at the current vertex. Note that also in the case where the gluon is soft, the leading eikonal contribution is similarly projected out. Therefore, we must go beyond the eikonal approximation.
The graph in fig. 2 is proportional to We define the transverse projector ω α β = δ α ⊥β and decompose every γ matrix in light cone components as We also define the effective vertex It enjoys the following properties and Then the factor in Eq. (8) becomes If only one gluon from thev-jet is connected to the hard scattering, as in fig. 2, then only the first term in (13) survives the massless Dirac equation. Recall that the leading spin structure , which is the part of the operator in Eq. (7) that sits to the left of the final state cut.
To show that the operator in Eq. (7) gives the leading fragmentation contributions for F L to all orders we proceed inductively by considering the coupling of twov-jet gluons to the hard scattering, fig. 3. Because the hard scattering subgraph is UV dominated, its Dirac structure is again γ µ hence for F L the same / p factor enters at the current vertex. Then the graph of fig. 2 is proportional to where terms that will vanish by the massless Dirac equation have been dropped. There is a hierarchy among the terms of the above expression once this factor is imbedded into a cross section diagram. The leading contribution comes from the first term and with its factor in square brackets eikonalized, i.e. from All the other terms in Eq. (14) are suppressed by at least one power of (1 − x) in the cross section, relative to the contribution that has been retained. We observe that Eq.
. This is how we show constructively, i.e. from analyzing the diagrams, that the new jet function J ′ , as defined by the operator of Eq. (7), will contain the leading, O((1 − x) 0 ), terms in the longitudinal cross section. Note that the operator in Eq. (7) contributes at leading twist level. It is also the only operator that generates O((1−x) 0 ) contributions from the fragmentation of the final state jet. It can be shown that any other operator with higher number of covariant derivatives than in J ′ will genarate contributions Once the jet function J ′ is defined in terms of the expectation value of an operator it can be calculated in perturbation theory. However, the momentum integrals that are obtain in this calculation will also contain contributions in which the gluons coupling to thev-jet are soft. It is a well known (and non-trivial) statement that such contributions can be expressed in terms of expectation values of Wilson lines. Such contributions are removed by the denominator in Eq. (7) so that J ′ is by definition free from soft contributions. The latter ones will be contained only in the function V , Eq. (2). In fixed order perturbation theory this amounts to performing IR subtractions from the expressions generated by the numerator in order to obtain a J ′ jet function that contains only enhancements from final state fragmentation.
There is also one more subtlety that arises when dealing with the function J ′ , which is again characteristic of F L . The UV renormalization of the effective operator that defines J ′ captures enhancements that are due to the near collinear fragmentation of the outgoing jet. These contributions in the full theory (i.e. before factorization) are generated by genuine collinear infrared singularities. Although these singularities cancel upon adding real and virtual contributions in a full calculation, their finite terms remain and these are captured by the factorization theorem. So, collinear enhancements in the full theory are associated with ultraviolet singlularities in the effective jet operator and its ensuing anomalous dimension. In addition to these singularities, though, there are genuine ultraviolet singularities in the full theory that lead to the renormalization of its parameters, i.e. the coupling constant. Such ultraviolet singularies can also be present in the effective operator and these are not related to collinear enhancements but to the renormalization of the coupling constant in the jet function. If one is not careful to subtract these genuine ultraviolet singularities from the renormalization of the operator then one is in danger of double counting them when the renormalization group analysis is performed. In other words, the UV singularities of J ′ generate both the collinear enhancements and the renormalization of the effective coupling of the jet. This effect is not seen in low orders in F 2 because the jet function J starts at O(α 0 s ), whereas in F 2 it starts at O(α s ). We shall return to this issue when we construct explicitly the anomalous dimension of J ′ .

The RG equation for the new jet function and exponentiation
The renormalization group equation for the jet and soft functions as well as the exponentiation of the logarithmic enhancements near the elastic limit are both consequences of the factorization formula. Therefore, we rewrite it here by making explicit the arguments that the participating functions depend on. Callp the total momentum of thev-jet and (1 − z) the longitudinal momentum fraction in the v direction that the jet carries into the final state. From the kinematics, x ≤ z ≤ 1, andp 2 = (1 −z)Q 2 /x is the small invariant mass of the jet. The jet has been defined with respect to thev direction hence the large scale on which it depends is (p · v) 2 . Then the F L factorization formula reads Recall that p f is the total momentum of the final state and that all functions are renormalized at scale µ 2 . The only singularities present are the collinear ones with respect to the initial v direction. In moment space the above convolution becomes a simple product: where the Mellin transformation has been defined as usual and similarly forṼ andφ. Note that we have not included the overall 1/N factor into the definition ofJ ′ . This function will contain only the logarithmic in N terms. Also note that the UV dominated hard scattering function |H L | 2 contains no dependence on N. Its presence is necessary, however, to render the right-hand side of Eq. (17) independent of µ 2 . Following the approach of Collins, Soper and Sterman, (for review and references see ref. [4]) we can write down the two equations that the jet function satisfies. The first is the covariance equation Since F L is invariant upon variation of the light cone directions with respect to which the jet has been defined, the jet variations K and G are cancelled by the corresponding variations of the soft functionṼ and the hard function |H L | 2 respectively. The dependence ofṼ on the light cone directions has been left implicit in Eq. (17). The quark distribution functionφ is by definition independent of the final state directionv so it does not play a rôle in Eq. (19). The function K is defined at scale µ 2 and depends only on the common arguments ofJ ′ and V . The sum K + G is renormalization group invariant. The renormalization group equation with the well known cusp anomalous dimension in MS scheme [14] γ The second equation that the jet function satisfies is its own renormalization group equation We employ MS scheme which is a massless scheme, so the anomalous dimensions γ K and γ J ′ can depend only on the renormalized coupling. The solution of the covariance equation (19) is where we have used Eq. (20) and have defined with β(α s (µ 2 )) = d d ln µ 2 α s (µ 2 ) = −α s The solution of both the covariance and the RG equation for the jet J ′ is and the boundary condition for the jet function is (setting (p · v) 2 = Q 2 ) The renormalization group invariance of F L imposes the condition where the anomalous dimensions of all the functions involved are defined just as in Eq. (22). The above relation allows us to writeF L , Eq. (17), in a variety of exponentiated forms depending on which anomalous dimensions we wish to include in the exponent and what boundary conditions we choose for the factors H L ,J ′ ,Ṽ andφ. The form that is suitable for our purposes is the one that involves the jet anomalous dimension γ J ′ . Explicitly, We shall be referring to this equation as the Sudakov exponentiated form ofF L . By inspecting the above expression it is easy to see where the logarithmic in N corrections come from and how they are organized. As usual, the leading ln N terms arise from the soft regions and are contained in the function A. The logarithmic corrections due to the fragmentation of the final state jet are contained in γ J ′ and the ones due to the renormalization of the coupling constant in the jet function are in beta function term. The beta function term arises from the evolution of the boundary condition from the low scaleJ ′ (1, 1, α s (Q 2 /N)) to the high scaleJ ′ (1, 1, α s (Q 2 )). It is worth empasizing again here that in F L even the lowest order contribution cannot be normalized unambiguously because it starts at O(α s ). Since the renormalizability of the theory determines the structure of the beta function term in the exponent, we can absorb it into the anomalous dimension γ J ′ as while keeping the boundary condition at the high scale, i.e.J ′ (1, 1, α s (Q 2 )). The practical advantage of this simple redefinition is that when we compute the anomalous dimension of the operator in Eq. (7) we do not have to subtract out the UV poles that renormalize the coupling as was suggested at the end of the previous section. In this fashion we compute γ ′ J ′ and not γ J ′ . Once the distinction between the two is made as in Eq. (31), there is no danger of overcounting. This is the approach that we follow in the next section. Not all logarithmic in N dependence is contained in the exponential of Eq. (30), but it resides also in theṼ ·φ factor. This can be captured too, if the corresponding anomalous dimensions are introduced in the exponent. We shall leave this factor as is, because it is common toF L andF 2 and contains no further information that is characteristic of the longitudinal structure function only. Our final objective is to connect the exponentiated form ofF L , Eq. (30), with the operator product expansion factorization as given by Eq. (8). In the MS scheme the coefficient function in moment space and at scale µ 2 = Q 2 is given bỹ The regularized quark distribution functionf (N, Q 2 , ǫ) can also be written as an exponential by simply integrating its DGLAP evolution equation in Q 2 , given one boundary condition. This is done in the presence of the regulator, i.e. in D = 4 − 2ǫ dimensions. It is useful to choose for boundary condition the normalization [4] f (N, Q 2 = 0, ǫ) = 1, which yields the solutioñ Here Γ qq is the Mellin transformation of the (singular) quark splitting function P qq (x, ǫ, α s (µ 2 )).
It is important to keep in mind that the expression in Eq. (34) is defined only in the presence of the regulator ǫ. The integration over the coupling constant down to µ 2 = 0 is also defined by using a dimensionally continued α s [4]. When we attempt to compute the coefficient functionC L from the ratio in Eq. (32) we encounter the following two problems. First, we cannot use the exponentiated form of Eq. (30) because Sudakov factorization organizes the exponent ofF L in terms of logarithmic in N dependence whereas the DGLAP evolution organizes the exponent off in terms of logarithms of Q 2 . In other words the ratioF L /f is not a simple ratio of two exponentials because parts of the prefactors in the Sudakov exponentiation ofF L are contained in the DGLAP exponent off . The consistent way to compute the coefficient function is to write alsoF L as a single exponential by solving its own DGLAP equation in D = 4 − 2ǫ dimensions. This is indeed what is done forF 2 . The boundary condition used in theF 2 DGLAP equation is the normalizatioñ and the solution is again a single exponential, like in Eq. (34). Then the coefficient functionC 2 follows simply from Eq. (32) with L → 2. It is here that we encounter the second problem for the longitudinal strucutre function.F L cannot be normalized as F 2 above, because it starts at O(α s ) 1 .
Although we cannot treat formallyF L likeF 2 , we observe that a large part in their respective Sudakov exponentiated forms is the same. Indeed, if we repeat the previous renormalization group analysis forF 2 , we get an answer forF 2 which can be obtained from Eq. (30) by simply substituting F L , H L → F 2 , H 2 and J → J ′ .F 2 can be normalized as in Eq. (35) because J = 1 + O(α s ). Therefore, if we know the coefficient functionC 2 to some order in α s we can predict the coefficient functionC L to one order higher by using the Sudakov factorization. All we need to do is to isolate the parts that are different between F 2 and F L in their exponentiated forms. The problem of normalization ofF L at Q 2 = 0 is solved by defining a functionF L as .
Since the jet has been moded out, we can normalize the ratio asF L (N, Q 2 = 0, ǫ) = 1. Then we compare the Sudakov exponentiation ofF L with the corresponding one forF 2 and we identify the parts in the exponent that are different. These parts are the corrections that must be included inC 2 to obtainC L . Once we do this we restoreF L via Eq. (36). The net result of this manipulation is that the coefficient function forC L is predicted to bẽ and the ratio of the two structure functions can be readily obtained from their Sudakov exponentiated forms. The net effect of the ratioF L /F 2 is to remove fromC 2 the contributions coming from the jet J of F 2 and replace them with the contributions from the J ′ jet of F L . To O(α 2 s ) the above relation yields

Explicit calculations to O(α 2 s )
In this section we prove via explicit calculations the F L factorization formula Eq. (6) to O(α 2 s ), or equivalently that Eq. (38) holds. This is done by comparing our results with the coefficient functionsC (n) L (x, 1), n = 1, 2 in Eqs. (9,12). The new object to be computed is the jet function J ′ defined in Eq. (7), since both the soft function V and the quark distribution φ are the same as in F 2 and can be taken from the literature. We show that the singular ln(1−x) fragmentation terms inC (2) L (x, 1) do arise from the UV renormalization of the operator defining J ′ and we construct its anomalous dimension γ ′ J ′ to O(α s ). The analysis of subsection 3.2 resulted in a set of Feynman rules for computing the logarithmic enhancements of the new jet function J ′ . Specifically, the coupling of gluons from thev-jet to the external lines is given by the effective non-local vertex O λ (l) defined in Eq. (10). In that expression, though, the transverse and light cone polarizations have been separated, whereas we find it more convenient for our calculations to compute traces with covariant objects. To this end we define the covariant vertex O λ (l) cov as Then tr pO λ (l)... = tr pO λ (l) cov ...
so we can use the covariant vertex for computing traces. Note that in terms of the operator definitions, O cov is generated by the operator in the numerator of Eq. (7) after performing a gauge transformation in the light cone direction. All functions entering the factorization formula have been defined in a gauge independent way and from now on we employ the Feynman gauge. Finally, in the MS scheme it is convenient for presenting intermediate results to define the ratio of scales

The O(α s ) contribution to J ′
To first order in the coupling F L has no logarithmic singularities, the answer being simply The F L factorization formula to this order reads We check that J ′(1) obtained from Eq. (7) comes with the correct normalization to reproduce the above lowest order result. There is only one diagram to compute, that of fig. 4. From Eqs. (2) and (6) we have We note that the whole contribution to F L comes from J ′(1) and no IR subtractions are needed to this order. For convenience we have included in J ′(1) all overall normalizing factors associated with F L . This short calculation shows that J ′ as defined by Eq. (7) is correctly normalized.

The O(α 2 s ) contribution to the jet function J ′
The non-trivial infrared behavior of F L starts at O(α 2 s ). We compute the jet function to this order, J ′(2) , and show that after summation over final state cuts and IR subtractions as required by the denominator of Eq. (7) the only surviving singularities for the operator diagrams are simple UV poles. In section 3.3 we argued that the J ′ operator can be multiplicatively renormalized in moment space as In a massless scheme such as MS, the anomalous dimension γ ′ J ′ can be constructed from the renormalization constant Z J ′ as Equivalently, to the order we compute, and because of the relations we can directly generate the O(ln(1 − x)) terms for J ′ in x-space. Their coefficient is the sign-inverted residue of the UV singular part of J ′ (2) . There are seven diagrams that determine J ′(2) and they are depicted in fig. 5. Since certain subtleties occur in the calculation arising from the presence of both linear (eikonal) and quadratic denominators in the loop integrations we proceed by analyzing each diagram in turn. In the full calculation to this order there are fourteen diagrams for F L . Whenever a diagram in our effective approach is free from IR divergences we provide direct comparison with the corresponding one of the full calculation taken from ref. [7]. At the practical level, the effective approach is vastly more economical in capturing the logarithmic enhancements.
Diagram (a): where a factor of 2 has been included from left↔right equal contributions. The numerator factor is and the loop integrals yields through the use of Eq. (2) in the appendix Then the diagram becomes By including the overall factor from Eq. (2) we obtain the contribution to the unrenormalized This leads to the anomalous dimension as well as to the contribution for the renormalized jet function Diagram (b): and the numerator factor is In the second step only terms that lead to UV singularities have been retained. There is a subtlety involved in the derivation of the above expression. A term proportional to l σ l ρ is to be replaced as l σ l ρ → (1/2)l 2 η σρ . Note the factor 1/2 instead of the usual 1/4. The latter case is valid for the UV part of a loop integral if the integrand contains only quadratic denominators. In our case, however, it also contains a linear (eikonal) denominator and one must be careful not to overcount the polarizations that contribute to the UV divergent part (2 and not 4). Then, through the use of Eq. (1) in the appendix we obtain Similarly to diagram (a), the corresponding contribution to the anomalous dimension is and to the renormalized jet function Eq. (22) reproduces the result for graph (10) in [7].
Diagram (c): and the numerator factor can be approximated in the UV as The same remark about counting the polarizations that contribute to the UV singularities as in diagram (b) applies here. Then, through the use of Eq. (1) of the appendix we obtain which contributes to the anomalous dimension as and to the renormalized jet function as Eq. (27) reproduces the result for graph (13) in [7]. So far, the three graphs we have computed are in one-to-one correspondence with the corresponding graphs of the full theory. They are finite in the full theory and become UV divergent in the effective theory. Now we proceed to calculate the rest of the graphs. These are not in one-to-one correspondence with graphs of the full theory, which is the default case in factorization. No more comparison for intermediate results will be possible, only for the final answer.
Diagrams (d) and (e): For these two diagrams no subtleties are involved. Standard manipulations yield and The last two diagrams (f ) and (g) are the most interesting to this order because they have both UV and IR singularities. Recall that the jet function J ′ is by definition free from IR singularities and its UV renormalization generates the fragmentation ln(1 − x) terms in the coefficient function. We must therefore show that the left over IR singularities in the graphs (f ) and (g), after summing over their two and three particle cuts, are precisely removed by the soft function V . The remaining UV poles will contribute to the J ′ anomalous dimension.
Diagram (f ): The two particle cut of this diagram (f 2 ) is We have separated UV from IR singularities by retaining the sign in front of the regulator.
The three particle cut (f 3 ) can similarly be obtained through the use of the three massless particle phase space given in Eq. (8) of the appendix in terms of convenient coordinates for our case. The result is Adding the two and three particle cuts we obtain for the full cut diagram We note the cancellation of the double IR poles which arise from soft gluon emission. Simple IR and UV poles and their combination do survive.
Diagram (g): Similar manipulations as in diagram (f ) yield the following results for the two and three particle cuts.
Adding the two and three particle cuts we obtain for the full cut diagram As in diagram (f ), we note the cancellation of the IR double poles. UV and IR simple poles do remain. Finally, we add the contributions from diagrams (f ) and (g) and multiply the result by (2πx 2 /α 2 s Q 2 ), see Eqs. (1,2), to obtain their contribution to the unrenormalized J ′ (2) . It is convenient to separate the color structures, since they behave quite differently. The result is We observe the absence of any IR poles in the C F C A color structure.
Total contribution to J ′(2) : To O(α 2 s ) the jet function J ′(2) is obtained from the sum of the diagrams (a)-(g) after subtracting the purely soft contributions to the same order. The IR subtraction term for F L starts at O(α 2 s ) and it is J ′(1) ⊗ V (1) . Explicitly, where the subscript u denotes that the jet function is unrenormalized and IR unsubtracted. The soft function V (1 − x) is universal and familiar. To O(α s ) it is defined in fig. 6. In the Feynman gauge it is The argument of V is the longitudinal momentum fraction that enters thev-jet. The first term in the square brackets is of UV origin and comes from the virtual contribution. Note the absence of IR double poles. We compute the convolution J ′(1) ⊗V (1) using the integral lemmata (9)(10)(11)(12)(13) in the appendix. The result is As only diagrams (f ) and (g) have non-trivial IR singularity structure, we consider only the (f ) + (g) part of J ′(2) and perform the IR subtraction there. Specifically, All IR singularities have been removed from J ′(2) as dictated by its definition in Eq. (7). The left over terms are simple UV poles. Now we are finally able to compute the contribution of the diagrams (f ) and (g) to the jet anomalous dimension and to the UV renormalized jet function in x space and J ′(2) It is true that we could have arrived at the above equation by just computing the UV simple poles of only the two particle cut graphs and converting them into ln(1 − x) contributions. We, nevertheless, took the pain of calculating also the three particle cuts and the subtraction terms to show that indeed only simple UV poles survive in the final expression, thus showing that an anomalous dimension can be defined to sum the singular terms.
The leading soft double logarithm is generated by the UV renormalized J ′(1) ⊗ V (1) , Eq. (41). It is We have organized the summation of the singular ln(1 − x) terms in such a way that all fragmentation O(ln(1 − x)) terms come from J ′(2) and the soft O(ln 2 (1 − x)) term comes from J ′(1) ⊗ V (1) . At this stage we can determine the coefficientC (2) L (x, 1) as given by the F L factorization formula. It isC (2) L (x, 1) = (J ′(1) ⊗ V (1) )(x) + J ′(2) (x), which, through Eqs. (45, 46), exactly reproduces the singular part of the result from the full calculation in Eq. (12). We have thus shown explicitly how the factorization formalism captures the logarithmic enhancements to O(α 2 s ). Finally, the jet anomalous dimension γ ′ J ′ is and the anomalous dimension γ J ′ that does not include the running of the coupling is obtained from Eq. (31). It is

Summary and conclusions
We have studied in detail the longitudinal structure function in deeply inelastic scattering near the elastic limit x → 1 and have shown that it can be Sudakov factorized in terms of two jet and one soft functions. The soft function and the target jet function are the same as in F 2 , whereas the current jet function is characteristic of F L and measures the transverse momentum of the struck parton. We have shown that the quasi-elastic logarithmic corrections from the fragmentation of the current jet can be resummed in moment space in terms of the anomalous dimension γ J ′ , which we calculated to O(α s ). We emphasize that our analysis is at the leading twist level but goes beyond the light-cone expansion. It resums terms that are of order ln k N/N, k > 0. This approach can be extended to the case of F 2 , where such terms are present, but subleading relative to the order ln k N, k ≥ 0 terms. It is quite remarkable that terms that are power suppressed in N are amenable to a similar systematic analysis as the leading ones, just through the introduction of a single new function, the jet J ′ . This leads us to the conjecture that, indeed, all power in N suppressed terms can be resummed in quasi-elastic processes, via the introduction of further jet functions defined by operators with higher powers of the transverse derivative. In Ref. [5] a systematic procedure has been presented for contructing these higher dimensional jet functions, and the present analysis is really a concrete example of that procedure.
Returning to F L , we note also that the logarithmic enhancements from the final state fragmentation are delayed by one order of α s relative to the ones in F 2 . As these enhancements are summed to all orders, part of their effect (the beta function term in Eq. (30)) is to change the overall normalization by shifting the scale of α s in the preexponent factor from Q 2 to Q 2 /N. Thus, in the region x → 1 higher order corrections are parametrically as important as lower order ones. This has implications for obtaining the power corrections to F L , using for instance renormalon methods, since the power corrections come predominantly from this near elastic region. This enhancement, which is characteristic of F L , can be thought of as an indication from perturbation theory of the large power corrections that this observable is known to receive [15]. These and applications of the formalism to other inclusive processes are currently under investigation.