{"id":"f7847fff-0a8d-4676-87fb-dbce7e1db8e5","arxiv_id":"2501.04298","paper_version":1,"verdict":"CONDITIONAL","confidence":"HIGH","novelty_score":4.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":4,"one_line_summary":"A coupled-channel quark model identifies the newly observed B_s(6064) and B_s(6158) states as the 1^3D_3 and 1^3D_1 bottom-strange mesons.","lead":"This paper calculates the masses and decay rates of bottom-strange mesons, particles made of a bottom quark and a strange antiquark, using a quark model that lets each particle temporarily become two lighter mesons. It argues that two recently discovered particles, B_s(6064) and B_s(6158), are the D-wave states 1^3D_3 and 1^3D_1.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Eq. (3) mass-shift integral for the above-threshold D-wave states has no stated pole prescription, so the -104 and -135 MeV shifts that carry the B_sJ assignments are not defined as written.","rationale":"The reader's conditional verdict already flags the quantitative reliability of the coupled-channel self-energy, specifically the unstated principal-value treatment of Eq. (3). The stress test sharpens this into the single most load-bearing concern: the headline D-wave assignments are entirely mediated by large self-energy shifts, and Eq. (3) is formally undefined for states above threshold unless a pole prescription is specified. This is an internal correctness risk, not just parameter sensitivity, because different acceptable readings of Eq. (3) can change the above-threshold contributions and therefore the predicted masses by tens of MeV. The concern does not by itself falsify the assignments; the framework is standard, the fitted C_bs is disclosed, and the final numbers agree well with LHCb, so a conditional verdict remains appropriate. The paper needs a stated pole prescription and a sensitivity scan before the assignments can be accepted as robust. Since the reader's verdict is already CONDITIONAL, no verdict change is needed.","tokens_in":13508,"tokens_out":4172,"duration_ms":44511,"concrete_test":"Recompute the 1^3D_1 and 1^3D_3 self-energies from Eq. (3) with an explicitly stated principal-value convention (e.g., split the integral at M−ε and M+ε, or take the real part of the iε-regulated integral), keeping all other parameters and wave functions identical. If either predicted mass moves by more than about 20 MeV relative to Table III, or if the sign of the BK contribution in the 1^3D_3 row changes, then the claimed assignments are prescription-dependent and should be presented as parameter-dependent candidates. As a secondary check, repeat with γ0 = 0.3 and 0.5 and r_q = 0.25 and 0.35 fm to bound the resulting mass predictions.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central claim that B_sJ(6064) = 1^3D_3 and B_sJ(6158) = 1^3D_1 depends on the coupled-channel mass shifts in Table III: -135 MeV for 1^3D_3 and -104 MeV for 1^3D_1 bring the bare quark-model masses to within 1-16 MeV of the LHCb masses. These shifts are computed from Eq. (3), ΔM = Σ ∫ p²dp |⟨BC|H_I|ψ0⟩|² / (M - E_BC). For every D-wave candidate the physical mass lies above the BK threshold (M_B + M_K ≈ 5861 MeV; M_B* + M_K ≈ 5909 MeV), so the integrand has a pole at E_BC = M. The paper never states whether the integral is a principal value, is regulated by an iε prescription, or is evaluated with some other subtraction; Eq. (3) is therefore not a defined quantity for the states that carry the paper's headline. The ambiguity is not academic: Table IV gives +12 and +5 MeV from BK and B*K for 1^3D_1 (above threshold), while for 1^3D_3 the BK entry is -19 MeV even though M = 6079 MeV > 5861 MeV; signs and magnitudes of such above-threshold pieces are convention-dependent. In addition, M appears on both sides of Eq. (3), and the text does not state how the self-consistent equation is solved or whether multiple solutions were checked. A 20-30% change in γ0, or moving r_q across the cited 0.25-0.35 fm range, rescales these shifts, so without the pole prescription the agreement with 6063.5/6158 MeV cannot be assessed. This is an internal omission in the argument, not merely a disagreement with other models.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper calculates the mass spectrum of bottom-strange mesons using a nonrelativistic quark model augmented with coupled-channel effects, and computes strong decay widths with the 3P0 model using numerically obtained wave functions. The authors assign the established states B_s, B_s*, B_s1(5830), and B_s2*(5840) to 1^1S0, 1^3S1, 1P', and 1^3P2, respectively, and propose that the newly observed LHCb states B_sJ(6064) and B_sJ(6158) are the 1^3D3 and 1^3D1 bottom-strange mesons, based on predicted masses and widths that agree with the experimental values within tens of MeV.","tokens_in":13716,"tokens_out":2702,"duration_ms":26706,"significance":"If the central assignments are correct, the paper would pin down the quantum numbers of two recently observed, as-yet-unassigned B_s states, which is a useful step for the spectroscopy program at LHCb and future facilities. The work has concrete strengths: it provides a transparent parameter set (Table II), compares with many other models (Tables III), and reproduces nontrivial decay ratios such as the B_s2*(5840) branching fraction ratio in Eq. (17). The main limitation is that the quantitative reliability of the D-wave predictions hinges on a coupled-channel self-energy integral whose definition is incomplete for the above-threshold states, so the headline agreement with the LHCb masses cannot be fully assessed as written.","major_comments":[{"comment":"The mass-shift integral in Eq. (3) has a pole at E_BC = M for every D-wave state, since the physical masses (6157 MeV and 6079 MeV for 1^3D1 and 1^3D3) lie above the principal thresholds (e.g., BK at ~5861 MeV and B*K at ~5909 MeV). The paper never states whether the integral is a principal value, regulated by an iε prescription, or treated with some subtraction. Without that prescription, the -104 MeV and -135 MeV shifts in Table III that carry the central assignments of B_sJ(6158) and B_sJ(6064) are not defined quantities. The text must specify the prescription and demonstrate that the predictions are stable under reasonable alternative choices.","section":"Section II.A, Eq. (3)"},{"comment":"The constant C_bs = 0.169 GeV is determined by reproducing the ground-state B_s mass, as stated in the opening paragraph of Section III, so the B_s row in Table III is a fit by construction rather than a prediction; the same applies to the other absolute masses, which all shift with C_bs. More importantly, no uncertainty or sensitivity analysis is provided for the D-wave predictions. The mass shifts scale as gamma_0^2 and depend on the Gaussian regulator r_q = 0.3 fm, which is itself only a middle value of a 0.25–0.35 fm range. A 20–30% change in gamma_0 would alter the D-wave masses by tens of MeV, comparable to the differences between the predictions and data, so the robustness of the assignments needs to be quantified.","section":"Section III, Table III and Table II"},{"comment":"Equation (3) is a self-consistent equation because the physical mass M appears on both sides, but the manuscript does not state how this equation is solved or whether multiple solutions were checked. Since Table III reports large coupled-channel shifts for all states, the iteration scheme or algebraic solution method should be described, and the absence of alternative solutions should be confirmed.","section":"Section II.A, Eq. (3) and Section III"}],"minor_comments":[{"comment":"The word \"botton-strange\" in the sentence before Eq. (8) is a typo and should read \"bottom-strange.\"","section":"Section II.B"},{"comment":"In the description of LHCb results, \"finial state\" should be \"final state.\"","section":"Introduction"},{"comment":"The phrase \"the the Bs(nL), andBs(nL′)\" contains a duplicated article and a missing space; it should read \"the B_s(nL) and B_s(nL′) states.\"","section":"Section II.B, above Eq. (14)"},{"comment":"The caption cites \"RPP [66]\" while the text and Table I use Ref. [2] for the Review of Particle Physics; consolidating to a single RPP reference would avoid confusion.","section":"Table III caption"},{"comment":"The PACS number line is left blank; either provide relevant PACS codes or remove the line.","section":"Abstract and Section I"}],"recommendation":"major_revision","confidential_remarks":"The manuscript addresses a timely topic and the qualitative picture is plausible, but the missing pole prescription in Eq. (3) directly affects the numerical values that motivate the two new quantum-number assignments. This is fixable within the manuscript's scope by specifying the prescription and adding a sensitivity study, so I recommend major revision rather than rejection."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Quick take: this is a solid, workmanlike coupled-channel quark-model paper from a group that has done this before. What is genuinely new is the 1^3D_1 assignment for B_s(6158) and the decay/prediction tables (channel probabilities, 2S/2P widths, the ratio 2.0 for B_s(2^3S_1)). The framework is standard, the one free parameter C_bs is disclosed, and the numerical wave functions are a real step up from a pure SHO basis.\n\nThe soft spots are real but mostly not fatal. The one I would make the referee focus on is Eq. (3). The mass shift integral has a pole for every D-wave candidate (all lie above BK threshold), and the paper never states whether the integral is principal value, iε-regulated, or something else. Table IV shows +12 and +5 MeV from BK and B*K for 1^3D_1 even though the state is above threshold — those signs are prescription-dependent. And M appears on both sides of Eq. (3); the text doesn't say how the self-consistent equation is solved. Without a stated pole prescription, the -104 and -135 MeV shifts that carry the headline assignments are not defined as written. That's an internal omission, not a model disagreement.\n\nAlso worth fixing: no uncertainty or sensitivity analysis (the shifts scale as gamma_0^2 and depend on r_q); the predicted 1^3D_1 and 1D' are nearly degenerate (6157 vs 6154 MeV) with widths that don't clearly discriminate (62 vs 35 MeV given the 72±18±25 experimental width); and the Summary says B_s(2^1S_0) main decay is BK, contradicting Table VI which shows B*K (91 MeV) and no BK. Minor, but sloppy.\n\nThe ground-state B_s mass is fit by construction via C_bs, so Table III's B_s row is not a prediction — the authors are upfront about this, so I don't hold it against them.\n\nNet: the central idea holds up in outline, and the 1^3D_3 assignment actually echoes the authors' earlier screened-potential paper. The real value is the 1^3D_1 claim plus the width-ratio predictions. Send it to a serious referee, but the referee should demand a stated pole prescription and a sensitivity scan over gamma_0 and r_q before the numbers can be trusted.","headline":"Competent coupled-channel calculation whose headline D-wave assignments rest on an unstated pole prescription in Eq. (3); worth refereeing but needs a stated regularization and a sensitivity analysis.","tokens_in":14559,"tokens_out":2323,"would_cite":false,"duration_ms":20883,"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":"Coupled-channel calculations identify the newly observed bottom-strange mesons BsJ(6064) and BsJ(6158) as D-wave quark-antiquark states.","keywords":["bottom-strange mesons","coupled channel effects","nonrelativistic quark model","3P0 model","D-wave mesons","strong decay widths","mass shifts"],"falsifier":"Measure the spin-parity of $B_{sJ}(6064)$ through the angular distribution of its $B^+K^-$ decay: the paper's assignment requires $J^P = 3^+$, whereas the alternative $B_s(2^3S_1)$ assignment gives a different angular dependence, so the observed distribution would settle the central claim.","tokens_in":13025,"feed_emoji":"⚛️","tokens_out":10897,"duration_ms":95372,"temperature":0.7,"pith_summary":"The paper sets out to show that the known and newly observed bottom-strange mesons can be understood as ordinary quark-antiquark states once the coupling to two-meson decay channels is included. It identifies the two newest states, $B_{sJ}(6064)$ and $B_{sJ}(6158)$, with the spin-triplet D-wave states $B_s(1^3D_3)$ and $B_s(1^3D_1)$, because coupled-channel effects shift the bare quark-model masses down to within about 15 MeV of the measured values and produce strong decay widths close to the data. It also assigns $B_{s1}(5830)$ and $B_{s2}^*(5840)$ to the $1P'$ and $1^3P_2$ states, reproducing their narrow widths and measured branching ratio. If these assignments are correct, the quantum numbers of the new states are fixed and the spectroscopy of this family becomes predictive up to D wave.","feed_headline":"New bottom-strange states identified as D-wave quark pairs","feed_subtitle":"Coupled-channel loops shift the predicted masses and widths into agreement with the observed states, fixing their quantum numbers.","key_machinery":"The central object is the coupled-channel self-energy, an integral over the momentum of each virtual two-meson pair with a denominator $M - E_{BC}$, where $E_{BC}$ is the total energy of the two mesons. The transition matrix element comes from the ${}^3P_0$ pair-creation operator, in which a quark-antiquark pair with vacuum quantum numbers is created and the pair has a Gaussian form factor of size $r_q = 0.3$ fm and strength $\\gamma_0 = 0.4$. The same matrix element, evaluated on shell, gives the strong decay widths. What carries the argument is that these self-energy corrections are large, about $-100$ to $-150$ MeV, and channel-dependent, so they move the bare quark-model states to the measured positions while the realistic numerical wave functions control the decay widths.","core_discovery":"The central claim is that a nonrelativistic potential model corrected by two-meson continuum channels reproduces the bottom-strange meson spectrum, and that this resolves the quantum numbers of the newly observed states. The paper computes bare masses from a nonrelativistic Hamiltonian, then adds self-energy shifts from ten coupled channels ($BK$, $B^*K$, $BK^*$, $B^*K^*$, $B_s\\eta$, $B_s\\eta'$, $B_s^*\\eta$, $B_s^*\\eta'$, $B_s\\phi$, $B_s^*\\phi$). The resulting masses for $B_s(1^3D_3)$ and $B_s(1^3D_1)$ are 6079 and 6157 MeV, with ${}^3P_0$ decay widths of 24 and 62 MeV, compared with measured masses $6063.5 \\pm 1.2 \\pm 0.8$ and $6158 \\pm 4 \\pm 5$ MeV and widths $26 \\pm 4 \\pm 4$ and $72 \\pm 18 \\pm 25$ MeV. The same calculation assigns $B_{s1}(5830)$ to $B_s(1P')$ and $B_{s2}^*(5840)$ to $B_s(1^3P_2)$, gives the ground states $B_s$ and $B_s^*$ as $1^1S_0$ and $1^3S_1$, and predicts the 2S and remaining 2P and 1D states.","pith_inferences":["Beyond the paper: the mass-shift calculation could be re-run with different pair-creation strengths and regulator sizes to map how stable the D-wave assignments are; the present paper quotes a single parameter set.","A further consequence the authors do not spell out is that the same unquenching mechanism should shift the corresponding bottom meson states by comparable amounts, so checking the analogous $B(5970)$ region would extend the framework to a neighboring family."],"forward_implications":["The quantum numbers of $B_{sJ}(6064)$ and $B_{sJ}(6158)$ would be fixed as $3^+$ and $1^+$ if the assignments hold, removing the present ambiguity.","The dominant strong decay modes of both new states are predicted to be $BK$ and $B^*K$, with widths near 24 and 62 MeV, so their observation in these channels is a direct check.","The two unobserved $1D$ partners, $B_s(1D)$ and $B_s(1D')$, are predicted at about 6077 and 6154 MeV with widths 84 and 35 MeV, giving concrete search targets.","The $2S$ states $B_s(2^1S_0)$ and $B_s(2^3S_1)$ are predicted near 5949 and 5992 MeV, with the ratio $\\Gamma(B_s(2^3S_1)\\to B^*K)/\\Gamma(B_s(2^3S_1)\\to BK)$ around 2.0.","The assignments of $B_{s1}(5830)$ to $B_s(1P')$ and $B_{s2}^*(5840)$ to $B_s(1^3P_2)$ fix the $1P$ mixing angle at about $-55.8$ degrees."],"supporting_citations":[{"why":"Provides the experimental masses and widths of the $B_{sJ}(6064)$ and $B_{sJ}(6158)$ states that the D-wave assignments must match.","marker":"[8]"},{"why":"Supplies the nonrelativistic quark model, the P-wave mixing treatment, and earlier assignments that the coupled-channel calculation starts from.","marker":"[14]"},{"why":"Establishes the coupled-channel framework with realistic wave functions used for the mass shifts.","marker":"[22]"},{"why":"Gives the ${}^3P_0$ transition operator and the self-energy formalism for meson-meson continuum coupling.","marker":"[31]"},{"why":"Applies the same coupled-channel method to bottomonium, providing the operational form of the mass-shift calculation.","marker":"[33]"},{"why":"Provides the constituent quark model whose parameters are used for the bare $B_s$ masses.","marker":"[45]"},{"why":"Supplies additional parameter validation for the nonrelativistic quark model in the heavy-light meson sector.","marker":"[46]"},{"why":"Determines the range for the Gaussian regulator $r_q$, whose midpoint is used to regulate the loop integrals.","marker":"[55]"}],"fun_headline_variants":["Coupled channels identify B_s(6064) and B_s(6158) as D-waves","Meson loops pin down D-wave nature of bottom-strange states","Coupled-channel shifts fix B_s(6064) and B_s(6158) identities","Bottom-strange mesons get D-wave tags from coupled channels","D-wave assignments for B_s(6064) and B_s(6158) from meson loops"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The quantitative reliability of the coupled-channel mass shift formula is load-bearing: its D-wave predictions are large shifts, about $-104$ and $-135$ MeV, that depend on the pair-creation strength $\\gamma_0 = 0.4$, the Gaussian regulator $r_q = 0.3$ fm, and an unspecified treatment of the pole when $M = E_{BC}$, so if any of those choices changes substantially, the match to the observed masses weakens.","fun_headline_variants_meta":{"raw":{"variants":["Coupled channels identify B_s(6064) and B_s(6158) as D-waves","Meson loops pin down D-wave nature of bottom-strange states","Coupled-channel shifts fix B_s(6064) and B_s(6158) identities","Bottom-strange mesons get D-wave tags from coupled channels","D-wave assignments for B_s(6064) and B_s(6158) from meson loops"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.001212,"raw_usage":{"total_tokens":5047,"prompt_tokens":1060,"completion_tokens":3987,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":676,"completion_tokens_details":{"reasoning_tokens":3875}},"tokens_in":676,"tokens_out":3987,"duration_ms":29339,"temperature":1.0,"reasoning_tokens":3875,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-10T21:37:38.975280+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Measure the spin-parity of $B_{sJ}(6064)$ through the angular distribution of its $B^+K^-$ decay: the paper's assignment requires $J^P = 3^+$, whereas the alternative $B_s(2^3S_1)$ assignment gives a different angular dependence, so the observed distribution would settle the central claim.","supporting_citations":[{"cited_title":"Observation of new excited B0 s states","cited_arxiv_id":null,"evidence_quote":"Provides the experimental masses and widths of the $B_{sJ}(6064)$ and $B_{sJ}(6158)$ states that the D-wave assignments must match."},{"cited_title":"Excited bottom and bottom-strange mesons in the quark model","cited_arxiv_id":null,"evidence_quote":"Supplies the nonrelativistic quark model, the P-wave mixing treatment, and earlier assignments that the coupled-channel calculation starts from."},{"cited_title":"Coupled-Channel Effects for the Bottomonium with Re- alistic Wave Functions","cited_arxiv_id":null,"evidence_quote":"Establishes the coupled-channel framework with realistic wave functions used for the mass shifts."},{"cited_title":"Ferretti, G","cited_arxiv_id":null,"evidence_quote":"Gives the ${}^3P_0$ transition operator and the self-energy formalism for meson-meson continuum coupling."},{"cited_title":"Ferretti and E","cited_arxiv_id":null,"evidence_quote":"Applies the same coupled-channel method to bottomonium, providing the operational form of the mass-shift calculation."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Provides the constituent quark model whose parameters are used for the bare $B_s$ masses."},{"cited_title":"The newly ob- served open-charm states in quark model","cited_arxiv_id":null,"evidence_quote":"Supplies additional parameter validation for the nonrelativistic quark model in the heavy-light meson sector."},{"cited_title":"Silvestre-Brac and C","cited_arxiv_id":null,"evidence_quote":"Determines the range for the Gaussian regulator $r_q$, whose midpoint is used to regulate the loop integrals."}],"review_version":1}