{"id":"f14ef17b-07b2-4d33-a083-ee5aa5caff1c","arxiv_id":"2411.12866","paper_version":1,"verdict":"ACCEPT","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":0,"one_line_summary":"The authors complete the O(alpha_s^2) perturbative correction to the q^2 spectrum of b -> X_c l nu_l, adding the triple-charm channel, and confirm the existing single-charm result.","lead":"This paper computes the missing next-to-next-to-leading order quantum chromodynamics corrections to the leptonic invariant mass spectrum in inclusive semileptonic bottom quark decays, including the triple-charm final state. The result sharpens the theoretical prediction used to extract the Cabibbo-Kobayashi-Maskawa matrix element |Vcb| from Belle and Belle II measurements.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Triple-charm completeness rests on an unverified cut-MI enumeration; the single-charm cross-check at 5e-4 precision is far looser than the ~1e-8 triple-charm signal and cannot validate that channel.","rationale":"The reader's weakest_assumption correctly identified the completeness of the cut-MI identification in Section 3.1 as the load-bearing point. The present pass adds the quantitative observation that the paper's single-charm cross-check cannot adjudicate this concern: the 5.2e-4 precision of the comparison with ref. [13] is far larger than the relative size of any plausible triple-charm error (~1e-8 to 1e-5). Thus the triple-charm result rests solely on the internal enumeration of 22 cut MIs. This does not imply the result is wrong; standard cut-identification algorithms are reliable, and the two topologies correspond to the two tree-level attachments of the splitting gluon (from the b line or from the c line). But the claim 'complete' would be technically compromised if a topology were missed, and the paper provides no way to verify this from the text alone. An independent direct phase-space computation of the tree-level b -> c c cbar l nu process is feasible and would settle the question. Recommendation: CONDITIONAL acceptance, pending this check or an explicit listing of the cut-MI enumeration.","tokens_in":15943,"tokens_out":15379,"duration_ms":160317,"concrete_test":"Compute the O(alpha_s^2) triple-charm spectrum directly from the tree-level process b -> c c cbar l nu (one gluon splitting, no loops), e.g., by integrating the 4-body phase space of c c cbar plus a vector boson of mass sqrt(q^2) with the effective weak vertex, using sector decomposition or a high-precision Monte Carlo. Evaluate dGamma_3c/dq^2 for a representative m_c (e.g., 0.25) across the allowed q^2 range and compare with the fit (3.1) and Tables 1-2. If the direct result differs by more than ~2% (i.e., beyond the stated 1.5% fit error), the 22 cut-MI enumeration is incomplete or the extraction is wrong; agreement to within the fit error validates the triple-charm channel.","verdict_should_be":"CONDITIONAL","load_bearing_attack":"The central claim is a complete O(alpha_s^2) q^2 spectrum, including the triple-charm channel. Section 3.1 isolates that channel by selecting from the 98 three-loop MIs of the full inclusive calculation those with an allowed cut across three c-quark lines, obtaining 22 cut MIs in two topologies (Fig. 4). The load-bearing assumption is that this enumeration is complete: if a triple-charm cut topology were missed, eq. (3.2) would mis-decompose the spectrum into single- and triple-charm parts, and the 'complete' claim would be inaccurate. The single-charm check against ref. [13] does not settle this: it is a 5.2e-4 relative comparison, while the triple-charm width is at most ~1e-5 of the single-charm width, and typically ~1e-8 or smaller (Fig. 6). A missing triple-charm topology would therefore shift the single-charm comparison by an amount far below its stated precision, so the agreement is consistent with both a complete and an incomplete triple-charm extraction. No other independent check of the triple-charm channel is provided.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper presents a calculation of the O(alpha_s^2) QCD correction to the leptonic invariant mass (q^2) spectrum in the partonic b -> X_c l nu_l decay. The authors use the optical theorem with an auxiliary vector boson of mass sqrt(q^2), generate the three-loop forward-scattering diagrams with QGRAF, reduce the scalar integrals with KIRA, and evaluate the master integrals with AMFlow. They identify 22 cut master integrals that isolate the triple-charm channel b -> c cbar c l nu, define the single-charm channel by subtracting the triple-charm contribution from the fully inclusive result, and provide numerical fits for both channels as functions of qhat^2 and mhat_c. They confirm the single-charm results of Fael and Herren [13] and show that the triple-charm channel has a numerically negligible effect on centralized q^2 moments for physical parameters.","tokens_in":16167,"tokens_out":8768,"duration_ms":93892,"significance":"If correct, this is the first complete NNLO partonic q^2 spectrum including the triple-charm channel, which is a needed ingredient for precision |Vcb| determinations from Belle and Belle-II q^2-moment measurements. The paper's strengths are its use of state-of-the-art multi-loop technology (QGRAF, KIRA, AMFlow), the independent confirmation of the dominant single-charm contribution of ref. [13], and the provision of reusable fit functions in an ancillary Mathematica-readable file. The triple-charm result is new, and the paper appropriately separates this channel from the better-controlled single-charm channel because the heavy-quark expansion is questionable near the triple-charm threshold.","major_comments":[{"comment":"The completeness of the triple-charm channel rests on the statement that, out of all 3-loop master integrals in the fully inclusive calculation, 22 cut master integrals with an allowed cut across three c-quark lines were extracted. The manuscript does not document the selection procedure, nor does it provide an independent check of this enumeration. This is load-bearing because eq. (3.2) subtracts the triple-charm contribution from the fully inclusive spectrum to define the single-charm channel; a missed cut topology would change both channels and invalidate the word 'complete' in the paper's central claim. The 5.2e-4 relative agreement with ref. [13] in section 3.2 cannot validate the triple-charm part, since the triple-charm width is roughly 1e-8 to 1e-5 of the single-charm width (figure 6), so an error in the triple-charm extraction would be far below the quoted comparison accuracy. Please provide a detailed description of how the 22 cut MIs were identified from the 98 MIs, or an independent cross-check such as a reverse-unitarity enumeration, or a comparison of the integrated total width (single-plus-triple charm) with known O(alpha_s^2) results, e.g. from ref. [16].","section":"3.1"}],"minor_comments":[{"comment":"The sentence 'We obtain the numerical results for the q^2 spectrum ... with 60 significant digits' refers to the AMFlow evaluation, but the final fit has a relative error as large as 1.5%. Please clarify that the 60-digit statement applies to the intermediate numerical integration, not to the fitted spectrum.","section":"3.1"},{"comment":"The step function Theta[10 L3c - 1] in the fit ansatz is not explained. Please define the step function and state its role in describing the high-q^2 region of the triple-charm spectrum.","section":"3.1, eq. (3.1)"},{"comment":"The comparison with ref. [13] is quoted as a relative error never exceeding 5.2e-4, while section 3.3 states that the moments agree 'to all digits given' in ref. [13]. Please clarify the relationship between these two statements, since a 5e-4 spectrum-level difference can still yield agreement at the level of the printed moment digits.","section":"3.2"},{"comment":"The fit coefficients in tables 1-6 are presented without any indication of their numerical precision or the resulting fit uncertainty. Since these fits are the main output of the paper, please state explicitly how many significant digits of the coefficients are meaningful and how the residuals quoted in sections 3.1 and 3.2 were computed.","section":"Appendix A"}],"recommendation":"major_revision","confidential_remarks":"The manuscript reports a solid, technically demanding calculation and is likely correct. My recommendation of major revision is driven solely by the need to substantiate the completeness of the triple-charm cut-master-integral enumeration, which is the load-bearing step for the paper's central claim. If the authors can document the selection from the 98 master integrals or provide an independent check of the triple-charm channel, I would support publication."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Mateusz, here's my read.\n\nThe genuinely new content is the first NNLO treatment of the b -> c cbar c l nu channel in the q^2 spectrum, computed via three-loop unitarity cuts, plus an independent confirmation of the Fael-Herren single-charm result. The single-charm moments agree to all digits quoted there; the triple-charm piece is small and carefully quantified.\n\nThe paper is technically solid. The pipeline (QGRAF, KIRA IBPs, AMFlow) is standard, and the authors work at 60-digit precision before fitting clean ansatze with documented errors. The fits are in an appendix and an ancillary file, so the results are directly usable. The replacement of the lepton pair by a massive auxiliary vector boson is a known trick, and there is no circularity: the spectrum is computed directly, and the fits are numerical representations rather than inputs.\n\nThe soft spot is the completeness of the triple-charm cut-MI enumeration. The claim of a complete O(alpha_s^2) correction rests on the statement that exactly 22 of the 98 three-loop master integrals belong to two topologies that admit an allowed cut across three c-quark lines. There is no independent cross-check of that enumeration for the triple-charm channel. The check against ref. [13] cannot provide one: it is a 5.2e-4 relative comparison of single-charm spectra, while the triple-charm width is five to nine orders of magnitude smaller. A missing triple-charm topology would shift that comparison by far less than its stated precision, so the agreement is compatible with both a complete and an incomplete triple-charm list. The authors are transparent about fit errors, but this structural assumption is unverified. I do not regard it as fatal: even if a topology were missed, the fully inclusive spectrum would be unaffected, and the triple-charm effect on the moments is at the 1e-9 to 1e-6 level for physical inputs. The practical impact on |Vcb| fits is negligible. Still, the word 'complete' is doing real work and deserves either a stronger argument or a small caveat.\n\nWho this is for: anyone doing inclusive semileptonic B decays, HQET, or |Vcb| extraction from Belle/Belle II q^2 moments. It belongs in JHEP.\n\nRecommendation: send it to peer review. I would accept, asking for a more systematic enumeration check of the triple-charm cuts, or at minimum an explicit statement of this caveat and why it is phenomenologically irrelevant.","headline":"Completes the NNLO q2 spectrum with the first triple-charm treatment and confirms the single-charm result; the unverified cut-MI completeness is a real but phenomenologically harmless caveat.","tokens_in":16736,"tokens_out":3440,"would_cite":true,"duration_ms":30975,"reading_group":"yes","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":["12.38.Bx","13.20.He","14.65.Fy"],"model":"deepseek-v4-flash","headline":"This paper computes the complete $O(\\\\alpha_s^2)$ correction to the partonic $q^2$ spectrum of $b\\\\to X_c l\\\\bar{\\\\nu}_l$, including the triple-charm channel, and confirms the earlier single-charm result.","keywords":["inclusive semileptonic B decays","q^2 moments","leptonic invariant mass spectrum","NNLO QCD corrections","triple-charm channel","|V_cb| extraction","heavy quark expansion","master integrals"],"falsifier":"Compute the $b\\\\to c\\\\bar{c}c\\\\,l\\\\bar{\\\\nu}_l$ differential rate at a representative kinematic point, for example $\\\\hat m_c = 0.27$ and $\\\\hat q^2 = 0.1$, by an independent method such as direct integration of the squared amplitude over the three-charm phase space, and compare it with the fit in eq. (3.1); a disagreement beyond the stated fit error of $1.5\\\\%$ would show that the cut-master-integral set or the subtraction is incomplete.","tokens_in":15722,"feed_emoji":"⚛️","tokens_out":11452,"duration_ms":104157,"temperature":0.7,"pith_summary":"This paper completes the next-to-next-to-leading-order ($O(\\\\alpha_s^2)$) QCD correction to the leptonic-invariant-mass squared ($q^2$) spectrum of the partonic $b\\\\to X_c l\\\\bar{\\\\nu}_l$ decay, the free-quark-level process underlying inclusive semileptonic $B$-meson decays. It adds the previously omitted triple-charm channel, whose final hadronic state is $c\\\\bar{c}c$, and separates it from the dominant single-charm channel. It also independently confirms the single-charm NNLO result already in the literature. The final spectrum is given as numerical fits in $q^2$, the charm-quark mass, and the renormalization scale, ready for computations of $q^2$ moments that are used to extract $|V_{cb}|$.","feed_headline":"NNLO q² spectrum is now complete for b→X_c lν decays","feed_subtitle":"New calculation adds the c c̄ c final state and confirms the single-charm result behind |V_cb| fits.","key_machinery":"The argument runs through an effective-width formula (eq. 2.13) that replaces the lepton pair by a light auxiliary vector boson, so the $q^2$ spectrum becomes the total width of $b\\\\to X_c W$ and can be obtained from imaginary parts of three-loop forward-scattering diagrams via the optical theorem. Integration-by-parts identities reduce those diagrams to 98 master integrals in five topologies; the triple-charm channel is isolated by keeping only unitarity cuts that separate three charm-quark lines, which leaves 22 master integrals in two topologies. All master integrals are evaluated numerically with 60-digit precision at hundreds of kinematic points, then fitted to the ansatz functions of eqs. (3.1) and (3.3), with renormalization performed in the on-shell scheme for quark masses and in the $\\\\overline{\\\\rm MS}$ scheme for $\\\\alpha_s$.","core_discovery":"The central claim is that the complete $O(\\\\alpha_s^2)$ correction to the partonic $q^2$ spectrum is now known, with the triple-charm $b\\\\to c\\\\bar{c}c\\\\,l\\\\bar{\\\\nu}_l$ contribution evaluated explicitly and subtracted from the fully inclusive correction to define the single-charm piece. The single-charm piece reproduces the result of ref. [13], including exact agreement with all digits of the centralized moments published there. The triple-charm piece, shown for the first time as a function of $q^2$, has a negligible effect on the first four normalized centralized moments for physical quark masses and realistic $q^2$ cuts, decreasing rapidly as either $\\\\hat m_c$ or the cut increases. This means that previous single-charm-only NNLO analyses were numerically adequate for current fits, while the present fits make the complete correction available for any future cut choice.","pith_inferences":["An independent evaluation of the triple-charm spectrum at one kinematic point, for example by reverse unitarity rather than cut master integrals, would test the completeness of the 22-integral set without requiring a full second computation.","Because the triple-charm effect on normalized centralized moments is so small, the limiting uncertainty for future $|V_{cb}|$ fits from $q^2$ moments will likely come from power-suppressed non-perturbative matrix elements and from renormalization-scale variation, not from this channel.","The same cut-based channel decomposition could be applied to other inclusive observables to separate final-state multiplicities whose heavy-quark-expansion behaviour differs, not just in semileptonic $B$ decays."],"forward_implications":["Any $q^2$ moment with an arbitrary lower cut can now be computed at NNLO by integrating the provided fits in $\\\\hat q^2$.","The centralized $q^2$ moments used in $|V_{cb}|$ fits can be quoted with the triple-charm channel included, eliminating a previously unquantified missing-channel assumption.","The single-charm NNLO spectrum of ref. [13] is independently confirmed, including exact agreement with its published moment values.","For physical quark masses and non-zero $q^2$ cuts, the relative triple-charm correction to the first four centralized moments is below $10^{-4}$ over most of the parameter range, so earlier single-charm-only analyses remain numerically sound.","The normalized spectrum shows better perturbative convergence than the unnormalized one, because the overall factor of $(m_b^{\\\\rm OS})^5$ and part of the renormalon ambiguity cancel in the width ratio."],"supporting_citations":[{"why":"Provides the single-charm NNLO $q^2$ spectrum that this paper confirms and extends; it is the central benchmark for comparison and moments.","marker":"[13]"},{"why":"Introduces the effective-width formula used throughout to turn the $q^2$ spectrum into a $b\\\\to X_c W$ width.","marker":"[19]"},{"why":"Generates the Feynman diagrams for the LO, NLO, and NNLO amplitudes.","marker":"[21]"},{"why":"Supplies the integration-by-parts reduction method used to organise the diagrams into master integrals.","marker":"[22, 23]"},{"why":"Supplies the numerical auxiliary-mass-flow method used to evaluate all NNLO master integrals at many kinematic points.","marker":"[27]"},{"why":"Provides the renormalization constants for $\\\\alpha_s$ and the quark masses used in the subtraction.","marker":"[28–30]"},{"why":"Defines the centralized moments whose triple-charm impact is analysed.","marker":"[31]"}],"fun_headline_variants":["Complete NNLO q² spectrum for b→X_c lν now available","Triple-charm channel added: NNLO b→X_c lν spectrum complete","Full NNLO corrections to b→X_c lν q² spectrum done","New triple-charm piece completes NNLO b→X_c lν spectrum","NNLO q² spectrum is complete with triple-charm for b→X_c lν"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The load-bearing assumption is that every possible $b\\\\to c\\\\bar{c}c\\\\,l\\\\bar{\\\\nu}_l$ contribution at this order is captured by the 22 master integrals from the two diagram families identified in figure 4; if a contributing family were missed, the separation into single- and triple-charm parts would be wrong even if the total inclusive spectrum were unaffected.","fun_headline_variants_meta":{"raw":{"variants":["Complete NNLO q² spectrum for b→X_c lν now available","Triple-charm channel added: NNLO b→X_c lν spectrum complete","Full NNLO corrections to b→X_c lν q² spectrum done","New triple-charm piece completes NNLO b→X_c lν spectrum","NNLO q² spectrum is complete with triple-charm for b→X_c lν"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000936,"raw_usage":{"total_tokens":3988,"prompt_tokens":917,"completion_tokens":3071,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":533,"completion_tokens_details":{"reasoning_tokens":2965}},"tokens_in":533,"tokens_out":3071,"duration_ms":23216,"temperature":1.0,"reasoning_tokens":2965,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-12T17:06:20.724835+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Compute the $b\\\\to c\\\\bar{c}c\\\\,l\\\\bar{\\\\nu}_l$ differential rate at a representative kinematic point, for example $\\\\hat m_c = 0.27$ and $\\\\hat q^2 = 0.1$, by an independent method such as direct integration of the squared amplitude over the three-charm phase space, and compare it with the fit in eq. (3.1); a disagreement beyond the stated fit error of $1.5\\\\%$ would show that the cut-master-integral set or the subtraction is incomplete.","supporting_citations":[],"review_version":1}