Gradient flow step-scaling function for SU(3) with Nf=6 or 4 fundamental flavors

Nonperturbative determinations of the renormalization group (RG) β function are crucial to understand properties of gauge-fermion systems at strong coupling and connect lattice simulations and the perturbative ultraviolet regime. Choosing well-understood, QCD-like systems with SU(3) gauge group and either six or four fundamental flavors, we investigate their step-scaling β function. In both cases we push the simulations to the boundary of chiral symmetry breaking and study the regime g_GF²≲8.2 with six, and g_GF²≲6.6 with four flavors. We carefully consider the lattice discretization errors by comparing three different gradient flows (GF), and for each flow three operators to estimate the renormalized finite volume coupling. We also consider the tree level improvement of the coupling. Noteworthy outcome is that nonperturbatively determined β functions sits between the universal 2-loop result and the 3-loop gradient flow scheme prediction. At strong coupling the renormalized coupling runs significantly slower than the 2-loop prediction.