{"id":"e7b89406-fae3-4ffb-8473-c9b84a70cd0c","arxiv_id":"2501.03147","paper_version":2,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":4.0,"correctness_risk":"high","formal_verification":"none","parameter_count":10,"one_line_summary":"A screened relativistic potential model reproduces many low-lying bottomonium masses and proposes S-D mixing assignments for Upsilon(10355), Upsilon(10580), Upsilon(10753), Upsilon(10860), and Upsilon(11020).","lead":"This paper calculates masses, decay widths, and radiative transitions for bottomonium using a screened relativistic potential model. It then assigns several measured Upsilon resonances to mixtures of S-wave and D-wave bottomonium states.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Table XXI reproduces Eqs. (35)–(36) only if the quoted mixing angles are read as radians, so the reported degree-valued angles and the S–D assignments are not reproducible.","rationale":"The reader's conditional verdict is appropriate, but the main load-bearing weakness is sharper than the one highlighted in the reader's weakest_assumption. The reader focused on the two-state mixing assumption and the abstract/body contradiction; my check reveals a concrete internal inconsistency: the table's mixing angles are reported in degrees yet reproduce the quoted masses and widths only when interpreted as radians. This is not merely a modeling assumption but a formal error in the quantitative core of the S–D mixing claim. It undermines every mixing fraction, mass shift, and leptonic width derived from those angles. The abstract/body contradiction on Υ(10753) adds further ambiguity. These issues do not require overturning the otherwise useful pure-spectrum and decay catalog, but they do require the mixing analysis to be recomputed and re-reported before the central claim can be trusted. Since the reader already assigned CONDITIONAL, the verdict remains CONDITIONAL, and the concrete test above should be run as part of the revision process.","tokens_in":30509,"tokens_out":14784,"duration_ms":129445,"concrete_test":"Substitute θ=19.28°, −28.82°, 44.55° into Eq. (35) and Eq. (36) using the pure-state masses and Table VIII widths. The tabulated Mφ and Γφ are recovered only when the same θ values are interpreted as radians (e.g., 19.28 rad ≡ 0.43 rad; 44.55 rad ≡ 0.57 rad). If the degree-interpreted values disagree with Table XXI but the radian-interpreted values agree, the angle column is mislabeled and every mixing fraction derived from it is wrong.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central S–D mixing results in Table XXI are internally inconsistent. Substituting the quoted angles (θ=19.28°, −28.82°, 44.55°) into Eqs. (35)–(36) does not give the tabulated masses or leptonic widths. For 3S–2D, Eq. (35) with θ=19.28° gives Mφ≈10384.0 MeV, not 10374.9 MeV; for 5S–4D, θ=44.55° gives cos2θ≈0.016 and a split of more than 1 GeV, whereas Table XXI shows a 50 MeV split. The table only matches Eqs. (35)–(36) if the θ column is interpreted as radians (e.g., 19.28 rad ≡ 0.43 rad, giving cos2θ≈0.65). Since the adjacent column θ[61] is in degrees, the authors appear to have used radians in the formulas while reporting degrees, so the physical mixing amplitudes (sinθ, cosθ) are misreported. The abstract/body contradiction on Υ(10753) (pure 3D vs 4S–3D mixed) compounds the ambiguity. Because the entire interpretation of Υ(10355), Υ(10580), Υ(10860), and Υ(11020) rests on these numbers, the central claim is unsupported as written.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"The manuscript develops a relativistic screened potential model for bottomonium and computes the mass spectrum, decay constants, radiative (E1, M1) transitions, and annihilation widths for S, P, D, F, and G waves. Its main interpretive claim is that several observed vector resonances — Υ(10355), Υ(10580), Υ(10753), Υ(10860), and Υ(11020) — can be described as S–D mixed bottomonium states, with the mixing angles determined from leptonic widths. The paper contains an extensive set of tables comparing the model with experiment and with other potential-model calculations.","tokens_in":30887,"tokens_out":7798,"duration_ms":73293,"significance":"If the S–D assignments were established, the paper would support a conventional bottomonium interpretation of these resonances and provide a useful compendium of decay predictions for future searches. The pure-spectrum part is a genuine eigenvalue calculation and the decay tables are a potentially valuable resource. However, the central S–D mixing conclusion is not currently supported: the mixing angles are fitted to the same leptonic widths used in the comparison, the reported angles and masses are mutually inconsistent with the stated formulas, and the abstract and body disagree on the nature of Υ(10753). These issues affect the main novelty of the paper, not merely its presentation.","major_comments":[{"comment":"The S–D mixing analysis is circular as presented. The mixing angle θ is obtained by fitting Eq. (36) to the experimental Γ_ee values, and that same θ is then inserted into Eq. (35) to compute the mixed masses in Table XXI. Consequently, the agreement of Mφ and Mφ′ with experiment is not an independent test of the S–D assignments; it is a consistency check at best. To support the claim that S–D mixing correctly reproduces the observed masses, the authors must either determine θ from a source that does not already enter the mass comparison (e.g., an independent coupled-channel or line-shape analysis) or explicitly reframe the analysis as a fit and propagate all correlations and uncertainties.","section":"Section IV, Eqs. (35)–(36) and Table XXI"},{"comment":"The quoted mixing angles and mixed masses are arithmetically inconsistent. For the 3S–2D row, inserting θ=19.28° into Eq. (35) with M_S=10394.2 MeV and M_D=10467.3 MeV gives Mφ≈10384.0 MeV, not the tabulated 10374.9 MeV. For the 5S–4D row, θ=44.55° gives cos2θ≈0.016, which would produce a mass split of more than 1 GeV, whereas Table XXI shows a split of about 50 MeV. The tabulated values are reproduced only if the θ entries are interpreted as radians (e.g., 19.28 rad ≡ 0.430 rad). Since the adjacent column from Ref. [61] is explicitly in degrees, the physical mixing amplitudes (sinθ, cosθ) used in Eqs. (35)–(36) are misreported. This is a load-bearing error: every S–D assignment in Table XXII depends on these numbers, and the central claim is not assessable until the tables and formulas are brought into mutual consistency.","section":"Table XXI vs. Eq. (35)"},{"comment":"The abstract states that Υ(10753) is obtained as a purely Υ1(3D) bottomonium state, whereas Section V, Table XXI, and Table XXII treat Υ(10753) as the 4S–3D mixed partner with Mφ′=10772.9 MeV. This is not a cosmetic inconsistency: the nature of Υ(10753) is one of the paper's four central interpretive claims, and the reader cannot tell which assignment is being advocated. The abstract and body must be reconciled.","section":"Abstract vs. Section V and Table XXII"},{"comment":"The conclusion that the computed masses exhibit strong agreement with experiment is not quantitatively supported for the states involved in the S–D analysis. The pure 4S1, 5S1, and 6S1 masses in Table II are 108.7, 53.7, and 160.9 MeV above the experimental values, and after mixing the 3S–2D and 4S–3D masses remain 19.8 and 77.0 MeV above the corresponding resonances. No uncertainties are propagated through Eqs. (35)–(36), so statements such as 'align closely' cannot be evaluated statistically. The authors should provide uncertainties on θ, Mφ, Mφ′, and Γφ, and should temper the agreement claim accordingly.","section":"Tables II and XXII; Conclusion"}],"minor_comments":[{"comment":"The χ² definition uses artificial 0.1% mass errors rather than experimental uncertainties. This choice should be justified or replaced with a robust fitting procedure; as written, the quoted χ²=14.1 does not have the usual statistical interpretation.","section":"Eq. (18)"},{"comment":"The table layout is confusing: the M_S/M_D and Mφ/Mφ′ entries are combined in single cells, and the negative values of θ[61] are not explained. A cleaner layout would separate the two masses and the two widths.","section":"Table XXI"},{"comment":"There are several typographical issues, including 'masess' (Section V), inconsistent spacing in 'E1' and 'M1' transitions, and the duplicated reference [30]/[13] for the Belle observation of Υ(10753). A careful proofreading pass is needed.","section":"Throughout"}],"recommendation":"major_revision","confidential_remarks":"The paper has a substantial body of spectral and decay calculations that could be useful, but the S–D mixing section — the main advertised result — is internally inconsistent and circular as written. If the authors can correct the unit inconsistency in Table XXI, recompute the mixing predictions accordingly, and reframe the S–D analysis as a fit rather than an independent prediction, the paper may become publishable. In its current form, I would not recommend acceptance."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Quick take: this is a broad bottomonium catalog from a relativistic screened potential model, and the pure-spectrum part is a usable reference. But the S-D mixing section, which is the paper's headline, has a units problem that makes Table XXI internally inconsistent, and the abstract contradicts the body on Upsilon(10753). The mixing claims should not be taken as predictions until that is fixed.\n\nWhat is new: the authors extend their charmonium model to bottomonium and give masses for S through G states, leptonic/radiative/annihilation widths, and mixing angles for the four Upsilon resonances. The eigenvalue calculation is a genuine solution of the spinless Salpeter equation, and the comparison tables with several other models are useful. For lower states the numbers are reasonable; the chi_b(1P,2P) splittings and several E1 widths track experiment well.\n\nSoft spots: the central S-D mixing section. Eq. (35) and (36) do not reproduce Table XXI when the quoted mixing angles are read as degrees. With theta = 19.28 deg for 3S-2D, Eq. (35) gives about 10384 MeV, not 10374.9; with theta = 44.55 deg for 5S-4D, the mass split would be over a GeV. The table only works if the angles are in radians, but they are reported in degrees and compared with degree-valued angles from Ref. [61]. So the physical mixing amplitudes are misreported. On top of that, the abstract calls Upsilon(10753) a pure 3D state while the body and Table XXI treat it as the 4S-3D mixed partner. The mixing angles are also fit to the leptonic widths and then used to compute the mixed masses, so the agreement in Table XXII is not an independent test. Add to that deviations of 50-160 MeV in the pure S-wave masses, no uncertainty propagation, and the unstated potential constant V0, and the model's predictive weight rests mostly on the lower states.\n\nWho it's for: hadron spectroscopists who want a broad set of bottomonium numbers for comparison; the mass and transition tables are worth having even if the mixing interpretation is not yet established.\n\nRecommendation: send it to a competent referee, but with a direct request to check Table XXI. As written, the mixing section needs major revision; the rest is publishable after standard polishing.","headline":"A useful bottomonium spectrum from a relativistic screened model, but the S-D mixing table is internally inconsistent (degrees vs radians) and the abstract contradicts the body on Upsilon(10753), so the headline interpretation needs major repair.","tokens_in":31335,"tokens_out":3711,"would_cite":false,"duration_ms":33177,"reading_group":"maybe","serious_thinker":"no","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":["12.39.Pn","14.40.Pq"],"model":"deepseek-v4-flash","headline":"A relativistic screened potential model with S-D mixing assigns Υ(10355), Υ(10580), Υ(10860), and Υ(11020) — plus Υ(10753) — as conventional b-bbar states, reproducing their masses and di-leptonic widths.","keywords":["bottomonium","relativistic potential model","screened potential","S-D mixing","Upsilon resonances","leptonic decay widths","radiative transitions","quarkonium spectroscopy"],"falsifier":"Measure the $e^+e^-\\to\\Upsilon(10753)$ di-leptonic width and its radiative transitions: the model predicts $\\Gamma_{ee}\\approx 0.129$ keV and a D-dominated wavefunction, so a width near the S-wave scale or an E1 pattern characteristic of a pure 4S state would falsify the 4S-3D assignment. A lattice QCD calculation of the vector bottomonium spectrum above 10.5 GeV that includes $B^{(*)} \\bar B$ thresholds would directly test whether these resonances contain significant non-$b\\bar b$ components.","tokens_in":30303,"feed_emoji":"🔄","tokens_out":14015,"duration_ms":110489,"temperature":0.7,"pith_summary":"This paper argues that five excited Upsilon resonances, several of which have been proposed as hybrids or tetraquarks, are in fact ordinary bottomonium (b-bbar) states once S-D mixing is included. Using a relativistic screened potential model, the authors compute the mass spectrum and decay properties and fix the two-state mixing angle by matching the measured di-leptonic widths. They assign Υ(10355) to 3S-2D, Υ(10580) and Υ(10753) to the two members of the 4S-3D pair, and Υ(10860) and Υ(11020) to the 5S-4D pair; the resulting masses and lepton-pair widths agree with experiment to within a few tens of MeV and a few percent, respectively. If this interpretation is correct, the conventional quark model, with screening and S-D mixing, is sufficient for these states, and the paper's extensive tables of radiative and annihilation widths become a practical guide for searching the still-unseen higher states.","feed_headline":"S-D mixing makes five Upsilon peaks ordinary bottomonium","feed_subtitle":"The model matches the masses and lepton-pair widths of the excited Upsilon states — no exotic matter required.","key_machinery":"The engine of the calculation is the spinless Salpeter equation with a screened Coulomb-plus-linear potential, $V(r) = -\\frac{4}{3}\\frac{\\alpha_s(r)}{r} + \\lambda\\frac{1-e^{-\\mu r}}{\\mu} + V_0$, whose screening parameter $\\mu$ flattens the confining term at large distances. The wavefunctions are obtained by expanding in spherical Bessel functions on a finite interval and solving the resulting matrix eigenvalue equation, with spin-dependent corrections added perturbatively. On top of this, the S-D mixing ansatz $|\\phi\\rangle = \\cos\\theta\\,|nS\\rangle + \\sin\\theta\\,|n'D\\rangle$ (and its orthogonal partner) converts the pure-state masses and radial wavefunctions into the physical masses and leptonic widths; the angle $\\theta$ is fixed by matching the measured $\\Gamma_{ee}$ values. This two-state mixing, not the potential alone, carries the central claim.","core_discovery":"On its own terms, the paper's central result is that the persistent overestimation of the excited vector bottomonium masses by potential models disappears when each physical resonance is written as a linear combination of one S-wave and one D-wave component, $|\\phi\\rangle = \\cos\\theta\\,|nS\\rangle + \\sin\\theta\\,|n'D\\rangle$. The mixing angle $\\theta$ is not chosen to fit masses: it is fixed by the measured di-leptonic widths $\\Gamma_{ee}$, and the mixed-state masses then fall much closer to experiment than the pure states do. The assignments are $\\Upsilon(10355)$ = 3S-2D, $\\Upsilon(10580)$ = the lower 4S-3D state, $\\Upsilon(10753)$ = the higher, D-dominated 4S-3D partner, and $\\Upsilon(10860)$/$\\Upsilon(11020)$ = the 5S-4D pair. The paper concludes that these resonances are conventional $b\\bar b$ states and that the S-D admixture, rather than exotic degrees of freedom, accounts for their masses and widths.","pith_inferences":["A decisive check is the di-leptonic width of Υ(10753), predicted here at 0.129 keV; because the state is mostly 3D, an experimental $\\Gamma_{ee}$ near the S-wave scale would rule out the 4S-3D assignment.","The mixing angle is inferred from $\\Gamma_{ee}$ alone; a coupled-channel calculation with explicit $B^{(*)} \\bar B$ thresholds could reveal whether the residual mass gaps of tens of MeV are absorbed by $\\theta$ or require further components.","The same screened-potential plus S-D mixing scheme, already applied to charmonium, should predict which higher states mix most strongly as a function of the screening parameter $\\mu$; that is testable in lattice or experimental data.","The model's radiative-transition predictions, such as $\\Upsilon(10753) \\to \\gamma \\chi_b(3P)$, could be searched for in Belle II data and would distinguish a D-dominated from an S-dominated wavefunction."],"forward_implications":["If the assignments hold, Υ(10355), Υ(10580), Υ(10753), Υ(10860), and Υ(11020) are all conventional $b\\bar b$ states, so hybrid and tetraquark explanations are not required for them.","The fitted di-leptonic widths match experiment closely: 0.440 keV vs 0.443 keV for Υ(10355), 0.272 vs 0.272 keV for Υ(10580), and 0.291 vs 0.31 keV and 0.142 vs 0.13 keV for the 5S-4D pair.","S-D mixing resolves part of the systematic overestimation of the pure 4S and 5S masses: the lower mixed state is pulled down by its D-wave admixture.","The E1 and M1 transition tables and annihilation widths give concrete guidance for Belle II, LHCb, and PANDA searches for the unobserved higher states such as 2D, 3P, F, and G waves.","For n ≥ 3 the di-leptonic width of the D-dominated partner becomes sensitive to even a small mixing angle, so the predicted $\\Gamma_{ee}$ values distinguish the S-D assignment from a pure-state interpretation."],"supporting_citations":[{"why":"Supplies the S-D mixing ansatz (two-state linear combinations and the mass/leptonic-width formulas) and the earlier prediction that Υ(10860) and Υ(11020) are 5S-4D mixtures.","marker":"[61]"},{"why":"Provides the experimental masses, di-leptonic widths, and branching ratios that the fit targets and compares against.","marker":"[67]"},{"why":"The higher-bottomonium 'zoo' paper whose mass and decay predictions serve as the main baseline this model improves.","marker":"[46]"},{"why":"Introduces the screened Q-Qbar potential (Coulomb plus $\\lambda(1-e^{-\\mu r})/\\mu$) that is the model's defining interaction.","marker":"[38]"},{"why":"Supplies the spectral-expansion method (spherical Bessel basis on a finite interval) used to solve the spinless Salpeter equation.","marker":"[39]"},{"why":"The authors' previous application of the same model to charmonium, establishing the framework's reliability before this bottomonium study.","marker":"[36]"},{"why":"The Van Royen-Weisskopf formula used to compute decay constants and the di-leptonic widths that fix the mixing angles.","marker":"[47]"},{"why":"The 3P0-model analysis arguing Υ(10580) has a significant meson-meson component, motivating S-D mixing above the $B^{(*)} \\bar B$ threshold.","marker":"[68]"}],"fun_headline_variants":["S-D mixing solves Upsilon mass puzzle without exotics","Bottomonium excited states explained by S-D admixture","Five Upsilon peaks are ordinary S-D mixed bottomonium","Relativistic screened potential fits Upsilon masses via S-D mixing","Upsilon(10355) to (11020) as conventional S-D states"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"Each physical Upsilon above threshold is exactly one S-wave plus one D-wave $b\\bar b$ component, with all coupled-channel, meson-loop, hybrid, and tetraquark effects ignored or absorbed into a single fitted angle $\\theta$; if any of those extra components are significant, the state assignments in Table XXII do not follow.","fun_headline_variants_meta":{"raw":{"variants":["S-D mixing solves Upsilon mass puzzle without exotics","Bottomonium excited states explained by S-D admixture","Five Upsilon peaks are ordinary S-D mixed bottomonium","Relativistic screened potential fits Upsilon masses via S-D mixing","Upsilon(10355) to (11020) as conventional S-D states"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000159,"raw_usage":{"total_tokens":1229,"prompt_tokens":948,"completion_tokens":281,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":564,"completion_tokens_details":{"reasoning_tokens":193}},"tokens_in":564,"tokens_out":281,"duration_ms":82517,"temperature":1.0,"reasoning_tokens":193,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-10T21:52:59.146785+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Measure the $e^+e^-\\to\\Upsilon(10753)$ di-leptonic width and its radiative transitions: the model predicts $\\Gamma_{ee}\\approx 0.129$ keV and a D-dominated wavefunction, so a width near the S-wave scale or an E1 pattern characteristic of a pure 4S state would falsify the 4S-3D assignment. A lattice QCD calculation of the vector bottomonium spectrum above 10.5 GeV that includes $B^{(*)} \\bar B$ thresholds would directly test whether these resonances contain significant non-$b\\bar b$ components.","supporting_citations":[{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Supplies the S-D mixing ansatz (two-state linear combinations and the mass/leptonic-width formulas) and the earlier prediction that Υ(10860) and Υ(11020) are 5S-4D mixtures."},{"cited_title":"Wang and X","cited_arxiv_id":null,"evidence_quote":"Provides the experimental masses, di-leptonic widths, and branching ratios that the fit targets and compares against."},{"cited_title":"Li and K.-T","cited_arxiv_id":null,"evidence_quote":"The higher-bottomonium 'zoo' paper whose mass and decay predictions serve as the main baseline this model improves."},{"cited_title":"Godfrey and N","cited_arxiv_id":null,"evidence_quote":"Introduces the screened Q-Qbar potential (Coulomb plus $\\lambda(1-e^{-\\mu r})/\\mu$) that is the model's defining interaction."},{"cited_title":"Yibing, C","cited_arxiv_id":null,"evidence_quote":"Supplies the spectral-expansion method (spherical Bessel basis on a finite interval) used to solve the spinless Salpeter equation."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"The authors' previous application of the same model to charmonium, establishing the framework's reliability before this bottomonium study."},{"cited_title":"Wang, Z.-F","cited_arxiv_id":null,"evidence_quote":"The Van Royen-Weisskopf formula used to compute decay constants and the di-leptonic widths that fix the mixing angles."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"The 3P0-model analysis arguing Υ(10580) has a significant meson-meson component, motivating S-D mixing above the $B^{(*)} \\bar B$ threshold."}],"review_version":1}