“…A perturbative one-loop calculation [16,46] of the matrix elements of the Wilson-line operators on the lattice has shown two nontrivial features of these operators: linear divergences (similar to those found in the continuum [44]), in addition to the logarithmic divergences, and mixing among certain pairs of the original operators under renormalization. Studies for the elimination of the linear divergences have been made using various methods, such as the static quark potential [17,47,48], the gradient flow [49][50][51], the nonperturbative bare matrix elements of the Wilson-line operators [16,46], and the auxiliary field formalism [52], [20,[53][54][55]. A complete nonperturbative renormalization prescription, which relies on nonperturbative matrix elements of Wilson-line operators, is described in Ref.…”