the drell - yan process @xcite , the production of a lepton pair in hadron - hadron collisions , has played an important role in establishing the parton picture underlying the description of hard interactions in qcd . in current experiments , studies of the drell - yan cross section as a function of the invariant mass of the lepton pair
are used to search for new heavy particles such as a hypothetical @xmath5 boson , while the differential distributions provide detailed information about the parton distribution functions ( pdfs ) , including in particular the sea - quark distributions .
drell - yan production also serves as a prototype for other collider processes , such as higgs production or the production of new particles .
a lot of effort has been put in obtaining accurate theoretical predictions for the drell - yan process in perturbative qcd .
the calculation of the cross section and rapidity distribution at next - to - leading order ( nlo ) in @xmath6 was accomplished in the pioneering work @xcite .
the first next - to - next - to - leading order ( nnlo ) result for the cross section was presented in @xcite and confirmed much later in @xcite .
the rapidity distribution at nnlo was derived in @xcite , while the fully differential cross section was obtained recently in @xcite .
as the invariant mass @xmath7 of the lepton pair approaches the center - of - mass energy of the collision , there is less and less phase space available for the emission of qcd radiation .
after the cancellation of virtual and real soft divergences large sudakov logarithms remain , because the scale associated with the soft radiation is much smaller than @xmath7 .
these  threshold logarithms " threaten the convergence of the perturbative expansion and need to be resummed to all orders .
for the inclusive cross section @xmath8 this was accomplished in the seminal papers @xcite based on the solution of certain evolution equations in mellin moment space ( see also @xcite ) .
the generalization of this method to the drell - yan rapidity distribution was obtained in @xcite .
recent analyses of threshold resummation for the rapidity distribution in drell - yan or electroweak gauge boson production can be found in @xcite in next - to - leading logarithmic ( nll ) approximation , while in @xcite the resummation is extended to the next - to - next - to - next - to - leading logarithmic ( n@xmath9ll ) order .    because the pdfs are strongly suppressed in the endpoint region @xmath10 , the cross section @xmath8 is a steeply falling function as @xmath7 approaches the kinematical endpoint @xmath11 .
in fact , in a typical experiment it will not be possible to observe drell - yan pairs with masses exceeding about one half of the center - of - mass energy . for instance , at the lhc one does not expect to discover new heavy particles with masses in the 10tev range . in practice ,
one is therefore never in a region where the ratio @xmath12 approaches 1 .
since threshold resummation deals with logarithms of the form @xmath13 , it is then not obvious why such terms should be treated on different footing than other higher - order terms . in view of this , it is surprising that large resummation effects were recently claimed to be important for the drell - yan rapidity distribution as measured by the e866/nusea collaboration @xcite .
specifically , ref .
@xcite claims that for @xmath14gev and @xmath15gev ( corresponding to @xmath16 ) the resummation of threshold logarithms at nnll order would lower the cross section by about 30% compared with the fixed - order nlo result , whereas the fixed - order nnlo corrections increase it by a small amount @xcite .
however , significantly smaller resummation effects were found by other authors for lhc kinematics ( @xmath17gev and @xmath18tev ) @xcite and for @xmath19-boson production at rhic ( @xmath20gev and @xmath21gev ) @xcite . an argument why threshold resummation effects could be important even if @xmath22 has been given in refs .
the idea is that the sharp fall - off of the parton luminosity at large @xmath0 dynamically enhances the contribution of the _ partonic _ threshold region @xmath23 , i.e. , the region where the center - of - mass energy @xmath24 of the initial - state partons is just sufficiently large to produce the drell - yan pair .
it could then be important to resum logarithms of the form @xmath25 in the hard partonic cross section .
however , since @xmath26 is not related to a small ratio of external physical scales , it is not obvious how to give a formal justification of this argument .    to study this question
quantitatively and to assess the importance of resummation effects for the drell - yan rapidity distribution were the main motivations for the present work .
to do so , we use a recent approach to sudakov resummation based on effective field theory @xcite .
contrary to the standard treatment in mellin moment space , this framework completely separates the effects associated with different scales in the problem , thereby avoiding the landau - pole ambiguities inherent in the standard approach .
it then uses renormalization - group ( rg ) evolution to resum logarithms of scale ratios .
the resummation is performed directly in momentum space , which makes it simpler to compare to and match onto fixed - order calculations .
our framework is particularly well suited to study the resummed rapidity distribution , for which we derive an exact analytic expression as a one - dimensional integral over pdfs .
using the convergence properties of the perturbative expansion after scale separation as the primary criterion , we study in detail how and under which circumstances an effective physical scale @xmath27 emerges , which is associated with the soft emission in the process .
our approach resums logarithms of the ratio @xmath28 to all orders in perturbation theory .
our main findings can be summarized as follows :    * in the true endpoint region @xmath29 , the effective soft scale @xmath30 is an order of magnitude smaller than the naive guess @xmath31 . for pdfs behaving like @xmath32 near @xmath10 , we find @xmath33 with @xmath34 .
this result provides a formal justification to the argument of a dynamical enhancement of the partonic threshold region due to the fall - off of parton densities . *
the dynamical enhancement of the threshold contributions remains effective down to moderate values @xmath35 , while at very small @xmath36 values the parameter @xmath37 decreases to about 2 .
this reflects the fact that for small @xmath0 values the fall - off of the pdfs is much weaker than for large @xmath0 . * even far away from the true threshold the drell - yan cross section receives its dominant contributions from those terms in the hard partonic cross section that are leading in the limit @xmath38 . assuming this is true for other processes as well , the evaluation of virtual corrections plus soft emissions provides a simple and efficient way to obtain useful approximations for higher - order perturbative corrections . * with the appropriate choice of the effective soft scale @xmath30 , the convergence of the perturbative expansion is greatly improved by the resummation .
however , for small drell - yan masses the terms beyond @xmath39 in the resummed expression for the cross section are numerically unimportant .
we thus do not confirm the large impact of threshold resummation on the drell - yan rapidity distribution reported in @xcite .
for larger masses the effects can be significant .
for instance , the experiment e866/nusea has reported data up to @xmath40gev ( corresponding to @xmath41 ) @xcite .
we find that at @xmath42gev resummation effects enhance the fixed - order predictions for the cross section by about 25% at nlo , and 7% at nnlo .
* for the case of the integrated cross - section @xmath8 , we perform a detailed comparison with the traditional resummation approach in moment space .
similar to the case of deep - inelastic scattering ( dis ) in the region @xmath10 , we find that the two approaches are equivalent up to power corrections , which turn out to be numerically small .
an important conceptual difference is that in the effective - theory approach the running coupling is evaluated at physical short - distance scales depending only on the external variables @xmath43 and @xmath7 ( and perhaps the rapidity @xmath44 ) . in this way ,
the landau - pole ambiguities inherent in the standard approach are avoided .
we begin our analysis discussing the structure of the hard - scattering kernels relevant for the drell - yan rapidity distribution in fixed - order perturbative qcd . in section  [ sec : scet ] we use the framework of soft - collinear effective theory ( scet ) @xcite to derive the standard factorization formula for the partonic cross section in the limit @xmath38 in terms of hard and soft functions , which we define in terms of wilson coefficients of operators in the effective theory .
the solutions to the rg equations obeyed by these coefficients are derived in section  [ sec : resummation ] . with these results at hand ,
we present exact analytic expressions for the resummed drell - yan cross section and rapidity distribution .
a detailed numerical analysis of our results is presented in section  [ sec : numerics ] . after choosing the hard and soft matching scales in the effective theory by analyzing the perturbative expansions of the wilson coefficient functions , we investigate the stability of the results under scale variations and discuss the impact of the resummation . before concluding , we discuss the connection with the conventional moment - space approach
we consider the production of a lepton pair with invariant mass @xmath7 in hadron - hadron collisions at center - of - mass energy @xmath45 ( drell - yan process ) , focusing for simplicity on the reaction @xmath46 followed by @xmath47 .
our goal is to calculate the double differential cross section in the variables @xmath48 and @xmath49 , where @xmath44 denotes the rapidity of the lepton pair in the center - of - mass frame . up to power corrections
this cross section can be calculated in perturbative qcd and expressed in terms of convolutions of short - distance partonic cross sections with pdfs : @xmath50 here @xmath51 is the probability of finding a parton @xmath52 with longitudinal momentum fraction @xmath0 inside the hadron @xmath53 , and @xmath54 is the factorization scale .
the hard - scattering kernels @xmath55 have an expansion in powers of the strong coupling @xmath6 .
the sum extends over all possible partonic channels contributing at a given order in this expansion . at leading order ( @xmath56 )
only the channels @xmath57 contribute , while at nlo ( @xmath58 ) one must include @xmath59 in the sum
.    it will be useful for our purposes to introduce the ratios @xmath60 where @xmath61 is the center - of - mass energy squared of the partonic subprocess that creates the lepton pair .
this determines the maximum energy transferred to the leptons and the maximum invariant mass @xmath7 they can have . in @xcite the coefficient functions @xmath55 are expressed in terms of the variable @xmath62 and a second quantity @xmath63 these variables take values on the intervals @xmath64 and @xmath65 subject to the condition that the parton momentum fractions @xmath66 do not exceed 1 .
the allowed range for the rapidity is such that @xmath67 .
we then define new kernels via @xmath68[1-(1-y)(1-z ) ] } \,.\ ] ] at nlo the explicit results for these functions can be written in the form ( with @xmath69 and @xmath70 denoting the electric charges of the quarks in units of @xmath71 ) @xcite @xmath72 \nonumber\\     & & \mbox{}+ \frac{c_f\alpha_s}{\pi}\,\bigg\ {      \frac{\delta(y ) + \delta(1-y)}{2}\,\bigg [       ( 1+z^2 ) \left [ \frac{1}{1-z } \ln\frac{m^2(1-z)^2}{\muf^2 z }       \right]_+ + 1 - z \bigg ] \nonumber\\     & & \quad \mbox{}+ \frac12       \left [ 1 + \frac{(1-z)^2}{z}\,y(1-y ) \right ]      \left [ \frac{1+z^2}{1-z } \left ( \left [ \frac{1}{y } \right]_+       + \left [ \frac{1}{1-y } \right]_+
\right ) - 2(1-z ) \right ]      \bigg\ } \ , , \nonumber\\     \frac{c_{qg}}{e_q^2 }     & = & \frac{t_f\alpha_s}{2\pi}\,\bigg\ { \delta(y ) \left [      \left ( z^2 + ( 1-z)^2 \right ) \ln\frac{m^2(1-z)^2}{\muf^2 z }       + 2z(1-z ) \right ] \nonumber\\     & & \quad \mbox{}+ \left [ 1 + \frac{(1-z)^2}{z}\,y(1-y ) \right ]      \left [ \left ( z^2 + ( 1-z)^2 \right ) \left [ \frac{1}{y } \right]_+      + 2z(1-z ) + ( 1-z)^2 y \right ] \bigg\ } \ , .
\quad\end{aligned}\ ] ] the @xmath54-dependent terms can be derived from the fact that the cross section in ( [ sig ] ) is scale independent , while the pdfs obey the dglap evolution equations @xcite .
the remaining functions follow from the symmetry relations @xmath73 the hard - scattering kernels at nnlo have been calculated in @xcite and are available in the form of a computer program @xcite .
the explicit expressions for the coefficient functions given above contain terms that are singular in the  partonic threshold region " @xmath38 , in which the center - of - mass energy of the parton subprocess is just large enough to create a lepton pair with invariant mass @xmath7 .
indeed , the arguments of the logarithms in ( [ cijres ] ) suggest the relevance of two mass scales : a  hard " scale @xmath74 , and a  soft " scale @xmath75 .
physically , the hard scale is set by the invariant mass of the lepton pair , while the soft scale is of the order of the energy of the remnant jet @xmath76 produced in the collision . in the region of parton kinematics
where @xmath38 these scales are separated , @xmath77 , in which case the coefficient functions contain large logarithms irrespective of the choice of the factorization scale @xmath54 .
threshold resummation for the drell - yan cross section @xcite aims at resumming these logarithms to all orders in perturbation theory .
let us return to the structure of the relations ( [ cijres ] ) and identify the leading singular terms in the partonic threshold region .
they are contained in @xmath78 and up to nlo multiply @xmath79-functions in the variable @xmath80 . beyond nlo
some of the leading singular terms in the expressions obtained in @xcite multiply nontrivial functions of @xmath80 , but since the @xmath80-dependence of the parton variables @xmath81 and @xmath82 in ( [ xi ] ) is subleading in the @xmath38 limit one can always rearrange the expressions in such a way that the leading singular terms multiply @xmath79-functions in @xmath80 .
explicitly , we then obtain @xmath83 where @xmath84_+ \right\ } \nonumber\\     & & \mbox{}+ c_f \left ( \frac{\alpha_s}{\pi } \right)^2       \big [ c_f p_f(z ) + c_a p_a(z ) + t_f n_f p_f(z ) \big ] \,,\end{aligned}\ ] ] and we have defined @xmath85 and @xmath86 $ ] .
the terms in the first line can be readily read off from ( [ cijres ] ) .
the two - loop coefficients @xmath87 are given in appendix  [ app : a ] .
note that the factor @xmath62 in the argument of the logarithm @xmath88 could be set to 1 at leading order , but it is correctly reproduced by our resummation formula below and so we will keep it .
the goal of this paper is to derive a formalism that resums these leading terms to all orders in perturbation theory .
upon performing the integration over @xmath80 , the leading singular terms in ( [ llterms ] ) give rise to the following contribution to the cross section : @xmath89 . \nonumber\end{aligned}\ ] ] the lower limit of the @xmath62 integral is @xmath90 , as appropriate for the two terms . at tree level
, we recover the parton - model result @xmath91 .\ ] ] we also note that integrating over rapidity , we obtain for the leading singular terms in the single - differential cross section @xmath92 , \ ] ] where @xmath93 , and the integration is restricted to the region where @xmath94 .
this result can be rewritten in the more convenient form @xmath95 where @xmath96\ ] ] denotes the mellin convolution of the pdfs .    in the following section we will derive a factorization formula for the coefficient @xmath97 using methods of effective field theory .
the result is @xmath98 where @xmath99 and @xmath100 will be referred to as hard and soft functions , respectively , and will be defined in terms of wilson coefficients of operators in scet .
the calculation of the components @xmath99 and @xmath100 at any order in perturbation theory is much simpler than the calculation of the drell - yan cross section at the same order .
( [ factform ] ) thus provides an approximation to the cross section that requires a minimal amount of calculational work .
the all - order resummation of the partonic threshold logarithms is then achieved by solving rg equations .
it must be emphasized at this point that the variable @xmath62 is not set by external kinematics , but instead is integrated over the interval between @xmath12 and 1 .
it is therefore necessary to specify under which conditions the partonic threshold region requires special attention .
for instance , also in inclusive processes such as @xmath101 there are kinematical situations where scales much smaller than the center - of - mass energy are important , such as the emission of a soft gluon into the final state .
it is a well - known fact , however , that upon integration over the entire phase space the perturbation series is insensitive to scales lower than @xmath11 at leading power .
there are two limits in which the threshold contributions are parametrically enhanced .
first , if the invariant mass of the lepton pair is near the kinematic limit set by the total center - of - mass energy of the hadron - hadron collision , then @xmath102 and hence @xmath103 is always near 1 .
threshold resummation is necessary in this case and proceeds in close analogy to the resummation of threshold logarithms for the dis structure functions for @xmath10 @xcite . in practice ,
the region @xmath102 is irrelevant for phenomenology , since the strong suppression of the pdfs near the endpoint implies a very low parton luminosity in this case .
[ ] [ ] @xmath80 [ b]@xmath104   fall - off of the parton luminosity function @xmath104 for @xmath105gev .
the dashed lines show the asymptotic behavior for small and large @xmath80.,title="fig:",scaledwidth=50.0% ]    a second way in which the threshold contributions can be enhanced arises dynamically , if the weight function multiplying the hard - scattering kernel under the @xmath62-integral is steeply falling with @xmath26 @xcite . in this case threshold resummation
can be justified even if @xmath36 is much less than 1 . in practice ,
such a behavior has to result from the fall - off of the parton densities with increasing @xmath0 . consider for simplicity the total cross section @xmath8 , for which the relevant combination of pdfs is given by the function @xmath104 in ( [ ffdef ] ) .
figure  [ fig : ffdrop ] shows that this function is indeed very steeply falling with @xmath80 .
taking @xmath105gev , one finds that @xmath106 for @xmath107 and @xmath108 for @xmath109 , where @xmath110 and @xmath111 .
the figure shows that the first form reasonably well describes the behavior for @xmath112 , while the second form holds for @xmath113 . using these asymptotic forms for the parton luminosity function
, we find that for @xmath114 @xmath115 while for @xmath116 @xmath117 in the first case the cross section is given by a low - order moment of the hard - scattering kernel @xmath118 , in which case the partonic threshold region @xmath38 is not parametrically enhanced . in the second case ,
on the other hand , in the limit where we treat the exponent @xmath119 as a large parameter , the @xmath62 integral receives important contributions only from the region where @xmath120 . even for @xmath36 values
not near 1 there is thus a parametric enhancement of the partonic threshold region , which turns the threshold logarithms into logarithms of the exponent @xmath119 .
the intermediate range @xmath121 is a transition region , in which the dynamical enhancement of the threshold region ceases to be effective as @xmath36 is lowered .
cc [ ] [ ] @xmath44 [ b]@xmath122[pb / gev ]   comparison of the complete fixed - order results ( solid lines ) and the contributions arising from the leading singular terms ( dashed lines ) to the drell - yan rapidity distribution at different orders in perturbation theory . on the left we use the pdf sets mrst01lo , mrst04nlo , and mrst04nnlo , as appropriate to the order of the calculation
; in the right plot mrst04nnlo is used throughout.,title="fig:",scaledwidth=45.0% ] & [ ] @xmath44 [ b]@xmath122[pb / gev ] [ ]    @xmath123   q & \!=\ ! & 38.76\,\mbox{gev } \nonumber\\[-0.25 cm ]   m & \!=\ ! & 8\,\mbox{gev } \nonumber\\[-0.25 cm ]   \muf & \!=\ !
& m \nonumber   \end{aligned}\ ] ]     comparison of the complete fixed - order results ( solid lines ) and the contributions arising from the leading singular terms ( dashed lines ) to the drell - yan rapidity distribution at different orders in perturbation theory .
on the left we use the pdf sets mrst01lo , mrst04nlo , and mrst04nnlo , as appropriate to the order of the calculation ; in the right plot mrst04nnlo is used throughout.,scaledwidth=45.0% ]    while , as we have just discussed , the dominance of the partonic threshold region can not be justified parametrically for @xmath114 , it nevertheless appears that even in this case the drell - yan cross section and rapidity distribution are well approximated by keeping only the leading singular terms in the hard - scattering kernel ( [ cdecomp ] ) .
the reason is an inherent property of the hard - scattering kernel , which appears to receive the largest radiative corrections from the region of phase space corresponding to born kinematics . in other words ,
the effects of hard real emissions appear to be suppressed compared with virtual corrections and soft emissions .
to illustrate this point , we show in figure  [ fig : fixedvsresummed ] the full fixed - order predictions for the drell - yan rapidity distribution at @xmath15gev and @xmath14gev and compare them to the results obtained by keeping only the threshold terms in the hard - scattering kernel .
we use the mrst sets of pdfs compiled in @xcite . throughout our work
we use the three - loop running coupling @xmath124 normalized to @xmath125 , and we take @xmath126 for the number of light quark flavors . in the right plot
we use the same nnlo parton densities for all curves , so that the size of the different perturbative corrections to the hard - scattering kernel can be read off directly . since our main focus in this work is on the behavior of the perturbative expansion of the hard - scattering kernel
, we will from now on always use the set mrst04nnlo for the pdfs . even though @xmath16 is very small in this example
, we observe that at central rapidity about 80% of the nlo correction and 63% of the nnlo correction arise from the leading singular terms .
for the total cross section @xmath8 about 93% of the nlo correction and 97% of the nnlo correction are accounted for by the leading singular terms .
if we lower the factorization scale to @xmath127 , then the leading singular terms come even closer to reproducing the fixed - order results .
the factorization theorem ( [ factform ] ) for the leading singular terms in the hard - scattering kernel @xmath118 has been established a long time ago @xcite . in particular , it has been understood that the soft function @xmath100 can be represented as the vacuum expectation value of a certain wilson loop of soft gluon fields .
nevertheless we find it useful to rederive this formula using methods of effective field theory .
the advantage of this approach is that we will relate the hard and soft functions , @xmath99 and @xmath100 , to wilson coefficients of operators in the effective theory , which obey certain rg equations . solving these equations we accomplish threshold resummation directly in momentum space .    in the framework of scet the standard factorization formula
( [ sig ] ) for the drell - yan cross section has been discussed in @xcite .
however , the issue of glauber gluons @xcite has not yet been addressed in the effective theory .
glauber - gluon interactions have presented an important complication in the factorization analysis of the drell - yan process @xcite , and it would be worthwhile to rederive their cancellation using scet . as discussed below , glauber gluons do not affect the factorization of the hard - scattering kernel @xmath128 in the threshold region .    the threshold resummation for the total cross section @xmath8 has been performed in @xcite using an effective - theory variant of the traditional mellin moment - space approach . in the present paper we go beyond these works by presenting a derivation of the factorization formula ( [ factform ] ) in scet and performing the threshold resummation directly in @xmath62 space .
this will provide us with the tools to discuss the relevance of resummation away from the true threshold , i.e. , for @xmath129 much smaller than @xmath43 .
we begin with the standard formula for the drell - yan cross section , @xmath130 where @xmath131 is the electromagnetic current . to derive the factorization theorem
, we match the product of currents onto operators in scet .
the matching proceeds in two steps , and the corresponding wilson coefficients are the hard function @xmath99 and the soft function @xmath100 , respectively . the remaining effective - theory matrix element can be identified with the pdfs .
a similar two - step matching procedure has been used in many scet applications . in particular
, the factorization theorem for dis in the region @xmath10 was derived in scet by proceeding this way @xcite .    in the present paper
we are solely interested in the factorization of the hard - scattering coefficient @xmath128 . for this purpose
, we can simplify the effective - theory analysis by considering ( [ textbook ] ) with partonic instead of hadronic matrix elements .
these partonic matrix elements correspond directly to the hard - scattering coefficients , and the hadronic cross sections are then obtained after convoluting the results with the pdfs .
an effective - theory analysis of the hadronic matrix elements would be more complicated because of their sensitivity to nonperturbative physics governed by the scale @xmath132 , which does not enter the hard - scattering kernels .
note that the product of current operators in the matrix element in ( [ textbook ] ) is not time - ordered .
this is in contrast to the case of dis , where the decay rate can be obtained from the discontinuity of a time - ordered product of currents .
it is not possible to rewrite ( [ textbook ] ) in the same way , because the imaginary part of the time - ordered product also gets contributions from virtual corrections to @xmath133 , where @xmath76 is a purely hadronic final state .
a path integral framework to analyze operator products that are not time ordered is the keldysh formalism @xcite , which we review in appendix  [ sec : keldysh ] . in our discussion below , the operator ordering does not lead to any complications , but it implies that the soft function is not given by feynman diagrams , but appropriate cuts of such diagrams , as explained in @xcite .      before entering the technicalities of the discussion of factorization in scet , it is useful to understand the physics behind the factorization theorem ( [ factform ] ) by means of a simple , intuitive argument .
consider the special kinematics of drell - yan production in the ( partonic ) threshold region , as illustrated in figure  [ fig : cartoon ] . because the partonic center - of - mass energy is just above the invariant mass of the produced drell - yan pair , only soft emissions from the initial - state partons are allowed .
as is well known , these soft emissions are described by eikonal interactions and exponentiate into wilson lines .
furthermore , to leading power they leave the incoming partons on the mass shell , so that the production of the drell - yan pair is described by an on - shell quark form factor .
[ b]@xmath134 [ b]@xmath135 [ b]@xmath136   soft gluon emissions from the initial - state partons in drell - yan production.,title="fig:",scaledwidth=32.0% ]    we will now formalize this argument using the language of effective field theory . the effective theory is constructed by introducing fields for the momentum regions that contribute to the matrix elements in the given kinematics .
as stressed above , for our purposes it is sufficient to analyze partonic matrix elements . to distinguish the different momentum modes
, we introduce the light - cone decomposition @xmath137 where @xmath138 and @xmath139 are two light - like reference vectors ( @xmath140 ) in the directions of the colliding partons .
we denote the small expansion parameter by @xmath141 and quote the components @xmath142 of parton momenta .
the relevant momentum regions are @xmath143 where the large momentum scale in the process is set by the partonic center - of - mass energy @xmath24 . in the first matching step ,
the contributions from the hard region are absorbed into wilson coefficients , while the remaining contributions are represented by fields in the effective theory .
the version of the effective theory with fields scaling in the above way is often called scet@xmath144 . note that , while in many applications the soft scale is associated with nonperturbative physics , in our case @xmath145 is of order the mass of the hadronic final - state jet and assumed to be in the perturbative domain . as explained above ,
the kinematics of the drell - yan process near threshold is such that the hadronic final state is made up of soft partons ; the virtuality of the hard - collinear partons is parametrically larger than the jet mass : @xmath146 .
the potential contribution from glauber gluons @xcite presents a challenge when deriving the standard factorization theorem of the drell - yan cross section into pdfs convoluted with hard partonic cross sections @xcite .
however , the glauber region does not present an additional difficulty for the factorization of the hard - scattering coefficient @xmath147 in the threshold region , which we will derive below .
the hard - scattering coefficient is given by the on - shell @xmath148 cross section , so that for its factorization only glauber exchanges involving active quarks need to be studied .
the momenta of the glauber gluons relevant for the threshold region scale as @xmath149 , and these gluons are thus not allowed in the hadronic final state , which is assumed to have @xmath150 .
writing down the loop integrals for glauber exchanges between the two inital - state quarks , one finds that they are scaleless and vanish in dimensional regularization . in more physical terms , this means that the glauber contribution is accounted for already by the standard momentum regions introduced above .    in order not to clutter the notation
, we consider for the moment a single quark flavor with charge @xmath151 .
we thus match the current @xmath152 onto an effective current operator in scet containing a hard - collinear quark and an anti - hard - collinear anti - quark @xcite . at leading power in @xmath153 ,
this yields @xmath154 where @xmath155 and @xmath156 are the usual hard - collinear ( along the @xmath157-direction ) and anti - hard - collinear ( along the @xmath158-direction ) wilson lines of scet . note that the matching relation for the ( transverse ) vector current does not include a two - gluon operator at leading power in the effective theory .
such an operator would arise , however , in the corresponding matching relation for a scalar current relevant for higgs production .
mixed quark - gluon operators corresponding to the @xmath159 , @xmath160 , @xmath161 , and @xmath162 scattering channels are power suppressed .
accordingly , the corresponding hard - scattering kernels @xmath163 do not contain terms that are singular in the limit @xmath38 , as is explicitly seen in ( [ cijres ] ) .
the fourier transform of the position - space wilson coefficient @xmath164 is a function of the product of the large light - cone momentum components carried by the quark fields . in our case , @xmath165 is equal to the hard scale set by the mass of the drell - yan pair . as a result
the wilson coefficient is evaluated at the time - like momentum transfer carried by the current , @xmath166 .
the @xmath167 prescription is required since this function has a branch cut along the positive real @xmath168 axis .
note that the same wilson coefficient @xmath169 appears in dis , but evaluated at space - like momentum transfer @xcite .
the coefficient @xmath169 can be determined by on - shell matching ; indeed , it is simply given by the on - shell massless form factor in qcd @xcite .
the drell - yan cross section involves the current squared , so that the hard function in ( [ factform ] ) is given by @xmath170 the expression for the coefficient @xmath169 up to @xmath39 was derived in @xcite and can be found in appendix  [ appendix : rgfunctions ] .    at leading power in @xmath153 , only the @xmath171 component of the soft gluon field couples to the hard - collinear fields .
these eikonal interactions can be represented by wilson lines . in the effective theory
this is achieved by redefining the hard - collinear fields as @xcite @xmath172 which implies @xmath173 . here
@xmath174 is a soft wilson line along the @xmath158 light - cone . the same redefinition , but with @xmath158 and @xmath157 interchanged , decouples the soft gluon field also from the anti - hard - collinear fields . as a result ,
the current operator @xmath175 splits into three parts , which no longer interact with each other . in the same way ,
the matrix element for the drell - yan process factorizes in the form @xmath176 to obtain this expression , we have fierz rearranged the fermion fields and have averaged over the color of the external states .
we have simplified the dirac algebra making use of the projection properties @xmath177 and @xmath178 of the hard - collinear fermion fields ( and likewise @xmath179 and @xmath180 for the anti - hard - collinear fields ) .
also , we have neglected the power - suppressed dependence of the hard - collinear matrix element on @xmath181 and @xmath182 ( @xmath183 and @xmath182 for the anti - hard - collinear matrix element ) , using the fact that up to power corrections the incoming partons fly along the beam axis . in more technical terms , we have multi - pole expanded the corresponding fields to leading power @xcite .
the soft matrix element @xmath184 ( not to be confused with the hard - collinear wilson lines ) is a closed wilson loop formed from the product of the soft wilson lines in the two currents , @xmath185\ ,     { \bf t } \big [ s_{\bar n}^\dagger(0 ) s_n(0 ) \big]|0 \rangle \,,\ ] ] where the trace is over color indices , and the operator @xmath186 makes explicit that the wilson lines in the complex conjugate current @xmath187 are anti - time - ordered . a detailed discussion of the operator ordering is given in appendix  [ sec : keldysh ] .
it will become evident in section  [ sec : kin ] that the function @xmath184 is closely related to the soft function for the drell - yan process .    in the second matching step we lower the renormalization scale below the soft scale and integrate out the soft fields .
because the soft scale is in the short - distance domain , this simply amounts to a perturbative calculation of the wilson loop ( [ softmatrixelement ] ) .    to make contact with the standard treatment , where the factorization theorem is derived within perturbative qcd
, we note that after the decoupling transformation ( [ decoupl ] ) the leading - power lagrangian in each of the three sectors of the effective theory is completely equivalent to the qcd lagrangian .
we can thus equally well evaluate the soft function with qcd wilson lines instead of soft wilson lines .
the same is true for the hard - collinear matrix elements , where one can replace @xmath188\ ,     \psi(0)|n_1\rangle \,,\ ] ] where @xmath189 is the qcd quark field and @xmath190 $ ] a qcd wilson line of @xmath191 gluons extending from @xmath192 to @xmath183 .
we have simplified the effective - theory treatment by restricting ourselves to on - shell partonic instead of hadronic matrix elements , which has the advantage that we do not need to discuss momentum regions that scale with powers of @xmath132 .
, a full factorization analysis ( including messenger fields @xcite ) was performed in  @xcite . ]
these momentum regions can not enter the hard or the soft functions , and given that the drell - yan cross section is known to factorize for @xmath193 , we are guaranteed that their presence will not spoil factorization . the restriction to partonic matrix elements
can thus safely be dropped .
the hadronic matrix elements of the form ( [ collinearmatrixelemt ] ) are the usual pdfs @xcite @xmath194\,\psi(0)\,|n(p)\rangle \,.\ ] ]    inserting the factorized form of the matrix element ( [ eq : factor ] ) into expression ( [ textbook ] ) for the cross section , we obtain @xmath195 \nonumber\\     & & \times \int\frac{d^3\vec{q}}{(2\pi)^3 2q^0}\ ,      \frac{1}{2\pi } \int d^4x\ ,      e^{i(x_1 p_{1-}+x_2 p_{2+}-q)\cdot x}\,\hat w_{\rm dy}(x ) \,.\end{aligned}\ ] ] from now on we will drop the @xmath167 prescription in the argument of @xmath169 .      to proceed , we need to study the kinematics in the partonic threshold region , where @xmath38 . with the expansion parameter @xmath141
, it is easy to show that in the parton center - of - mass system @xmath196 it follows that @xmath197 is parametrically larger than the spatial components of @xmath198 .
this observation implies that the term @xmath197 in the denominator in ( [ sigma1 ] ) is independent of @xmath199 at leading power , in which case performing the integral over @xmath200 yields @xmath201 .
it also follows that the rapidity of the drell - yan pair in the parton center - of - mass system vanishes at leading power , @xmath202 which explains why the rapidity in the laboratory frame , @xmath203 , is determined by the ratio @xmath204 at leading power in @xmath26 , as can be seen from ( [ xi ] ) .
finally , we need that @xmath205    using these results , we obtain at leading power in the @xmath38 limit the expression @xmath206 \nonumber\\     & & \times \frac{1}{2\sqrt{\hat s } }      \int\frac{dx^0}{2\pi}\,e^{i\sqrt{\hat s}(1-z)x^0/2}\ ,
\hat w_{\rm dy}(x^0,\vec{x}=0,\muf ) \,,\end{aligned}\ ] ] where the integration region over the parton momentum fractions is such that @xmath94 , so that @xmath207 . in the final step
we introduce the fourier transform of the position - space wilson loop at time - like separation via @xmath208 this wilson loop plays the role of the jet function in dis in the sense that it describes the properties of the hadronic final state .
it has been introduced previously in @xcite .
the drell - yan cross section now takes the form @xmath209
\nonumber\\     & & \times \sqrt{\hat{s}}\,w_{\rm
dy}(\sqrt{\hat s}\,(1-z),\muf ) \,,\end{aligned}\ ] ] and from comparison with ( [ dsigdm2 ] ) and ( [ factform ] ) we identify the soft function as @xmath210 we should mention that at leading power in @xmath26 the argument of the soft function could be simplified as @xmath211 ; however , since the exact expressions ( [ cijres ] ) for the hard - scattering kernels at nlo contain the logarithm @xmath212 , we prefer to keep the argument of the soft function in the form written above .
at this point we have identified the two components in the factorized expression ( [ factform ] ) for the hard - scattering coefficient @xmath128 with field - theoretic objects defined in terms of operator matrix elements .
the resummation of threshold logarithms arising in the @xmath38 region can now be accomplished by solving the rg evolution equations obeyed by these quantities .
the evolution equation for the hard matching coefficient @xmath169 evaluated at time - like momentum transfer and its solution can be obtained from the corresponding results valid for space - like momentum transfer @xcite by analytic continuation .
this leads to @xmath213 c_v(-m^2-i\epsilon,\mu ) \,.\ ] ] we have reinserted the @xmath167 regulator , which determines the sign of the imaginary part of the anomalous dimension .
the appearance of the logarithm and its coefficient , the cusp anomalous dimension @xmath214 @xcite , can be explained using arguments presented in @xcite .
this term in the evolution equation is associated with sudakov double logarithms .
the remaining term , @xmath215 , accounts for single - logarithmic evolution .
the exact solution to ( [ gammav ] ) is @xmath216     \left ( \frac{m^2}{\muh^2 } \right)^{-a_\gamma(\muh,\muf)}\ ,     c_v(-m^2,\muh ) \,,\ ] ] where @xmath74 is a hard matching scale , at which the value of @xmath169 is calculated using fixed - order perturbation theory .
note that the wilson coefficient at time - like momentum transfer is a complex quantity .
the sudakov exponent @xmath100 and the exponents @xmath217 are given by @xcite @xmath218 and similarly for the function @xmath219 .
the perturbative expansions of the anomalous dimensions and the resulting expressions for the evolution functions valid at nnlo in rg - improved perturbation theory are collected in appendix  [ appendix : rgfunctions ] .      in order to derive the evolution equation for the soft function
it is important to have a consistent definition of the threshold region . to this end
we consider the limit @xmath220 in ( [ sigma3 ] ) , which implies @xmath38 , since @xmath103 .
the condition @xmath94 then implies that @xmath221 , in which case the dglap evolution for the pdfs can be written in the simplified form @xmath222_+ } + \dots \,,\ ] ] which is obtained by expanding the altarelli - parisi splitting function as @xmath223_+ }     + 2\gamma^\phi(\alpha_s)\,\delta(1-z ) + \dots \,.\ ] ] this asymptotic form of the splitting function holds to all orders in perturbation theory @xcite . the cusp anomalous dimension and the coefficient @xmath224
have been calculated at three - loop order @xcite .
they are given in appendix  [ appendix : rgfunctions ] .
requiring that the drell - yan cross section in the threshold region be independent of the arbitrary factorization scale @xmath54 , we find that the momentum - space wilson loop obeys the integro - differential evolution equation @xmath225 w_{\rm dy}(\omega,\mu ) \nonumber\\     & & \mbox{}- 4\gamma_{\rm cusp}(\alpha_s ) \int_0^{\omega}\!d\omega'\ ,      \frac{w_{\rm dy}(\omega',\mu)-w_{\rm dy}(\omega,\mu ) }          { \omega-\omega ' } \,,\end{aligned}\ ] ] where @xmath226 curiously , this quantity starts only at two - loop order .
the analogous relation in the case of dis is @xmath227 , where the anomalous dimension @xmath228 of the dis jet function is defined in analogy with @xmath229 in ( [ jrge ] ) . combining the two relations we obtain @xmath230 , which is an exact relation that links the anomalous dimensions of the soft function in drell - yan production and the dis jet function to the anomalous dimension of the vector current .    the exact solution to the evolution equation ( [ jrge ] )
can be written in the form @xcite
@xmath231     \widetilde s_{\rm dy}(\partial_\eta,\mui)\ ,     \frac{1}{\omega }     \left ( \frac{\omega}{\mui } \right)^{2\eta}\ ,     \frac{e^{-2\gamma_e\eta}}{\gamma(2\eta ) } \,,\ ] ] where @xmath232 denotes a derivative with respect to an auxiliary parameter @xmath233 , which is then identified with @xmath234 .
this result is well defined for @xmath235 .
the solution for negative @xmath233 is obtained by analytic continuation .
for instance , to obtain the result for @xmath236 we use the identity @xmath237 where @xmath238 is a smooth test function . for @xmath239
additional subtractions are required .
the function @xmath240 is obtained from the momentum - space wilson loop by the laplace transformation @xcite @xmath241 the position - space wilson loop on the right - hand side of ( [ wdydef ] ) can be shown to have the functional form @xmath242 @xcite , and a straightforward calculation shows that the function @xmath243 can be expressed in terms of @xmath244 as @xmath245 we can then use the explicit two - loop expression for the position - space wilson loop obtained in @xcite to compute the function @xmath240 at two - loop order .
the result is presented in appendix  [ appendix : rgfunctions ] .
we are now in a position to derive the rg - resummed expression for the hard - scattering coefficient @xmath128 in ( [ factform ] ) , which is given by the product of the solutions ( [ cvsol ] ) and ( [ sonice ] ) for the hard and soft functions . the result can be simplified by eliminating the anomalous dimension @xmath229 using ( [ gw ] ) , and combining the sudakov exponents using the relation @xmath246 we find @xmath247 where @xmath234 , and we have defined the evolution function @xmath248 .\ ] ] as before , equation  ( [ cfinal ] ) is valid for @xmath235 ( @xmath249 ) . for negative @xmath233 ( @xmath250 ) , integrals of @xmath251 with test functions @xmath252 must be defined using a subtraction at @xmath253 and analytic continuation in @xmath233 .
we emphasize that the result ( [ cfinal ] ) is formally independent of the scales @xmath254 and @xmath30 , at which the matching conditions for the hard and soft functions are evaluated . on the other hand ,
the hard - scattering kernel @xmath118 does depend on the factorization scale @xmath54 , at which the pdfs are renormalized . in practice , a residual dependence on the matching scales
arises when the perturbative expansions of the matching coefficients and anomalous dimensions are truncated , and this dependence can be used to estimate the remaining perturbative uncertainties .
setting the three scales @xmath254 , @xmath30 , and @xmath54 equal to each other in the resummed expression ( [ cfinal ] ) , one can readily reproduce the leading singular terms for @xmath38 in the fixed - order perturbative qcd expression for the hard - scattering kernel . in this way we have obtained the two - loop corrections in ( [ llterms ] ) .
ccrccc rg - impr .
pt & log .  approx .
& accuracy @xmath255 & @xmath214 & @xmath215 , @xmath224 & @xmath169 , @xmath240 +   +  & ll & @xmath256 & 1-loop & tree - level & tree - level + lo & nll & @xmath257 & 2-loop & 1-loop & tree - level + nlo & nnll & @xmath258 & 3-loop & 2-loop & 1-loop + nnlo & nnnll & @xmath259 & 4-loop & 3-loop & 2-loop +    the final expression ( [ cfinal ] ) for the hard - scattering kernel can be evaluated at any desired order in resummed perturbation theory .
table  [ tab : counting ] shows what is required to obtain different levels of accuracy . in this work
we adopt the counting scheme of rg - improved perturbation theory , where at lo one includes all @xmath260 terms , at nlo one includes all @xmath261 terms , etc .
the large logarithm @xmath262 is counted like @xmath263 . in the literature on threshold resummation the alternative notation n@xmath264ll
is often used instead of n@xmath265lo .
the leading logarithmic ( ll ) approximation is listed only for completeness , as it misses some @xmath260 terms .    in the following section we will perform a detailed numerical analysis of the drell - yan cross section and rapidity distribution . in most cases of phenomenological relevance
the invariant mass of the drell - yan pair will be small compared with the center - of - mass energy , i.e.  @xmath266 . nevertheless , it is interesting to briefly consider the limit @xmath29 , in which the need for threshold resummation is justified parametrically ( see the discussion in section  [ sec : pqcd ] ) . in this case
the convolution integrals in formula ( [ dsigdm2 ] ) for the drell - yan cross section can be performed analytically if a reasonably simple model for the pdfs near the endpoint is adopted .
we parameterize the behavior near @xmath267 as @xmath268 \,,\ ] ] and similarly for the anti - quark distribution .
it then follows that at leading power in @xmath269 the parton luminosity function defined in ( [ ffdef ] ) is given by @xmath270 furthermore , the sum over flavors is dominated by the terms with the smallest exponent @xmath271 .
for example , in @xmath1 collisions at one finds @xmath272 , @xmath273 , and @xmath274 at @xmath275gev @xcite .
these exponents increase by a flavor - independent amount @xmath276 when @xmath54 is raised to larger values , e.g. , by about 0.4 at @xmath105gev @xcite .
it follows that the leading behavior near the endpoint is due to the down - quark contribution .
introducing the drell - yan @xmath277-factor as the ratio @xmath278 we then obtain at leading power @xmath279 this may be compared with the @xmath277-factor for dis at large bjorken @xmath0 and momentum transfer    @xmath280 , which has been derived in @xcite and reads @xmath281 \nonumber\\     & & \times ( 1-x)^\eta\,\ ,      \widetilde j_{\rm dis}\big ( \ln\frac{q^2(1-x)}{\mui^2 }       + \partial_\eta,\mui \big)\ ,      \frac{e^{-\gamma_e\eta}\,\gamma(1+b_u)}{\gamma(1+b_u+\eta ) } \,,\end{aligned}\ ] ] where in this case the dominant contribution comes from the valence up - quark , which has the smallest @xmath282 parameter , @xmath283 at @xmath284gev @xcite .
obviously the structure of threshold logarithms is very similar in the two cases once we consider equal hard scales ( @xmath285 ) and compare the small parameter @xmath286 in dis near the endpoint with the small parameter @xmath287 in drell - yan production .
differences arise from the following facts : ( i ) drell - yan production has time - like kinematics , whereas dis probes the nucleon at space - like momentum transfer .
this gives rise to a difference in the hard matching coefficients @xmath288 starting at one - loop order .
( ii ) the soft / jet functions in the two cases are different starting from one - loop order .
( iii ) the drell - yan cross section involves a convolution with two pdfs , whereas in dis a single parton density appears .
this explains the different coefficients in front of @xmath289 and @xmath233 in the expressions for the @xmath277-factors .
( iv ) for the same reason , the resulting convolution integrals over the pdfs give rise to different expressions involving the @xmath282 and @xmath290 exponents .
there is one more important piece of information that we can extract from the result ( [ kdy ] ) for the drell - yan @xmath277-factor . as we have seen , the exponents @xmath282 and @xmath290 take rather large values .
therefore , the arguments of the @xmath291-functions in ( [ kdy ] ) contain the large quantity @xmath292 .
it is straightforward to show that the derivative with respect to @xmath233 in the argument of the soft function has the effect of changing the argument of the logarithm as follows : @xmath293 up to @xmath260 factors .
it follows that a proper choice for the soft matching scale near @xmath29 is @xmath294 which is an order of magnitude less than the naive choice @xmath31 .
the fact that the fall - off of the parton densities strongly favors the large-@xmath62 region leads to a strong additional suppression of the effective soft scale .
in this section we perform a detailed numerical analysis of our results .
one of our goals is to study to what extent threshold resummation is important ( or even justified ) in processes where the invariant mass of the drell - yan pair is not very close to the center - of - mass energy .
we have seen in section  [ sec : pqcd ] that for very small values of the ratio @xmath12 the threshold contributions are not parametrically enhanced . even though empirically these terms still give rise to the dominant contributions to the cross section
, there is no need to perform a resummation of the threshold logarithms . on the other hand ,
the result ( [ magicmui ] ) derived in the previous section shows that near the true endpoint the threshold logarithms are enhanced by two effects : the kinematic restriction of the @xmath62-integral to the small interval between @xmath36 and 1 , as well as the dynamical enhancement from the strong fall - off of the parton densities . as a result ,
the appropriate value of the soft matching scale is an order of magnitude smaller than the naive choice @xmath31 . from our discussion
so far it is not obvious how to interpolate between these two extreme cases .
we now perform a detailed numerical study to assess the importance of resummation at intermediate values of @xmath36 .
we begin with a discussion of the proper choice of the matching scales @xmath254 and @xmath30 , using the convergence of the perturbative expansions of the matching coefficients @xmath169 and @xmath240 in the resummed hard - scattering kernel ( [ final ] ) as the primary guiding principle .
while it is obvious that the hard matching scale should be chosen of order @xmath7 , the choice of the soft scale is more problematic .
naively , based on the structure of the result ( [ cfinal ] ) one would expect that @xmath295 should in some sense be identified with the  average " value of @xmath296 $ ] . unfortunately , however , the distribution in the variable @xmath62 is both singular at @xmath253 and not positive definite , so that it does not lend itself to a probabilistic interpretation .
a simple way to avoid large logarithmic contributions would be to make the scale choice @xmath297 $ ] inside the @xmath62-integral .
however , this would lead to landau - pole singularities in the integrand and hence would upset the proper scale separation that is at the heart of our approach .    [ ]
@xmath298 [ ] @xmath299 [ ] @xmath300 [ ] @xmath301   dependence of the first two expansion coefficients in the perturbative series for the hard function on the matching scale @xmath254.,title="fig:",scaledwidth=50.0% ]    our approach to threshold resummation based on effective field - theory methods applied directly in momentum space provides a natural way of resolving this question .
the matching scales @xmath254 and @xmath30 should be chosen such that the perturbative expansions of the wilson coefficient functions @xmath169 and @xmath240 in ( [ cfinal ] ) are well behaved . for the case of the soft function
, this criterion should be applied _
after _ the integration over @xmath62 in ( [ final ] ) has been performed .
this is the essence of effective field theory : one performs matching calculations at scales where these calculations can be done in fixed - order perturbation theory , and use the renormalization group to perform the evolution (  running " ) between the different matching scales .
[ ] @xmath30[gev ] [ b]@xmath6 correction   relative contributions to the drell - yan cross section @xmath8 at @xmath302gev arising from the one - loop ( left ) and two - loop ( right ) corrections to the soft function @xmath240 , as a function of the soft matching scale @xmath30 .
the curves are labeled by the corresponding values of @xmath12.,title="fig:",scaledwidth=45.0% ] [ ] @xmath303 correction   relative contributions to the drell - yan cross section @xmath8 at @xmath302gev arising from the one - loop ( left ) and two - loop ( right ) corrections to the soft function @xmath240 , as a function of the soft matching scale @xmath30 .
the curves are labeled by the corresponding values of @xmath12.,title="fig:",scaledwidth=46.0% ]    we begin by applying this criterion to the hard function @xmath304^n .\ ] ] figure  [ fig : muhchoice ] shows the dependence of the expansion coefficients @xmath300 and @xmath301 on the ratio @xmath298 .
the one - loop coefficient @xmath300 vanishes for @xmath305 and @xmath306 .
the second solution is in a region of very small @xmath254 , where the expansion coefficients vary strongly and where the @xmath307 and @xmath308 terms have opposite sign
. we will thus discard it . in the region around the first solution
the two - loop coefficient @xmath301 is stable and positive . in our numerical analysis
we will vary @xmath254 between @xmath7 and @xmath309 , taking @xmath310 as the default choice .
the matching scale @xmath30 must be determined separately for each process , since it is sensitive to the integration range of the @xmath62 variable ( which depends on @xmath36 and @xmath44 ) and to the shape of the pdfs . in figure
[ fig : muichoice1 ] , we plot the relative contributions to the cross section @xmath8 ( normalized to the total cross section ) arising from the one- and two - loop terms in the soft function @xmath240 as a function of @xmath30 .
we choose @xmath302gev and consider different values of @xmath12 between 0.01 and 0.7 .
the plots have been obtained by setting the factorization scale equal to @xmath7 and using mrst2004nnlo parton densities @xcite .
we have checked that virtually indistinguishable results are obtained when @xmath54 is varied by a factor of  2 .
notice that with increasing @xmath36 values the regions where the one- and two - loop contributions are of modest size shift toward lower @xmath30 values .
to be specific , we consider two criteria for a good convergence of the perturbative expansion ( see the left plot in the figure ) :    * starting from a high scale , we determine the value of @xmath30 at which the one - loop correction drops below 15% .
* we choose the value of @xmath30 for which the one - loop contribution is minimal .
note that with either choice the two - loop corrections at the corresponding @xmath30 values are very small , indicating that the first two terms in the perturbation series for the soft function are well behaved for the same choice of scale .
the same analysis can be repeated for different masses of the drell - yan pair .
[ ] @xmath12 [ b]@xmath311   scale - setting results for the soft matching scale @xmath30 for different values of @xmath12 and drell - yan masses @xmath312 , 50 , 20 , 10gev .
the lighter curves correspond to lower scales .
the upper set of curves corresponds to convergence criterion i , the lower one to criterion ii.,title="fig:",scaledwidth=50.0% ]    the resulting values for the soft scale @xmath30 determined using these convergence criteria are shown in figure  [ fig : scaling ] , where we consider the choices @xmath313 , 20 , 50 , 100gev . to a good approximation the curves for the ratio @xmath311 as a function of @xmath36 exhibit scaling , i.e. , they are almost independent of @xmath7 .
we thus obtain a relation of the form @xmath314 .
small scaling violations arise from the scale dependence of the pdfs and of the running coupling in the perturbative expansion of the soft function .
our numerical results are well reproduced by the empirical functions @xmath315 for the two criteria .
their form should not be taken too seriously except to note that for @xmath29 both functions approach @xmath316 , as required by the resummation formula ( [ kdy ] ) for the @xmath277-factor valid near the true endpoint .
indeed , the smallness of the constants in the two cases ( 0.125 and 0.080 , respectively ) is in good agreement with our estimate in ( [ magicmui ] ) .
it results from the dynamical suppression provided by the pdfs . in the opposite limit @xmath317 ,
both forms yield @xmath318 with an @xmath260 constant , as required by the fact that there is only a single physical scale in this case . below we will vary the soft scale between @xmath319 and @xmath320 and , somewhat arbitrarily , use the average of the two scales as the default choice .
the same analysis can be carried out for the double differential decay rate @xmath321 , and it leads to similar results . in general , one finds that at central rapidity ( @xmath322 ) the resulting values for the soft matching scale are higher by about 20% than in the case of @xmath8 , while at larger rapidity they are lower , so that on average , after integration over @xmath44 , the results for the total cross section are reproduced .
we now proceed to study the stability of our resummed expression for the drell - yan cross section , using scale dependence as an estimator of yet unknown higher - order perturbative effects .
we vary the matching scales @xmath254 and @xmath30 about the default values @xmath323 and @xmath324 determined in the previous section .
to simplify comparisons with the literature we adopt the conventional choice @xmath325 for the factorization scale , at which the pdfs are renormalized . from the point of view of effective field theory
it would be more natural to choose a lower value for @xmath54 , given that the cross section is sensitive to physics at scales much below the hard scale @xmath7 .
we will see , however , that our results are very stable with respect to variations of the factorization scale , so that the choice of the default value is not particularly important .
since our focus is on the behavior of different perturbative approximations to the hard - scattering kernel , we use the same set of pdfs ( mrst04nnlo ) throughout the analysis .
our results are shown in the first three plots in figure  [ fig : scales ] , in which we study the drell - yan @xmath277-factor defined in ( [ kdef ] ) for @xmath302gev and various values of @xmath12 . in the calculation of the @xmath277-factor we keep the factorization scale
@xmath325 fixed in the leading - order expression in the denominator , even when @xmath54 is varied in the numerator . for the time being we only include the leading terms in the @xmath38 limit , corresponding to the result ( [ cfinal ] ) . adding the small power - suppressed corrections would not change any of our conclusions .
we observe an excellent convergence of the perturbative expansion for the cross section after resummation .
the bands corresponding to the lo , nlo , and nnlo approximations overlap , and the dependence on the matching scales @xmath254 and @xmath30 becomes negligible beyond lo , indicating that the residual perturbative uncertainty is very small .
[ l ]  @xmath326 [ b]@xmath327 [ ]    [ cols="^,^ " , ]     the fact that the moment - space resummation approach differs from the momentum - space approach by first - order corrections in @xmath26 leads to visible numerical effects . in order to illustrate this
we take the inverse mellin transformation of ( [ sign ] ) to compute the drell - yan @xmath277-factor . for simplicity
, we fit the polynomial form @xmath328 to the parton luminosity function defined in ( [ ffdef ] ) , so that the product of the moments of the pdfs in ( [ sign ] ) can be evaluated in closed form .
the first plot in figure  [ fig : compmoment ] shows that for the same choices of matching scales the moment - space approach leads to somewhat larger results for the @xmath277-factor than the momentum - space approach , and that the convergence is slightly better in the second case .    in the traditional moment - space approach based on expression ( [ gn ] )
one implicitly makes the scale choice @xmath329 , which leads to a landau pole at @xmath330 in the mellin inversion . to see how the results obtained with this choice compare with our default results we adopt the so - called minimal prescription @xcite , which amounts to choosing an inversion contour in the @xmath53-plane that does not include this pole .
when this is done , we find that the moment - space formula ( [ momentspacecn ] ) evaluated with @xmath331 and @xmath329 gives numerical results similar to those found using our default scale choices once we go beyond leading order .
this is illustrated in the right plot in figure  [ fig : compmoment ] .
we have presented a detailed analysis of threshold resummation for drell - yan production . instead of the conventional mellin moment formalism , we have performed the resummation directly in momentum space , using rg evolution .
we have obtained analytic expressions for the resummed hard - scattering kernels that are free of the landau - pole ambiguities inherent in the traditional approach . after deriving the necessary matching coefficients and anomalous dimensions from known perturbative results in the literature
, we have performed the resummation for the cross section and rapidity distribution at nnlo in rg - improved perturbation theory , corresponding to n@xmath9ll accuracy .
the main goal of this work was to address the question to what extent resummation can be phenomenologically relevant given that in all practical applications the true threshold region is experimentally not accessible , because the parton distributions fall off very steeply at large @xmath0 .
it has been argued in the literature that precisely this fall - off forces the drell - yan process to the partonic threshold region , in which large logarithmic corrections arise .
our analysis confirms the existence of this effect and quantifies its importance . in the true endpoint region
, we find that the scale of soft radiation is an order of magnitude smaller than the naive estimate .
the extra suppression factor is given by the sum of the exponents of the fall - off of the quark and anti - quark distributions near @xmath267 .
on the other hand , for very small lepton - pair masses this effect becomes inoperative because the fall - off of the parton distributions is not very steep at low @xmath0 .
an analysis of the convergence of perturbation theory in the intermediate region shows that resummation effects become relevant for @xmath7 larger than about @xmath332 . at the largest measured drell - yan masses resummation effects
can thus be significant . in the case of the experiment e866/nusea ,
a fixed - target experiment with @xmath1 collisions at @xmath14gev , these effects increase the nnlo ( nlo ) cross section by 7% ( 25% ) at @xmath42gev .
we do not confirm the recent claim of large negative resummation effects for the rapidity distribution measured in the e866/nusea experiment at @xmath15gev @xcite .
we find that the fixed - order nnlo threshold contribution is positive and a good approximation to the full nnlo correction .
higher - order terms beyond nnlo turn out to be negligible , so that it is unnecessary to resum them . to check whether the large effect seen in @xcite could be an artifact of the moment - space formalism ,
we have applied our approach also in moment space and showed that it encompasses the traditional moment - space approach .
we find that the two methods give compatible results , and that in both schemes resummation increases the cross section .
the drell - yan cross section is dominated by the contribution of the threshold region even at small masses of the lepton pair .
for example , for @xmath15gev the threshold terms give rise to @xmath333 of the nlo correction to the total cross section and @xmath334 of the nnlo contribution .
if this is a more generic feature of hard cross sections , then threshold resummation techniques provide an efficient way to obtain the dominant part of the higher - order corrections .
perhaps the most important outcome of our analysis is a quantitative understanding of the emergence of an effective physical scale characterizing the soft radiation in the drell - yan process .
this scale is generated through an intricate interplay of dynamical and kinematical effects .
it will be interesting to explore how the same mechanism affects other collider processes such as higgs production at the tevatron and lhc .
furthermore , it would be useful to analyze resummation for @xmath335-production , which is phenomenologically relevant , since the tevatron now produces top - quark pairs with quite high invariant masses .
there was some controversy on how to best perform the resummmation in this case @xcite .
it would be interesting to revisit the issue using effective field - theory methods .
we are grateful to lance dixon , keith ellis , frank petriello , and werner vogelsang for useful discussions .
m.n .  likes to thank the fermilab theory group for hospitality and support as a frontier fellow .
the research of t.b .  was supported by the department of energy under grant de - ac02 - 76ch03000 .
the research of g.x .  was supported by the national science foundation under grant phy-0355005 .
fermilab is operated by the fermi research alliance under contract with the u.s .
department of energy .
the leading singular terms in the fixed - order perturbative expressions for the hard - scattering kernels @xmath163 can be derived by using the perturbative expansions of the hard and soft functions in ( [ factform ] ) . in this way
, we obtain for the coefficients @xmath87 in ( [ llterms ] ) @xmath336 \nonumber\\     & & \mbox{}+ 16\zeta_3 \left [ \frac{1}{1-z }
\right]_+      + \left ( - l^2 + 3l - 8 - \frac{2\pi^2}{3 } \right )      \left [ \frac{l_z}{1-z } \right]_+      + \left [ \frac{l_z^3}{1-z } \right]_+ \ , , \nonumber\\     p_a(z ) & = & \delta(1-z ) \left [ - \frac{11}{16}\,l^2       + \left ( \frac{193}{48 } - \frac{11\pi^2}{36 }       - \frac{3\zeta_3}{2 } \right ) l      - \frac{1535}{192 } + \frac{47\pi^2}{36 } - \frac{23\pi^4}{720 }       + \frac{43\zeta_3}{12 } \right ] \nonumber\\     & & \mbox{}+ \left ( - \frac{101}{27 } + \frac{11\pi^2}{18 }       + \frac{7\zeta_3}{2 } \right ) \left [ \frac{1}{1-z } \right]_+      + \left ( \frac{67}{18 } - \frac{\pi^2}{6 } \right )      \left [ \frac{l_z}{1-z } \right]_+      - \frac{11}{12 } \left [ \frac{l_z^2}{1-z } \right]_+ \ , , \nonumber\\     p_f(z ) & = & \delta(1-z ) \left [ \frac14\,l^2      + \left ( - \frac{17}{12 } + \frac{\pi^2}{9 } \right ) l       + \frac{127}{48 } - \frac{4\pi^2}{9 } + \frac{\zeta_3}{3 } \right ]      \nonumber\\     & & \mbox{}+ \left ( \frac{28}{27 } - \frac{2\pi^2}{9 } \right )      \left [ \frac{1}{1-z } \right]_+      - \frac{10}{9 } \left [ \frac{l_z}{1-z } \right]_+      + \frac13 \left [ \frac{l_z^2}{1-z } \right]_+ \,.\end{aligned}\ ] ] as before , we use the abbreviations @xmath85 and @xmath86 $ ] .
if desired , the plus distributions involving powers of @xmath88 can be reduced to distributions of the form @xmath337_+$ ] using the identity @xmath338_+
= f(z ) \left [ g(z ) \right]_+ - \delta(1-z ) \int_0 ^ 1\!dz'\,f(z ' )     \left [ g(z ' ) \right]_+ \,.\ ] ]
for completeness we list the perturbative expansions of the various matching coefficients and anomalous dimensions required to evaluate our rg - improved result ( [ cfinal ] ) at nnlo .
the matching conditions for the wilson coefficient @xmath169 evaluated at time - like momentum transfer can be obtained by analytic continuation from the corresponding expression valid at space - like momentum transfer . using the known two - loop result for the on - shell qcd form factor @xcite
, we find @xcite @xmath339 , \ ] ] where @xmath340 and @xmath341 .
this result agrees with the corresponding expression given in @xcite
.    the matching condition for the function @xmath240 can be derived most easily using relation ( [ wonderful ] ) , which relates it to the perturbative expansion of the position - space wilson loop @xmath342 at time - like separation .
we can then use the explicit two - loop expression for the position - space wilson loop obtained in @xcite . to this end , we have obtained the fourth - order coefficient in the expansion of a certain appell hypergeometric function , @xmath343              2+\epsilon,1 - 2\epsilon \end{array}\,\big|\,1,1 \big )     = - \frac{1+\epsilon}{2\epsilon } \left (       1 - \frac{\pi^2}{3}\,\epsilon^2 - 14\zeta_3\epsilon^3       - \frac{5\pi^4}{18}\,\epsilon^4 + \dots \right ) .\ ] ] the resulting two - loop expression for the soft function reads @xmath344 , \ ] ] where @xmath345 note that the @xmath346 terms exponentiate , which is a consequence of the non - abelian exponentiation theorem for wilson loops @xcite .
our result for the function @xmath240 agrees with a corresponding expression entering in the moment - space resummation approach studied in @xcite .
here we list expressions for the anomalous dimensions and the qcd @xmath347-function , quoting all results in the @xmath348 renormalization scheme .
we define the expansion coefficients of the anomalous dimensions and the qcd @xmath347-function as @xmath349 , \end{aligned}\ ] ] and similarly for the other anomalous dimensions .
the expansion of the cusp anomalous dimension @xmath214 to two - loop order was obtained some time ago @xcite , while recently the three - loop coefficient has been calculated in @xcite .
for the four - loop coefficient @xmath350 we use the pad approximant @xmath351 .
the results are @xmath352 ,      \nonumber\\     \gamma_2 & = & 4 c_f \bigg [ c_a^2 \left ( \frac{245}{6 }       - \frac{134\pi^2}{27 }      + \frac{11\pi^4}{45 } + \frac{22}{3}\,\zeta_3 \right )       + c_a t_f n_f   \left ( - \frac{418}{27 } + \frac{40\pi^2}{27 }      - \frac{56}{3}\,\zeta_3 \right ) \nonumber\\     & & \mbox{}+ c_f t_f
n_f \left ( - \frac{55}{3 } + 16\zeta_3 \right )       - \frac{16}{27}\,t_f^2 n_f^2 \bigg ] \ , , \nonumber\\     \gamma_3 & \approx & 7849,~ 4313 , ~ 1553 \quad \mbox{for } \quad      n_f = 3,\,4,\,5 \,.\end{aligned}\ ] ] the anomalous dimension @xmath215 can be determined up to three - loop order from the partial three - loop expression for the on - shell quark form factor in qcd , which has recently been obtained in @xcite .
we find @xmath353 the anomalous dimension @xmath224 is know to three - loop order from the nnlo calculation of the altarelli - parisi splitting functions @xcite .
the expansion coefficients are @xmath354 using these results , one can compute the expansion coefficients for the anomalous dimension @xmath229 of the drell - yan soft function from the relation @xmath355 .
this yields @xmath356 we do not list the three - loop coefficient .
finally , the expansion coefficients for the qcd @xmath347-function to four - loop order are @xmath357 the value of @xmath358 is taken from @xcite and corresponds to @xmath359 and @xmath360 .
we now give the perturbative expansions of the functions @xmath100 and @xmath361 defined in ( [ rgesols ] ) , working consistently at nnlo in rg - improved perturbation theory . at this order
we need to keep terms through @xmath39 in the final expressions .
the resulting expression for @xmath361 is given by @xmath362      \frac{\alpha_s^2(\mu ) - \alpha_s^2(\nu)}{32\pi^2 } + \dots \bigg\ }      \,.\end{aligned}\ ] ] similar expressions with the @xmath363 replaced by the coefficients @xmath364 or @xmath365 hold for the functions @xmath219 and @xmath289 , respectively .
the nnlo expression for the sudakov exponent @xmath100 is more complicated .
it reads @xcite @xmath366      \nonumber\\     & & \mbox{}+ \left ( \frac{\alpha_s(\nu)}{4\pi } \right)^2 \bigg [      \left ( \frac{\beta_1\beta_2}{\beta_0 ^ 2 }       - \frac{\beta_1 ^ 3}{2\beta_0 ^ 3 }      - \frac{\beta_3}{2\beta_0 } + \frac{\beta_1}{\beta_0 }      \left ( \frac{\gamma_2}{\gamma_0 } - \frac{\beta_2}{\beta_0 }      + \frac{\beta_1 ^ 2}{\beta_0 ^ 2 }       - \frac{\beta_1\gamma_1}{\beta_0\gamma_0 } \right )       \frac{r^2}{2 } \right ) \ln r \nonumber\\     & & \hspace{1.0 cm }      \mbox{}+ \left ( \frac{\gamma_3}{\gamma_0 }       - \frac{\beta_3}{\beta_0 }      + \frac{2\beta_1\beta_2}{\beta_0 ^ 2 }       + \frac{\beta_1 ^ 2}{\beta_0 ^ 2 }      \left ( \frac{\gamma_1}{\gamma_0 } - \frac{\beta_1}{\beta_0 } \right )      - \frac{\beta_2\gamma_1}{\beta_0\gamma_0 }      - \frac{\beta_1\gamma_2}{\beta_0\gamma_0 } \right )       \frac{(1-r)^3}{3 } \nonumber\\     & & \hspace{1.0 cm }      \mbox{}+ \left ( \frac{3\beta_3}{4\beta_0 }       - \frac{\gamma_3}{2\gamma_0 } + \frac{\beta_1 ^ 3}{\beta_0 ^ 3 }      - \frac{3\beta_1 ^ 2\gamma_1}{4\beta_0 ^ 2\gamma_0 }      + \frac{\beta_2\gamma_1}{\beta_0\gamma_0 }      + \frac{\beta_1\gamma_2}{4\beta_0\gamma_0 }      - \frac{7\beta_1\beta_2}{4\beta_0 ^ 2 } \right ) ( 1-r)^2 \nonumber\\     & & \hspace{1.0 cm }      \mbox{}+ \left ( \frac{\beta_1\beta_2}{\beta_0 ^ 2 }       - \frac{\beta_3}{\beta_0 }      - \frac{\beta_1 ^ 2\gamma_1}{\beta_0 ^ 2\gamma_0 }      + \frac{\beta_1\gamma_2}{\beta_0\gamma_0 } \right ) \frac{1-r}{2 }      \bigg ] + \dots \bigg\ } \,,\end{aligned}\ ] ] where @xmath367 . whereas the three - loop anomalous dimensions and @xmath347-function are required in ( [ asol ] ) , the expression for @xmath100 also involves the four - loop coefficients @xmath350 and @xmath358 .
to perform the effective theory analysis of the process , we would like to have a path - integral definition of the hadronic quantity of interest . in the case of dis
( and also for inclusive @xmath368-decays ) , one considers the discontinuity of forward matrix elements of time - ordered products of the electroweak currents and studies their factorization properties . the drell - yan cross section can not be written in a similar form .
the reason is that not all cuts of the relevant feynman graphs correspond to the same physical process .
for example , in addition to the contribution from the cut indicated in figure  [ fig : cutdiagram ] , the discontinuity of the same diagram also gets a contribution from a cut through the triangle loop on the left , which describes @xmath369 with a virtual lepton pair .
a path - integral method for the direct evaluation of cut diagrams is the keldysh ( or time - loop ) formalism @xcite .
let us illustrate the method for a scalar field theory with field @xmath372 . instead of the usual action , where one integrates over time from @xmath373
, one considers a path integral with an action @xmath374 where the contour @xmath375 first runs from @xmath376 and then back , with a negative imaginary part , along the contour shown in figure  [ fig : contour ] . when evaluating expectation values of fields , one needs to specify on which part of the contour the fields reside .
we denote the fields living on the first half of the contour by @xmath377 and the ones on the second half by @xmath378 .
the path - integral expectation values of fields correspond to path - ordered vacuum expectation values . since the fields @xmath379 are positioned at earlier points along the contour , they are to the right of all fields @xmath380 .
furthermore path - ordering translates into anti - time - ordering on the second part of the contour , so that we find @xmath381 here @xmath382 denotes anti - time - ordering . note that under a field redefinition @xmath383 , the plus and minus fields transform separately @xmath384
this is relevant when one uses a field redefinition to decouple the soft gluon fields from the quark fields in the current operator .
the keldysh formalism is useful , because it gives a path integral formulation of squared amplitudes , e.g. @xmath385 which is then rewritten as an expectation value of @xmath386 fields using ( [ expectation ] ) .    for our perturbative analysis
, we need the feynman rules to calculate matrix elements of the form ( [ expectation ] ) .
the rules are simple : those for the field @xmath379 are the usual feynman rules and those for @xmath380 are the complex conjugate of the usual rules
. since all interactions are local , there are no vertices involving both @xmath379 and @xmath380 fields .
the only connection is the cut propagator @xmath387 this expression is familiar from the cutkosky rules @xcite used to extract the contribution of a given cut to an ordinary feynman diagram .
[ t]@xmath388 -plane along which the lagrangian is integrated to obtain the action in the keldysh formalism .
( to dampen the oscillatory behavior of the path integrand , the entire contour should further be rotated clockwise by a small angle . ) [ fig : contour],title="fig:",scaledwidth=40.0% ]        s.  d.  drell and t.  m.  yan , phys .
lett .   * 25 * , 316 ( 1970 ) [ erratum - ibid .   * 25 * , 902 ( 1970 ) ] .
g.  altarelli , r.  k.  ellis and g.  martinelli , nucl .
b * 157 * , 461 ( 1979 ) .
r.  hamberg , w.  l.  van neerven and t.  matsuura , nucl .
b * 359 * , 343 ( 1991 ) [ erratum - ibid .
b * 644 * , 403 ( 2002 ) ] .
r.  v.  harlander and w.  b.  kilgore , phys .
lett .   * 88 * , 201801 ( 2002 ) [ hep - ph/0201206 ] . c.  anastasiou , l.  j.  dixon , k.  melnikov and f.  petriello , phys .
lett .   * 91 * , 182002 ( 2003 ) [ hep - ph/0306192 ] . c.  anastasiou , l.  j.  dixon , k.  melnikov and f.  petriello , phys .
d * 69 * , 094008 ( 2004 ) [ hep - ph/0312266 ] .
k.  melnikov and f.  petriello , phys .
lett .   * 96 * , 231803 ( 2006 ) [ hep - ph/0603182 ] .
k.  melnikov and f.  petriello , phys .
d * 74 * , 114017 ( 2006 ) [ hep - ph/0609070 ] .
g.  sterman , nucl .
b * 281 * , 310 ( 1987 ) .
s.  catani and l.  trentadue , nucl .
b * 327 * , 323 ( 1989 ) .
l.  magnea , nucl .
b * 349 * , 703 ( 1991 ) .
g.  p.  korchemsky and g.  marchesini , phys .
b * 313 * , 433 ( 1993 ) .
g.  sterman and w.  vogelsang , jhep * 0102 * , 016 ( 2001 ) [ hep - ph/0011289 ] .
a.  mukherjee and w.  vogelsang , phys .
d * 73 * , 074005 ( 2006 ) [ hep - ph/0601162 ] .
p.  bolzoni , phys .
b * 643 * , 325 ( 2006 ) [ hep - ph/0609073 ] .
v.  ravindran , j.  smith and w.  l.  van neerven , hep - ph/0608308 . v.  ravindran and j.  smith , arxiv:0708.1689 [ hep - ph ] .
j.  c.  webb _ et al .
_ [ nusea collaboration ] , hep - ex/0302019 .
d.  appell , g.  sterman and p.  b.  mackenzie , nucl .
b * 309 * , 259 ( 1988 ) .
s.  catani , m.  l.  mangano and p.  nason , jhep * 9807 * , 024 ( 1998 ) [ hep - ph/9806484 ] .
t.  becher and m.  neubert , phys .
lett .   * 97 * , 082001 ( 2006 ) [ hep - ph/0605050 ] .
t.  becher , m.  neubert and b.  d.  pecjak , jhep * 0701 * , 076 ( 2007 ) [ hep - ph/0607228 ] . c.  w.  bauer , s.  fleming , d.  pirjol and i.  w.  stewart , phys .
d * 63 * , 114020 ( 2001 ) [ hep - ph/0011336 ] . c.  w.  bauer , d.  pirjol and i.  w.  stewart , phys .
d * 65 * , 054022 ( 2002 ) [ hep - ph/0109045 ] .
m.  beneke , a.  p.  chapovsky , m.  diehl and t.  feldmann , nucl .
b * 643 * , 431 ( 2002 ) [ hep - ph/0206152 ] .
v.  n.  gribov and l.  n.  lipatov , sov .
j.  nucl .
phys .   * 15 * , 438 and 675 ( 1972 ) [ yad .
15 * , 781 and 1218 ( 1972 ) ] . y.  l.  dokshitzer , sov .
jetp * 46 * , 641 ( 1977 ) [ zh .  eksp .
73 * , 1216 ( 1977 ) ] .
g.  altarelli and g.  parisi , nucl .
b * 126 * , 298 ( 1977 ) .
the code dyrap we used is available at http://www.slac.stanford.edu/~lance / vrap/.    a.  d.  martin , r.  g.  roberts , w.  j.  stirling and r.  s.  thorne , phys .
b * 604 * , 61 ( 2004 ) [ hep - ph/0410230 ] . c.  w.  bauer , s.  fleming , d.  pirjol , i.  z.  rothstein and i.  w.  stewart , phys .
d * 66 * , 014017 ( 2002 ) [ hep - ph/0202088 ] .
g.  t.  bodwin , s.  j.  brodsky and g.  p.  lepage , phys .
lett .   *
47 * , 1799 ( 1981 ) .
j.  c.  collins , d.  e.  soper and g.  sterman , phys .
b * 109 * , 388 ( 1982 ) .
j.  c.  collins , d.  e.  soper and g.  sterman , nucl .
b * 261 * , 104 ( 1985 ) .
j.  c.  collins , d.  e.  soper and g.  sterman , nucl .
b * 308 * , 833 ( 1988 ) .
a.  idilbi , x.  d.  ji and f.  yuan , nucl .
b * 753 * , 42 ( 2006 ) [ hep - ph/0605068 ] .
j.  s.  schwinger , j.  math .
* 2 * , 407 ( 1961 ) .
l.  v.  keldysh , zh .  eksp .
47 * , 1515 ( 1964 ) [ sov .
jetp * 20 * , 1018 ( 1965 ) ] .
a.  v.  belitsky , phys .
b * 442 * , 307 ( 1998 ) [ hep - ph/9808389 ] .
t.  becher , r.  j.  hill , b.  o.  lange and m.  neubert , phys .
d * 69 * , 034013 ( 2004 ) [ hep - ph/0309227 ] .
a.  v.  manohar , phys .
d * 68 * , 114019 ( 2003 ) [ hep - ph/0309176 ] .
t.  becher , r.  j.  hill and m.  neubert , phys .
d * 69 * , 054017 ( 2004 ) [ hep - ph/0308122 ] .
m.  beneke and t.  feldmann , phys .
b * 553 * , 267 ( 2003 ) [ hep - ph/0211358 ] .
j.  c.  collins and d.  e.  soper , nucl .
b * 194 * , 445 ( 1982 ) .
g.  p.  korchemsky and a.  v.  radyushkin , nucl .
b * 283 * , 342 ( 1987 ) .
i.  a.  korchemskaya and g.  p.  korchemsky , phys .
b * 287 * , 169 ( 1992 ) .
m.  neubert , eur .
j.   c * 40 * , 165 ( 2005 ) [ hep - ph/0408179 ] .
g.  p.  korchemsky and g.  marchesini , nucl .
b * 406 * , 225 ( 1993 ) [ hep - ph/9210281 ] .
s.  moch , j.  a.  m.  vermaseren and a.  vogt , nucl .
b * 688 * , 101 ( 2004 ) [ hep - ph/0403192 ] .
m.  neubert , phys .
d * 72 * , 074025 ( 2005 ) [ hep - ph/0506245 ] .
s.  alekhin , jetp lett .   * 82 * , 628 ( 2005 ) [ pisma zh .
* 82 * , 710 ( 2005 ) ] [ hep - ph/0508248 ] .
a.  d.  martin , r.  g.  roberts , w.  j.  stirling and r.  s.  thorne , phys .
b * 531 * , 216 ( 2002 ) [ hep - ph/0201127 ] .
w.  k.  tung , h.  l.  lai , a.  belyaev , j.  pumplin , d.  stump and c.  p.  yuan , jhep * 0702 * , 053 ( 2007 ) [ hep - ph/0611254 ] .
t.  becher and m.  neubert , phys .
b * 637 * , 251 ( 2006 ) [ hep - ph/0603140 ] .
s.  moch and a.  vogt , phys .
b * 631 * , 48 ( 2005 ) [ hep - ph/0508265 ] .
e.  laenen and l.  magnea , phys .
b * 632 * , 270 ( 2006 ) [ hep - ph/0508284 ] .
s.  moch , j.  a.  m.  vermaseren and a.  vogt , nucl .
b * 726 * , 317 ( 2005 ) [ hep - ph/0506288 ] .
v.  ravindran , nucl .
b * 746 * , 58 ( 2006 ) [ hep - ph/0512249 ] .
v.  ravindran , nucl .
b * 752 * , 173 ( 2006 ) [ hep - ph/0603041 ] .
m.  beneke and v.  m.  braun , nucl .
b * 454 * , 253 ( 1995 ) [ hep - ph/9506452 ] .
s.  catani , m.  l.  mangano , p.  nason and l.  trentadue , nucl .
b * 478 * , 273 ( 1996 ) [ hep - ph/9604351 ] .
t.  matsuura and w.  l.  van neerven , z.  phys .
c * 38 * , 623 ( 1988 ) .
t.  matsuura , s.  c.  van der marck and w.  l.  van neerven , nucl .
b * 319 * , 570 ( 1989 ) .
t.  gehrmann , t.  huber and d.  maitre , phys .
b * 622 * , 295 ( 2005 ) [ hep - ph/0507061 ]
. s.  moch , j.  a.  m.  vermaseren and a.  vogt , jhep * 0508 * , 049 ( 2005 ) [ hep - ph/0507039 ] .
j.  g.  m.  gatheral , phys .
b * 133 * , 90 ( 1983 ) .
j.  frenkel and j.  c.  taylor , nucl .
b * 246 * , 231 ( 1984 ) .
t.  van ritbergen , j.  a.  m.  vermaseren and s.  a.  larin , phys .
b * 400 * , 379 ( 1997 ) [ hep - ph/9701390 ] .
r.  e.  cutkosky , j.  math .
* 1 * , 429 ( 1960 ) . s.  catani , m.  l.  mangano , p.  nason and l.  trentadue , phys
.  lett .
b * 378 * , 329 ( 1996 ) [ hep - ph/9602208 ] .
e.  l.  berger and h.  contopanagos , phys .
d * 54 * , 3085 ( 1996 ) [ hep - ph/9603326 ] .
e.  l.  berger and h.  contopanagos , phys .
d * 57 * , 253 ( 1998 ) [ hep - ph/9706206 ] .