# Two-loop divergences of scattering amplitudes with massive partons

## Abstract

We complete the study of two-loop infrared singularities of scattering amplitudes with an arbitrary number of massive and massless partons in non-abelian gauge theories. To this end, we calculate the universal functions and , which completely specify the structure of three-parton correlations in the soft anomalous-dimension matrix, at two-loop order in closed analytic form. Both functions are found to be suppressed like in the limit of small parton masses, in accordance with mass factorization theorems proposed in the literature. On the other hand, they are unsuppressed and diverge logarithmically near the threshold for pair production of two heavy particles. As an application, we calculate the two-loop anomalous-dimension matrix for near threshold and show that it is not diagonal in the -channel singlet-octet basis.

###### pacs:

11.15.Bt,12.38.Bx,12.38.Cy,13.87.-a^{1}

## I Introduction

Recently, much progress has been achieved in the understanding of the infrared (IR) singularities of massless scattering amplitudes in non-abelian gauge theories. While factorization proofs guarantee the absence of IR divergences in inclusive observables (1), in many cases large Sudakov logarithms remain after this cancellation. A detailed control over the structure of IR poles in the virtual corrections to scattering amplitudes is a prerequisite for the resummation of these logarithms beyond the leading order (2); (3); (4); (5). Catani was the first to predict the singularities of two-loop scattering amplitudes apart from the pole term (6), whose general form was only understood much later in (7); (8); (9). In recent work (10), it was shown that the IR singularities of on-shell amplitudes in massless QCD can be derived from the ultraviolet (UV) poles of operator matrix elements in soft-collinear effective theory (SCET). They can be subtracted by means of a multiplicative renormalization factor, whose structure is constrained by the renormalization group. It was proposed in this paper that the simplicity of the corresponding anomalous-dimension matrix holds not only at one- and two-loop order, but may in fact be an exact result of perturbation theory. This possibility was raised independently in (11). Detailed theoretical arguments supporting this conjecture were presented in (12), where constraints derived from soft-collinear factorization, the non-abelian exponentiation theorem, and the behavior of scattering amplitudes in two-parton collinear limits were studied.

It is relevant for many physical applications to generalize these results to the case of massive partons. The IR singularities of one-loop amplitudes containing massive partons were obtained some time ago in (13), but until very recently little was known about higher-loop results. In the limit where the parton masses are small compared with the typical momentum transfer among the partons, mass logarithms can be predicted based on collinear factorization theorems (14); (15). This allows one to obtain massive amplitudes from massless ones with a minimal amount of calculational effort. A major step toward solving the problem of finding the IR divergences of generic two-loop scattering processes with both massive and massless partons has been taken in (16); (17). One finds that the simplicity of the anomalous-dimension matrix observed in the massless case no longer persists in the presence of massive partons. Important constraints from soft-collinear factorization and two-parton collinear limits are lost, and only the non-abelian exponentiation theorem restricts the allowed color structures in the anomalous-dimension matrix. At two-loop order, two different types of three-parton color and momentum correlations appear, whose effects can be parameterized in terms of two universal, process-independent functions and (17). Apart from some symmetry properties, the precise form of these functions was left unspecified. In this Letter we calculate these functions at two-loop order and study their properties in some detail.

## Ii Anomalous-dimension matrix

We denote by a UV-renormalized, on-shell -parton scattering amplitude with IR singularities regularized in dimensions. Here and denote the momenta and masses of the external partons. The amplitude is a function of the Lorentz invariants and , where the sign factor if the momenta and are both incoming or outgoing, and otherwise. For massive partons we define 4-velocities with and . We further define the recoil variables . We use the color-space formalism (18), in which -particle amplitudes are treated as -dimensional vectors in color space. is the color generator associated with the -th parton and acts as a matrix on its color index. The product is summed over . Generators associated with different particles commute, and is given in terms of the eigenvalue of the quadratic Casimir operator of the corresponding color representation, i.e., for quarks and for gluons. Below, we will label massive partons with capital indices () and massless ones with lower-case indices ().

It was shown in (10); (12); (17) that the IR poles of such amplitudes can be removed by a multiplicative renormalization factor , which acts as a matrix on the color indices of the partons. More precisely, the product is finite for after the coupling constant used in the calculation of the scattering amplitude is properly matched onto the coupling in the effective theory, in which the heavy partons are integrated out (17). The relation

(1) |

links the renormalization factor to a universal anomalous-dimension matrix , which governs the scale dependence of effective-theory operators built out of collinear SCET fields for the massless partons and soft heavy-quark effective theory fields for the massive ones. For the case of massless partons, the anomalous dimension has been calculated at two-loop order in (8); (9) and was found to contain only two-parton color-dipole correlations. It has recently been conjectured that this result may hold to all orders of perturbation theory (10); (11); (12). On the other hand, when massive partons are involved in the scattering process, then starting at two-loop order correlations involving more than two partons appear (16). At two-loop order, the general structure of the anomalous-dimension matrix is (17)

(2) | |||||

The one- and two-parton terms depicted in the first three lines start at one-loop order, while the three-parton terms in the last two lines appear at . The notation etc. means that the corresponding sum extends over tuples of distinct parton indices. The cusp angles are defined via

(3) |

They are associated with the hyperbolic angles formed by the time-like Wilson lines of two heavy partons. The physically allowed values are (one parton incoming and one outgoing), corresponding to , or (both partons incoming or outgoing), corresponding to with real . These possibilities correspond to space-like and time-like kinematics, respectively. Since in a three-parton configuration there is always at least one pair of partons either incoming or outgoing, at least one of the or variables must be time-like, and hence the functions and have non-zero imaginary parts.

The anomalous-dimension coefficients and (for ) in (2) have been determined to three-loop order in (12) by considering the case of the massless quark and gluon form factors. Similarly, the coefficients for massive quarks and color-octet partons such as gluinos have been extracted at two-loop order in (17) by analyzing the anomalous dimension of heavy-light currents in SCET. In addition, the velocity-dependent function has been derived from the known two-loop anomalous dimension of a current composed of two heavy quarks moving at different velocity (19); (20).

Here we complete the calculation of the two-loop anomalous-dimension matrix by deriving closed analytic expressions for the universal functions and , which parameterize the three-parton correlations in (2).

## Iii Calculation of and

To calculate the function we compute the two-loop vacuum matrix element of the operator , which consists of three soft Wilson lines along the directions of the velocities of three massive partons, without imposing color conservation. The anomalous dimension of this operator contains a three-parton term given by . The function follows from the coefficient of the pole in the bare matrix element of . We will then obtain from a limiting procedure.

The operator is renormalized multiplicatively, so that is UV finite, where is linked to the anomalous dimension in the same way as shown in (1). In order to calculate the two-loop factor, we have evaluated the two-loop non-planar and planar graphs shown in the first row of Figure 1, as well as the one-loop counterterm diagrams depicted in the second row. Contrary to a statement made in (16), we find that receives contributions from all five diagrams, not just from the non-planar graph. The most challenging technical aspect of the analysis is the calculation of the diagram involving the triple-gluon vertex. We have computed this diagram using a Mellin-Barnes representation and checked the answer numerically using sector decomposition (21). We have also checked that for Euclidean velocities our result for the triple-gluon diagram agrees numerically with a position-space based integral representation derived in (16). Combining all contributions, we find

(4) |

where we have introduced the functions

(5) | |||||

The function can be obtained from the above result by writing , and taking the limit at fixed . In that way, we obtain

(6) |

Whether a factorization of the three-parton terms into two functions depending on only a single cusp angle persists at higher orders in is an open question.

It is interesting to expand the two functions and about the threshold point with . We find

(7) | ||||

Based on the symmetry properties of and , it was concluded in (16); (17) that these functions vanish whenever two of the velocities of the massive partons coincide. Indeed, this seems to be an obvious consequence of the fact that is totally anti-symmetric in its arguments, while is odd in its second argument. This reasoning implicitly assumes that the limit of equal velocities is non-singular, but is invalidated by the presence of Coulomb singularities in and near threshold. Consider the limit where two massive partons 1 and 2 are produced near threshold, with a small relative 3-velocity defined in their rest frame. It is then straightforward to derive that

(8) |

where is the scattering angle formed by the 3-momenta of particles 1 and 3. A similar formula holds for . This result is anti-symmetric in the parton indices 1 and 2 as required, but it does not vanish at threshold.

Another interesting limit is that of large recoil, where all the scalar products become large in magnitude. In that case, both and are suppressed like , because for large

(9) |

Note that the non-planar contribution from the first graph in Figure 1, which was studied in the Euclidean region in (16), contains the leading-power term

(10) |

and is unsuppressed in this limit. However, this contribution cancels against a leading-power term in the planar and counterterm contributions.

Using that , we see that the large-recoil limit corresponds to at fixed . It follows that the three-parton correlation terms described by and vanish like in the small-mass limit. This observation is in accordance with a factorization theorem proposed in (14); (15), which states that massive amplitudes in the small-mass limit can be obtained from massless ones by a simple rescaling prescription for the massive external legs.

## Iv Anomalous dimension for near threshold

As a sample application, we apply our formalism to the calculation of the two-loop anomalous-dimension matrices for top-quark pair production near threshold in the channel. This matrix (along with the corresponding one in the channel) forms the basis for soft-gluon resummation at the next-to-next-to-leading logarithmic (NNLL) order. We adopt the -channel singlet-octet basis, in which the pair is either in a color-singlet or color-octet state. For the quark-antiquark annihilation process , we thus choose the independent color structures as and . In the threshold limit it is convenient to define the quantity , which is related to the relative 3-velocity between the top-quark pair in the center-of-mass frame by . We find that in the limit the two-loop anomalous-dimension matrices reduces to

(11) |

where the two-loop expressions for the anomalous dimensions , , and can be found in (17), and

(12) |

arises from the threshold expansion of the two-loop coefficient of the velocity-dependent cusp anomalous dimension .

We stress that, as a consequence of the Coulomb singularities, the three-parton correlation term does not vanish near threshold. Instead, it gives rise to scattering-angle dependent, off-diagonal contribution in (IV). The off-diagonal terms were omitted in two recent papers (22); (23), where threshold resummation for top-quark pair production was studied at NNLL order. We leave it to future work to explore if and how the results obtained by these authors need to be modified in light of our findings.

## V Conclusions

The IR divergences of scattering amplitudes in non-abelian gauge theories can be absorbed into a multiplicative renormalization factor, whose form is determined by an anomalous-dimension matrix in color space. At two-loop order this anomalous-dimension matrix contains pieces related to color and momentum correlations between three partons, as long as at least two of them are massive. This information is encoded in two universal functions: , describing correlations between three massive partons, and , describing correlations between two massive and one massless parton. In this Letter we have calculated these functions at two-loop order. Using the exact analytic expressions, we studied the properties of the three-parton correlations in the small-mass and threshold limits. We found that the functions and vanish as in the small-mass limit, in accordance with existing factorization theoreoms for massive scattering amplitudes (14); (15). On the other hand, and contrary to naive expectations, the two functions do not vanish in the threshold limit, because Coulomb singularities compensate a zero resulting from the anti-symmetry under exchange of two velocity vectors. This fact has been overlooked in the recent papers (22); (23), where the three-parton correlations were neglected near threshold.

Our results allow for the calculation of the IR poles in an arbitrary on-shell, -particle scattering amplitude to two-loop order, where any number of the partons can be massive. As an application, we have derived the anomalous-dimension matrix for top-quark pair production in the channel. We will explore in future work to what extent the new off-diagonal entries, arising from three-parton correlation terms, affect the numerical results for the threshold-resummed production cross sections at the Tevatron and LHC.

This Letter completes the study of IR divergences of two-loop scattering amplitudes with an arbitrary number of massive and massless external particles, and in arbitrary non-abelian (or abelian) gauge theories with massless gauge bosons. Details of our calculations will be presented in a forthcoming article.

Acknowledgments: We are grateful to Thomas Becher for collaboration during the early stages of this work.

### Footnotes

- preprint: MZ-TH/09-25

### References

- J. C. Collins, D. E. Soper and G. Sterman, Adv. Ser. Direct. High Energy Phys. 5, 1 (1988).
- G. Sterman, Nucl. Phys. B 281, 310 (1987).
- S. Catani and L. Trentadue, Nucl. Phys. B 327, 323 (1989).
- H. Contopanagos, E. Laenen and G. Sterman, Nucl. Phys. B 484, 303 (1997).
- N. Kidonakis and G. Sterman, Nucl. Phys. B 505, 321 (1997).
- S. Catani, Phys. Lett. B 427, 161 (1998).
- G. Sterman and M. E. Tejeda-Yeomans, Phys. Lett. B 552, 48 (2003).
- S. M. Aybat, L. J. Dixon and G. Sterman, Phys. Rev. Lett. 97, 072001 (2006).
- S. M. Aybat, L. J. Dixon and G. Sterman, Phys. Rev. D 74, 074004 (2006).
- T. Becher and M. Neubert, Phys. Rev. Lett. 102, 162001 (2009).
- E. Gardi and L. Magnea, JHEP 0903, 079 (2009).
- T. Becher and M. Neubert, JHEP 0906, 081 (2009).
- S. Catani, S. Dittmaier and Z. Trocsanyi, Phys. Lett. B 500, 149 (2001).
- A. Mitov and S. Moch, JHEP 0705, 001 (2007).
- T. Becher and K. Melnikov, JHEP 0706, 084 (2007).
- A. Mitov, G. Sterman and I. Sung, Phys. Rev. D 79, 094015 (2009).
- T. Becher and M. Neubert, Phys. Rev. D 79, 125004 (2009).
- S. Catani and M. H. Seymour, Phys. Lett. B 378, 287 (1996).
- G. P. Korchemsky and A. V. Radyushkin, Nucl. Phys. B 283, 342 (1987).
- N. Kidonakis, Phys. Rev. Lett. 102, 232003 (2009).
- A. V. Smirnov and M. N. Tentyukov, Comput. Phys. Commun. 180 (2009) 735.
- M. Beneke, P. Falgari and C. Schwinn, arXiv:0907.1443 [hep-ph].
- M. Czakon, A. Mitov and G. Sterman, arXiv:0907.1790 [hep-ph].