Over the past few years, significant progress has been made in computing parton distributions from lattice QCD within the framework of large momentum effective theory (LaMET). The light-cone parton distributions can be derived from gauge-invariant equal-time corrrelators with large momentum boost. In this talk, we introduce a new approach to calculate parton distributions from correlations of boosted quarks and gluons in the Coulomb gauge (CG), without Wilson lines. This new method can substantially improve statistical precision and simplify renormalization, thus providing a more efficient way to calculate parton distributions in lattice QCD, particularly for 3D distributions. Through a CG quasi-TMD, we derived the Collins-Soper (rapidity) evolution kernel of transverse-momentum-dependent PDFs (TMDs) up to $b_T$=4.2 GeV^{-1}, demonstrating its effectiveness and potential for future lattice calculations.