Impact of gauge fixing precision on the continuum limit of nonlocal quark-bilinear lattice operators

We analyze the gauge fixing precision dependence of some nonlocal quark-bilinear lattice operators interesting in computing parton physics for several measurements, using five lattice spacings ranging from 0.032 to 0.121 fm. Our results show that gauge-dependent nonlocal measurements are significantly more sensitive to the precision of gauge fixing than anticipated. The impact of imprecise gauge fixing is significant for fine lattices and long distances. For instance, even with the typically defined precision of Landau gauge fixing of 10 − 8 , the deviation caused by imprecise gauge fixing can reach 12%, when calculating the trace of Wilson lines at 1.2 fm with a lattice spacing of approximately 0.03 fm. Similar behavior has been observed in ξ gauge and Coulomb gauge as well. For both quasi-parton-distribution functions (quasi-PDFs) and quasi-transverse-momentum-dependent PDFs operators renormalized using the regularization independent momentum subtraction (RI/MOM) scheme, convergence for different lattice spacings at long distance is only observed when the precision of Landau gauge fixing is sufficiently high. To describe these findings quantitatively, we propose an empirical formula to estimate the required precision.