{"id":"3b88787e-e6d6-4839-8963-2723381fa548","arxiv_id":"1908.09131","paper_version":1,"verdict":"REJECT","confidence":"LOW","novelty_score":3.0,"correctness_risk":"high","formal_verification":"none","parameter_count":8,"one_line_summary":"Mass spectra are calculated from an extended Cornell potential with fitted parameters, and the reported agreement with experiment is partly guaranteed by fitting to the same states.","lead":"The paper solves the N-dimensional Schrödinger equation with a Cornell potential extended by quadratic and inverse-square terms, then reports masses for charmonium, bottomonium, and heavy-light mesons. It claims better agreement with experiment than several earlier works, but the agreement is produced by fitting potential parameters to the experimental masses it then compares against.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The inverse-square term d is re-derived per state via Eq. (18), so each tabulated state is solved with a different Hamiltonian; the reported agreement is a per-state fitting artifact.","rationale":"I read the paper as claiming that a single extended Cornell potential, with parameters fixed once per meson, predicts the masses of many states better than previous models. For that claim to hold, d must be the same in every state of a given meson. The paper never states this; instead, d is computed from Eq. (11)/(18), which depends on n and l. This makes the potential state-dependent, so each tabulated mass is the eigenvalue of a different Hamiltonian. This is not a minor technicality: it is exactly the mechanism the authors use to avoid the constraints that limited Ref. [6]. The reported improvement over previous work therefore does not demonstrate that the extended Cornell potential is physically better; it demonstrates that a per-state free parameter can lower the total error. I agree with the reader's weakest assumption, and I see no way to reinterpret the text as using a fixed d. A fixed-d recomputation would settle the matter. The reader's REJECT verdict is appropriate; no adjustment is needed.","tokens_in":7596,"tokens_out":12552,"duration_ms":126541,"concrete_test":"For charmonium, keep a=0.058 GeV and b=0.3366 GeV from Table 1, fix d to the value obtained from Eq. (18) for the 1S state (n=0, l=0, N=3), and solve the N=3 radial Schrödinger equation numerically for V(r)=a+b/r+c r+d/r^2 with c held at the 1S value. If the resulting 1P, 1D, 2S, 2P, 3S, 4S eigenvalues do not match the paper's P.W. column to within the quoted total error, the tabulated masses are not eigenvalues of one potential. Also check whether Eq. (18), with this fixed d, admits any solution for 1P or 2S; if not, the series terminates only for 1S.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central claim requires one fixed potential V(r)=a+b/r+c r+d/r^2 whose Schrödinger eigenvalues are the meson masses. Section 3 states that 'the value of parameter d was calculated by using Eq. (11)', and Eq. (11), restated as Eq. (18), explicitly contains n and l: (3n+2l)(3n+2l+2N-4)-4l(l+N-2)-8μd=0. Therefore d changes for every tabulated state. With a fixed d, Eq. (18) would restrict n and l, reproducing the very defect the authors attribute to Ref. [6]; the constraint is removed only by allowing d to be re-fit per state. The P.W. columns in Tables 1-6 are thus not eigenvalues of a single Hamiltonian. The quoted total errors (0.162, 0.074, 0.00003, 0.00001 GeV) measure the flexibility of per-state tuning, not the predictive power of a single extended Cornell potential.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The manuscript studies the mass spectra of heavy and heavy-light mesons using an extended Cornell potential V(r)=a+b/r+cr+d/r^2. The authors solve the N-dimensional radial Schrödinger equation by a power-series ansatz and derive an energy eigenvalue formula, Eq. (14). They then compute charmonium, bottomonium, D_s, and related heavy-light meson masses by relating the meson mass to the quark masses plus the eigenvalue, Eq. (16). The parameters a and b are fixed by fitting Eq. (16) to selected experimental masses, and the parameter d is stated to be calculated from Eq. (11). The paper reports total errors in Tables 1-6, compares with several earlier works, and claims improved agreement with experiment; it also discusses the dependence of the masses on the spatial dimensionality N.","tokens_in":7884,"tokens_out":6083,"duration_ms":67567,"significance":"If the model were shown to describe all states of a meson family with one fixed potential and without fitting the same states that are compared, the result could be a useful phenomenological contribution to quarkonium spectroscopy. The paper also correctly identifies a technical defect in the polynomial-termination condition of Ref. [6] and extends the calculation to heavy-light systems, which is a legitimate goal. However, as presented, the central claim is not established: the state-dependent determination of d and the fitting of potential parameters to the same experimental masses mean that the reported agreement is largely a fitting artifact. The manuscript therefore does not provide a credible prediction of the meson spectra from a single extended Cornell potential.","major_comments":[{"comment":"The central issue is the state dependence of the parameter d. Eq. (11), restated as Eq. (18), reads (3n+2l)(3n+2l+2N-4)-4l(l+N-2)-8μd=0, and Section 3 states that the value of d is calculated using Eq. (11). Because this relation contains n and l, d takes a different value for every state, so each entry in the P.W. columns of Tables 1-6 is computed with a different Hamiltonian V(r)=a+b/r+cr+d_{n,l}/r^2. With a fixed d, Eq. (18) would constrain n and l, reproducing the very defect the authors attribute to Ref. [6]; the constraint is removed only by re-fitting d per state. The reported spectra are therefore not eigenvalues of a single extended Cornell potential, and the comparison with experiment is not a test of the model's predictive power.","section":"§3, Eqs. (11), (18)"},{"comment":"The total errors quoted in Tables 1-6 are not defined, and the fitting procedure makes them uninformative. Section 3 states that a and b are determined by inserting experimental masses for selected states into Eq. (16); for example, charmonium uses the 2S and 2P states and bottomonium uses the 1S and 2S states. Those fitted states contribute zero or near-zero error by construction, so the total errors (0.162 GeV for charmonium, 0.074 GeV for bottomonium, and the much smaller values for the heavy-light mesons) measure fitting flexibility rather than predictive accuracy. The improvement claim in the abstract requires an error evaluation on states that were not used in the fit, or at minimum a clearly defined leave-one-out procedure.","section":"§4, Tables 1-6"},{"comment":"The derivation of the energy eigenvalue is not reproducible from the manuscript. After substituting the ansatz Eq. (3) and the series Eq. (5), the displayed Eq. (6) is garbled and the explicit recurrence relation for the coefficients is not written. The relations Eqs. (7)-(11) are stated as consequences of equating coefficients, but without the recurrence and the polynomial-termination condition the reader cannot verify the energy formula Eq. (14). The authors should supply the full recurrence, the termination condition, and the intermediate steps leading to Eqs. (7)-(11) and Eq. (14).","section":"§2, Eqs. (4)-(14)"},{"comment":"The manuscript does not clearly specify which experimental masses are used to fix a and b for each meson family. For charmonium the text says the 2S and 2P states are used, but Table 1 lists no experimental value for 2P; for other families several entries in the Experimental columns are also marked '-'. Without this information, the reader cannot determine which states are fitted and which are genuine predictions, and the error totals in Tables 1-6 cannot be audited.","section":"§3-§4, experimental input states"}],"minor_comments":[{"comment":"Eq. (15) is garbled as M=mq+qm; it should be written as M=m_q+m_{\\bar q}+E with the quark and antiquark masses clearly identified.","section":"Eq. (15)"},{"comment":"The units of the potential parameters are not stated consistently; for example, Table 1 lists a=0.058 GeV and b=0.3366 GeV, but in Eq. (1) the terms a, b/r, cr, and d/r^2 require different dimensions. The units of a, b, c, and d should be specified explicitly.","section":"Tables 1-6, parameter units"},{"comment":"The manuscript contains numerous typographical and formatting errors, including missing symbols in equations, inconsistent notation for the Schrödinger equation, and incomplete reference entries; these issues make the paper difficult to read and should be corrected throughout.","section":"General presentation"}],"recommendation":"reject","confidential_remarks":"I see no evidence of scientific misconduct, but the central physical claim is undermined by the state-dependent parameter d and by the fitting of the same experimental states that are later compared. These are not merely presentational issues: they affect the validity of the central conclusion. If the authors can reformulate the model so that one fixed potential describes all states of each meson family and can define and report a proper error metric on non-fitted states, a resubmission might be reconsidered; in the present form the manuscript is not suitable for publication."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Bottom line: this paper's central claim doesn't survive contact with its own equations. The coefficient d of the inverse-square term is computed separately for each state via Eq. (11)/(18), so every tabulated mass is obtained with a different potential. The reported agreement with experiment is a per-state tuning artifact, not a prediction from a single extended Cornell potential. That is a load-bearing flaw, not a fixable omission.\n\nWhat's genuinely here: the authors correctly identify a real defect in earlier work — Kumar and Chand's Eq. (17) unreasonably constrains n and l. They also extend the series-solution calculation to heavy-light mesons, which was not in Ref. [6]. The paper is organized clearly and compares against a reasonable set of references.\n\nThe problem is that their proposed remedy just moves the constraint into the inverse-square term. With a fixed d, Eq. (18) would again restrict n and l; the authors avoid that only by allowing d to change from state to state. That means Tables 1–6 are not spectra of a single potential. The total errors quoted (0.162 GeV for charmonium, 0.00001 for Bc) measure the flexibility that per-state d provides, not how well a fixed potential works.\n\nThere are also smaller, real omissions: the recurrence relation from Eq. (6) is never shown, so the derivation cannot be checked; the error definition is absent; potential parameters are only partially reported (table headers give a and b, but not c and d). The OCR is garbled in places, but these gaps are not OCR artifacts.\n\nWho gets value from this? Possibly someone wanting a quick mass formula, but the formula is not actually fixed, so the numbers are not trustworthy. The authors would need to treat d as a global parameter and face the constraint, or present this as a fitting exercise with the per-state freedom made explicit and physically justified. As written, the central argument fails.\n\nMy recommendation: desk reject. The load-bearing flaw is clear and the claimed improvement is an artifact. If the authors can reformulate with a fixed potential, there might be a modest phenomenological study worth looking at, but this version is not it.","headline":"The paper's central claim collapses: the inverse-square coefficient d is re-fitted per state, so each tabulated mass comes from a different Hamiltonian.","tokens_in":8338,"tokens_out":3874,"would_cite":false,"duration_ms":39744,"reading_group":"no","serious_thinker":"no","would_accept_peer_review":false},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":[],"model":"deepseek-v4-flash","headline":"The paper claims that extending the Cornell potential with an inverse-square term and solving the N-dimensional radial Schrödinger equation by power series yields improved masses for charmonium, bottomonium, and heavy-light mesons.","keywords":["Schrödinger equation","Cornell potential","heavy mesons","heavy-light mesons","power series method","meson mass spectra","quarkonium","N-dimensional radial equation"],"falsifier":"Compute every charmonium state with $a$ and $b$ fixed to the values in Table 1 and with $d$ set once from Eq. (11) using the 1S state, then compare the resulting total error with the reported 0.162 GeV; if the error grows substantially, the claimed improvement depends on per-state refitting.","tokens_in":7411,"feed_emoji":"⚛️","tokens_out":11251,"duration_ms":104230,"temperature":0.7,"pith_summary":"This paper tries to establish that an extended Cornell potential, which adds a quadratic term and an inverse-square term to the usual Coulomb-plus-linear potential, can describe the masses of heavy and heavy-light mesons more accurately than earlier potential-model calculations. The authors solve the N-dimensional radial Schrödinger equation with a power-series ansatz, obtain a closed-form energy formula, fix two potential parameters from selected experimental masses, and then compute the remaining states of the spectrum. They report total errors of 0.162 GeV for charmonium, 0.074 GeV for bottomonium, and as low as 0.00001 GeV for heavy-light systems, and they state that these results improve on recent Cornell-potential treatments. A reader would care because a simple potential model that produces a whole meson ladder from a few experimental inputs is a practical way to estimate states that have not been measured yet. The paper also claims that going to five spatial dimensions raises every computed mass, which it interprets as stronger binding in higher-dimensional space.","feed_headline":"Inverse-square term brings meson mass predictions closer to data","feed_subtitle":"A power-series solution with an inverse-square term reproduces charmonium, bottomonium, and heavy-light spectra.","key_machinery":"The machinery is the power-series solution of the N-dimensional radial Schrödinger equation. The radial wavefunction is written as $R(r)=e^{-\\alpha r^2-\\beta r}F(r)$, with $F(r)=\\sum_k c_k r^{k+\\sigma}$ and the offset $\\sigma$ chosen to avoid degeneracies. Substituting this ansatz converts the differential equation into a recursion for the coefficients, and requiring the series to terminate produces the energy eigenvalue formula (Eq. (14)) together with the quantization condition (Eq. (18)). The inverse-square coefficient $d$ enters through the term $8\\mu d$, and it is precisely that contribution which removes the earlier restriction to $n=0$. Meson masses are then computed from $M=m_q+m_{\\bar q}+E_{n,l}$, with quark masses and two potential parameters taken from experiment.","core_discovery":"The paper's central claim is that the quark-antiquark interaction $V(r)=a r^2+b r-\\frac{c}{r}+\\frac{d}{r^2}$, with positive parameters $a,b,c,d$, reproduces quarkonium and heavy-light meson spectra when the N-dimensional radial Schrödinger equation is solved by a terminating power series. The added inverse-square term is the load-bearing novelty: it changes the quantization condition from $(3n+2l)(3n+2l+2N-4)-4l(l+N-2)=0$, which the paper says leaves $n=0$ as the only acceptable radial state, to $(3n+2l)(3n+2l+2N-4)-4l(l+N-2)-8\\mu d=0$, which allows $n$ and $l$ to vary freely. With $a$ and $b$ fixed from two experimental states per meson and $d$ set by the quantization condition, the model generates the full S, P, D ladder and, according to the paper, lowers the total deviation from experiment to 0.162 GeV for charmonium, 0.074 GeV for bottomonium, and 0.00001-0.00003 GeV for the heavy-light systems. The paper further claims that in $N=5$ dimensions all computed masses increase, corresponding to stronger binding in higher dimensions.","pith_inferences":["The paper leaves open whether $d$ can be fixed once per meson; if it can, the model becomes a genuine predictive spectral tool, and testing that is the first follow-up.","Because the inverse-square term is what lifts the quantization constraint, the series method may extend to other power-law combinations whose coefficient recursion closes in a similar way.","The reported rise of masses with dimensionality suggests a concrete, testable consequence: higher-dimensional quarkonium should bind more deeply and could dissociate at different temperatures, which bears on extra-dimension searches.","The comparisons in the paper use energy levels only; a sharper test would compute hyperfine splittings or leptonic widths from the same wavefunctions, observables the current paper does not address."],"forward_implications":["If the extended potential is right, the same series solution can generate radial excitations for any quarkonium or heavy-light system without numerical integration.","The corrected quantization condition (Eq. (18)) removes the $n=0$ obstruction in the earlier series treatment, so a full ladder of S, P, and D states becomes available.","At $d=0$ the potential reduces to the extended Cornell potential treated by earlier Nikiforov-Uvarov calculations, so those results are included as a limiting case.","The reported total errors, if accepted, mean the model predicts unmeasured excited states, such as the 2P charmonium state listed as unknown in Table 1, with an estimated precision comparable to the measured low-lying states."],"supporting_citations":[{"why":"Supplies the series-solution method and the restrictive quantization condition (Eq. (17)) that the inverse-square term is designed to fix, and gives the bottomonium total error of 0.090 that the present work claims to improve.","marker":"[6]"},{"why":"Baseline asymptotic-iteration solution of the Cornell potential; the paper compares its charmonium and bottomonium total errors of 0.441 and 0.172 GeV against the present results.","marker":"[5]"},{"why":"Earlier Nikiforov-Uvarov solution of the extended Cornell potential at d=0; the paper states this is a limiting case of its potential and compares errors of 0.905 and 0.485 GeV.","marker":"[23]"},{"why":"Experimental charmonium and bottomonium masses used to fix the parameters a and b and to judge the agreement of the computed spectra.","marker":"[33]"},{"why":"Comparison source for heavy-light mesons whose total errors of 0.019 and 0.011 GeV the paper says the present calculation reduces.","marker":"[29]"},{"why":"Quasipotential-model calculation for B_c with total error 0.0008 GeV, one of the heavy-light benchmarks the paper says it improves.","marker":"[30]"}],"fun_headline_variants":["Inverse-square term brings meson masses closer to data","Power-series solution with r^-2 term refines meson spectra","Extended Cornell potential with inverse-square term improves fits","New term in Cornell potential yields precise heavy-light masses","Series method with r^-2 term tightens quarkonium predictions"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The load-bearing premise is that $d$ in the potential is a fixed physical parameter, but the paper appears to evaluate $d$ separately for each state through Eq. (11), which would mean the potential is refit state by state rather than being one potential for the whole meson.","fun_headline_variants_meta":{"raw":{"variants":["Inverse-square term brings meson masses closer to data","Power-series solution with r^-2 term refines meson spectra","Extended Cornell potential with inverse-square term improves fits","New term in Cornell potential yields precise heavy-light masses","Series method with r^-2 term tightens quarkonium predictions"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000197,"raw_usage":{"total_tokens":1342,"prompt_tokens":904,"completion_tokens":438,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":520,"completion_tokens_details":{"reasoning_tokens":357}},"tokens_in":520,"tokens_out":438,"duration_ms":5074,"temperature":1.0,"reasoning_tokens":357,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-14T11:20:31.217897+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Compute every charmonium state with $a$ and $b$ fixed to the values in Table 1 and with $d$ set once from Eq. (11) using the 1S state, then compare the resulting total error with the reported 0.162 GeV; if the error grows substantially, the claimed improvement depends on per-state refitting.","supporting_citations":[{"cited_title":"Abu-Shady, T","cited_arxiv_id":null,"evidence_quote":"Supplies the series-solution method and the restrictive quantization condition (Eq. (17)) that the inverse-square term is designed to fix, and gives the bottomonium total error of 0.090 that the present work claims to improve."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Baseline asymptotic-iteration solution of the Cornell potential; the paper compares its charmonium and bottomonium total errors of 0.441 and 0.172 GeV against the present results."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Earlier Nikiforov-Uvarov solution of the extended Cornell potential at d=0; the paper states this is a limiting case of its potential and compares errors of 0.905 and 0.485 GeV."},{"cited_title":"Gupta and I","cited_arxiv_id":null,"evidence_quote":"Comparison source for heavy-light mesons whose total errors of 0.019 and 0.011 GeV the paper says the present calculation reduces."},{"cited_title":"Godfrey and K","cited_arxiv_id":null,"evidence_quote":"Quasipotential-model calculation for B_c with total error 0.0008 GeV, one of the heavy-light benchmarks the paper says it improves."}],"review_version":1}