many of the most accurate heavy - quark mass determinations are related to the spectral functions of the heavy - quark vector current , which can be measured in @xmath0 collisions .
the energy region near the pair production threshold is particularly sensitive to the mass value .
the bottom and charm masses are usually inferred from suitable averages of the pair production cross section .
looking to the future , the measurement of the top quark pair production cross section in the threshold region at the planned international linear collider ( ilc ) would lead to a very precise knowledge of the top mass directly from the energy dependence of the cross section , even though the toponium resonances are smeared out due to the large top quark width  @xcite . to put this into perspective
, the current top mass value from direct production at the fermilab tevatron is @xmath5gev @xcite , and @xmath6gev from the large hadron collider ( lhc ) run at @xmath7tev centre - of - mass energy  @xcite .
a less precise value can be obtained from the total production cross section as illustrated for example in @xcite .
further reduction of the uncertainty below the one of the tevatron result is very difficult due to the complicated theoretical systematics of top jet mass reconstruction at hadron colliders .
the above numbers should be compared to a precision of @xmath8mev that can be achieved experimentally  @xcite from the @xmath9 threshold scan at the ilc . in particular , while it is not evident which renormalized top mass parameter is determined with the quoted accuracy at the tevatron and lhc , the threshold cross section in @xmath0 collisions provides an observable that can be unambiguously related to a particular top mass definition .    aside from determining a fundamental parameter of the standard model ( sm ) ,
accurate top mass measurements are of interest for extrapolating the sm or its tev scale extension to higher energies , either in the context of `` electroweak precision tests '' or a theory of the yukawa couplings .
furthermore , the recent discovery of a higgs boson @xcite with a mass of about 125@xmath10gev if interpreted as the sm higgs boson determines the higgs self - coupling and leads to the conclusion that the sm vacuum becomes metastable at scales above @xmath11gev  @xcite .
the precise value of this scale turns out to depend very sensitively on the value of the top quark mass . finally , the absence of any hint on physics beyond the sm in high - energy collisions at lhc has renewed the interest in performing precision measurements of properties of sm particles , including the top quark mass , width and yukawa coupling .
measurements at the top pair production threshold are uniquely suited for this purpose .
the challenge is thus to calculate the heavy - quark to zero in the top - quark expressions . ]
spectral functions precisely in the threshold region .
this kinematic region is characterized by two features , which make the theoretical calculation of qcd corrections rather different from standard loop calculations : the small three - velocity @xmath12 of the heavy quarks , which allows to expand feynman diagrams in @xmath12 rather than calculate them exactly , and the strong colour - coulomb force , which on the other hand requires certain diagrams to be summed to all orders in the strong coupling @xmath13 similar to bound - state calculations in quantum electrodynamics .
the expansion of the cross section relative to the ultra - relativistic point - particle cross section @xmath14 is then organized as @xmath15 where the overall factor of @xmath12 arises from the phase - space of the two produced massive particles , and the order of the various terms is indicated explicitly . to perform the expansion and the required summation of perturbation theory to all orders beyond the next - to - leading order ( nlo ) , non - relativistic effective field theory  @xcite and
the threshold expansion of feynman diagrams  @xcite are the methods of choice .
the cross section is then obtained from the expression - boson exchange . see section  [ sec : eft ] for the full expression .
] @xmath16 g(e ) + \ldots\right ) \right],\ ] ] where @xmath17 , @xmath18 denote certain relativistic matching corrections , @xmath19 , @xmath20 the number of colours , and @xmath21 represents a two - point function of heavy - quark currents in the non - relativistic effective theory .
the purpose of this paper is to present results of the part of the third - order ( nnnlo ) qcd corrections to the heavy quark anti - quark pair production cross section near threshold related to the correlation function @xmath21 , which contains the all - order summation . since this concludes the non - relativistic third - order calculation , we also present details of the methods and calculations that have been used but not given in earlier publications .    for the top pair production threshold the leading order ( lo ) and next - to - leading order ( nlo ) approximations to the cross section ( [ systematics ] )
have been examined long ago  @xcite .
several other aspects of the threshold such as top momentum distributions and polarization have been computed at this order  @xcite . beyond nlo , matching the non - relativistic approximation to qcd is non - trivial , because the separation between relativistic and non - relativistic physics is no longer unambiguous . a consistent field theoretical approach based on non - relativistic effective qcd is now required . the second - order ( nnlo ) qcd corrections to the total pair production cross section
have been computed in this framework about ten years ago  @xcite and turned out to be surprisingly large even for top quarks . while some of the large corrections can be understood as being due to mass renormalization  @xcite , and can be avoided by a change of renormalization convention , there remains an apparently slow convergence of successive approximations to the normalization of the cross section , which necessitates the calculation of the nnnlo term . an alternative approach that sums logarithms of @xmath12
has also been pursued , and an improvement of convergence has been found in a ( still partially incomplete ) next - to - next - to - leading logarithmic ( nnll ) approximation  @xcite .
nonetheless , the nnnlo non - logarithmic terms not included in the nnll approximation are required to be certain that the theoretical calculation is sufficiently accurate for the proposed mass measurement at the ilc .
this is the main motivation for the present work . over the past years
a significant number of results relevant to a nnnlo calculation or partial results for third - order quarkonium energy levels and wave - functions at the origin have already appeared  @xcite . in the following we
summarize the status of the nnnlo top cross section calculation , the contribution of the present paper , and what remains to be done to complete the nnnlo calculation .
while the resummed cross section is calculated in an effective field theory ( eft ) , a number of matching calculations has to be performed to guarantee that the eft reproduces qcd to the required accuracy .
this is done in two steps .
hard matching ( scale @xmath22 ) yields the coefficients of the non - relativistic qcd ( nrqcd ) interactions and heavy - quark currents ; soft matching ( scale @xmath23 ) the quark - anti - quark potentials . at nnnlo
the coefficients of several subleading nrqcd interactions must be determined with one - loop precision .
this calculation has been performed in  @xcite . however , as shall be explained , the calculation of the cross section requires the @xmath24 terms of the coefficient functions .
we therefore repeated the nrqcd matching calculation , confirm the results of  @xcite and provide the expressions for the @xmath24 terms in this paper ( see also @xcite ) .
hard matching of the non - relativistic currents at nnnlo is needed at the one - loop level for the sub - leading currents and at the three - loop level for the leading current .
the former are known  @xcite , and will be rederived in the present paper , but only the fermionic contributions to the three - loop matching coefficient @xmath25 have been calculated up to now  @xcite . as discussed in  @xcite , the incomplete result for @xmath25 presently constitutes the most important missing nnnlo contribution to the cross section .    as concerns soft matching , the potentials of order @xmath26 must be determined with one - loop precision , the @xmath27 potential with two - loop precision , and the @xmath28 coulomb potential at three loops , since @xmath29 counts as the bohr radius @xmath30 . except for the coulomb potential ,
the coefficients of the other potentials have been calculated in  @xcite , but again these results are not sufficient for the cross section calculations , since the @xmath24 terms of all these potentials are needed .
we repeated the calculation of the one - loop @xmath26 potentials , confirm the results of  @xcite and provide the expressions for the @xmath24 terms in this paper ( see also  @xcite ) .
the @xmath24 term of the two - loop @xmath27 potential coefficient @xmath31 remains presently unknown . as concerns the three - loop coulomb potential , the fermionic contributions to the three - loop coefficient @xmath32
have been calculated first in  @xcite and the full result is now also known  @xcite .    to summarize : the matching coefficients required for the nnnlo calculation of the heavy - quark production cross section near threshold are known except for the non - fermionic contributions to @xmath25 , and the @xmath24 terms of @xmath31
. out of these , only @xmath25 is expected to have a major impact on the final result  @xcite .
work on these missing contributions is in progress .      after matching the heavy - quark currents and lagrangians the cross section calculation
is mapped to the calculation of the imaginary part of the two - point function @xmath21 of non - relativistic currents in ( [ r1 ] ) .
the leading - order colour coulomb potential is now part of the leading - order effective lagrangian , since the coulomb interaction is strong near threshold .
the propagators to be used in the perturbative calculation of the two - point function are the coulomb green functions , making this part of the computation similar to qed bound - state problems .
the calculation can be divided into three parts : contributions up to the third order involving only the coulomb potential , already completed in  @xcite ; the ultrasoft contribution , appearing first at nnnlo , which has been computed in  @xcite ; and finally , contributions involving at least one potential other than the coulomb potential ( `` non - coulomb potential contribution '' ) , which are not yet known .
the main result of this paper is the missing non - coulomb potential contribution .
compared to the coulomb contributions the major complication is the singular nature of the potential insertions .
the ultraviolet divergences must be regulated in dimensional regularization in a scheme consistent with the calculation of the matching coefficients order by order in the strong coupling , while retaining the resummation of infinitely many coulomb gluon exchanges by the use of coulomb green functions , whose @xmath33-dimensional expression is unknown .
the techniques we apply are an extension of those used in the nnlo calculation  @xcite .
since the method of that calculation was never written up ( though some results are scattered in @xcite ) , we devote some effort to presenting the third - order calculation in some technical detail .
in addition to the dominant production of the top - quark pair though a virtual photon there is also a @xmath34-boson contribution .
the contribution from the vector - coupling of the @xmath34 is trivially inferred from the photon - mediated cross section , while the axial - vector coupling contribution is suppressed near threshold and begins only at nnlo .
thus , only the first - order correction to the axial - vector in non - relativistic perturbation theory is needed .
some results at this order are available  @xcite , but none of these results are given in dimensional regularization .
we recalculated these small contributions , and will present them elsewhere , since in this paper we focus on the missing third - order terms in the vector current contribution .    to summarize : together with the results of this paper , the matrix element calculation is complete to nnnlo .
much less work has been done on electroweak and electromagnetic corrections .
counting the electromagnetic and electroweak coupling as two powers of the strong coupling , electromagnetic corrections contribute from nlo through the electromagnetic coulomb potential . this effect is easily included and has been discussed in  @xcite .
electroweak contributions to the matching coefficients of nrqcd currents and production operators have been calculated in  @xcite .
the formalism for calculating initial - state radiation , and soft and collinear photon corrections in general , simultaneous with summing coulomb exchange has been worked out in  @xcite for @xmath35-boson pair production .
its application to top quarks is straightforward .
it remains to be done to put all these electroweak effects together .    for top quarks
the sizeable decay width @xmath36 introduces further complications . in leading order the width is correctly included by evaluating the current two - point functions in pnrqcd at complex energy argument @xmath37  @xcite , where @xmath38 .
we adopt this prescription as the _ definition _ of the pure qcd contributions to the cross section .
beginning at nlo there exist further contributions related to the finite width .
since the physical final state is @xmath39 , boson may be regarded as stable .
] there exist irreducible `` backgrounds '' related to off - shell and single or non - resonant top - quark pair production .
in fact , the qcd contribution as defined above can not be unambiguously separated from electroweak effects at this order  perhaps not surprisingly , since the top quark width itself is such an effect  and the fact that a physical scattering cross section should refer only to stable ( or sufficiently long - lived ) particles in the final state , must be taken into account .
the incompleteness of the qcd cross section is signaled explicitly by the presence of uncanceled singularities in dimensional regularization with coefficients proportional to @xmath36 starting at nnlo .
the origin and consistent cancellation of these singularities is discussed in  @xcite and the corresponding calculations of electroweak effects and non - resonant contributions to a physical final state such as @xmath40 can in principle be performed in the framework of unstable - particle effective field theory  @xcite as already done for @xmath35 pair production  @xcite .
the corresponding calculation of the non - resonant nlo correction for top production has been performed in  @xcite and confirmed by a different method in @xcite . rather than embarking on the rather difficult computation of non - resonant contributions up to nnnlo
, however , a more promising approach suppresses them by appropriate invariant mass cuts  @xcite , which will anyway be applied experimentally .
calculations of top quark pair production near threshold with cuts on the final - state @xmath41 invariant masses have appeared recently  @xcite in the non - relativistic qcd and unstable particle effective theory frameworks .
as should be expected , the non - resonant contributions to the @xmath40 are sizeable below the nominal top pair production threshold , and hence can change the shape of the threshold cross section in the region of interest for the top - quark mass determination . in the present paper ,
we focus only on the qcd part of the problem .
we shall , however , make the finite - width @xmath42 poles explicit , so that they can be canceled analytically with future computations of electroweak effects .
indeed , the pole parts of the nnlo non - resonant contribution have recently been computed  @xcite and the cancelation has been verified .    to summarize : while the formalism for computing electroweak and finite - width effects consistently is in place , and bits and pieces have already been calculated , significant work remains to be done on electroweak effects at nnnlo . with the third - order qcd calculation now being nearly complete , improvements in the prediction of realistic cross sections ( including electroweak effects and cuts ) must be made .
0.3 cm since our intention to present the concepts , techniques and calculations in some detail resulted in a rather lengthy text , we have split it into two parts .
part i presents the effective field theory set - up , the nrqcd and pnrqcd matching coefficients and ends with a master formula for the third - order heavy - quark pair production cross sections .
this part could also be read as a review of non - relativistic effective theory in the weak - coupling regime complementary to  @xcite .
part ii contains the actual pnrqcd matrix element calculation together with a numerical estimate of the new contributions to the top cross section .
( a preliminary version has been presented in  @xcite . ) the small p - wave contributions to the cross section from the axial - vector coupling of the @xmath34-boson will be presented separately in  @xcite .
the outline of the paper i is as follows : in section  [ sec : eft ] we review the effective field theory framework and discuss the power - counting arguments that lead to the identification of the matching coefficients needed for the nnnlo calculation .
the subsequent two sections deal with matching qcd to a sequence of two non - relativistic effective theories , nrqcd and pnrqcd .
section  [ sec : nrqcd ] discusses the nrqcd aspects of the calculation .
in particular , we calculate the relevant one - loop matching coefficients including the new @xmath24 terms and collect all other results that feed into the cross section calculation .
the second matching step from nrqcd to pnrqcd is discussed at length in section  [ sec : pnrqcd ] , since a coherent summary is not yet available in the literature . among the new results of this section
are the path - integral derivation of the pnrqcd lagrangian ( neglecting ultrasoft gluons ) , the @xmath24 terms of the one - loop potentials , and a discussion on the non - renormalization of currents in the nrqcd to pnrqcd matching .
pnrqcd perturbation theory in the pole - mass scheme provides a poor approximation to the top - quark pair production cross section near threshold .
the top quark decay width is obviously an important effect , as is the conversion from the pole mass renormalization scheme that is employed in the primary calculations to renormalization schemes that absorb large corrections into the mass counterterm , which is a prerequisite for reliable perturbative calculations  @xcite .
furthermore , a resummation of pnrqcd perturbation theory for the green function is necessary in the vicinity of the bound state poles despite the sizeable top quark width .
these refinements will be explained in paper ii . in section  [ sec : master ] we conclude paper i by providing the master formula for the computation of the third - order cross section in pnrqcd .
in this section we present the relation of the pair production cross section to correlation functions of heavy quark currents together with the arguments , why this relation holds true at nnnlo .
we review the scales and momentum regions relevant to the problem both of which are central to the systematics of the effective theory approach .
the basic top pair production mechanisms in @xmath0 annihilation are shown in the upper part of figure  [ fig : eett ] . since we work to lowest order in the electromagnetic and electroweak couplings , the optical theorem allows us to relate the total cross section @xmath43 of the process @xmath44 to the two point functions of the vector and axial - vector heavy quark current .
we define @xmath45 for the vector current @xmath46 and the axial vector current @xmath47 .
the cross section is then given by @xmath48,\end{aligned}\ ] ] where @xmath49 is the high - energy limit of the @xmath50 production cross section , @xmath51 the center - of - mass energy squared , and @xmath52 the @xmath34-boson mass .
@xmath53 denotes the top quark electric charge in units of positron charge and @xmath54 is the electromagnetic coupling . the vector and axial - vector couplings of fermion @xmath55 to the @xmath34-boson
are given by @xmath56 with @xmath57 the weak mixing angle , @xmath58 the electric charge of fermion @xmath55 and @xmath59 its third component of the weak isospin .
the dominant production mechanism is through the coupling to the virtual photon .
the vector coupling of the @xmath34-boson increases the photon - mediated cross section by only about 8% in the threshold region @xmath60 .
the axial - vector contribution is even smaller , since the axial coupling is suppressed near threshold by the small velocity of the top quarks .
@xmath61 contributes to ( [ eq : sigtt ] ) only at nnlo relative @xmath62 .
( [ eq : sigtt ] ) which relates the inclusive top cross section to the spectral functions of heavy - quark currents is not exact .
there exist top production mechanisms , shown in the lower part of figure [ fig : eett ] , which are not captured by the heavy - quark current correlation functions , since the photon or @xmath34-boson couples to light quarks .
vice versa , there exist cuts contributing to @xmath63 related to annihilation subdiagrams , see figure  [ fig : annihilation ] below , which do not contain top quarks and hence should be excluded .
we shall discuss in the next subsection that these contributions are either highly suppressed and not relevant at third order , or can easily be included .
the characteristic non - relativistic energy in threshold production is much smaller than @xmath64 .
we define @xmath19 and the top quark `` velocity '' @xmath65 .
this is related to another often used velocity parameter @xmath66 by @xmath67 real values of @xmath68 and @xmath69 should be interpreted with a @xmath70 prescription .
we will extensively use another variable @xmath71 defined by @xmath72 here and below @xmath13 without any argument denotes the strong coupling in the @xmath73 scheme @xcite at the renormalization scale @xmath74 , and @xmath75 .
the value of @xmath71 determines when resummation to all orders in @xmath13 is necessary . above threshold ,
the variable @xmath71 is purely positive - imaginary , below threshold it is real and positive .
the threshold region is characterized by an absolute value of @xmath71 of order 1 or larger . in particular , the coulomb bound state poles are found at positive integer values of @xmath71 .
conventional fixed - order perturbation theory can be used only when @xmath76 . as will be discussed later the width of the top quark
is accounted for by substituting @xmath77 .
thus , as @xmath69 varies from @xmath78 to @xmath79 the variable @xmath71 sweeps through a curve in the complex plane that begins at the origin , moves out into the first quadrant into the direction of the positive real axis and returns to the origin from above near the imaginary axis .
the absolute value of @xmath71 along this curve is always smaller than @xmath80 , which is about 1.2 for top quarks , with @xmath81 , but since it reaches order one in the threshold region , the perturbation expansion in @xmath13 breaks down and resummation is necessary .      near the heavy - quark pair production threshold only a small kinetic energy @xmath82 is available to the final state . in the natural frame where @xmath83
this implies that the typical three momentum of a heavy quark is of order @xmath84 ( about 20  gev for top quarks ) , while the energy and momentum of any other nearly massless particle can at most be @xmath85 ( about 2  gev for tops ) .
the presence of several small scales propagates into the loop diagrams that contribute to the spectral functions and causes a breakdown of the standard perturbation expansion in the strong coupling @xmath13 .
however , since @xmath12 is small one does not have to compute the loop integrals exactly  an expansion in @xmath12 suffices .
this leads to a reorganized expansion as shown in ( [ systematics ] ) , in which @xmath13 and @xmath12 are expansion parameters but @xmath86 or , equivalently , @xmath71 is of order one .    for a given feynman diagram
the expansion in @xmath12 can be constructed without first computing the full expression using the threshold expansion @xcite .
the method uses that every diagram is the sum of terms , for which each loop momentum is in one of the following four regions : @xmath87 when on - shell , only massless particles ( gluons , light quarks and ghosts ) can be ultrasoft , and only the heavy quarks can be potential . in each region , the loop integrand is expanded in the terms which are small in the corresponding region and the loop integration of the expanded integrand is carried out over the complete @xmath33-dimensional space - time volume .
the expansion generates ultraviolet and infrared divergences which are regulated dimensionally ( @xmath88 ) and subtracted according to the @xmath73 prescription .
however , the divergences generated by the separation of the diagram into regions cancel in the sum over all terms .
the procedure just described is largely equivalent to constructing appropriate effective lagrangians within dimensional regularization , but it clarifies the correct matching procedure which is subtle in dimensional regularization if the effective theory contains more than one scale as is the case in non - relativistic qcd @xcite .
the threshold expansion also synthesizes the non - relativistic velocity power counting rules developed for the different modes ( momentum regions ) in @xcite . in the construction of the non - relativistic and resummed expansion of the pair production cross section we derive effective lagrangians in two steps by integrating out the large momentum modes according to the following scheme :    @xmath89 \quad & \mu >
m \\[0.3 cm ] \big\downarrow & \\[0.3 cm ] { \cal l}_{\rm nrqcd}\,[q(s , p),\,g(s , p , us ) ] \quad &   m v < \mu < m \\[0.3 cm ] \big\downarrow & \\[0.3 cm ] { \cal l}_{\rm pnrqcd}\,[q(p),\,g(us ) ] \quad & \mu
< m v\\[0.3 cm ] \end{array}\ ] ]    in square brackets we display the modes of the heavy quarks ( @xmath90 ) and massless particles ( @xmath91 ) which are still contained in the effective lagrangian ; the others are integrated out when the energy cut - off @xmath74 is lowered as indicated on the right .
the first step leads to nrqcd  @xcite , in which all interactions are local , since only the short - distance hard modes have been eliminated .
the expansion rules of the threshold expansion _ define _ the dimensionally regularized nrqcd lagrangian .
the second step whereby soft modes and potential massless modes are integrated out was suggested in @xcite in the context of the effective lagrangian and threshold expansion method .
the result is the potential nrqcd ( pnrqcd ) lagrangian  @xcite .
the pnrqcd lagrangian is not local .
it contains spatially non - local but temporally local , i.e. instantaneous interactions of the heavy quarks , since the three - momentum of the potential heavy quark field still present in pnrqcd is of the same order as the one of the modes integrated out .
these interactions provide a precise definition of the concept of `` heavy - quark potentials '' .
perturbation theory in pnrqcd resembles quantum - mechanical perturbation theory closely , since the leading colour - coulomb interaction is part of the unperturbed theory .
thus , the propagator of pnrqcd includes the leading coulomb interaction exactly , which effects the required resummation of conventional perturbation theory to all orders .    to illustrate the velocity scaling of feynman diagrams , we consider the power counting for the loop integrand .
( [ eq : scales ] ) implies that the integration measure @xmath92 scales as @xmath93 , @xmath94 , @xmath95 and @xmath96 , when @xmath97 is hard , soft , potential and ultrasoft , respectively .
the denominator of a gluon ( massless ) propagator with momentum @xmath98 is approximated at leading order in a given region by : @xmath99 \displaystyle -\frac{1}{\bm{\ell}^2}+\ldots & \mbox {   potential } ( v^{-2 } ) \end{array } \right.\ ] ] the velocity scaling is given in brackets .
for the heavy - quark propagator with momentum @xmath100 the denominators are : @xmath101 \displaystyle \frac{1}{2 m}\,\frac{1}{\ell^0 } + \ldots & \mbox {   soft } ( v^{-1 } ) \\[0.5 cm ] \displaystyle
\frac{1}{2 m}\,\frac{1}{e/2+\ell^0-\bm{\ell}^2/(2 m ) } + \ldots & \mbox {   potential } ( v^{-2 } ) \end{array } \right.\ ] ]     ladder diagrams.,scaledwidth=90.0%,height=68 ]    0.2 cm using these scaling rules , it is easy to see why an all - order resummation of feynman diagrams is required in the threshold region .
it will become clear from the later systematic derivation that the relevant diagrams are the ladder diagrams shown in figure  [ fig : ladder ] , and that the dominant term in the velocity expansion arises from the loop momentum region when all loop momenta are in the potential region .
adding an additional rung to the ladder adds one potential gluon ( @xmath102 ) and two potential heavy - quark propagators ( @xmath103 ) to the diagram .
the numerator of the diagram contains no velocity suppression factors , hence accounting for the potential loop measure ( @xmath95 ) and strong coupling from the two additional vertices ( @xmath104 ) , we find that each rung provides a factor of order @xmath86 , which is unsuppressed in the threshold region .
it will be seen below that only potential gluon exchange generates this @xmath105 enhancement , which is equivalent to the statement that only the leading coulomb interaction must be included in the unperturbed effective lagrangian .
we now return to the discussion of heavy - quark production mechanisms not captured by ( [ eq : sigtt ] ) , which expresses the cross section in terms of the heavy - quark current spectral functions . in the case of heavy - quark radiation ( lower left in figure  [ fig : eett ] ) the final state consists of @xmath106 , and since the available energy at threshold is limited to @xmath85 , the light quarks must be ultrasoft . in the three - loop diagram that represents the square of the heavy - quark radiation amplitude ,
the @xmath107 loop must be potential and the two other loops ultrasoft , which leads to a factor of @xmath108 from the loop integration measure .
the intermediate gluon and light - quark propagators in the amplitude must be hard to produce the @xmath107 pair and hence do not contribute inverse powers of @xmath12 .
the two potential heavy quark and the two ultrasoft light quark propagators ( @xmath109 ) supply a factor of @xmath102 each , so the heavy - quark radiation contribution to the cross section scales at least as @xmath110 which should be compared to @xmath12 for the leading term .
inspection of the analytic expression @xcite confirms this result , hence this contribution can be safely neglected . in the case of singlet production
( lower right in figure  [ fig : eett ] ) through the coupling of the virtual photon or @xmath34-boson to light quarks the dominant term comes from three hard loops , leading to the counting @xmath111 , which represents a third - order correction to the cross section . while not part of @xmath112 this mechanism can be included in the three - loop short - distance coefficient @xmath113 of the non - relativistic heavy - quark current , which is discussed below , although it is not known at present .
note that this contribution to @xmath113 is complex , but the imaginary part should be discarded , since it corresponds to the three - gluon and light - quark cut , which is not part of the heavy - quark production cross section .
a similar singlet - production diagram exists for the axial - vector coupling with only two gluons coupling to the light - quark triangle , but due to velocity suppression this contribution begins only at fourth order .    we have thus argued that the production mechanisms not included in @xmath114 are either suppressed or easily included at third order .
consider now figure  [ fig : annihilation ] , which shows a diagram contained in @xmath112 , but whose three - gluon cut should not be part of the heavy - quark cross section .
the possible loop momentum regions for this diagram are h - h - h - h , p - h - h - p , p - h - h - h and h - h - h - p , where the first and last letter refers to the left and right heavy quark loop , respectively . in the all - hard configuration
only the three - gluon cut contributes to @xmath112 , so the correct prescription is to simply not include this configuration .
the p - h - h - p configuration may be interpreted as heavy - quark production followed by rescattering through annihilation .
annihilation is suppressed by @xmath115 relative to coulomb exchange as can be seen by counting loop integration and propagator factors , so this configuration is relevant only from fourth order as a contribution to the @xmath116 potential .
the remaining two configurations with only one potential heavy quark loop are analogous to the singlet production mechanisms .
that is , one drops the imaginary part of the h - h - h subdiagram , which comes from the three - gluon cut and associates its real part to @xmath117 .
there exists a three - loop diagram similar to figure  [ fig : annihilation ] but with two gluon lines only for the axial - vector coupling , but as in the singlet diagram with a light - quark triangle discussed above , the axial - vector coupling implies another factor of @xmath118 , so this diagram is never relevant at third order .
in this section we discuss the matching of the vector current correlation function for @xmath60 to its equivalent representation in non - relativistic qcd . this amounts to integrating out the hard modes , which correspond to `` relativistic effects '' involving the scale of the heavy - quark mass .
non - relativistic qcd is expressed in terms of a two - component quark field @xmath119 and the corresponding anti - quark field @xmath120 to represent the remaining soft and potential fluctuations of the original quark field .
the effective gluon field @xmath121 can be soft , potential and ultrasoft .    before going into the details of the lagrangian and power counting we briefly sketch the result . as will be shown below the expansion of the vector current @xmath122 in terms of the non - relativistic fields is given by @xmath123 where the hard matching coefficients @xmath17 , @xmath18 have perturbative expansions in @xmath13 .
in the `` rest frame '' @xmath124 , eq .
( [ eq : pi ] ) implies @xmath125 , so @xmath126 g(e ) + \ldots,\ ] ] where the neglected terms on the right - hand side include a subtraction term that does not contribute to the imaginary part of @xmath62 as well as terms beyond the third order ( nnnlo ) .
the important quantity is the two - point function of the non - relativistic current @xmath127(x)\ , [ \psi^{\dag}\sigma^i\chi](0 ) )    evaluated in non - relativistic qcd ( nrqcd ) .
the terms proportional to @xmath69 in ( [ eq : pitonrqcd ] ) arise from expanding the prefactor @xmath128 and from the @xmath129 suppressed current in ( [ eq : qcdvectorcurrent ] ) , whose matrix element can be reduced to the one of the leading current by an equation - of - motion relation derived later .
thus the main ingredients to the non - relativistic representation are the calculation of @xmath21 and the current matching coefficients .
similar relations hold for the axial - vector contribution to the cross section ( [ eq : sigtt ] ) , which arises from @xmath34-boson exchange .
the axial - vector current @xmath130 is represented in nrqcd by the expansion
@xmath131\chi + \ldots,\end{aligned}\ ] ] with hard matching coefficient @xmath132 .
as is the case for the vector current , only the spatial components of the current contribute to the cross section , since the lepton tensor from the @xmath0 initial state is transverse to both initial state momenta when the electron mass is neglected . only the leading term in the @xmath133 expansion
is needed for nnnlo accuracy , since the derivative in the leading current implies the well - known p - wave velocity suppression . the qcd correlation function
is then given by the expression @xmath134^\dagger(x)\ , [ \psi^{\dag}\gamma^i\chi](0 ) )    + \ldots,\qquad \nonumber\end{aligned}\ ] ] where @xmath135 $ ] .
for the present purpose the non - relativistic effective lagrangian can be divided into five parts , @xmath136     nrqcd feynman rules for two - quark vertices up to order @xmath129 . dashed ( curly )
lines denote the @xmath137 ( @xmath138 ) gluon field .
@xmath139.,title="fig:",height=748 ] 0.3 cm    -0.3cmthe gluon field is contained in the gauge - covariant derivative @xmath140 , field strength tensor @xmath141 $ ] , and the chromoelectric and chromomagnetic fields defined as @xmath142 , \nonumber \\ \bm{\sigma}\cdot{{\bf { b } } } & \equiv & -\frac{1}{2}\,\sigma^{ij } g^{ij } , \label{eq : magneticfield}\end{aligned}\ ] ] where @xmath143 $ ] and @xmath144 . with these definitions
the bilinear heavy - quark lagrangian is given by .
the @xmath145 term is misprinted in eq .
( 8) of @xcite , where @xmath146\psi$ ] should read @xmath147 . ]
@xmath148   & & + \ , \psi^{\dag}\bigg(\frac{d_2 g_s}{8 m^2}[d^i , e^i ] +   i \ , \frac{d_3 g_s}{8 m^2}\sigma^{ij}\{d^i , e^j\}\bigg ) \psi \nonumber \\[0.1 cm ]   \nonumber & & + \,\psi^{\dag}\bigg(-d_w g_s\frac{\{{\bf d}^2,{\bm { \sigma } }   \cdot { \bf b}\}}{8m^3 } + d_a g_s^2 \frac{{\bf b}^2-{\bf e}^2}{8m^3}+ d_b g_s^2\frac{\sigma^{ij}(b^i b^j- e^i e^j)}{8m^3}\,\bigg ) \psi \nonumber\\[0.1 cm ] & & + \,o\bigg(\frac{1}{m^4}\bigg ) , \\[0.2 cm ]
{ \cal l}_{\chi}&= & - { \cal l}_{\psi}\quad\mbox{with   $ \psi\to\chi , i d^0\to -i d^0 ,   e^i\to - e^i$}\,.\end{aligned}\ ] ] note that in @xmath33 dimensions we can not define the @xmath149 individually , since they do not represent the components of a @xmath150 dimensional vector
. however , all we need are scalars such as ( [ eq : magneticfield ] ) and @xmath151 , @xmath152 , which can be consistently defined through the @xmath33 dimensional field strength tensor : @xmath153\,.\ ] ]    the four - fermion quark - antiquark terms in the effective lagrangian read @xmath154 \psi \ , \chi^{\dag}[\sigma^i,\sigma^j ]   \chi + \nonumber \\ \label{eq : nrqcdqqqq } & & + \,\frac{d_{vs}}{m^2 } \psi^{\dag } t^a \psi \
, \chi^{\dag } t^a \chi - \frac{d_{vv}}{8 m^2 } \psi^{\dag } t^a [ \sigma^i,\sigma^j ] \psi \ , \chi^{\dag } t^a [ \sigma^i,\sigma^j ] \chi+o\bigg(\frac{1}{m^3}\bigg),\end{aligned}\ ] ] where the factor @xmath155 in the definition of the spin - triplet operators has been inserted , since in four dimensions @xmath156\otimes [ \sigma^i,\sigma^j ] =   -8 \,\sigma^i\otimes\sigma^i$ ] .
the pure gluon lagrangian takes the form @xmath157 finally , @xmath158 is the same as the light - quark lagrangian in full qcd .
the feynman rules for the @xmath133 and @xmath129 terms in @xmath159 are given in figure  [ fig : nrqcd_fr ] .
in the fifth row can be eliminated using the heavy - quark equation of motion .
this generates @xmath160 terms not shown in the figure and modifies the @xmath129 four - point vertices in the last row of the figure as follows : drop the @xmath161 and @xmath162 terms and multiply the remaining terms by two .
] the hard matching coefficients @xmath163 equal one at tree level , while the others vanish .
all coefficients obtain one - loop corrections , which must be determined by matching the qcd diagrams to the nrqcd diagrams order by order in the non - relativistic expansion . for reasons which will be explained later , the matching coefficients
are needed to order @xmath164 in the dimensional regularization parameter .
in writing ( [ eq : nrqcdqq ] ) we did not include @xmath160 interactions with vanishing tree - level coefficients as well as mixed heavy - light quark operators of the form @xmath165 in @xmath166 , since they are irrelevant for the nnnlo calculation , as we discuss now .
the velocity scaling of the fields depends on the momentum region .
heavy quark fields can be potential or soft . from the scaling rules of the propagator ( compare ( [ eq : quarkprop ] ) ) and the integration measure
, it follows that in both cases the field scales as @xmath167 .
the same power counting can be done with the gluon field , which is either soft , potential or ultrasoft and one finds that @xmath168 scales as @xmath167 , @xmath169 and @xmath170 respectively , which includes a factor of @xmath171 from the coupling constant @xmath172 .
turning to the effective lagrangian , we first consider the bilinear terms in the heavy lagrangian @xmath166 . for potential quarks , the bilinear terms in the kinetic term @xmath173
are both of order @xmath95 .
the relativistic correction @xmath174 scales as @xmath175 .
being suppressed by @xmath118 relative to the leading terms it contributes from nnlo .
the next term in the expansion of the relativistic energy - momentum relation would be of order @xmath176 and is beyond nnnlo .
for soft quarks , the quadratic kinetic energy term is an order @xmath12 correction to the leading - order static lagrangian @xmath177 , which scales as @xmath94 .
this explains why static heavy - quark propagators can be used in the calculation of the heavy - quark potentials . for soft heavy quarks
the quartic kinetic energy correction is a nnnlo effect .
consider now the interactions of the heavy quark with the gauge field , i.e. terms of the form @xmath178 , potentially with derivatives on the quark and gluon fields .
the @xmath179 interaction that arises from @xmath180 scales as @xmath95 when all fields are potential .
therefore it is not suppressed relative to the bilinear terms that persist as @xmath181 . in the potential region
, this interaction has to be treated non - perturbatively ; this is why the cross section near threshold requires a summation of some loop momentum contributions to all orders in @xmath13 . when the gluon field is soft , the velocity scaling of the @xmath179 interaction is @xmath182 , but since in this case the leading term @xmath183 scales as @xmath94 , this interaction is now a perturbation .
it follows that all soft interactions can be treated in conventional perturbation theory .
further three - point interactions of the form above carry derivatives .
each derivative gives at least a suppression of @xmath12 so interactions with up to three derivatives in @xmath184 may contribute at nnnlo .
the requirement of gauge invariance restricts the possible interaction terms to the so - called chromomagnetic interaction @xmath185 at order @xmath186 ( nlo ) , and the darwin and spin - orbit interactions at order @xmath175 ( nnlo ) , multiplied by the short - distance coefficients @xmath187 , @xmath188 and @xmath145 , respectively
. however , a single chromomagnetic interaction contributes only in connection with a @xmath12 suppressed quark - gluon vertex from the @xmath189 interaction , and not at all to the current correlation function , since the trace over an odd number of pauli matrices vanishes .
thus the chromomagnetic , darwin and spin - orbit interaction all start to contribute at nnlo ( with two insertions of the chromomagnetic term ) .
hence , at nnnlo one needs the coefficient functions @xmath187 , @xmath188 and @xmath145 in the one - loop approximation . beyond order
@xmath175 only the interactions with non - vanishing tree - level coefficient functions given in ( [ eq : nrqcdqq ] ) can potentially contribute to nnnlo .
none of them does , however , since single insertions of interactions with pauli matrices vanish , as discussed above , while the terms with two electric or magnetic field strengths are of the form @xmath190 with two derivatives , which is smaller than nnnlo in both , the potential and soft region .
bilinear heavy quark operators in conjunction with light quarks , @xmath165 , are of order @xmath186 . for
these operators to contribute to the heavy - quark current correlation function at least one interaction of the form @xmath191 is required , which costs a factor of @xmath12 or , equivalently @xmath13 .
thus , a @xmath165 operator is relevant at nnnlo , if its short - distance coefficient function is of order @xmath13 .
the only operator that may have a tree - level coefficient is @xmath192 , since in this case the intermediate potential gluon propagator can be cancelled , making the operator local .
this operator is generated at tree level from the darwin term in the lagrangian ( [ eq : nrqcdqq ] ) by the use of the field equation for the chromo - electric field .
our convention is that we do not eliminate the darwin term by the field equation , hence @xmath165 operators must be added to the lagrangian only with coefficients of order @xmath193 producing corrections to the heavy - quark current correlation function beyond nnnlo .    we therefore conclude that only the terms in the first two lines of ( [ eq : nrqcdqq ] ) are needed for the nnnlo calculation , the same as in nnlo .
the only difference to nnlo is that the short - distance coefficients @xmath187 , @xmath188 and @xmath145 are required at the one - loop order .
we also see that only the @xmath194 interaction in the potential region is non - perturbative , and this explains why only ladder diagrams of coulomb gluons must be summed to all orders .
the four - quark operator lagrangian @xmath195 ( [ eq : nrqcdqqqq ] ) is generated by hard scattering of quarks and anti - quarks , or by quark anti - quark annihilation .
hard scattering with momentum exchange of order @xmath196 requires the exchange of at least two gluons , corresponding to one - loop diagrams so the coefficient functions are of order @xmath193 .
the four - quark operator counts as @xmath186 . including the one - loop coefficient function gives the counting @xmath197 , which is a nnnlo effect relative to the leading - order lagrangian of order @xmath95 for the scattering of potential quarks .
the annihilation contribution is present already at order @xmath13 ( tree - level ) , but the operator has the colour structure @xmath198 , which does not contribute to the current correlation function due to @xmath199 .
similarly , annihilation into two gluons does not contribute to the vector correlation function , since it leads to fermion loops with three vector couplings that vanish by charge conjugation .
thus , at nnnlo , we can restrict ourselves to the four - fermion operators generated by hard quark anti - quark scattering .
for this reason we write the operators in the `` scattering ordering '' @xmath200 rather than the `` annihilation ordering '' @xmath201 .
although the two orderings are related by a fierz transformation in four space - time dimensions , the two are inequivalent in dimensional regularization . in general
, we would have to introduce the difference of the two as evanescent operators .
this complication is avoided here , since there are no annihilation contributions at nnnlo . adopting the `` scattering ordering '' in the lagrangian
, we do not need to perform any fierz transformations .
the necessity to avoid relations that hold only in four dimensions is also the reason for introducing the definition ( [ eq : magneticfield ] ) that does not make use of the three - dimensional @xmath202 symbol , which is not defined in @xmath150 dimensions . in particular , the commutator @xmath203\ , , \ ] ] must be considered as an independent element of the @xmath33-dimensional algebra of pauli matrices .
this poses no difficulty in the calculation of the vector current correlation function , since in the end all expressions can be evaluated using the @xmath33-dimensional identities @xmath204 where @xmath88 and @xmath205 .
the pure gauge field lagrangian ( [ eq : nrqcdg ] ) follows from integrating out heavy - quark loops with small momentum gluon lines attached .
the renormalization of the standard kinetic term by the coefficient @xmath206 is well - known to be related to the matching of the strong coupling from the @xmath207 flavour theory including the heavy quark to the theory with @xmath208 light flavours . in the following
we express all results in terms of the strong coupling in the @xmath73 scheme in the @xmath208 flavour theory , which is the appropriate coupling for calculations in nrqcd , where the heavy quark short - distance fluctuations have been integrated out . after redefining the strong coupling , @xmath206
should be set to one in ( [ eq : nrqcdg ] ) , so that the kinetic term is canonically normalized .
the next term , @xmath209 , in the gauge field lagrangian involves two derivatives and a coefficient function @xmath210 . therefore it must be included at nnnlo . on the other hand , the term involving three gluon field strengths can be neglected at this order .
having collected the relevant terms in the effective lagrangian , we are now in the position to discuss the matching calculations .
most of the results required at nnnlo are available in the literature .
however , many of the matching coefficients multiply nrqcd correlation functions , which exhibit @xmath42 poles . thus , as will be explained in section  [ sec : matchcurrent ] below , we also need the @xmath24 terms of the matching coefficients , which have not been calculated or presented up to now .
we therefore had to repeat these matching calculations and extend them to the @xmath24 terms .
the matching calculation is performed in the center - of - mass frame , so that the three - momenta of the heavy quark and anti - quark are of opposite sign .
the external heavy quark spinors in qcd are given by @xmath211 \bm{\sigma}\cdot{{\bf { p}}}\,\xi \end{array}\right ) , \quad v(p ) = \frac{1}{(e_p+m)^{1/2}}\left(\begin{array}{c } \bm{\sigma}\cdot{{\bf { p}}}\,\eta \\[0.2 cm ]   ( e_p+m)\ , \eta \end{array}\right),\ ] ] for external momentum @xmath212 with @xmath213 .
the variables @xmath214 and @xmath215 denote the quark and anti - quark two - spinors , respectively .
they are normalized according to @xmath216 .
the coefficient functions of the interactions of heavy quarks with a single gauge field in @xmath166 can be deduced from the heavy quark form factors in background field gauge .
the finite part of the one - loop form factors was calculated in @xcite .
however , the order @xmath164 coefficients were not computed there .
the on - shell form factors can be brought into the general form @xmath217u(p),\end{aligned}\ ] ] where @xmath218 .
for the one - loop diagrams shown in figure [ fig : formfactors ] we obtain the following result for the expansion of the ultraviolet - renormalized form factors in @xmath13 and @xmath219 up to order @xmath220 : @xmath221 , \\
\nonumber f_2 & = & \frac{\alpha_s}{\pi}\bigg(\frac{\mu}{m}\bigg)^{2\epsilon } \frac{\gamma(\epsilon)e^{\gamma_e\epsilon}}{24(4\epsilon^2 - 1 ) } \bigg[c_a\bigg(6(2\epsilon+1)(2\epsilon^2 - 1)+\frac { q^2}{m^2}(4\epsilon^4 + 8\epsilon^3 + 5\epsilon^2 - 2\epsilon-6)\bigg ) \\   & &   + \,c_f\bigg(-12\epsilon(2\epsilon+1)^2-\frac { q^2}{m^2}2\epsilon(\epsilon+1)(2\epsilon+1)^2\bigg)\bigg].\end{aligned}\ ] ] here and below we use the standard colour factors @xmath222 , @xmath223 , @xmath224 .
the remaining divergences are infrared divergences of the on - shell form factors . by calculation of the corresponding form factors in the effective theory
we obtain the relations between the coefficients @xmath225 and the form factors : @xmath226 where @xmath227 and @xmath228 . since @xmath229 exactly , this implies the well - known relation @xmath230 .
the @xmath73 renormalized coefficient functions follow by subtracting the @xmath42 poles from the above expressions .
the @xmath24 terms of the above expressions are in agreement with @xcite .
one - loop form factor diagrams : wavefunction renormalization and vertex corrections.,height=113 ]      to obtain the bilinear pure gauge field lagrangian @xmath231 , the gluon self energy has to be matched . in the one - loop order
the relevant diagram is the heavy - quark loop , which gives : @xmath232 the results agree for @xmath233 with the ones in @xcite .
as mentioned above , the operator with coefficient @xmath234 is not needed , because it can contribute to the heavy quark correlation function only with an additional loop , which is beyond nnnlo .
recall that @xmath235 should be used after having normalized the fields canonically .
qcd diagrams for the four - fermion operators at one loop.,height=68 ]    the remaining part of the lagrangian to be matched to qcd is the one containing the four - fermion operators . as discussed in section  [ sec : hqcf ] we do not need the annihilation contributions and therefore restrict ourselves to the scattering diagrams shown in figure  [ fig:4quark ] . the results in @xmath233 can be obtained from the equal mass limit of the unequal mass case given in @xcite . here
we present the @xmath33-dimensional matching coefficients : @xmath236 this agrees with the finite part of the unequal mass case in @xcite .
the @xmath24 terms of the above expressions are in agreement with @xcite .
and @xmath237 defined in @xcite are @xmath238 times those above . ]
we finally need an expression for the heavy quark vector current @xmath239 to nnnlo accuracy in the effective theory .
the perturbative matching coefficients of the nrqcd currents come from diagrams where the hard loop connects to one of the external current vertices .
since the zero component of the vector current is irrelevant , we focus on matching the operator @xmath240 .    at leading order in the velocity expansion
the unique nrqcd vector current is @xmath241 with coefficient function @xmath17 as given by the first term on the right - hand side of ( [ eq : qcdvectorcurrent ] ) .
the precise definition of the matching coefficient is @xcite @xmath242 where @xmath243 are the on - shell wave function renormalization constants in qcd and nrqcd , respectively .
@xmath244 represents the amputated , bare electromagnetic current vertex function evaluated for on - shell heavy quarks directly at threshold , i.e.  with zero relative momentum , expressed in terms of the renormalized qcd coupling and pole mass . in dimensional regularization ,
@xmath245 , and @xmath246 , the corresponding nrqcd vertex function , equals its tree - level expression @xmath247 , since the nrqcd integrals for zero external relative momentum are scaleless . here
it is important that the threshold expansion is employed to define nrqcd in dimensional regularization .
thus , @xmath17 equals the uv renormalized on - shell qcd vertex directly at threshold with infrared divergences subtracted recursively by the nrqcd renormalization factor @xmath248 .
this definition is equivalent at order @xmath249 to extracting the purely hard ( h - h-
-h ) momentum regions in the threshold expansion of the @xmath250-loop vertex function with external heavy quark momenta in the potential region .
the coefficient @xmath17 is needed to three - loop accuracy to achieve nnnlo precision .
while the two - loop expression has been known for some time @xcite , the three - loop correction is not yet completed and constitutes one of the missing pieces in the complete third - order calculation of the cross section .
the three - loop diagrams involving at least one fermion loop are , however , already known ) . ]
@xcite , as are the logarithmic terms related to the anomalous dimension of the current and strong coupling renormalization  @xcite .
defining @xmath251 the coefficients of perturbative expansions of any quantity @xmath252 in @xmath253 , through @xmath254 and the coefficients of the qcd @xmath255 function in the @xmath73 scheme with @xmath208 light flavours ( not including the heavy quark ) @xcite @xmath256 the known results for the vector current matching coefficients are given by : @xmath257 c_v^{(2)}(\mu ) & = & 2
\beta_0 l_m   c_v^{(1)}(m )   + l_m   c_f \pi^2 \bigg[- \frac{16}{3 } c_f-8 c_a \bigg ] + c_v^{(2)}(m ) , \\[0.2 cm ] c_v^{(3)}(\mu ) & = & \left(4\beta_0 ^ 2 l_m^2 + 2\beta_1 l_m\right)c_v^{(1)}(m ) \nonumber \\ & & + 4\beta_0 l_m   \bigg\ {   l_m c_f \pi^2 \bigg[- \frac{16}{3 } c_f-8 c_a \bigg ] + c_v^{(2)}(m ) \bigg\ } \nonumber \\ & & + l_m^2 c_f \pi^2 \bigg [ -40c_f^2 -\frac{188}{9 } c_f c_a + \frac{128}{3 } c_a^2 \nonumber \\ & & \hspace{2 cm } -\frac{128}{9 } c_f t_f n_f -\frac{64}{3 } c_a   t_f n_f \bigg ] \nonumber \\ & & + l_m c_f \pi^2 \bigg [   \left(-72 + 192\ln{2}\right ) c_f^2 + \left(-\frac{1888}{27}-96\ln{2}\right ) c_f c_a \nonumber   \\ & & \hspace{2 cm } + \left ( -\frac{256}{9}-96\ln{2}\right )   c_a^2 + \frac{800}{27 } c_f t_f n_f \nonumber   \\ & & \hspace{2 cm } + \frac{296}{9 } c_a t_f n_f -\frac{32}{5 } c_f t_f \bigg ] + c_v^{(3)}(m),\end{aligned}\ ] ] where @xmath258 is the matching coefficient evaluated at @xmath259 .
the two - loop non - logarithmic terms are fully known but the three - loop ones only for the fermionic contributions .
they read : @xmath260 , \\[0.2 cm ] c_{v}^{(3)}(m ) & = & 64\,\bigg [ c_f t_f n_f   \bigg [    c_f \left(46.692(1)\right )   + c_a \left(39.623(1)\right ) \nonumber \\ & & \hspace{1.5 cm } + t_f n_f \left(-\frac{163}{162}-\frac{4}{27}\pi^2\right ) + t_f \left(-\frac{557}{162}+\frac{26}{81}\pi^2\right )   \bigg ] \nonumber \\ & &   + c_f t_f   \bigg [    c_f \left(-0.840(2)\right )   + c_a \left(-0.09(2)\right ) \nonumber \\ & & \hspace{1.5 cm } + t_f \left(-\frac{427}{162}+\frac{158}{2835}\pi^2+\frac{16}{9}\zeta_3\right )   \bigg]\,\bigg ] + c_{v , g}^{(3)}(m).\end{aligned}\ ] ] where the last term @xmath261 is the unknown gluonic and singlet contribution and @xmath262 is a short - hand for the riemann zeta function value @xmath263 .
we recall our convention that @xmath13 denotes the strong coupling in the @xmath73 scheme with @xmath208 light flavours .    turning to the next orders in the velocity expansion , we find the operators @xmath264 which are suppressed by @xmath265 relative to the leading nrqcd current .
further operators of dimension five contain the ultrasoft gauge field strength @xmath266 of order @xmath182 .
thus up to nnnlo all production vertices contain only the quark - antiquark pair .
the on - shell heavy quark - antiquark production vertex in full qcd can be decomposed into the expression @xmath267 v(p_2),\ ] ] where now @xmath268 , @xmath269 , and @xmath270 . inserting ( [ eq : spinors ] ) for the external spinors ,
we obtain the exact expression @xmath271 \xi^\dagger\sigma^i\eta -2 e_p \left[\frac{m}{e_p } \hat f_1(q^2)-\hat f_2(q^2)\right ]   \frac{p^i}{m ( e_p+m ) }   \,\xi^\dagger{{\bf { \sigma}}}\cdot{{\bf { p}}}\,\eta.\quad\ ] ] this shows explicitly that only quark - antiquark operators with a sigma matrix can appear as assumed in ( [ eq : currentsv2 ] ) .
expanding this expression in @xmath272 we find @xmath273_{\,|q^2=4 m^2 , \,\mbox{\scriptsize hard } } , \nonumber\\[0.2 cm ]   d_{va } & = &   [ \hat f_1-\hat f_2]_{\,|q^2=4 m^2 , \,\mbox{\scriptsize hard } } , \nonumber\\   d_{vb } & = & ( -4 ) \,\frac{d}{d(q^2/m^2 ) }   [ \hat f_1+\hat f_2]_{\,|q^2=4 m^2 , \,\mbox{\scriptsize hard}}.\end{aligned}\ ] ] the subscript `` hard '' means that only the hard regions should be included in the computation .
the one - loop hard form factors can be extracted from @xcite , by dropping the non - analytic terms in the expansion in @xmath274 , which originate from the potential region .
we obtain @xmath275 given in ( [ eq : cv1 ] ) above and @xmath276 + o(\alpha_s^2)\ , .
\label{eq : dvab}\end{aligned}\ ] ] the logarithm in @xmath277 arises as the consequence of mixing with the leading order current through sub - leading nrqcd interactions , see ( [ eq : usdiv ] ) below .
these results agree with the computation of the one - loop corrected matching coefficients of the sub - leading current operators through explicit nrqcd matching @xcite .    at nnnlo
the correlation functions of velocity - suppressed currents will be evaluated only with the leading and next - to - leading order coulomb potential , which is spin - independent . hence , only the traces @xmath278 appear .
this allows us to combine @xmath279 such that the qcd vector current is now represented by @xmath280 as anticipated in ( [ eq : qcdvectorcurrent ] ) . from ( [ eq : dvab ] ) and ( [ eq : currentspinprojection ] ) we obtain @xmath281 + o(\alpha_s^2)\ , .
\label{eq : dv1loop}\end{aligned}\ ] ]    the explicit scale dependence of the matching coefficients is due to evolution of the strong coupling and the factorization of the hard scale .
it must cancel when all contributions to the cross section are combined .
we have checked explicitly that this is indeed the case .
note that we do not need the @xmath24 terms of the coefficient functions @xmath17 and @xmath18 to compute the heavy - quark current correlation function , since they multiply the finite , renormalized nrqcd correlation function in ( [ eq : pitonrqcd ] ) .
one may wonder then what is the difference between the nrqcd current and the nrqcd lagrangian matching coefficients @xmath187 , @xmath188 etc .
, since for the latter we need the @xmath24 terms as stated and given above .
the reason is the particular definition ( [ eq : newmatch ] ) of the current matching coefficient .
imagine that we calculate the qcd and nrqcd vertex functions @xmath244 with non - vanishing external relative momentum .
then the nrqcd diagrams are no longer scaleless and @xmath246 is the sum of potential , soft and ultrasoft loop momentum contributions . because of the @xmath105 factors from potential gluon exchange , the higher - dimensional nrqcd interactions contribute to the leading current matching equation at some order in perturbation theory . as a result @xmath246 will be different whether one uses the nrqcd lagrangian with four - dimensional or with @xmath33-dimensional short - distance coefficients .
the difference is , however , a local term that can be absorbed into the matching coefficients of the external currents .
which definition does ( [ eq : newmatch ] ) correspond to ?
suppose we follow the more conventional path to define the renormalized effective lagrangian with @xmath233 short - distance coefficients . in this case
@xmath246 does _ not _ represent the sum of all potential , soft and ultrasoft loop momentum regions plus the hard ones not connected to the external vertex ( encapsulated in the lagrangian matching coefficients @xmath187 etc . ) , since one misses some @xmath24 terms from hard subgraphs that multiply divergent soft , potential or ultrasoft loops .
these missing local contributions can be and must be added back by adapting the external current matching coefficient .
hence the matching coefficient @xmath17 defined by this prescription does not correspond to ( [ eq : newmatch ] ) . on the other hand
, the nrqcd lagrangian with @xmath33-dimensional short - distance coefficients reproduces these missing terms directly , so the matching coefficient corresponding to this case is simply the contribution from the purely hard ( h - h- ... -h ) regions as it was defined in ( [ eq : newmatch ] ) .
moreover , the purely hard regions can now be computed directly at zero external relative momentum , which simplifies the calculation .
the same discussion applies to the matching of the potentials in pnrqcd , to which we turn next .
however , while the @xmath24 terms of the nrqcd lagrangian are relevant only at nnnlo , the difference between four- and @xmath33-dimensional potentials in the lagrangian matters already at nnlo .
the @xmath33-dimensional ones must be used in conjunction with ( [ eq : newmatch ] ) as was done in @xcite .      due to the @xmath118 suppression of the p - wave correlation function ( [ eq : piatonrqcd ] ) relative to the s - wave case ( [ eq : g ] ) , the hard matching coefficient @xmath132 of the axial - vector current is needed only with one - loop accuracy for the nnnlo calculation of the top - quark pair production cross section near threshold .
the relevant expression is @xmath282 the two - loop correction is also known @xcite .
as discussed in section  [ sec : momentumregions ] to perform the all - order resummation , a second matching procedure is required , by which the soft region and potential light fields ( gluons and light quarks ) are integrated out .
this results in the potential nrqcd ( pnrqcd ) effective field theory  @xcite . in pnrqcd
the light fields are purely ultrasoft and the heavy quarks are potential , hence the terms in the effective lagrangian can be assigned a unique scaling in the velocity expansion .
the effective lagrangian relevant to the third - order calculation takes the simple form @xmath283(x+{\bf r } ) \ , v_{ab;cd } ( r,{{\mbox{\boldmath $ { \partial}$}}})\ , \big[\chi^\dag_c \chi^{}_d\big](x ) \\ & & -g_s\psi^\dagger(x){{\bf { x}}}\cdot{{\bf { e}}}(t,{{\bf { 0}}})\psi(x)- g_s\chi^\dagger(x){{\bf { x}}}\cdot{{\bf { e}}}(t,{{\bf { 0}}})\chi(x ) , \label{eq : pnrqcd}\end{aligned}\ ] ] where @xmath284 with @xmath285 the tree - level colour coulomb potential . the pnrqcd lagrangian consists of kinetic terms ( first line ; including the relativistic corrections proportional to @xmath286 ) , heavy - quark potential interactions ( second line ) and an ultrasoft interaction that contributes first at third order .
the heavy - quark potentials generated in the the matching to pnrqcd should be considered as short - distance coefficients of the interactions .
they are split into the tree - level coulomb potential , which must be treated non - perturbatively and a remainder @xmath287 , which represents a perturbation . to achieve a homogeneous velocity scaling the position argument of ultrasoft fields
should be multipole - expanded in interactions with heavy quarks @xcite , which explains the space - time argument of @xmath288 and the chromo - electric field in the ultrasoft interaction terms .    as will be discussed below no further matching of the non - relativistic vector current
is needed , that is @xmath289 to the required accuracy .
thus , instead of ( [ eq : g ] ) , we have to calculate @xmath290(x)\ , [ \psi^{\dag}\sigma^i\chi](0 ) )    evaluated to third - order in pnrqcd perturbation theory .
the dimensionally regulated pnrqcd lagrangian required for second - order calculations of heavy - quark pair production near threshold was provided in @xcite , and the explicit derivation of the ultrasoft interaction in the third line of ( [ eq : pnrqcd ] ) from nrqcd was given in @xcite .
the only new piece that is needed is the third - order heavy - quark potential in @xmath287 . in the remainder of this subsection
we first give the pnrqcd feynman rules ( when the ultrasoft interactions are neglected ) and then sketch several ways of deriving these rules and the form of the propagator .
subsequently , we summarize the heavy - quark potentials .
we also derive equation - of - motion relations that allow us to reduce the number of potential insertions to be calculated and briefly discuss the ultrasoft contribution already calculated in  @xcite .
we begin by summarizing the rules for calculating pnrqcd diagrams with insertions of potential interactions , but no ultrasoft interactions .
since the leading - order lagrangian includes the coulomb potential @xmath291 , the propagator is the one for a heavy quark anti - quark _
pair_. we draw the propagator as in the left diagram of figure  [ fig : pnrqcd ] , where the blob stands for the sum of all potential ( coulomb ) ladder diagrams , which is included in the propagator . for a pair in a colour - singlet state
each propagator gives a factor @xmath292 where @xmath19 is the non - relativistic energy of the pair and @xmath161 ( @xmath293 ) the three - momentum of the in - coming ( out - going ) quark .
for the colour - octet state the propagator is @xmath294 .
the function @xmath295 is the solution to the @xmath33-dimensional lippmann - schwinger equation for a pair in irreducible su(3 ) colour representation  @xmath296 , @xmath297 where the @xmath298 and @xmath299 for the colour - singlet and colour - octet representation , respectively .
explicit expressions will be given below .
the scale @xmath300^{1/2}$ ] is defined such that minimal subtraction of @xmath42 poles corresponds to the @xmath73 rather than ms scheme @xcite .
the fourier transform @xmath301 satisfies in @xmath233 dimensions the schrdinger equation @xmath302 in the following , when the superscript is left out , the propagator refers to the colour - singlet representation .
the vertex associated with the insertion of a perturbation potential @xmath303 in momentum space is given by @xmath304 and internal relative momenta @xmath305 are integrated over with measure @xmath306 .
note that the insertion of a potential does not change the colour - state of the quark anti - quark pair , when it is in an irreducible representation , that is the colour - singlet or the colour - octet state .
the reason for this is that @xmath307 and @xmath308 constitute a complete set of orthogonal projectors .
if the in - coming pair is , for example , in the colour - singlet state , all propagators will be colour - singlet propagators and the potential insertions are effectively projected to the colour - singlet potential @xmath309 this is different when an ultrasoft gluon is emitted , in which case the pair changes its colour state , as explicitly seen in ( [ eq : gfcorr ] ) below .
a general pnrqcd diagram with multiple insertions of perturbation potentials is therefore an expression of the form @xmath310
i{g_0}({\bf p}_1,{\bf p}_2;e ) i\delta{v}_1({\bf p}_2,{\bf p}_3 ) i{g_0}({\bf p}_3,{\bf p}_4;e)\ ,   i\delta{v}_2({\bf p}_4,{\bf p}_5 )   i{g_0}({\bf p}_5,{\bf p}_6;e)\ ,   \ldots \label{multipleinsertions}\ ] ] ( for the case of a colour - singlet state ) . for convenience of notation here and below
we often leave out factors of @xmath311 required in dimensional regularization to restore the proper dimension of the given expression .      in this section
we sketch three derivations of the rules for pnrqcd perturbation theory : diagrammatic , quantum - mechanical , and by path - integral methods . for simplicity
we assume the colour - singlet representation , but the derivation is easily generalized to an arbitrary irreducible representation .
consider the amputated amplitude of the heavy - quark scattering process @xmath312 with non - relativistic external momenta @xmath313 , @xmath314 and @xmath315 , @xmath316 in the rest frame of the @xmath107 pair .
the sum of all ( ladder ) diagrams with any number ( greater than zero ) of leading - order potential insertions @xmath291 is given in momentum space by @xmath317 \frac{(i g_s)^2 i}{({{\bf { k}}}_1-{{\bf { k}}}_0)^2}\ , \frac{(i g_s)^2 i}{({{\bf { k}}}_2-{{\bf { k}}}_1)^2}\ , \ldots \frac{(i g_s)^2 i}{({{\bf { k}}}_{n+1}-{{\bf { k}}}_n)^2}\ , \nonumber\\ & & \cdot \,\prod_{i=1}^n\ , \frac{i}{\frac{e}{2}+k_i^0-\frac{({{\bf { p}}}+{{\bf { k}}}_i)^2}{2m}+ i\epsilon}\ , \frac{-i}{\frac{e}{2}-k_i^0-\frac{({{\bf { p}}}+{{\bf { k}}}_i)^2}{2m}+ i\epsilon } \,,\end{aligned}\ ] ] where we define @xmath318 and @xmath319 .
we perform the integrations over the loop momentum zero components @xmath320 by closing the contour in the upper half plane , and pick up the residues from the poles at @xmath321 , which results in @xmath322 \,\frac{1}{{{\bf { k}}}_1 ^ 2}\ ; \nonumber\\ & & \hspace*{0.0cm}\times   \prod_{i=1}^n\ , \frac{1}{({{\bf { k}}}_{i+1}-{{\bf { k}}}_i)^2   ( e-\frac{({{\bf { p}}}+{{\bf { k}}}_i)^2}{2m}+i\epsilon)}\ , .
\qquad \label{eq : hc}\end{aligned}\ ] ] the @xmath323 term in this and the previous sum is understood as @xmath324 , which is the expression for the exchange of a single potential ( coulomb ) gluon .
@xmath325 sums the leading p - p- ... -p region to all orders .
next we multiply the propagator factors @xmath326 for the external pairs of lines and add the zero - coulomb exchange graph .
multiplying by @xmath327 this defines @xmath328 it is straightforward to show that this expression satisfies the @xmath33-dimensional lippmann - schwinger equation ( [ eq : lippmann - schwinger ] ) , and hence represents the colour - singlet coulomb green function .
the summation of ladder diagrams is therefore accomplished by associating the quantity @xmath329 with the propagator of the quark anti - quark pair and the vertex @xmath330 with the interaction potentials .
it is understood that the integrations over the zero - components of loop momenta are already done .
note that while no closed expression for the green function is known in @xmath33 dimensions , it is important that the above expression is defined in dimensional regularization , and that it can be expanded perturbatively in @xmath104 in @xmath33 dimensions .
this guarantees the consistency of the dimensional regularization procedure , which requires subtracting terms with a finite number of coulomb exchanges from @xmath331 in @xmath33 dimensions as will be seen in part ii of the paper .
the pnrqcd correlation function ( [ eq : gpnrqcd ] ) describes a quark anti - quark pair created in a colour - singlet , spin - triplet state at point 0 , which propagates and is destroyed locally at point @xmath332 . in terms of the momentum - space propagator
it is given by @xmath333\,,\qquad \label{eq : gexp}\end{aligned}\ ] ] where the potential refers to the colour - singlet potential ( [ eq : singlet - potential ] ) , and the ellipses to terms with multiple potential insertions .
both the propagator and the potential carry spin - indices in general .
however , the unperturbed pnrqcd lagrangian is spin - independent , so the propagator is diagonal in the spin indices and we drop the @xmath334 spin factor .
the spin - dependence of the perturbation potential enters only at nnlo , hence up to n@xmath335lo only one of the @xmath336 insertions can carry a non - trivial spin - dependence .
when this insertion appears in ( [ eq : gexp ] ) , @xmath336 is understood as @xmath337 and the trace must be carried out in @xmath33 dimensions .
greek indices of the potential refer to spin ( rather than colour ) .
the normalization factor and pauli matrices result from the definition of ( [ eq : gpnrqcd ] ) and correspond to the spin - triplet projection of the potential .
for the insertions of spin - independent potentials the spin - factor @xmath338 is included in ( [ eq : gexp ] )      with only potential interactions the pnrqcd lagrangian can be projected onto the fock states with a single quark and a single anti - quark without loss of content , since potential interactions do not change particle number .
we define the centre - of - mass wave function of a quark anti - quark state @xmath339 in the position - space representation via @xmath340^\dagger     simplicity , we assume again a projection on the colour - singlet representation . by reversing the steps that lead from the schrdinger equation to a second - quantized schrdinger field theory , making use of the field equation and the canonical commutation relations , we find that @xmath341    \psi(t,{{\bf { r}}})\,.\ ] ] the green function of the schrdinger operator is given by @xmath342^{-1}|{{\bf { r}}}'\rangle$ ] with @xmath343 a quark anti - quark separation eigenstate with eigenvalue @xmath344 . in operator notation
the green function is @xmath345^{-1}$ ] such that @xmath346 where @xmath347^{-1}$ ] and @xmath348 . since @xmath349 , the previous equation verifies the pnrqcd feynman rules . in this notation
the correlation function ( [ eq : gpnrqcd ] ) is given by @xmath350 which is equivalent to ( [ eq : gexp ] ) upon using ( [ eq : goperatorexp ] ) and inserting complete sets of momentum eigenstates , @xmath351 leaving the spin - average being implicit .      in @xcite the pnrqcd lagrangian is expressed in terms of composite colour - singlet and colour - octet fields . here
we provide a path - integral derivation of this formulation .
we focus on the colour - singlet field and drop the colour and spin indices of the composite field @xmath352_{x^0=y^0 } =   \left [ \psi^\dag(x)\chi(y ) \right]_{x^0=y^0}\ , , \qquad   \left [ s^\dag(y , x)\right]_{x^0=y^0 } =   \left [ \chi^\dag(y)\psi(x ) \right]_{x^0=y^0}\,.\ ] ]    the partition function @xmath353 of pnrqcd is defined as @xmath354 where we use the leading - order pnrqcd lagrangian @xmath355\psi(x ) + \chi^\dag(x)\bigg[i\partial_0 - \frac{{\mbox{\boldmath $ { \partial}$}}^2}{2m}\bigg]\chi(x ) \nonumber\\ & & -\int d^4y \left[\psi^\dag(x)\chi(y)\right ] v(x , y )   \left[\chi^\dag(y ) \psi(x)\right],\end{aligned}\ ] ] with the quark anti - quark potential @xmath356 .
the derivation remains however valid with an arbitrary potential , not just the leading - order coulomb potential .
the different sign of the potential term relative to ( [ eq : pnrqcd ] ) is due to the different order of fermion fields .
the composite field is introduced by means of @xmath357_{x^0=y^0 } -   \left[\psi^\dag(x)\chi(y)\right]_{x^0=y^0}\right )   \nonumber\\ & = & \int { \cal d}s \,{\cal d}\sigma\ , \exp\bigg\ { \,i \!\int d^4x \int d^4y \ , \sigma(x , y ) \delta(x^0-y^0 )   \left ( s(x , y)- \left[\psi^\dag(x)\chi(y)\right]\right ) \bigg\}\ , , \qquad\end{aligned}\ ] ] where the second line is the fourier representation of the functional delta - function which introduces the auxiliary field @xmath358 . with a similar formula for the complex conjugate field we can rewrite the partition function as @xmath359 where @xmath360\right ) \nonumber\\ & & + \,\sigma^\dag(y , x)\left ( s^\dag(y , x)- \left[\chi^\dag(y)\psi(x)\right ] \right)\,.\end{aligned}\ ] ] to make equations concise we introduced @xmath361 and @xmath362    our task is to integrate over @xmath363 and their conjugates to obtain the lagrangian for the @xmath252 field . in the following
we do this step by step obtaining the sequence of lagrangians @xmath364 we drop field - independent factors which can be absorbed into the path - integral measure .
first we integrate out @xmath119 and its conjugate by completing squares in the exponent : @xmath365_{\psi=0 } \nonumber\\ & & \qquad + \ , \left ( \psi^\dag- \chi^\dag\cdot\sigma^\dag\cdot k_q^{-1 } \right )   \cdot k_q\cdot \left(\psi- k_q^{-1 } \cdot\sigma\cdot \chi \right ) -\chi^\dag\cdot\sigma^\dag\cdot k_q^{-1 } \cdot \sigma \cdot \chi\ , .
\qquad\end{aligned}\ ] ] in the second line the abbreviation @xmath366 stands for the integration of the adjoining variables , such that , for example , the last term reads explicitly @xmath367 the inverse operators are given by @xmath368 changing variables to @xmath369 ( similarly for @xmath370 ) and integrating over @xmath119 , @xmath370 we obtain @xmath371_{\psi=0 } -\chi^\dag\cdot\sigma^\dag\cdot k_q^{-1 } \cdot \sigma \cdot \chi . \qquad\end{aligned}\ ] ] the integration over @xmath372 is done analogously resulting in the partition function @xmath373 the determinant contains the @xmath358 field . to write it as a term in the lagrangian we use @xmath374 and keep only the bilinear term in the @xmath375 field in the expansion of the logarithm .
we comment on the other terms below .
the lagrangian after this procedure is @xmath376 now we integrate over the delta - functions in the time coordinates implicit in the definition of the @xmath375 field and obtain @xmath377 where @xmath378 represents the coordinates of the quark and anti - quark at coincident time @xmath379 .
the fields with argument @xmath380 are defined as @xmath381_{x^0=y^0 } \nonumber\\ \sigma(z)&=&\left [ \sigma(x , y)\right]_{x^0=y^0 } , \nonumber\\ k_\sigma(z',z)&=&\left[k_{q}^{-1}(x',x ) k_{\bar{q}}^{-1}(y , y ' ) \right]_{t = x^0=y^0 , \,t'=x^{0'}=y^{0 ' } } \label{eq : ksig}\end{aligned}\ ] ] the last step consists of performing the gaussian integral over @xmath358 .
the inverse of @xmath382 is defined by @xmath383 resulting in @xmath384 to compute @xmath385 we use ( [ eq : ksig ] ) and the definitions ( [ eq : inversekq ] ) .
the integrals over the zero - components of the two momenta can be written as integrals over relative momentum @xmath386 and total momentum @xmath387 . since
the exponentials are independent of @xmath387 , the @xmath387 integral can be performed by contour integration which gives @xmath388 where @xmath389 .
it follows that @xmath390\,,\ ] ] and therefore @xmath391 with @xmath392 after separating the free centre - of - mass motion this represents the pnrqcd lagrangian expressed in terms of the composite quark anti - quark field .    when one keeps the higher - order terms in the expansion of the logarithm in ( [ eq : expdet ] ) the path - integral over @xmath358
can no longer be done exactly .
expanding the quartic and higher - order terms in the exponential , we obtain vertices involving four and more @xmath252 fields , which describe scattering of composite fields
. these terms are clearly not relevant to the threshold dynamics of a single quark anti - quark pair .      in four dimensions
explicit solutions for the the schrdinger equation ( [ eq : schroedinger ] ) can be found , equivalent to the sum of diagrams ( [ eq : coulombfn ] ) .
we quote the results for the colour - singlet green function .
the general case is obtained by substituting @xmath393 everywhere .
the momentum space pnrqcd propagator ( coulomb green function ) can be expressed in the form @xcite @xmath394 which closely resembles ( [ eq : coulombfn ] ) and shows that the sum from @xmath395 to infinity in ( [ eq : hc ] ) can be transformed into a remarkably simple integral . at this point
we omit the @xmath70 prescription on @xmath69 and regard @xmath331 as a function of a complex energy variable , which has a cut for @xmath396 and isolated poles on the negative real axis .
the variable @xmath71 equals @xmath397 as defined in ( [ eq : deflambda ] ) .
the first line of ( [ eq : schwingerrep ] ) separates the zero- and one - coulomb gluon exchange terms . in practice
, we find it simpler to perform the all - order summation in the position space representation , where the potential insertions take a simple multiplicative ( rather than convolutive ) form , and therefore we do not make use of the above representation in the calculation in paper  ii .
an integral representation for the position space coulomb green function is @xmath398^\lambda [ t(s-1)]^{-\lambda}\ , \nonumber\\ & & \times   \frac{\partial^2}{\partial t\partial s }   \left(\frac{ts}{|s{{\bf { r}}}-t{{\bf { r}}}^\prime| } \ , e^{-\sqrt{-me}\,((1-t)r^\prime+(s-1)r+|s{{\bf { r}}}-t{{\bf { r}}}^\prime| ) } \right),\end{aligned}\ ] ] valid for @xmath399 , where @xmath400 , @xmath401  @xcite . for @xmath402 exchange @xmath403 in the above expression .
putting one of the arguments to zero , this simplifies to @xmath404 which depends only on @xmath400 .
we use this form of the coulomb green function mainly for propagators connecting to the external current vertex , in which case ( [ eq : greenint ] ) applies .    for the general case of a propagator in between two potential insertions the representation of the position - space green function in terms of laguerre polynomials @xmath405 @xcite turns out to be most useful . in this representation one
first performs a partial wave expansion @xmath406}(r , r^\prime;e)\ , , \label{eq : gplexpand}\ ] ] where @xmath407 are the legendre polynomials .
the partial - wave green functions read @xmath408}(r , r^\prime;e ) =   \frac{m p}{2\pi } \ , ( 2 p r)^l ( 2 p r^\prime)^l\ , e^{-p(r+r^\prime ) }   \sum_{s=0}^{\infty}\ , \frac{s!\,l_s^{(2l+1)}(2pr ) l_s^{(2l+1)}(2 p r^\prime ) } { ( s+2l+1 ) ! ( s+l+1-\lambda)}\ , , \label{eq : gpartial}\ ] ] where @xmath409 , and the laguerre polynomials are defined by @xmath410.\ ] ] since at nnnlo accuracy the potential insertions can not change the angular momentum of the quark anti - quark pair and since the production current @xmath241 creates an s - wave state , we only need the @xmath411 green function to compute the potential contributions to the pnrqcd correlation function ( [ eq : gpnrqcd ] ) .
the p - wave green function is required to compute the ultrasoft contribution  @xcite and the contribution ( [ eq : qcdaxialvectorcurrent ] ) from the p - wave production current .
a general property of the coulomb interaction is that the ultraviolet behaviour of the ladder diagrams improves with the number of exchanges .
thus , when the external current or potential insertions cause uv divergences , it is necessary to subtract only the first few terms in the sum of ladder diagrams .
the divergent diagrams must be done in @xmath33 dimensions using standard methods , while for the convergent remainder one of the above expressions , properly subtracted , can be used .
we therefore use the notation @xmath412 for example , the three terms in ( [ eq : schwingerrep ] ) correspond to @xmath413 .    from ( [ eq : gexp ] ) it follows that the leading term in @xmath21 equals @xmath414 , which is , however , divergent as can be seen from ( [ eq : greenint ] ) . to compute @xmath415 in dimensional regularization
, we note that the zero - coulomb exchange term is linearly divergent , the one - coulomb exchange logarithmically , and the remainder is convergent .
we therefore compute @xmath416 from the first line of ( [ eq : schwingerrep ] ) , which yields @xmath417 .
\label{eq : g01ex}\end{aligned}\ ] ] the remaining terms can be calculated from a modified version of ( [ eq : greenint ] ) with @xmath418 and integrand @xmath419 . after subtracting the first two terms of the @xmath13 expansion
the result is finite as the regulator @xmath420 and gives the remaining contribution @xmath421 .
the final result for the @xmath73 subtracted zero - distance green function @xcite is @xmath422\ , .
\label{eq : g00msbar}\ ] ] the poles of the euler psi - function at positive integer @xmath71 correspond to the s - wave quark anti - quark bound states . near the bound - state poles the green function takes the form @xmath423 where @xmath424 is the wave - function at the origin of the @xmath250th bound state with energy @xmath425 .
the imaginary part of the green function for @xmath396 is known as the sommerfeld factor  @xcite .
explicitly , the imaginary part below and above threshold is given by @xmath426 for real energies @xmath427 .
this expression can be used in ( [ r1 ] ) to obtain the leading - order approximation to the resummed top pair production cross section in the threshold region for zero width of the top quark .
we now summarize the momentum - space potentials required for the nnnlo calculation with the pnrqcd lagrangian ( [ eq : pnrqcd ] ) .
only the colour - singlet projection @xmath428 of the general quark anti - quark potential @xmath429 is relevant to top - quark pair production through the electromagnetic and electroweak current , and this will be given in the following .
the various potential terms can be ordered in a @xmath133 expansion , beginning with the coulomb potential of order @xmath430 . allowing for the spin - dependence from order @xmath129
, we write the singlet - potential in the general form : @xmath431\otimes[\sigma_i,\sigma_j ]   \nonumber   \\   & &    -\,{\cal v}_{p}(\alpha_s)\frac{2\pi c_f   \alpha_s({\bf p}^2+{\bf p}^{\prime 2})}{m^2{\bf{q}}^2 } + { \cal v}_{hf}(\alpha_s)\frac{\pi c_f \alpha_s}{4m^2{\bf q}^2 }   [ \sigma_i,\sigma_j]q_j\otimes [ \sigma_i,\sigma_k]q_k \nonumber   \\   & &   - { \cal v}_{so}(\alpha_s)\frac{3\pi c_f \alpha_s}{2m^2 { \bf q}^2 }   \bigg([\sigma_i,\sigma_j]q_i p_j   \otimes 1   -1\otimes [ \sigma_i,\sigma_j]q_i p_j\bigg ) + \ldots .
\label{eq : potentialbeforespin } \end{aligned}\ ] ] the coefficients @xmath432 of the potentials are @xmath13 dependent : @xmath433 the above representation of the potential therefore corresponds to an expansion in @xmath13 and @xmath12 , since the potential momenta @xmath161 , @xmath293 and @xmath434 are of order @xmath435 . with @xmath436
the leading term is of order @xmath437 .
the ellipses denote terms of order @xmath438 , which would contribute from n@xmath439lo .
a notation has been used which is valid in @xmath33 dimensions by avoiding the use of vector products or the totally antisymmetric @xmath164 tensor that would arise from using the three - dimensional identity for the commutator of pauli matrices .
the coefficients of the potentials are chosen such that the leading - order coefficients are either one or zero .
the tensor products @xmath440 refer to the spin matrices on the quark ( @xmath441 ) and anti - quark line ( @xmath442 ) . for the first three and the fifth terms of ( [ eq : potentialbeforespin ] ) , which are spin - independent
, we omitted the trivial @xmath443 factor .
the on - shell matching calculation of the potential coefficients @xmath444 , @xmath445 , @xmath446 , @xmath447 results in infrared ( ir ) divergences @xcite , which are related to ultraviolet divergences in the calculation of the ultrasoft correction .
it is convenient to subtract these divergences in the results given below and add them back to the ultrasoft calculation ( see @xcite and section [ sec : ultrasoft ] ) .
the subtraction term , which is _ added _ to the potential , is @xmath448 .
\label{eq : counter}\end{aligned}\ ] ]    in addition to potential insertions , the relativistic correction to the kinetic energy term @xmath449 in ( [ eq : pnrqcd ] ) needs to be included in pnrqcd perturbation theory .
formally , this can be done by adding @xmath450 to the potential . since @xmath451
, it counts as a nnlo potential .
the delta function eliminates the momentum integration that is associated with a potential insertion , but which is not present for the kinetic energy correction to a ( anti-)quark propagator .
we now present the results for the potential coefficients , starting with the coulomb potential .
the coefficient @xmath452 encodes the quantum corrections to the coulomb potential , which are needed up to the three - loop order .
the insertions of coulomb potentials are finite , so we do not need the @xmath33-dimensional expression of the potential , as long as only coulomb potential insertions are considered .
this reflects the fact that the schrdinger equation with the @xmath28 potential is non - singular and could be solved exactly , without referring to pnrqcd perturbation theory , as was done , for instance in @xcite . however , in the third - order computation of the top anti - top production cross section also the double insertion of the nlo coulomb potential together with the singular insertion of a nnlo non - coulomb potential has to be taken into account ; hence the order @xmath164 part of the one - loop coulomb potential multiplies a divergent quantity and contributes to the final result .
the coefficient @xmath453 is therefore given with the full @xmath164 dependence .
the first four terms in the expansion of the coulomb potential can be represented in the form @xmath454\ , \frac{\beta_0}{\epsilon}\ , + \bigg(\frac{\mu^2}{{{\bf { q}}}^2 } \bigg)^{\!\epsilon}\ , a_{1}(\epsilon ) , \label{eq : vcoulombnlo}\\ { \cal v}_{c}^{\,(2 ) } & = & a_2 +   \left(2a_1\beta_0+\beta_1\right)\ln\frac{\mu^2}{{{\bf { q}}}^2 } + \beta_0 ^ 2\ln^2\frac{\mu^2}{{{\bf { q}}}^2 } , \\ { \cal v}_{c}^{\,(3 ) } & = & a_3 +   \left(2a_1\beta_1+\beta_2 + 3a_2\beta_0 + 8\pi^2c_a^3\right ) \ln\frac{\mu^2}{{{\bf { q}}}^2 } \nonumber\\   & & + \left(\frac{5}{2}\beta_0\beta_1 + 3a_1\beta_0 ^ 2\right ) \ln^2\frac{\mu^2}{{{\bf { q}}}^2}+\beta_0 ^
3\ln^3\frac{\mu^2}{{{\bf { q}}}^2 } , \label{eq : vcoulombn3lo}\end{aligned}\ ] ] with @xmath455 - 4\,t_fn_f\bigg)\ , \frac{e^{\gamma_e\epsilon}\,\gamma(1-\epsilon)\,\gamma(2-\epsilon)\ , \gamma(\epsilon)\ , } { ( 3 - 2\epsilon)\,\gamma(2 - 2\epsilon ) } -\frac{\beta_0}{\epsilon } , \\ & & a_{2 } = \left(\frac{4343}{162}+4\pi -\frac{\pi^{2}}{4 } + \frac{22}{3}\zeta_{3}\right)c_{a}^{2 } -\left(\frac{1798}{81}+\frac{56}{3}\zeta_{3}\right ) c_{a}t_{f}n_f \nonumber\\ & & \hspace*{1 cm }   -\,\left(\frac{55}{3}-16\zeta_{3}\right ) c_f t_{f}n_f +   \frac{400}{81}\,(t_f n_f)^2.\end{aligned}\ ] ] the one - loop correction @xmath456 has been known in @xmath233 for some time @xcite .
we computed the @xmath33-dimensional expression and confirmed the result first shown in  @xcite .
the term @xmath457 without a momentum factor arises from the charge renormalization counterterm and the square bracket in ( [ eq : vcoulombnlo ] ) produces the associated logarithm .
however , the expansion in @xmath164 can only be done after the momentum integrals of the potential insertion are performed , if these integrals are divergent .
the two - loop coefficient @xmath458 has first been calculated in @xcite and correctly in @xcite .
the three - loop coefficient @xmath32 is also known  @xcite , though not all colour coefficients have been computed analytically : @xmath459 ( n_f t_f)^2   \nonumber \\ & & + \bigg[\ ,    \,\left(-709.717\right ) c_a^{\,2 }   + \left(-\frac{71281}{162}+264\zeta_3 + 80\zeta_5\right ) c_a c_f \nonumber\\ & & \hspace{1 cm }   + \left(\frac{286}{9}+\frac{296}{3}\zeta_3 - 160\zeta_5\right ) c_f^{\,2 } \bigg]\,n_f t_f \nonumber \\ & & + \bigg [   -56.83(1)\ , \biggl(\frac{d_f^{abcd}d_f^{abcd}}{2n_a\,t_f}\biggr)\ , \bigg ] n_f \nonumber \\ & & + \bigg[\ , 502.24(1)\ , c_a^{\,3 } + \,\left(-136.39(12)\right)\ ,   \biggl(\frac{d_a^{abcd}d_f^{abcd}}{2n_a\,t_f}\,\biggr ) \bigg]\ , .
\label{eq : a3}\end{aligned}\ ] ]
the numerical coefficients are quoted from  @xcite , and the color factors @xmath460 for @xmath461 are given by @xmath462 parts of the nnnnlo coulomb potential are also known @xcite , but not needed for the third - order cross section calculation .
the third - order coulomb potential has an ir divergence  @xcite , which cancels against a divergence in the calculation of the nnnlo ultrasoft calculation .
the corresponding @xmath42 pole is subtracted and therefore does not appear in ( [ eq : vcoulombn3lo ] ) , but the logarithmic part multiplied by @xmath463 in this equation comes from this divergence .
note that the coefficient of this logarithm agrees with @xcite , but is three times larger than the one in @xcite .
the reason for this is that here as in @xcite all three loops are computed in @xmath33 dimensions not just the divergent one , as is required for consistency with the ultrasoft calculation . hence the divergence related to the coulomb potential in ( [ eq : counter ] ) is multiplied by @xmath464 rather than @xmath465 .
the coefficient of the @xmath466 potential is generated first at the one - loop order , where it is suppressed by @xmath467 relative to the leading coulomb potential .
hence the two - loop coefficient is required for the nnnlo calculation of the cross section .
the insertions of this potential cause ultraviolet divergences such that we need the one - loop coefficient to @xmath468 and the two - loop one to @xmath24 .
up to the two - loop order we can represent the @xmath469 potential in the form @xmath470\left(-\frac{8}{3\epsilon}\right )             \bigg(2 c_f c_a + c_a^2 \bigg)\ , \nonumber \\ & & + \ , \bigg[\bigg(\frac{\mu^2}{{{\bf { q}}}^2}\bigg)^{\!2\epsilon }      - \bigg ( \frac{{\mu } ^2}{{{\bf { q}}}^2}\bigg)^{\!\epsilon}\ ,    \bigg]\ , \frac{2\beta_0}{\epsilon } \ , b_{1}(\epsilon)\ , + \bigg(\frac{\mu^2}{{{\bf { q}}}^2 } \bigg)^{\ ! 2\epsilon } \ , 4 b_{2}(\epsilon)\ , , \label{eq : vb2}\end{aligned}\ ] ] with @xmath471        -c_a\,[1-\epsilon ] \right)\ , \frac { e^{\gamma_e\epsilon}\,\gamma(\frac{1}{2}-\epsilon)^{2 } \gamma(\frac{1}{2}+\epsilon)}{\pi^{\frac{3}{2}}\,\gamma(1 - 2\epsilon)}\ , , \label{b1eps}\\ b_{2}(\epsilon ) & = &   \left[\frac{65}{18}-\frac{8 } { 3}\ln{2}\right ] c_ac_f -\left[\frac{101}{36}+\frac{4 } { 3}\ln{2}\right ] c_a^{2 } \nonumber \\ & & + \left[\frac{49}{36}\,c_a -\frac{2}{9 } c_f\right ] t_f n_f + \epsilon \,b_{2}^{(\epsilon)}+o(\epsilon^2)\,.\end{aligned}\ ] ] once again we subtracted the @xmath42 ir pole that remains after charge renormalization by adding the relevant part of @xmath472 from ( [ eq : counter ] ) .
hence ( [ eq : vb2 ] ) is finite ; however , the expansion in @xmath164 must be performed only after the computation of the potential insertion .
the one - loop expression @xmath473 has been computed in @xmath33 dimensions  @xcite , and the four - dimensional value @xmath474 of the two - loop coefficient is quoted from @xcite .
the @xmath24 term of @xmath475 has been parameterized above by @xmath476 .
while known , it has not yet been published @xcite . in the numerical results in part
ii of the paper we use the estimate @xmath477 , which is significantly larger but more realistic than our previous estimate @xmath478 ( @xmath479 ) @xcite .
the coefficients of the @xmath480 potentials are generated at tree level , where they are suppressed by @xmath118 relative to the leading coulomb potential .
hence the one - loop coefficients are required for the nnnlo calculation of the cross section .
the tree - level coefficients are : @xmath481 as can be seen from ( [ eq : potentialbeforespin ] ) , spin - dependence arises first within the @xmath480 potentials .
the insertions of these potentials are again ultraviolet divergent .
we therefore need the @xmath24 term of the one - loop coefficients .
these are available only from @xcite .
we computed the @xmath33-dimensional expressions and confirmed the previous result .
the spin - projected expression has already been given in our previous work @xcite .
nrqcd tree level diagrams of order @xmath129 . dashed ( curly )
lines denote the @xmath137 ( @xmath138 ) gluon field .
the number @xmath482 at the vertex refers to the nrqcd interaction with coefficient function @xmath225 , see also figure  [ fig : nrqcd_fr ] and ( [ eq : nrqcdqq ] ) , ( [ eq : nrqcdqqqq ] ) , ( [ eq : nrqcdg ] ) .
symmetric diagrams are not shown.,height=132 ]    the complete one - loop coefficients consist of two different contributions .
the first arises from the one - loop correction to the nrqcd couplings @xmath225 , already discussed in the previous section .
we call this the `` hard '' contribution
. it can be extracted from the tree diagrams in figure [ fig : nrqcdtree ] with one - loop corrected nrqcd vertices .
the second contribution arises from explicitly integrating out soft loops .
potential as discussed in section  [ subsec : vectormatch ] below . however , there is no such contribution to the @xmath129 potentials at one loop .
the reason for this is that in the case of the @xmath133 potential the relevant contribution arises from the box integral with all vertices of the leading - order @xmath483 type , expanded to subleading order in the potential region .
the next correction is always suppressed by _ two _ powers of @xmath12 ( for instance , from replacing one of the vertices by a @xmath265 vertex from the nrqcd lagrangian ) , and hence can contribute only at order @xmath160 . ]
the corresponding one - loop nrqcd diagrams are shown in figure  [ fig : nrqcdloop ] .
these contributions will be called `` soft '' . the total potential including tree and one - loop correction is then @xmath484 .
the hard one - loop contributions can be easily calculated with the feynman rules presented in figure [ fig : nrqcd_fr ] , since all required @xmath33-dimensional nrqcd matching coefficients are already known at one - loop order from the previous section .
the result reads : @xmath485 + o(\alpha_s^2 ) , \nonumber\\ { \cal v}_{p}^{(hard)}(\alpha_s ) & = & 1+o(\alpha_s^2 ) , \\ { \cal v}_{so}^{(hard)}(\alpha_s ) & = & \frac{1}{3}(2d_1+d_3)+o(\alpha_s^2 ) \\ & = & \nonumber 1+\frac{\alpha_s}{\pi}\bigg(\frac{\mu^2}{m^2}e^{\gamma_e}\bigg)^{\epsilon } \gamma(\epsilon ) \bigg[\frac{c_a(2\epsilon^2 - 1)-2c_f\epsilon(2\epsilon+1)}{6\epsilon-3}\bigg ] + o(\alpha_s^2 ) , \\ { \cal v}_{hf}^{(hard)}(\alpha_s ) & = & d_1 ^ 2+o(\alpha_s^2 ) \\ & = & \nonumber 1+\frac{\alpha_s}{\pi}\bigg(\frac{\mu^2}{m^2}e^{\gamma_e}\bigg)^{\epsilon } \gamma(\epsilon ) \bigg[\frac{c_a(2\epsilon^2 - 1)-2c_f\epsilon(2\epsilon+1)}{4\epsilon-2}\bigg ] + o(\alpha_s^2 ) , \\ { \cal v}_{s}^{(hard)}(\alpha_s ) & = & \frac{1}{2\pi c_f \alpha_s}(d_{sv}+c_f d_{vv})+o(\alpha_s^2 ) \label{eq : vm2hlast } \\ & = & \nonumber \frac{\alpha_s}{\pi}\bigg(\frac{\mu^2}{m^2}e^{\gamma_e}\bigg)^{\epsilon } \gamma(\epsilon ) \bigg[-\frac{c_a}{8}+\frac{\epsilon}{4\epsilon+2}c_f\bigg]+o(\alpha_s^2)\ , .\end{aligned}\ ] ] in this notation the tree - level value of the potential coefficient is included and assigned to the hard contribution .
nrqcd one - loop diagrams of order @xmath129 ( non - vanishing diagrams shown only , symmetric diagrams not displayed).,height=377 ]    the soft contributions come from diagrams of order @xmath129 shown in figure [ fig : nrqcdloop ] , from diagrams of order @xmath430 , where the denominator of a propagator has been expanded to higher orders as appropriate to the soft region , from soft self - energy insertions containing gluons and light quarks into the tree diagrams , and from charge renormalization counterterms .
the final result is : @xmath486+o(\alpha_s^2 ) , \label{eq : vm2s}\\ { \cal v}_{p}^{(soft)}(\alpha_s ) & = & \frac{\alpha_s}{4\pi } \bigg[\bigg(\frac{\mu^2}{{\bf q}^2}\bigg)^{\!\epsilon } \frac{(-\epsilon)\gamma(-\epsilon)^2\gamma(\epsilon ) e^{\gamma_e\epsilon}\big(c_a(56\epsilon^2 - 121\epsilon+57 ) + 12n_f t_f(\epsilon-1)\big ) } { 6(4\epsilon^2 - 8\epsilon+3)\gamma(-2\epsilon ) }   \nonumber \\    & & -\frac{\beta_0}{\epsilon}\bigg]+o(\alpha_s^2 ) , \\ { \cal v}_{so}^{(soft)}(\alpha_s ) & = & \frac{\alpha_s}{4\pi } \bigg[\bigg(\frac{\mu^2}{{\bf q}^2}\bigg)^{\!\epsilon } \frac{-\epsilon\gamma(-\epsilon)^2\gamma(\epsilon)e^{\gamma_e\epsilon } } { 6(4\epsilon^2 - 8\epsilon+3)\gamma(-2\epsilon ) } \big(c_a(d_3(8\epsilon^2 - 19\epsilon+11 ) \nonumber \\ & & + 2d_1(8\epsilon^2 - 15\epsilon+5))+ 4(2d_1+d_3)n_f t_f(\epsilon-1)\big)-\frac{2d_1+d_3}{3}\frac{\beta_0}{\epsilon}\bigg ] \nonumber \\   & & + o(\alpha_s^2 ) , \\
{ \cal v}_{hf}^{(soft)}(\alpha_s ) & = & \frac{\alpha_s}{4\pi } \bigg[\bigg(\frac{\mu^2}{{\bf q}^2}\bigg)^{\!\epsilon } \frac{(-d_1 ^ 2)\epsilon(\epsilon-1 ) \gamma(-\epsilon)^2\gamma(\epsilon)e^{\gamma_e\epsilon } \big(c_a(4\epsilon-5)+ 4n_f t_f\big ) } { ( 8\epsilon^2 - 16\epsilon+6)\gamma(-2\epsilon)}- d_1 ^ 2\frac{\beta_0}{\epsilon}\bigg ] \nonumber \\   & & + o(\alpha_s^2 ) , \\ { \cal v}_{s}^{(soft)}(\alpha_s ) & = & \frac{\alpha_s}{4\pi } \bigg(\frac{\mu^2}{{\bf q}^2}\bigg)^{\!\epsilon } \frac{d_1 ^ 2\epsilon\gamma(-\epsilon)^2\gamma(\epsilon)e^{\gamma_e\epsilon } c_a } { ( 8\epsilon-4)\gamma(-2\epsilon)}+o(\alpha_s^2 ) .
\label{eq : vm2slast}\end{aligned}\ ] ] here the @xmath225 can be set to their tree - level value 1 , since the entire expressions are already of order @xmath13 .
an alternative method of calculating the soft contribution is to extract the soft region directly from the qcd diagrams . in this way , there are fewer diagrams than in the nrqcd calculation ( but more terms from the expansion ) .
this calculation has also been performed and we checked that the results are exactly the same . note that with the second method we can not get the result in terms of the hard matching coefficients , but we obtain ( [ eq : vm2s ] ) to ( [ eq : vm2slast ] ) with @xmath487 .
some remarks are in order on the pole structure of the hard and soft @xmath129 potential coefficients .
we first note that we did not yet add the relevant contribution from the counterterm ( [ eq : counter ] ) , so the expressions given are still divergent .
the @xmath42 pole in the hard contribution is infrared , while the soft contribution contains both infrared and ultraviolet divergences .
it is instructive to separate the two : @xmath488_{\rm ir } \nonumber\\ & & + \,\bigg(\frac{\mu^2}{{\bf q}^2}\bigg)^{\!\epsilon } \left ( \left[\frac{13}{6 } c_a+\frac{2}{3 } c_f\right]_{\rm uv } + \left[\frac{8}{3 } c_a-\frac{16}{3 } c_f\right]_{\rm ir }   \right ) \bigg\}+o(\alpha_s^2 ) , \\ { \cal v}_{p}(\alpha_s ) & = & 1 + \frac{\alpha_s}{4\pi \epsilon }   \bigg(\frac{\mu^2}{{\bf q}^2}\bigg)^{\!\epsilon } \left[\frac{8}{3}c_a \right]_{\rm ir }   + o(\alpha_s^2 ) , \\ { \cal v}_{so}(\alpha_s ) & = & 1 + \frac{\alpha_s}{4\pi \epsilon } \left\ { \bigg(\frac{\mu^2}{m^2}\bigg)^{\!\epsilon }   \left[\frac{4}{3 } c_a\right]_{\rm ir } + \bigg(\frac{\mu^2}{{\bf q}^2}\bigg)^{\!\epsilon } \left[-\frac{4}{3 } c_a\right]_{\rm uv } \right\}+o(\alpha_s^2 ) , \\ { \cal v}_{hf}(\alpha_s ) & = & 1 + \frac{\alpha_s}{4\pi \epsilon } \left\ { \bigg(\frac{\mu^2}{m^2}\bigg)^{\!\epsilon }   \left[2 c_a\right]_{\rm ir } + \bigg(\frac{\mu^2}{{\bf q}^2}\bigg)^{\!\epsilon } \left[-2 c_a\right]_{\rm uv }   \right\}+o(\alpha_s^2 ) , \\
{ \cal v}_{s}(\alpha_s ) & = & \frac{\alpha_s}{4\pi \epsilon } \left\ { \bigg(\frac{\mu^2}{m^2}\bigg)^{\!\epsilon }   \left[-\frac{c_a}{2}\right]_{\rm ir } + \bigg(\frac{\mu^2}{{\bf q}^2}\bigg)^{\!\epsilon } \left[\frac{c_a}{2 } \right]_{\rm uv }   \right\}+o(\alpha_s^2 ) .
\label{eq : vm2polelast}\end{aligned}\ ] ] a term @xmath489 that was included in ( [ eq : vm2s ] ) to ( [ eq : vm2slast ] ) , which is related to the logarithms from charge renormalization , has been omitted from the expression in curly brackets for all potentials except the last one , which has no tree - level term .
we see that the ir poles from the hard region cancel the uv poles from the soft region , as it should be , since these singularities arise from hard - soft factorization .
the remaining ir singularities in the soft contribution appear only in the spin - independent potentials @xmath490 , @xmath491 .
they are precisely of the form of the remaining terms in the subtraction term ( [ eq : counter ] ) and therefore cancel with uv divergences in the ultrasoft calculation .
again , the structure of divergences is as required by consistency , since the ultrasoft contribution is spin - independent at nnnlo .
since the spin - dependent potentials appear first at nnlo , their double insertion is of higher order than nnnlo .
hence , when one computes the correlation function of the spin - triplet current ( [ eq : gpnrqcd ] ) , or the corresponding spin - singlet one , the spin - algebra can effectively be performed before the computation by working with spin - projected potentials . given a potential with spin - dependence @xmath440 , where @xmath441 ( @xmath442 ) refers to the spin - matrix on the quark ( anti - quark ) line , we replace @xmath492 note that the traces must be performed in @xmath150 space dimensions .
only the spin - triplet projection is relevant to the third - order top production cross section in @xmath0 collisions .
for the three spin - dependent terms in the general potential ( [ eq : potentialbeforespin ] ) , the projections result in @xmath493q_i p_j \otimes 1   -1\otimes [ \sigma_i,\sigma_j]q_i p_j\rightarrow 0 , \\   & & [ \sigma_i,\sigma_j]q_j\otimes [ \sigma_i,\sigma_k]q_k \rightarrow    \frac{10 - 7d+d^2}{1-d}\ , 4{\bf q}^2 , \\   & & [ \sigma_i,\sigma_j]\otimes[\sigma_i,\sigma_j ] \rightarrow ( -4)(10 - 7d+d^2)\,,\end{aligned}\ ] ] where we have omitted the trivial @xmath443 dependence as done earlier .
after the projection the four potentials @xmath494 , @xmath495 , @xmath496 , @xmath497 in ( [ eq : potentialbeforespin ] ) can be merged into a single expression @xmath498 , and we arrive at the spin - triplet , colour - singlet potential already presented in  @xcite : @xmath499.\end{aligned}\ ] ] the coulomb and @xmath133 potentials are as given earlier .
the two @xmath129 terms read , up to the one - loop order : @xmath500\ , \frac{1}{\epsilon}\ ,     \bigg(\frac{8}{3}c_a+\beta_0\bigg )   + \bigg ( \frac{\mu^2}{\bf{q}^2 } \bigg)^{\!\epsilon}\ ,          v^{(1)}_p(\epsilon)\ , , \label{eq : vp1 } \\[0.2 cm ]   { \cal v}_{1/m^2}^{\,(0)}&= & v_0(\epsilon ) = -\frac{4-\epsilon-2\,\epsilon^2}{6 - 4\epsilon}\ , , \\
\label{eq : vc1 } { \cal v}_{1/m^2}^{(1 ) } & = &     \bigg[\bigg(\frac{\mu^2}{\bf{q}^2 } \bigg)^{\!\epsilon } -1     \bigg]\ ,     \frac{1}{\epsilon}\ ,     \bigg(\,\frac{7}{3}c_f-\frac{11}{6}c_a\ , + \beta_0\ , v _ {       0}(\epsilon)\bigg ) \nonumber \\ & & + \,\bigg[\left(\frac{{\mu } ^2}{m^2 } \right)^{\!\epsilon}-1     \bigg]\ , \frac{1}{\epsilon}\ ,     \bigg(\,\frac{c_f}{3}+\frac{c_a}{2}\,\bigg )    + \bigg(\frac{\mu^2}{\bf{q}^2 } \bigg)^{\!\epsilon}\ ,          v^{(1)}_{q}(\epsilon )    + \bigg(\frac{{\mu } ^2}{m^2 } \bigg)^{\!\epsilon}\ ,          v^{(1)}_{m}(\epsilon)\,.\quad\end{aligned}\ ] ] the one - loop coefficients ( expanded up to @xmath24 ) are given by @xmath501 v^{(1)}_{q}(\epsilon ) & = & -\frac{c_f}{3 } -\frac{11}{27 } c_a + \frac{40}{27 } \,t_f n_f   \label{eq : epsterms}\\   & & + \ , \epsilon \left\{\bigg(-\frac{419}{81}+\frac{77\pi^2}{216}\bigg)c_a + \bigg(2 -\frac{7\pi^2}{36 } \bigg)c_f+\bigg(\frac{274}{81 } -\frac{2\pi^2}{27}\bigg )
t_f n_f \right\}+o(\epsilon^2 ) \ , , \nonumber \\[0.2 cm ] v^{(1)}_{m}(\epsilon ) & = & -\frac{c_f}{3 } -\frac{29}{9 } c_a + \frac{4}{15}\,t_f + \epsilon \left\{\bigg(\frac{379}{54 } + \frac{\pi^2}{24}\bigg ) c_a+\bigg(-10 + \frac{\pi^2}{36}\bigg ) c_f \right\ } \nonumber\\ & & + \,o(\epsilon^2)\,,\end{aligned}\ ] ] where now the @xmath129 pieces of the subtraction term ( [ eq : counter ] ) have been added so that there are no @xmath42 poles .
the four - dimensional expressions @xmath502 for @xmath503 agree with those obtained from @xcite , and the @xmath24 term agrees with @xcite .
having determined the matching coefficients of the pnrqcd lagrangian we now return to the question whether the nrqcd spin - triplet current is renormalized when it is matched to pnrqcd .
in general , we may write , in analogy with ( [ eq : qcdvectorcurrent ] ) , @xmath504 \psi^\dag\sigma^i\,{\bf d^2}\chi_{|\rm nrqcd } & = &   \tilde{d}_{v2 } \,\psi^\dag\sigma^i\,{\bf d^2}\chi_{|\rm pnrqcd } + \ldots.\end{aligned}\ ] ] non - trivial ( @xmath505 , @xmath506 ) pnrqcd matching coefficients of the currents can arise from three sources : ( 1 ) soft loops not accounted in the matching of the lagrangian .
this implies that the soft loop momentum must flow through the external current vertex .
( 2 ) off - shell effects .
since the pnrqcd lagrangian is matched on - shell , off - shell effects that are not reproduced by the lagrangian interactions must be absorbed into a renormalization of the external currents .
( 3 ) @xmath507 terms of soft loops contributing to the matching of the lagrangian that multiply @xmath42 poles of pnrqcd loops are local and must be absorbed into a renormalization of the external currents , when the matching coefficients ( the potentials ) are defined in four dimensions . as discussed before
, we choose to work with @xmath33-dimensional potentials , hence these contributions are included in the pnrqcd calculation without a modification of the external current .
we shall now prove that there is also no further renormalization of the currents from ( 1 ) and ( 2 ) , that is , @xmath508 to three loops , and @xmath509 , @xmath510 at one loop .
we first consider the issue ( 1 ) of soft renormalization of the nrqcd currents .
the relevant vertex diagram at the one - loop order is shown in figure  [ fig : vertex1loop ] with external momenta @xmath511 , @xmath512 and @xmath513 , and @xmath514 .
the nrqcd expression is @xmath515 where the @xmath70 prescription for the propagators is left implicit .
the unspecified polynomial factor arises from derivatives in the subleading nrqcd interactions or from the subleading external current .
power - counting for soft loop momentum @xmath516 shows that this integral gives rise to a @xmath517 correction to @xmath518 .
however , in the soft region we pick the pole at @xmath519 of the gluon propagator and expand the two quark propagators in @xmath69 and @xmath520 . the resulting integral @xmath521 is scaleless and
vanishes in dimensional regularization .
this holds to any order in the expansion in the soft region  @xcite , hence there is no soft one - loop correction to the matching of the pnrqcd currents .    moving to the two - loop level
, we consider as an example the planar vertex diagram in figure  [ fig : vertexplanar2loop ] .
the momentum regions of interest are s - s and s - p , where the first letter refers to the inner vertex subgraph and loop momentum @xmath522 , the second to the box subgraph and loop momentum @xmath523 .
the p - s combination is not relevant , since the soft loop does not flow through the external vertex .
such contributions are included in the one - loop potentials .
the inner vertex subdiagram in the s - s and s - p regions is an expression similar to ( [ vertex1loopnrqcd ] ) with @xmath524 in the quark propagator and @xmath525 in the gluon propagator . picking up the gluon propagator pole results in @xmath526 for the inner integral
now , in the s - p region @xmath523 is potential , and @xmath527 must be expanded , in which case the integral is scaleless as before .
if , as in the s - s region , @xmath523 is also soft , then the @xmath528-integration picks up the pole of the second gluon propagator at @xmath519 and the entire 2-loop integral is scaleless .
this discussion evidently applies to all 2-loop vertex diagrams . since the vanishing of the integrals is due to the analytic structure of the propagator denominators , it generalizes to higher orders in the velocity expansion , which contains higher powers of the same propagators and numerator polynomials .
hence , we conclude that there are no s - p and s - s contributions to two - loop vertex diagrams in any order in the threshold expansion , in accordance with the results of  @xcite , and hence no soft renormalization of the currents at two loops .
the structure of the analysis extends to higher loop orders .
either one of the inner sub - graphs is scaleless , because an outer loop momentum is potential . or the entire diagram is soft and scaleless , because the external quark momenta are potential .
we therefore conclude that there is no contribution to the matching of the nrqcd currents to pnrqcd from soft loops ( item ( 1 ) ) .
we now turn to the discussion of off - shell effects , item ( 2 ) , and start again with the one - loop diagram shown in figure  [ fig : vertex1loop ] .
now , however , the loop momentum is potential , and we have to compare the nrqcd potential contribution contained in ( [ vertex1loopnrqcd ] ) with the pnrqcd expression @xmath529 which arises from the single insertion of the tree - level pnrqcd potential and the @xmath517 term of the coulomb green function .
the factor @xmath530 corresponds to the coulomb potential insertion , while higher - order potentials as well as derivative factors from the external current are contained in the unspecified polynomial .
the potential contribution in nrqcd is defined as the contribution from the quark - propagator pole at @xmath531 in ( [ vertex1loopnrqcd ] ) .
the difference @xmath532 between @xmath533 and @xmath534 contributes to the matching of the external current and arises as follows : when the tree - level pnrqcd potential is derived from the one - gluon exchange diagram , the external quark lines are assumed on - shell , which implies @xmath535 and @xmath536 with the momentum assignment as in figure  [ fig : vertex1loop ] .
however , no such restrictions are imposed on the loop momentum @xmath523 in the calculation of the nrqcd diagram , figure  [ fig : vertex1loop ] .
thus the difference between @xmath533 and @xmath534 comes from the expansion of the gluon propagator in the potential region @xmath537
^ 2}{{{{\bf { k}}}}^4 } + { \cal o}(v^2 ) , \label{offshell1}\ ] ] and the difference of polynomial factors , which after a short computation can be determined to be @xmath538 note that the leading contribution to @xmath533 and @xmath534 is of order @xmath86 .
the terms neglected in ( [ offshell1 ] ) and ( [ offshell2 ] ) are therefore of order @xmath539 .
these are fourth - order corrections to the cross section beyond the accuracy we aim at . in total
, we arrive at @xmath540 ^ 2}{{{{\bf { k}}}}^2}\right\}_{\!|k^0 = \frac{e}{2}-\frac{({{{\bf { p}}}}+{{{\bf { k}}}})^2}{2 m } } \nonumber\\ & = &    ( i g_s)^2 c_f\tilde\mu^{2\epsilon }   \!\int \frac{d^{d-1 } { { { \bf { k}}}}}{(2\pi)^{d-1}}\ , \frac{1}{{{{\bf { k}}}}^2}\ , \frac{1}{2 k^0 } \times   \left\{\frac{k^0}{m } - \frac{[k^0]^2}{{{{\bf { k}}}}^2}\right\}_{\!|k^0 = \frac{e}{2}-\frac{({{{\bf { p}}}}+{{{\bf { k}}}})^2}{2 m}}\,.\end{aligned}\ ] ] to arrive at the second line we used that the on - shell condition implies @xmath541 .
the two terms in curly brackets each cancel the heavy - quark propagator denominator @xmath542 , after which the integral is scaleless and vanishes .
this must be so , since a non - zero contribution to @xmath532 at this order would have scaled as @xmath467 , but there is no @xmath543 production current .
that off - shell effects from the potential region _ are _ relevant in general can be seen from the calculation of the @xmath133-potential at @xmath544 .
the potential region of the one - loop box diagram shown in figure  [ fig : box1loop ] is not completely reproduced by pnrqcd . the difference is a contribution to the @xmath133 potential , that is , a pnrqcd matching coefficient , which is crucial to obtain the gauge - invariant result given in  @xcite and ( [ b1eps ] ) .
the difference between the box and the vertex diagram discussed above is that the box loop integral is not scaleless after the cancellation of the quark propagator by the off - shell terms , since there is a second gluon propagator @xmath545 .
returning to the vertex diagrams , we now consider the planar two - loop diagram of figure  [ fig : vertexplanar2loop ] in the p - p and p - s region , where the first letter refers again to the inner vertex subgraph .
the corresponding pnrqcd diagrams are the two - loop vertex diagram with tree - level potentials and a one - loop vertex diagram with insertion of a one - loop potential , respectively .
the off - shell terms of the nrqcd diagram in the p - p region are precisely the ones that contribute to the @xmath133 potential discussed in the previous paragraph ; they are correctly reproduced by the pnrqcd diagram with the insertion of the @xmath133 potential .
the off - shell terms in the p - s region have a similar origin as in the one - loop vertex diagram .
the soft box graph gives rise to a one - loop potential , but since the potentials are matched on - shell , the soft box graph is not completely reproduced when it appears as a subgraph in a larger diagram .
since the leading p - s region is @xmath546 , if all the off - shell corrections were of order @xmath118 relative to the leading term as in ( [ offshell1 ] ) , ( [ offshell2 ] ) , we could immediately dismiss them , since there is no @xmath547 hard vertex correction into which it could be absorbed .
this is indeed true for the planar diagram but not in general .    as an example
, we consider the non - planar nrqcd two - loop diagram shown in figure  [ fig : vertexnonplanar2loop ] , which ( neglecting constant factors ) is given by the expression @xmath548 in the p - s region the @xmath522-integral for the soft crossed - box subdiagram is exactly the same that appears in the computation of the one - loop coulomb and @xmath133 potentials except for the additional @xmath528 in the last quark propagator that is absent when the quark lines of the inner vertex subgraph are on - shell ( see discussion above ) . since @xmath523 is potential , @xmath549 must be expanded relative to @xmath550 , which results in a series of corrections beginning at @xmath543 .
if this off - shell correction were non - zero , it would result in a nnlo @xmath544 contribution to the coefficient function @xmath551 .
however , as @xmath523 is potential , the @xmath528-integral is the contribution from the pole of the quark propagator in the first line of ( [ nonplanar2loopnrqcd ] ) .
thus , in complete analogy with the discussion of the one - loop vertex diagram , the expansion in @xmath528 cancels the remaining quark propagator and renders the @xmath552-integral scaleless .
this cancellation is generic for all two - loop vertex diagrams in the p - s region .
the structure of the argument extends to the 3-loop order .
the possible off - shell terms are either already accounted in the matching of subleading potentials ; or an inner vertex subgraph becomes scaleless due to a cancellation of the remaining quark - propagator of a potential loop .
hence , we conclude that there is no contribution to the matching of the nrqcd currents to pnrqcd from off - shell terms ( item ( 2 ) ) , at least up to the nnnlo order .
the integrals for pnrqcd diagrams with insertions of potentials or external currents such as ( [ multipleinsertions ] ) can sometimes be simplified by the equation of motion for the pnrqcd quark - antiquark propagator .
we will make use of this to reduce the number of independent insertions that need to be calculated ( in part ii of the paper ) , and provide the relevant identities here .    in the present context
the equation of motion is the lippmann - schwinger equation ( [ eq : lippmann - schwinger ] ) for the colour - singlet coulomb green function written in the form @xmath553 in the remainder of this subsection we drop the energy argument of the green function , which is always @xmath69 , and set @xmath554 to simplify the notation .    as our first example we consider the insertion of the subleading derivative current ( [ eq : qcdvectorcurrent2 ] ) into one of the vertices .
the relevant integral is @xmath555 the second term on the right - hand side is a scaleless integral .
the third one is seen to contain @xmath556 after shifting @xmath557 .
hence , the insertion of @xmath558 can be replaced by the factor @xmath559 .
this holds true in expressions of the form @xmath310 \,\frac{{{\bf { p}}}^2_1}{m^2}\ ,   i{g_0}({\bf p}_1,{\bf p}_2 ) i\delta{v}_1({\bf p}_2,{\bf p}_3 ) i{g_0}({\bf p}_3,{\bf p}_4)\ ,   i\delta{v}_2({\bf p}_4,{\bf p}_5 )   i{g_0}({\bf p}_5,{\bf p}_6)\ ,   \ldots\ ] ] that contain multiple potential insertions .
this shows that the _
unrenormalized , dimensionally regulated _ matrix element of the subleading external current operator @xmath560 is @xmath561 times the unrenormalized matrix element of the leading - order current @xmath562 , which has been used in ( [ r1 ] ) .
as our second example we consider the insertion of the @xmath563 potential .
with @xmath564 we find @xmath565^{\frac{1}{2}+\epsilon } } \bigg ] g_0({\bf p}_3,{\bf p}_4)\,.\end{aligned}\ ] ] the final expression follows , since the second term on the left - hand side of the first equation contains the scaleless @xmath566-integral , while in the third term the integration over @xmath567 can be performed with the help of the integral @xmath568^{1+u } ( { { \bf { q}}}-{{\bf { k}}})^2 }   \equiv \frac{1}{[{{\bf { q}}}^2]^{\frac{1}{2}+u+\epsilon}}\ , \frac{k(u)}{8}\end{aligned}\ ] ] with @xmath569 whose general form for @xmath570 will be needed below .
this shows that the insertion of the @xmath563 potential can be eliminated in favour of the insertion of the coulomb and the @xmath33-dimensional @xmath571 potential .    in the general case ,
after applying the spin projection , only six different types of insertions are needed for the nnnlo calculation : @xmath572 the first four come from the potentials in ( [ eq : potentialafterspin ] ) , the fifth is the kinetic energy correction and the last might be used for the conversion of the pole scheme to threshold mass scheme as discussed in paper  ii .
the identities given below show that the last three types of insertions can be reduced by using the equation of motion to the first three and the delta - function potential @xmath573 without factors of @xmath274 . at third order , both , single insertions and double insertions with an additional coulomb potential insertion have to be considered . for the single insertions ,
the equation of motion relations read : @xmath574^{\frac{1}{2}+\epsilon } } \frac{k(0)}{8 } g_0({\bf p}_3,{\bf p}_4 ) , \\
\frac{{\bf p}_2 ^ 2+{\bf p}_3 ^ 2}{2m^2 { \bf q}^2}:&&\int \prod_{i=1}^4\frac{d^{d-1}{\bf p}_i}{(2\pi)^{d-1}}g_0({\bf p}_1,{\bf p}_2 ) \frac{{\bf p}_2 ^ 2+{\bf p}_3 ^ 2}{2m^2 { \bf q}^2 } g_0({\bf p}_3,{\bf
p}_4 ) \label{eompoverq } \\ & = & \int \prod_{i=1}^4\frac{d^{d-1}{\bf p}_i}{(2\pi)^{d-1}}g_0({\bf p}_1,{\bf p}_2 ) \bigg[\frac{e}{m } \frac{1}{{\bf q}^2 } + \frac{\pi c_f\alpha_s}{2 m } \frac{k(0)}{[{\bf q}^2]^{\frac{1}{2}+\epsilon } } \bigg ] g_0({\bf p}_3,{\bf p}_4).\nonumber\end{aligned}\ ] ] the delta - function potential @xmath573 appears implicitly in the first two relations in integrands of the form @xmath575 , where the delta - function has been eliminated to set two arguments equal .
the equation of motion identities for the double insertions with a coulomb potential ( with @xmath576 and @xmath577 ) read : @xmath578^{1+u}}\nonumber \\
\nonumber   & = & e\int \prod_{i=1,3 - 6}\frac{d^{d-1}{\bf p}_i}{(2\pi)^{d-1 } } g_0({\bf p}_1,{\bf p}_3 ) \frac{g_0({\bf
p}_3,{\bf p}_4)g_0({\bf p}_5,{\bf p}_6)}{[{\bf q}_2 ^ 2]^{1+u } } \\ \nonumber & + & \int \prod_{i=1}^4\frac{d^{d-1}{\bf p}_i}{(2\pi)^{d-1 } } \frac{g_0({\bf p}_1,{\bf p}_2)g_0({\bf p}_3,{\bf p}_4)}{[{\bf q}_1 ^ 2]^{1+u } } \\   & + & 4\pi c_f \alpha_s\int \prod_{i=1}^6\frac{d^{d-1}{\bf p}_i}{(2\pi)^{d-1}}g_0({\bf p}_1,{\bf p}_2)\frac{1}{{\bf q}_1 ^ 2 } \frac{g_0({\bf p}_3,{\bf p}_4)g_0({\bf p}_5,{\bf p}_6)}{[{\bf q}_2 ^ 2]^{1+u } } , \\
\delta^{(d-1)}({\bf q}_1)\frac{{\bf p}_2 ^ 4}{m^3}:&&\int \prod_{i=1}^6\frac{d^{d-1}{\bf p}_i}{(2\pi)^{d-1}}g_0({\bf p}_1,{\bf p}_2 ) ( 2\pi)^{d-1}\delta^{(d-1)}({\bf q}_1)\frac{{\bf p}_2 ^ 4}{m^3 } \frac{g_0({\bf p}_3,{\bf p}_4)g_0({\bf p}_5,{\bf p}_6)}{[{\bf q}_2 ^ 2]^{1+u}}\nonumber \\
\nonumber   & = & \frac{e^2}{m}\int \prod_{i=1,3 - 6}\frac{d^{d-1}{\bf p}_i}{(2\pi)^{d-1 } } g_0({\bf p}_1,{\bf p}_3 ) \frac{g_0({\bf
p}_3,{\bf p}_4)g_0({\bf p}_5,{\bf p}_6)}{[{\bf q}_2 ^ 2]^{1+u } } \\
\nonumber&+&\int \prod_{i=1}^4\frac{d^{d-1}{\bf p}_i}{(2\pi)^{d-1}}g_0({\bf p}_1,{\bf p}_2 ) \left[\frac{2e}{m}\frac{1}{[{\bf q}_1 ^
2]^{1+u}}+   \frac{\pi c_f \alpha_s}{2 m}\frac{k(u)}{[{\bf q}_1 ^ 2]^{\frac{1}{2}+u+\epsilon}}\right ] g_0({\bf p}_3,{\bf p}_4 ) \\ \nonumber&+&4\pi c_f \alpha_s\int \prod_{i=1}^6\frac{d^{d-1}{\bf p}_i}{(2\pi)^{d-1}}g_0({\bf p}_1,{\bf p}_2 ) \left[\frac{2e}{m}\frac{1}{{\bf q}_1 ^ 2}+ \frac{\pi c_f \alpha_s}{2m}\frac{k(0)}{[{\bf q}_1 ^ 2]^{\frac{1}{2}+\epsilon}}\right ] \\ & & \times\frac{g_0({\bf p}_3,{\bf p}_4)g_0({\bf p}_5,{\bf p}_6)}{[{\bf q}_2 ^ 2]^{1+u } } ,
\\ \frac{{\bf p}_2 ^ 2+{\bf p}_3 ^ 2}{2 m^2 { \bf q}_1 ^ 2}:&&\int \prod_{i=1}^6\frac{d^{d-1}{\bf p}_i}{(2\pi)^{d-1}}g_0({\bf p}_1,{\bf p}_2 ) \frac{{\bf
p}_2 ^ 2+{\bf p}_3 ^ 2}{2 m^2 { \bf q}_1 ^ 2 } \frac{g_0({\bf p}_3,{\bf p}_4)g_0({\bf p}_5,{\bf p}_6)}{[{\bf q}_2 ^ 2]^{1+u}}\nonumber \\ & & \hspace*{-1 cm } = \,\int \prod_{i=1}^6\frac{d^{d-1}{\bf p}_i}{(2\pi)^{d-1}}g_0({\bf p}_1,{\bf p}_2 ) \left[\frac{e}{m } \frac{1}{{\bf q}_1 ^ 2 } + \frac{\pi c_f\alpha_s}{2 m } \frac{k(0)}{[{\bf q}_1 ^ 2]^{\frac{1}{2}+\epsilon}}\right ] \frac{g_0({\bf p}_3,{\bf p}_4)g_0({\bf p}_5,{\bf p}_6)}{[{\bf q}_2 ^ 2]^{1+u}}\nonumber \\    & & \hspace*{-1 cm }   + \ , \frac{k(u)}{16m}\int\prod_{i=1}^4\frac{d^{d-1}{\bf p}}{(2\pi)^{d-1 } }   \frac{g_0({\bf p}_1,{\bf p}_2)g_0({\bf p}_3,{\bf p}_4)}{[{\bf q}_1 ^ 2]^{\frac{1}{2}+u+\epsilon } } \,.\end{aligned}\ ] ] we note that we have kept the insertion of the coulomb potential in the more general form @xmath579^{1+u}$ ] , since in the double insertions we may also need the @xmath33-dimensional coulomb potential , implying @xmath580 .      at third order there is for the first time a contribution from an ultrasoft loop momentum region .
the ultrasoft correction to the heavy - quark correlation function has been computed separately in @xcite and will be incorporated in the results shown in paper  ii . here
we provide a short overview of the ultrasoft calculation and discuss some issues of the factorization , which are important to understand the splitting of various divergent parts .
the relevant ultrasoft interaction terms in the pnrqcd lagrangian ( [ eq : pnrqcd ] ) are given by @xmath581\psi(x ) +   g_s\chi^\dagger(x)\big[a_0(t,{{\bf { 0 } } } ) -{{\bf { x}}}\cdot { { \bf { e}}}(t,{{\bf { 0}}})\big]\chi(x ) .
\label{usterms}\ ] ] the derivation of the chromoelectric dipole interaction from the multipole expansion of the nrqcd lagrangian can be found in @xcite .
the interaction with @xmath582 can be removed by a field redefinition involving a time - like wilson line .
this modifies the external current that creates the heavy - quark pair , as discussed in @xcite . in the present case of colour - singlet production in @xmath0 collisions
the wilson lines cancel , and the @xmath582 terms in ( [ usterms ] ) can be dropped . with @xmath583 , and @xmath584 for ultrasoft gluon fields , it follows that the chromoelectric dipole interaction is suppressed by @xmath167 relative to the kinetic term in the action .
two ultrasoft interaction vertices are required to form a loop , from which it follows that the leading ultrasoft contribution arises first at the third order .
the ultrasoft correction can be expressed in the form @xmath585,\ ,   \label{eq : gfcorr}\end{aligned}\ ] ] with the understanding that one picks up only the pole at @xmath586 in the gluon propagator . here
@xmath587 is the colour - singlet and @xmath588 the colour - octet coulomb green function ( [ eq : lippmann - schwinger ] ) .
however , as explained in  @xcite , this expression can not be used in practice , because it is ultraviolet ( uv ) divergent .
the regularization and subtraction of divergences must be done consistently with the calculation of potential insertions and hard matching coefficients , which have been done in dimensional regularization . in order to apply dimensional regularization to the ultrasoft contribution , ( [ eq : gfcorr ] ) is transformed to momentum space .
it also turns out to be convenient to revert the derivation of the pnrqcd ultrasoft interaction ( [ usterms ] ) and to instead use the nrqcd vertices .
the reason for this is that the derivation of ( [ usterms ] ) uses the pnrqcd equation of motion , which reshuffles the loop expansion , and employs four - dimensional identities @xcite .
the correspondence between uv divergences in the ultrasoft calculation and ir divergences in the potential and hard matching calculations is more directly seen at the level of nrqcd diagrams , and the correct evaluation of the finite terms requires the consistent use of dimensional regularization in every loop order .
the uv divergences arise from the integral over the three - momentum @xmath589 of the ultrasoft gluon , and from the subsequent potential loop integrations .
the former divergence is related to the factorization of the ultrasoft scale from the other scales , and cancels when all pieces of the calculation are added .
the uv - divergent part of the ultrasoft integral has the form of a single insertion of a third - order potential and of a one - loop correction to the coefficient @xmath18 of the derivative current in ( [ eq : dv1loop ] ) .
we therefore define the ultrasoft correction by adding counterterms that cancel these ultrasoft subdivergences . with these subtractions ,
the ultrasoft correction reads  @xcite @xmath590 ^ 2\int \frac{d^{d-1}{\mbox{\boldmath $ { \ell}$}}}{(2\pi)^{d-1 } } \frac{d^{d-1}{\mbox{\boldmath $ { \ell ^\prime}$}}}{(2\pi)^{d-1 } } \bigg\{\delta d_v^{\rm div}\,(-1 ) \frac{{\mbox{\boldmath $ { \ell}$}}^2+{\mbox{\boldmath $ { \ell}$}}^{\prime\,2 } } { 6 m^2}\ , g^{(1)}_0({\mbox{\boldmath $ { \ell}$}},{\mbox{\boldmath $ { \ell}$}}^\prime;e ) \nonumber\\[0.2 cm ] & & \hspace*{-1.5 cm }   + \,\big[\tilde\mu^{2\epsilon}\big]^2 \int \frac{d^{d-1}{\bf p}}{(2\pi)^{d-1 } } \frac{d^{d-1}{\bf p^\prime}}{(2\pi)^{d-1 } } \,g^{(1)}_0({\mbox{\boldmath $ { \ell}$}},{\bf p};e ) \,i\big[\delta u -\delta   v_{c.t.}\big ] \,ig^{(1)}_0({\bf p^\prime},{\mbox{\boldmath $ { \ell^\prime}$}};e ) \bigg\}.\qquad \label{eq : defus}\end{aligned}\ ] ] here @xmath472 represents the potential subtraction ( [ eq : counter ] ) , and @xmath591 is the ultrasoft insertion ( containing the octet green function ) . with an opposite sign compared to  @xcite . ]
the first line of ( [ eq : defus ] ) is related to the renormalization of the @xmath129 suppressed vector current @xmath592 . if @xmath593_{\rm ren } = z_{1/m^2 }
[ j_{1/m^2}]_{\rm bare}$ ] , the one - loop counterterm @xmath594 equals the infrared divergence @xmath595 that was subtracted to obtain the finite expression ( [ eq : dv1loop ] ) .
the explicit expression is @xmath596 the remaining divergences are associated with the three - loop hard matching coefficient @xmath25 .
note that this divergence structure implies that the nrqcd current mixes with the @xmath129 suppressed current through ultrasoft interactions .
schematically , @xmath597^{2}\,\right)_{| \rm us }   = \mbox{const}\times\frac{\alpha_s}{\epsilon } \ , \frac{1}{6 m^2}\,\psi^\dagger \sigma^i { { \bf { d^2}}}\chi + \ldots \ , , \label{eq : usdiv}\ ] ] where the ellipses denote the remaining divergences related to time - ordered products with potentials and the leading - power current itself .
the subtracted expression is then simplified and reduced to a number of integrations that can mostly be done only numerically .
the code that computes the ultrasoft correction was developed in conjunction with @xcite , and is implemented in our program for the third - order cross section .
we have now collected all prerequisites to write down the expansion of the non - relativistic correlation function @xmath598(x)\ , [ \psi^{\dag}\sigma^i\chi](0 ) )    non - relativistic ( pnrqcd ) perturbation theory . adopting the operator notation from section  [ sec : qm ] ,
the expansion is given by @xmath599 with @xmath600 as given in ( [ eq : g00msbar ] ) , and @xmath601 \delta_2 g(e ) & = & \langle { { \bf { 0}}}| \hat g_0(e)i \delta v_1   i\hat g_0(e ) i\delta v_1 i\hat g_0(e)|{{\bf { 0}}}\rangle +    \langle { { \bf { 0}}}| \hat g_0(e)i\delta v_2   i\hat g_0(e)|{{\bf { 0}}}\rangle,\\[0.2 cm ] \delta_3 g(e ) & = & \langle { { \bf { 0}}}| \hat g_0(e)i \delta v_1   i\hat g_0(e ) i \delta v_1   i\hat g_0(e ) i\delta v_1 i\hat g_0(e)|{{\bf { 0}}}\rangle   \nonumber\\ & & + \ ,   2 \langle { { \bf { 0}}}| \hat g_0(e)i \delta v_1   i\hat g_0(e ) i\delta v_2 i\hat g_0(e)|{{\bf { 0}}}\rangle +    \langle { { \bf { 0}}}| \hat g_0(e)i\delta v_3   i\hat g_0(e)|{{\bf { 0}}}\rangle \nonumber\\   & & + \ , \delta^{us}g(e)\ ,
. \label{masterthirdorder}\end{aligned}\ ] ] in momentum space these expressions are of the from of single and multiple insertions of potentials as defined in ( [ multipleinsertions ] ) plus the ultrasoft correction .
@xmath602 denotes a potential correction of order @xmath250 .
the first - order potential consists only of the one - loop correction to the coulomb potential : @xmath603 at second order , we have the two - loop coulomb potential , the one - loop @xmath604 potential and the @xmath118-suppressed potentials at tree - level . together with the kinetic - energy correction , we obtain @xmath605   \nonumber\\[-0.1 cm ]   & & - \,2\times \frac{{{\bf { p}}}^4}{8 m^3}\,(2\pi)^{d-1}\delta^{(d-1)}({{\bf { q}}}).\end{aligned}\ ] ] the notation for the potentials is defined in ( [ eq : potentialafterspin ] ) , where also the explicit expressions are given .
note that different from ( [ eq : calvexp ] ) the @xmath250-loop potential @xmath606 includes the coupling constant factor , i.e. @xmath607 . there are no new potentials appearing at third order , hence @xmath608.\end{aligned}\ ] ] note that there is no kinetic - energy term in @xmath609 , because the kinetic energy term in the lagrangian is not renormalized .
these potentials appear as single insertions in @xmath610 .
in addition , @xmath610 receives contributions from multiple insertions of the lower - order potentials . to complete the perturbative expansion of the third - order cross section
, the expansion ( [ greenexpand ] ) of the green function is inserted into ( [ r1 ] ) , @xmath611 g(e ) + \ldots\right ) \right].\ ] ] the coefficient functions @xmath17 and @xmath18 are likewise expanded , and product terms of order higher than three in non - relativistic perturbation theory are dropped . note that @xmath612 counts as second order in this expansion .
this work concerns the third - order correction @xmath613 .
the triple insertion of the first - order coulomb potential is algebraically complicated but has no uv or ir divergences and can therefore be computed numerically as done in  @xcite .
the ultrasoft correction was obtained in  @xcite . in part
ii of the present work we will give the details of the calculation of the remaining single and double insertions in the second line of ( [ masterthirdorder ] ) in dimensional regularization as is necessary for a consistent combination with the matching calculations performed and summarized in part i. a more precise master formula for the third - order correction to the green function that accounts for the pole structure of the @xmath33-dimensional potentials will be given in part ii .
we thank a.  maier , j.  piclum and m.  steinhauser for comments on the text .
this work has been supported by the dfg sonderforschungsbereich / transregio  9 `` computergesttzte theoretische teilchenphysik '' , the dfg graduiertenkolleg `` elementarteilchenphysik an der tev - skala '' and the dfg cluster of excellence `` origin and structure of the universe '' .
i.  i.  y. bigi , y.  l. dokshitzer , v.  a. khoze , j.  h. khn and p.  m. zerwas , _ production and decay properties of ultraheavy quarks _ , _ phys .
_ * b181 * ( 1986 ) 157 . and * d0 * collaboration , t.  aaltonen _ et  al .
_ , _ combination of the top - quark mass measurements from the tevatron collider _ , _ phys.rev . _
* d86 * ( 2012 ) 092003 [ http://arxiv.org/abs/1207.1069[arxiv:1207.1069 ] ] . and
* cms * collaboration , _ combination of atlas and cms results on the mass of the top quark using up to 4.9 fb@xmath614 of data _ , atlas - conf-2012 - 095 , cms - pas - top-12 - 001 .
m.  beneke , p.  falgari , s.  klein , j.  piclum , c.  schwinn _ et  al .
_ , _ inclusive top - pair production phenomenology with topixs _ , _ jhep _ * 1207 * ( 2012 ) 194 [ http://arxiv.org/abs/1206.2454[arxiv:1206.2454 ] ] . m.  martinez and r.  miquel , _ multi - parameter fits to the t anti - t threshold observables at a future @xmath0 linear collider _ , _ eur .
phys . j. _ * c27 * ( 2003 ) 4955 [ http://arxiv.org/abs/hep-ph/0207315 [ hep - ph/0207315 ] ] .
k.  seidel , f.  simon , m.  tesar and s.  poss , _ top quark mass measurements at and above threshold at clic _ , _ eur.phys.j . _
* c73 * ( 2013 ) 2530 [ http://arxiv.org/abs/1303.3758[arxiv:1303.3758 ] ] .
t.  horiguchi , a.  ishikawa , t.  suehara , k.  fujii , y.  sumino _ et  al .
_ , _ study of top quark pair production near threshold at the ilc _ , [ http://arxiv.org/abs/1310.0563[arxiv:1310.0563 ] ] .
collaboration , g.  aad _
_ , _ observation of a new particle in the search for the standard model higgs boson with the atlas detector at the lhc _ , _ phys.lett . _ * b716 * ( 2012 ) 129 [ http://arxiv.org/abs/1207.7214[arxiv:1207.7214 ] ] .
collaboration , s.  chatrchyan _ et  al .
_ , _ observation of a new boson at a mass of 125 gev with the cms experiment at the lhc _ , _ phys.lett . _ * b716 * ( 2012 ) 3061 [ http://arxiv.org/abs/1207.7235[arxiv:1207.7235 ] ] .
g.  degrassi , s.  di  vita , j.  elias - miro , j.  r. espinosa , g.  f. giudice _
_ , _ higgs mass and vacuum stability in the standard model at nnlo _ , _ jhep _ * 1208 * ( 2012 ) 098 [ http://arxiv.org/abs/1205.6497[arxiv:1205.6497 ] ] . b.  a. thacker and g.  p. lepage , _ heavy quark bound states in lattice qcd _ , _ phys .
* d43 * ( 1991 ) 196208 .
g.  p. lepage , l.  magnea , c.  nakhleh , u.  magnea and k.  hornbostel , _ improved nonrelativistic qcd for heavy quark physics _ , _ phys .
* d46 * ( 1992 ) 40524067 [ http://arxiv.org/abs/hep-lat/9205007 [ hep - lat/9205007 ] ] .
g.  t. bodwin , e.  braaten and g.  p. lepage , _ rigorous qcd analysis of inclusive annihilation and production of heavy quarkonium _ , _ phys .
* d51 * ( 1995 ) 11251171 [ http://arxiv.org/abs/hep-ph/9407339 [ hep - ph/9407339 ] ] .
m.  beneke and v.  a. smirnov , _
asymptotic expansion of feynman integrals near threshold _ , _ nucl .
phys . _ * b522 * ( 1998 ) 321344 [ http://arxiv.org/abs/hep-ph/9711391[hep-ph/9711391 ] ]
. v.  s. fadin and v.  a. khoze , _ threshold behavior of heavy top production in @xmath0 collisions _ , _ jetp lett .
_ * 46 * ( 1987 ) 525529 .
v.  s. fadin and v.  a. khoze , _ production of a pair of heavy quarks in @xmath615 annihilation in the threshold region _ , _ sov .
_ * 48 * ( 1988 ) 309313 .
m.  j. strassler and m.  e. peskin , _ the heavy top quark threshold : qcd and the higgs _ , _ phys .
* d43 * ( 1991 ) 15001514 .
m.  jezabek , j.  h. khn and t.  teubner , _ momentum distributions in t anti - t production and decay near threshold _ , _ z. phys . _
* c56 * ( 1992 ) 653660 .
y.  sumino , k.  fujii , k.  hagiwara , h.  murayama and c.  k. ng , _ top quark pair production near threshold _ , _ phys .
* d47 * ( 1993 ) 5681 .
k.  fujii , t.  matsui and y.  sumino , _ physics at @xmath9 threshold in @xmath0 collisions _ , _ phys .
_ * d50 * ( 1994 ) 43414362 .
r.  harlander , m.  jezabek , j.  h. khn and t.  teubner , _ polarization in top quark pair production near threshold _ , _ phys .
* b346 * ( 1995 ) 137142 [ http://arxiv.org/abs/hep-ph/9411395[hep-ph/9411395 ] ] .
a.  h. hoang and t.  teubner , _ top quark pair production at threshold : complete next - to - next - to - leading order relativistic corrections _
* d58 * ( 1998 ) 114023 [ http://arxiv.org/abs/hep-ph/9801397[hep-ph/9801397 ] ] .
k.  melnikov and a.  yelkhovsky , _ top quark production at threshold with o(@xmath616 ) accuracy _ , _ nucl .
* b528 * ( 1998 ) 5972 [ http://arxiv.org/abs/hep-ph/9802379[hep-ph/9802379 ] ] .
m.  beneke , a.  signer and v.  a. smirnov , _ top quark production near threshold and the top quark mass _ , _ phys .
* b454 * ( 1999 ) 137146 [ http://arxiv.org/abs/hep-ph/9903260[hep-ph/9903260 ] ] .
a.  h. hoang and t.  teubner , _ top quark pair production close to threshold : top mass , width and momentum distribution _ , _ phys .
* d60 * ( 1999 ) 114027 [ http://arxiv.org/abs/hep-ph/9904468 [ hep - ph/9904468 ] ] .
o.  i. yakovlev , _ top quark production near threshold : nnlo qcd correction _
lett . _ * b457 * ( 1999 ) 170176 [ http://arxiv.org/abs/hep-ph/9808463[hep-ph/9808463 ] ] .
t.  nagano , a.  ota and y.  sumino , _
@xmath617 corrections to @xmath618 total and differential cross sections near threshold _ , _ phys .
rev . _ * d60 * ( 1999 ) 114014 [ http://arxiv.org/abs/hep-ph/9903498[hep-ph/9903498 ] ] .
a.  a. penin and a.  a. pivovarov , _ analytical results for @xmath618 and @xmath619 observables near the threshold up to the next- to - next - to - leading order of nrqcd _ , _ phys . atom .
* 64 * ( 2001 ) 275293 [ http://arxiv.org/abs/hep-ph/9904278[hep-ph/9904278 ] ] .
m.  beneke , _ a quark mass definition adequate for threshold problems _ , _ phys .
lett . _ * b434 * ( 1998 ) 115125 [ http://arxiv.org/abs/hep-ph/9804241[hep-ph/9804241 ] ] .
a.  h. hoang , a.  v. manohar , i.  w. stewart and t.  teubner , _ a renormalization group improved calculation of top quark production near threshold _
* 86 * ( 2001 ) 19511954 [ http://arxiv.org/abs/hep-ph/0011254[hep-ph/0011254 ] ] .
a.  h. hoang , a.  v. manohar , i.  w. stewart and t.  teubner , _ the threshold @xmath9 cross section at nnll order _ , _ phys .
_ * d65 * ( 2002 ) 014014 [ http://arxiv.org/abs/hep-ph/0107144[hep-ph/0107144 ] ] .
a.  pineda and a.  signer , _ heavy quark pair production near threshold with potential non - relativistic qcd _
* b762 * ( 2007 ) 6794 [ http://arxiv.org/abs/hep-ph/0607239[hep-ph/0607239 ] ] . a.  hoang and m.  stahlhofen , _ the top - antitop threshold at the ilc : nnll qcd uncertainties _ , [ http://arxiv.org/abs/1309.6323[arxiv:1309.6323 ] ] .
n.  brambilla , a.  pineda , j.  soto and a.  vairo , _ the infrared behaviour of the static potential in perturbative qcd _ , _ phys .
rev . _ * d60 * ( 1999 ) 091502 [ http://arxiv.org/abs/hep-ph/9903355 [ hep - ph/9903355 ] ] .
b.  a. kniehl and a.  a. penin , _ ultrasoft effects in heavy quarkonium physics _ , _ nucl .
_ * b563 * ( 1999 ) 200210 [ http://arxiv.org/abs/hep-ph/9907489[hep-ph/9907489 ] ] .
n.  brambilla , a.  pineda , j.  soto and a.  vairo , _ the heavy quarkonium spectrum at order @xmath620 _ , _ phys .
_ * b470 * ( 1999 ) 215 [ http://arxiv.org/abs/hep-ph/9910238 [ hep - ph/9910238 ] ] . b.  a. kniehl and a.  a. penin , _
order @xmath621 corrections to heavy- quarkonium creation and annihilation _ , _ nucl .
phys . _ * b577 * ( 2000 ) 197208 [ http://arxiv.org/abs/hep-ph/9911414[hep-ph/9911414 ] ] .
a.  v. manohar and i.  w. stewart , _ running of the heavy quark production current and @xmath105 potential in qcd _
* d63 * ( 2001 ) 054004 [ http://arxiv.org/abs/hep-ph/0003107[hep-ph/0003107 ] ] .
y.  kiyo and y.  sumino , _
@xmath622 quarkonium 1s spectrum in large @xmath623 approximation and renormalon cancellation _ , _ phys
. lett . _ * b496 * ( 2000 ) 8388 [ http://arxiv.org/abs/hep-ph/0007251[hep-ph/0007251 ] ] . b.  a. kniehl , a.  a. penin , m.  steinhauser and v.  a. smirnov , _
nonabelian @xmath624 heavy - quark - antiquark potential _ , _ phys .
* d65 * ( 2002 ) 091503 [ http://arxiv.org/abs/hep-ph/0106135 [ hep - ph/0106135 ] ] .
b.  a. kniehl , a.  a. penin , v.  a. smirnov and m.  steinhauser , _ potential nrqcd and heavy - quarkonium spectrum at next - to - next - to - next - to - leading order _ , _ nucl .
phys . _ * b635 * ( 2002 ) 357383 [ http://arxiv.org/abs/hep-ph/0203166[hep-ph/0203166 ] ] .
a.  a. penin and m.  steinhauser , _ heavy quarkonium spectrum at @xmath625 and bottom / top quark mass determination _ , _ phys .
_ * b538 * ( 2002 ) 335345 [ http://arxiv.org/abs/hep-ph/0204290[hep-ph/0204290 ] ] .
b.  a. kniehl , a.  a. penin , m.  steinhauser and v.  a. smirnov , _ heavy - quarkonium creation and annihilation with @xmath626 accuracy .
_ , _ phys .
* 90 * ( 2003 ) 212001 [ http://arxiv.org/abs/hep-ph/0210161[hep-ph/0210161 ] ] .
a.  h. hoang , _ three - loop anomalous dimension of the heavy quark pair production current in non - relativistic qcd _ , _ phys .
* d69 * ( 2004 ) 034009 [ http://arxiv.org/abs/hep-ph/0307376 [ hep - ph/0307376 ] ] . m.  beneke , y.  kiyo and k.  schuller , _ third - order coulomb corrections to the s - wave green function , energy levels and wave functions at the origin _ , _ nucl .
_ * b714 * ( 2005 ) 6790 [ http://arxiv.org/abs/hep-ph/0501289[hep-ph/0501289 ] ] .
a.  a. penin , v.  a. smirnov and m.  steinhauser , _ heavy quarkonium spectrum and production / annihilation rates to order @xmath627 _ , _ nucl .
phys . _ * b716 * ( 2005 ) 303318 [ http://arxiv.org/abs/hep-ph/0501042[hep-ph/0501042 ] ] .
p.  marquard , j.  h. piclum , d.  seidel and m.  steinhauser , _ fermionic corrections to the three - loop matching coefficient of the vector current _ , _ nucl .
* b758 * ( 2006 ) 144160 [ http://arxiv.org/abs/hep-ph/0607168[hep-ph/0607168 ] ] .
m.  beneke , y.  kiyo and k.  schuller , _ third - order non - coulomb correction to the s - wave quarkonium wave functions at the origin _ , _ phys .
_ * b658 * ( 2008 ) 222229 [ http://arxiv.org/abs/arxiv:0705.4518    [ hep - ph][arxiv:0705.4518 ] ] .
m.  beneke , y.  kiyo and a.  a. penin , _ ultrasoft contribution to quarkonium production and annihilation _ , _ phys .
lett . _ * b653 * ( 2007 ) 5359 [ http://arxiv.org/abs/arxiv:0706.2733 [ hep - ph][arxiv:0706.2733 ] ] . m.  beneke , y.  kiyo , a.  penin and k.  schuller , _ nnnlo correction to the toponium and bottomonium wave - functions at the origin _ ,
econf * c0705302 * ( 2007 ) top01 [ http://arxiv.org/abs/arxiv:0710.4236 [ hep - ph][arxiv:0710.4236 ] ] .
m.  beneke , y.  kiyo and k.  schuller , _ nnnlo results on top - quark pair production near threshold _ , pos radcor * 2007 * ( 2007 ) 051 [ http://arxiv.org/abs/arxiv:0801.3464    [ hep - ph][arxiv:0801.3464 ] ] .
m.  beneke and y.  kiyo , _ ultrasoft contribution to heavy - quark pair production near threshold _ , _ phys .
* b668 * ( 2008 ) 143147 [ http://arxiv.org/abs/0804.4004 [ arxiv:0804.4004 ] ] .
a.  v. smirnov , v.  a. smirnov and m.  steinhauser , _ fermionic contributions to the three - loop static potential _ , _ phys .
lett . _ * b668 * ( 2008 ) 293298 [ http://arxiv.org/abs/0809.1927[arxiv:0809.1927 ] ] .
p.  marquard , j.  piclum , d.  seidel and m.  steinhauser , _ completely automated computation of the heavy - fermion corrections to the three - loop matching coefficient of the vector current _ , _ phys.lett .
_ * b678 * ( 2009 ) 269275 [ http://arxiv.org/abs/arxiv:0904.0920 [ arxiv:0904.0920 ] ] .    c.  anzai , y.  kiyo and y.  sumino , _
static qcd potential at three - loop order _ , _ phys .
* 104 * ( 2010 ) 112003 [ http://arxiv.org/abs/0911.4335[arxiv:0911.4335 ] ] .
a.  v. smirnov , v.  a. smirnov and m.  steinhauser , _ three - loop static potential _ , _ phys .
* 104 * ( 2010 ) 112002 [ http://arxiv.org/abs/0911.4742[arxiv:0911.4742 ] ] .
y.  kiyo and y.  sumino , _ perturbative heavy quarkonium spectrum at next - to - next - to - next - to - leading order _ , [ arxiv:1309.6571 ] .
a.  v. manohar , _ the hqet / nrqcd lagrangian to order @xmath628 _ , _ phys .
* d56 * ( 1997 ) 230237 [ http://arxiv.org/abs/hep-ph/9701294[hep-ph/9701294 ] ] .
s.  wster , _ heavy quark potential at order @xmath629 [ in german ] _ , _ diploma thesis _
rwth aachen university ( 2003 ) .
m.  e. luke and m.  j. savage , _
power counting in dimensionally regularized nrqcd _ , _ phys .
_ * d57 * ( 1998 ) 413423 [ http://arxiv.org/abs/hep-ph/9707313[hep-ph/9707313 ] ] .
a.  a. penin and a.  a. pivovarov , _ top quark threshold production in gamma gamma collision in the next - to - leading order _ , _ nucl .
* b550 * ( 1999 ) 375396 [ http://arxiv.org/abs/hep-ph/9810496 [ hep - ph/9810496 ] ] .
j.  h. khn and t.  teubner , _ axial contributions at the top threshold _ , _ eur .
phys . j. _
* c9 * ( 1999 ) 221228 [ http://arxiv.org/abs/hep-ph/9903322[hep-ph/9903322 ] ] .
a.  h. hoang , c.  j. reisser and p.  ruiz - femenia , _ phase space matching and finite lifetime effects for top- pair production close to threshold _
* d82 * ( 2010 ) 014005 [ http://arxiv.org/abs/1002.3223[arxiv:1002.3223 ] ] .
m.  beneke , b.  jantzen and p.  ruiz - femenia , _ electroweak non - resonant nlo corrections to @xmath630 in the @xmath631 resonance region _ , _ nucl .
* b840 * ( 2010 ) 186213 [ http://arxiv.org/abs/1004.2188[arxiv:1004.2188 ] ] .
r.  j. guth and j.  h. khn , _ top quark threshold and radiative corrections _
phys . _ * b368 * ( 1992 ) 3856 .
a.  h. hoang and c.  j. reisser , _ electroweak absorptive parts in nrqcd matching conditions _ , _ phys .
* d71 * ( 2005 ) 074022 [ http://arxiv.org/abs/hep-ph/0412258[hep-ph/0412258 ] ] .
a.  h. hoang and c.  j. reisser , _ on electroweak matching conditions for top pair production at threshold _ , _ phys .
* d74 * ( 2006 ) 034002 [ http://arxiv.org/abs/hep-ph/0604104[hep-ph/0604104 ] ] .
d.  eiras and m.  steinhauser , _ complete higgs mass dependence of top quark pair threshold production to order @xmath632 _ , _ nucl .
phys . _ * b757 * ( 2006 ) 197210 [ http://arxiv.org/abs/hep-ph/0605227 [ hep - ph/0605227 ] ] .
y.  kiyo , d.  seidel and m.  steinhauser , _
@xmath633 corrections to the @xmath634 vertex at the top quark threshold _ , http://arxiv.org/abs/0810.1597[arxiv:0810.1597 ] .
s.  actis , m.  beneke , p.  falgari and c.  schwinn , _ dominant nnlo corrections to four - fermion production near the w - pair production threshold _ , _ nucl .
* b807 * ( 2009 ) 132 [ http://arxiv.org/abs/0807.0102[arxiv:0807.0102 ] ] .
m.  beneke , a.  p. chapovsky , a.  signer and g.  zanderighi , _ effective theory approach to unstable particle production _ , _ phys .
* 93 * ( 2004 ) 011602 [ http://arxiv.org/abs/hep-ph/0312331 [ hep - ph/0312331 ] ] .
m.  beneke , a.  p. chapovsky , a.  signer and g.  zanderighi , _ effective theory calculation of resonant high - energy scattering _ , _ nucl .
* b686 * ( 2004 ) 205247 [ http://arxiv.org/abs/hep-ph/0401002 [ hep - ph/0401002 ] ] .
m.  beneke , p.  falgari , c.  schwinn , a.  signer and g.  zanderighi , _ four - fermion production near the w pair production threshold _ , _ nucl .
* b792 * ( 2008 ) 89135 [ http://arxiv.org/abs/0707.0773 [ arxiv:0707.0773 ] ] .
a.  a. penin and j.  h. piclum , _ threshold production of unstable top _ , _ jhep _ * 1201 * ( 2012 ) 034 [ http://arxiv.org/abs/1110.1970 [ arxiv:1110.1970 ] ] .
a.  h. hoang , c.  j. reisser and p.  ruiz - femenia , _ implementing invariant mass cuts and finite lifetime effects in top - antitop production at threshold _ , _ nucl .
* 186 * ( 2009 ) 403406 [ http://arxiv.org/abs/0810.2934[arxiv:0810.2934 ] ] .
b.  jantzen and p.  ruiz - femenia , _ nnlo non - resonant corrections to threshold top - pair production from @xmath0 : endpoint - singular terms _ , _ phys .
_ * d88 * ( 2013 ) 054011 [ http://arxiv.org/abs/1307.4337[arxiv:1307.4337 ] ] .
a.  pineda , _ review of heavy quarkonium at weak coupling _ , _
prog.part.nucl.phys . _
* 67 * ( 2012 ) 735785 [ http://arxiv.org/abs/1111.0165[arxiv:1111.0165 ] ] .
m.  beneke , j.  piclum and t.  rauh , _ in preparation_.    m.  beneke , _ new results on heavy quarks near threshold _
, in : proceedings of the 33rd rencontres de moriond : electroweak interactions and unified theories , les arcs , france , 14 - 21 mar 1998 , [ http://arxiv.org/abs/hep-ph/9806429[hep-ph/9806429 ] ] . w.  a. bardeen , a.  buras , d.  duke and t.  muta , _ deep inelastic scattering beyond the leading order in asymptotically free gauge theories _ ,
_ phys.rev . _ * d18 * ( 1978 ) 3998 . m.  beneke , _ nonrelativistic effective theory for quarkonium production in hadron collisions _ , in : ` the strong interaction , from hadrons to partons ' , stanford 1996 , p.  549574
[ http://arxiv.org/abs/hep-ph/9703429 [ hep - ph/9703429 ] ] .
p.  labelle , _ effective field theories for qed bound states : extending nonrelativistic qed to study retardation effects _ , _ phys . rev .
_ * d58 * ( 1998 ) 093013 [ http://arxiv.org/abs/hep-ph/9608491 [ hep - ph/9608491 ] ] .
m.  e. luke and a.  v. manohar , _ bound states and power counting in effective field theories _ , _ phys .
* d55 * ( 1997 ) 41294140 [ http://arxiv.org/abs/hep-ph/9610534[hep-ph/9610534 ] ] . b.  grinstein and i.  z. rothstein , _ effective field theory and matching in non - relativistic gauge theories _ ,
* d57 * ( 1998 ) 7882 [ http://arxiv.org/abs/hep-ph/9703298[hep-ph/9703298 ] ] .
a.  pineda and j.  soto , _ effective field theory for ultrasoft momenta in nrqcd and nrqed _ , _ nucl .
_ * 64 * ( 1998 ) 428432 [ http://arxiv.org/abs/hep-ph/9707481[hep-ph/9707481 ] ] .
a.  pineda and j.  soto , _ the lamb shift in dimensional regularization _
lett . _ * b420 * ( 1998 ) 391396 [ http://arxiv.org/abs/hep-ph/9711292[hep-ph/9711292 ] ] . n.  brambilla , a.  pineda , j.  soto and a.  vairo , _ potential nrqcd : an effective theory for heavy quarkonium _ , _ nucl .
phys . _ * b566 * ( 2000 ) 275 [ http://arxiv.org/abs/hep-ph/9907240[hep-ph/9907240 ] ] .
a.  h. hoang , m.  jezabek , j.  h. khn and t.  teubner , _ radiation of heavy quarks _
lett . _ * b338 * ( 1994 ) 330335 [ http://arxiv.org/abs/hep-ph/9407338[hep-ph/9407338 ] ] .
a.  pineda and j.  soto , _ matching at one loop for the four - quark operators in nrqcd _ , _ phys .
_ * d58 * ( 1998 ) 114011 [ http://arxiv.org/abs/hep-ph/9802365[hep-ph/9802365 ] ] .
m.  beneke , a.  signer and v.  a. smirnov , _ two - loop correction to the leptonic decay of quarkonium _ , _ phys .
* 80 * ( 1998 ) 25352538 [ http://arxiv.org/abs/hep-ph/9712302 [ hep - ph/9712302 ] ] .
a.  czarnecki and k.  melnikov , _ two - loop qcd corrections to the heavy quark pair production cross section in @xmath0 annihilation near the threshold _
_ * 80 * ( 1998 ) 25312534 [ http://arxiv.org/abs/hep-ph/9712222[hep-ph/9712222 ] ] .
o.  v. tarasov , a.  a. vladimirov and a.  y. zharkov , _ the gell - mann - low function of qcd in the three loop approximation _ , _ phys .
* b93 * ( 1980 ) 429432 .
w.  bernreuther , r.  bonciani , t.  gehrmann , r.  heinesch , t.  leineweber _
_ , _ two - loop qcd corrections to the heavy quark form - factors : the vector contributions _ , _ nucl.phys .
_ * b706 * ( 2005 ) 245324 [ http://arxiv.org/abs/hep-ph/0406046[hep-ph/0406046 ] ] .
b.  kniehl , a.  onishchenko , j.  piclum and m.  steinhauser , _ two - loop matching coefficients for heavy quark currents _ , _ phys.lett . _
* b638 * ( 2006 ) 209213 [ http://arxiv.org/abs/hep-ph/0604072 [ hep - ph/0604072 ] ] .
m.  beneke , _ perturbative heavy quark - antiquark systems _ , in : proceedings of the 8th international symposium on heavy flavor physics ( heavy flavors 8) , 25 - 29 july 1999 , southampton , england , [ http://arxiv.org/abs/hep-ph/9911490[hep-ph/9911490 ] ]
. m.  beneke , p.  falgari and c.  schwinn , _ soft radiation in heavy - particle pair production : all - order colour structure and two - loop anomalous dimension _
_ nucl.phys .
_ * b828 * ( 2010 ) 69101 [ http://arxiv.org/abs/0907.1443[arxiv:0907.1443 ] ] .
m.  beneke , p.  falgari and c.  schwinn , _ threshold resummation for pair production of coloured heavy ( s)particles at hadron colliders _ , _ nucl.phys . _
* b842 * ( 2011 ) 414474 [ http://arxiv.org/abs/1007.5414[arxiv:1007.5414 ] ] .
m.  b. voloshin , _ nonperturbative effects in hadronic annihilation of heavy quarkonium _ , _ sov .
* 40 * ( 1984 ) 662667 .
m.  b. voloshin , _ precoulombic asymptotics for energy levels of heavy quarkonium _ , _ sov
. j. nucl .
* 36 * ( 1982 ) 143 .
d.  eiras and j.  soto , _ effective field theory approach to pionium _ ,
_ phys.rev . _ * d61 * ( 2000 ) 114027 [ http://arxiv.org/abs/hep-ph/9905543[hep-ph/9905543 ] ] .
a.  sommerfeld , _ ber die beugung und bremsung der elektronen _ , _ annalen der physik _ * 403 * ( 1931 ) 257330 .
t.  appelquist , m.  dine and i.  muzinich , _ the static limit of quantum chromodynamics _ , _ phys.rev . _
* d17 * ( 1978 ) 2074 .
w.  fischler , _ quark - anti - quark potential in qcd _
phys . _ * b129 * ( 1977 ) 157174 .
a.  billoire , _ how heavy must be quarks in order to build coulombic q anti - q bound states _ , _ phys .
lett . _ * b92 * ( 1980 ) 343 .
y.  schrder , _ the static potential in qcd _ , _ ph.d .
thesis _ ( 1999 ) .
desy - thesis-1999 - 021 .
m.  peter , _ the static quark - antiquark potential in qcd to three loops _
* 78 * ( 1997 ) 602605 [ http://arxiv.org/abs/hep-ph/9610209[hep-ph/9610209 ] ] .
m.  peter , _ the static potential in qcd : a full two - loop calculation _ , _ nucl .
* b501 * ( 1997 ) 471494 [ http://arxiv.org/abs/hep-ph/9702245[hep-ph/9702245 ] ] .
y.  schrder , _ the static potential in qcd to two loops _ , _ phys .
lett . _ * b447 * ( 1999 ) 321326 [ http://arxiv.org/abs/hep-ph/9812205[hep-ph/9812205 ] ] .
a.  v. smirnov , v.  a. smirnov and m.  steinhauser , _ full result for the three - loop static quark potential _ , _ pos _ * radcor2009 * ( 2010 ) 075 [ http://arxiv.org/abs/1001.2668[arxiv:1001.2668 ] ] .
n.  brambilla , x.  garcia  i tormo , j.  soto and a.  vairo , _ the logarithmic contribution to the qcd static energy at nnnnlo _ , _ phys .
lett . _ * b647 * ( 2007 ) 185 [ http://arxiv.org/abs/hep-ph/0610143 [ hep - ph/0610143 ] ] . m.  steinhauser , _ private communication _ , _ _ ( 2010 ) .