  We present a high-statistics lattice QCD determination of the valence parton
distribution function (PDF) of the pion, with a mass of 300 MeV, using two very
fine lattice spacings of $a=0.06$ fm and 0.04 fm. We reconstruct the
$x$-dependent PDF, as well as infer the first few even moments of the PDF using
leading-twist 1-loop perturbative matching framework. Our analyses use both
RI-MOM and ratio-based schemes to renormalize the equal-time bi-local
quark-bilinear matrix elements of pions boosted up to 2.4 GeV momenta. We use
various model-independent and model-dependent analyses to infer the large-$x$
behavior of the valence PDF. We also present technical studies on lattice
spacing and higher-twist corrections present in the boosted pion matrix
elements.
