text:
# x Hard Scattering Factorization from Effective Field Theory

## I Introduction

The principle of factorization underlies all theoretical predictions for hadronic processes. Simply put, factorization is the statement that short and long distances contributions to physical processes can be separated, up to corrections suppressed by powers of the relevant large scale in the process. The predictive power gained from this result stems from the fact that the incalculable long distance effects are universal, defined in an unambiguous way in terms of matrix elements. As a consequence, the non-perturbative long distance effects can be extracted in one process and then used in another. In general, proving factorization is a difficult task Ellis:ty . The proof of factorization in Drell Yan processes, for instance, took several years to sort out DYrefs  (for reviews on factorization see Collins:1987pm ; Sterman:1995fz ; Jaffe:zw ). Indeed, there are still some processes such as $`B\pi \pi `$ where a proof of factorization only exists at one-loop Bpipi .

Given that we would like to retain our predictive power over the largest possible range of energies, we are compelled to understand power corrections to the factorized rates. These corrections are not necessarily universal, and as such, the relevant size of the power corrections are process dependent. In processes for which there exists an operator product expansion (OPE), there is a systematic way in which to include power corrections. However, for a majority of observables we do not have an OPE at our disposal, and the nature of the power correction is not always known. For instance, in the case of shape variables there is still some on going discussion about the form of subleading corrections nason ; Korchemsky:1999kt .

The purpose of this paper is to show that an effective theory framework can be used to simplify proofs of factorization and describe processes with an operator formalism. To do this we extend the soft-collinear effective theory (SCET) developed in Refs. bfl ; bfps ; cbis ; bpssoft ; bps , to high energy processes. It should be emphasized that there are several other useful advantages in using an effective field theory (EFT). For instance, the EFT makes any symmetries which emerge in the $`Q\mathrm{}`$ limit manifest in the Lagrangian and operators, and allow statements to be made to all orders in perturbation theory. The calculation of hard coefficients reduces to simple matching calculations, where subtracting the EFT graphs automatically removes all infrared divergences from the QCD calculation. Perhaps most importantly, it provides a framework for systematically investigating power corrections. Finally, the EFT framework allows standard renormalization group techniques to be used for the resummation of logarithms that are often necessary in calculating rates for certain high energy scattering events bfl ; bfps ; resum ; resum2 . The factorization formulae that we prove in this paper are not new, but serve to illustrate our approach in familiar settings. The results are valid to all orders in $`\alpha _s`$ and leading order in the power expansion. The simplicity of our approach lies in the fact that factorization occurs at the level of the SCET Lagrangian and operators, and is facilitated by gauge symmetry in the EFT. This provides the advantage that our proofs do not rely on making use of Ward identities and induction, or on specifying a particular gauge.666In fact our factorization proofs rely heavily on the gauge symmetry structure of SCET. When a gauge fixing term is required for explicit calculations we use general covariant gauges. Furthermore, it becomes rather simple to derive factorization formulae for a myriad of processes, since many results are universal. The examples given here serve to illustrate these simplifications. Developments on the issues of power corrections and resummations are left to future publications.777For recent work on subleading corrections in SCET for heavy-to-light transitions see Ref. Chay:2002vy .

In section II we review the construction of the SCET. The formalism developed in Refs. bfl ; bfps ; cbis ; bpssoft  is extended to include two types of collinear particles moving in opposite directions in section II.2, and factorization for $`\gamma ^{}`$ to two collinear states is discussed as an example. In section II.3 we define the non-perturbative matrix elements such as the parton distribution functions that will be needed for the processes presented in the paper, and in section II.4 we discuss some of the symmetries in SCET that may be used to place restrictions on matrix elements. In the remaining sections we give various examples on how factorization theorems emerge in the effective theory language. In section III we prove factorization theorems for two exclusive processes, namely the $`\pi `$-$`\gamma `$ form factor, and meson form factors ($`\gamma ^{}MM`$) for arbitrary spin and isospin structure. In section IV two inclusive processes are treated, namely DIS ($`e^{}pe^{}X`$) and Drell-Yan ($`p\overline{p}X\mathrm{}^+\mathrm{}^{}`$), and we also give results for deeply virtual Compton scattering ($`\gamma ^{}p\gamma ^{()}p`$). In these processes we include all leading power contributions in the factorization proofs (even if the operators are only matched onto at higher orders in perturbation theory such as for the gluon distribution functions). Our conclusions are given in section V. In appendix A we show how auxiliary fields can be used to prove the simultaneous factorization of soft fields from collinear fields for particles in back-to-back directions.

## II Formalism

Effective field theories provide a simple and elegant way of organizing physics in processes containing widely disparate energy scales. In constructing an EFT, some degrees of freedom are eliminated, and the remaining degrees of freedom must reproduce all the infrared physics of the full theory in the domain where the EFT is valid. The EFT is organized by an expansion in $`\lambda `$, defined as the ratio of small to large energy scales. As a useful guideline the following steps are used to identify the infrared degrees of freedom: 1) Determine the relevant scales in a problem from the size of the momenta and masses of all particles that can make up the initial and final states, 2) Construct all momenta from these scales whose components correspond to propagating degrees of freedom, and which have offshellness less than the large scale, i.e. $`p^2m^2Q^2`$. Effective theory fields are then constructed for each unique set of these momenta.

We will be interested in an EFT with particles of energy $`Q`$ much greater than their mass. The dynamics of these particles can be described by constructing a soft-collinear effective theory (SCET). This theory is organized as an expansion in powers of $`\lambda p_{}/Q`$, and offshell fluctuations with $`p^2(Q\lambda )^2`$ are integrated out. In section II.1 we begin by describing this procedure and comparing the construction to other EFT’s. We then give a brief review of the soft-collinear effective theory developed in Refs. bfl ; bfps ; cbis ; bpssoft ; bps . We do not attempt to give a comprehensive treatment, but instead emphasize the main results and refer the reader to the literature for details. In section II.2 we extend the formulation of SCET to describe processes with collinear particles moving in back to back directions, and prove the factorization formula for $`\gamma ^{}`$ to two collinear states as an example. In section II.3 we define the non-perturbative matrix elements that are needed for our examples, then in section II.4 we discuss some of the symmetry properties of collinear fields and currents.

### II.1 Soft-Collinear Effective Theory

In the standard construction of an EFT one removes the short distance scales and massive fields by integrating them out one at a time. A classical example is integrating out the $`W`$ boson to obtain the effective electroweak Hamiltonian with 4-fermion operators. However, in some situations we are interested in integrating out large momentum fluctuations without fully removing the corresponding field. The simplest example of this is Heavy Quark Effective Theory (HQET) bbook , which is constructed to describe the low energy properties of mesons with a heavy quark. Here the heavy anti-quarks are integrated out and only heavy quarks with fluctuations close to their mass-shell are retained. This is accomplished by removing fluctuations of order the heavy quark mass $`m_Q`$ with a field redefinitiongeorgi

$$\psi (x)=\underset{v}{}e^{im_Qvx}h_v(x),$$
(1)
where $`v`$ is the heavy quark velocity and $`h_v`$ is the field in the EFT. While $`^\mu \psi (x)m_Q\psi (x)`$, the effective field has $`^\mu h_v(x)\mathrm{\Lambda }_{\mathrm{QCD}}h_v(x)`$, indicating that it no longer describes short-distance fluctuations about the perturbative scale $`m_Q`$. Instead these effects are encoded in calculable Wilson coefficients. The HQET degrees of freedom with offshellness $`p^2\mathrm{\Lambda }_{\mathrm{QCD}}^2`$ are the heavy quarks, soft gluons, and soft quarks.

Similarly, for collinear particles with energy $`Qm`$, one needs to remove momentum fluctuations $`Q`$ while retaining effective theory fields to describe smaller momenta. However, unlike heavy quarks the collinear particles have two low energy scales. Consider the light-cone momenta, $`p^+=np`$ and $`p^{}=\overline{n}p`$ where $`n^2=\overline{n}^2=0`$ and $`n\overline{n}=2`$. Here $`n`$ parameterizes a light-cone direction close to that of the collinear particle and $`\overline{n}`$ the opposite direction (eg. for motion in the $`z`$ direction $`n_\mu =(1,0,0,1)`$ and $`\overline{n}_\mu =(1,0,0,1)`$). For a particle of mass $`mp_{}Q`$, we have $`p^{}Q`$, and a small parameter $`\lambda p_{}/Q`$. The scaling of the $`p^+`$ component is then fixed by the equations of motion $`p^+p^{}+p_{}^2=m^2`$, so that $`(p^+,p^{},p_{})Q(\lambda ^2,1,\lambda )`$.

The appearance of two small scales, $`Q\lambda ^2Q\lambda Q`$, is similar to the situation in non-relativistic QCD (NRQCD), which is an EFT for systems of two heavy quarks with an expansion in their relative velocity $`\beta `$. In a non-relativistic bound state the momentum of a heavy quark is $`𝐩m_Q\beta `$, but the equations of motion $`E=𝐩^2/(2m_Q)`$ make the energy $`Em_Q\beta ^2`$, giving scales $`m_Q\beta ^2m_Q\beta m_Q`$. The two low energy scales can be distinguished by following Eq. (1) with a further field redefinition LMR  $`h_v(x)=_𝐩e^{i𝐩𝐱}\psi _𝐩(x)`$, so that derivatives on $`\psi _𝐩`$ only pick out the $`m\beta ^2`$ scale. The on-shell degrees of freedom are then the heavy quarks, soft quarks and gluons with $`p^2(m_Qv)^2`$, and ultrasoft quarks and gluons with $`p^2(m_Qv^2)^2`$.

SCET fields:
For collinear particles the analogous field redefinitions are bfps ; cbis

$$\varphi (x)=\underset{n}{}\underset{p}{}e^{ipx}\varphi _{n,p}(x),$$
(2)
where the collinear fields $`\varphi _{n,p}`$ are labelled by light-cone vectors $`n`$ and label momentum $`p`$. Here $`p`$ contains the $`\overline{n}pQ`$ and $`p_{}Q\lambda `$ momenta so that $`^\mu \varphi _{n,p}(Q\lambda ^2)\varphi _{n,p}`$. The field $`\varphi _p`$ can either be a quark or gluon field. Similar to $`\psi _𝐩`$, the missing $`Q`$ fluctuations are described by Wilson coefficients and the $`Q\lambda `$ labels simplify the power counting by distinguishing the $`Q\lambda `$ and $`Q\lambda ^2`$ scales. Now

$`\varphi _{n,p}\varphi _{n,p}^++\varphi _{n,p}^{},`$ (3)
so collinear particles and antiparticles are contained in the same effective theory field, but have momentum labels with the opposite sign. In the large energy limit the four component fermion spinors contain two large and two small components. One therefore defines collinear quark fields $`\xi _{n,p}`$ which only retain the large components for motion in the $`n`$ direction and satisfies $`n/\xi _{n,p}=0`$. For these fields $`\xi _{n,p}^+/\xi _{n,p}^{}`$ destroy/create the particles/antiparticles with large momentum $`\overline{n}p>0`$ cbis . For collinear gluons $`A_{n,q}^\mu =A_{n,q}^\mu `$, and $`(A_{n,q}^\mu )^+`$/$`(A_{n,q}^\mu )^{}`$ destroy/create gluons with $`\overline{n}q>0`$.

For simplicity we will ignore quark masses and only consider massless $`u`$ and $`d`$ quarks. For the processes considered here SCET then requires three types of degrees of freedom: collinear, soft, and ultrasoft (usoft) fields. These are distinguished by the scaling of the light cone components ($`p^+`$,$`p^{}`$,$`p^{}`$) of their momenta: $`(\lambda ^2,1,\lambda )`$ for collinear modes in the $`n`$ direction ($`A_{n,q}`$, $`\xi _{n,p}`$), $`(\lambda ,\lambda ,\lambda )`$ for the soft modes ($`A_q^s`$, $`q_p^s`$), and $`(\lambda ^2,\lambda ^2,\lambda ^2)`$ for the usoft modes ($`A_{us}`$, $`q_{us}`$). The soft modes are labelled by their order $`Q\lambda `$ momenta, so $`A_q^s`$ and $`q_p^s`$ are essentially just momentum space fields. The usoft fields have no labels and depend only on the coordinate $`x`$. The fields are assigned a scaling with $`\lambda `$ to make the action for their kinetic terms order $`\lambda ^0`$ bfl ; bfps ; bpssoft . For instance $`\xi _{n,p}\lambda `$, $`A_{n,q}^\mu (\lambda ^2,1,\lambda )`$, $`A_q^s\lambda `$, and $`A_{us}^\mu \lambda ^2`$. At leading order only order $`\lambda ^0`$ vertices are necessary to correctly account for all order $`\lambda ^0`$ Feynman diagrams.

In HQET only external currents with momenta of order $`m_b`$ can change the label $`v`$. Thus the Lagrangian has a superselection rule forbidding changes in the four-velocity of the heavy quark bbook ; georgi . In NRQCD the $`v`$ labels are also conserved, but the smaller momentum labels $`𝐩`$ are changed by operators in the effective theory such as the Coulomb potential. A novel feature of SCET is that interactions in the leading action can change both the large and small parts of the momentum labels $`p^\mu `$. However, only external currents can change the direction $`n`$ of a collinear particle, so this label is conserved. Thus, for each distinct direction $`n`$ a separate set of collinear fields are needed. In the remainder of this section we will restrict ourselves to collinear particles with a single $`n`$. We will generalize the discussion to the case of two back-to-back directions and discuss the factorization of collinear particles with different $`n`$’s in section II.2.

Since in SCET interactions can change the order $`Q`$ label momenta it turns out to be very useful to introduce a label operator, $`𝒫^\mu `$ cbis , for which the collinear fields satisfy $`𝒫^\mu \xi _{n,p}=p^\mu \xi _{n,p}`$. More generally, $`𝒫^\mu `$ acts on a product of labelled fields as

$`f(𝒫^\mu )\left(\varphi _{q_1}^{}\mathrm{}\varphi _{q_m}^{}\varphi _{p_1}\mathrm{}\varphi _{p_n}\right)=f(p_1^\mu +\mathrm{}+p_n^\mu q_1^\mu \mathrm{}q_m^\mu )\left(\varphi _{q_1}^{}\mathrm{}\varphi _{q_m}^{}\varphi _{p_1}\mathrm{}\varphi _{p_n}\right).`$ (4)
so conjugate field labels come with a minus sign. The operator $`𝒫_\mu `$ acts to the right, while the conjugate operator $`𝒫_\mu ^{}`$ acts to the left. As explained in Ref. cbis  the label operator allows all large phases to be moved to the front of operators with a factor $`\mathrm{exp}(ix𝒫)`$. This phase and the label sums can then be suppressed if we impose that interactions conserve label momenta and that the momentum indices on fields are implicitly summed over. Basically, for labels $`p`$ and $`p^{}`$ and residual momenta $`k`$ and $`k^{}`$

$`{\displaystyle d^4xe^{i(p^{}p+k^{}k)x}}=\delta (pp^{}){\displaystyle d^4xe^{i(k^{}k)x}},`$ (5)
so that the label and residual momenta are individually conserved. (Although technically the label momenta are discrete we abuse notation and use $`\delta (pp^{})`$ rather than $`\delta _{p,p^{}}`$ because it makes the subscripts easier to read.) For convenience we define the operator $`\overline{𝒫}`$ to pick out only the order $`\lambda ^0`$ labels on collinear fields, and the operator $`𝒫^\mu `$ to pick out only the order $`\lambda `$ labels. For the matrix element of any collinear operator $`𝒪`$, momentum conservation constrains the sum of field labels cbis , giving

$`M_{n,p_1}\left|\left[f(\overline{𝒫})𝒪\right]\right|M_{n,p_2}`$ $`=`$ $`f(\overline{n}(p_2p_1))M_{n,p_1}\left|𝒪\right|M_{n,p_2},`$ (6)
for any function $`f`$.

For a single $`n`$ the Lagrangian can be broken up into three sectors: collinear, usoft, and soft. We therefore write

$`=_{c,n}[\xi _{n,p},A_{n,q}^\mu ,A_{us}^\mu ]+_{us}[q_{us},A_{us}^\mu ]+_s[q_{s,p},A_{s,q}^\mu ],`$ (7)
where we have made the field content of each sector explicit. We will discuss each of these terms separately.

Collinear sector:
As explained in detail in Ref. bpssoft , gauge invariance in SCET restricts the Lagrangian and allowed form of operators. Only local gauge transformations whose action is closed on the effective theory fields need to be considered. These include collinear, soft, and usoft transformations. Each of these vary over different distance scales, with collinear gauge transformations satisfying $`^\mu U_n(x)Q(\lambda ^2,1,\lambda )U_n(x)`$, soft satisfying $`^\mu V_s(x)Q\lambda V_s(x)`$, and usoft transformations with $`^\mu V_{us}(x)Q\lambda ^2V_{us}(x)`$. All particles transform under $`V_{us}(x)`$ and usoft gluons act like background fields for collinear particles. Invariance under $`U_n(x)`$ requires a collinear Wilson line built out of the order $`\lambda ^0`$ gluon fields bfps ; cbis

$`W_n(x)=\left[\text{ }\stackrel{{\scriptscriptstyle }}{\text{perms}}\text{ }\mathrm{exp}\left(g{\displaystyle \frac{1}{\overline{𝒫}}}\overline{n}A_{n,q}(x)\right)\right].`$ (8)
Here the operator $`\overline{𝒫}`$ acts only inside the square brackets, the $`n`$ on $`W_n`$ refers to the direction of the collinear quanta, and $`W_n`$ is local with respect to $`x`$ (corresponding to the residual momenta). Taking the Fourier transform of $`\delta (\omega \overline{𝒫})W_n(0)`$ with respect to $`\omega `$ gives the more familiar path-ordered Wilson line $`W_n(y,\mathrm{})=\mathrm{P}\mathrm{exp}\left[ig_{\mathrm{}}^yds\overline{n}A_n(s\overline{n})\right]`$. Under a collinear gauge transformation $`W_n`$ transforms as $`W_nU_nW_n`$. An invariant under collinear gauge transformations can therefore be formed by combining a collinear fermion $`\xi _{n,p}`$ and the Wilson line $`W_n`$ in the form

$`W_n^{}(x)\xi _{n,p}(x).`$ (9)
This combination still transforms under an usoft gauge transformation, $`W_n^{}\xi _{n,p}V_{us}(x)W_n^{}\xi _{n,p}`$. We will often suppress the $`x`$ dependence of the combination $`W_n^{}\xi _{n,p}`$.

Integrating out hard fluctuations gives Wilson coefficients in the effective theory that are functions of the large $`\overline{n}p_i`$ collinear momenta, $`C(\overline{n}p_i)`$. However, collinear gauge invariance restricts these coefficients to only depend on the linear combination of momenta picked out by the order $`\lambda ^0`$ operator $`\overline{𝒫}`$ cbis . In general the Wilson coefficients are then functions $`C(\overline{𝒫},\overline{𝒫}^{})`$ which must be inserted between gauge invariant products of collinear fields. In general the Wilson coefficients also depend on the large momentum scales in a process such as $`Q`$ and the renormalization scale $`\mu `$.

To construct the collinear Lagrangian one can match full QCD onto operators with collinear fields that are invariant under usoft and collinear gauge transformations. The collinear Lagrangian at order $`\lambda ^0`$ is bfps ; cbis ; bpssoft

$`_{c,n}`$ $`=`$ $`\overline{\xi }_{n,p^{}}\{inD+gnA_{n,q}+(𝒫/_{}+gA/_{n,q}^{})W{\displaystyle \frac{1}{\overline{𝒫}}}W^{}(𝒫/_{}+gA/_{n,q^{}}^{})\}{\displaystyle \frac{\overline{n}/}{2}}\xi _{n,p}`$ (10)
$`+{\displaystyle \frac{1}{2g^2}}\mathrm{tr}\left\{[i𝒟^\mu +gA_{n,q}^\mu ,i𝒟^\nu +gA_{n,q^{}}^\nu ]^{\mathrm{\hspace{0.17em}2}}\right\}+_c^{g.f.},`$
where $`_c^{g.f.}`$ are gauge fixing terms, $`iD^\mu =i^\mu +gA_{us}^\mu `$, and

$`i𝒟^\mu ={\displaystyle \frac{n^\mu }{2}}\overline{𝒫}+𝒫_{}^\mu +{\displaystyle \frac{\overline{n}^\mu }{2}}inD.`$ (11)
Since usoft gluons act as background fields in the collinear gauge transformation the couplings, $`g(\mu )`$, for both types of gluons must be identical.

Usoft and Soft sectors:
The usoft and soft Lagrangians for gluons and massless quarks are the same as those in QCD. From Eq. (7) we see that collinear quarks and gluons interact with usoft gluons, however at order $`\lambda ^0`$ only the $`nA_{us}`$ component appears in Eq. (10). In order to prove factorization formulae it is essential to disentangle the collinear and usoft modes. This can be done by introducing an usoft Wilson line

$`Y_n(x)`$ $`=`$ $`\mathrm{P}\mathrm{exp}\left(ig{\displaystyle _{\mathrm{}}^x}dsnA_{us}(sn)\right),`$ (12)
where the subscript $`n`$ on $`Y_n`$ labels the direction of the Wilson line (we emphasize that this is different from the meaning of the subscript on $`W_n`$ in Eq. (8)). An usoft gauge transformation takes $`Y_nV_{us}Y_n`$. In Ref. bpssoft  it was shown that the field redefinitions

$`\xi _{n,p}=Y_n\xi _{n,p}^{(0)},A_{n,p}^\mu =Y_nA_{n,p}^{(0)\mu }Y_n^{},`$ (13)
imply $`W_n=Y_nW_n^{(0)}Y_n^{}`$ and decouple the usoft gluons from the collinear particles in the leading order Lagrangian

$`_{c,n}[\xi _{n,p},A_{n,q}^\mu ,nA_{us}]=_{c,n}[\xi _{n,p}^{(0)},A_{n,q}^{(0)\mu },0].`$ (14)
Thus, the new collinear fields with superscript $`(0)`$ no longer interact with usoft gluons or transform under an usoft gauge transformation. Since the field redefinitions do not change physical $`S`$ matrix elements, the new fields give an equally valid parameterization of the collinear modes. The leading SCET Lagrangian therefore factors into separate collinear and usoft sectors. This alone does not guarantee factorization in operators and currents, since after the field redefinition these operators may still contain both usoft and collinear fields. However, the field redefinition makes factorization transparent since identities such as $`Y_n^{}Y_n=1`$ may be applied directly to the operators. This will become clear in the examples in sections III and IV.

The coupling of soft gluons to collinear particles differs from the usoft-collinear interactions. Interactions of a soft gluon with a collinear particle results in a particle with momentum $`pQ(\lambda ,1,\lambda )`$, so soft gluons can not appear in the collinear Lagrangian. These offshell particles have $`p^2Q^2\lambda `$ and since $`Q^2\lambda (Q\lambda )^2`$ these offshell quarks and gluons can be integrated out. At leading order in $`\lambda `$ it was shown in Ref. bpssoft  that in operators with collinear fields this simply builds up factors of a soft Wilson line $`S_n`$ involving the $`nA_s`$ component of the soft gluon field,

$`S_n`$ $`=`$ $`\left[\text{ }\stackrel{{\scriptscriptstyle }}{\text{perms}}\text{ }\mathrm{exp}\left(g{\displaystyle \frac{1}{n𝒫}}nA_{s,q}\right)\right].`$ (15)
The factors of $`S_n`$ appear outside gauge invariant products of collinear fields, and their location is restricted by soft gauge invariance.

### II.2 SCET for $`n`$ and $`\overline{n}`$ collinear fields

In this section we extend SCET to include the possibility of collinear fields moving in different light-cone directions: $`n_1`$, $`n_2`$, $`n_3,\mathrm{}`$. These directions can be considered to be distinct provided that $`n_in_j\lambda ^2`$ for $`ij`$. This follows from the fact that if $`n_1n_2\lambda ^2`$ then the directions $`n_1`$ and $`n_2`$ are too close to be distinguished. For example, a momentum $`p_2=Qn_2`$ can be considered to be collinear in the $`n_1`$ direction if $`n_1p_2=Qn_1n_2Q\lambda ^2`$, since this is the correct scaling for the small momentum component of an $`n_1`$-collinear particle.

For simplicity we will only consider the case of back-to-back jets corresponding to collinear particles moving in the $`n`$ and $`\overline{n}`$ directions. These are clearly distinct since $`n\overline{n}=2`$. Collinear particles in the $`\overline{n}`$ direction have $`(+,,)`$ momenta $`Q(1,\lambda ^2,\lambda )`$, and the $`np1`$ and $`p_{}\lambda `$ momenta appear as labels on the corresponding fields: $`\xi _{\overline{n},p}`$ and $`A_{\overline{n},p}^\mu `$. Emission of a collinear particle moving in the $`n`$ direction from a collinear particle in the $`\overline{n}`$ direction results in a particle with momentum $`kQ(1,1,\lambda )`$ and offshellness $`k^2Q^2`$. These offshell modes are integrated out to construct the SCET, so collinear modes in the $`n`$ direction do not directly couple to collinear modes in the $`\overline{n}`$ direction. A distinct set of collinear gauge transformations is associated with each of $`n`$ and $`\overline{n}`$, and fields in one direction do not transform under the gauge symmetry associated with the opposite direction. Two Wilson lines $`W_n(x)`$ and $`W_{\overline{n}}(x)`$ are necessary (defined as in Eq. (8)), and they appear in a way that makes collinear operators gauge invariant. For instance the combinations

$`W_n^{}\xi _{n,p},W_{\overline{n}}^{}\xi _{\overline{n},p}`$ (16)
are invariant under collinear gauge transformations in the $`n`$ and $`\overline{n}`$ directions, respectively. We also require two types of label operators, $`\overline{𝒫}`$ as before, and an operator $`𝒫`$ to pick out $`np`$ labels that are order $`\lambda ^0`$. Thus, $`\overline{𝒫}`$ and $`𝒫`$ act only on the $`n`$ and $`\overline{n}`$ collinear fields respectively. (The label operator $`𝒫^\mu `$ still picks out order $`\lambda `$ momentum components and therefore acts on both $`n`$ and $`\overline{n}`$ fields.) With two collinear directions, decoupling usoft gluons requires introducing both $`Y_n`$ and $`Y_{\overline{n}}`$ Wilson lines, defined as in Eq. (12), but along the $`n`$ or $`\overline{n}`$ paths respectively. Finally, integrating out $`Q^2\lambda `$ fluctuations at leading order induces both $`S_n`$ and $`S_{\overline{n}}`$ soft Wilson lines defined analogous to Eq. (15). This is discussed in greater detail in Appendix A where we show explicitly to all orders in $`g`$ that integrating out the $`Q^2\lambda `$ fluctuations causes

$`W_n^{}\xi _{n,p}S_nW_n^{}\xi _{n,p},\overline{\xi }_{n,p}W_n\overline{\xi }_{n,p}W_nS_n^{},`$
$`W_{\overline{n}}^{}\xi _{\overline{n},p}S_{\overline{n}}W_{\overline{n}}^{}\xi _{\overline{n},p},\overline{\xi }_{\overline{n},p}W_{\overline{n}}\overline{\xi }_{\overline{n},p}W_{\overline{n}}S_{\overline{n}}^{}.`$ (17)
Relations for operators with collinear gluon fields are also derived in Appendix A.

Note that we have not included ”Glauber gluons” with momenta $`p^\mu (\lambda ^2,\lambda ^2,\lambda )`$, which are kinematically allowed in $`t`$-channel Coulomb exchange between $`n`$ and $`\overline{n}`$ collinear quarks. In determining the relevant degrees of freedom we have assumed that Glauber gluons are not necessary to describe the infrared for the processes considered in this paper. Intuitively, this can be seen from the fact these gluons are instantaneous in both time and longitudinal separation, and only could contribute when the $`n`$ and $`\overline{n}`$ jets overlap for a duration of order $`1/(Q\lambda ^2)`$ in a space-time diagram. In processes with a hard interaction the overlap scale is always much shorter than this (however this need not be the case in processes such as forward scattering). For the Drell-Yan process more quantitative arguments can be found in Refs. DYrefs ; CSSnew .

At order $`\lambda ^0`$ it is not possible to construct a gauge invariant kinetic Lagrangian with terms that involve both $`n`$ and $`\overline{n}`$ fields. Thus, the $`n`$ and $`\overline{n}`$ collinear modes are described by independent Lagrangians (however $`n`$ and $`\overline{n}`$ modes may still both appear in an external operator). The collinear sector of the SCET Lagrangian is therefore

$`_{c,n}[\xi _{n,p},A_{n,q}^\mu ,nA_{us}]+_{c,\overline{n}}[\xi _{\overline{n},p},A_{\overline{n},q}^\mu ,\overline{n}A_{us}].`$ (18)
Making the field redefinitions

$`\xi _{n,p}`$ $`=`$ $`Y_n\xi _{n,p}^{(0)},A_{n,p}^\mu =Y_nA_{n,p}^{(0)\mu }Y_n^{},`$
$`\xi _{\overline{n},p}`$ $`=`$ $`Y_{\overline{n}}\xi _{\overline{n},p}^{(0)},A_{\overline{n},p}^\mu =Y_{\overline{n}}A_{\overline{n},p}^{(0)\mu }Y_{\overline{n}}^{},`$ (19)
gives $`W_n=Y_nW_n^{(0)}Y_n^{}`$, $`W_{\overline{n}}=Y_{\overline{n}}W_{\overline{n}}^{(0)}Y_{\overline{n}}^{}`$, and usoft degrees of freedom once again decouple from the collinear modes since

$`_{c,n}[\xi _{n,p},A_{n,q}^\mu ,nA_{us}]+_{c,\overline{n}}[\xi _{\overline{n},p},A_{\overline{n},q}^\mu ,\overline{n}A_{us}]=_{c,n}[\xi _{n,p}^{(0)},A_{n,q}^{(0)\mu },0]+_{c,\overline{n}}[\xi _{\overline{n},p}^{(0)},A_{\overline{n},q}^{(0)\mu },0].`$ (20)
Thus, usoft gluons are removed from the collinear Lagrangian at the expense of inducing $`Y_n`$ and $`Y_{\overline{n}}`$ factors in operators with collinear fields. In certain cases the identities $`Y_n^{}Y_n=1`$ and $`Y_{\overline{n}}^{}Y_{\overline{n}}=1`$ can be used in these operators to cancel usoft gluon interactions. Perturbatively these cancellations would occur by adding an infinite set of Feynman diagrams.

To see in more detail how this works consider the simple example of the $`\gamma ^{}`$-production of back-to-back collinear states $`X_n`$ and $`X_{\overline{n}}`$. The full theory current $`\overline{\psi }(x)\mathrm{\Gamma }\psi (x)`$ matches onto an effective theory operator $`O_{n\overline{n}}`$. Naively one might guess that the SCET operator mediating this process is

$`O_{n\overline{n}}=\overline{\xi }_{n,p_1}\mathrm{\Gamma }\xi _{\overline{n},p_2}.`$ (21)
However, this operator is not invariant under the collinear gauge transformations $`U_n`$ and $`U_{\overline{n}}`$, so the process is instead mediated by the invariant operator

$`O_{n\overline{n}}=\overline{\xi }_{n,p_1}W_n\mathrm{\Gamma }W_{\overline{n}}^{}\xi _{\overline{n},p_2}.`$ (22)
A hard matching coefficient $`C(\overline{𝒫},\overline{𝒫}^{},𝒫,𝒫^{})`$ can be inserted in any location in the operator that does not break apart the gauge invariant combinations of fields in Eq. (16). The operators $`\overline{𝒫}`$ and $`𝒫`$ in the coefficient only pick out momenta that are order $`\lambda ^0`$ in the power counting. Thus, $`𝒫`$ does not act on fields in the $`n`$ direction and $`\overline{𝒫}`$ does not act on fields in the $`\overline{n}`$ direction, and the most general result is

$`O_{n\overline{n}}=\overline{\xi }_{n,p_1}W_n\mathrm{\Gamma }C(\overline{𝒫}^{},𝒫)W_{\overline{n}}^{}\xi _{\overline{n},p_2}.`$ (23)
Next, we integrate out the offshell $`Q^2\lambda `$ fluctuations which induces additional soft Wilson lines in the operator. This is discussed in detail in Appendix A and from Eq. (II.2) gives

$`O_{n\overline{n}}=\overline{\xi }_{n,p_1}W_nS_n^{}\mathrm{\Gamma }C(\overline{𝒫}^{},𝒫)S_{\overline{n}}W_{\overline{n}}^{}\xi _{\overline{n},p_2}.`$ (24)
Note that $`\overline{𝒫}`$ and $`𝒫`$ do not act on the fields in the soft Wilson lines since soft gluons carry only order $`\lambda `$ momenta. Finally, we can make the usoft gluon couplings explicit by switching to the $`(0)`$ fields using Eq. (II.2):

$`O_{n\overline{n}}=\overline{\xi }_{n,p_1}^{(0)}W_n^{(0)}Y_n^{}S_n^{}\mathrm{\Gamma }C(\overline{𝒫}^{},𝒫)S_{\overline{n}}Y_{\overline{n}}W_{\overline{n}}^{(0)}\xi _{\overline{n},p_2}^{(0)}.`$ (25)
This operator is manifestly invariant under collinear gauge transformation in the $`n`$ and $`\overline{n}`$ direction, as well as under soft and usoft gauge transformations.

To separate the short distance Wilson coefficient from the long-distance operator one introduces convolution variables $`\omega `$ and $`\omega ^{}`$ to give

$`O_{n\overline{n}}`$ $`=`$ $`{\displaystyle 𝑑\omega 𝑑\omega ^{}C(\omega ,\omega ^{})O_{n\overline{n}}(\omega ,\omega ^{})},`$ (26)
$`O_{n\overline{n}}(\omega ,\omega ^{})`$ $`=`$ $`\left[\overline{\xi }_{n,p_1}^{(0)}W_n^{(0)}\delta (\overline{𝒫}^{}\omega )Y_n^{}S_n^{}\mathrm{\Gamma }S_{\overline{n}}Y_{\overline{n}}\delta (𝒫\omega ^{})W_{\overline{n}}^{(0)}\xi _{\overline{n},p_2}^{(0)}\right].`$
The function $`C(\omega ,\omega ^{})`$ contains all the short distance physics and is determined by matching the full theory onto this effective theory operator. $`O_{n\overline{n}}(\omega ,\omega ^{})`$ contains all the infrared long-distance QCD contributions at leading order in $`\lambda `$.

Now consider the matrix element of the production current between $`X_nX_{\overline{n}}|`$ and the vacuum. Taking the $`\gamma ^{}`$ to have large time-like momentum $`q^\mu =(Q,0,0,0)`$ (and zero residual momentum) we have

$`{\displaystyle d^4xe^{iqx}X_nX_{\overline{n}}\left|\overline{\psi }(x)\mathrm{\Gamma }\psi (x)\right|0}={\displaystyle d^4xX_nX_{\overline{n}}\left|O_{n\overline{n}}(x)\right|0}`$
$`={\displaystyle d^4xe^{ikx}(X_nX_{\overline{n}})(k)\left|O_{n\overline{n}}(x=0)\right|0}=(X_nX_{\overline{n}})(0)\left|O_{n\overline{n}}(0)\right|0.`$ (27)
In the first step the conservation of the large label momentum $`q`$ was made implicit in the matrix element (c.f. Eq. (5)). Since we are in the center-of-mass frame the $`X_n`$ has large momentum $`\overline{n}p=Q`$, and the $`X_{\overline{n}}`$ has momentum $`np^{}=Q`$. Now using translation invariance, we see that the remaining $`x`$ integral forces the $`|X_nX_{\overline{n}}`$ state to have zero residual momentum. Using Eq. (26) this matrix element is equal to

$`{\displaystyle 𝑑\omega 𝑑\omega ^{}C(\omega ,\omega ^{})X_n\overline{X}_{\overline{n}}\left|O_{n\overline{n}}(\omega ,\omega ^{})\right|0}={\displaystyle 𝑑\omega 𝑑\omega ^{}C(\omega ,\omega ^{})J_n(\omega )\mathrm{\Gamma }S_{n\overline{n}}J_{\overline{n}}(\omega ^{})},`$ (28)
where the functions $`C`$, $`J_n`$, $`J_{\overline{n}}`$, and $`S`$ also depend on the renormalization point $`\mu `$. Here we have used the fact that both $`X_n`$ and $`X_{\overline{n}}`$ can be described entirely by collinear particles in the $`n`$ and $`\overline{n}`$ directions respectively. Since the Lagrangians for the collinear, soft, and usoft fields are factorized the remaining matrix element splits into distinct matrix elements for each class of modes. These matrix elements are

$`J_n(\omega )`$ $`=`$ $`X_n(Q/2)\left|\overline{\xi }_{n,p_1}^{(0)}W_n^{(0)}\delta (\overline{𝒫}^{}\omega )\right|0,`$
$`J_{\overline{n}}(\omega ^{})`$ $`=`$ $`X_{\overline{n}}(Q/2)\left|\delta (𝒫\omega ^{})W_{\overline{n}}^{(0)}\xi _{\overline{n},p_2}^{(0)}\right|0,`$
$`S_{n\overline{n}}`$ $`=`$ $`0\left|Y_n^{}S_n^{}S_{\overline{n}}Y_{\overline{n}}\right|0,`$ (29)
(and are matrices whose color, spin, and flavor indices are suppressed). Note that $`J_n`$, $`J_{\overline{n}}`$, and $`S_{n\overline{n}}`$ are explicitly invariant under the collinear, soft, and usoft gauge transformations bpssoft  of SCET, but still transform globally under a color rotation. Now using the momentum conservation identity in Eq. (6), the large momentum of the $`X_n`$ and $`X_{\overline{n}}`$ states set $`𝒫Q`$ and $`\overline{𝒫}^{}Q`$. Label conservation also implies that the total perpendicular momentum of each of $`J_n`$, $`J_{\overline{n}}`$, and $`S_{n\overline{n}}`$ is zero. The sum over $`\omega `$ and $`\omega ^{}`$ can then be performed to give the final factorized form

$`C(Q,Q)J_n(Q)\mathrm{\Gamma }S_{n\overline{n}}J_{\overline{n}}(Q).`$ (30)
Although rather idealized, the above example illustrates the main steps needed to derive a factorization formula. Taking $`X_n`$ and $`X_{\overline{n}}`$ to be single quark states the result in Eq. (30) also agrees with the factorization formula for the onshell production form factor for $`q\overline{q}`$ Korchemskyff ; CollinsP .888In this case depending on the choice of infrared regulator(s), it may not be possible to distinguish the $`Y_n`$ and $`S_n`$ Wilson lines in $`S_{n\overline{n}}(\mu )`$. For instance if one chooses $`\mathrm{\Lambda }_{\mathrm{IR}}Q\lambda `$ then the usoft gluons give scaleless loop integrals and can be dropped, so that $`Y_n^{}S_n^{}S_{\overline{n}}Y_{\overline{n}}S_n^{}S_{\overline{n}}`$. If instead one chooses $`\mathrm{\Lambda }_{\mathrm{IR}}Q\lambda ^2`$ then the soft gluons give scaleless loop integrals (they simply act to pull-up the ultraviolet divergences in the usoft integrals to the hard scale amis3 ; hms1 ), so the soft Wilson lines can be suppressed. This is why one only finds $`S_n^{}S_{\overline{n}}`$ for this operator in the literature. For typical regulator choices the other gluons are simply not required to reproduce the infrared structure of the full theory result. In the above example the factors of $`S`$ and $`Y`$ in the operator in Eq. (25) do not cancel. In the examples we will consider in sections III and IV there are several operators at leading order, however the factors of $`S`$ and $`Y`$ cancel in observable matrix elements. The collinear matrix elements of long-distance operators, such as those in Eq. (II.2), are the ones that have interpretations as structure functions or wave functions. In the next section we give the operator definitions for these functions that will be needed in the remainder of the paper.

### II.3 Non-perturbative Matrix Elements

Predictions for hadronic processes depend on universal matrix elements that are not computable in perturbation theory. For exclusive processes these include light-cone wavefunctions and form factors, while for inclusive processes they include parton distribution functions and fragmentation functions. In this section we define matrix elements in SCET that are needed for our examples. All the collinear operators considered here decouple from usoft gluons since they are local in the residual coordinate $`x`$ and because $`Y^{}(x)Y(x)=1`$. Thus

$`\overline{\xi }_{n,p^{}}W\mathrm{\Gamma }W^{}\xi _{n,p}=\overline{\xi }_{n,p^{}}^{(0)}W^{(0)}\mathrm{\Gamma }W^{(0)}\xi _{n,p}^{(0)},`$ (31)
and expressions with and without the $`(0)`$ superscript are equal. For convenience we will write the collinear fields without the superscript in the remainder of this section.

Consider first the light-cone wavefunctions. For the pion isotriplet $`\pi ^a`$, the wavefunction $`\varphi _\pi (x)`$ is conventionally defined by brodskylepage

$`\pi ^a(p)\left|\overline{\psi }(y)\gamma ^\mu \gamma ^5{\displaystyle \frac{\tau ^b}{\sqrt{2}}}Y(y,x)\psi (x)\right|0=if_\pi \delta ^{ab}p^\mu {\displaystyle _0^1}𝑑ze^{i[zpy+(1z)px]}\varphi _\pi (\mu ,z),`$ (32)
Here $`f_\pi 131\mathrm{MeV}`$ and the QCD field $`\psi `$ denotes the isospin doublet $`\{\psi ^{(u)},\psi ^{(d)}\}`$. The coordinates satisfy $`(yx)^2=0`$, which ensures that the path from $`y^\mu `$ to $`x^\mu `$ is along the light-cone, and $`Y(y,x)`$ is a Wilson line along this path. In SCET we require the matrix elements of highly energetic pions, which therefore have collinear constituents. Boosting the matrix element in Eq. (32) and letting $`y^\mu =y\overline{n}^\mu `$, $`x^\mu =x\overline{n}^\mu `$ we have

$`\pi _{n,p}^a\left|\overline{\xi }_{n,y}\mathrm{\Gamma }_\pi ^bW(y,x)\xi _{n,x}\right|0=if_\pi \overline{n}p\delta ^{ab}{\displaystyle _0^1}𝑑ze^{i\overline{n}p[zy+(1z)x]}\varphi _\pi (\mu ,z),`$ (33)
where $`\mathrm{\Gamma }_\pi ^b=\overline{n}/\gamma _5\tau ^b/\sqrt{2}`$, and $`\xi _{n,z}`$ is a collinear field with position space label $`z\overline{n}^\mu `$. For our purposes it is more useful to use the operator with momentum space labels bps

$`\pi _{n,p}^a\left|\overline{\xi }_{n,p_1}W\mathrm{\Gamma }_\pi ^b\delta (\omega \overline{𝒫}_+)W^{}\xi _{n,p_2}\right|0`$ $`=`$ $`{\displaystyle \frac{dy}{2\pi }e^{i\omega y}\pi _{n,p}^a\left|\overline{\xi }_{n,y}\mathrm{\Gamma }_\pi ^bW(y,y)\xi _{n,y}\right|0}`$
$`=`$ $`if_\pi \overline{n}p\delta ^{ab}{\displaystyle _0^1}𝑑z\delta [\omega (2z1)\overline{n}p]\varphi _\pi (\mu ,z),`$
where $`\overline{𝒫}_\pm =\overline{𝒫}^{}\pm \overline{𝒫}`$. In Eq. (II.3) the delta function fixes $`\omega `$ to the sum of labels picked out by the $`\overline{𝒫}_+`$ operator. The combination picked out by $`\overline{𝒫}_{}`$ is equivalent to ($`\overline{𝒫}`$) acting on the entire operator, and using Eq. (6) is fixed to the $`\overline{n}p`$ momentum of the pion state.

In some situations it is convenient to have delta functions which fix the labels of both $`W^{}\xi _{n,p}`$ and $`\overline{\xi }_{n,p}W`$. In this case a useful field is

$`\chi _{n,\omega }^{(i)}`$ $``$ $`\left[\delta (\omega \overline{𝒫})W_n^{}\xi _{n,p}^{(i)}\right].`$ (35)
Here $`i`$ is the flavor index and will be omitted if the flavor doublet field is implied. Note that unlike the $`p`$ in $`\xi _{n,p}`$ the label $`\omega `$ on $`\chi _{n,\omega }`$ is not summed over. A matrix element with $`\chi _{n,\omega }`$ fields is related to a matrix element like the one in Eq. (II.3) through

$`M_{n,p}\left|\overline{\chi }_{n,\omega }\mathrm{\Gamma }\chi _{n,\omega ^{}}\right|M_{n,p^{}}`$ $`=`$ $`2\delta (\omega _{}\overline{n}p_{})M_{n,p}\left|\overline{\xi }_{n,p_1}W\mathrm{\Gamma }\delta (\omega _+\overline{𝒫}_+)W^{}\xi _{n,p_2}\right|M_{n,p^{}},`$
where $`\omega _\pm =\omega \pm \omega ^{}`$ and $`p_{}=pp^{}`$. Thus, with the $`\chi `$ notation the momentum conserving delta functions become explicit. The factor of two appears from treating the $`\omega `$’s as continuous variables, and in the final results cancels with a factor of $`1/2`$ from a Jacobian.

For inclusive processes such as DIS it is the proton parton distribution functions for quarks of flavor $`i`$, $`f_{i/p}(z)`$, antiquarks $`\overline{f}_{i/p}(z)`$, and gluons, $`f_{g/p}(z)`$ which are needed. The standard coordinate space definitions Soper  are ($`y^\mu =y\overline{n}^\mu `$)

$`f_{i/p}(z)`$ $`=`$ $`{\displaystyle \frac{dy}{2\pi }e^{i\mathrm{\hspace{0.17em}2}z\overline{n}py}p\left|\overline{\psi }^{(i)}(y)Y(y,y)\overline{n}/\psi ^{(i)}(y)\right|p}|_{\text{spin avg.}},`$ (37)
$`f_{g/p}(z)`$ $`=`$ $`{\displaystyle \frac{2}{z\overline{n}p}}{\displaystyle \frac{dy}{2\pi }e^{i\mathrm{\hspace{0.17em}2}z\overline{n}py}\overline{n}_\mu \overline{n}^\nu p\left|G_a^{\mu \lambda }(y)Y^{ab}(y,y)G_{\lambda \nu }^b(y)\right|p}|_{\text{spin avg.}},`$
and $`\overline{f}_{i/p}(z)=f_{i/p}(z)`$. Here $`G_{\mu \lambda }^a(y)`$ is the gluon field strength, $`Y(y,y)`$ and $`Y^{ab}(y,y)`$ are path-ordered Wilson lines in the fundamental and adjoint representations, and $`|p`$ is the proton state with momentum $`p`$. In SCET these distribution functions can be defined by the matrix elements of collinear fields with collinear proton states

$`{\displaystyle \frac{1}{2}}{\displaystyle \underset{\mathrm{spin}}{}}p_n\left|\overline{\chi }_{n,\omega }^{(i)}\overline{n}/\chi _{n,\omega ^{}}^{(i)}\right|p_n`$ $`=`$ $`4\overline{n}p{\displaystyle _0^1}𝑑z\delta (\omega _{})\delta (\omega _+2z\overline{n}p)f_{i/p}(z)`$
$``$ $`4\overline{n}p{\displaystyle _0^1}𝑑z\delta (\omega _{})\delta (\omega _++2z\overline{n}p)\overline{f}_{i/p}(z),`$
$`{\displaystyle \frac{1}{2}}{\displaystyle \underset{\mathrm{spin}}{}}p_n\left|\mathrm{Tr}\left[B_{n,\omega }^\mu B_\mu ^{n,\omega ^{}}\right]\right|p_n`$ $`=`$ $`{\displaystyle \frac{\omega _+\overline{n}p}{2}}{\displaystyle _0^1}𝑑z\delta (\omega _{})\delta (\omega _+2z\overline{n}p)f_{g/p}(z),`$
where $`\omega _\pm =\omega \pm \omega ^{}`$, and $`B_{n,\omega }^\mu \overline{n}_\nu (𝒢_{n,\omega })^{\nu \mu }`$ with the collinear gauge invariant field strength

$`(𝒢_{n,\omega })^{\mu \nu }={\displaystyle \frac{i}{g}}\left[\delta (\omega \overline{𝒫})W^{}[i𝒟_n^\mu +gA_{n,q}^\mu ,i𝒟_n^\nu +gA_{n,q^{}}^\nu ]W\right].`$ (39)
Both operators in Eq. (II.3) are order $`\lambda ^2`$ since $`\xi _{n,\omega }B_{n,\omega }^{}\lambda `$. Note that the matrix element of a single operator ($`\overline{\chi }_{n,\omega }^{(i)}\overline{n}/\chi _{n,\omega ^{}}^{(i)}`$) contains both the quark and antiquark distributions. This is due to Eq. (3), from which we see that for $`\omega =\omega ^{}>0`$ ($`\omega =\omega ^{}<0`$) this operator reduces to the number operator for collinear quarks (antiquarks) with momentum $`\omega `$.

Processes other than DIS sometimes depend on more complicated distribution functions. In deeply virtual Compton scattering (DVCS) we will need to parameterize the matrix element of an operator between proton states with different momenta. In terms of QCD fields $`\psi `$ the nonforward parton distribution function (NFPDF) defined by Radyushkin in Eq.(4.1) of Ref. Radyushkin:1997ki  are (up to a trivial translation)

$`p^{},\sigma ^{}\left|\overline{\psi }^{(i)}(y)Y(y,y)\overline{n}/\psi ^{(i)}(y)\right|p,\sigma `$ (40)
$`=e(\sigma ^{},\sigma ){\displaystyle _0^1}𝑑z\left[e^{i\overline{n}p(2z\zeta )y}_\zeta ^{(i)}(z;t)e^{i\overline{n}p(2z\zeta )y}\overline{}_\zeta ^{(i)}(z;t)\right]`$
$`+h(\sigma ^{},\sigma ){\displaystyle _0^1}𝑑z\left[e^{i\overline{n}p(2z\zeta )y}𝒦_\zeta ^{(i)}(z;t)e^{i\overline{n}p(2z\zeta )y}\overline{𝒦}_\zeta ^{(i)}(z;t)\right],`$
where $`t=(pp^{})^2`$, and $`\zeta =1\overline{n}p^{}/\overline{n}p`$. Here $`e(\sigma ,\sigma ^{})`$ and $`h(\sigma ,\sigma ^{})`$ are matrix elements which respectively preserve or flip the proton spin. They are defined in terms of the proton spinors

$`e(\sigma ^{},\sigma )`$ $`=`$ $`\overline{u}(p^{},\sigma ^{})\overline{n}/u(p,\sigma ),h(\sigma ^{},\sigma )={\displaystyle \frac{1}{2m_p}}\overline{u}(p^{},\sigma ^{})[\overline{n}/,p/p/^{}]u(p,\sigma ),`$ (41)
where $`m_p`$ is the proton mass. The NFPDF for gluons is similarly given by Radyushkin:1997ki

$`\overline{n}_\mu \overline{n}^\nu p^{},\sigma ^{}\left|G_a^{\mu \lambda }(y)Y^{ab}(y,y)G_{\lambda \nu }^b(y)\right|p,\sigma `$ (42)
$`={\displaystyle \frac{\overline{n}p}{2}}e(\sigma ^{},\sigma ){\displaystyle _0^1}𝑑z\left[e^{i\overline{n}p(2z\zeta )y}+e^{i\overline{n}p(2z\zeta )y}\right]_\zeta ^g(z;t)`$
$`+{\displaystyle \frac{\overline{n}p}{2}}h(\sigma ^{},\sigma ){\displaystyle _0^1}𝑑z\left[e^{i\overline{n}p(2z\zeta )y}+e^{i\overline{n}p(2z\zeta )y}\right]𝒦_\zeta ^g(z;t).`$
In SCET the definition of the NFPDFs in terms of collinear fields are

$`p_n^{},\sigma ^{}\left|\overline{\chi }_{n,\omega }^{(i)}\overline{n}/\chi _{n,\omega ^{}}^{(i)}\right|p_n,\sigma =2\delta (\omega _{}+\overline{n}p\zeta ){\displaystyle _0^1}𝑑z`$ (43)
$`\times \{e(\sigma ^{},\sigma )[\delta (\omega _+(2z\zeta )\overline{n}p)_\zeta ^{(i)}(z;t)\delta (\omega _++(2z\zeta )\overline{n}p)\overline{}_\zeta ^{(i)}(z;t)]`$
$`+h(\sigma ^{},\sigma )[\delta (\omega _+(2z\zeta )\overline{n}p)𝒦_\zeta ^{(i)}(z;t)\delta (\omega _++(2z\zeta )\overline{n}p)\overline{𝒦}_\zeta ^{(i)}(z;t)]\},`$
$`p_n^{},\sigma ^{}\left|\mathrm{Tr}\left[B_{n,\omega }^\mu B_\mu ^{n,\omega ^{}}\right]\right|p_n,\sigma ={\displaystyle \frac{\overline{n}p}{2}}\delta (\omega _{}+\overline{n}p\zeta ){\displaystyle _0^1}𝑑z`$
$`\times \left\{e(\sigma ^{},\sigma )\delta (\omega _+(2z\zeta )\overline{n}p)_\zeta ^g(z;t)+h(\sigma ^{},\sigma )\delta (\omega _+(2z\zeta )\overline{n}p)𝒦_\zeta ^g(z,t)\right\},`$
where the spinors in $`e(\sigma ^{},\sigma )`$ and $`h(\sigma ^{},\sigma )`$ are two component effective theory spinors, so that $`uu_n`$ where $`n/u_n=0`$. Note that for $`p^{}p`$ both $`\zeta 0`$ and $`t0`$. In this limit the NFPDFs reduce to the standard PDFs:

$`\underset{pp^{}}{lim}_\zeta ^{(i)}`$ $`=`$ $`f_{i/p}(z),\underset{pp^{}}{lim}\overline{}_\zeta ^{(i)}=\overline{f}_{i/p}(z),\underset{pp^{}}{lim}_\zeta ^{(g)}=zf_{g/p}(z).`$ (44)
### II.4 Symmetries for collinear fields

In this section we discuss spin and discrete symmetry constraints on operators involving collinear fields.

The possible spin structures of currents with $`\xi _{n,p}`$ fields is restricted by the fact that they have only two components, $`n/\xi _{n,p}=0`$. The four most general spin structures for currents with two collinear particles moving in the same or opposite directions are

$`\overline{\xi }_{n,p^{}}\mathrm{\Gamma }_1\xi _{n,p}\mathrm{\Gamma }_1=\{\overline{n}/,\overline{n}/\gamma _5,\overline{n}/\gamma _{}^\mu \},`$
$`\overline{\xi }_{\overline{n},p^{}}\mathrm{\Gamma }_2\xi _{n,p}\mathrm{\Gamma }_2=\{1,\gamma _5,\gamma _{}^\mu \}.`$ (45)
Other choices for $`\mathrm{\Gamma }_1`$ and $`\mathrm{\Gamma }_2`$ either vanish between the fields or are related to those in Eq. (II.4). This result can be expressed in a compact way by the trace formulae

$`\overline{\xi }_{n,p^{}}\mathrm{\Gamma }\xi _{n,p},\mathrm{\Gamma }={\displaystyle \frac{\overline{n}/}{8}}\mathrm{Tr}[n/\mathrm{\Gamma }]{\displaystyle \frac{\overline{n}/\gamma _5}{8}}\mathrm{Tr}[n/\gamma _5\mathrm{\Gamma }]{\displaystyle \frac{\overline{n}/\gamma _{}^\mu }{8}}\mathrm{Tr}[n/\gamma _\mu ^{}\mathrm{\Gamma }],`$
$`\overline{\xi }_{\overline{n},p^{}}\mathrm{\Gamma }\xi _{n,p},\mathrm{\Gamma }={\displaystyle \frac{1}{8}}\mathrm{Tr}[n/\overline{n}/\mathrm{\Gamma }]+{\displaystyle \frac{\gamma _5}{8}}\mathrm{Tr}[\gamma _5n/\overline{n}/\mathrm{\Gamma }]+{\displaystyle \frac{\gamma _{}^\mu }{8}}\mathrm{Tr}[\gamma _\mu ^{}n/\overline{n}/\mathrm{\Gamma }],`$ (46)
which reduce a general $`\mathrm{\Gamma }`$ to a linear combination of the terms in Eq. (II.4). For instance, it implies that $`2i\overline{\xi }_n\sigma ^{\mu \nu }\xi _n=n^\nu \overline{\xi }_n\overline{n}/\gamma _{}^\mu \xi _nn^\mu \overline{\xi }_n\overline{n}/\gamma _{}^\nu \xi _n`$, and $`\overline{\xi }_{\overline{n}}\gamma _{}^\mu \gamma _5\xi _n=iϵ_{}^{\mu \nu }\overline{\xi }_{\overline{n}}\gamma _\nu ^{}\xi _n`$ where $`ϵ_{}^{\mu \nu }=ϵ^{\mu \nu \alpha \beta }\overline{n}_\alpha n_\beta /2`$. Furthermore, each of the two components of $`\xi _n`$ and also $`\xi _{\overline{n}}`$ can be chosen to be eigenstates of their helicity operators, $`h=\widehat{p}\stackrel{}{S}`$ with eigenvalues $`\pm 1/2`$. For these fields $`h`$ is equivalent to the chiral rotation, $`h=\gamma _5/2`$. The structures in Eq. (II.4) split into two classes depending on whether they conserve or flip the helicity

chiral even: $`\text{ }\overline{\xi }_{n,p^{}}\{\overline{n}/,\overline{n}/\gamma _5\}\xi _{n,p},\overline{\xi }_{\overline{n},p^{}}\gamma _{}^\mu \xi _{n,p},`$ (47)
chiral odd: $`\text{ }\overline{\xi }_{n,p^{}}\overline{n}/\gamma _{}^\mu \xi _{n,p},\overline{\xi }_{\overline{n},p^{}}\{1,\gamma _5\}\xi _{n,p}.`$
Since gluon interactions in QCD preserve helicity, integrating out hard QCD fluctuations results in effective theory operators with the same helicity structure as the original operators at leading order in $`\lambda `$.

The presence of labels on the effective theory fields makes their transformation properties under the discrete symmetries $`C`$, $`P`$, and $`T`$ slightly different than in QCD. For example under $`P`$ or $`T`$ we have $`n\overline{n}`$ so these transformations relate collinear fields for different directions. Under charge conjugation, parity, and time-reversal the collinear fields transform as

$`\begin{array}{cccc}& C^1\xi _{n,p}(x)C=\left[\overline{\xi }_{n,p}(x)𝒞\right]^T,\hfill & & C^1A_{n,p}^\mu (x)C=[A_{n,p}^\mu (x)]^T,\hfill \\ & P^1\xi _{n,p}(x)P=\gamma _0\xi _{\overline{n},\stackrel{~}{p}}(x_P),\hfill & & P^1A_{n,p}^\mu (x)P=g_{\mu \nu }A_{\overline{n},\stackrel{~}{p}}^\nu (x_P),\hfill \\ & T^1\xi _{n,p}(x)T=𝒯\xi _{\overline{n},\stackrel{~}{p}}(x_T),\hfill & & T^1A_{n,p}^\mu (x)T=g_{\mu \nu }A_{\overline{n},\stackrel{~}{p}}^\nu (x_T),\hfill \end{array}`$ (51)
where $`𝒞^1\gamma _\mu 𝒞=\gamma _\mu ^T`$ and $`𝒯=\gamma ^5𝒞`$, while if $`x^\mu =(x^+,x^{},x^{})`$ and $`p^\mu =(p^+,p^{},p^{})`$ then $`\stackrel{~}{p}^\mu (p^{},p^+,p_{})`$, $`x_P^\mu (x^{},x^+,\stackrel{}{x}^{})`$, and $`x_T^\mu (x^{},x^+,\stackrel{}{x}^{})`$. The transformation properties of $`W_n`$ can be worked out using Eq. (51), for instance $`C^1W_nC=[W_n^{}]^T`$.

The collinear effective Lagrangian (10) is invariant under the transformations in Eq. (51) (adding the $`\overline{n}n`$ terms). These symmetries also constrain the form of non-perturbative matrix elements. As an example, for a meson which is an eigenstate of $`C`$ one finds

$`M_n\left|\overline{\xi }_{n,p}W_n\overline{n}/\gamma _5\delta (\omega \overline{𝒫}_+)W_n^{}\xi _{n,p^{}}\right|0`$
$`=(1)^CM_n\left|(𝒞W_n^{}\xi _{n,p})^T\overline{n}/\gamma _5\delta (\omega \overline{𝒫}_+)(\overline{\xi }_{n,p^{}}W_n𝒞)^T\right|0`$
$`=(1)^CM_n\left|\overline{\xi }_{n,p^{}}W_n\overline{n}/\gamma _5\delta (\omega +\overline{𝒫}_+)W_n^{}\xi _{n,p}\right|0.`$ (52)
For the isotriplet pion state $`(1)^C=+1`$ so combining Eq. (II.4) with Eq. (II.3) gives

$`\pi _{n,p}^a\left|\overline{\xi }_{n,p_1}W\mathrm{\Gamma }_\pi ^b\delta (\omega \overline{𝒫}_+)W^{}\xi _{n,p_2}\right|0`$
$`=if_\pi \delta ^{ab}\overline{n}p{\displaystyle _0^1}\text{d}x\delta [\omega (2x1)\overline{n}p]\varphi _\pi (x)`$
$`=if_\pi \delta ^{ab}\overline{n}p{\displaystyle _0^1}\text{d}x\delta [\omega (2x1)\overline{n}p]\varphi _\pi (1x).`$ (53)
Together with Eq. (II.3), charge conjugation therefore implies that $`\varphi _\pi (1x)=\varphi _\pi (x)`$.

## III Exclusive Processes

### III.1 $`\pi `$-$`\gamma `$ Form Factor

The pion-photon form factor $`F_{\pi \gamma }(Q^2)`$ is perhaps the simplest setting for factorization since there is only one hadron in the external state. The form factor is measurable in single-tagged two photon $`e^{}e^{}e^{}e^{}\pi ^0`$ reactions. This process involves the scattering of a highly virtual photon and a quark-anti-quark constituent pair off an on-shell photon. The photon scatters the quark pair away from the incoming photon into a pion, so that $`\gamma ^{}\gamma \pi ^0`$. The matrix element for this transition defines the $`\pi `$-$`\gamma `$ form factor

$`\pi ^0(p_\pi )\left|J_\mu (0)\right|\gamma (p_\gamma ,ϵ)`$ $`=`$ $`ieϵ^\nu {\displaystyle d^4ze^{ip_\gamma z}\pi ^0(p_\pi )\left|\mathrm{T}J_\mu (0)J_\nu (z)\right|0}`$ (54)
$`=`$ $`ieF_{\pi \gamma }(Q^2)ϵ_{\mu \nu \rho \sigma }p_\pi ^\nu ϵ^\rho q^\sigma .`$
Here $`J_\mu =\overline{\psi }\widehat{𝒬}\gamma _\mu \psi `$ is the full theory electromagnetic current with isodoublet field $`\psi `$ and charge matrix $`\widehat{𝒬}=\tau _3/2+\mathrm{𝟏}/6`$, and $`q^2=Q^2\mathrm{\Lambda }_{\mathrm{QCD}}^2`$ where $`q=p_\pi p_\gamma `$ is the virtual photon momentum. It has been shown that the form factor can be written as a one-dimensional convolution of a hard coefficient with the light-cone pion wavefunction Pink . Here we show how this factorization takes place in the SCET.

In the Breit frame $`q^\mu =Q(n^\mu \overline{n}^\mu )/2`$, the real photon’s momentum is $`p_\gamma ^\mu =E\overline{n}^\mu Q\overline{n}^\mu /2`$, and the pion is made up of collinear particles with momenta $`\overline{n}p_iQ`$. The particles exchanged between the two currents in Eq. (54) have hard momenta and can be integrated out. At leading order in $`\lambda `$ the time ordered product of the two currents in Eq. (54) matches onto a single operator in the effective theory. For simplicity we restrict ourselves to the tensor and spin structures that are relevant when the meson is a pion 999Note that a pure glue operator would not have the same isospin as the pion state.

$`O`$ $`=`$ $`{\displaystyle \frac{i}{Q}}ϵ_{\mu \nu }^{}[\overline{\xi }_{n,p}W]\mathrm{\Gamma }C(\overline{𝒫},\overline{𝒫}^{},\mu )[W^{}\xi _{n,p^{}}],`$ (55)
where $`\xi _{n,p}`$ is an isodoublet collinear quark field, and $`2ϵ_{\mu \nu }^{}=ϵ_{\mu \nu \rho \beta }\overline{n}^\rho n^\beta `$. $`O_1`$ is of dimension two, just like the time-ordered product in Eq. (54), and a power of $`1/Q`$ is included to make $`C(\mu ,\overline{𝒫},\overline{𝒫}^{})`$ dimensionless. The time-ordered product in Eq. (54) is even under charge conjugation, so the operators in Eq. (55) must also be even. This implies $`C_{\pi \gamma }(\mu ,\overline{𝒫},\overline{𝒫}^{})=C_{\pi \gamma }(\mu ,\overline{𝒫}^{},\overline{𝒫})`$. The location of the $`W`$’s in Eq. (55) is fixed by gauge invariance, and $`\mathrm{\Gamma }`$ contains the spin and flavor structure

$`\mathrm{\Gamma }`$ $`=`$ $`\left(\overline{n}/\gamma _5\right)\left(3\sqrt{2}\widehat{𝒬}^2\right).`$ (56)
Since the offshellness of the collinear particles in the pion is $`p^2\mathrm{\Lambda }_{\mathrm{QCD}}^2`$ we can also integrate out offshell modes with $`p^2Q\mathrm{\Lambda }_{\mathrm{QCD}}`$ which come from soft-collinear interactions. For the collinear operators $`O_j`$, Eq. (II.2) implies that factors of the soft Wilson line $`S_n`$ are induced. However, the location is such that $`S_n^{}S_n=1`$, so no coupling to soft gluons occur at leading order. The coupling of the collinear fields to usoft gluons can be simplified with the field redefinitions in Eq. (13). As discussed in section II.1 this moves all couplings into $`O_j`$, and using $`Y_n^{}Y_n=1`$ gives

$`O`$ $`=`$ $`{\displaystyle \frac{i}{Q}}ϵ_{\mu \nu }^{}[\overline{\xi }_{n,p}^{(0)}W^{(0)}]\mathrm{\Gamma }C(\overline{𝒫},\overline{𝒫}^{},\mu )[W^{(0)}\xi _{n,p^{}}^{(0)}].`$ (57)
Thus, usoft gluons also decouple.

In the Breit frame the pion momentum satisfies $`p_\pi ^\mu =E_\pi n^\mu +𝒪(\lambda )`$, and comparing Eq. (54) with the SCET matrix element $`i\pi _{n,p_\pi }^0|O_{\pi \gamma }|0`$, gives

$`{\displaystyle \frac{Q^2}{2}}F_{\pi \gamma }(Q^2)`$ $`=`$ $`{\displaystyle \frac{i}{Q}}\pi _n^0\left|\overline{\xi }_{n,p}^{(0)}W^{(0)}\mathrm{\Gamma }C(\overline{𝒫},\overline{𝒫}^{},\mu )W^{(0)}\xi _{n,p^{}}^{(0)}\right|0.`$ (58)
Defining $`\overline{𝒫}_\pm =\overline{𝒫}^{}\pm \overline{𝒫}`$, the operator $`\overline{𝒫}_{}`$ is related to $`\overline{𝒫}`$ acting from the outside on the fields. Using Eq. (6) it can therefore be set equal to the momentum label of the state, $`\overline{𝒫}_{}=\overline{n}p_\pi =Q`$. Suppressing this dependence we write $`C(\overline{𝒫},\overline{𝒫}^{},\mu )=C_1(\overline{𝒫}_+,\mu )`$ leaving

$`F_{\pi \gamma }(Q^2)`$ $`=`$ $`{\displaystyle \frac{2i}{Q^3}}\pi _n^0\left|\overline{\xi }_{n,p}^{(0)}W^{(0)}\mathrm{\Gamma }C_1(\overline{𝒫}_+,\mu )W^{(0)}\xi _{n,p}^{(0)}\right|0`$ (59)
$`=`$ $`{\displaystyle \frac{2i}{Q^3}}{\displaystyle 𝑑\omega C_1(\omega ,\mu )\pi _n^0\left|\overline{\xi }_{n,p}^{(0)}W^{(0)}\mathrm{\Gamma }\delta (\omega \overline{𝒫}_+)W^{(0)}\xi _{n,p}^{(0)}\right|0}.`$
Using Eq. (II.3) the remaining matrix identity in Eq. (59) can be written in terms of the light-cone pion wavefunction

$`F_{\pi \gamma }(Q^2)`$ $`=`$ $`{\displaystyle \frac{2f_\pi }{Q^2}}{\displaystyle 𝑑\omega _0^1𝑑x\delta (\omega (2x1)2E_\pi )C_1(\omega ,\mu )\varphi _\pi (x,\mu )}`$ (60)
$`=`$ $`{\displaystyle \frac{2f_\pi }{Q^2}}{\displaystyle _0^1}𝑑xC_1((2x1)Q,\mu )\varphi _\pi (x,\mu ).`$
This is the final result and is valid to leading order in $`\lambda `$ and all orders in $`\alpha _s`$. From Eq. (II.4) charge conjugation implies that $`\varphi _\pi (x)=\varphi _\pi (1x)`$ and $`C_1(\omega )=C_1(\omega )`$. Eq. (60) agrees with the Brodsky-Lepage brodskylepage  result that the form factor can be written as the convolution of a short distance function with the light-cone pion wavefunction. The SCET formalism gives a concise derivation of this result and defines the short distance function in terms of the Wilson coefficient of an effective theory operator.

As an illustrative example consider the tree level matching onto $`C`$ illustrated in Fig. 1. Since the location of the $`W`$’s in $`O`$ are fixed by gauge invariance, $`C(\mu ,\overline{𝒫},\overline{𝒫}^{})`$ can be determined by matching with $`W=1`$. Expanding the full theory graphs to leading order gives

$`i(\text{Fig.}1)={\displaystyle \frac{ie}{2}}ϵ_{\mu \nu \rho \beta }ϵ^\nu \overline{n}^\rho n^\beta \left({\displaystyle \frac{\overline{n}/}{2}}\gamma _5\right)\left(\widehat{𝒬}^2\right)\left({\displaystyle \frac{1}{\overline{n}p}}{\displaystyle \frac{1}{\overline{n}p^{}}}\right),`$ (61)
where we have dropped isosinglet terms, contributions with opposite parity to the pion, as well as those proportional to $`n/\gamma _5`$ since $`n/\xi _{n,p}=0`$. Comparing Eq. (61) to Eq. (55) gives

$`C(\mu ,\overline{𝒫},\overline{𝒫}^{})={\displaystyle \frac{Q}{6\sqrt{2}}}\left({\displaystyle \frac{1}{\overline{𝒫}^{}}}{\displaystyle \frac{1}{\overline{𝒫}}}\right)+𝒪\left(\alpha _s(Q)\right),`$ (62)
so that

$`C_1(\mu ,\omega =(2x1)Q)={\displaystyle \frac{1}{6\sqrt{2}}}\left({\displaystyle \frac{1}{x}}+{\displaystyle \frac{1}{1x}}\right)+𝒪\left(\alpha _s(Q)\right).`$ (63)
This result is again in agreement with Ref. brodskylepage , and the order $`\alpha _s(Q)`$ corrections to this Wilson coefficient can be read off from the results in Ref. delAguila:1981nk ; Braaten:yp . An identical analysis applies for operators with different spin structures such as the ones contributing to $`\gamma ^{}\gamma \rho ^0`$.

### III.2 The large $`Q^2`$ meson form factor

Another example of an exclusive process which can be treated in the effective theory is the classic case of the electromagnetic pion form factor at large space-like momentum transfer. For generality we consider in this section the electromagnetic form factors for arbitrary mesons (pseudoscalar $`P`$ or vector $`V`$), defined as

$`P^{}(p^{})\left|J_\mu \right|P(p)=F_P(Q^2)(p_\mu +p_\mu ^{}),`$ (64)
$`V^{}(p^{},\epsilon ^{})\left|J_\mu \right|P(p)=G(Q^2)iϵ_{\mu \nu \alpha \beta }p^\nu p^\alpha \epsilon ^\beta ,`$
$`V^{}(p^{},\epsilon ^{})\left|J_\mu \right|V(p,\epsilon )=F_1(Q^2)(\epsilon ^{}\epsilon )(p+p^{})_\mu +F_2(Q^2)[(\epsilon ^{}p)\epsilon _\mu +(\epsilon p^{})\epsilon _\mu ^{}],`$
where $`q^2=Q^2,q=pp^{}`$. For simplicity we suppress the dependence of the form factors on the isospin of the two mesons. We will restrict ourselves in the following to the case of hadrons made up only of $`u,d`$ quarks. The electromagnetic current is defined as usual by $`J_\mu =\overline{q}\widehat{𝒬}\gamma _\mu q`$, with charge matrix $`\widehat{𝒬}=\text{diag}(2/3,1/3)`$, which can be written in terms of the up and down quark charges as $`\widehat{𝒬}=(Q_uQ_d)\tau _3/2+(Q_u+Q_d)\mathrm{𝟏}/2`$.

We will be interested in the asymptotic form of the form factors in the region with $`Q^2m_P^2`$, where it can be expanded in a power series in $`1/Q^2`$ brodskylepage . It is convenient to work in the Breit frame, where the momentum transfer has the light-cone components $`q=(q^+,q^{},q_{})=(Q,Q,0)`$. In this frame the meson momenta are $`p=(Q,m_P^2/Q,0)`$, $`p^{}=(m_P^2/Q,Q,0)`$, so the partons in the incoming/outgoing meson are collinear along the $`\overline{n}_\mu `$/$`n_\mu `$ direction.

The electromagnetic current in Eq. (64) is matched in the effective theory onto the most general combination of operators constructed from collinear fields which are compatible with collinear gauge invariance. Operators such as the dimension three current

$`[\overline{\xi }_{\overline{n}}W_{\overline{n}}]\mathrm{\Gamma }C(\mu ,𝒫^{},\overline{𝒫})[W_n^{}\xi _n],`$ (65)
can contribute, but only overlap with the asymmetric meson states with one energetic collinear quark and one usoft or soft quark. Often this overlap is referred to as the tail of the wavefunction contribution or the Feynman mechanism of generating the form factor brodskyfeynman ; isgurfeynman . There are other operators with significant overlap with more symmetric meson states (where all the constituents are allowed to be energetic). The leading such operators have the form 101010There are also gluon operators that can contribute when one or more of the mesons is a neutral isosinglet, however for simplicity these are not discussed here.

$`{\displaystyle \frac{1}{Q^3}}\left[\overline{\xi }_{n,p_1}W_n\mathrm{\Gamma }W_{\overline{n}}^{}\xi _{\overline{n},p_2}\right]C(\mu ,\overline{𝒫},\overline{𝒫}^{},𝒫,𝒫^{})\left[\overline{\xi }_{\overline{n},p_3}W_{\overline{n}}\mathrm{\Gamma }^{}W_n^{}\xi _{n,p_4}\right],`$ (66)
with $`C`$ a dimensionless Wilson coefficient. As usual, collinear gauge invariance is enforced by the location of the $`W`$’s in Eqs. (65) and (66). There is some argument about the relative size of Eqs. (65) and (66) in the literature brodskyfeynman ; isgurfeynman . Often it is argued that the tail of the wavefunction is suppressed by an extra $`\mathrm{\Lambda }_{\mathrm{QCD}}^2/Q^2`$ brodskyfeynman , in which case the operator in Eq. (66) dominates by two powers of $`Q`$. An analysis of the tail of wavefunction contributions has not yet been performed in the effective theory framework. Therefore, we choose to ignore the operator in Eq. (65), and below only analyze the operator in Eq. (66). We emphasize that we do not claim to have shown that this is justified by the effective theory power counting.

There are two different structures possible for the operator in Eq. (66), and we write the general matching for the electromagnetic current as

$`J^\nu {\displaystyle \frac{1}{Q^3}}{\displaystyle 𝑑\omega _j\left[C_0(\mu ,\omega _j)𝒥_0^\nu (\omega _j,\mu )+C_8(\mu ,\omega _j)𝒥_8^\nu (\mu ,\omega _j)\right]},`$ (67)
where $`j=1,2,3,4`$. The SCET currents are dimension-6 operators

$`𝒥_0^\nu `$ $`=`$ $`\overline{\chi }_{n,\omega _1}\mathrm{\Gamma }\chi _{\overline{n},\omega _2}\overline{\chi }_{\overline{n},\omega _3}\mathrm{\Gamma }^{}\chi _{n,\omega _4}(\mathrm{\Gamma }\mathrm{\Gamma }^{},\omega _{1,2}\omega _{4,3}),`$ (68)
$`𝒥_8^\nu `$ $`=`$ $`\overline{\chi }_{n,\omega _1}\mathrm{\Gamma }T^a\chi _{\overline{n},\omega _2}\overline{\chi }_{\overline{n},\omega _3}\mathrm{\Gamma }^{}T^a\chi _{n,\omega _4}(\mathrm{\Gamma }\mathrm{\Gamma }^{},\omega _{1,2}\omega _{4,3}),`$
where the $`\chi `$ fields are defined in Eq. (35). In terms of the charge matrix $`\widehat{𝒬}`$, the spin and flavor structure is

$`\mathrm{\Gamma }\mathrm{\Gamma }^{}=(n^\nu +\overline{n}^\nu )\left(\gamma _\alpha ^{}\widehat{𝒬}\gamma _{}^\alpha \mathrm{𝟏}\right).`$ (69)
The Wilson coefficients $`C_{0,8}`$ can be computed in a power series in $`\alpha _s(Q)`$. They are functions of $`\mu `$, $`Q`$, and the $`\omega _j`$ which are the sum of momentum labels for gauge invariant products of collinear fields in the SCET currents.

The currents operators in Eq. (68) are the most general allowed operators which are gauge invariant, transform the same way as $`J^\mu `$ under charge conjugation and satisfy current and helicity conservation. To see how these properties constrain the form of the allowed operators, we begin by noting that Eq. (II.4) implies that $`\mathrm{\Gamma },\mathrm{\Gamma }^{}=\{1,\gamma _5,\gamma _{}^\mu \}`$ are the most general allowed spin structures. For massless quarks the electromagnetic and QCD couplings preserve helicity, whereas $`\overline{\xi }_n\{1,\gamma _5\}\xi _{\overline{n}}`$ cause the helicity to flip. Thus, only the structure $`\overline{\xi }_n\gamma _{}^\mu \xi _{\overline{n}}`$ is allowed. Current conservation $`q^\nu 𝒥_\nu =0`$, together with $`q^\nu =Q(\overline{n}^\nu n^\nu )/2`$ implies $`𝒥_\nu (n_\nu +\overline{n}_\nu )`$. Under charge conjugation $`J^\mu J^\mu `$ so the same must be true for the SCET currents. In the current operators, charge conjugation switches $`\omega _1\omega _4`$, $`\omega _2\omega _3`$, and $`\mathrm{\Gamma }\mathrm{\Gamma }^{}`$, as can be seen from Eq. (51). Thus, the second term in $`𝒥_{0,8}^\nu `$ is required to make these operators odd under charge conjugation. The operators $`𝒥_{0,8}^\nu `$ and the full electromagnetic current are invariant under a combined $`PT`$ transformation. This requires that the Wilson coefficients are real.

The operators $`𝒥_{0,8}`$ are responsible for the $`P_{\overline{n}}P_n`$ transition, while the reverse transition $`P_nP_{\overline{n}}`$ is described by similar operators with $`\overline{n}n`$. Parity invariance requires the Wilson coefficients of these operators to be identical to $`C_{0,8}(\omega _i)`$. Demanding hermiticity of the electromagnetic current in the effective theory then gives the relation $`C_{0,8}(\omega _1,\omega _2,\omega _3,\omega _4)=C_{0,8}^{}(\omega _2,\omega _1,\omega _4,\omega _3)`$. Since the coefficients are real they must therefore satisfy $`C_{0,8}(\omega _1,\omega _2,\omega _3,\omega _4)=C_{0,8}(\omega _2,\omega _1,\omega _4,\omega _3)`$.

To compute the matrix elements in the effective theory, it is convenient to Fierz transform the four-quark operators in Eq. (68). This gives

$`C_0𝒥_0+C_8𝒥_8=(n^\nu +\overline{n}^\nu ){\displaystyle \underset{j=1}{\overset{4}{}}}C_j𝒥_j,`$ (70)
where

$`𝒥_j`$ $`=`$ $`\left[\overline{\chi }_{n,\omega _1}\mathrm{\Gamma }_j\chi _{n,\omega _4}\right]\left[\overline{\chi }_{\overline{n},\omega _3}\mathrm{\Gamma }_j^{}\chi _{\overline{n},\omega _2}\right].`$ (71)
The spin, flavor, and color structures are

$`\mathrm{\Gamma }_1\mathrm{\Gamma }_1^{}`$ $`=`$ $`{\displaystyle \frac{1}{4}}(Q_uQ_d)iϵ^{3bc}(\tau ^b\tau ^c)[\overline{n}/n/+\overline{n}/\gamma _5n/\gamma _5],`$
$`\mathrm{\Gamma }_2\mathrm{\Gamma }_2^{}`$ $`=`$ $`{\displaystyle \frac{1}{4}}\left[(Q_u+Q_d)(\mathrm{𝟏}\mathrm{𝟏}+\tau ^a\tau ^a)+(Q_uQ_d)(\mathrm{𝟏}\tau ^3+\tau ^3\mathrm{𝟏})\right]`$ (72)
$`\times [\overline{n}/n/+\overline{n}/\gamma _5n/\gamma ^5],`$
while $`\mathrm{\Gamma }_{3,4}=T^a\mathrm{\Gamma }_{1,2}`$ and $`\mathrm{\Gamma }_{3,4}^{}=T^a\mathrm{\Gamma }_{1,2}^{}`$. The new Wilson coefficients are

$`C_1(\mu ,\omega _j)`$ $`=`$ $`[{\displaystyle \frac{1}{8}}(1{\displaystyle \frac{1}{N_c^2}})C_8(\mu ,\omega _j)+{\displaystyle \frac{1}{4N_c}}C_0(\mu ,\omega _j)]+(\omega _{1,2}\omega _{4,3}),`$
$`C_2(\mu ,\omega _j)`$ $`=`$ $`[{\displaystyle \frac{1}{8}}(1{\displaystyle \frac{1}{N_c^2}})C_8(\mu ,\omega _j)+{\displaystyle \frac{1}{4N_c}}C_0(\mu ,\omega _j)](\omega _{1,2}\omega _{4,3}),`$ (73)
with similar relations for $`C_{3,4}`$ which are also in terms of $`C_{0,8}`$.

A few general predictions follow from the form of the operators in Eq. (71).111111These predictions depend on the dominance of the operators in Eq. (66) over those in Eq. (65). For mesons with spin, only helicity conserving form factors appear, and furthermore no off-diagonal (e.g., $`PV`$) matrix elements are present at leading order in $`1/Q^2`$. These results agree with Ref. Pink . We also see that the form factors between arbitrary meson states are determined at leading power by only two hard coefficients, $`C_0`$ and $`C_8`$.

Now consider what factorization tells us about the matrix element of the operators in Eq. (71). For the decoupling of usoft and soft gluons we will follow section II.1. Integrating out offshell modes with $`p^2Q\mathrm{\Lambda }_{\mathrm{QCD}}`$ induces soft Wilson lines $`S_n`$ and $`S_{\overline{n}}`$, while the field redefinitions in Eq. (13) make all couplings to usoft gluons explicit in the operators. Together these give

$`𝒥_j^\nu `$ $`=`$ $`\left[\overline{\chi }_{n,\omega _1}^{(0)}Y_n^{}S_n^{}\mathrm{\Gamma }_jS_nY_n\chi _{n,\omega _4}^{(0)}\right]\left[\overline{\chi }_{\overline{n},\omega _3}^{(0)}Y_{\overline{n}}^{}S_{\overline{n}}^{}\mathrm{\Gamma }_j^{}S_{\overline{n}}Y_{\overline{n}}\chi _{\overline{n},\omega _2}^{(0)}\right].`$ (74)
Consider first the color singlet currents $`j=1,2`$. Here the $`Y`$’s and $`S`$’s all cancel using unitarity of the Wilson lines. Since the $`A_{n,q}^\mu `$ and $`A_{\overline{n},q}^\mu `$ gluons only interact with fields in the $`n`$ and $`\overline{n}`$ directions respectively, collinear gluons are not exchanged between the $`n`$ and $`\overline{n}`$ quark bilinears. Thus, the matrix element between states with particles moving in the $`n`$ and $`\overline{n}`$ directions factors

$`n|𝒥_{1,2}|\overline{n}=n|\overline{\chi }_{n,\omega _1}^{(0)}\mathrm{\Gamma }_{1,2}\chi _{n,\omega _4}^{(0)}|00|\overline{\chi }_{\overline{n},\omega _3}^{(0)}\mathrm{\Gamma }_{1,2}^{}\chi _{\overline{n},\omega _2}^{(0)}|\overline{n}.`$ (75)
Next consider the currents $`𝒥_{3,4}`$, which have color structure $`T^aT^a`$ in $`\mathrm{\Gamma }_j\mathrm{\Gamma }_j^{}`$. In this case the usoft and soft gluons do not cancel, but can all be moved into one quark bilinear using the color identity $`Y_n^{}S_n^{}T^aS_nY_nY_{\overline{n}}^{}S_{\overline{n}}^{}T^aS_{\overline{n}}Y_{\overline{n}}=T^aY_{\overline{n}}^{}S_{\overline{n}}^{}S_nY_nT^aY_n^{}S_n^{}S_{\overline{n}}Y_{\overline{n}}`$. After this rearrangement it is clear that the (u)soft gluons and $`A_{n,q}^\mu `$ and $`A_{\overline{n},q}^\mu `$ gluons only interact with the fields in one of the quark bilinears. Thus, the matrix element $`n|𝒥_{3,4}|\overline{n}`$ factors, similar to Eq. (75). For color singlet states, however, the matrix element of an octet operator vanishes identically since

$`n|\overline{\chi }_{n,\omega _1}^{(0)}T^a\chi _{n,\omega _4}^{(0)}|0=0.`$ (76)
Thus, the effective theory currents $`𝒥_{3,4}`$ do not contribute to the form factors at any order in perturbation theory.

Eq. (75) shows that for arbitrary meson states factorization occurs. It remains to show that the matrix elements in Eq. (75) are given by a two-dimensional convolution with the light-cone meson wavefunctions. To do this we consider the simple example of the $`0^+0^+`$ form factor for the charged pion. It should be obvious that the same steps go through for other meson states.

The symmetry of the pion wavefunction $`\varphi _\pi (x)`$ under charge conjugation ($`x1x`$) implies that only the $`𝒥_1`$ current contributes. Thus,

$`F_{\pi ^\pm }(Q^2)={\displaystyle \frac{2}{Q^4}}{\displaystyle 𝑑\omega _jC_1(\mu ,\omega _j)\pi _n^\pm (p^{})\left|\overline{\chi }_{n,\omega _1}^{(0)}\mathrm{\Gamma }_1\chi _{n,\omega _4}^{(0)}\right|00\left|\overline{\chi }_{\overline{n},\omega _3}^{(0)}\mathrm{\Gamma }_1^{}\chi _{\overline{n},\omega _2}^{(0)}\right|\pi _{\overline{n}}^\pm (p)}.`$ (77)
The required matrix elements can be obtained from Eqs. (II.3) and (II.3) with $`|\pi ^\pm =(\pi ^1\pm i\pi ^2)/\sqrt{2}`$. The momentum conserving delta functions fix $`\omega _1\omega _4=\overline{n}p^{}=Q`$ and $`\omega _2\omega _3=\overline{n}p=Q`$, while the $`\omega =\omega _1+\omega _4`$ and $`\omega ^{}=\omega _3+\omega _2`$ integrations can be done with the delta functions. This leaves

$`F_{\pi ^\pm }`$ $`=`$ $`\pm (Q_uQ_d){\displaystyle \frac{f_\pi ^2}{Q^2}}{\displaystyle _0^1}\text{d}x{\displaystyle _0^1}\text{d}yT_1(x,y,\mu )\varphi _\pi (x,\mu )\varphi _\pi (y,\mu ),`$ (78)
where $`T_1(x,y)`$ is defined in terms of $`C_1(\omega _1,\omega _2,\omega _3,\omega _4)`$ as

$`T_1(x,y)=C_1(xQ,yQ,(y1)Q,(x1)Q).`$ (79)
The coefficients $`C_j(\mu ,\omega _j)`$, and therefore also $`T_j(x,y)`$, can be obtained at the scale $`\mu =Q`$ by a matching calculation, as illustrated in Fig. 2. For this purpose, it is sufficient to compute the matrix element of the currents with free collinear quarks. To lowest order in $`\alpha _s(Q)`$, only $`C_8(\omega _j,\mu =Q)`$ is nonvanishing

$`C_0(\omega _j,\mu =Q)=0,C_8(\omega _j,\mu =Q)=4\pi \alpha _s(Q){\displaystyle \frac{Q^2}{\omega _3\omega _4}}.`$ (80)
This implies

$`T_1(x,y,\mu =Q)={\displaystyle \frac{4\pi \alpha _s(Q)}{9}}\left[{\displaystyle \frac{1}{xy}}+{\displaystyle \frac{1}{(1x)(1y)}}\right],`$ (81)
Using the asymptotic light-cone pion wavefunction $`\varphi _\pi (x)=6x(1x)`$ we find agreement with Ref. brodskylepage ,

$`F_{\pi ^\pm }(Q^2)=\pm {\displaystyle \frac{8\pi f_\pi ^2\alpha _s(Q)}{9Q^2}}\left[{\displaystyle _0^1}\text{d}x{\displaystyle \frac{\varphi _\pi (x)}{x}}\right]^2\pm {\displaystyle \frac{8\pi f_\pi ^2\alpha _s(Q)}{Q^2}}.`$ (82)
The order $`\alpha _s^2(Q)`$ corrections to Eq. (80) can be found in Refs. Field:wx ; Dittes:aw ; Braaten:yy .

## IV Inclusive Processes

### IV.1 Deep inelastic scattering

DIS is a process which is both simple and rich in physics. As such it provides an ideal introduction to inclusive factorization in QCD, which we study from an effective field theory point of view in this section. The aim is to prove that to all orders in $`\alpha _s`$ and leading order in $`\lambda `$ the DIS forward scattering amplitude can be written as an integral over hard coefficients times the parton distribution functions. This is done by matching onto local operators in SCET. The properties of SCET are used to show that matrix elements of the leading local operator can be written as a convolution of a hard coefficient with the parton distribution functions for the proton.

The first step is to understand the kinematics of the process. The hard scale $`Q^2=q^2`$ is set by the invariant mass of the photon, and $`x=Q^2/(2pq)`$ is the Bjorken variable. In the Breit frame the momentum of the virtual photon is $`q^\mu =Q(\overline{n}^\mu n^\mu )/2`$, and the incoming proton momentum is $`p^\mu =n^\mu \overline{n}p/2+\overline{n}^\mu m_p^2/(2\overline{n}p)n^\mu Q/(2x)+\overline{n}^\mu xm_p^2/(2Q)`$ up to terms $`m_p^2/Q^2`$, where $`m_p`$ is the proton mass. By momentum conservation the final state momentum is $`P_X^\mu =q^\mu +p^\mu `$, which gives an invariant mass $`P_X^2=(Q^2/x)(1x)+m_p^2`$. Values $`1x\mathrm{\Lambda }_{\mathrm{QCD}}/Q`$ correspond to the endpoint region where the particles in $`X`$ are collimated into a jet, while values $`1x\mathrm{\Lambda }_{\mathrm{QCD}}^2/Q^2`$ correspond to the resonance region. We will consider the standard OPE region where $`1x\mathrm{\Lambda }_{\mathrm{QCD}}/Q`$ so that the final state has virtuality of order $`Q^2`$ and can be integrated out. In contrast, although the incoming proton has a large momentum component in the $`n^\mu `$ direction it has a small invariant mass $`p^2=m_p^2\mathrm{\Lambda }_{\mathrm{QCD}}^2`$, and therefore is described by collinear fields in the effective theory.

Consider the spin-averaged cross section for DIS which can be written as

$`d\sigma ={\displaystyle \frac{d^3𝐤^{}}{2|𝐤^{}|(2\pi )^3}}{\displaystyle \frac{\pi e^4}{sQ^4}}L^{\mu \nu }(k,k^{})W_{\mu \nu }(p,q),`$ (83)
where $`k`$ and $`k^{}`$ are the incoming and outgoing lepton momenta with $`q=k^{}k`$, $`L^{\mu \nu }`$ is the lepton tensor, and $`s=(p+k)^2`$. The hadronic tensor $`W^{\mu \nu }`$ can be related to the imaginary part of the DIS forward scattering amplitude:

$`W_{\mu \nu }(p,q)`$ $`=`$ $`{\displaystyle \frac{1}{\pi }}\mathrm{Im}T_{\mu \nu }(p,q),`$ (84)
$`T_{\mu \nu }(p,q)`$ $`=`$ $`{\displaystyle \frac{1}{2}}{\displaystyle \underset{\mathrm{spin}}{}}p\left|\widehat{T}_{\mu \nu }(q)\right|p,\widehat{T}_{\mu \nu }(q)=i{\displaystyle d^4ze^{iqz}\mathrm{T}[J_\mu (z)J_\nu (0)]},`$
where for an electromagnetic current $`J_\mu `$ we can write

$`T_{\mu \nu }(p,q)=\left(g_{\mu \nu }+{\displaystyle \frac{q_\mu q_\nu }{q^2}}\right)T_1(x,Q^2)+\left(p_\mu +{\displaystyle \frac{q_\mu }{2x}}\right)\left(p_\nu +{\displaystyle \frac{q_\nu }{2x}}\right)T_2(x,Q^2).`$ (85)
As explained above, the intermediate hadronic state has invariant mass $`P_X^2Q^2`$. Therefore, one can perform an OPE and match $`\widehat{T}^{\mu \nu }(q)`$ onto operators in SCET. All fields in the resulting operators are evaluated at the same residual space time point, however, the presence of Wilson lines and label momenta make the operators nonlocal along a particular light cone direction. These nonlocal operators sum the infinite set of purely local operators of a given twist, however this is built into the formalism automatically.

To match we write down the most general leading operator in SCET which contains collinear fields moving in the $`n^\mu `$ direction, and enforce the condition from current conservation $`q^\mu \widehat{T}_{\mu \nu }=0`$. This leads to

$`\widehat{T}^{\mu \nu }`$ $``$ $`{\displaystyle \frac{g_{}^{\mu \nu }}{Q}}\left({\displaystyle \underset{i}{}}O_1^{(i)}+{\displaystyle \frac{O_1^g}{Q}}\right)+{\displaystyle \frac{(n^\mu +\overline{n}^\mu )(n^\nu +\overline{n}^\nu )}{Q}}\left({\displaystyle \underset{i}{}}O_2^{(i)}+{\displaystyle \frac{O_2^g}{Q}}\right),`$ (86)
where

$`O_j^{(i)}`$ $`=`$ $`\overline{\xi }_{n,p^{}}^{(i)}W{\displaystyle \frac{\overline{n}/}{2}}C_j^{(i)}(\overline{𝒫}_+,\overline{𝒫}_{})W^{}\xi _{n,p}^{(i)},`$
$`O_j^g`$ $`=`$ $`\overline{n}_\mu \overline{n}_\nu \mathrm{tr}\left[W^{}(G_n)^{\mu \lambda }WC_j^g(\overline{𝒫}_+,\overline{𝒫}_{})W^{}(G_n)_\lambda ^\nu W\right],`$ (87)
where $`i`$ labels the flavor of the fermions and $`igG_n^{\mu \lambda }=[i𝒟_n^\mu +gA_{n,q}^\mu ,i𝒟_n^\lambda +gA_{n,q^{}}^\lambda ]`$. The Wilson coefficients are dimensionless functions of $`\overline{𝒫}_+`$, $`\overline{𝒫}_{}`$, $`Q`$, and $`\mu `$. As in previous sections we can separate the hard coefficients from the long distance operators by introducing trivial convolutions. This gives

$`O_j^{(i)}`$ $`=`$ $`{\displaystyle 𝑑\omega _1𝑑\omega _2C_j^{(i)}(\omega _+,\omega _{})\left[\overline{\chi }_{n,\omega _1}^{(i)}\frac{\overline{n}/}{2}\chi _{n,\omega _2}^{(i)}\right]},`$
$`O_j^g`$ $`=`$ $`{\displaystyle 𝑑\omega _1𝑑\omega _2C_j^g(\omega _+,\omega _{})\mathrm{tr}\left[B_{n,\omega _1}^\mu B_\mu ^{n,\omega _2}\right]},`$ (88)
where $`\omega _\pm =\omega _1\pm \omega _2`$, and $`B_{n,\omega }^\mu \overline{n}_\nu (𝒢_{n,\omega })^{\nu \mu }`$ with $`(𝒢_{n,\omega })^{\mu \lambda }`$ defined in Eq. (39). Next we factor the coupling of usoft gluons from the collinear fields using the field redefinitions in Eq. (13). The operator $`O_j^{(i)}`$ has the structure in Eq. (31) so the $`Y`$’s cancel trivially, while for $`O_j^g`$ we find

$`B_n^\mu =Y_nB_n^{\mu (0)}Y_n^{},`$ (89)
and the factors of $`Y`$ cancel in the trace. It is easy to see that soft gluons also decouple using Eq. (II.2) or by noting that there is no non-trivial soft gauge invariant way of adding soft Wilson lines $`S_n`$ to $`O_j^{(i)}`$ or $`O_j^g`$. Under charge conjugation the full theory electromagnetic current $`J_\mu J_\mu `$ and therefore the operator $`\widehat{T}_{\mu \nu }\widehat{T}_{\mu \nu }`$. This implies relations for the effective theory Wilson coefficients since the operators $`O_j^{(i)}`$ must also respect this symmetry. Thus charge conjugation gives

$`{\displaystyle 𝑑\omega _1𝑑\omega _2C_j^{(i)}(\omega _+,\omega _{})\overline{\chi }_{n,\omega _1}^{(i)}\overline{n}/\chi _{n,\omega _2}^{(i)}}`$ $`\stackrel{\text{C}}{=}`$ $`{\displaystyle 𝑑\omega _1𝑑\omega _2C_j^{(i)}(\omega _+,\omega _{})\overline{\chi }_{n,\omega _2}^{(i)}\overline{n}/\chi _{n,\omega _1}^{(i)}}`$ (90)
$`=`$ $`{\displaystyle 𝑑\omega _1𝑑\omega _2\left[C_j^{(i)}(\omega _+,\omega _{})\right]\overline{\chi }_{n,\omega _1}^{(i)}\overline{n}/\chi _{n,\omega _2}^{(i)}}.`$
In the second line we changed variable $`\omega _1\omega _2`$ and $`\omega _2\omega _1`$ which takes $`\omega _+\omega _+`$ and $`\omega _{}\omega _{}`$. Thus, to all orders in perturbation theory $`C_j^{(i)}(\omega _+,\omega _{})=C_j^{(i)}(\omega _+,\omega _{})`$. This relates the Wilson coefficients for quarks and anti-quarks. Note that the above results are all independent of the collinear hadron on which DIS is performed.

Next we take the matrix element between proton states. Using the definitions of the nonperturbative matrix elements given in Section II.3, and picking out the coefficients of the tensor structures we find that the delta functions in Eq. (II.3) set $`\omega _+=\pm 2Q\xi /x`$ and $`\omega _{}=0`$. Since charge conjugation relates negative and positive values of $`\omega _+`$ only coefficients, $`C_j(\omega _+,0)`$, with positive $`\omega _+`$ are needed in the formulae for DIS. Therefore we define

$`H_j(z)C_j(2Qz,0,Q,\mu ),`$ (91)
where here we have made the dependence on $`Q`$ and $`\mu `$ explicit. Combining this with Eqs. (85), (86), (IV.1), and (II.3), gives the final result

$`T_1(x,Q^2)`$ $`=`$ $`{\displaystyle \frac{1}{x}}{\displaystyle _0^1}𝑑\xi \left\{H_1^{(i)}\left({\displaystyle \frac{\xi }{x}}\right)\left[f_{i/p}(\xi )+\overline{f}_{i/p}(\xi )\right]+{\displaystyle \frac{\xi }{2x}}H_1^g\left({\displaystyle \frac{\xi }{x}}\right)f_{g/p}(\xi )\right\},`$
$`T_2(x,Q^2)`$ $`=`$ $`{\displaystyle \frac{4x}{Q^2}}{\displaystyle _0^1}d\xi \{[4H_2^{(i)}\left({\displaystyle \frac{\xi }{x}}\right)H_1^{(i)}\left({\displaystyle \frac{\xi }{x}}\right)][f_{i/p}(\xi )+\overline{f}_{i/p}(\xi )]`$ (92)
$`+{\displaystyle \frac{\xi }{2x}}[4H_2^g\left({\displaystyle \frac{\xi }{x}}\right)H_1^g\left({\displaystyle \frac{\xi }{x}}\right)]f_{g/p}(\xi )\},`$
where a sum over $`i`$ is implicit. The hadronic tensor components $`W_{1,2}(x,Q^2)=\mathrm{Im}T_{1,2}(x,Q^2)/\pi `$ and therefore are determined by the imaginary part of the Wilson coefficients. The Wilson coefficients are dimensionless and therefore can only have $`\alpha _s(Q)\mathrm{ln}(\mu /Q)`$ dependence on $`Q`$. This reproduces the Bjorken scaling of the structure functions.

Finally, consider the tree level matching onto the Wilson coefficients shown in Fig. 3. From these graphs only the quark coefficient functions $`C_j^{(i)}`$ can be non-zero and we find

$`\mathrm{Im}H_1^{(i)}(z)=Q_i^2\pi \delta (z1),\mathrm{Im}H_2^{(i)}=0,`$ (93)
where $`Q_i`$ is the charge of parton $`i`$. The vanishing of $`\mathrm{Im}H_2^{(i)}`$ at tree level reproduces the Callan-Gross relation $`W_1/W_2=Q^2/(4x^2)`$.

### IV.2 Drell-Yan, $`p\overline{p}\mathrm{}^+\mathrm{}^{}X`$

Next we will extend the DIS analysis to the Drell-Yan (DY) process: $`p\overline{p}\mathrm{}^+\mathrm{}^{}X`$. Specifically we consider the $`Q^2`$ distribution, where $`Q^2`$ is the invariant mass of the lepton pair. Drell-Yan is more complicated than DIS because one has two hadrons in the initial state. In the center-of-mass frame the incoming proton and anti-proton move in opposite lightlike directions, and to prove factorization we use the fact that collinear modes in different lightlike directions can only couple to each other in external operators in SCET. We take the incoming proton to move in the $`n^\mu `$ direction and the incoming antiproton to move in the $`\overline{n}^\mu `$ direction. The hard scales in DY are $`Q^2`$ and the invariant mass of the colliding $`p\overline{p}`$ pair $`s=(p+\overline{p})^2`$. The lepton pair has an invariant mass $`Q^2`$, and the invariant mass of the final hadronic state is

$`p_X^2=Q^2\left(1+{\displaystyle \frac{1}{\tau }}{\displaystyle \frac{1}{x_1}}{\displaystyle \frac{1}{x_2}}\right),`$ (94)
where

$`\tau ={\displaystyle \frac{Q^2}{s}},x_1={\displaystyle \frac{Q^2}{2pq}},x_2={\displaystyle \frac{Q^2}{2\overline{p}q}}.`$ (95)
We are interested in the kinematic region where $`p_X^2Q^2`$, which implies that both $`x_1`$ and $`x_2`$ are far away from one. As $`\tau `$ approaches one the invariant mass becomes too small for the treatment given here to apply. However, the effective theory can be used to deal with this region as well. It is also possible to study the $`q_{}`$ distribution, but this again requires a generalization of the discussion given below.

The spin averaged cross section for Drell-Yan is

$`d\sigma ={\displaystyle \frac{32\pi ^2\alpha ^2}{Q^4s}}L_{\mu \nu }W^{\mu \nu }{\displaystyle \frac{d^3k_1}{(2\pi )^3(2k_1^0)}}{\displaystyle \frac{d^3k_2}{(2\pi )^3(2k_2^0)}},`$ (96)
where

$`W^{\mu \nu }`$ $`=`$ $`{\displaystyle \frac{1}{4}}{\displaystyle \underset{\mathrm{spins}}{}}{\displaystyle \underset{X}{}}(2\pi )^4\delta ^{(4)}(p+\overline{p}qp_X)p\overline{p}|J^\mu (0)|XX|J^\nu (0)|p\overline{p}`$ (97)
$`=`$ $`{\displaystyle \frac{1}{4}}{\displaystyle \underset{\mathrm{spins}}{}}{\displaystyle d^4xe^{iqx}p\overline{p}|J^\mu (x)J^\nu (0)|p\overline{p}}.`$
The sum over spins refers to the initial hadron spins (the sum over final hadron spins is included in the sum over $`X`$). Integrating Eq. (96) over the emission angles of the final leptons one obtains

$`{\displaystyle \frac{d\sigma }{dQ^2}}={\displaystyle \frac{2\alpha ^2}{3Q^2s}}{\displaystyle \frac{1}{4}}{\displaystyle \underset{\mathrm{spins}}{}}p\overline{p}|\widehat{W}|p\overline{p},`$ (98)
where we have neglected the lepton masses and defined the operator

$`\widehat{W}(\tau ,Q^2)={\displaystyle \frac{d^4q}{(2\pi )^3}\theta (q^0)\delta (q^2Q^2)d^4xe^{iqx}J^\mu (x)J_\mu (0)}.`$ (99)
As we discussed above in the region of phase space under consideration $`p_X^2Q^2`$, so these hard fluctuations can be integrated out. Operationally this means we match $`\widehat{W}`$ onto local operators in the effective theory. We would like to show that the minimal set of order $`\lambda ^4`$ operators that contribute to Drell-Yan are

$`\widehat{W}`$ $`{\displaystyle \frac{1}{Q^2}}{\displaystyle 𝑑\omega _iC_{qq}(\omega _i,Q)\left[\overline{\chi }_{n,\omega _1}^{(i)}\overline{n}/\chi _{n,\omega _2}^{(i)}\right]\left[\overline{\chi }_{\overline{n},\omega _3}^{(i)}n/\chi _{\overline{n},\omega _4}^{(i)}\right]}`$ (100)
$``$ $`{\displaystyle \frac{1}{Q^3}}{\displaystyle 𝑑\omega _iC_{qg}(\omega _i,Q)\left[\overline{\chi }_{n,\omega _1}^{(i)}\overline{n}/\chi _{n,\omega _2}^{(i)}\right]\mathrm{Tr}\left[_{\overline{n},\omega _3}^\beta _\beta ^{\overline{n},\omega _4}\right]}`$
$``$ $`{\displaystyle \frac{1}{Q^3}}{\displaystyle 𝑑\omega _iC_{gq}(\omega _i,Q)\mathrm{Tr}\left[_{n,\omega _1}^\nu _\nu ^{n,\omega _2}\right]\left[\overline{\chi }_{\overline{n},\omega _3}^{(i)}n/\chi _{\overline{n},\omega _4}^{(i)}\right]}`$
$`+`$ $`{\displaystyle \frac{1}{Q^4}}{\displaystyle 𝑑\omega _iC_{gg}(\omega _i,Q)\mathrm{Tr}\left[B_{n,\omega _1}^\nu B_\nu ^{n,\omega _2}\right]\mathrm{Tr}\left[_{\overline{n},\omega _3}^\beta _\beta ^{\overline{n},\omega _4}\right]},`$
where the powers of $`Q`$ are included to make the coefficients dimensionless. The operators displayed in Eq. (100) are just products of the operators that occured in DIS, so for these terms the decoupling of soft and usoft gluons occurs in a straightforward manner. To show that the operators in Eq. (100) are the most general set needed we must show that all other operators that are order $`\lambda ^4`$ either reduce to these or vanish between the matrix elements in Eq. (98). For instance, $`\lambda ^4`$ operators also exist where a $`B_{n,\omega }^\nu `$ field is contracted with a $`B_{\overline{n},\omega ^{}}^\mu `$ field, or the color structures of the operators in Eq. (100) could be arranged in a different way.

We now give a general argument for why we can always rewrite an arbitrary operator in the form of Eq. (100) or show that it does not contribute to DY. All operators relevant for DY contain four order $`\lambda `$ collinear fields chosen from $`\xi _{n,p}`$, $`\xi _{\overline{n},p}`$, $`B_{n,p}^\mu `$, or $`B_{\overline{n},p}^\mu `$. Furthemore, two must move in direction $`n`$ and two in the direction $`\overline{n}`$ (other possibilities end up vanishing by baryon number conservation or because they involve a set of fields between physical states that can not possibly form a color singlet operator). For operators with 4 quark fields, Fierz transformations can always be made to arrange the fields such that those in the same direction sit in the same bilinear. Using as an example the operator with two collinear quarks in the $`n`$ directions and two gluons in the $`\overline{n}`$ direction and leaving out the soft Wilson lines for the moment, the most general matrix element is

$`p_n\overline{p}_{\overline{n}}|\overline{\chi }_{n,\omega _1}^{a,\alpha }\chi _{n,\omega _2}^{b,\beta }_{\overline{n},\omega _3}^{A,\mu }_{\overline{n},\omega _4}^{B,\nu }|p_n\overline{p}_{\overline{n}}\mathrm{\Delta }_{\mu \nu ;\alpha \beta }^{ab;AB},`$ (101)
where $`a,b`$ are quark colors, $`A,B`$ are gluon colors, and $`\alpha ,\beta `$ are spinor indices for the quarks. $`\mathrm{\Delta }_{\mu \nu ;\alpha \beta }^{ab;AB}`$ is some tensor that connects the indices in an arbitrary way. In the contraction of $`a,b`$ and $`A,B`$ there are two possible ways to make an overall color singlet, one where both the quarks and gluons are in a color singlet, and another where both the quarks and gluons are in a color octet. We will discuss both of these possibilities in turn.

In the color singlet case, including the soft and ultrasoft Wilson lines is trivial, since using Eqs. (II.2), (II.2), and (A) we see that they cancel due to unitarity/orthogonality of the various Wilson lines in the fundamental/adjoint representations. Thus, there are no soft, usoft, or collinear interactions that connect the $`n`$ and the $`\overline{n}`$ fields. As in previous sections this leads to a factorization of the matrix element in Eq. (101), namely

$`p_n|\overline{\chi }_{n,\omega _1}^{(0)a,\alpha }\chi _{n,\omega _2}^{(0)b,\beta }|p_n\overline{p}_{\overline{n}}|_{\overline{n},\omega _3}^{(0)A,\mu }_{\overline{n},\omega _4}^{(0)B,\nu }|\overline{p}_{\overline{n}}\mathrm{\Delta }_{\mu \nu ;\alpha \beta }^{ab;AB}.`$ (102)
Since the proton spins are summed over, we can write (with the help of Eq. (II.4))

$`p_n|\overline{\chi }_{n,\omega _1}^{(0)a,\alpha }\chi _{n,\omega _2}^{(0)b,\beta }|p_n\delta ^{ab}(n/)_{\alpha \beta }p_n|\overline{\chi }_{n,\omega _1}^{(0)c,\gamma }\overline{n}/\chi _{n,\omega _2}^{(0)c,\gamma }|p_n,`$ (103)
so that spin and color are summed over in the matrix element. Similarly the antiproton matrix element can be simplified to

$`\overline{p}_{\overline{n}}|_{\overline{n},\omega _3}^{(0)A,\mu }_{\overline{n},\omega _4}^{(0)B,\nu }|\overline{p}_{\overline{n}}\delta ^{AB}g_{}^{\mu \nu }\overline{p}_{\overline{n}}|\mathrm{Tr}\left[_{\overline{n},\omega _3}^{(0)\alpha }_\alpha ^{(0)\overline{n},\omega _4}\right]|\overline{p}_{\overline{n}}.`$ (104)
Here we used the fact that the matrix element is symmetric in $`\mu `$ and $`\nu `$, and that only the perpendicular index $`\mu `$ of the field $`B^\mu `$ is order $`\lambda `$. Using Eqs. (103) and (104) the original matrix element in Eq. (101) can be written as

$`p_n\overline{p}_{\overline{n}}|\overline{\chi }_{n,\omega _1}^{a,\alpha }\chi _{n,\omega _2}^{b,\beta }_{\overline{n},\omega _3}^{A,\mu }_{\overline{n},\omega _4}^{B,\nu }|p_n\overline{p}_{\overline{n}}\mathrm{\Delta }_{\mu \nu ;\alpha \beta }^{ab;AB}`$ (105)
$`\mathrm{Tr}[\mathrm{\Delta }_\mu ^\mu n/]\overline{p}_{\overline{n}}|\mathrm{Tr}\left[_{\overline{n},\omega _3}^{(0)\nu }_\nu ^{(0)\overline{n},\omega _4}\right]|\overline{p}_{\overline{n}}p_n|\overline{\chi }_{n,\omega _1}^{(0)}\overline{n}/\chi _{n,\omega _2}^{(0)}|p_n,`$
where the trace of $`\mathrm{\Delta }`$ is over spin and color, and just gives an overall constant. The final result in Eq. (105) is identical to the matrix element of the second operator in Eq. (100).

If each of the $`n`$ and $`\overline{n}`$ field bilinears involve color octet structures, then the soft and usoft Wilson lines don’t cancel, since they don’t commute with the SU(3) generators. However, one can use the color identity

$`Y_n^{}S_n^{}𝒯^xS_{\overline{n}}Y_{\overline{n}}Y_{\overline{n}}^{}S_{\overline{n}}^{}𝒯^xS_nY_n=𝒯^xS_nY_nY_{\overline{n}}^{}S_{\overline{n}}^{}𝒯^xS_nY_nY_{\overline{n}}^{}S_{\overline{n}}^{}`$ (106)
where each $`𝒯^x`$, $`S`$, and $`Y`$ factor is in the appropriate representation of the color group. Eq. (106) moves all the soft and ultrasoft interactions between either the $`n`$ or the $`\overline{n}`$ collinear fields. Thus, again the fields in one bilinear can not be contracted with fields in the other bilinear and the matrix element factors. However this time the factored matrix element vanishes. For the example discussed above,

$`p_n|\overline{\chi }_{n,\omega }^{(0)\alpha }T^C\chi _{n,\omega ^{}}^{(0)\beta }|p_n=0,`$ (107)
since an color octet operator vanishes between color singlet states. The same holds true for the matrix element of an octet gluon operator.

An identical proof of decoupling goes through for the case of 4 quarks, where we again either have two color singlet or two color octet $`n`$ and $`\overline{n}`$ bilinears. With 4 gluon fields we can either have the two $`n`$ and two $`\overline{n}`$ fields coupled as singlets, or coupled in the same higher representation (an 8, {10, $`\overline{\mathrm{𝟏𝟎}}`$}, or 27). In the latter case the matrix element between color singlet states still vanishes so the proof for the 4 gluon operators also goes through in an identical way.

Thus we have shown that the matrix element of an operator with an arbitrary contraction of indices either vanishes or can be written in terms of a product of a matrix element which is related to a proton pdf and a matrix element which is related to an antiproton pdf as in the example in Eq. (105). This is the result we want. To see how the final formulae are derived note that we can write the matrix element of Eq. (100) in the form of a convolution

$`{\displaystyle \frac{1}{4}}{\displaystyle \underset{\mathrm{spins}}{}}p_n\overline{p}_{\overline{n}}|\widehat{W}|p_n\overline{p}_{\overline{n}}`$ $`=`$ $`{\displaystyle \underset{a,b}{}}{\displaystyle 𝑑\omega _iC_{a,b}(\omega _i)p_n|O_n^a(\omega _+,\omega _{})|p_n\overline{p}_{\overline{n}}|O_{\overline{n}}^b(\omega _+^{},\omega _{}^{})|\overline{p}_{\overline{n}}},`$
where $`\omega _\pm =\omega _1\pm \omega _2`$ and $`\omega _\pm ^{}=\omega _3\pm \omega _4`$. The operators here are the same as in DIS, with $`a=(i)`$ for the quark operator, and $`a=g`$ for the gluon operator

$`O_n^{(i)}(\omega _+,\omega _{})`$ $`=`$ $`{\displaystyle \frac{1}{Q}}\left[\overline{\chi }_{n,\omega _1}^{(i)}{\displaystyle \frac{\overline{n}/}{2}}\chi _{n,\omega _2}^{(i)}\right],`$
$`O_n^g(\omega _+,\omega _{})`$ $`=`$ $`{\displaystyle \frac{1}{Q^2}}\mathrm{tr}\left[B_{n,\omega _1}^\mu B_\mu ^{n,\omega _2}\right].`$ (109)
Apart from the dependence on the labels, the Wilson coefficients in Eq. (IV.2) can also depend on the renormalization point $`\mu `$ and the kinematic variable $`Q`$. Using Eq. (II.3) we see that the matrix elements in Eq. (IV.2) set $`\omega _{}=\omega _{}^{}=0`$ and $`\omega _+=2\sqrt{s}z_1`$, $`\omega _+^{}=2\sqrt{s}z_2`$ where $`z_1`$ and $`z_2`$ are the convolution variables. Since all kinematic variables aside from $`Q^2`$ are integrated over in Eq. (99) the only other variable that can appear in the Wilson coefficient is the center of mass energy which produces the $`\mathrm{}^+\mathrm{}^{}`$ pair, namely $`4\widehat{s}=\omega _+\omega _+^{}`$. Thus, the Wilson coefficients only depend on $`\omega _+\omega _+^{}=4sz_1z_2`$. Defining new coefficients

$`H^{ab}(z_1z_2)`$ $`=`$ $`C^{ab}(\omega _+\omega _+^{}=4sz_1z_2,Q,\mu )`$ (110)
we can replace the matrix elements in Eq. (IV.2) with parton distribution functions using Eq. (II.3) to obtain:

$`{\displaystyle \frac{1}{4}}{\displaystyle \underset{\mathrm{spins}}{}}p_n\overline{p}_{\overline{n}}|\widehat{W}|p_n\overline{p}_{\overline{n}}`$ $`=`$ $`{\displaystyle \frac{1}{\tau }}{\displaystyle _0^1}dz_1dz_2\{H^{(i)(j)}(z_1z_2)[f_{i/p}(z_1)\overline{f}_{j/\overline{p}}(z_2)+\overline{f}_{i/p}(z_1)f_{j/\overline{p}}(z_2)]`$ (111)
$`+H^{(i)(j)}(z_1z_2)\left[f_{i/p}(z_1)f_{j/\overline{p}}(z_2)+\overline{f}_{i/p}(z_1)\overline{f}_{j/\overline{p}}(z_2)\right]`$
$`+{\displaystyle \frac{z_2}{2\sqrt{\tau }}}H^{(i)g}(z_1z_2)f_{i/p}(z_1)f_{g/\overline{p}}(z_2){\displaystyle \frac{z_2}{2\sqrt{\tau }}}H^{(i)g}(z_1z_2)\overline{f}_{i/p}(z_1)f_{g/\overline{p}}(z_2)`$
$`+{\displaystyle \frac{z_1}{2\sqrt{\tau }}}H^{g(j)}(z_1z_2)f_{g/p}(z_1)f_{j/\overline{p}}(z_2){\displaystyle \frac{z_1}{2\sqrt{\tau }}}H^{g(j)}(z_1z_2)f_{g/p}(z_1)\overline{f}_{j/\overline{p}}(z_2)`$
$`+{\displaystyle \frac{z_1z_2}{4\tau }}H^{gg}(z_1z_2)f_{g/p}(z_1)f_{g/\overline{p}}(z_2)\}.`$
This is the final convolution formulae for Drell-Yan and is valid to all orders in $`\alpha _s`$ and leading order in the power expansion. At tree level the matching calculation shown in Fig. 4 yields zero for all the Wilson coefficients except

$`H^{(i)(i)}(z_1z_2)={\displaystyle \frac{2\pi \tau }{3}}Q_i^2\delta (\tau z_1z_2),`$ (112)
where $`Q_i`$ is the charge of parton $`i`$. The coefficients $`H^{(i)g}(\pm z_1z_2)`$ and $`H^{g(j)}(\pm z_1z_2)`$ start at order $`\alpha _s(Q)`$, while $`H^{(i)(i)}(z_1z_2)`$ and $`H^{gg}(z_1z_2)`$ start at order $`\alpha _s^2(Q)`$.

### IV.3 Deeply Virtual Compton Scattering, $`\gamma ^{}p\gamma ^{()}p^{}`$

Next we examine deeply virtual Compton scattering (DVCS). To be more precise we examine the exclusive reaction $`\gamma ^{}p\gamma ^{()}p`$, where the incoming photon is highly virtual, the final photon is either off-shell or real, and the incoming and outgoing protons have different momenta. The reason we have included this process in the inclusive section is that DVCS has the remarkable property that the nonperturbative physics is described by a so called non-forward parton distribution function (NFPDF). The NFPDF is a more general distribution function that reduces to the standard parton distribution functions (familiar from DIS) for some values of the momentum fraction, and behaves like a lightcone wavefunction (familiar from the pion examples) for other values. Deeply virtual Compton scattering was first studied in perturbative QCD in Refs Ji:1996ek ; Ji:1996nm ; Radyushkin:1996nd , and proofs of factorization to all orders in perturbation theory were later presented in Refs Collins:1998be ; Ji:1998xh . In addition properties of NFPDFs were studied in Ref. Radyushkin:1997ki . Here we present a proof of factorization for DVCS based on SCET.

As with the previous proofs it is important to understand the kinematics of the process. We take the incoming proton and photon momenta to be $`p`$ and $`q`$ respectively, with $`x=Q^2/(2pq)`$ and $`q^2=Q^2\mathrm{\Lambda }_{\mathrm{QCD}}^2`$. The outgoing proton and photon momentum are $`p^{}`$ and $`q^{}`$ respectively, with $`0q^2Q^2`$. It is convenient to define a parameter $`\zeta 1\overline{n}p^{}/\overline{n}p`$, which measures the change to the proton’s large momentum. Working in the Breit frame and neglecting contributions that are $`\mathrm{\Lambda }_{\mathrm{QCD}}^2/Q`$ we have:

(119)
Here the label momenta are order $`Q`$ or $`Q\lambda `$, while the residual momenta are order $`Q\lambda ^2`$ and depend on $`p^2=m_p^2`$ and $`t=(p^{}p)^2=(p_{}^2\zeta ^2m_p^2)/(1\zeta )`$, which are both $`\mathrm{\Lambda }_{\mathrm{QCD}}^2`$. The invariant mass of the intermediate hadronic state is $`(p+q)^2Q^2(1x)/x`$ just like DIS, so for $`1x\mathrm{\Lambda }_{\mathrm{QCD}}/Q`$ the intermediate state can be integrated out.

We will proceed in a manner analogous to the analysis for DIS. The amplitude (up to an overall momentum conserving $`\delta `$-function) is given by a time ordered product of currents:

$`T_{\mu \nu }(p,q,q^{})`$ $`=`$ $`p^{},\sigma ^{}|\widehat{T}_{\mu \nu }(q,q^{})|p,\sigma `$
$`\widehat{T}_{\mu \nu }(q,q^{})`$ $`=`$ $`i{\displaystyle d^4ze^{i(q+q^{})z/2}T\left[J_\mu (z/2)J_\nu (z/2)\right]}.`$ (120)
This time ordered product is contracted with a lepton tensor to obtain the amplitude. Now current conservation requires $`q^\mu T_{\mu \nu }=q^\nu T_{\mu \nu }=0`$, however the DVCS $`T_{\mu \nu }`$ is not symmetric under $`\mu \nu `$. For electromagnetic currents $`J_\mu `$ we have

$`T_{\mu \nu }`$ $`=`$ $`\left(g_{\mu \nu }{\displaystyle \frac{q_\mu ^{}q_\nu }{qq^{}}}\right)T_1+\left(p_\mu {\displaystyle \frac{q_\mu ^{}pq}{qq^{}}}\right)\left(p_\nu {\displaystyle \frac{q_\nu pq^{}}{qq^{}}}\right)T_2`$ (121)
$`+\mathrm{}_\mu \mathrm{}_\nu ^{}T_3+\mathrm{}_\mu \left(p_\nu {\displaystyle \frac{q_\nu pq^{}}{qq^{}}}\right)T_4+\left(p_\mu {\displaystyle \frac{q_\mu ^{}pq}{qq^{}}}\right)\mathrm{}_\nu T_5+\mathrm{},`$
where the functions $`T_i=T_i(x,\zeta ,Q^2,t)`$, and the vectors $`\mathrm{}_\mu q_\mu ^{}q_\mu +p_\mu (q^2q^2+t)/(2pq)`$ and $`\mathrm{}_\mu ^{}q_\mu q_\mu ^{}+p_\mu (q^2q^2+t)/(2pq^{})`$ are defined so that $`q\mathrm{}=q^{}\mathrm{}^{}=0`$. In Eq. (121) and below the ellipses denote spin dependent terms. For simplicity we will show how factorization is achieved for the spin independent contributions shown in Eq. (121) with the understanding that it is no more difficult to also include the other terms.

It is convenient to define a parameter $`0\alpha 1`$, by $`q^2\alpha Q^2`$. The DIS hadronic time-ordered product is obtained in the limit $`p^{}p`$, where $`\alpha 1`$ and $`\zeta 0`$. From Eq. (119) we see that

$`\zeta =x(1\alpha )+𝒪({\displaystyle \frac{t}{Q^2}},{\displaystyle \frac{m_p^2}{Q^2}}),`$ (122)
so these parameters are not independent. Since the intermediate hadronic state has invariant mass $`𝒪(Q^2)`$ we can match $`\widehat{T}_{\mu \nu }`$ onto operators in SCET. Requiring $`q^\mu \widehat{T}_{\mu \nu }=0`$ and $`q^\nu \widehat{T}_{\mu \nu }=0`$ for the order $`Q`$ label momenta leads to

$`\widehat{T}_{\mu \nu }{\displaystyle \frac{g_{}^{\mu \nu }}{Q}}\left(O_1^{(i)}+{\displaystyle \frac{O_1^g}{Q}}\right)+{\displaystyle \frac{1}{Q}}\left(n^\mu +\overline{n}^\mu \right)\left(\alpha n^\nu +\overline{n}^\nu \right)\left(O_2^{(i)}+{\displaystyle \frac{O_2^g}{Q}}\right)+\mathrm{},`$ (123)
where the ellipses are spin dependent terms and the displayed operators are

$`O_j^{(i)}`$ $`=`$ $`\overline{\xi }_{n,p^{}}^{(i)}W{\displaystyle \frac{\overline{n}/}{2}}C_j^{(i)}(\overline{𝒫}_+,\overline{𝒫}_{})W^{}\xi _{n,p}^{(i)},`$
$`O_j^g`$ $`=`$ $`\overline{n}_\mu \overline{n}_\nu \mathrm{tr}\left[W^{}(G_n)^{\mu \lambda }WC_j^g(\overline{𝒫}_+,\overline{𝒫}_{})W^{}(G_n)_\lambda ^\nu W\right].`$ (124)
We have suppressed the dependence of the Wilson coefficients $`C(\overline{𝒫}_+,\overline{𝒫}_{})`$ on $`Q`$, $`\alpha `$, and $`\mu `$. The form of the operators in Eq. (IV.3) looks the same as the DIS operators given in Eq. (IV.1), however the operators here are more general because the Wilson coefficients depend on $`\alpha `$. In the limit $`\alpha 1`$ the DVCS operators reduce to the DIS operators. However, since the field structure of the DVCS operators is identical to DIS several results follow immediately. For instance, the steps which factorize soft and usoft gluons and leave fields with superscript $`(0)`$ are the same and are not repeated here,

$`O_j^{(i)}`$ $`=`$ $`\overline{\chi }_{n,\omega }^{(0)(i)}{\displaystyle \frac{\overline{n}/}{2}}C_j^{(i)}(\overline{𝒫}_+,\overline{𝒫}_{})\chi _{n,\omega ^{}}^{(0)(i)},`$
$`O_j^g`$ $`=`$ $`\mathrm{tr}\left[B_{n,\omega }^{(0)\mu }C_j^g(\overline{𝒫}_+,\overline{𝒫}_{})(B_{n,\omega ^{}}^{(0)})_\mu \right].`$ (125)
The restrictions on the DVCS Wilson coefficients from charge conjugation are the same as in Eq. (90), $`C_j(\overline{𝒫}_+,\overline{𝒫}_{})=C_j(\overline{𝒫}_+,\overline{𝒫}_{})`$, however because $`pp^{}`$ this is not simply a relation between quark and anti-quark Wilson coefficients. The way in which DVCS is unique is that the matrix elements involve nucleon states with different momenta. This is what leads to results in terms of non-forward parton distribution functions.

The definition of the NFPDFs are given in Eq. (43), and can be used along with the relations above to obtain expressions for the $`T_i`$ in terms of the NFPDFs. Before we give this result we note that the Wilson coefficients depend on the operators $`\overline{𝒫}_\pm `$ which become the variables $`\omega _\pm `$ after introducing trivial convolutions and the $`\chi _{n,\omega _i}`$ fields in Eq. (35). The delta functions in Eq. (43) then set $`\omega _{}=Q\zeta /x`$ and $`\omega _+=\pm Q(2\xi \zeta )/x`$, where $`\xi `$ is the convolution variable. Note that $`\zeta /x=1\alpha `$, and just like DIS it is the combination $`\xi /x`$ which appears. Since charge conjugation relates the Wilson coefficients for $`\omega _+>0`$ and $`\omega _+<0`$ it is convenient to define

$`_j(\xi /x)C_j(Q(2\xi /x1+\alpha ),Q(\alpha 1),Q,\alpha ,\mu ),`$ (126)
where in the last three arguements we have made the dependence on $`Q`$, $`\alpha `$, and $`\mu `$ explicit. Combining Eqs. (43), (IV.3), and (126) then gives

$`T_1`$ $`=`$ $`{\displaystyle \frac{e(\sigma ^{},\sigma )}{2Q}}{\displaystyle _0^1}𝑑\xi \left\{_1^{(i)}\left({\displaystyle \frac{\xi }{x}}\right)\left[_\zeta ^i(\xi ;t)+\overline{}_\zeta ^i(\xi ;t)\right]+{\displaystyle \frac{1}{2x}}_1^g\left({\displaystyle \frac{\xi }{x}}\right)_\zeta ^g(\xi ;t)\right\}+\mathrm{},`$
$`T_2`$ $`=`$ $`{\displaystyle \frac{x^2(1+\alpha )}{Q^3}}e(\sigma ^{},\sigma ){\displaystyle _0^1}d\xi \{[2(1+\alpha )_2^{(i)}\left({\displaystyle \frac{\xi }{x}}\right)_1^{(i)}\left({\displaystyle \frac{\xi }{x}}\right)][_\zeta ^i(\xi ;t)+\overline{}_\zeta ^i(\xi ;t)]`$
$`+{\displaystyle \frac{1}{2x}}[2(1+\alpha )_2^g\left({\displaystyle \frac{\xi }{x}}\right)_1^g\left({\displaystyle \frac{\xi }{x}}\right)]_\zeta ^g(\xi ;t)\}+\mathrm{},`$
$`T_3`$ $`=`$ $`0,T_4=0,T_5=0,`$ (127)
which are the final convolution results valid to all orders in $`\alpha _s`$ and leading order in the power expansion. The structure functions $`T_{3,4,5}`$ vanish since the vectors $`\mathrm{}^\mu =\mathrm{}^\mu =0`$ at leading order in the power expansion. The terms with ellipses are for the spin flip terms and involve the NFPDF $`𝒦`$ defined in Eq. (43). The results for these terms have a similar form to those in Eq. (IV.3).

Finally we match at tree level. The tree level diagram in QCD is the same as in Fig. 3 except the outgoing photon and proton have momenta $`q^{}`$ and $`p^{}`$ respectively. Only the quark Wilson coefficients are nonzero at tree level. We find

$`C_1^{(i)}(\omega _+,\omega _{},Q,\alpha )`$ $`=`$ $`e^2Q_i^2\left({\displaystyle \frac{2Q}{2Q+\omega _+\omega _{}}}{\displaystyle \frac{2Q}{2Q\omega _+\omega _{}}}\right)`$ (128)
$`C_2^{(i)}(\omega _+,\omega _{},Q,\alpha )`$ $`=`$ $`0,`$
which gives

$`_1^{(i)}\left({\displaystyle \frac{\xi }{x}}\right)`$ $`=`$ $`e^2Q_i^2\left({\displaystyle \frac{1}{1\xi /x}}{\displaystyle \frac{1}{1+(\xi \zeta )/x}}\right)`$ (129)
$`_2^{(i)}\left({\displaystyle \frac{\xi }{x}}\right)`$ $`=`$ $`0.`$
Since $`_2^{(i)}=0`$ at tree level, DVCS also obeys a Callan-Gross relation.

## V Conclusion

What we hope we have demonstrated here is the power of effective field techniques in the context of factorization for hard scattering processes. The explicit separation of modes and the implementation of gauge invariance for these modes greatly simplifies seemingly complex problems. What is normally accomplished by diagramatic Ward identities and induction techniques now falls out as a consequence of the gauge symmetry of operators in a low energy soft-collinear effective theory.

As we have emphasized the factorization formulae derived in this paper are not new. The purpose here was simply to extend the formalism introduced in bfl ; bfps ; cbis ; bpssoft ; bps  to cases with back-to-back collinear particles, and apply these ideas to more general processes than previously considered. Furthermore, the factorization proofs presented are perhaps simpler than those previously given (certainly they are more concise). We believe that within the confines of the SCET more difficult, and unresolved problems can be addressed, such as power corrections in cases without an OPE, and proofs of factorization for more complex processes.

###### Acknowledgements.

This work was supported in part by the Department of Energy under the grant DOE-FG03-97ER40546.

## Appendix A Factorization of soft and collinear $`n`$ & $`\overline{n}`$ modes

This Appendix discusses the simultaneous factorization of the soft $`(\lambda ,\lambda ,\lambda )`$ modes, $`n`$-collinear $`(\lambda ^2,1,\lambda )`$ modes, and $`\overline{n}`$-collinear $`(1,\lambda ^2,\lambda )`$ modes. These three classes of modes can not interact with each other in a local manner and therefore do not couple through the SCET Lagrangian. However, they can couple in a gauge invariant way through external operators and currents. These interactions in currents are built up by integrating out offshell fluctuations with $`p^2(Q\lambda )^2`$. For the special case of factorization of soft from $`n`$-collinear modes this was shown in detail in the Appendix of Ref. bpssoft . There it was shown that integrating out certain modes with offshellness $`p^2Q^2\lambda `$ causes the Wilson lines $`W_n`$ and $`S_n`$ to appear in operators in a gauge invariant way. Here we will extend this approach to the factorization of modes for cases involving two classes of collinear particles. For simplicity we restrict ourselves to the case where the original operators involve only collinear quark or gluon fields. This type of factorization was used for the pion form factor example discussed in Sec. III.B and the Drell-Yan process presented in Sec. IV.B.

The basic idea is to first match onto a Lagrangian with couplings between onshell and offshell modes that give all order $`\lambda ^0`$ diagrams. The offshell modes (with $`p^2(Q\lambda )^2`$) can then be integrated out, so that all operators are expressed entirely in terms of the onshell degrees of freedom. In table 1 a summary is given of the three types of offshell momenta that are induced by adding soft, $`n`$-collinear, and $`\overline{n}`$-collinear momenta. For each type auxiliary quark and gluon fields are defined, and for convenience momentum labels are suppressed in this section. For example, the $`\psi _Q`$ quarks are created by the interaction of a $`n`$-collinear quark with an $`\overline{n}`$-collinear gluon, whereas the $`\psi _n^L`$ quarks are created when a collinear quark $`\xi _n`$ is knocked offshell by a soft gluon. For the field $`\psi _Q`$ we write $`\psi _Q=\psi _n^Q+\psi _{\overline{n}}^Q`$, where $`\psi _n^Q=\frac{1}{4}n/\overline{n}/\psi _Q`$ and $`\psi _{\overline{n}}^Q=\frac{1}{4}\overline{n}/n/\psi _Q`$. Then we have $`n/\psi _n^Q=\overline{n}/\psi _{\overline{n}}^Q=0`$ and $`n/\psi _n^L=\overline{n}/\psi _{\overline{n}}^L=0`$. Various Wilson lines are also required and are listed in the table. It is convenient to define a generic Wilson line $`L[a,A]`$ along direction $`a`$ with field $`A`$ by the solution of

$`\left(a𝒫+gaA\right)L[a,A]=0.`$ (130)
With this notation the on-shell Wilson lines are $`W_n=L[\overline{n},A_n]`$, $`W_{\overline{n}}=L[n,A_{\overline{n}}]`$, $`S_n=L[n,A_s]`$, and $`S_{\overline{n}}=L[\overline{n},A_s]`$. (Recall that the subscripts on $`W`$ and $`S`$ mean different things.) The Wilson lines involving offshell fields that we will require are

$`X_n`$ $`=`$ $`L[\overline{n},A_Q+A_n^X+A_n],X_{\overline{n}}=L[n,A_Q+A_{\overline{n}}^X+A_{\overline{n}}],`$ (131)
$`W_n^X`$ $`=`$ $`L[\overline{n},A_n^X+A_n],W_{\overline{n}}^X=L[n,A_{\overline{n}}^X+A_{\overline{n}}],`$
$`S_n^X`$ $`=`$ $`L[n,A_n^X+A_s],S_{\overline{n}}^X=L[\overline{n},A_{\overline{n}}^X+A_s].`$
Below we discuss the results which allow us to integrate out offshell fluctuations. The structure of the auxiliary Lagrangians and construction of their solutions are very similar to the case presented in Ref. bpssoft , to which we refer for a more detailed presentation.

From table 1 we see that adding $`n`$ and $`\overline{n}`$-collinear momenta gives $`p^2Q^2`$, whereas adding soft and collinear momenta gives $`p^2Q^2\lambda `$. Loops that are dominated by offshell momenta only modify Wilson coefficients and not the infrared physics. Therefore, to determine the structure of SCET fields in an operator it sufficient to integrate out the offshell fields at tree level. For convenience we can integrate out the fluctuations starting with those with the largest offshellness. Recall that we only wish to consider offshell propagators connected to external operators. A subtlety for quarks is that distinct auxiliary fields are needed for the incoming and outgoing offshell propagators. However, the solution for the outgoing field turns out to be the conjugate of the incoming field, so to avoid a proliferation of notation we simply denote the outgoing terms in the Lagrangian by $`+\text{h.c.}`$, and present a solution for the incoming fields. Finally, note that for the gluon field $`A_Q`$ the fields $`A_n`$, $`A_{\overline{n}}`$, $`A_n^X`$, $`A_{\overline{n}}^X`$, and $`A_s`$ appear as background fields while for the fields $`A_n^X`$ and $`A_{\overline{n}}^X`$ it is $`A_n`$, $`A_{\overline{n}}`$, and $`A_s`$ that appear as background fields.

The terms in the auxiliary Lagrangian involving the $`p^2Q^2`$ fields are

$`_{\mathrm{aux}}^Q`$ $`=`$ $`\overline{\psi }_n^Qgn(A_Q+A_{\overline{n}}^X+A_{\overline{n}}){\displaystyle \frac{\overline{n}/}{2}}(\psi _n^L+\xi _n)+\overline{\psi }_n^Q\left[n𝒫+gn(A_Q+A_{\overline{n}}^X+A_{\overline{n}})\right]{\displaystyle \frac{\overline{n}/}{2}}\psi _n^Q`$ (132)
$`+(n\overline{n})+\text{h.c.}`$
$`+{\displaystyle \frac{1}{2g^2}}\text{tr}\left\{[iD_Q^\mu +gA_Q^\mu ,iD_Q^\nu +gA_Q^\nu ]^2\right\}+{\displaystyle \frac{1}{\alpha _Q}}\text{tr}\{[iD_Q^\mu ,A_{Q\mu }]^2\}.`$
where $`iD_Q^\mu =\frac{1}{2}n^\mu [\overline{𝒫}+g\overline{n}(A_n^X+A_n)]+\frac{1}{2}\overline{n}^\mu [𝒫+gn(A_{\overline{n}}^X+A_{\overline{n}})]`$. The solution of the equations of motion for these modes are

$`\psi _n^Q=(X_{\overline{n}}1)(\psi _n^L+\xi _n),\psi _{\overline{n}}^Q=(X_n1)(\psi _{\overline{n}}^L+\xi _{\overline{n}}),`$
$`X_{\overline{n}}^{}X_n=W_n^XW_{\overline{n}}^X.`$ (133)
(In addition to the last equation a constraint on the components $`nA_Q`$ and $`\overline{n}A_Q`$ also comes from the gauge fixing term, but will not be needed.) The terms in the auxiliary Lagrangian involving the $`p^2Q^2\lambda `$ fields are bpssoft

$`_{\mathrm{aux}}^X`$ $`=`$ $`\overline{\psi }_n^Lgn(A_n^X+A_s){\displaystyle \frac{\overline{n}/}{2}}\xi _n+\overline{\psi }_n^L[n𝒫+gn(A_n^X+A_s)]{\displaystyle \frac{\overline{n}/}{2}}\psi _n^L+(n\overline{n})+\text{h.c.}`$
$`+`$ $`{\displaystyle \frac{1}{2g^2}}\text{tr}\left\{[iD_{nX}^\mu +gA_n^{X\mu },iD_{nX}^\nu +gA_n^{X\nu }]^2\right\}+{\displaystyle \frac{1}{\alpha _n}}\text{tr}\left\{[iD_{nX}^\mu ,A_{n\mu }^X]^2\right\}+(n\overline{n}),`$
where $`iD_{nX}^\mu =\frac{1}{2}n^\mu [\overline{𝒫}+g\overline{n}A_n]+\frac{1}{2}\overline{n}^\mu [n𝒫+gnA_s]`$. The solutions for these modes are

$`\psi _n^L=(S_n^X1)\xi _n,S_n^XW_n^X=W_nS_n^{},`$
$`\psi _{\overline{n}}^L=(S_{\overline{n}}^X1)\xi _{\overline{n}},S_{\overline{n}}^XW_{\overline{n}}^X=W_{\overline{n}}S_{\overline{n}}^{}.`$ (135)
Together Eqs. (A) and (A) can be used at leading order to eliminate the fields representing offshell fluctuations with $`p^2(Q\lambda )^2`$. For collinear quarks this leads to the rules in Eq. (II.2). Note that we did not need to use the gauge fixing term to resolve the ambiguity in the implicit solution for the $`\overline{n}A`$ and $`nA`$ auxillary fields.

As an illustration of these results we discuss the soft-collinear factorization for the production of a $`q_n\overline{q}_{\overline{n}}`$ pair with a large invariant mass $`Q^2`$. This process is mediated in the full theory by the electromagnetic current $`J=\overline{\psi }\mathrm{\Gamma }\psi `$ ($`\mathrm{\Gamma }`$ a color singlet). This current will match onto a current in SCET that is built entirely out of onshell fields. Using the results in this appendix this current can be systematically derived. To start the quark field in $`J`$ matches onto $`\xi _n`$ plus all possible fields which the auxiliary Lagrangian can create starting from $`\xi _n`$, so

$`J(\overline{\xi }_n+\overline{\psi }_n^L+\overline{\psi }_n^Q)\mathrm{\Gamma }(\xi _{\overline{n}}+\psi _{\overline{n}}^L+\psi _{\overline{n}}^Q).`$ (136)
Integrating out the $`p^2Q^2`$ fluctuations with Eq. (A) and inserting a hard Wilson coefficient $`C`$ which depends on label operators turns Eq. (136) into

$`(\overline{\xi }_n+\overline{\psi }_n^L)X_{\overline{n}}^{}C\mathrm{\Gamma }X_n(\xi _{\overline{n}}+\psi _{\overline{n}}^L)=(\overline{\xi }_n+\overline{\psi }_n^L)W_n^XC\mathrm{\Gamma }W_{\overline{n}}^X(\xi _{\overline{n}}+\psi _{\overline{n}}^L).`$ (137)
To construct the first operator we used the equations of motion for $`\psi _n^Q`$ and $`\psi _{\overline{n}}^Q`$, and in the second operator we used the equation of motion identity for the gluons in $`X_n`$ and $`X_{\overline{n}}`$. In a similar fashion we can now integrate out the $`p^2Q^2\lambda `$ fluctuations with Eq. (A) to give

$`\overline{\xi }_nS_n^XW_n^XC\mathrm{\Gamma }W_{\overline{n}}^XS_{\overline{n}}^X\xi _{\overline{n}}=\overline{\xi }_nW_nS_n^{}C\mathrm{\Gamma }S_{\overline{n}}W_{\overline{n}}^{}\xi _{\overline{n}}.`$ (138)
The operator on the right is the final result used in Eq. (24), and is soft, collinear, and usoft gauge invariant. It should be obvious from this example how the equations of motion in Eqs. (A) and (A) can be used to determine the factorized form of a general leading order operator.