  In applying large-momentum effective theory, renormalization of the Euclidean
correlators in lattice regularization is a challenge due to linear divergences
in the self-energy of Wilson lines. Based on lattice QCD matrix elements of the
quasi-PDF operator at lattice spacing $a$= 0.03 fm $\sim$ 0.12 fm with clover
and overlap valence quarks on staggered and domain-wall sea, we design a
strategy to disentangle the divergent renormalization factors from finite
physics matrix elements, which can be matched to a continuum scheme at short
distance such as dimensional regularization and minimal subtraction. Our
results indicate that the renormalization factors are universal in the hadron
state matrix elements. Moreover, the physical matrix elements appear
independent of the valence fermion formulations. These conclusions remain valid
even with HYP smearing which reduces the statistical errors albeit reducing
control of the renormalization procedure. Moreover, we find a large
non-perturbative effect in the popular RI/MOM and ratio renormalization scheme,
suggesting favor of the hybrid renormalization procedure proposed recently.
