The parton distribution functions (PDFs) are crucial to the understanding of the internal structure of hadrons, and their precise determination is necessary for searches of BSM physics in collider experiments. Directly accessing PDFs in Lattice QCD calculations is impossible due to the Euclidean space-time geometry. The moments of PDFs are defined in terms of local operators that can be computed in Lattice QCD, but suffer from power divergent mixing beyond a certain order due to the hypercubic symmetry of the lattice. Recently, a method was proposed [1] to use gradient flow in order to circumvent this mixing and access moments of PDFs of any order. In this talk, I will summarize this method, and present preliminary results on moments of the unpolarized isovector PDF of the pion using four stabilized Wilson fermion ensembles generated by the OpenLat [2] initiative.
[1] A. Shindler 2311.18704