  Transverse-momentum-dependent parton distribution functions and wave
functions (TMDPDFs/TMDWFs) can be extracted from lattice calculations of
appropriate Euclidean matrix elements of staple-shaped Wilson line operators.
We investigate the mixing pattern of such operators under lattice
renormalization using symmetry considerations. We perform an analysis for
operators with all Dirac structures, which reveals mixings that are not present
in one-loop lattice perturbation theory calculations. We also present the
relevant one-loop matching in a renormalization scheme that does not introduce
extra non-perturbative effects at large distances, both for the TMDPDFs and for
the TMDWFs. Our results have the potential to greatly facilitate numerical
calculations of TMDPDFs and TMDWFs on the lattice.
