Single-quark electromagnetic form factors of charmonium up to $J=2$

We calculate the single-quark electromagnetic form factors of a broad subset of charmonium, including $η_c(1S)$, $η_c(2S)$, $χ_{c0}(1P)$, $χ_{c0}(2P)$, $J/ψ(1S)$, $J/ψ(2S)$, $χ_{c1}(1P)$, $χ_{c1}(2P)$, $h_c(1P)$, $h_c(2P)$, $χ_{c2}(1P)$ and $χ_{c2}(2P)$, via a relativized quark model. The reference frame dependence of the results is estimated as the computational error. We compare our results with those of the lattice quantum chromodynamics (LQCD), the Dyson-Schwinger equation (DSE) and the basis light front quantization (BLFQ) approaches where available and we find that most of our results agree with the other results. We also predict the single-quark electromagnetic form factors of $χ_{c0}(2P)$, $χ_{c1}(2P)$, $h_c(1P)$, $h_c(2P)$, $χ_{c2}(1P)$ and $χ_{c2}(2P)$, where no direct comparisons are available.