where $\ell$ is the total orbital angular momentum, $n$ the radial quantum number. We also cross-check the state identification with the decay constants and the wave functions themselves (see Sect. \[sect:lfwfs\]).

Spectroscopy {#sect:spectroscopy}
------------

![The reconstructed charmonium (*left panel*) and bottomonium (*right panel*) spectra at $N_{\max}=L_{\max}=32$. The horizontal and vertical axises are $j^{\mathsf P\mathsf C}$ and invariant mass in GeV, respectively. Model parameters are listed in Table \[tab:model\_parametersI\]. Calculated states are marked by boxes to represent the spread of the mass eigenvalues in $m_j$ owing to violation of the rotational symmetry (see text). The mean mass spreads, i.e. the average heights of the boxes, are 17 MeV and 8 MeV for charmonium and bottomonium, respectively. The r.m.s. deviations of the masses from the PDG values are 31 MeV and 38 MeV for charmonium and bottomonium, respectively. See text for details. []{data-label="fig:spectroscopy"}](Figures/charm_spectrum_Runalf0p6_kap0p97_mc1p60_NL32_style8.pdf "fig:"){width="45.00000%"} ![The reconstructed charmonium (*left panel*) and bottomonium (*right panel*) spectra at $N_{\max}=L_{\max}=32$. The horizontal and vertical axises are $j^{\mathsf P\mathsf C}$ and invariant mass in GeV, respectively. Model parameters are listed in Table \[tab:model\_parametersI\]. Calculated states are marked by boxes to represent the spread of the mass eigenvalues in $m_j$ owing to violation of the rotational symmetry (see 