{"id":"11e5edc1-33f3-43af-b304-3c7406a9cf93","arxiv_id":"2412.02997","paper_version":3,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":13,"one_line_summary":"A QCD-sum-rule plus rescattering calculation finds that X0(2900) as a D*bar K* molecule reproduces the measured B+ to D+X0 branching fraction and X0 to DK width, though a compact tetraquark remains allowed within errors.","lead":"This paper calculates how often B+ decays produce the exotic particle X0(2900) and how fast X0 decays, using QCD sum rules and a rescattering model. The D*bar K* molecule picture matches experiment within large uncertainties, while the compact tetraquark picture is not ruled out.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Dominant FSI loop integrals sample vertex couplings at Q²≈0.04–1 GeV², far below the sum-rule window; fit shape parameters are quoted with zero uncertainty, so the claimed LHCb consistency is not robust.","rationale":"The reader correctly identified the off-shell extrapolation of QCD sum rule couplings as the weakest assumption. My concern sharpens it: the FSI loop integrals do not simply evaluate the couplings at the physical K pole; they sample Q² ≈ 0.04–8 GeV², with the propagator emphasizing the lowest Q² values. The sum rules are only computed at Q² ≥ 3.5 GeV², and the fitted functional forms are used far outside that region with zero uncertainty quoted for the shape parameters. The paper itself notes that the form factor does not always behave as a monopole, which underscores the model dependence. A change in the vertex shapes at low Q² would alter the dominant amplitudes and could move the branching fraction and width away from the experimental values. Since the reader's verdict is already CONDITIONAL and this concern is a concrete instance of the identified weakness, I recommend keeping the verdict unchanged. The proposed test—re-fitting with alternative functional forms and floating the shape parameters—would determine whether the consistency with LHCb survives a reasonable variation of the extrapolation.","tokens_in":26050,"tokens_out":17918,"duration_ms":156414,"concrete_test":"Recompute the dominant amplitudes by replacing the monopole/exponential fits for g_DKX0, g_D*K*X0, and g_D*sDK with alternative forms (e.g., dipole g = a/(b+Q²)² or a Padé form) that pass through the same Q² ∈ [3.5, 6] GeV² sum-rule points, and allow the shape parameter b to float so its uncertainty is propagated through Eqs. (11)–(20) and Eq. (68). If the resulting Br(B⁺→D⁺X0) and Γ(X0→DK) move outside the LHCb windows (2.4 ± 1.0) × 10⁻⁵ and 57 ± 13 MeV, the consistency claim depends on an unvalidated extrapolation.","verdict_should_be":"UNCHANGED","load_bearing_attack":"In the FSI loop integrals Eqs. (11) and (13), the exchanged kaon momentum squared t is negative throughout the physical integration region (for diagram (a), t ∈ [-7.98, -0.04] GeV², so Q² = -t ∈ [0.04, 7.98] GeV²). The propagator denominator (Q² + m_K²) is smallest near Q² ≈ 0, making the loop dominated by the low-Q² part of the vertex functions g(Q²). However, the QCD sum rules are evaluated only for Q² ≳ 3.5–4 GeV² (Figs. 4, 5, 8), and the paper extrapolates to small positive Q² using ad hoc monopole/exponential fits (Eqs. 44, 48; Table III) whose shape parameters b are quoted with zero uncertainty (e.g., b = 57.43 ± 0.00 and b = 0.22 ± 0.00). The sum-rule expression itself, Eq. (42), contains a factor (Q² + m_K²)/Q², so the OPE cannot be trusted near Q² = 0. The dominant amplitudes A_a, A_b, and A_d receive large contributions from Q² ≲ 1 GeV², so an unconstrained change in the vertex shape there can shift the branching fraction by a large factor. The decay width Γ(X0 → DK) likewise relies on extrapolating g_DKX0 to Q² = -m_K², where the same sum rule is invalid. Because the quoted uncertainties omit this shape error, the central claim of consistency with LHCb is not yet established.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper studies the production and decay of the open-flavor tetraquark candidate X0(2900) using a final-state-interaction rescattering model. Strong vertices such as g_{\\bar D K X_0}, g_{\\bar D^* K^* X_0}, and the D_s^(**) D K^(**) vertices are computed with three-point QCD sum rules, and the resulting couplings are used in hadronic loop amplitudes for B^+ -> D^+ X_0. For the \\bar D^* K^* molecular interpretation, the authors obtain Br(B^+ -> D^+ X_0) = (2.47 +2.07 -1.55) x 10^-5 and \\Gamma(X_0 -> \\bar D K) = 64.07 ± 18.23 MeV, both consistent with the LHCb values of about (2.4 ± 1.0) x 10^-5 and 57 ± 12 ± 4 MeV. The paper concludes that the molecular interpretation is favored but that the compact tetraquark interpretations cannot be excluded within uncertainties.","tokens_in":26610,"tokens_out":4746,"duration_ms":50876,"significance":"If the extrapolation procedure is robust, the paper provides a valuable comparison of three X0(2900) interpretations within one framework, replacing the ad hoc monopole cutoff of earlier FSI studies with sum-rule-derived vertex functions. It also presents a large set of three-point sum rules and on-shell couplings that can be checked against other approaches. The main significance is the claim of quantitative consistency with LHCb for the molecular scenario, together with an explicit statement that the compact tetraquark is not excluded. This claim, however, currently rests on long and untested off-shell extrapolations of the sum-rule couplings, so the significance is conditional on those extrapolations being validated.","major_comments":[{"comment":"The couplings g_{\\bar D K X_0} and g_{\\bar D^* K^* X_0} are extracted only for Q^2 ≳ 3.5-4 GeV^2 and then extrapolated to the physical pole and to the low-Q^2 region using monopole or exponential fits. The FSI loop integrals in Eqs. (11) and (13) sample exchanged momenta with Q^2 = -t ∈ [0.04, 7.98] GeV^2 for diagram (a), and the integrand is largest near Q^2 ≈ 0 where the propagator is smallest. In that region Eq. (42) contains the factor (Q^2 + m_K^2)/Q^2, so the OPE cannot be trusted. Since Table V shows that A_a, A_b, and A_d dominate the branching fraction, an unconstrained change in the low-Q^2 behavior of the vertex functions can shift Br by a large factor. The paper provides no independent check of the monopole/exponential extrapolations, no alternative functional form, and no estimate of this shape uncertainty. This is a load-bearing uncertainty for the central consistency claim.","section":"Sec. III A and Sec. IV, Eqs. (42)-(49), (11), (13), and Table V"},{"comment":"The continuum threshold s0 = 10 GeV^2 for the X0 vertices is taken from Ref. [38], which shares an author with this paper. The sum-rule results for g_{\\bar D K X_0} and g_{\\bar D^* K^* X_0} are quoted only for this single s0 value, and no sensitivity study of s0 in the three-point sum rules is reported. Because the spectral integrals in Eqs. (42) and (47) depend explicitly on s0, the quoted uncertainties likely underestimate the model dependence. I request a variation of s0 (for example ±0.5-1 GeV^2) or a stability test showing that the extracted couplings and the final Br and width are insensitive to this choice.","section":"Sec. III A, Eqs. (42) and (47), and s0 input"},{"comment":"Several fit parameters are quoted with zero uncertainty, e.g., b = 57.43 ± 0.00 GeV^2 in Eq. (45), b = 0.22 ± 0.00 GeV^-2 in Eq. (49), and b = 0.22 ± 0.00 in Table III. This is not credible for fits to data points that themselves carry errors, and it hides the dominant systematic of the off-shell extrapolation. The authors should either propagate the fit-parameter covariance or assign a realistic uncertainty to the shape parameters; the present error budget for Br and \\Gamma omits this uncertainty.","section":"Tables II and III, Eqs. (45), (49), and (67)"}],"minor_comments":[{"comment":"Eq. (67) reads e^{-(-0.22 ± 0.01) Q^2}, which has a double negative in the exponent and would grow with Q^2. The decreasing curve in Fig. 8 indicates that e^{-0.22 Q^2} is intended; this typo should be corrected because a reader using Eq. (67) at face value would obtain the wrong sign of the slope.","section":"Sec. III B, Eq. (67)"},{"comment":"The text says the decay width is derived from Eq. (52), but Eq. (52) is the transition matrix element for \\bar D_1 K_1 X_0. The correct reference for the \\bar D K X_0 matrix element is Eq. (39).","section":"Sec. IV, Eq. (68)"},{"comment":"The manuscript contains several typographical artifacts, such as 'resoncance', 'e fforts', 'of fective', and 'compansate', which should be cleaned up in a final version.","section":"Throughout"},{"comment":"The comparison with experimental Br(B^+ -> D^+ X_0) uses the isospin-inferred value (2.4 ± 1.0) x 10^-5. It would be helpful to state explicitly that this value inherits the LHCb uncertainty and that the quoted consistency does not include the experimental errors in the theory error budget.","section":"Sec. I, Eq. (3)"}],"recommendation":"major_revision","confidential_remarks":"The paper's use of s0 = 10 GeV^2 from Ref. [38] is not a direct fit to the X0 observables, but Ref. [38] is co-authored by one of the current authors, so the input is not fully independent. The central claim of consistency with LHCb is interesting but is dominated by the unquantified low-Q^2 extrapolation of the vertex functions; the requested robustness tests are essential before the claim can be accepted. The manuscript is within the scope of the journal."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Short version: this is a real calculation, not a toy. It computes a batch of strong vertices with three-point QCD sum rules and feeds them into an eight-diagram FSI model for B+ -> D+ X0(2900). For the D*bar K* molecular interpretation it gets Br = (2.47 +2.07 -1.55) x 10^-5 and Gamma(X0 -> DK) = 64.07 +- 18.23 MeV, which overlap the LHCb values. That is a legitimate data point for a difficult four-quark candidate. The paper is also honest: it says explicitly that the compact tetraquark interpretation cannot be excluded within the uncertainties.\n\nWhat is genuinely new: the vertex list is broader than in earlier rescattering studies, and the finding that g_D*barK*X0 behaves exponentially in Q^2 rather than as a monopole is new and worth knowing. The OPE side is detailed enough to reproduce, and the comparisons with earlier determinations in Table III show the calculation sits in the expected range. The citation pattern is appropriate; the literature on X0(2900) is well covered.\n\nThe soft spot is the extrapolation. The two largest amplitudes A_a and A_b are driven by vertex couplings at Q^2 near and below 1 GeV^2, while the sum rules are evaluated only around Q^2 = 3.5-6 GeV^2. The extrapolation through monopole or exponential fits is long, and the fit shape parameters are quoted with zero uncertainty (e.g., b = 57.43 +- 0.00 and b = 0.22 +- 0.00). Since the sum-rule expression itself contains (Q^2 + m_K^2)/Q^2, the OPE cannot be trusted near Q^2 = 0, where the loop integrand peaks. No independent check of the extrapolated vertex is provided. That means the agreement with LHCb, while real, is not robust: a different but still plausible low-Q^2 shape could shift the branching fraction by a large factor. The decay width also requires g_DKX0 at Q^2 = -m_K^2, another long extrapolation. The continuum threshold s0 = 10 GeV^2 is imported from a paper co-authored by W. Chen with no sensitivity scan; this is a mild point but should be addressed.\n\nWho it is for: hadron phenomenologists working on open-flavor exotics, especially the X0(2900) question. It will not settle the molecule-versus-tetraquark debate, but it is one of the more complete production-side calculations and deserves referee time. A serious referee should push for a sensitivity analysis of the low-Q^2 extrapolation and a scan over s0, but the core calculation is coherent and the claims are stated honestly.","headline":"A useful QCD sum rule + FSI calculation of X0(2900) production that moderately favors the D*K* molecule, but the central consistency claim leans heavily on long, unvalidated vertex extrapolations.","tokens_in":27152,"tokens_out":2585,"would_cite":true,"duration_ms":26544,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":["12.39.Mk","12.38.Lg","14.40.Ev","14.40.Rt"],"model":"deepseek-v4-flash","headline":"The X0(2900) state's production and decay both match a D* K* molecule interpretation, computed from QCD sum rules and final-state rescattering.","keywords":["X0(2900)","open-flavor tetraquark","hadronic molecule","D* K* molecular interpretation","final-state interaction","QCD sum rules","B meson production","exotic hadron"],"falsifier":"A lattice-QCD computation of the off-shell $\\bar D K X_0$ and $\\bar D^* K^* X_0$ couplings at the kaon and $K^*$ poles would settle the extrapolation; alternatively, measuring the ratio of $B^0\\to D^0 X_0\\to D^0\\bar D^0 K^0$ to $B^+\\to D^+ X_0\\to D^+D^-K^+$ would distinguish the molecular-rescattering prediction of equal rates from the triangle-singularity prediction of strong suppression.","tokens_in":25840,"feed_emoji":"⚛️","tokens_out":10884,"duration_ms":86079,"temperature":0.7,"pith_summary":"The paper asks what the X0(2900), a resonance seen in $B^+\\to D^+D^-K^+$ decays, is made of. It computes the branching fraction of $B^+\\to D^+ X_0$ and the width of $X_0\\to DK$ under three interpretations: a $\\bar D^*K^*$ hadronic molecule, a scalar diquark-scalar antidiquark compact tetraquark, and an axial-vector diquark-axial-vector antidiquark tetraquark. The calculation uses QCD sum rules for the strong vertices and a final-state rescattering mechanism for the production. For the $\\bar D^* K^*$ molecule picture, the branching fraction $(2.47^{+2.07}_{-1.55})\\times 10^{-5}$ and width $(64.07\\pm 18.23)\\,\\mathrm{MeV}$ agree with the measured values, while the compact tetraquark pictures are not ruled out. The point of the work is that the molecular interpretation is a viable, quantitative description of the observed state.","feed_headline":"A D* K* molecule fits X0(2900) production and decay","feed_subtitle":"QCD sum rules plus final-state rescattering reproduce the measured rate and width of the exotic state.","key_machinery":"The load-bearing object is the set of strong coupling constants $g_{ABC}(Q^2)$ at hadronic vertices, computed from three-point QCD sum rules and then extrapolated from the Euclidean momentum region where the sum rule is valid ($Q^2\\sim 3.5$ to $6$ GeV$^2$) to the physical meson pole ($Q^2=-m_K^2$ or $-m_{K^*}^2$) with fitted monopole or exponential forms. These $Q^2$-dependent couplings replace the common monopole form factor, because the sum-rule results do not always follow a monopole curve; the $\\bar D^* K^* X_0$ coupling in particular is better described by an exponential. The production amplitude is assembled from eight final-state-interaction rescattering diagrams, with the weak $b\\to c(c\\bar s)$ transition factorized into known $B\\to D^{(*)}$ form factors and meson decay constants. The two largest contributions come from the $D_s^*\\bar D\\to DX_0$ ($K$ exchange) and $D_{s1}\\bar D^*\\to DX_0$ ($K^*$ exchange) rescattering processes.","core_discovery":"On the paper's own terms, the central result is that the $\\bar D^* K^*$ molecular interpretation of $X_0(2900)$ survives a quantitative comparison with experiment. The authors build the production amplitude $B^+\\to D^+X_0$ from eight rescattering diagrams in which the weakly produced $c\\bar s$ pair hadronizes into $D_s$ or $D_{s1}$ and rescatters off $D$ or $D^*$ mesons by exchanging $K$, $K^*$, $K_1$, or $K_0^*$ mesons. The strong couplings at these vertices are obtained from three-point QCD sum rules, including $g_{\\bar D K X_0}$ and $g_{\\bar D^* K^* X_0}$. Summing the amplitudes gives $\\mathrm{Br}(B^+\\to D^+ X_0)=(2.47^{+2.07}_{-1.55})\\times 10^{-5}$ and $\\Gamma(X_0\\to DK)=(64.07\\pm 18.23)\\,\\mathrm{MeV}$ for the molecule, consistent with the measured values of roughly $(2.4\\pm 1.0)\\times 10^{-5}$ and $(57\\pm 12\\pm 4)\\,\\mathrm{MeV}$. The authors also state that the compact tetraquark interpretation cannot be excluded within the uncertainties.","pith_inferences":["The extrapolation of the sum-rule vertices from positive $Q^2$ to the meson pole is the step the numerical agreement hinges on; a lattice-QCD determination of $g_{\\bar D K X_0}$ at the kaon pole would provide a direct check.","The equal-branching prediction for the $B^+$ and $B^0$ cascade modes is a sharper discriminant between a genuine resonance and a triangle-singularity interpretation than the width alone, because the triangle model is argued in the literature to suppress some of these channels.","The same QCD-sum-rule-plus-rescattering machinery can be applied to other open-flavor candidates near a $\\bar D^{(*)}K^{(*)}$ threshold, making the approach a general test of molecular interpretations.","A similar analysis of the companion $X_1(2900)$ state, whose larger width remains unexplained here, would show whether the molecular picture extends to the full observed doublet."],"forward_implications":["If the molecular interpretation is correct, the observed $X_0(2900)$ is a hadronic molecule rather than a compact four-quark state, and its dominant strong decay into $DK$ follows naturally from the $\\bar D^* K^*$ constituents.","The production rate in $B$ decays is set largely by long-distance rescattering, so measuring $B\\to DX_0$ channels probes how the weak-produced $c\\bar s$ pair hadronizes before final-state interactions.","The four cascade channels $B^+\\to D^+X_0\\to D^+(D^-K^+)$, $B^+\\to D^+X_0\\to D^+(\\bar D^0 K^0)$, $B^0\\to D^0 X_0\\to D^0(\\bar D^0K^0)$, and $B^0\\to D^0X_0\\to D^0(D^-K^+)$ are predicted to have equal branching fractions, a consequence of factorization combined with $X_0\\to DK$ dominance.","Tighter future measurements of the width and production rate would separate the three interpretations, since the molecular values sit closest to the current central data."],"supporting_citations":[{"why":"Supplies the experimental mass, width, and branching-fraction data for $X_0(2900)$ that the calculation is compared with.","marker":"[22, 23]"},{"why":"Provides the $D_{s1}\\bar D^*\\to DX_0$ rescattering process with $K^*$ exchange, one of the two dominant amplitudes.","marker":"[24]"},{"why":"Provides the $D_s^*\\bar D\\to DX_0$ rescattering process with $K$ exchange, the other dominant amplitude, and the monopole form-factor treatment replaced here by sum-rule couplings.","marker":"[48]"},{"why":"Two-point QCD sum rule for the $\\bar D^*K^*$ molecule, supplying the continuum threshold $s_0=10$ GeV$^2$ in the vertex sum rules.","marker":"[38]"},{"why":"Earlier light-cone sum-rule calculation of the molecular decay width, giving the comparison value $49.6\\pm 9.3$ MeV.","marker":"[39]"},{"why":"Effective-Lagrangian triangle-diagram calculation of $X_0\\to DK$ that gives a reference width of $59.37^{+24.94}_{-17.96}$ MeV for the molecular picture.","marker":"[44]"},{"why":"Parameterizes the $B\\to D^{(*)}$ weak transition form factors used in the factorized production amplitudes.","marker":"[57]"}],"fun_headline_variants":["D* K* molecule explains X0(2900) production and decay","Molecular X0(2900) consistent with measured rates","X0(2900) data consistent with D* K* molecule","QCD sum rules support D* K* molecule for X0(2900)","D* K* molecule fits X0(2900) production and decay"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"Everything rests on extrapolating the QCD sum-rule vertex couplings, which are computed only at positive momentum transfers around $3.5$ to $6$ GeV$^2$, down to the negative momentum transfer where the exchanged kaon or $K^*$ is on shell; if that extrapolation is not the true behaviour of the vertices, the two largest rescattering amplitudes change and the agreement with experiment is no longer guaranteed.","fun_headline_variants_meta":{"raw":{"variants":["D* K* molecule explains X0(2900) production and decay","Molecular X0(2900) consistent with measured rates","X0(2900) data consistent with D* K* molecule","QCD sum rules support D* K* molecule for X0(2900)","D* K* molecule fits X0(2900) production and decay"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000341,"raw_usage":{"total_tokens":1933,"prompt_tokens":1057,"completion_tokens":876,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":673,"completion_tokens_details":{"reasoning_tokens":779}},"tokens_in":673,"tokens_out":876,"duration_ms":7481,"temperature":1.0,"reasoning_tokens":779,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-11T22:53:26.860263+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"A lattice-QCD computation of the off-shell $\\bar D K X_0$ and $\\bar D^* K^* X_0$ couplings at the kaon and $K^*$ poles would settle the extrapolation; alternatively, measuring the ratio of $B^0\\to D^0 X_0\\to D^0\\bar D^0 K^0$ to $B^+\\to D^+ X_0\\to D^+D^-K^+$ would distinguish the molecular-rescattering prediction of equal rates from the triangle-singularity prediction of strong suppression.","supporting_citations":[{"cited_title":"Aaij et al.[LHCb], Phys","cited_arxiv_id":null,"evidence_quote":"Provides the $D_{s1}\\bar D^*\\to DX_0$ rescattering process with $K^*$ exchange, one of the two dominant amplitudes."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Provides the $D_s^*\\bar D\\to DX_0$ rescattering process with $K$ exchange, the other dominant amplitude, and the monopole form-factor treatment replaced here by sum-rule couplings."},{"cited_title":"Liu, J.-J","cited_arxiv_id":null,"evidence_quote":"Two-point QCD sum rule for the $\\bar D^*K^*$ molecule, supplying the continuum threshold $s_0=10$ GeV$^2$ in the vertex sum rules."},{"cited_title":"Wang, and S.-L","cited_arxiv_id":null,"evidence_quote":"Earlier light-cone sum-rule calculation of the molecular decay width, giving the comparison value $49.6\\pm 9.3$ MeV."},{"cited_title":"Buchalla, A.J","cited_arxiv_id":null,"evidence_quote":"Parameterizes the $B\\to D^{(*)}$ weak transition form factors used in the factorized production amplitudes."}],"review_version":1}