Large Transverse Momentum in Semi-Inclusive Deeply Inelastic Scattering Beyond Lowest Order

Motivated by recently observed tension between $O\left(\alpha_s^2\right)$ calculations of very large transverse momentum dependence in both semi-inclusive deep inelastic scattering and Drell-Yan scattering, we repeat the details of the calculation through $O\left(\alpha_s^2\right)$ transversely differential cross section. The results confirm earlier calculations, and provide further support to the observation that tension exists with current parton distribution and fragmentation functions.


I. INTRODUCTION
In a previous article [1], we discussed the semi-inclusive deep inelastic scattering (SIDIS) process: l(l) + Proton(P ) → l (l ) + Hadron(P H ) + X, and we highlighted the challenge of finding agreement between O (α 2 s ) calculations and existing SIDIS data in the very large transverse momentum transverse momentum limit where standard collinear factorization is expected to be valid.One motivation is that obtaining a description of the small transverse momentum behavior associated with nucleon structure requires a good understanding of the matching to large transverse momentum where a transverse momentum dependent (TMD) factorization description fails.Given the current focus on using deeply inelastic hadro-production to access nucleon structure sensitivity sensitivity, it is imperative to examine theoretical framework for the full q T -range in greater detail.In the language of Ref. [1], we are interested in this paper in what was there called "region 3" behavior, corresponding to where q T is so large that the small q T /Q approximations associated with TMDs are not reliable, but where ordinary collinear factorization should be applicable and reliable.Calculations to O (α 2 s ) have existed for some time [2,3].The main observation of Ref. [1] was that, while the O (α 2 s ) correction gives an order of magnitude increase over the leading order, that still is not sufficient to achieve reasonable agreement with data for Q in the region of one to several GeVs, transverse momentum of order Q, and for moderate Bjorken-x.Combined with similar observations concerning the Drell-Yan pointed out in [4], this points to general tension between transverse momentum dependent cross sections at q T ∼ Q and collinear factorization.
There are a number of potential explanations or solutions, including a direct re-tuning of collinear parton distribution and/or fragmentation functions to transversely differential SIDIS cross sections.(A proposal to constrain gluon PDFs in transversely differential Drell-Yan cross sections was made already in [5].)But before proceeding to consider these directions, it is important to validate the O(α 2 S ) in [2,3] that led to our conclusions in [1].Therefore, we have in this paper repeated the large transverse momentum order O (α 2 s ) calculation following a slightly different formal framework.We are able to reproduce the results in [2] very closely, thus bolstering the observations we made earlier in [1].
The specific purposes of this article are as follows: i.) to lay out the logical steps of the calculation with enough detail, we hope, to lead to clues as to how to improve phenomenological agreement, ii.) to make available a convenient numerical implementation for the large q T region and, iii.) to present some quantitative results relevant to current experimental programs such as those at COMPASS and Jefferson Lab 12 GeV.Overall, our results support the general observations made in [1].
In our calculations we use qgraf [6] to generate Feynman graphs and FORM [7] to carry out spinor and color traces for the amplitudes.The renormalization counter terms are computed in MATHEMATICA with the packages FeynArts [8,9] and FeynCalc [10,11].
Further analytic manipulations were performed in MATHEMATICA.The code to compute the cross sections are publically available at [12].
In Sec.II, we explain our setup and notation.In Sec.III we summarize the organization of our calculations.This includes a classification of the partonic Feynman graphs needed at order α s (Q) 2 and a discussion of the phase space integrals for multiparton final states.We discuss the results of the calculation in Sec.IV, and give concluding remarks in Sec.V.

II. NOTATION AND CONVENTIONS
We will express quantities in terms of the conventional kinematical variable z ≡ P H •P/(P • q).P H,T is the Breit frame transverse momentum of the produced hadron, and P and q are the four-momenta of the incoming target hadron and the virtual photon respectively.We will focus on the unpolarized and azimuthally-independent cross section since this is the most straight forward observable to calculate.(Although the results have potential implications also for polarization dependent observables.)Many general treatments of SIDIS as a process are available [13][14][15][16][17][18][19][20][21][22][23][24].Our notation is consistent with trends among these, with modifications as needed for our current purposes.
The unpolarized differential cross section is where σ H is the unpolarized hadronic cross section.S lP = (l + P ) 2 .The hadron transverse momentum P H,T is defined in a frame where the photon and incoming hadron are backto-back (a "photon" frame).The Bjorken x and the y variable are the usual definitions x ≡ Q 2 /2P • q and y ≡ P • q/P • l.L µν is the usual leptonic tensor and W µν is the hadronic tensor for SIDIS: where we have omitted spin and azimuthal angle dependent terms since we do not consider these in this paper.Also, we assume that x and 1/Q are small enough that both the proton and lepton mass can be dropped in kinematical and phase space factors.The normalization convention in Eq. ( 4) is so that the prefactor on the right hand side of Eq. ( 2) is similar to the unpolarized case.
The unpolarized structure functions F 1 and F 2 are defined by the usual gauge invariant decomposition W µν = −g µν + q µ q ν q 2 F 1 + (P µ − q µ P • q/q 2 )(P ν − q ν P • q/q 2 ) P • q F 2 . (5) Then, the cross section is In calculations, it is convenient to work with Lorentz invariant structure function extraction tensors P µν Γ with Γ ∈ {g, P P } where P µν g = g µν , P µν P P = P µ P ν . Then, where in 4 − 2 dimensions, (We we always use the massless target approximation -see discussion in [25].)It is useful to express transverse momentum in terms of In a frame where the incoming and outgoing hadrons are back-to-back, q T is the transverse momentum of the virtual photon (assuming, as always in this paper, that external hadron masses are negligible).
Our treatment of factorization will follow the general style of [26].The factorization theorem that relates the hadronic and partonic differential cross sections in SIDIS is The 1/ξ is from the partonic flux factor, and the 1/ζ 2 is from the conversion between k 1 and P H .The indices i and j denote, respectively, the flavors of the parton in the proton (with a momentum fraction ξ) and of the outgoing parton that fragments into hadron H, whose momentum is a fraction ζ of parton j momentum.The incoming and outgoing parton momenta p and k 1 satisfy p = ξP and k 1 = P H /ζ. (Indices i and j for incoming and outgoing partons p i and k 1,j are not shown explicitly but are understood).f i/P (ξ) and d H/j (ζ) are the collinear parton distribution and fragmentation functions respectively.It is also useful The differential partonic hard part in Eq. ( 11) is finite and well-behaved, and at large transverse momentum it starts at O (α s ).
The unpolarized partonic structure tensor for scattering off parton i into parton j is defined in exact analogy with the hadronic tensor: with, Thus, from Eq. ( 11), The partonic structure function decomposition is Then, III. ORGANIZATION

A. Basic Setup
For a process with N final state partons, P µν g W µν and P µν P P W µν are squared amplitudes integrated over the N -particle phase space for the outgoing partons, Here |M g | 2 or |M P P | 2 is the squared amplitude for the process γ * (q) + parton(p) → parton(k 1 ) + (N − 1) spectator partons, with the polarization sum of the virtual photon replaced by g µν or P µ P ν .The subtraction terms in Eq. ( 19) are needed to remove double counting with lower orders of perturbation theory and cancel singularities in the first term.The form of the subtraction term will be discussed in Sec.III B. It is assumed in Eq. ( 19) that all integrals that allow kinematical δ-functions to be evaluated have been performed.Also, the phase space factors associated with k 1 are excluded from the partonic phase space since they give the ẑ and k 2 1T dependence of the differential partonic cross section.The 1/(2π) 4 on the right hand side of Eq. ( 19) is the same factor in Eq. ( 4).The dΠ (N ) represents a generic phase space factor for 2 → N scattering (see Eq. (35) and Eq. ( 39)).
Thus it is convenient to express the structure functions in the form where the N is a sum over all possible partonic final states.
We will express the phase space in terms of Mandelstam variables: where i, j = 1, 2, 3, with the labeling in Fig. 1.For simplicity, u 1 and t 1 will be abbreviated as u and t from here on.All Mandelstam variables will refer to the partonic cross sections.
Occasionally it will be useful to change kinematical variables, For example, The O(α 0 s ) contribution to Eq. ( 19) is kinematically constrained to k 1T = 0.At O(α 1 s ), only tree level processes contribute, and no singularities appear.At O(α 2 s ), soft, collinear and UV singularities arise as 1/ and 1/ 2 poles in dimensional regularization with space-time dimension n ≡ 4 − 2 .
These singularities cancel in the sum of real, virtual and counterterm graphs and after applying collinear factorization.More discussion of the singularity structure at O(α 2 S ) is available in [2,27].The O(α s ) partonic cross section is described in details in [17], and our results match that calculation.At O(α 2 s ), details below can also be found in, e.g., [27][28][29][30][31].The change of variables in Eq. ( 20) from momentum fractions to Mandelstam variables is where ζ is replaced by s 23 , the virtuality of the spectator parton system.(s 23 is k 2 X in the notation of [1]).Overall momentum conservation s + u The graphical structures needed at order α s (Q) 2 can be classified as in Fig. 2. Each diagram corresponds to a contribution to the squared amplitude, and the blobs include all possible attachments.For the 2 → 2 subprocesses, the blobs include up to one loop corrections in one side of the cut to give O(α 2 S ) accuracy to the hard parts.The black dot marks the fragmenting parton and the other lines are integrated over all phase space.
The following are assumed and not shown explicitly in the graphs: • Subtractions, consistent with factorization, for all collinear divergences.
• One QCD loop at all possible positions inside the blob in virtual processes.
• Associated UV counter term graphs for virtual processes.
Define the hard parts H for individual graphs as where Γ ∈ {g, pp} and i, j, k, l ∈ {q, q, g}.As before, i and j label incoming and fragmenting parton flavors respectively while k and l label the unobserved parton flavors.The graph is virtual when only one flavor index appears after the ";".Figures 2(V-A) through 2(V-F) correspond to the O (α s ) graphs if all virtual loops in the blobs are removed.
So, for example, H g qq;gg is represented by Fig. 2(R-A), when the partonic tensor is contracted with g µν .Fig. 2(R-E) includes both H g qq;q q and H g qq ;q q contributions where the prime indicates a different flavor.Note that graphs like Fig. 2(R-G) give contributions like H g q q;qq or H g q q ;q q .We work in the approximation that all particles are massless, except for the photon which is highly virtual, q 2 = −Q 2 .In Eq. ( 19), for 2 → 2 scattering, 2. Structure of graphs needed at order α s (Q) 2 .The last six correspond to virtual corrections to α s (Q) order graphs.
Thus the 2 → 2 phase space factor is For the 2 → 3 case, we follow the strategy in [29,30] where we work in the rest frame of the k 2 + k 3 system, in which where kn−3 denotes unit vectors for the first n − 3 components of the n − 1 dimensional unit spatial vector in spherical coordinates.See Fig. 3.
In the center-of-mass of p and q, the spatial vectors p, q, k 1 , and k 2 + k 3 are in the same plane.We then boost to the k 2 + k 3 rest frame p, where q, and k 1 are still in the same plane -see Fig. 3.We choose the spatial orientations of this coordinate system so that p, q, and k 1 are in the plane created by the last two spatial axes.For these vectors, the first n − 3 spatial components are zero.Then the scattering amplitudes do not depend on the first n − 3 spatial components of k 2 or k 3 , and the 3-body phase space simplifies (see Appendix A for useful identities relating to the 3-body final state phase space).In this frame, where in the third equality we use that the scattering amplitudes are independent of the 3. Schematic of SIDIS with 2 → 3 partonic scattering, using the orientation of the vectors defined in Eqs. ( 36)-( 37) in 4-dimensions.A 2 → 2 scattering configuration is recovered if (k 2 +k 3 ) 2 becomes much less than Q 2 , or if k 2 or k 3 align with one of the other hadrons.The three frames defined in this section are obtain by rotating p, q and k 1 inside the plane; frame 1 corresponds to orienting the z-axis along k 1 (Eq.(B4)), frame 2 orients the z-axis along p (Eq. (B5)), frame 3 orients the z-axis along q (Eq.(B9)).

From overall momentum conservation
Since s 23 can become zero in certain regions of x and ẑ, dimensional regularization is needed for the s 23 → 0 behavior.Further details of the 2 → 3 phase space are provided in the next subsection.
To summarize the basic steps, integrals for the unobserved final state phase space are first performed over angular variables.Then, integrals over momentum fractions are expressed via Eq.( 29) as integrals over ξ and s 23 .In the s 23 → 0 limit there are soft and collinear singularities, proportional to δ(s 23 )/ or δ(s 23 )/ Once all real and virtual Feynman graphs have been computed in dimensional regularization, they must be combined in an algorithm consistent with a factorization theorem.
According to the collinear factorization theorem, hard factors are independent of the species of target and final state particles, so hard scattering calculations can be performed for free and massless partons as the target and final state without loss of generality.(In other words, we will consider perturbatively calculable parton-in-parton PDFs and and parton-toparton FFs.)This simplifies calculations, so we start our calculations by rewriting Eq. ( 15) as a factorization theorem for a partonic initial and final state: All manipulations will be done in n = 4 − 2 dimensions until the very end.The left side is now an argument of p , q, and k 1 , and the PDFs and FFs have subscripts i/i and j /j.The momenta p and k 1 play the roles that P and P H,T played earlier.The initial p has a flavor i and the final k 1 has a flavor j .The second line in Eq. ( 41) defines the "⊗" notation as the usual shorthand for the convolution integrals in ξ and ζ.On the last line, W µν,i j and W (NLO) µν,i j are hadronic tensors calculated to LO and NLO in the hard parts, i.e. what one normally means when one speaks of a "leading order" or "next-to-leading order" calculation in collinear pQCD.
To get Ŵ (NLO) µν,ij (q, x/ξ, z/ζ), first write The term in braces is the correction to the leading order so it vanishes by construction at order α s in the hard part.Also, it contains subtractions for the overlap with the LO, so it is infrared safe through order α 2 s in hard scattering.Thus, we may rewrite it as: µν,i j (q, x/ξ, z/ζ) The term in braces on the second line is equal to the first line evaluated with hatted partonic variables and expanded only to order α 2 s in the coupling.Comparing with Eq. ( 42) and Eq.(44) shows that the factor in braces in Eq. ( 45) is Ŵ (NLO) µν,ij (q, x/ξ, z/ζ).We use Eq. ( 43) to evaluate it explicitly in terms of low order PDFs and FFs: Here, the W (NLO) µν (q, x, ẑ) unsub is the sum of order O(α 2 s ) graphs contributing to W µν (q, x, ẑ) without any subtractions (i.e., graphs in Fig. 2(R-A) through Fig. 2(V-F) integrated over phase space, but absent the subtractions).The second two terms are obtained from Eq. ( 43), expanded to order α 2 s .The last line in Eq. ( 46) uses the known results for PDFs and FFs calculated to O (α s ) for massless quarks and gluons: Note that, consistent with the MS scheme, the renormalized order α s PDFs and FFs in dimensional regularization for massless pQCD are just S / poles.There are no µ factors or logarithms of µ.Thus the last line of Eq. ( 46) involves no factors of µ .For completeness, the one-loop splitting functions [32] are Our calculations need Eq. ( 46), contracted with specific P µν Γ tensors.The substraction scheme has been carry out without contracting the Lorentz indices of the partonic tensor.
However the method is equally applicable for any kind of contraction with external momenta.
In our calculations we carry out the substractions separately for the two extraction tensors in Eq.( 7) and verified analytically the cancellation of all infrared and collinear singularities.

A. Combining real and virtual contributions
After the hard real and virtual contributions are calculated, they must be combined into infrared safe squared amplitudes.Table I shows graphs from corresponding real and virtual processes.Note that the correspondence between real and virtual processes is not one-toone.In particular, infrared (soft and collinear) singularities in H qq;g are canceled by three real processes.A subtlety arises when a quark loop appears on the gluon leg of H qq;g .This creates a collinear pole term proportional to N f , the number of massless quark flavors in the loop.This pole term is canceled after adding the real processes H qq;q q and H qq;q q with all massless flavors q (other than q) [29].For processes with a non-spectator gluon leg, the N f dependence of the collinear pole is removed by factorization.Also notice that some real processes have no corresponding virtual ones, and in these cases factorization subtractions are sufficient to remove all infrared poles.Since many graphs need to be combined, the results are presented as the six scattering hard parts listed in the last column of Table I.
We also need the hard parts for processes with reversed quark number flow in one or two open fermion lines in the graph.That is, in one or two open quark lines, quarks and their corresponding anti-quarks are interchanged.These are easily related to the hard parts already listed in Tabel I by using the fact that QED vertices acquire a minus sign under charge conjugation.The results are summerized in Tabel II.Note that when a quark line links an outgoing spectator quark anti-quark pair, there is no need for the reversed quark flow.This is because the spectator momenta are integrated over and interchanging the quark and anti-quark will double count the contribution.

B. Comparison with existing results
After combining the graphs in each of the six scattering channels in Table I, we have verified explicitly that all single and double poles cancel.The terms left are finite in the limit → 0, and constitute the infrared safe hard parts in the last column of Table I.As a check, we have compared with a (privately obtained) computer calculation that appears to reproduce the results of Ref. [2] but is modified to be consistent with the kinematics of the + H qq;q q H qq H qg;q H qg;qg H qg N/A H gg;q q H gg N/A H q q;qq H q q N/A H qq ;q q H qq TABLE I. Correspondence between real and virtual graphs at order α s (Q) 2 .The superscript of the hard part Γ ∈ {g, pp} is supressed.For virtual graphs, the additional gluon is inside one of the blobs.
current experimental data.For most kinematics, the calculations agree within experimental uncertainties, but we found several possible differences.
• In the last row of Table I, our calculation of terms from the P µν g projection with the charge structure e 2 q and e q e q differs from the corresponding terms calculated in the previously existing code by a minus sign.(By contrast, terms with e 2 q agree.) • As is discussed in subsection IV A, adding the real processes H qq;q q and H qq;q q with all massless flavors q gives terms proportional to N f .The pole parts of these terms are canceled by the corresponding poles from H qq;g .However, a finite part remains since the real processes above with N f final state quark pairs give identical contribution even when the quark pairs are not collinear.These terms contribute to the third row of Table I.In the previously existing code, we did not find an explicit dependence on N f in this channel.
The numerical result of our calculation can in some cases differ from the previously existing code by as much as ∼ 100% for certain individual subchannels.In the kinematics of the COMPASS experiment, this translates into a discrepancy of up to ∼ 20% for the overall cross section.

C. Phenomenological Results
We examine the impact of the O(α 2 S ) corrections in SIDIS by plotting the NLO to LO ratio K-factor for the F 1 structure function.We will consider values of z between 0.2 and 0.8 since this is a region where factorization theorems based on current fragmentation in SIDIS is conventionally expected to apply.For x, Q, q T , we choose values of Q = 2, 20 GeV, q T = Q, 2Q and 0.001 ≤ x ≤ 1, which correspond to kinematics ranges accessible by existing experimental facilities such as COMPASS and JLab 12. Fig. 4 shows the K-factor across the aforementioned kinematical range, and it shows a clear concave up shape in its dependence on x with values that decrease as Q varies from Q = 2 to 20 GeV.
It is notable that even for Q as large as 2 GeV, the K-factor is greater than 2, even at its minimum value.Note also that this minimum is approximately in the region of 0.01 ≤ x ≤ 0.1, close to the valence region relevant to many hadron structure studies.
The K-factor increases at smaller values of x, indicating that the NLO corrections becomes increasingly important for describing regions with large (P + q) 2 .It also increases at large x until it reaches the kinematic boundary where the phase space for the hadron production vanishes.
The enhancements at large x can be traced to the logarithms of B (see Eq. 32 and appendix D) which are associated with soft gluon effects near the kinematical threshold.
Thus it is likely that threshold resummation becomes important in this region.Point-bypoint in q T , effects near the kinematical boundaries are more apparent than in q T -inclusive cross sections.It is interesting to compare with observations found in [33] where threshold resummation effects for q T integrated SIDIS were found to give sizable corrections, but not the factor of 5 enhancement needed to explain the largest q T present in Fig. 4.
The minimum in the K-factor is helpful for identifying regions in x where ordinary fixed order treatments are most likely to be sufficient.Given a set of values for z, q T , and Q, let x 0 denote the x value corresponding to the minimum K-factor.Then bins in x, Q, z ,q T can be classified according to x < x 0 and x > x 0 .We show this in Fig. 5 for kinematical bins corresponding to recent COMPASS q T -dependent SIDIS multiplicities h ± [34] for q T /Q > 1.
The values of x the for which the K-factor reaches its minimum for q T ≈ Q can be read from the plots.Away from these regions, additional resummation techniques may be needed.
It may also be necessary to update collinear PDFs outside these regions in order to fully describe the large q T behavior.As mentioned in [1], the large q T tails of SIDIS are sensitive to the large ξ PDFs and large ζ FFs.These could potentially open up new opportunities to constrain collinear PDFs at large momentum fractions.
At the same time, factorization theorems for the full q T -dependent SIDIS spectrum, including q T ≈ 0 with non-perturbative TMD PDFs and TMD FFs, need the large q T ∼ Q component in order to have a completely reliable tests of the factorization formalism.
Ideally, such validation will take place when q T -dependent SIDIS data are included in the simultaneous extraction of collinear PDFs and FFs as well as non-perturbative TMD PDFs and FFs in QCD global analysis.

V. CONCLUSION
The computational tools necessary to reproduce Fig. 4, Fig. 5 and other similar calculations are available at [12], along with documentation.As indicated in the introduction, these generally confirm earlier calculations (e.g., [2,3]) of the overall cross section to within about 20%, well within present experimental uncertainties, and thus strengthens our earlier position [1] that significant tension exists between existing SIDIS data and collinear factorization calculations at large q T Q.
We have identified region where collinear factorization appears most reliable as the regions with minimal K-factors, with current sets of collinear PDFs and FFs.As discussed in Sec.IV C, minimal values for the K-factor lie in the region 0.01 < x < 0.1, that is, ap-virtual real H gq;q H g q;q = H gq;q H gq;qg H g q;qg = H gq;qg H qq;g H q q;g = H qq;g H qq;gg H q q;gg = H qq;gg H qq;q q H q q;q q = H qq;q q H qq;q q H e 2 q q q;q q = H e 2 q qq;q q ; H e 2 q q q;q q = H e 2 q qq;q q H qg;q H qg;q = H qg;q H qg;qg H qg;qg = H qg;qg N/A H gg;q q N/A N/A H q q;qq H qq;q q = H q q;qq N/A H qq ;q q H q q ;qq = H qq ;q q H e 2 q qq ;q q = H e 2 q qq ;q q ; H e 2 q qq ;q q = H e 2 q qq ;q q ; H eqe q qq ;q q = −H eqe q qq ;q q H e 2 q q q ;qq = H e 2 q qq ;q q ; H e 2 q q q ;qq = H e 2 q qq ;q q ; H eqe q q q ;qq = −H eqe q qq ;q q TABLE II.The first column of real and virtual are the hard parts as in Table I.The second column gives the hard parts needed other than those in the first column.The superscript of the hard part Γ ∈ {g, pp} is supressed.e 2 q , e 2 q , and e q e q are various combinations of quark electric charges.A charge combination as a superscript denotes terms proportional to that particular charge structure in a hard part.
4. K-factor ratio.We use CJ15nlo PDFs [35] and DSS09 FFs [36].The ratio is bounded at large x due to limiting phase space for hadron production.The small noise in the lower left panel is due to the oscillatory behavior in the interpolation of FFs tables.x > x 0 x ≤ x 0 FIG. 5. SIDIS kinematics of h ± at COMPASS [34].x 0 is the location in x where the minumum of K-factor is attained at a given (Q, q T , z).Samples with x > x 0 indicates sensitivity to large threshold corrections.
proximately in valence kinematics.The increase of the K-factor at both smaller and larger values of x, may point at the importance of resummation effects outside the valence region, or the need to refit collinear PDFs and FFs.This last interpretation is consistent with that of our previous study [1].It is also important to note that the 2 → 2 kinematics of the O (α s ) contribution places severe kinematical constraints on the relationship between the initial and final state partons.Therefore, it is likely that the generally large K-factors are at least partly due simply to a kinematical suppression of the O (α s ) contribution, and are not a fundamental problem with the convergence of the perturbation series.We note that Reference [37] has addressed somewhat similar issues, but in the region where q T is still small enough that a q T /Q power expansion is still meaningful, and also in the limit of Q 2 s.
As a next step, we plan to refit collinear functions in the large q T region, using SIDIS and e + e − -annihilation to back-to-back, and Drell-Yan scattering to explore the possibility that this allows the region of minimum K-factor in Fig. 4 to be fully accommodated with minimal modification to existing fits.
Frame 3: q = (q 0 , 0 n−3 , 0, q ), (B9) where sin sin sin There are three cases (see Appendix C for a proof) for the coefficients in the integral in Eq. (B1): (1) Integrals in case (3) can be transformed to case (2) by switching to another frame.With multiple frames, the number of integrals to compute can be reduced.In the following we show how the integrals in case ( 1) and ( 2) are obtained.
Case (1): In this case the integral in Eq. (B1) has a closed form [71] π where for simplicity an overall factor is not shown and we choose a = A = −b = 1, B = cos χ, and C = sin χ.In the frames we choose above, χ can be α 1 or θ 2 .Here the signs of the trigonometric function terms in the denominator are chosen to be negative.All other choices can be transformed to this one by applying substitutions , and/or χ → π + χ.The derivation of this integral can be found in appendix A of [71].

Case (2):
In this case the integral no longer has a closed form.Following steps similar to appendix A of [71], we arrive at where D > 1.The integral representation of hypergeometric function is used.Note the I j,l notation for integrals with different denominator structures.The result of evaluating Eq. (B19) for specific j and l is given in Appendix F.Here χ can be α 3 or θ 3 .The remaining z integral has to be performed by expanding the integrand into series of , and computing the integral order by order.The subtlety comes from the treatment of the factor (1 − z) n/2−l−2 , which can produce poles for l ≥ 1.The integral has the form where f (z) is a regular function at z = 1.When l = 1, the pole term can be made explicit using the identity where the "plus" functions are defined, for a function p(ξ), through When l > 1, the integral in Eq. (B21) is divergent.Nonetheless, it can be analytically continued.To do this, we write Note that, in the first term on the right-hand side of Eq. (B24), the power of 1 − z is increased by 1, with a regular function f 1 (z).The second term is divergent, but can be analytically continued using Γ functions: The manipulation in Eq. (B24) can be repeated until the power of 1 − z in the integral becomes − − 1, in which case Eq. (B22) can be used.To expand the hypergeometric function in the integrand in Eq. (B19), we note that if either j or l is less than 1, the hypergeometric series terminates and the function reduces to a polynomial.For j = 1 and l = 1 we use the expansion For the case that one of j and l is larger than 1 and the other is at least 1, we use Gauss's contiguous relations to reduce j or l.For example, for F (1, 2, 1 − ) and Three ADMV relations (3ARs): Another type of relation is between three ADMVs, each from a different ADMV type.We can get such relations for any three ADMVs of different types.To do this notice that exchanging the positions of k 2 or k 3 with k 1 in Eq. (D4) and squaring will give one of such relation.For instance, squaring p + q − k 2 = k 1 + k 3 : a relations between s 13 , t 2 , and u 2 .After one such relation is obtained, the others are obtained by replacing the ADMVs in that relation using three 2ARs.
Now the partial fraction is implemented in two steps: Step 1: Consider a term T 1 in the squared amplitude.The first step is to reduce the number of ADMVs in the denominator to two or fewer, and to ensure that no two ADMVs of the same type appear in the denominator.This is done in two substeps: Substep 1: First, separate any ADMVs of the same type in the denominator.
If two ADMVs of the same type appear in the denominator of T 1 , use a 2AR for the two ADMVs to write a factor of unity in the form and write with c j = ±1.Expanding the right side of Eq. (D10) produces a sum of new terms.Each will have one power lower of these two ADMVs in its denominator.
Further multiplication by unity factors like Eq. (D9) can be repeated until in all terms only one of the two ADMVs appears in any denominator.As an example, consider t 2 and t 3 both appear in the denominator, so we will use Eq.(D7) to write Eq. (D9): (D12) Then Substep 2: After Substep 1 is done for all three types, the leftover terms have denominators with at most three ADMVs, all from different types.For terms with two or fewer ADMVs in the denominator, nothing more needs to be done in this step.For terms with three ADMVs in the denominator, one may write a 3AR in the form of Eq. (D9) and reduce the number of ADMVs in the denominator to two or fewer in a way similar to Substep 1.For example, say that Then it is possible to eliminate s 13 , u 2 and t 2 by using Eq.(D8) to construct Eq. (D9).Then, Step 2: After Step 1 is completed, all terms have two or fewer ADMVs in the denominator.To get to the final form in Eq. (D3), we need to make sure no third ADMV appears in the numerator.This is done by noticing that the ADMVs in the numerator can be written in terms of the ADMVs of the denominator using 2ARs and 3ARs.
For instance, if one term has t 2 u 3 in the denominator, then any t-type ADMV in the numerator can be replaced by t 2 , and u-type by u 3 , using 2ARs.Any s-type ADMV can be replaced by a linear combination of t 2 and u 3 using the corresponding 3AR involving t 2 and u 3 .This puts T 1 in the form of Eq. (D3 At order O(α 2 s ), virtual corrections to the partonic cross section involve only one-loop integrals.Comprehensive reviews of the methods for this type of calculation can be found in [72][73][74].In this article, we use the traditional Passarino-Veltman (PV) approach [75], for which we closely follow the notation of [74].
Our O(α 2 s ) virtual contributions are the interference terms between tree-level 2 → 2 amplitudes and 2 → 2 amplitudes with an additional virtual loop, averaged(summed) over initial(final) states.The structure of each such interference term is The denominators contain massless propagators of the form d i = (l + ∆ i ) 2 + i0, where ∆ i depends only on external momenta, which we will call v i .That is, v i ∈ {p, q, k 1 } or combinations thereof.The trace in (E1) contains Dirac gamma matrices contracted with either the external momenta v i , or with the loop momentum l.Both IR and UV divergences are handled by standard dimensional regularization techniques.After these steps, the numerator in Eq (E1) is a collection of terms with Lorentz invariant products of loop and external momenta.All integrals needed to calculate the virtual corrections in Figure 2 can be written in one of the following ways: where higher rank integrals cannot appear since, in every case, virtual correction diagrams at one-loop involve at least one gluon propagator.We have left the Feynman prescription implicit in the denominators in Eq. (E2), and partonic cross section calculations are done in massless QCD.These tensor integrals can be written in terms the scalar box (D 0 ), triangle (C 0 ) and bubble (B 0 ) integrals using the Passarino-Veltman reduction procedure [75].
As an example, consider the steps to reduce C αβ .Symmetry properties allow C αβ , C β and B β to be written in terms of form factors B Defining the Graham matrix G i,j = v i • v j and the row vectors C = C i2 , the contraction of both sides of Eq. (E3c) with v iα and g αβ respectively gives where R is defined by the relations Here we have used the notation u = v 1 + v 2 .The results for R ij above are obtained by contracting v iα with the right-hand side of Eq. (E3c) and using the property that the products v i • l can always be written in terms of the denominators This allows v iα C αβ to be reduced to rank-1 integrals, whereupon expansion according to Eqs. (E3b), lead to Eqs. (E5).Solving Eq. (E4) gives C ij and C 00 in terms of form factors corresponding to lower rank tensor integrals.In a similar way, one may reduce residual C i and B 1 form factors so that only the scalar integrals C 0 and B 0 appear at the end of the reduction procedure.Equations analogous to (E4) and (E5) for all cases in Eq. (E2) are provided in Appendix A of [74].
After the reduction procedure, all virtual contributions from Figure 2 are in terms of scalar integrals, whose values depend only on Lorentz invariants constructed from momenta that appear in the denominators.The complete set of scalar integrals are in [76].In our calculation, all of these contain single or double pole singularities.
The singular behavior of D 0 and C 0 corresponds to soft and collinear divergences that exactly cancel the soft and collinear singularities in the 2 → 3 calculation of the corresponding channel.
For our computations, we only need expressions for a reduced number of cases for C 0 and D 0 , following the notation of [76]: C 0 divergent type 1 and type 2, D 0 divergent type 2 The usual UV singularities introduced by virtual corrections ultimately are all produced by B 0 -type integrals.UV divergent tadpole integrals A 0 also appear but are zero in a massless theory.Thus, keeping track of B 0 integrals accounts for all UV singularities, which cancel in the sum of virtual graphs.Self-energy diagrams on external legs enter for each channel via the corresponding factors of field strength normalization Z Certain order terms are not needed and may be dropped in calculations at O (α 2 s ).

FIG. 1 .
FIG. 1. Momentum labels in amplitudes for (a) 2 → 2 and (b) 2 → 3 partonic scattering.The dashed lines represent partons of unspecified flavor.The dot on the end of k 1 indicates this is the fragmenting parton.The other momenta are integrated in SIDIS.
2, that cancel after combining 2 → 3 and the corresponding virtual 2 → 2 processes.There are additional collinear singularities in 2 B. Factorization Order-By-Order ). (Recall that α 1 and α 2 can be negative or zero.)This example had two ADMVs in the denominator.If the number of ADMVs in the denominator is 1, one can express the ADMVs in the numerator in terms of two ADMVs of different types, where one is the ADMV in the denominator and the other is chosen randomly.If the number of ADMVs in the denominator is 0, then the ADMVs in the numerator can be expressed by any two randomly chosen ADMVs of different types.There are six ADMVs, but the four relations from Eqs. (D5)-(D8) eliminate four in terms of the other two.