{"id":"0b6752a4-3fcd-4dfe-878a-fd5a104dc1f0","arxiv_id":"2505.02742","paper_version":2,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":6,"one_line_summary":"A coupled-channel meson-exchange calculation predicts P-wave B(*)Bbar(*) states Y(10600) and Y(10650), with Y(10650) matching threshold anomalies in Belle and Belle II data.","lead":"This paper predicts two new bottomoniumlike particles, Y(10600) and Y(10650), formed from P-wave interactions of bottom mesons near their thresholds. It argues that existing Belle and Belle II data already show signs of Y(10650), making it a testable candidate for a new class of exotic hadrons.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Flatté-style line-shape formula for the P-wave Y(10650) uses an S-wave self-energy (√ΔE ∼ k); the data-support comparison may therefore be invalid.","rationale":"The reader correctly identified the fragility of the data comparison: a single normalization point, no background, and no interference terms make the claimed Belle/Belle II support weak. My concern is more specific and more damaging to the paper's central claim as stated in the abstract, because it concerns the internal consistency of the line-shape formula used to produce Fig. 4. The formula is an S-wave Flatté-type propagator, but Y(10650) is asserted to be a P-wave pole. For a P-wave resonance, the open-channel width near threshold is proportional to k^3, and the self-energy should reflect that; using k instead changes the line shape, especially near threshold where the data point used for normalization sits. This is not a matter of taste or background modeling: it is a discrepancy between the claimed angular momentum and the equation used to fit the data. The theoretical pole existence, supported by the coupled-channel calculation and uncertainty scan, is not directly challenged by this concern; the issue is that the experimental 'support' for Y(10650) may be an artifact of an S-wave line-shape formula applied to a P-wave pole. A direct recomputation with the correct P-wave threshold behavior would settle the issue. If the corrected line shapes no longer reproduce the threshold enhancement and the B B̄* dip, the abstract's claim that data 'naturally account' for the pole would need to be retracted, while the prediction itself might survive as a model prediction. Because the paper's headline result combines the dynamical pole with experimental support, and the support hinges on this formula, the verdict should remain conditional, with the explicit additional condition that the line-shape formula be corrected and the comparison repeated. The reader's verdict is therefore not wrong, but it under-specifies the most concrete technical check needed.","tokens_in":15073,"tokens_out":10280,"duration_ms":117692,"concrete_test":"Recompute the Fig. 4 curves with the correct P-wave self-energy: replace √(−2μ_i(E−m_i_th)) by (−2μ_i(E−m_i_th))^{3/2} with the appropriate sign and reduced-mass prefactors, or, better, extract the line shape directly from the coupled-channel T-matrix of the same model. Refit B(Λ) to the 3.7-MeV B*B̄* data point and compare the predicted e+e−→B B̄* cross section to the Belle/Belle II data; if the dip position or magnitude shifts beyond the experimental uncertainties, the claimed data support rests on the incorrect S-wave formula.","verdict_should_be":"CONDITIONAL","load_bearing_attack":"In the section 'THE INDICATION OF Y(10650) IN THE OPEN-BOTTOM CROSS SECTIONS', the cross section is written as dσ_j/dE = B(Λ) μ_j g_j k_j(E)/(4π^2) divided by |E−E_Y − Σ_i μ_i/(8π^2) g_i √(−2μ_i(E−m_i_th))|^2. The self-energy term √(−2μ(E−m_th)) is the S-wave loop function, proportional to the c.m. momentum k. For the claimed P-wave B*B̄* pole with L=1, unitarity requires the width (imaginary part of the self-energy) to scale as k^{2L+1} = k^3, equivalently (−2μ(E−m_th))^{3/2}, not √ΔE. The residue g_i defined in Eq. (21) as |⟨k_R|V|ϕ⟩|^2 already contains k_R^2 for a P-wave, so the missing k^2 energy dependence cannot be absorbed into a redefinition of g_i without also altering the energy dependence away from the pole. Using the linear-k form changes the threshold line shape and the normalization at the single 3.7-MeV data point, and can alter the predicted B B̄* dip that is claimed to match Belle/Belle II data. The paper does not justify applying an S-wave Flatté formula to a P-wave pole; this is an internal-consistency concern, distinct from the background/interference question raised by the reader.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper studies B(*)\\bar B(*) scattering in a coupled-channel meson-exchange framework with complex scaling, reporting two isoscalar I(J^PC)=0(1^{--}) P-wave poles near the B\\bar B* and B*\\bar B* thresholds, named Y(10600) and Y(10650). The authors scan coupling constants and form factors to argue that both poles persist under theoretical uncertainties. They then construct a Flatté-like line-shape formula for e+e- -> B\\bar B* and e+e- -> B*\\bar B*, fix the production strength B(Λ) to a single data point of the B*\\bar B* threshold enhancement, and claim that the resulting prediction for a dip in e+e- -> B\\bar B* agrees with Belle and Belle II data. The abstract additionally promises a hidden-bottom calculation of Y(10650) -> \\Upsilon(2S)\\eta with a width estimate, but that calculation does not appear in the body of the manuscript.","tokens_in":15368,"tokens_out":9398,"duration_ms":106587,"significance":"If the two poles are confirmed, this would be the first evidence for P-wave dimeson states in the bottomonium sector, providing a strong test of the mechanism previously proposed for G(3900). A genuine strength is the systematic parameter scan: 1250 pole searches over four coupling-constant rescaling factors and a comparison of non-local and local form factors, which indicates that the near-threshold nature of Y(10650) is not fine-tuned. The paper also clearly separates the pole search from the data normalization, so the main cross-check is not circular. However, the quantitative experimental support is currently fragile: the line-shape comparison is built on a one-point fit without background, uncertainties, or goodness-of-fit, and, more seriously, the Flatté formula used for the P-wave pole has an S-wave threshold behavior. The advertised hidden-bottom analysis is absent, which further weakens the paper in its present form.","major_comments":[{"comment":"The line-shape formula for dσ_j/dE contains the self-energy term Σ_i (μ_i/8π^2) g_i sqrt(-2μ_i(E-m_i^th)). This is the nonrelativistic S-wave loop function, whose imaginary part is proportional to k_i. For the claimed P-wave B*\\bar B* pole, unitarity requires the imaginary part of the self-energy to scale as k_i^3, i.e. (E-m_i^th)^{3/2}, and the production vertex should carry an additional factor k^2. The residue g_i defined in Eq. (21) is evaluated at the pole position and cannot absorb this energy dependence without changing the line shape away from the pole. Since the fit anchors the normalization at the single point 3.7 MeV above the B*\\bar B* threshold and the predicted B\\bar B* dip depends on the interference between channels, the use of the S-wave Flatté form is an internal-consistency problem for the data comparison in Fig. 4. The authors should either replace the self-energy with the proper P-wave loop function or justify why the S-wave form is an adequate approximation in this energy region.","section":"THE INDICATION OF Y(10650) IN THE OPEN-BOTTOM CROSS SECTIONS"},{"comment":"The experimental support is built by fitting the production strength B(Λ) to exactly one data point, with no background model, no interference with non-resonant contributions, no propagated experimental uncertainties, and no goodness-of-fit measure. The text states: 'we fit the theoretical line shape of Y(10650) to match the enhancement data point at the threshold above 3.7 MeV in the B*\\bar B* channel, we then predict the associated cross section in the e+e- -> B\\bar B*.' The subsequent agreement is judged visually. In addition, Fig. 4 shows predictions only for Λ=0.45, 0.50, 0.55, while Fig. 3 and Tables I-II consider the wider range 0.40-0.60; the omission of Λ=0.4 and 0.6 is not explained. A quantitative fit including background and errors, or a clear statement that only a qualitative comparison is claimed, is needed before the abstract's 'support from Belle and Belle II data' can be accepted.","section":"THE INDICATION OF Y(10650) IN THE OPEN-BOTTOM CROSS SECTIONS and Fig. 4"},{"comment":"The abstract states: 'We further study the hidden-bottom transition Y(10650) -> \\Upsilon(2S)\\eta through a near-threshold B*\\bar B* loop mechanism. The resulting O(10~100 keV) width ... is sufficient to account for the corresponding cross sections measured by Belle II.' No such analysis appears anywhere in the body or the appendix. This is a load-bearing advertised result in the abstract's case for Y(10650) appearing in both open- and hidden-bottom channels. The manuscript must either include the hidden-bottom calculation with enough detail to be checked, or the abstract must be revised to remove the claim.","section":"Abstract and full text"}],"minor_comments":[{"comment":"The caption of both tables defines δE = E_pole - m^0_th with m^0_th = m_B + m_B*, but for Y(10650) the relevant threshold is m_B* + m_B*. The quoted Y(10650) δE values (a few MeV) are inconsistent with m_B + m_B* ≈ 10604 MeV and only make sense relative to the B*\\bar B* threshold. The caption should be corrected.","section":"Tables I and II"},{"comment":"The channel probabilities P_i defined by the c-product are complex and in some rows take unphysical negative values, e.g. Table II at Λ=0.4 gives P(\\bar B*B + c.c.) = -1.2-1.4i %. Calling these quantities 'channel probabilities' or 'roughly reflecting the ratio of the channel' is misleading; the authors should either use |P_i| or another positive-definite measure, or explicitly explain how to interpret complex and negative values.","section":"Tables I and II"},{"comment":"The y-axis labels in Fig. 4 appear garbled as 'σ d ressed' and the legend entry 'Hint from e+ e- -> b b  Belle II fit' is unclear. The figure should be cleaned up for readability.","section":"Figure 4"},{"comment":"The title in the body text, 'Emergence of new heavy quarkoniumlike states: Y (10600) and Y (10650)', differs from the abstract title 'Emergence of the exotic bottomoniumlike state Y(10650) and support from Belle and Belle II data'. The two should be aligned, especially since the body presents both Y(10600) and Y(10650) as the main predictions.","section":"Title and abstract"}],"recommendation":"major_revision","confidential_remarks":"The central pole calculation is interesting and the parameter scan is a commendable robustness check, but the advertised experimental support and the hidden-bottom analysis are not in a publishable state. The missing hidden-bottom section in particular suggests the submission may be a version mismatch; the editor may wish to check the intended version before review proceeds."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"The paper's real content is a prediction, not a discovery claim. In a coupled-channel meson-exchange framework, the authors find two P-wave J^PC=1^-- poles near the B Bbar* and B* Bbar* thresholds, which they call Y(10600) and Y(10650). That part is new, concrete, and worth taking seriously.\n\nWhat's good: the framework is the same one that produced G(3900) in the charm sector, with a track record and external support from two global fits. They do a serious uncertainty scan—1250 coupling combinations plus a local form-factor check—and the poles survive. That is diligent and reproducible in spirit even if no code is shipped.\n\nThe soft spots are mostly in the data-support section. First, and most worrying, the line-shape formula they use has a self-energy term proportional to sqrt(-2μ(E-m_th)), which is linear in the c.m. momentum k. That is the S-wave loop function. For a claimed P-wave resonance, unitarity requires the imaginary part to scale as k^{2L+1}=k^3. The residue g_i defined in Eq. (21) already contains k_R^2, so you cannot absorb the missing k^2 into a redefinition of g_i without altering the energy dependence. This is an internal inconsistency, not just an external objection. It changes the threshold line shape, the normalization at the single 3.7-MeV data point, and therefore the predicted dip in e+e- → B Bbar* that is claimed to match Belle/Belle II. The authors never justify using an S-wave form for a P-wave pole.\n\nSecond, the data comparison is built on one data point. B(Λ) is fixed by the 3.7-MeV enhancement; there is no background, no uncertainties, no fit quality. The dip looks plausible visually, but between the background question and the k-vs-k^3 issue, I'd call the experimental support suggestive at best.\n\nMinor: the abstract mentions only Y(10650) while the full text and title discuss Y(10600) as well. That looks like a version mismatch and should be cleaned up.\n\nBottom line: the pole prediction deserves attention—it gives Belle II a concrete measurement program—but the data-support claim does not hold as presented. I would not cite the data comparison as evidence; I would cite the pole prediction with a caveat. Send to peer review: a serious referee can push for a proper P-wave unitarized line shape and a background-inclusive fit, and the core result is important enough to warrant that revision.","headline":"A real and testable P-wave bottomonium prediction buried under an over-claimed data comparison and an S-wave line-shape mistake.","tokens_in":15979,"tokens_out":3333,"would_cite":true,"duration_ms":39564,"reading_group":"yes","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":[],"model":"deepseek-v4-flash","headline":"A near-threshold $P$-wave $B^*\\bar B^*$ pole, $Y(10650)$, explains Belle and Belle II open-bottom line shapes and would fill the $\\Upsilon(4S)$-$\\Upsilon(5S)$ gap.","keywords":["Y(10650)","P-wave dimeson","exotic bottomonium","complex scaling method","meson exchange model","B*Bbar* interaction","Belle II data","near-threshold resonance"],"falsifier":"A high-statistics measurement of $e^+e^-\\to B^*\\bar B^*$ and $e^+e^-\\to B\\bar B^*$ across the full 10.60--10.70 GeV region: if the observed bump and dip cannot both be reproduced by a single $Y(10650)$ state with the relative strength fixed by the pole residues, the claim is ruled out.","tokens_in":14836,"feed_emoji":"⚛️","tokens_out":14515,"duration_ms":148913,"temperature":0.7,"pith_summary":"This paper argues that a state called $Y(10650)$, sitting essentially at the threshold where a $B^*$ meson and its antiparticle can be produced, is not a conventional quark-antiquark meson but a $P$-wave \"dimeson\" object generated by meson exchange between two heavy mesons. Using a coupled-channel meson-exchange model solved with the complex scaling method, the authors find a $J^{PC}=1^{--}$ pole very close to the $B^*\\bar B^*$ threshold and show that its production line shape reproduces both the sharp enhancement just above that threshold in $e^+e^-\\to B^*\\bar B^*$ and the dip in $e^+e^-\\to B\\bar B^*$ seen by Belle and Belle II. If the claim is right, $Y(10650)$ becomes the first neutral isoscalar exotic bottomoniumlike state in the gap between $\\Upsilon(4S)$ and $\\Upsilon(5S)$, and it would establish that higher partial waves, normally suppressed by the centrifugal barrier, can still generate near-threshold exotic hadrons.","feed_headline":"Belle II bumps point to new exotic bottomonium state Y(10650)","feed_subtitle":"If real, a P-wave B* meson pair would fill the gap between Upsilon(4S) and Upsilon(5S).","key_machinery":"The central object is a coupled-channel, partial-wave projected potential for $B^{(*)}\\bar B^{(*)}$ scattering built from heavy meson effective Lagrangians, solved by the complex scaling method (CSM), which rotates coordinates into the complex plane so that both bound states and resonances appear as stable poles of the complex-scaled Schrodinger equation. The load-bearing identity is the pole expansion of the $T$-matrix near the pole: the cross section for $e^+e^-\\to Y\\to B^{(*)}\\bar B^{(*)}$ is expressed in terms of the pole residues $g_i$ (effective couplings of the $Y$ to each open-bottom channel) and phase-space factors, with a single overall production-strength parameter $B(\\Lambda)$. This machinery converts the dynamical pole into concrete line-shape predictions whose relative normalization across channels is fixed by the residues.","core_discovery":"The paper's central claim is that the $P$-wave $B^*\\bar B^*$ interaction generates a $J^{PC}=1^{--}$ (vector) pole, $Y(10650)$, within a few MeV of the $B^*\\bar B^*$ threshold, and that this pole is responsible for the observed threshold enhancement in $e^+e^-\\to B^*\\bar B^*$ and the associated dip in $e^+e^-\\to B\\bar B^*$. The pole is dynamically generated by coupled-channel meson exchange ($\\pi$, $\\eta$, $\\rho$, $\\omega$, $\\sigma$) rather than by a quark-antiquark pair, and the centrifugal barrier plays a constructive role: for weakly attractive potentials, the barrier turns a would-be bound state into a narrow near-threshold resonance. Once the production strength is fixed by the $B^*\\bar B^*$ enhancement data point, the pole residues fix the relative strength in $B\\bar B^*$, and the predicted dip location and magnitude match the Belle and Belle II measurements. A companion pole, $Y(10600)$, is also found near the $B\\bar B^*$ threshold and is identified as the bottom analogue of the charmoniumlike $G(3900)$.","pith_inferences":["If confirmed, the same mechanism suggests searching for analogous $P$-wave near-threshold states in other heavy-flavor systems, such as $B_s^{(*)}\\bar B_s^{(*)}$ or $B_c^{(*)}\\bar B^{(*)}$, where the reduced kinetic energy should make the effect even more pronounced.","The dip-peak correlation between the two open-bottom channels is a generic single-pole signature; a combined fit to both channels with an explicit background model could be turned into a model-independent pole-extraction test.","Because the production strength is anchored to a single data point, the full energy dependence of the cross sections provides a stricter test than the one shown; the predicted line shapes can be checked point-by-point with more Belle II data."],"forward_implications":["$Y(10650)$ would be the first neutral isoscalar exotic bottomoniumlike state in the $\\Upsilon(4S)$--$\\Upsilon(5S)$ gap, where conventional quark models predict no vector states.","The state would confirm that $P$-wave dimeson dynamics can produce near-threshold exotics despite the centrifugal barrier, supporting the $P$-wave interpretation of $G(3900)$.","The predicted hidden-bottom transition $Y(10650)\\to\\Upsilon(2S)\\eta$, with a width of order 10--100 keV, should be observable in Belle II cross-section data.","A companion pole $Y(10600)$ near the $B\\bar B^*$ threshold predicts a corresponding threshold enhancement around 10.60 GeV, giving a second search target."],"supporting_citations":[{"why":"Supplies the unified meson-exchange framework and the $P$-wave resonance mechanism that this paper extends to the bottom sector.","marker":"[15]"},{"why":"Provides the Belle data for $e^+e^-\\to B\\bar B^*$ that show a dip near the $B^*\\bar B^*$ threshold.","marker":"[19]"},{"why":"Provides the Belle II data for $e^+e^-\\to B^*\\bar B^*$ that show the threshold enhancement just above 10.65 GeV.","marker":"[20]"},{"why":"Provides the high-precision BESIII confirmation of $G(3900)$, the charm-sector state that motivates the search for bottom $P$-wave partners.","marker":"[14]"},{"why":"Provides the pole-expansion formula used to compute the cross-section line shapes from the pole residues.","marker":"[44]"},{"why":"Introduces the complex scaling method used to extract bound-state and resonance poles.","marker":"[31]"},{"why":"Companion to [31], establishing the complex scaling method's mathematical foundations.","marker":"[32]"},{"why":"Provides the inclusive $e^+e^-\\to b\\bar b$ data used to set the dip minimum in $B\\bar B^*$ to zero.","marker":"[21]"}],"fun_headline_variants":["P-wave B* pairs spawn exotic Y(10650)","Y(10650) emerges from P-wave B*Bbar*","Belle data support P-wave origin of Y(10650)","Exotic Y(10650) born from P-wave B*Bbar*","Belle and Belle II back Y(10650) P-wave origin"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The whole experimental comparison depends on the assumption that the measured cross sections are dominated by the single $Y(10650)$ state, with no important background from other processes.","fun_headline_variants_meta":{"raw":{"variants":["P-wave B* pairs spawn exotic Y(10650)","Y(10650) emerges from P-wave B*Bbar*","Belle data support P-wave origin of Y(10650)","Exotic Y(10650) born from P-wave B*Bbar*","Belle and Belle II back Y(10650) P-wave origin"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000983,"raw_usage":{"total_tokens":4289,"prompt_tokens":1178,"completion_tokens":3111,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":794,"completion_tokens_details":{"reasoning_tokens":3018}},"tokens_in":794,"tokens_out":3111,"duration_ms":28418,"temperature":1.0,"reasoning_tokens":3018,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-16T00:42:49.090604+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"A high-statistics measurement of $e^+e^-\\to B^*\\bar B^*$ and $e^+e^-\\to B\\bar B^*$ across the full 10.60--10.70 GeV region: if the observed bump and dip cannot both be reproduced by a single $Y(10650)$ state with the relative strength fixed by the pole residues, the claim is ruled out.","supporting_citations":[{"cited_title":"Double pole structures of $X_1(2900)$ as the $P$-wave $\\bar{D}^*K^*$ resonances","cited_arxiv_id":"2408.08965","evidence_quote":"Provides the Belle data for $e^+e^-\\to B\\bar B^*$ that show a dip near the $B^*\\bar B^*$ threshold."}],"review_version":1}