Four loop renormalization of QCD in the Curci-Ferrari gauge

We renormalize quantum chromodynamics (QCD) when gauge fixed in the nonlinear Curci-Ferrari gauge to four loops in the modified minimal subtraction (MS) scheme. We reproduce the four loop QCD MS /3-function from the Slavnov-Taylor identity for this gauge which relates the coupling constant renormalization to the gluon, Faddeev-Popov ghost and gauge parameter anomalous dimensions. This is carried out for a nonzero gauge parameter, without having to evaluate a vertex function. The anomalous dimension of the BRST invariant dimension two gluon and ghost mass term is deduced from a similar Slavnov-Taylor identity for this gauge. Consequently we construct the renormalization group functions in the minimal momentum subtraction scheme to four loops. As a corollary we deduce the five loop /3-function and quark mass anomalous dimensions in the same scheme. We also outline the pros and cons of employing the Curci-Ferrari gauge to access the six loop QCD /3-function in the MS scheme.