{"id":"b9c0bb4b-66fe-4235-8f87-ab538f69f927","arxiv_id":"2611.15628","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":5.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":10,"one_line_summary":"A non-relativistic model with Coulomb-plus-linear and spin-dependent potentials reproduces the masses of conventional charmonium and bottomonium states to within about 1%.","lead":"This paper models heavy quarkonium particles (like charmonium and bottomonium) using a simple quantum mechanics approach with a potential that mimics the strong force. Its main value lies in being a clear, pedagogical demonstration of how matrix methods, variational methods, and perturbation theory can reproduce real particle masses.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek/deepseek-v4-flash","headline":"Fitted parameters lack error bars; chi-squared and degrees of freedom not reported, so 'good agreement' cannot be statistically evaluated.","rationale":"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.","tokens_in":17762,"tokens_out":606,"duration_ms":41770,"concrete_test":"Re-perform a chi-squared fit of the model to the experimental masses, including experimental uncertainties from the PDG. Compute the reduced chi-squared (χ²/d.o.f.) for the full fit (Tables III, IV) and for the 'predictive' fit (Tables VII, VIII). Report the number of parameters, degrees of freedom, and the P-value. If χ²/d.o.f. is significantly less than 1, overfitting is likely; if >2, the fit is poor. Propagation of parameter errors to the mass predictions would further quantify the reliability.","verdict_should_be":"UNCHANGED","load_bearing_attack":"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.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek/deepseek-v4-flash","summary":"The paper develops a non-relativistic potential model for heavy quarkonium (charmonium and bottomonium) using a Gaussian expansion method (GEM) to solve the radial Schrödinger equation with a Coulomb-plus-linear confining potential and spin-dependent terms (hyperfine and spin–orbit). The model parameters are fitted separately to experimental spectra of charmonium and bottomonium via chi-squared minimization. The resulting masses are compared to data in Tables III–IV, with most percent errors below 1% for the full fit. A subset fit (low-lying states only) is used to claim predictive power for higher states, though deviations increase (average ~70 MeV for charmonium). The paper acknowledges convergence limitations for S-waves, the empirical nature of the basis extrapolation, and that the full fit has 'no apparent predictive power' (§IVA).","tokens_in":17510,"tokens_out":2683,"duration_ms":21358,"significance":"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.","major_comments":[{"comment":"No chi-squared, reduced chi-squared, or degrees of freedom are reported for the fits. The paper claims 'good agreement' based on percent errors, but without a statistical goodness-of-fit measure it is impossible to assess whether the model overfits the data given 10 free parameters (5 per sector) fitted to a limited number of states. The number of data points used in each fit should be stated explicitly, and the fit quality must be quantified.","section":"§IV, Tables III–IV","authors_acknowledge":null,"authors_acknowledge_p":null},{"comment":"Parameter uncertainties are absent. The parameters are obtained by minimization, but no error bars are provided. This omission makes it impossible to propagate uncertainties to the predicted masses, to compare with other potential models, or to assess the sensitivity of the results to the parameter choices. Confidence intervals should be reported (e.g., via bootstrap or covariance matrix from the chi-squared fit).","section":"§III.A, Table I","authors_acknowledge":null,"authors_acknowledge_p":null},{"comment":"The S-wave basis convergence is acknowledged to be incomplete at N_max=11 (the hyperfine matrix element for ηb(1S) changes by 1.6% from N_max=10 to 11, and the infinite-basis extrapolation shifts masses by 7–10 MeV). Since the fit is performed at N_max=11, the systematic uncertainty from basis truncation is not propagated into the fitted parameters or the reported masses. The authors should estimate how the extrapolated limit affects the parameter values and the spectral agreement, or explain why this uncertainty is negligible relative to other model errors.","section":"§III.B, Fig. 2, Table II, Table V","authors_acknowledge":"However, the figure also shows that theS−wave does not converge rapidly enough in the r <0.6 GeV−1 regime, where Coulombic and hyperfine effects are strongest.","authors_acknowledge_p":0.69},{"comment":"The 'predictive power' section uses a subset of low-lying states to fit parameters, then predicts the remaining masses. The average deviation for charmonium is 69.8 MeV, which is large (≈2% for several states). The authors interpret this as identifying exotic candidates, but without error bars on the predictions, it is unclear which discrepancies are significant. Moreover, the subset fit produces a different parameter set (Table VI), so the model does not provide a unique prediction. The authors should provide estimated prediction uncertainties (e.g., from parameter covariance) and a statistical test for whether the outliers are indeed inconsistent with the model.","section":"§IV.A, Tables VII–VIII","authors_acknowledge":null,"authors_acknowledge_p":null}],"minor_comments":[{"comment":"The statement in §IV.A that the model 'has no apparent predictive power' (sentence b174) contradicts the abstract and conclusion, which assert 'good agreement'. The authors should reconcile these statements and clearly define the model's status as a fit rather than a prediction.","section":"Abstract and §V","authors_acknowledge":null,"authors_acknowledge_p":null},{"comment":"The spin–orbit regulator r0 = 0.1 GeV^{-1} is introduced ad hoc. The text says it was chosen by decreasing until results were insensitive, but no quantitative sensitivity study is shown. A brief plot or table demonstrating insensitivity would strengthen the justification.","section":"§III, Eq. (8)","authors_acknowledge":null,"authors_acknowledge_p":null},{"comment":"The smearing parameter σ for the hyperfine interaction is listed in Table I but its definition in Eq. (9) uses the same symbol; clarify that σ is the Gaussian width parameter and that its values in Table I are in GeV.","section":"§III, Eq. (9)","authors_acknowledge":null,"authors_acknowledge_p":null},{"comment":"The expressions for matrix elements are given assuming ℏ=1, but earlier equations include ℏ explicitly. Ensure consistent notation throughout (e.g., Eq. (1) vs. Eq. (8)).","section":"§VII, Appendix","authors_acknowledge":null,"authors_acknowledge_p":null},{"comment":"The extrapolation formula M(N_max) = M_∞ + a/N_max + b/N_max^2 is empirical. The authors note it is not derived, but they should also discuss its reliability for other states beyond the two ground-state S-waves examined.","section":"§III.B, Eq. (24)","authors_acknowledge":null,"authors_acknowledge_p":null},{"comment":"The residuals plots are too small to read. Consider enlarging or presenting them as separate panels. Also label axes clearly (ΔM in MeV).","section":"Fig. 4","authors_acknowledge":null,"authors_acknowledge_p":null}],"recommendation":"major_revision"},"author_rebuttal":null,"desk_editor":{"model":"deepseek/deepseek-v4-flash","letter":"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.\n\nThe 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.\n\nOverall, 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.","headline":"A careful pedagogical calculation of heavy quarkonium spectra with standard potentials; honest about its limitations but lacking statistical rigor on the fit.","tokens_in":17928,"tokens_out":1054,"duration_ms":9302},"pith_extraction":{"msc":[],"pacs":[],"model":"deepseek/deepseek-v4-flash","headline":"A non-relativistic Schrödinger model with spin-dependent forces predicts heavy quarkonium masses with most errors below 1%, matching experiment for both charm and bottom mesons.","keywords":["heavy quarkonium","charmonium","bottomonium","non-relativistic potential model","Gaussian expansion method","spin-dependent interactions","mass spectrum","constituent quark model"],"falsifier":"If a future experiment discovers a conventional charmonium or bottomonium state whose mass deviates from the model prediction by more than 3-4% (well outside the average deviation), or if the predicted mass of a known conventional state (e.g., ηb(2S)) is later revised and moves away from the model value by more than 10 MeV, then the model's parameterization or its neglect of the tensor interaction would be disfavored.","tokens_in":18591,"feed_emoji":"🔬","tokens_out":4093,"duration_ms":32220,"temperature":0.7,"pith_summary":"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).","feed_headline":"Non-relativistic model hits quarkonium masses to 1%","feed_subtitle":"Schrödinger equation with spin-dependent forces reproduces charm and bottom meson spectra within measurement precision.","key_machinery":"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_discovery":"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.","pith_inferences":["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."],"forward_implications":["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."],"supporting_citations":[{"why":"Supplies lattice QCD results showing the static quark–antiquark potential has short-distance Coulombic and long-distance linear behavior, motivating the Cornell-potential choice.","marker":"[12,13]"},{"why":"Provides the Breit–Fermi reduction of one-gluon exchange from which the spin-dependent hyperfine and spin–orbit terms are derived.","marker":"[14]"},{"why":"Justifies the sign of the spin–orbit term by assuming the confining potential is scalar, which is needed to reproduce the observed fine-structure splitting pattern.","marker":"[15]"},{"why":"Introduces the Gaussian Expansion Method (GEM) used to expand the radial wavefunctions, yielding analytic matrix elements and fast convergence.","marker":"[16]"},{"why":"Provides the experimental masses for all quarkonium states used both in the fit and in the comparison tables, as well as the review of constituent quark models.","marker":"[7]"},{"why":"Previous work applying matrix mechanics to charmonium only; the present paper extends it to bottomonium and includes the spin–orbit interaction.","marker":"[5]"}],"fun_headline_variants":["Schrödinger model predicts quarkonium within 1%","Non-relativistic approach matches charm and bottom meson spectra","Variational solution yields heavy quarkonium masses to 1%","Matrix mechanics plus spin terms reproduces quarkonium data","Student-friendly quarkonium model hits 1% accuracy"],"cache_read_input_tokens":0,"weakest_assumption_plain":"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).","fun_headline_variants_meta":{"raw":{"variants":["Schrödinger model predicts quarkonium within 1%","Non-relativistic approach matches charm and bottom meson spectra","Variational solution yields heavy quarkonium masses to 1%","Matrix mechanics plus spin terms reproduces quarkonium data","Student-friendly quarkonium model hits 1% accuracy"]},"model":"deepseek/deepseek-v4-flash","effort":"low","cost_usd":0.000177,"raw_usage":{"cost":0.000176568,"is_byok":false,"cost_details":{"upstream_inference_cost":0.000176568,"upstream_inference_prompt_cost":8.6968e-05,"upstream_inference_completions_cost":8.96e-05},"total_tokens":1146,"prompt_tokens":826,"completion_tokens":320,"prompt_tokens_details":{"audio_tokens":0,"video_tokens":0,"cached_tokens":256,"cache_write_tokens":0},"completion_tokens_details":{"audio_tokens":0,"image_tokens":0,"reasoning_tokens":238}},"tokens_in":826,"tokens_out":320,"cost_source":"openrouter_usage","duration_ms":3505,"image_pages":[],"temperature":1.0,"billed_attempts":1,"reasoning_tokens":238,"transport_provider":"openrouter","cache_read_input_tokens":0,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-11-13T18:47:30.971437+00:00","model_set":{"reader":"deepseek/deepseek-v4-flash"},"pipeline_version":"pith-pipeline@v1.4.1-daily-deepseek-or-medium","input_kind":null,"input_chars":null,"input_total_chars":null,"passes_missing":null,"published_standing":null,"published_confidence":null,"objection_audit":null,"run_record":null,"basis":{"kind":"not_recorded","source":null,"chars_given":null,"chars_total":null},"receipt":null,"falsifier":"If a future experiment discovers a conventional charmonium or bottomonium state whose mass deviates from the model prediction by more than 3-4% (well outside the average deviation), or if the predicted mass of a known conventional state (e.g., ηb(2S)) is later revised and moves away from the model value by more than 10 MeV, then the model's parameterization or its neglect of the tensor interaction would be disfavored.","supporting_citations":[{"cited_title":"``Spin-dependent forces in QCD,'' Phys. Rev. D 23, 2724--2744","cited_arxiv_id":null,"evidence_quote":"Provides the Breit–Fermi reduction of one-gluon exchange from which the spin-dependent hyperfine and spin–orbit terms are derived."},{"cited_title":"``Dynamics of light and heavy bound quarks,'' Phys. Rev. D 12, 841--854","cited_arxiv_id":null,"evidence_quote":"Justifies the sign of the spin–orbit term by assuming the confining potential is scalar, which is needed to reproduce the observed fine-structure splitting pattern."},{"cited_title":"Hiyama, Y","cited_arxiv_id":null,"evidence_quote":"Introduces the Gaussian Expansion Method (GEM) used to expand the radial wavefunctions, yielding analytic matrix elements and fast convergence."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Provides the experimental masses for all quarkonium states used both in the fit and in the comparison tables, as well as the review of constituent quark models."},{"cited_title":"``The mass spectrum of quarkonium using matrix mechanics,'' Am.\\ J.\\ Phys.\\ 91, 34--41","cited_arxiv_id":null,"evidence_quote":"Previous work applying matrix mechanics to charmonium only; the present paper extends it to bottomonium and includes the spin–orbit interaction."}],"lane_accuracy":{"status":"not yet measured","lane":{"model":"deepseek/deepseek-v4-flash","input_kind":null},"url":"https://pith.science/accuracy"},"review_version":1}