Parton distribution functions of ground state mesons composed of $c$ or $b$ quarks

The valence quark parton distribution functions (PDFs) of all ground state heavy mesons that composed of $b$ or $c$ quarks, are discussed; namely, the pseudoscalar $\eta_c(1S)$, $\eta_b(1S)$ and $B_c$, together with the corresponding vector ones, $J/\psi$, $\Upsilon(1S)$ and $B_c^\ast$. We use a QCD-inspired constituent quark model, which has been applied with success to conventional heavy mesons, so that one advantage here is that all parameters have already been fixed by previous studies. The wave functions of the heavy mesons in the rest frame are obtained by solving the Schr\"odinger equation, then boosted to its light-front based on Susskind's Lorentz transformation. The PDFs at the hadron scale, are then obtained by integrating out the transverse momenta of the modulus square of the light-front wave function. Our study shows how the valence quark distributions differ between pseudoscalar and vector mesons, as well as among charmonia, bottomonia and bottom-charmed mesons. Comparisons with other theoretical calculations demonstrate that the PDFs obtained herein are in general narrower but align well with the expected patterns. Moreover, each PDF's point-wise behavior is squeezed with respect to the scale-free parton-like PDF.