Recently years have seen significant progress in the first-principles calculation of TMD physics from lattice QCD. We will describe the theoretical method for calculating both quark and gluon TMDs, which has been developed under the framework of large-momentum effective theory. Then we review its most recent applications to the non-perturbative quark TMDs and their rapidity evolution anomalous dimension, i.e., the Collins-Soper kernel, and discuss the control and improvement of systematic uncertainties in such calculations.