A Non-Relativistic Model of Heavy Quarkonium

This paper does exactly what it says: it builds a non-relativistic quarkonium model using a Cornell potential plus spin-dependent terms, solves it with a Gaussian expansion, and compares to data. The computation is clearly explained, the matrix elements are given explicitly, and the code availability (implied by the C++ mention) is a plus. For a student or a teacher looking for a self-contained project that ties together variational methods, matrix diagonalization, and perturbation theory, this is a solid resource. The inclusion of bottomonium and spin-orbit terms extends the earlier Belhouari paper, and the convergence study (Table II, Figure 4) is a nice touch. The authors are upfront about the model's limitations: they note that the tensor term is neglected, that the Gaussian basis struggles with S-wave cusps, and that the full fit has “no apparent predictive power.” That honesty deserves credit.

The soft spots are real but not fatal for the intended audience. The main one—flagged by the stress test—is the lack of any error bars on the fitted parameters. Without chi-squared values or degrees of freedom, the “good agreement” in Tables III and IV cannot be statistically evaluated, and it is hard to judge whether the 10-parameter fit (5 per sector) is overfitting the roughly 20 conventional states per sector. The authors partly address this with the predictive-power section, but even there no parameter uncertainties are given. For a research paper this would be a serious gap; for a pedagogical paper it is a weakness that a sufficiently advanced instructor could use as a classroom discussion point. The basis extrapolation is empirical and examined only for the ground S-waves, which the authors acknowledge. Also, the charmonium fit in the predictive section shows a 69.8 MeV average deviation—compared to 26.8 MeV in the full fit—which suggests that the full fit massages the parameters quite a bit.

Overall, this is not a research paper that advances the theory of quarkonium. It is a well-written teaching tool that shows how standard quantum mechanics techniques can produce believable numbers. I would recommend it for peer review in a pedagogy-focused journal like the American Journal of Physics, where the transparency and the pedagogical structure are valued. The authors should be asked to add a brief statement about the absence of parameter error bars and to consider reporting chi-squared per degree of freedom for the fits, but the work is otherwise sound enough to warrant a referee round.

The central claim—that the theoretical mass spectrum agrees with experimental values—rests on percent errors shown in Tables III and IV. However, the paper does not report the chi-squared value, reduced chi-squared, or degrees of freedom from the fit. Without these, it is impossible to assess whether the agreement is statistically significant or whether the model overfits the data by using 10 fitted parameters (5 per sector) for a limited number of states. The paper itself acknowledges that the model has 'no apparent predictive power' (sentence b174) when fitted to the full spectrum, and the subsequent 'predictive power' section shows larger deviations (average 69.8 MeV for charmonium). This raises a concern: the good agreement in the full fit may be a consequence of overfitting rather than the validity of the model. Parameter uncertainties are absent, so confidence intervals for the masses cannot be computed, and the sensitivity of the results to small changes in parameters is unknown. Furthermore, the extrapolation to infinite basis (Fig. 4) shows a systematic shift of ~7–10 MeV below the quoted masses, yet parameters were fit at Nmax=11; the impact of this on the fit quality is not propagated.

The reader noted that the main results are fits with many parameters, so agreement is expected. My concern extends this by highlighting the absence of any statistical measure of fit quality (chi-squared, d.o.f., parameter uncertainties). This makes it impossible to evaluate the claim quantitatively. However, the reader's verdict already sets the bar at CONDITIONAL, and this additional concern does not change that assessment; it merely strengthens the rationale for conditionality. The paper's own limitations are well acknowledged, but the missing statistical rigor is a load-bearing gap for the central claim. Hence I recommend the verdict remain unchanged.

If the model were genuinely predictive, it would provide a simple, pedagogical tool for heavy quarkonium spectroscopy. However, the paper's central claim of 'good agreement' is weakened by the absence of statistical fit diagnostics and the self-admitted lack of predictive power in the full fit. The methods (GEM, variational diagonalization, first-order perturbation) are standard; the novelty lies in combining them for both quarkonium systems with explicit wavefunction convergence studies. The work may be suitable for an educational journal, but for a research journal in hep-ph, the statistical reporting is insufficient and the circularity of fitting to the same data used to claim agreement undermines the scientific contribution.

This paper develops a non-relativistic, spin-dependent phenomenological model of heavy quarkonium (charmonium and bottomonium). The Hamiltonian includes a Coulomb-plus-linear confining potential plus hyperfine and spin-orbit interactions, derived from one-gluon exchange. The radial Schrödinger equation is solved by expanding in a Gaussian basis, diagonalizing the resulting matrix, and varying a width parameter to minimize the energy. The model reproduces the experimental masses of nearly all well-established quarkonium states; the average deviation is 26.8 MeV for charmonium and 14.6 MeV for bottomonium, with most percent errors below 1%. By fitting only a few low-lying states, the model can predict higher excitations and identify states that likely have non-conventional structure, such as the X(3872) and Z_c(3900).

What carries the argument. The central object is the non-relativistic Hamiltonian $H = 2m_q - \frac{1}{m_q}\big(\frac{d^2}{dr^2}+\frac{2}{r}\frac{d}{dr}-\frac{\ell(\ell+1)}{r^2}\big) + V_{\text{SI}}(r) + V_{\text{SD}}(r)$ with $V_{\text{SI}} = -\frac{4}{3}\frac{\alpha_s}{r}+br+C$ and $V_{\text{SD}} = \frac{32\pi\alpha_s}{9m_q^2}\tilde\delta_\sigma(r)\mathbf{S}_q\cdot\mathbf{S}_{\bar q} + \frac{1}{m_q^2}\big(\frac{2\alpha_s}{3r^3}-\frac{b}{2r}\big)\mathbf{S}\cdot\mathbf{L}}$. The radial wavefunctions are expanded in a Gaussian basis $\psi_n(r)\propto r^{\ell}e^{-n\beta^2 r^2/2}$, leading to analytic matrix elements; the generalized eigenvalue problem is solved by matrix diagonalization, the variational parameter $\beta$ is optimized for each state, and the spin–orbit term is added perturbatively after the diagonalization.

Core claim. The central claim is that a Hamiltonian containing a static, spin-independent Cornell potential (Coulomb + linear) together with a regularized hyperfine contact term and a spin–orbit term, when solved by the Gaussian expansion method with variational optimization and first-order perturbation theory for the spin–orbit piece, yields a spectrum of heavy quarkonium states that agrees with experimental data to within roughly 1% for most states. The bottomonium predictions are more accurate than those for charmonium, consistent with the heavier bottom quark making the non-relativistic approximation more reliable. When the model is fitted to only the lowest S-, P-, and D-wave multiplets and used to predict the remaining states, it still achieves agreement within about 1% for conventional states while flagging candidates for exotic hadrons (e.g., X(3872), Z_c(3900), ψ(4230)) as outliers with errors of several percent.

The premise it rests on. The non-relativistic Schrödinger equation with a static potential is an adequate description of heavy quarkonium; this is least reliable for charmonium (where relativistic corrections are non-negligible) and for states with orbital angular momentum L>0 (where the neglected tensor interaction may contribute).

The model can be used to predict masses of yet-unobserved high-excitation states of charmonium and bottomonium, providing targets for experimental searches.

It naturally flags states with large deviations (e.g., X(3872), Z_c(3900), ψ(4230)) as candidates for exotic hadrons (tetraquarks, hybrid mesons, molecules), thereby serving as a filter for identifying non-conventional structures.

The close values of the confining parameter b for charm and bottom suggest that the linear confining potential may be flavor-independent, while the different α_s values argue for a running coupling that depends on the heavy-quark scale.

The good agreement for S-waves but larger errors for P- and D-waves indicates that the omitted tensor interaction and relativistic corrections become more important for higher orbital angular momentum.

The paper offers a pedagogical bridge: a student who understands quantum mechanics can, by implementing the described matrix and variational methods, reproduce state-of-the-art hadron spectrum results and explore the boundaries of the constituent quark model.

The systematic deviation for χc0(1P) (3.36% error) and ψ3(3842) (1.53%) suggests that the tensor interaction, though neglected here, may be required to achieve sub-percent accuracy for all P- and D-wave states.

If the model were extended to include a running coupling constant α_s(Q^2) and the tensor term, it might be possible to fit both charmonium and bottomonium with a single universal parameter set, further testing the flavor independence of confinement.

This review was written by a machine. No person has reviewed this paper on Pith yet. A scientist listed on Pith can review it for a fee, and their review is published here whatever it concludes.

Request a human review