Correlator with tensor currents and two masses at two loops

We calculate the vacuum-to-vacuum correlator of two quark tensor currents with two massive quarks, retaining full momentum dependence. For the first time, we include perturbative corrections up to next-to-leading order. Our fully analytical expressions are provided in machine-readable form. Furthermore, we present numerical results for various input parameters, including an estimate of the scale uncertainties. Our results are essential input for applications of dispersive methods, including unitarity bounds and QCD sum rules.