Two loop integrals and QCD scattering

We present the techniques for the calculation of one- and two-loop integrals contributing to the virtual corrections to 2→2 scattering of massless particles. First, tensor integrals are related to scalar integrals with extra powers of propagators and higher dimension using the Schwinger representation. Integration By Parts and Lorentz Invariance recurrence relations reduce the number of independent scalar integrals to a set of master integrals for which their expansion in є = 2 — D/2 is calculated using a combination of Feynman parameters, the Negative Dimension Integration Method, the Differential Equations Method, and Mellin-Barnes integral representations. The two-loop matrix-elements for light-quark scattering are calculated in Conventional Dimensional Regularisation by direct evaluation of the Feynman diagrams. The ultraviolet divergences are removed by renormalising with the MS scheme. Finally, the infrared singular behavior is shown to be in agreement with the one anticipated by the application of Catani's formalism for the infrared divergences of generic QCD two-loop amplitudes.