Pion and kaon distribution amplitudes in the continuum limit

We present a lattice-QCD calculation of the pion, kaon, and η s distribution amplitudes using large-momentum effective theory. Our calculation is carried out using three ensembles with 2 + 1 + 1 flavors of highly improved staggered quarks, generated by MILC Collaboration, at 310-MeV pion mass with 0.06, 0.09, and 0.12 fm lattice spacings. We use clover fermion action for the valence quarks and tune the quark mass to match the lightest light and strange masses in the sea. The resulting lattice matrix elements are nonperturbatively renormalized in regularization-independent momentum-subtraction scheme and extrapolated to the continuum. We use two approaches to extract the x dependence of the meson distribution amplitudes: (i) we fit the renormalized matrix elements in coordinate space to an assumed distribution form through a one-loop matching kernel and (ii) we use a machine-learning algorithm trained on pseudo lattice-QCD data to make predictions on the lattice data. We found the results are consistent between these methods with the latter method giving a less smooth shape. Both approaches suggest that as the quark mass increases, the distribution amplitude becomes narrower. Our pion distribution amplitude has broader distribution than predicted by light-front constituent-quark model, and the moments of our pion distributions agree with previous lattice-QCD results using the operator production expansion.