Complete two-loop QCD amplitudes for tW production at hadron colliders

We calculate the complete two-loop QCD amplitudes for hadronic $tW$ production by combining analytical and numerical techniques. The amplitudes have been first reduced to master integrals of eight planar and seven non-planar families, which can contain at most four massive propagators. Then a rational transformation of the master integrals is found to obtain a good basis so that the dimensional parameter decouples from the kinematic variables in the denominators of reduction coefficients. The master integrals are computed by solving their differential equations numerically. We find that the finite part of the two-loop squared amplitude is stable in the bulk of the phase space. After phase space integration and convolution with the parton distributions, it increases the LO cross section by about $3\%$.