  We present a lattice-QCD calculation of the pion, kaon and $\eta_s$
distribution amplitudes using large-momentum effective theory (LaMET). Our
calculation is carried out using three ensembles with 2+1+1 flavors of highly
improved staggered quarks (HISQ), 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 (RI/MOM) scheme and extrapolated to the continuum. We use
two approaches to extract the $x$-dependence of the meson distribution
amplitudes: 1) we fit the renormalized matrix elements in coordinate space to
an assumed distribution form through a one-loop matching kernel; 2) 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.
