Leading-twist parton distribution amplitudes of S-wave heavy-quarkonia

The leading-twist parton distribution amplitudes (PDAs) of ground-state $^1S_0$ and $^3S_1$ $c\bar c$- and $b\bar b$-quarkonia are calculated using a symmetry-preserving continuum treatment of the meson bound-state problem which unifies the properties of these heavy-quark systems with those of light-quark bound-states, including QCD's Goldstone modes. Analysing the evolution of $^1S_0$ and $^3S_1$ PDAs with current-quark mass, $\hat m_q$, increasing away from the chiral limit, it is found that in all cases there is a value of $\hat m_q$ for which the PDA matches the asymptotic form appropriate to QCD's conformal limit and hence is insensitive to changes in renormalisation scale, $\zeta$. This mass lies just above that associated with the $s$-quark. At current-quark masses associated with heavy-quarkonia, on the other hand, the PDAs are piecewise convex-concave-convex. They are much narrower than the asymptotic distribution on a large domain of $\zeta$; but nonetheless deviate noticeably from $\varphi_{Q\bar Q}(x) = \delta(x-1/2)$, which is the result in the static-quark limit. There are also material differences between $^1S_0$ and $^3S_1$ PDAs, and between the PDAs for different vector-meson polarisations, which vanish slowly with increasing $\zeta$. An analysis of moments of the root-mean-square relative-velocity, $\langle v^{2m}\rangle$, in $^1S_0$ and $^3S_1$ systems reveals that $\langle v^4\rangle$-contributions may be needed in order to obtain a reliable estimate of matrix elements using such an expansion, especially for processes involving heavy pseudoscalar quarkonia.