The perturbative QCD gradient flow to three loops

The gradient flow in QCD is treated perturbatively through next-to-next-to-leading order in the strong coupling constant. The evaluation of the relevant momentum and flow-time integrals is described, including various means of validation. For the vacuum expectation value of the action density, which turns out to be a useful quantity in lattice calculations, we find a very well-behaved perturbative series through NNLO. Quark mass effects are taken into account through NLO. The theoretical uncertainty due to renormalization-scale variation is significantly reduced with respect to LO and NLO, as long as the flow time is smaller than about 0.1 fm.