{"id":"bd4ee315-8764-400b-ba38-cad67aa72481","arxiv_id":"2608.08515","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":4.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":6,"one_line_summary":"A fitted Killingbeck potential reproduces B, B_s, and B_c masses to under 0.5% for selected parameters, but its wavefunction at the origin is unreliable, inflating the B decay constant by up to 84%.","lead":"This paper uses a non-relativistic potential model with the Killingbeck potential to calculate the masses, decay constants, and other properties of B, B_s, and B_c mesons, comparing the matrix Numerov and variational methods. It finds the model reproduces meson masses within about 0.5% for selected parameters, but overestimates the B meson decay constant, exposing the model's short-distance limitations.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The <0.5% mass claim is only true for selected parameter sets; the B_s scan averages deviate by 1.2–1.5% from PDG.","rationale":"The reader's strongest_claim already identifies the mass-claim overstatement, and my stress-test agrees with that. However, the reader's weakest_assumption focuses on the Killingbeck fit at short distance and its effect on |psi(0)|^2 and the decay constant. That is a genuine model limitation, and the authors themselves acknowledge it, but it does not directly undermine the headline mass claim as sharply as the representativeness of the reported masses. I therefore treat the mass-claim representativeness as the single most load-bearing concern. The issue is not an internal inconsistency in the Numerov or variational solvers; it is a mismatch between the paper's broad '<0.5%' language and its own summary tables. This does not move the reader's conditional verdict, because the reader already recommended conditioning acceptance on clarifying the range of validity of the claim. The concrete check is straightforward and uses only data already in the paper.","tokens_in":19452,"tokens_out":6787,"duration_ms":75010,"concrete_test":"Recompute, directly from Tables II–VII, the mean mass and the fraction of scan rows within 0.5% of the PDG mass for each meson and each method, and compare with Table VIII. In particular, verify the B_s Numerov and variational averages: if they remain 5.300 and 5.287 GeV, the deviations are 1.2% and 1.5%, respectively, and the abstract/conclusion must be revised to attribute the <0.5% claim only to the specifically identified closest parameter sets, with the selection criterion defined.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central claim is the abstract and conclusion statement that the ground-state masses of B, B_s, and B_c deviate by less than 0.5% from PDG. The paper's own Table VIII contradicts this for B_s: the Numerov average is 5.300 ± 0.057 GeV and the variational average is 5.287 ± 0.057 GeV, versus the PDG value 5.366 GeV, i.e., deviations of about 1.2% and 1.5%. Moreover, many individual rows in Tables II–VII fall outside the 0.5% band (e.g., B masses below 5.253 GeV, B_s masses below 5.326 GeV, and B_c masses below 6.243 GeV). The <0.5% statement therefore holds only for the explicitly selected 'closest obtained' rows, and no selection rule is stated. This is load-bearing because mass accuracy is the paper's headline contribution; if the claim is restricted to the best-fit parameter sets, that restriction and its criterion must be stated, and the averages must be reported as the model's typical predictive output.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"The manuscript computes ground-state masses, leptonic decay constants, B-meson oscillation frequencies, and Isgur-Wise slope/curvature parameters for B, B_s, and B_c mesons in a non-relativistic potential model. The interaction is the Killingbeck potential, whose three parameters are obtained by an ordinary-least-squares fit to the Cornell potential over r in [0, 2.0] GeV^{-1}. The Schrödinger equation is solved with both the matrix Numerov method and a variational Gaussian trial wave function. The paper's central quantitative claim is that the ground-state masses deviate by less than 0.5% from PDG values, and the results are compared extensively with lattice QCD, QCD sum rules, and other potential-model calculations.","tokens_in":19674,"tokens_out":3958,"duration_ms":46176,"significance":"If the <0.5% mass claim were supported by the paper's own tables, the manuscript would be a useful cross-check of two independent numerical methods in a standard phenomenological framework, and the extensive comparison tables would make it a convenient reference for this class of models. The authors are also candid in the conclusion about the decay-constant overestimation and about the need for relativistic corrections. However, the headline accuracy claim is not supported by the reported averages, and the parameter-selection procedure is not sufficiently specified. The contribution is therefore incremental and needs revision before the central claims can be accepted.","major_comments":[{"comment":"The claim that the ground-state masses of B, B_s, and B_c deviate by less than 0.5% from PDG values is contradicted by the paper's own Table VIII. For B_s, the PDG mass is 5.366 GeV; the Numerov average is 5.300 +/- 0.057 GeV, a deviation of about 1.2%, and the variational average is 5.287 +/- 0.057 GeV, a deviation of about 1.5%. Many individual rows in Tables II-VII also fall outside the 0.5% band. Since mass accuracy is the headline result, the claim must be restricted to the explicitly selected 'closest obtained' rows, the selection rule must be stated, and the averages must be reported as the model's typical predictive output.","section":"Abstract and Section IV"},{"comment":"The OLS fit of the Killingbeck potential to the Cornell potential is performed over r in [0, 2.0] GeV^{-1}, but the Cornell potential contains a 1/r term and diverges at r = 0. The manuscript does not state whether r_min = 0 is actually used in Eq. (3), nor how the singular point is handled in the configuration matrix X of Eq. (6). This is not a purely technical detail: the fitted c' coefficient controls the short-distance behaviour, and |psi(0)|^2 enters both the spin-spin mass term in Eq. (28) and the decay constant in Eqs. (29)-(31). The short-distance sensitivity is visible in the large decay-constant overestimates, so the fit grid and the quality of the fit should be documented explicitly.","section":"Section II.A"},{"comment":"The selection of parameter sets is not described in a way that supports the stated accuracy claim. Section III.A says the tables list parameter sets 'that yield mass estimates in closest agreement with experimental results', but no algorithm or a priori criterion is given for choosing these sets among the allowed ranges of alpha_s and C, and b is fixed at 0.195 without justification in this paper. Because rows outside the 0.5% band appear throughout Tables II-VII, the headline claim could reflect selection bias unless the selection procedure is predefined and transparent. The authors should specify how the reported sets were chosen and present the full scan statistics as the model's predictive range.","section":"Section III.A and Tables II-VII"}],"minor_comments":[{"comment":"Reference [2] contains an editorial placeholder ('check for published proceedings version and DOI') and should be completed.","section":"References"},{"comment":"The section ordering is inconsistent: the oscillation-frequency subsection is labeled F but appears before the Isgur-Wise subsection labeled E, and the text refers to 'Equations (16) and (17)' for the oscillation frequency, although that quantity is defined in Eqs. (32) and (34).","section":"Sections II.D-F and III.D"},{"comment":"The column headers and signs in Tables IV-VII are inconsistent: Table IV labels Cpot in GeV^2 while Table II labels it in GeV, and several Cpot entries in Tables V-VII switch sign (e.g., positive 0.998, 0.978, 0.958) without explanation. Please check whether these entries are the Cornell constant C or another quantity.","section":"Tables IV-VII"},{"comment":"The text says the variational results are listed in 'tables 6, 7 and 8', but the actual variational tables are numbered V, VI, and VII; renumber or re-reference consistently.","section":"Section III.B"},{"comment":"The caption says solid and dashed lines represent fitted curves from Table II and Table V, but only two curves per panel are described; please clarify which datasets correspond to which line style.","section":"Figure 1"}],"recommendation":"major_revision","confidential_remarks":"The paper is within the scope of the journal and the numerical machinery appears workable, but the central '<0.5%' mass claim is not supported by the paper's own averages, and the parameter-selection procedure needs to be made explicit. I would not accept the manuscript in its present form; a revised version that honestly reports scan averages and qualifies the headline claim could be reconsidered. I also noted a heavy reliance on the authors' own prior work (Ref. [6] and related), which should be checked for appropriate contextualization."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Quick take: this is a legitimate, internally consistent potential-model benchmark study, not a new-physics result. The authors solve the Killingbeck-potential Schrödinger equation with two independent methods — matrix Numerov and a Gaussian variational ansatz — fit the potential parameters to the Cornell potential, and tabulate B, B_s, and B_c masses, decay constants, mixing frequencies, and Isgur–Wise slope and curvature. The two methods agree well with each other on masses, and the best-fit rows do land within 0.5% of the PDG values. The conclusion is also refreshingly explicit about the model's limitations: it admits the B decay constant is overestimated by 58–84% and that the Isgur–Wise slope falls below the Uraltsev bound. That honesty earns credit.\n\nThe main problem is packaging. The abstract and conclusion say the masses deviate by less than 0.5%, without saying this holds only for the closest-obtained parameter rows. Table VIII shows the Numerov average for B_s is 5.300 ± 0.057 GeV, about 1.2% away from the PDG 5.366 GeV, and the variational average is similar. Many individual rows in Tables II–VII sit outside the 0.5% band. This is not fatal — the model is being used as an interpolation device — but the claim is load-bearing and needs to be restated as best-fit results, not typical predictive output. The B_s oscillation frequency also shows large method-to-method disagreement (Numerov average 9.0 ps⁻¹, variational average 23.5 ps⁻¹, versus experimental 17.8 ps⁻¹), which deserves a sentence of explanation.\n\nSofter issues: the tables have sign and convention glitches, most notably the shifted and mis-signed C_pot entries in Table V; that looks like copy-paste and should be cleaned up. There is no code or detailed Numerov discretization supplied, so reproducibility is incomplete. The parameter fitting to the Cornell potential is sensible, but it does mean the Killingbeck parameters are not independently derived.\n\nWho is this for? Someone doing non-relativistic heavy-meson phenomenology who wants a quick cross-check of two solving methods and a compact comparison table. It is not a high-significance result, but it is honest and useful. I would send it to a serious referee with a clear request for revision: restate the 0.5% claim, fix the tables, add numerical details or code. It is below the threshold for acceptance as-is, but it is the kind of work a referee can materially improve.","headline":"A workmanlike potential-model scan with honest limitations: the <0.5% mass claim is true only for selected parameter rows, not the scan average, and the B-meson decay constant is off by 58–84%.","tokens_in":20241,"tokens_out":3603,"would_cite":true,"duration_ms":37517,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":["12.39.Pn","14.40.Nd"],"model":"deepseek-v4-flash","headline":"By fitting the Killingbeck potential to the Cornell potential and solving the Schrödinger equation with two independent methods, this paper reproduces the ground-state masses of B, B_s, and B_c mesons to within 0.5% of the experimental…","keywords":["heavy meson spectroscopy","Killingbeck potential","Cornell potential","matrix Numerov method","variational method","decay constants","Isgur-Wise function","quarkonium"],"falsifier":"Recompute the $B$ meson decay constant after replacing the OLS-fitted Killingbeck potential with one that matches the Cornell potential more faithfully at small $r$—for instance by extending the fit window, adding a logarithmic short-distance term, or sampling $r$ densely near zero—and check whether $f_B$ moves from the paper's 0.31–0.37 GeV toward the experimental and lattice values around 0.19–0.22 GeV. In parallel, a direct measurement of the $B$-meson Isgur–Wise slope exceeding 0.75 would contradict the paper's non-relativistic prediction of $\\rho^2\\approx0.34$.","tokens_in":19198,"feed_emoji":"⚛️","tokens_out":13254,"duration_ms":117042,"temperature":0.7,"pith_summary":"The paper asks whether a single, simple interaction—the three-parameter Killingbeck potential—can describe the ground state of the heavy mesons $B$, $B_s$, and $B_c$ when solved non-relativistically. The parameters are fixed once by an ordinary least-squares fit to the Cornell potential over the window $r\\in[0,2]\\,\\mathrm{GeV}^{-1}$, and the Schrödinger equation is then solved independently by the matrix Numerov method and by a variational Gaussian trial wavefunction. The headline finding is that the closest obtained masses sit within 0.5% of the experimental averages, and the two methods agree with each other across all computed observables. The paper also reports that decay constants, oscillation frequencies, and Isgur–Wise slope and curvature parameters follow, and that the model's accuracy improves as the lighter quark gets heavier: $B_s$ and $B_c$ decay constants land near experiment or theory, while the $B$ decay constant overshoots by 58% to 84%. A reader should care because, if the mass claim holds, this provides a nearly parameter-free, analytic cross-check for lattice QCD and sum-rule calculations.","feed_headline":"B, B_s, B_c masses land within 0.5% of experiment","feed_subtitle":"Two independent numerical methods agree on heavy-meson observables from one fitted potential.","key_machinery":"The machinery is a two-step pipeline. First, the three-knob Killingbeck potential $V(r)=a'r^2+b'r+c'/r$ is fitted to the Cornell potential $V(r)=-\\frac{4\\alpha_s}{3r}+br+C$ by ordinary least squares over $[0,2]\\,\\mathrm{GeV}^{-1}$, so the potential parameters come from a regression rather than from spectroscopy. Second, the radial Schrödinger equation is solved with two independent methods: the matrix Numerov method, which discretizes the kinetic and potential operators into $N\\times N$ matrices and solves a generalized eigenvalue problem with $O(\\Delta r^4)$ accuracy, and a variational method using a Gaussian trial wavefunction, chosen because the $a'r^2$ term makes a Gaussian the natural ground-state ansatz. All observables then follow from standard formulas: the pseudoscalar mass is $M_P=M_Q+M_{\\bar Q}+\\langle E\\rangle+\\langle H_{SS}\\rangle$ with $\\langle H_{SS}\\rangle\\propto \\alpha_s|\\psi(0)|^2/(M_Q M_{\\bar Q})$; the decay constant is $f_P=\\sqrt{12|\\psi(0)|^2/M_P}\\,\\bar C$ with a QCD correction; the oscillation frequency $\\Delta m_B$ is proportional to $f_{B_q}^2$; and the Isgur–Wise slope $\\rho^2$ and curvature $C$ are the first two moments of the ground-state density. The quantity that carries the argument is the wavefunction at the origin $|\\psi(0)|^2$: it sets the hyperfine splitting and the decay constant, and the paper identifies its short-distance behavior as the main source of error.","core_discovery":"The central discovery is that the Killingbeck potential $V_{\\rm kill}(r)=a'r^2+b'r+c'/r$, calibrated by ordinary least squares to the Cornell potential $V_{\\rm Cornell}(r)=-\\frac{4\\alpha_s}{3r}+br+C$ on $0\\le r\\le 2\\,\\mathrm{GeV}^{-1}$, reproduces the experimental ground-state masses of $B$, $B_s$, and $B_c$ to better than 0.5% for the best-fitting parameter sets, with the matrix Numerov and variational Gaussian methods giving statistically consistent results (mass standard deviations near 0.03 GeV). The same wavefunctions yield decay constants, mixing frequencies, and Isgur–Wise parameters. The paper's own comparison shows the method's limits: the $B$-meson decay constant is overestimated by roughly 58% (variational) and 84% (Numerov) because the fitted potential is too rigid at short distances and inflates $|\\psi(0)|^2$, while the $B_s$ and $B_c$ values are much closer to experiment and to other theoretical frameworks. The Isgur–Wise slope for $B$ comes out at $\\rho^2\\approx0.34$, above the model-independent lower bound of $1/4$ but below the stronger heavy-quark-limit bound of $3/4$, which the authors interpret as the expected signature of a purely non-relativistic treatment.","pith_inferences":["The fitted potential's short-distance failure is arguably baked into the ordinary least-squares procedure: over the window $[0,2]\\,\\mathrm{GeV}^{-1}$ the fit is dominated by the linear and quadratic terms, so the singular $c'/r$ term is under-constrained near $r=0$; a weighted regression that emphasizes small $r$ would probably fix $f_B$ at the cost of some mass accuracy.","Because the Gaussian trial state is exact for a pure $a'r^2$ potential, the close variational–Numerov mass agreement indirectly measures how harmonic the fitted potential is in the energy-dominating region; computing $\\langle b'r\\rangle$ would quantify that directly.","The paper's remark that a temperature-dependent quadratic term connects the Killingbeck form to screened potentials in quark–gluon plasma implies a concrete extension: recomputing ground-state masses and decay constants as $a'(T)$ grows would yield meson dissociation temperatures, which the authors mention but do not calculate."],"forward_implications":["If the sub-0.5% mass agreement holds, the Killingbeck-plus-OLS recipe becomes a cheap, parameter-light estimator for ground-state masses of heavy-heavy mesons, requiring only quark masses, $\\alpha_s$, $b$, and $C$ as inputs.","The close agreement between the matrix Numerov and variational Gaussian methods (mass spreads of about 0.03 GeV) supports the use of a simple Gaussian ansatz for these ground states in other heavy-heavy systems.","Because $\\Delta m_B$ grows with $f_{B_q}^2$, the model's overestimated $B$ decay constant forces an overestimated oscillation frequency: the paper obtains $0.89$–$0.98\\,\\mathrm{ps}^{-1}$ against the experimental $0.507\\,\\mathrm{ps}^{-1}$, so mixing observables are not reliable for $B$ mesons.","For $B_s$, the Numerov and variational decay constants bracket the experimental value, suggesting the framework can be useful for $B_s$ mixing phenomenology once the short-distance part is improved.","The Isgur–Wise slopes place the model in the non-relativistic regime; improving $\\rho^2$ to satisfy the heavy-quark-limit bound would require adding relativistic corrections, not changing the confining potential."],"supporting_citations":[{"why":"Defines the Killingbeck potential used as the paper's central interaction model.","marker":"[3]"},{"why":"Provides the matrix Numerov discretization used to solve the radial Schrödinger equation.","marker":"[5]"},{"why":"Introduces the Cornell potential form that the Killingbeck potential is fitted to.","marker":"[8]"},{"why":"Supplies the Cornell parameters $\\alpha_s$ and $C$ and their allowed ranges for the OLS fit.","marker":"[9]"},{"why":"Supplies the constituent quark masses and mass-relation inputs used to compute meson masses.","marker":"[13]"},{"why":"Gives the mass formula, spin–spin correction, decay-constant, and oscillation-frequency expressions adopted here.","marker":"[16]"},{"why":"Provides Isgur–Wise slope and curvature reference values and the relativistic-correction comparison.","marker":"[27]"},{"why":"Another source of Isgur–Wise parameters used as comparison for slope and curvature.","marker":"[28]"},{"why":"The experimental meson masses used as the benchmark for the sub-0.5% deviation claim.","marker":"[30]"}],"fun_headline_variants":["Two methods, one potential: heavy meson masses within 0.5%","B-family masses predicted to 0.5% by fitted potential","Matrix Numerov and variational agree on B meson masses","Killingbeck potential yields sub-0.5% mass errors for B mesons","Heavy meson mass accuracy: 0.5% with OLS-fitted potential"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The load-bearing premise is that the Killingbeck potential obtained by an ordinary least-squares fit to the Cornell potential on the single interval $0\\le r\\le 2\\,\\mathrm{GeV}^{-1}$ is accurate enough at short distances to give the wavefunction at the origin, $|\\psi(0)|^2$; the paper's own $B$-meson decay-constant excess (58–84%) is evidence that this premise is the weakest point.","fun_headline_variants_meta":{"raw":{"variants":["Two methods, one potential: heavy meson masses within 0.5%","B-family masses predicted to 0.5% by fitted potential","Matrix Numerov and variational agree on B meson masses","Killingbeck potential yields sub-0.5% mass errors for B mesons","Heavy meson mass accuracy: 0.5% with OLS-fitted potential"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000354,"raw_usage":{"total_tokens":1984,"prompt_tokens":1062,"completion_tokens":922,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":678,"completion_tokens_details":{"reasoning_tokens":821}},"tokens_in":678,"tokens_out":922,"duration_ms":9593,"temperature":1.0,"reasoning_tokens":821,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-14T04:33:51.739417+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Recompute the $B$ meson decay constant after replacing the OLS-fitted Killingbeck potential with one that matches the Cornell potential more faithfully at small $r$—for instance by extending the fit window, adding a logarithmic short-distance term, or sampling $r$ densely near zero—and check whether $f_B$ moves from the paper's 0.31–0.37 GeV toward the experimental and lattice values around 0.19–0.22 GeV. In parallel, a direct measurement of the $B$-meson Isgur–Wise slope exceeding 0.75 would contradict the paper's non-relativistic prediction of $\\rho^2\\approx0.34$.","supporting_citations":[{"cited_title":"Killingbeck, Phys","cited_arxiv_id":null,"evidence_quote":"Defines the Killingbeck potential used as the paper's central interaction model."},{"cited_title":"Pillai, J","cited_arxiv_id":null,"evidence_quote":"Provides the matrix Numerov discretization used to solve the radial Schrödinger equation."},{"cited_title":"Eichten, K","cited_arxiv_id":null,"evidence_quote":"Introduces the Cornell potential form that the Killingbeck potential is fitted to."},{"cited_title":"Kingkar, S","cited_arxiv_id":null,"evidence_quote":"Supplies the Cornell parameters $\\alpha_s$ and $C$ and their allowed ranges for the OLS fit."},{"cited_title":"Hoque, B","cited_arxiv_id":null,"evidence_quote":"Supplies the constituent quark masses and mass-relation inputs used to compute meson masses."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Gives the mass formula, spin–spin correction, decay-constant, and oscillation-frequency expressions adopted here."},{"cited_title":"Isgur-Wise function in a QCD inspired potential model with WKB Approximation","cited_arxiv_id":"1112.2800","evidence_quote":"Provides Isgur–Wise slope and curvature reference values and the relativistic-correction comparison."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Another source of Isgur–Wise parameters used as comparison for slope and curvature."}],"review_version":1}