Gradient flow step-scaling function for SU(3) with ten fundamental flavors

We calculate the step scaling function, the lattice analog of the renormalization group $\beta$-function, for an SU(3) gauge theory with ten fundamental flavors. We present a detailed analysis including the study of systematic effects of our extensive data set generated with ten dynamical flavors using the Symanzik gauge action and three times stout smeared M\"obius domain wall fermions. Using up to $32^4$ volumes, we calculate renormalized couplings for different gradient flow schemes and determine the step-scaling $\beta$ function for a scale change $s=2$ on up to five different lattice volume pairs. In an accompanying paper we discuss that gradient flow can promote lattice dislocations to instanton-like objects, introducing nonperturbative lattice artifacts to the step scaling function. Motivated by the observation that Wilson flow sufficiently suppresses these artifacts, we choose Wilson flow with the Symanzik operator as our preferred analysis. We study systematic effects by calculating the step-scaling function based on alternative flows (Zeuthen or Symanzik), alternative operators (Wilson plaquette, clover), and also explore the effects of the perturbative tree-level improvement. Further we investigate the effects due to the finite value of $L_s$.