{"id":"cf5bc61c-3d30-4169-ba39-e359fcd0ffe4","arxiv_id":"1908.07707","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":4.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":4,"one_line_summary":"Heavy quarkonia mass spectra and root mean square radii are reproduced by a hard-wall trigonometric Rosen-Morse potential with a small fitted conformal-symmetry-breaking parameter d.","lead":"A hard-wall potential made from cotangent and cosecant functions is used to model heavy quark-antiquark pairs. The model reproduces many charmonium and bottomonium masses with a small fitted parameter that the author interprets as a mild breaking of conformal symmetry.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The csc^2 replacement of the centrifugal term (Eq. 16) is uncontrolled for the fitted wave functions (r/a ~ 1), so the 'exact' spectra and the conformal-symmetry interpretation are for a different potential; the empirical agreement is also overstated (2S c-c off by 315 MeV).","rationale":"The reader's weakest_assumption correctly identifies the uncontrolled centrifugal approximation in Eq. (16) as the central technical flaw. My read agrees: the replacement 1/r^2 → csc^2(r/a)/a^2 is only justified for r/a << 1, but the fitted parameters and the explicit wave functions show the wave functions extend to regions where r/a ∼ O(1) or larger, making the approximation quantitatively unreliable. Because the 'exact' solution and the SO(4) conformal symmetry interpretation both depend on this replacement, the central claim that a small d describes heavy quarkonia via a conformally perturbed potential is not supported. Additionally, the empirical support is weaker than claimed: the tables show deviations of several hundred MeV for excited states, and the paper itself concedes that level splittings are strongly underestimated. These issues are addressable (by solving the original potential numerically or by reframing the model as the csc^2 potential and improving the fit), so a conditional verdict is appropriate. The proposed numerical test would settle whether the approximation changes the predictions substantially.","tokens_in":12507,"tokens_out":9589,"duration_ms":82873,"concrete_test":"Numerically integrate the radial Schrödinger equation for the original potential (11) including the true centrifugal term ℓ(ℓ+1)/r^2 on (0, a*π) with the parameters in Tables 1 and 2 (a, d, V0, quark masses), using a shooting method. Compare the resulting 1S, 2S, 1P, 3S masses and r.m.s. radii with Tables 1–3. If any level shifts by more than ~50 MeV relative to the approximate solution, the replacement in Eq. (16) is not valid and the paper's spectra and conformal-symmetry conclusions do not apply to the claimed potential. An even simpler check: compute the expectation value of a^2[1/r^2 − csc^2(r/a)/a^2] in the reported wave functions; if the ratio of this deviation to the kinetic energy is not small, the approximation fails.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The paper solves the Schrödinger equation after replacing the centrifugal term 1/r^2 by csc^2(r/a)/a^2 (Eq. 16), with the stated condition r/a << 1. However, the reported wave functions (Eqs. 49–54) and the fitted parameters (a = 0.56 fm for c-c, a = 0.29 fm for b-b) yield rms radii of 0.85–1.0 fm and 0.44–0.52 fm, respectively, so the wave functions have substantial support at r/a ∼ 1.5–1.8. At r/a = 1, the csc^2 term is about 1.41 times the true 1/r^2 term, and at r/a = π/2 the ratio is about 2.47. Thus the approximation is not valid over the integration domain, and the 'exact' energy formula (29) is the spectrum of a different potential (cot + [d(d+1)+ℓ(ℓ+1)] csc^2), not of Eq. (11) with the true centrifugal barrier. The conformal-symmetry interpretation rests on the csc^2 form that results from this replacement, so it does not carry over to the flat-space potential. Even setting this aside, the fit quality is poor: Table 1 shows the 2S c-c mass off by −315 MeV and 3S by −231 MeV, and Table 2 shows 2S b-b off by −234 MeV, directly contradicting the claim of 'good agreement with data.'","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper solves the radial Schrödinger equation for a hard-wall trigonometric Rosen-Morse potential, after replacing the centrifugal term 1/r^2 by csc^2(r/a)/a^2, using the Nikiforov-Uvarov method. It obtains closed-form energy eigenvalues and wave functions in terms of Romanovski polynomials, fits the potential parameters a and d to the 1S and 2P masses of charmonium and bottomonium, and presents predicted mass spectra and root-mean-square radii. The central claim is that a small conformal-symmetry-breaking parameter d suffices to describe heavy quarkonia and that conformal symmetry remains a viable dynamical symmetry in the heavy-flavor sector.","tokens_in":12939,"tokens_out":9000,"duration_ms":83204,"significance":"If the central claims were established, the paper would provide a simple analytic potential model for heavy quarkonia with closed-form wave functions and an explicit dynamical-symmetry interpretation. Strengths include the self-contained NU derivation, explicit polynomial wave functions, and detailed comparison tables for both c-cbar and b-bbar. However, the two load-bearing issues identified below—the inconsistent printed eigenvalue formula and the uncontrolled centrifugal approximation—currently prevent the results from supporting the conformal-symmetry interpretation. The empirical agreement is also more modest than the text claims, and the smallness of d is in part a fit outcome.","major_comments":[{"comment":"The printed energy formula is inconsistent with the derivation. Solving the second condition in Eq. (22) for the dimensionless energy gives E = (hbar^2/(2 mu a^2))[(alpha-1)^2 - V0^2/(4(alpha-1)^2)], i.e. Eq. (29) should contain (alpha-1)^2, not (-alpha-1)^2, in both places. For d=0, alpha = -(n+ell), so (alpha-1)^2 = (n+ell+1)^2, which reproduces Eq. (30). The printed version gives (-alpha-1)^2 = (n+ell-1)^2; for the 1P state (n=0, ell=1, d=0) it yields zero energy, contradicting the value 3.333 GeV in Table 1. Because every numerical result in the paper is generated from this formula, the inconsistency is load-bearing and must be corrected.","section":"Sec. 2.2, Eq. (29)"},{"comment":"The replacement 1/r^2 ≈ csc^2(r/a)/a^2 is introduced with the stated condition r/a << 1, but the wave functions are supported on the full interval [0, a*pi]. With the fitted values a=0.56 fm for c-cbar and a=0.29 fm for b-bbar, the rms radii in Table 3 (0.85-1.03 fm and 0.44-0.52 fm) correspond to r/a values around 1.5-1.8, where the csc^2 approximation differs from the true 1/r^2 term by factors of 1.4 at r/a=1 and 2.5 at r/a=pi/2. Thus Eq. (29) is the exact spectrum of a different potential, one with a csc^2 barrier, rather than of Eq. (11) with the true centrifugal term. The SO(4)/conformal interpretation follows directly from the csc^2 form and therefore does not automatically apply to the flat-space potential. The authors should quantify the error by solving the flat-space equation with the true 1/r^2 term numerically and comparing the resulting spectra.","section":"Sec. 2.2, Eq. (16)"},{"comment":"The claim of 'pretty good agreement' with data is not supported by the numbers. In Table 1 the 2S charmonium mass is off by -315 MeV and the 3S by -231 MeV; in Table 2 the 2S bottomonium mass is off by -234 MeV. These deviations are large compared with the typical splittings being modeled. Since the 1S and 2P masses are fitted by construction, the predictive content of the model is limited to the remaining states, and those show systematic deviations. The text should report a quantitative measure of fit quality and should not use 'good agreement' without qualification, especially in light of the paper's own statement in Sec. 4 that the level splittings are 'strongly underestimated.'","section":"Sec. 3.1, Tables 1-2"},{"comment":"The smallness of d is presented as evidence that conformal symmetry is only mildly broken, but d is one of two parameters fitted to the 1S and 2P masses. The fitted value of d is therefore an output of the fit, not a prediction, and cannot by itself support the claim that conformal symmetry is a viable dynamical symmetry. The paper should state this explicitly and identify a test of the symmetry-breaking parameter that does not rely on the fitted states, such as the rms radii or the relative splittings among the higher states.","section":"Sec. 3.1 and Sec. 4"}],"minor_comments":[{"comment":"The condition r/a << 1 is stated but never quantified; please specify the range over which the approximation is intended to hold and compare it with the actual support of the wave functions used in the fits.","section":"Sec. 2.2, Eq. (16)"},{"comment":"The ratio in Eq. (59) is computed from the model values in Table 3, yet the text describes the closed-form result in Eq. (63) as 'pretty close to the experimental one given in (59).' Equation (59) is not an experimental quantity; this wording should be corrected.","section":"Sec. 3.2, Eqs. (59) and (63)"},{"comment":"The closed-form rms expression uses the small-angle approximation sin^nu(r/a) ≈ (r/a)^nu and extends the integration to infinity, although the original problem has a hard wall at r=a*pi. The full integral is already computed for Table 3, so the approximate closed form is unnecessary and its use should be justified or removed.","section":"Sec. 3.2, Eqs. (60)-(62)"},{"comment":"The experimental values are cited collectively to Refs. [18,19], but individual states would benefit from explicit PDG names and uncertainties; some rows in the comparison column are left as '—' without explanation.","section":"Sec. 3.1, Tables 1-2"},{"comment":"There are numerous typographical and notation issues, including 'Such\"odinger,' inconsistent use of H_tRM(r/a) versus H(r/a), and the undefined symbols V_c and U0 in the text around Eq. (18); these should be cleaned up in a revision.","section":"General"}],"recommendation":"major_revision","confidential_remarks":"The manuscript is within the scope of a quarkonium potential-model paper, but the conformal-symmetry conclusion is currently overstated relative to the evidence. The inconsistency between Eq. (29) and Eqs. (22)/(30) is likely typographical, but the uncontrolled nature of the csc^2 centrifugal approximation is a more fundamental concern. I recommend major revision rather than rejection because the model could be reframed as an approximate description on S^3, or the flat-space connection could be defended with numerical checks; however, as it stands, the central claim is not supported."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Dear colleague,\n\nRead the Al-Jamel paper on heavy quarkonia with the trigonometric Rosen-Morse potential. My bottom line: it is a technically competent re-derivation of a known exact solution, but the central approximation undermines the physics, and the agreement with data is overstated.\n\nWhat is actually new: the paper takes the d≠0 trigonometric Rosen-Morse potential, solves it via the Nikiforov-Uvarov method, and applies it to charmonium and bottomonium. The wave functions and spectra are given explicitly, and r.m.s. radii are computed. The d parameter is fitted to the 1S and 2P masses, and the rest of the spectrum is predicted.\n\nThe math for d=0 is already in the literature (Kirchbach and Compean; Ciftci et al.), and the d extension is just a reparametrization of the csc^2 coefficient. So the novelty is mainly the application, not the solution.\n\nThe main problem is Eq. (16), where 1/r^2 is replaced by csc^2(r/a)/a^2 with the stated condition r/a << 1. But the fitted parameters give rms radii of about 1.5a for charmonium, so the wave functions have substantial support at r/a ~ 1.5–1.8. The approximation is not small there. This means the \"exact\" energy formula is not for the original potential with the true centrifugal barrier. The conformal symmetry interpretation rests entirely on the csc^2 form, which appears only after this replacement. So the conclusion that conformal symmetry remains viable in heavy-flavor QCD does not follow from the flat-space problem.\n\nThe fit quality is another soft spot. Table 1 has 2S c-c off by 315 MeV and 3S off by 231 MeV; Table 2 has 2S b-b off by 234 MeV. That is not \"pretty good agreement.\" The 1S and 2P states are fitted, so agreement there is by construction. The r.m.s. radii are compared to another theory, not to experiment, and the celebrated 2:1 ratio mostly reflects the fitted a values.\n\nThere is also an inconsistency in Eq. (29) relative to Eq. (22) and the d=0 limit in Eq. (30), likely a sign error in writing (-α-1) instead of (α-1). The n=0 case also gives α=0, violating the negativity condition stated earlier. These are fixable but need correction.\n\nWhat the paper does well: the NU method is applied carefully, the wave functions are explicit, and the authors honestly admit that level splittings are underestimated. But the central approximation is uncontrolled, and the physics claim overreaches.\n\nWho is this for? Someone interested in exact solutions of trigonometric potentials might find the technique useful, but the quarkonium phenomenology is not reliable. I would not cite it for the mass predictions. A serious editor should send it to peer review because the mathematical part is checkable and the application is a legitimate attempt, but I would expect major revision or rejection unless the approximation is justified or the model is reframed as genuinely living on S3 rather than flat space.\n\nRecommendation: do not desk-reject, but do not trust the results as they stand.\n\nYours,\n[Your name]","headline":"Competent NU-method solution of a known potential, but the centrifugal approximation invalidates the connection to the original potential, and the fit to quarkonia is significantly worse than claimed.","tokens_in":13392,"tokens_out":5134,"would_cite":false,"duration_ms":115589,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":["12.39.Pn","03.65.Ge","02.30.Ik","14.40.Lb","14.40.Nd"],"model":"deepseek-v4-flash","headline":"One small symmetry-breaking parameter reproduces heavy quarkonium masses and radii.","keywords":["heavy quarkonia","charmonium","bottomonium","hard-wall confinement","conformal symmetry","exact solutions","root mean square radii","dynamical symmetry"],"falsifier":"Numerically solve the radial wave equation for the same potential with the true centrifugal term $1/r^2$ and compare the excited-state masses to the closed-form values in Tables 1 and 2: if the shifts are larger than the experimental uncertainties, the reported ``exact'' spectrum belongs to the approximated potential, not the physical one.","tokens_in":12327,"feed_emoji":"⚛️","tokens_out":9254,"duration_ms":243839,"temperature":0.7,"pith_summary":"This paper's central claim is that heavy quark–antiquark bound states—charmonium and bottomonium—can be described by a hard-walled potential built from a cotangent term plus a squared-cosecant term, with a single fitted parameter $d$ that measures how strongly the heavy-quark masses break the conformal symmetry of the strong interaction. Solving the radial wave equation in closed form, the author obtains energy levels and wave functions, then fits $d$ and the wall size $a$ to the $1S$ and $2P$ masses of each meson family and predicts the remaining spectrum and the root-mean-square radii. The results match the measured masses and radii fairly well, while $d$ stays small ($0.109$ for $c\\bar c$ and $0.131$ for $b\\bar b$). A sympathetic reader would take this as evidence that a simple analytic potential with one symmetry-breaking parameter captures much of heavy quarkonium spectroscopy, and that conformal symmetry can survive as a dynamical symmetry in the heavy-flavor regime.","feed_headline":"One fitted parameter reproduces heavy quarkonium masses","feed_subtitle":"A hard-wall confining potential needs just one small symmetry-breaking parameter to match quarkonium masses and radii.","key_machinery":"The load-bearing object is the hard-wall potential $V(r/a) = -V_0 \\cot(r/a) + d(d+1)\\csc^2(r/a)$, solved for nonzero angular momentum by writing the centrifugal term as $1/r^2 \\approx \\csc^2(r/a)/a^2$. That identity converts the radial equation into the angular equation on a three-sphere $S^3$, whose isometry $SO(4)$ lies inside the conformal group $SO(2,4)$; the $d$-term then plays the role of a conformal-symmetry perturbation. The closed-form energies follow from a generalized hypergeometric solution method, and the wave functions are expressed through finite orthogonal polynomials with an orthogonality relation inherited from the sphere.","core_discovery":"The central discovery is that the exact solution of the hard-wall cotangent-plus-cosecant-squared potential, under the approximation that replaces $1/r^2$ by $\\csc^2(r/a)/a^2$, turns the radial problem into motion on a three-dimensional sphere, whose isometry group $SO(4)$ is the maximal compact subgroup of the conformal group $SO(2,4)$. In the $d=0$ limit the energy formula $E_{n\\ell} = \\frac{\\hbar^2}{2\\mu a^2}(n+\\ell+1)^2 - \\frac{\\hbar^2 V_0^2}{8\\mu a^2 (n+\\ell+1)^2}$ reproduces a spectrum degenerate in $N=n+\\ell+1$, and the small $d$ term removes those degeneracies within each $SO(4)$ multiplet but does not mix different multiplets. Fitting the two parameters $a$ and $d$ to the observed $1S$ and $2P$ masses of charmonium and bottomonium produces predicted masses and root-mean-square radii in reasonable agreement with experiment, including the observed $c\\bar c$ to $b\\bar b$ radius ratio near $2:1$. The author concludes that conformal symmetry remains a viable dynamical symmetry in all regimes of QCD, with $d$ quantifying its violation in the heavy-flavor sector.","pith_inferences":["Inference: because the cosecant substitution for the centrifugal term is accurate only for $r/a \\ll 1$, the ``exact'' formula may be an approximate description of the physical potential; a numerical solution with the true $1/r^2$ term would show how much the excited states change.","Inference: the clean mass dependence of $d$ suggests a future model could parameterize $d$ as a function of quark mass and extend the same potential to other flavors, such as strangeonium or toponium.","Inference: if conformal symmetry is genuinely dynamical here, the same wave functions should control transition rates and decay constants; computing those observables would test the model beyond masses and radii."],"forward_implications":["One parameter set $(a,d,V_0,\\mu)$ yields both the mass spectra and the root-mean-square radii for each quarkonium family.","At $d=0$, states with the same $N=n+\\ell+1$ are degenerate, explaining near-degeneracies such as $2S$–$1P$; the small fitted $d$ lifts these degeneracies without mixing multiplets.","The model reproduces the signs of the level splittings but underestimates their magnitudes, so improved splittings are expected from relativistic kinematics.","The charmonium-to-bottomonium radius ratio comes out close to the experimental $2:1$, driven mainly by the ratio of the fitted wall sizes $a$, indicating stronger localization of the $b\\bar b$ system.","The fitted $d$ increases from charm to bottom by a factor of about $1.6$, consistent with the larger heavy-quark mass causing a stronger conformal-symmetry violation."],"supporting_citations":[{"why":"Supplies the conformal-symmetry parameterization of the cotangent-plus-cosecant potential and its prior application to meson spectra.","marker":"[7]"},{"why":"Provides the exact and approximate solutions for this class of trigonometric potentials, including the $d=0$ limit reproduced here.","marker":"[12]"},{"why":"Supplies the generalized hypergeometric method used to obtain the closed-form energies and wave functions.","marker":"[15]"},{"why":"Gives the alternative theoretical mass spectra used as a comparison in Tables 1 and 2.","marker":"[17]"},{"why":"Provides the experimental quarkonium masses used for fitting and comparison.","marker":"[18]"},{"why":"Provides experimental data for the quarkonium masses in the tables.","marker":"[19]"},{"why":"Defines dynamical symmetry as degeneracy removal without multiplet mixing, the criterion used to conclude conformal symmetry remains viable.","marker":"[20]"}],"fun_headline_variants":["Tiny symmetry-breaking term captures quarkonium spectra","Small conformal-breaking tweak reproduces heavy quark data","Exact solution reveals quarkonia with minimal symmetry breaking","Hard-wall potential with tiny tweak matches quarkonium data","Quarkonium masses solved exactly via symmetry-breaking term"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The whole exact-solution construction rests on replacing the true centrifugal barrier $1/r^2$ by $\\csc^2(r/a)/a^2$, which is accurate only for $r/a \\ll 1$, while the wave functions extend all the way to the hard wall at $r = a\\pi$.","fun_headline_variants_meta":{"raw":{"variants":["Tiny symmetry-breaking term captures quarkonium spectra","Small conformal-breaking tweak reproduces heavy quark data","Exact solution reveals quarkonia with minimal symmetry breaking","Hard-wall potential with tiny tweak matches quarkonium data","Quarkonium masses solved exactly via symmetry-breaking term"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.001139,"raw_usage":{"total_tokens":4758,"prompt_tokens":1001,"completion_tokens":3757,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":617,"completion_tokens_details":{"reasoning_tokens":3680}},"tokens_in":617,"tokens_out":3757,"duration_ms":26476,"temperature":1.0,"reasoning_tokens":3680,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-14T11:59:34.587817+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Numerically solve the radial wave equation for the same potential with the true centrifugal term $1/r^2$ and compare the excited-state masses to the closed-form values in Tables 1 and 2: if the shifts are larger than the experimental uncertainties, the reported ``exact'' spectrum belongs to the approximated potential, not the physical one.","supporting_citations":[{"cited_title":"Kirchbach and C","cited_arxiv_id":null,"evidence_quote":"Supplies the conformal-symmetry parameterization of the cotangent-plus-cosecant potential and its prior application to meson spectra."},{"cited_title":"Hall, Nasser Saad, Exact and ap proximate solutions of Schr¨ odingers equation for a class of trigonometric potentials, Cent","cited_arxiv_id":null,"evidence_quote":"Provides the exact and approximate solutions for this class of trigonometric potentials, including the $d=0$ limit reproduced here."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Supplies the generalized hypergeometric method used to obtain the closed-form energies and wave functions."},{"cited_title":"Bhardwaj, and Fakir Chand, Mass Spectr a of Heavy and Light Mesons Using Asymptotic Iteration Method, Commun","cited_arxiv_id":null,"evidence_quote":"Gives the alternative theoretical mass spectra used as a comparison in Tables 1 and 2."},{"cited_title":"Tanabashi et al","cited_arxiv_id":null,"evidence_quote":"Provides the experimental quarkonium masses used for fitting and comparison."},{"cited_title":"Olive et","cited_arxiv_id":null,"evidence_quote":"Provides experimental data for the quarkonium masses in the tables."},{"cited_title":"Group Theory","cited_arxiv_id":null,"evidence_quote":"Defines dynamical symmetry as degeneracy removal without multiplet mixing, the criterion used to conclude conformal symmetry remains viable."}],"review_version":1}