{"id":"04697e59-e25f-4667-91ca-670d85c33cdb","arxiv_id":"2501.15110","paper_version":2,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":7,"one_line_summary":"A coupled-channel quark model reproduces bottomonium data and explains the Y(10753) bump as a B* anti-B* threshold effect, not a new resonance.","lead":"This paper models bottomonium, a bound state of a bottom quark and its antiquark, while allowing it to fluctuate into pairs of B mesons. It predicts that the Y(10753) bump seen in electron-positron collisions is a threshold effect rather than a genuine new particle, and it gives a full spectrum of higher bottomonium states to search for.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The claimed Y(10753) threshold-effect bump is not present in Eq. (41), an incoherent sum of constant-width Breit-Wigners; the threshold effect is asserted, not computed.","rationale":"The reader's weakest assumption—neglect of continuum-continuum interactions—is a valid model limitation, but it is explicitly acknowledged by the authors and does not by itself invalidate the reported calculation. The more load-bearing issue is that the central claim about Y(10753) is extracted from a cross-section formula (Eq. 41) that cannot produce threshold cusps. The paper's Fig. 4 and Sec. VII describe bumps 'due to threshold effects', yet Eq. (41) is a sum of constant-width Breit-Wigner shapes, which contain no threshold singularities. Therefore the model, as presented, does not actually compute what the abstract claims. This is an internal inconsistency rather than a disagreement with external models, so it goes to the correctness risk of the headline result. The manuscript otherwise gives a systematic and mostly self-consistent treatment of masses, widths, S-D mixing, and dielectron decays, with reasonable agreement for several well-established states. The recommended CONDITIONAL verdict stands, but the conditions must include recomputing the cross section from the spectral density functions or otherwise demonstrating that the threshold effect is genuinely present in the model output. The reader's FSI concern remains worth pursuing, but it is secondary to the fact that the reported calculation does not contain the claimed effect.","tokens_in":48649,"tokens_out":10010,"duration_ms":92322,"concrete_test":"Replace the constant-width Breit-Wigner sum in Eq. (41) with the spectral density functions ωR(M) computed from the coupled-channel propagator (Eqs. 22-24, using g(M) from Eq. 13), and plot the resulting e+e-→B*B* cross section. If the bump at 10.71 GeV disappears or is a pure Lorentzian centered at the Υ1(3D) mass (10670 MeV), the threshold-effect interpretation is refuted; if a cusp-like asymmetry appears, the claim is supported.","verdict_should_be":"CONDITIONAL","load_bearing_attack":"The paper's central claim is that the Y(10753) bump arises from threshold effects of the B*B* channel coupling to Υ(4S). The only cross-section calculation reported in Sec. VII is Eq. (41): an incoherent sum over vector resonances of constant-width Breit-Wigner shapes, with no energy-dependent widths, no interference, and no explicit threshold singularities. Such a sum produces Lorentzian peaks at each resonance mass, possibly overlapping tails, but it cannot generate a cusp at a threshold. The authors nonetheless attribute the bump near 10.71 GeV to a 'threshold effect' and conclude Y(10753) is not a genuine state. This is an internal inconsistency: the conclusion is not a consequence of the calculation but an interpretation imposed on it. The acknowledged neglect of continuum-continuum interactions (Sec. II.B) is a separate model uncertainty; even if FSI were included, the current calculation as written does not contain the claimed effect. Thus the central claim is unsupported by the reported computation.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper develops an unquenched quark model for bottomonium in which bare quark-model states are coupled to two-body open-bottom meson channels via a 3P0 hadronic-loop vertex, with a once-subtracted self-energy and a Gaussian high-momentum cutoff. The model parameters are fitted to the low-lying bottomonium masses (chi2 = 61.5 for 14 states) and to the width of Upsilon(10580). Using this framework the authors compute mass shifts, bbbar-core probabilities, strong decay widths, hadronic-loop-induced S-D mixing, dielectron widths, and an e+e- -> bbbar cross section. They conclude that most high-lying states have large continuum components, that Upsilon1(3D,5D,6D) mix strongly with nearby S-wave states, and that the Upsilon(10753) bump may be a B*B* threshold effect rather than a genuine resonance.","tokens_in":48933,"tokens_out":3754,"duration_ms":43185,"significance":"If the calculations were fully substantiated, the paper would provide a broad, unified description of bottomonium above the open-bottom threshold, including masses, widths, mixing, and line shapes, with falsifiable predictions for missing 5P, 6P, 5D, and 6D states. The systematic inclusion of excited B-meson channels and the use of the same 3P0 vertex for both real and imaginary parts of the self-energy are valuable strengths, as is the effort to connect the model to Belle II observations of e+e- -> B(*)B(*) cross sections. However, the central interpretive claim about Upsilon(10753) is currently not a consequence of the reported cross-section calculation, and the absence of any propagated theoretical uncertainties makes it difficult to assess the precision of the many numerical predictions.","major_comments":[{"comment":"The cross-section calculation in Eq. (41) is an incoherent sum of constant-width Breit-Wigner resonances and therefore cannot generate threshold cusps or threshold-induced bumps. No energy-dependent widths, no threshold phase-space factors, and no self-energy line shapes from Sec. III enter Eq. (41). The two structures near 10.63 and 10.71 GeV shown in Fig. 4 are merely overlapping tails of constant-width Lorentzians, yet the text attributes them to 'threshold effects' of the Bbar B* and B* B* channels and uses this to claim that Upsilon(10753) may arise from the B* B* threshold. As written, the abstract and Sec. VII conclusions are not supported by the computation. The authors should instead evaluate the cross section using the model's own energy-dependent spectral functions or coupled-channel amplitudes, for example by including the complex self-energy of Eq. (13)-(16) in the resonance propagators; otherwise the threshold-effect interpretation should be removed or explicitly labeled as an inference rather than a result of the calculation.","section":"Sec. VII, Eq. (41)"},{"comment":"The neglect of continuum-continuum interactions is acknowledged, but it is load-bearing for the threshold interpretation. Near the B* B* threshold, where the claimed Upsilon(10753) effect resides, meson-meson final-state interactions can alter the cusp shape, shift the apparent peak, and change the extracted widths. Since the manuscript uses the absence of a genuine Upsilon(10753) state as a conclusion, it needs at least an estimate of the FSI uncertainty (e.g., by comparing with unitarized or K-matrix approaches such as Ref. [68]) or a clear statement that the threshold interpretation is model-dependent and could be modified by FSI. This issue is separate from the missing energy dependence in Eq. (41), but both must be addressed before the central claim can be accepted.","section":"Sec. II.B, Eq. (8)"},{"comment":"No theoretical uncertainties are propagated, although the model has seven parameters (alpha_s, sigma, C0, rc, gamma, Lambda, alpha'_s) and several procedural choices (the once-subtracted zero point, the 100 MeV virtual-channel window, the integration interval Delta in Eq. (25)). Many predictions are quoted to the MeV or sub-keV level, and comparisons such as the 30 MeV discrepancy for Upsilon(10860) or the claimed agreement for Upsilon(11020) cannot be evaluated without an error estimate. At minimum, the authors should quantify the sensitivity of the central results to variations of Lambda and rc within their stated ranges (0.8 +/- 0.2 GeV and the fitted rc) and to the choice of M0; otherwise the precision implied by the tables is not justified.","section":"Sec. III, Tables II-XI"}],"minor_comments":[{"comment":"The fit yields chi2 = 61.5 for 14 states with Merr = 5 MeV, but the paper does not discuss the quality of this fit or the number of degrees of freedom; a sentence stating the corresponding rms deviation would help readers calibrate the model uncertainty.","section":"Sec. II.C, Table I"},{"comment":"Entries such as '33D3 1-- 10702 -52/+16 10650/10718 30/60' are ambiguous; the caption should explicitly state that the two numbers correspond to the two physical solutions and which quantity (mass shift, width, or core probability) is being split.","section":"Sec. III.C, Table II"},{"comment":"The off-diagonal terms are written as DeltaM_SD(M) and DeltaM_DS(M) but the text immediately below writes 'DeltaMSD(M) and DeltaMSD(M)'; please use consistent notation and state whether the matrix is symmetric.","section":"Sec. V, Eq. (32)"},{"comment":"Equation (37) uses alpha'_s for the QCD radiative correction while Eqs. (39) and (40) use alpha_s; clarify whether these are the same quantity or intended to be different, and define both.","section":"Sec. VI, Eqs. (37) and (39)"},{"comment":"The phrase 'The interpolations of Upsilon(10753) as a hybrid or a tetraquark state are also possible' should read 'interpretations' or 'explanations'.","section":"Sec. VII"},{"comment":"In the manuscript version provided, the figure captions and axis labels contain unreadable encoded text (e.g., '/s57/s46/s52'); the figures must be regenerated with properly rendered labels and captions before publication.","section":"Captions of Figs. 5-8"}],"recommendation":"major_revision","confidential_remarks":"The paper is a solid extension of the authors' earlier charmonium and heavy-light unquenched model work, and many of the spectral and width predictions are valuable. The main problem is that the flagship Upsilon(10753) threshold-effect claim is not actually computed in Sec. VII: Eq. (41) is an incoherent Breit-Wigner sum, and the energy-dependent self-energy developed in Secs. II-III is not used in the cross section. This is fixable within the scope of the paper by recalculating the cross section with the model's own line shapes, but until then the abstract's strongest claim should not be stated as a model result. I would not reject the manuscript, because the bulk of the work is defensible and the missing calculation is clearly identified; I recommend major revision."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"You should know two things about arXiv:2501.15110. First, it is a serious, mostly self-consistent unquenched quark model calculation for bottomonium up to 11.3 GeV, including excited open-bottom channels. It produces useful predictions: mass shifts of tens of MeV, significant non-b bbar components, two-pole structures for 3D and 5D states, hadronic-loop S-D mixing angles around -14 to -18 degrees, and suppressed dielectron widths. The low-lying spectrum is described reasonably (chi2 = 61.5 for 14 states with 5 MeV scale errors), the Y(11020) width comes out 16.5 MeV versus 24(+8,-6) MeV, and the lower Y dielectron widths agree well. That part is worth taking seriously.\\n\\nNow the soft spots. The framework is not new—it was developed for charmonium and heavy-light mesons in refs [55,69]—and this application carries the same caveats: no propagated parameter uncertainties, hand-chosen cutoffs (Lambda = 0.78 GeV, rc = 0.101 fm), and neglect of continuum-continuum interactions. The Y(10860) width is a clear failure: predicted 8 MeV versus 37 +/- 4 MeV, and the paper resorts to speculative channels to patch it. Those are ordinary phenomenological limitations.\\n\\nThe real problem is the cross-section section. Equation (41) is an incoherent sum of constant-width Breit-Wigners. There are no energy-dependent widths and no explicit threshold singularities. A sum of Breit-Wigners cannot produce a cusp at the B*B* threshold. Yet the text says the bumps near 10.63 and 10.71 GeV are threshold effects and uses the latter to explain Y(10753) as not a genuine state. That conclusion is not a consequence of the reported calculation; the bump in Fig. 4 is just overlapping tails of the vector resonances. Threshold effects do appear in the self-energy mass-shift functions (the dips in Fig. 8), which generate the two-pole structures, but those are not what goes into Eq. (41). So the abstract's claim is stronger than what is computed. This is fixable—one should compute the cross section from the spectral densities with energy-dependent widths, or use a proper coupled-channel unitarized amplitude—but as written it is an internal inconsistency.\\n\\nWho is this for? Hadron spectroscopists, especially people working on bottomonium exotics and Belle II final states. There are concrete search predictions for 4P and 5D states and their decay channels that are worth having on the record. It deserves a serious referee because the spectroscopy is solid enough to be useful even if the Y(10753) interpretation needs reworking. I would send it to review and ask for a proper cross-section calculation and uncertainty estimates.","headline":"Solid systematic unquenched bottomonium calculation, but the headline Y(10753) threshold explanation is not actually computed in the reported cross-section formula.","tokens_in":49456,"tokens_out":4066,"would_cite":true,"duration_ms":37714,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":[],"model":"deepseek-v4-flash","headline":"An unquenched quark model that lets bottomonium states mix with open-bottom meson pairs reproduces the observed spectrum and predicts that the $\\Upsilon(10753)$ bump in $e^+e^-\\to b\\bar b$ is a threshold effect of the $\\bar B^*B^*$…","keywords":["bottomonium","unquenched quark model","coupled-channel effects","3P0 hadronic loops","threshold effects","Y(10753)","S-D mixing","dielectron decay widths"],"falsifier":"A high-statistics measurement of $e^+e^-\\to \\bar B^*B^*$ across $\\sqrt{s}=10.70$--$10.75$ GeV would settle it: if the enhancement at 10.71 GeV is a genuine Breit-Wigner resonance, its pole position and width will remain stable when the $B^*B^*$ threshold is treated with different final-state-interaction assumptions, whereas a threshold cusp will change shape and shift as those assumptions vary. A direct check of the $\\Upsilon_1(3D)$ interpretation is the predicted $\\Gamma_{ee}\\simeq 0.028$ keV and $\\sim 92\\%$ $\\bar B^*B^*$ branching; an experimental upper bound well below either would falsify the model's assignment.","tokens_in":48420,"feed_emoji":"⚛️","tokens_out":9427,"duration_ms":79677,"temperature":0.7,"pith_summary":"The paper tries to establish that the bottomonium spectrum, including the hard-to-classify resonances above the open-bottom threshold, can be described by an unquenched quark model: a conventional $b\\bar b$ core that is allowed to mix with the two-meson continuum through $^3P_0$ quark-pair-creation loops. Using all OZI-allowed channels near each bare state and a once-subtracted self-energy for distant channels, the authors reproduce the masses and widths of the well-established states and give predictions for the $S$-, $P$-, and $D$-wave spectrum up to about 11.3 GeV. The outcome is a set of concrete claims: coupled-channel effects shift high-lying masses by tens of MeV, most high-lying resonances carry large non-$b\\bar b$ components, hadronic loops induce $S$-$D$ mixing angles of roughly $-14^\\circ$ to $-18^\\circ$ for the D-wave vector states, and the non-$b\\bar b$ components suppress dielectron widths. Their most striking claim is that $\\Upsilon(10753)$ need not be a genuine state at all: the bump seen near 10.71 GeV in $e^+e^-\\to b\\bar b$ can arise as a threshold effect of the strong $\\Upsilon(4S)$--$\\bar B^*B^*$ coupling.","feed_headline":"Y(10753) is a threshold bump, not a particle, model says","feed_subtitle":"Bottomonium masses, decay widths, and the 10.71 GeV bump all fall out of one quark model with meson-pair loops.","key_machinery":"The load-bearing object is the hadronic self-energy $g(M)=\\sum_{BC} \\int |\\mathcal{M}_{A\\to BC}(p)|^2/(M-E_{BC}+i\\epsilon)\\,p^2 dp$, built from $^3P_0$ model transition amplitudes between the bare bottomonium state and two-meson channels. Its real part, computed with a once-subtracted dispersion integral, gives the mass shift; its imaginary part gives each partial width. A momentum cutoff $\\Lambda=0.78$ GeV damps unphysical high-momentum loop contributions. For resonances, the spectral density function $\\omega_R(M)$ supplies the $\\bar b\\bar b$ core probability, and the off-diagonal self-energy between $n^3S_1$ and $m^3D_1$ states, computed from the same loops, produces the $S$-$D$ mixing angles without free mixing parameters.","core_discovery":"The central claim is that the unquenched coupled-channel mechanism simultaneously explains the masses and strong widths of the known $\\Upsilon$ states and predicts the missing higher states. Mass shifts from hadronic loops are tens of MeV and non-monotonic: $\\Upsilon(4S)$ and $\\Upsilon(6S)$ shift down by about $-40$ and $-56$ MeV, bringing them close to $\\Upsilon(10580)$ and $\\Upsilon(11020)$, while $\\Upsilon(5S)$ shifts by only about $-2$ MeV. The model yields $\\bar b\\bar b$ core fractions as low as roughly 30 percent for the lower $\\Upsilon_3(3D)$ solution and about 37 percent for $\\Upsilon(11020)$; gives hadronic-loop-induced mixing angles $\\phi_D\\simeq -(14\\text{--}18)^\\circ$ for $\\Upsilon_1(3D,5D,6D)$; and produces two-pole structures for $\\Upsilon_3(3D)$, $\\Upsilon_3(5D)$, and $\\eta_{b2}(5D)$ from strong coupling to $\\bar B^*B^*$ and $\\bar B^*B_2^*(5747)$. In the computed $e^+e^-\\to b\\bar b$ cross section, the $\\bar BB^*$ and $\\bar B^*B^*$ thresholds generate bumps near 10.63 and 10.71 GeV, and the 10.71 GeV bump is the paper's explanation for $\\Upsilon(10753)$.","pith_inferences":["Neglect of continuum-continuum interactions is the main uncontrolled piece; including them could sharpen or shift the $\\bar B^*B^*$ cusp, so the $\\Upsilon(10753)$ line shape is also a probe of near-threshold meson-meson scattering.","The same once-subtracted loop machinery could be applied to other heavy-quark systems; if threshold artifacts explain $\\Upsilon(10753)$, analogous 'states' in charmonium or bottom-strange spectra deserve re-examination before exotic interpretations are adopted.","A direct test: measure the phase of the $e^+e^-\\to \\bar B^*B^*$ amplitude across 10.71 GeV; a threshold cusp produces a different phase motion than a resonance pole, so the data can distinguish the two explanations without model dependence."],"forward_implications":["If the model is right, $\\Upsilon(10753)$ is a threshold bump rather than a quark-antiquark state, so its line shape should not follow a stable Breit-Wigner form.","The predicted $\\Upsilon_1(3D)$ state at about 10.67 GeV, with $\\Gamma_{ee}\\simeq 0.028$ keV and about 92% branching to $\\bar B^*B^*$, should appear as an enhancement in $e^+e^-\\to \\bar B^*B^*$ near threshold.","$\\Upsilon(11020)$ should decay dominantly to $\\bar B B(1P_1)$ (about 75%), making that mode the best place to search for the missing excited $B(1P_1)$ meson.","The higher D-wave vectors $\\Upsilon_1(5D)$ and $\\Upsilon_1(6D)$ get dielectron widths enhanced to about 0.019 and 0.023 keV by loop-induced $S$-$D$ mixing, while the S-wave states' dielectron widths are suppressed by their continuum components.","The predicted masses, widths, and dominant channels for the $7S$, $5P$, $6P$, $5D$, and $6D$ states give concrete search targets in open-bottom final states."],"supporting_citations":[{"why":"Establishes the unquenched quark model machinery and once-subtracted self-energy treatment that this paper adapts from charmonium to bottomonium.","marker":"[55]"},{"why":"Supplies the heavy-light meson wave functions, masses, and parameter choices used for the final-state B and Bs mesons in decays and loops.","marker":"[69]"},{"why":"Provides the once-subtracted dispersion relation method that absorbs far-above-threshold channels into a redefined bare mass.","marker":"[70]"},{"why":"Introduces the 3P0 quark-pair-creation operator adopted as the bare-to-continuum interaction.","marker":"[85]"},{"why":"Gives the formula relating the derivative of the self-energy to the bare-state probability for bound states.","marker":"[98]"},{"why":"Supplies the spectral density function method used to estimate b-bbar core components of resonances above threshold.","marker":"[99]"},{"why":"A recent independent coupled-channels quark-model analysis that also attributes the 10.71 GeV structure to threshold effects and is used for comparison.","marker":"[113]"},{"why":"Provides the measured energy-dependent e+e- -> B*B* cross sections that motivate the predicted Y1(3D) signal near threshold.","marker":"[119]"},{"why":"Supplies the experimental masses, widths, and dielectron widths used to fix parameters and compare predictions.","marker":"[7]"}],"fun_headline_variants":["Upsilon(10753) is just a hadronic-loop bump","Y(10753) from meson-pair thresholds in quark model","Unquenched quark model predicts Y(10753) as threshold bump","Hadronic loops explain bottomonium masses and Y(10753)"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The load-bearing premise is that the two mesons in the continuum do not interact after being created, so $H_c$ contains only kinetic energy; near the $\\bar B^*B^*$ threshold, unmodeled meson-meson final-state interactions could alter the cusp that the paper identifies with $\\Upsilon(10753)$ and could shift the mass integrals that generate the two-pole structures.","fun_headline_variants_meta":{"raw":{"variants":["Upsilon(10753) is just a hadronic-loop bump","Y(10753) from meson-pair thresholds in quark model","Unquenched quark model predicts Y(10753) as threshold bump","Hadronic loops explain bottomonium masses and Y(10753)"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000644,"raw_usage":{"total_tokens":3080,"prompt_tokens":1182,"completion_tokens":1898,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":798,"completion_tokens_details":{"reasoning_tokens":1822}},"tokens_in":798,"tokens_out":1898,"duration_ms":11764,"temperature":1.0,"reasoning_tokens":1822,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-10T14:37:41.276556+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"A high-statistics measurement of $e^+e^-\\to \\bar B^*B^*$ across $\\sqrt{s}=10.70$--$10.75$ GeV would settle it: if the enhancement at 10.71 GeV is a genuine Breit-Wigner resonance, its pole position and width will remain stable when the $B^*B^*$ threshold is treated with different final-state-interaction assumptions, whereas a threshold cusp will change shape and shift as those assumptions vary. A direct check of the $\\Upsilon_1(3D)$ interpretation is the predicted $\\Gamma_{ee}\\simeq 0.028$ keV and $\\sim 92\\%$ $\\bar B^*B^*$ branching; an experimental upper bound well below either would falsify the model's assignment.","supporting_citations":[{"cited_title":"Vijande, F","cited_arxiv_id":null,"evidence_quote":"Supplies the heavy-light meson wave functions, masses, and parameter choices used for the final-state B and Bs mesons in decays and loops."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Provides the once-subtracted dispersion relation method that absorbs far-above-threshold channels into a redefined bare mass."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Introduces the 3P0 quark-pair-creation operator adopted as the bare-to-continuum interaction."},{"cited_title":"Barnes, S","cited_arxiv_id":null,"evidence_quote":"Gives the formula relating the derivative of the self-energy to the bare-state probability for bound states."},{"cited_title":"Weinberg, Elementary particle theory of composite p arti- cles, Phys","cited_arxiv_id":null,"evidence_quote":"A recent independent coupled-channels quark-model analysis that also attributes the 10.71 GeV structure to threshold effects and is used for comparison."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Provides the measured energy-dependent e+e- -> B*B* cross sections that motivate the predicted Y1(3D) signal near threshold."}],"review_version":1}