  We compute the master integrals for two-loop QCD corrections to quasi parton
distribution functions (PDFs) in large momentum effective theory. Analytical
results of the master integrals are derived using the method of differential
equations, along with a proper choice of canonical basis. The results of master
integrals are expressed in terms of Goncharov polylogarithms. These integrals
allow to extract the two-loop short-distance matching coefficients between
quasi and light cone PDFs in large momentum effective theory, and are helpful
to extract the nucleon PDFs from first principles.
