{"id":"161d7e87-07a9-4d13-ba84-6d37183d9889","arxiv_id":"2505.22740","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"low","formal_verification":"none","parameter_count":3,"one_line_summary":"The full charm-mass dependence of the three-loop current-current contribution to the B-meson decay matrix is computed, yielding Delta_Gamma_s = (0.077 +/- 0.016) ps^-1 with the leading-term perturbative error reduced to the experimental level.","lead":"This paper computes the three-loop quantum chromodynamics (QCD) corrections to the width difference of neutral B mesons, updating Standard Model predictions such as Delta_Gamma_s = (0.077 +/- 0.016) ps^-1. A generalist should care because these precision predictions are the yardstick for new-physics searches in B-meson mixing, and the remaining uncertainty is now largely non-perturbative rather than perturbative.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Subleading 1/N_c evanescent-scheme dependence of the NNLO matching coefficients is unverified; it could shift ΔΓ_s beyond the quoted scale uncertainty.","rationale":"The paper presents a highly nontrivial three-loop calculation with an unusually strong verification structure: two independent computational setups, analytic master integrals checked numerically with pySecDec, reproduction of NLO and fermionic NNLO benchmarks, and regulator/gauge independence checks for α1 and α2. The phenomenological prediction ΔΓ_s = (0.077 ± 0.016) ps^-1 is reasonably robust against the mass-scheme choice (MS vs. PS). However, the construction of the evanescent operator basis in the |ΔB|=2 theory is the conceptual foundation of the matching: the chosen basis is designed to be compatible with four-dimensional lattice matrix elements by enforcing Fierz symmetry. The paper explicitly verifies scheme independence only for the leading O(N_c^2) term, leaving the subleading 1/N_c terms unexamined. Since the final numerical result for N_c=3 includes contributions at all orders in 1/N_c, an unquantified scheme dependence of the NNLO matching coefficients could in principle alter the central value or the error budget. This is a genuine soft spot, not an internal inconsistency; it is exactly the kind of subtlety that the conditions in Sect. 2.3 are meant to control but do not fully pin down. A direct recomputation with an alternative allowed evanescent scheme would settle the matter. Until then, the reader's CONDITIONAL verdict is appropriate: the claim is likely correct but rests on an unverified subleading scheme independence.","tokens_in":45093,"tokens_out":9091,"duration_ms":105561,"concrete_test":"Recompute the NNLO matching coefficients H^(2) and \tilde H_S^(2) (and hence ΔΓ_s in Eq. (84)) using a different allowed solution for the second-generation evanescent constants, for instance by swapping the regular/Fierz assignments in Eqs. (32)–(34) or by setting k and \tilde k to non-zero values while still satisfying conditions 1–3 of Sect. 2.3. Compare the resulting physical matching coefficients with those published. If they are identical (or the shift in ΔΓ_s is below 0.001 ps^-1), the scheme ambiguity is numerically irrelevant; if they differ by more than the largest scale uncertainty in Eq. (83), the NNLO claim requires an additional scheme-ambiguity uncertainty and the verdict should be CONDITIONAL or REJECT.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The Fierz-symmetric renormalisation scheme of the |ΔB|=2 theory (Sect. 2.3) is fixed by conditions 1–3 only up to a solution space for the second-generation evanescent coefficients. The paper selects one solution (Eqs. (31)–(34), App. D) and notes in App. D that the leading O(N_c^2) term of the renormalised physical matrix elements is independent of the remaining undetermined constant. No equivalent check is shown for the subleading 1/N_c terms. For N_c=3 these subleading terms are not numerically negligible, and the physical NNLO matching coefficients H^(2) and \\tilde H_S^(2) entering Eq. (84) could in principle depend on the arbitrary choice of the remaining constants. Because the lattice bag parameters used for the non-perturbative input are four-dimensional (Fierz-symmetric) quantities, they do not automatically cancel a perturbative scheme dependence of the Wilson coefficients at order α_s^2. If a different allowed choice of the evanescent constants shifts ΔΓ_s by more than the quoted scale uncertainty (about 0.004–0.008 ps^-1 from Eq. (83)), the central claim that the NNLO perturbative uncertainty is reduced to the experimental level would need revision. The paper's strong internal cross-checks (two independent setups, IR pole cancellation, reproduction of literature benchmarks) do not address this residual scheme ambiguity, making it the most load-bearing untested assumption.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"This paper presents the first complete three-loop (NNLO) QCD calculation of the current-current operator contribution to the off-diagonal decay matrix element Γ_12 in B_q–\\bar B_q mixing, retaining the full dependence on z = m_c^2/m_b^2 through a semi-analytic expansion of the three-loop master integrals. The authors construct a Fierz-symmetric renormalisation scheme for the |ΔB|=2 effective theory, compute the matching coefficients H^(2) and \\tilde H_S^(2), and perform a phenomenological analysis of ΔΓ_q, ΔΓ_q/ΔM_q, and a_fs^q for q = d, s. Their final result is ΔΓ_s = (0.077 ± 0.016) ps^-1, with the scale uncertainty of the leading 1/m_b term reduced to a level comparable to the current experimental error; the remaining uncertainty is dominated by the 1/m_b-suppressed hadronic matrix elements.","tokens_in":45306,"tokens_out":12282,"duration_ms":138942,"significance":"If the result stands, this is a substantial step forward in B-meson mixing phenomenology. The NNLO correction stabilises the renormalisation-scale dependence of the leading-power term, and the deep semi-analytic expansion in z makes the result usable over the full physical range. The paper is distinguished by strong internal cross-checks: two independent automated setups, reproduction of the NLO benchmark of Ref. [11] and the fermionic NNLO benchmark of Ref. [17], cancellation of IR poles in the matching (Eq. (67)), gauge-parameter independence of the α_i constants, and analytic master integrals checked with pySecDec. The authors also provide computer-readable ancillary files for the renormalisation constants and matching coefficients, which is a valuable resource. The principal reservation is the unverified independence of the physical NNLO matching coefficients from the residual freedom in the second-generation evanescent operator definitions at subleading order in 1/N_c.","major_comments":[{"comment":"The construction of the Fierz-symmetric renormalisation scheme in Section 2.3 leaves a solution space for the second-generation evanescent coefficients; the paper selects one member of this space (Eqs. (31)–(34)) after imposing conditions 1–3. Appendix D states that only the leading O(N_c^2) term of the renormalised physical matrix elements is independent of the remaining undetermined constant. The NNLO matching coefficients H^(2) and \\tilde H_S^(2) that enter Eq. (84) are computed in this chosen scheme, and the paper does not demonstrate that the subleading 1/N_c parts of these coefficients are also independent of the arbitrary constants (e.g., k and \\tilde k set to zero in Eq. (33)). For N_c = 3 the subleading terms are not parametrically negligible, and the lattice bag parameters used in Eq. (77) are four-dimensional quantities that cannot cancel a perturbative scheme dependence of the Wilson coefficients at order α_s^2. The authors should either prove the full independence of the physical matching coefficients from all allowed evanescent constants, or quantify the numerical shift in ΔΓ_s obtained by varying those constants within the allowed solution space. Without this check, the central claim that the NNLO perturbative uncertainty of the leading 1/m_b term is reduced to the experimental level is not fully established.","section":"Section 2.3, Eqs. (31)–(34), Appendix D"}],"minor_comments":[{"comment":"The phrase \"The calculated NNLO correction reduce\" contains a subject-verb disagreement; it should be \"The calculated NNLO correction reduces\".","section":"Abstract"},{"comment":"The entry for f_Bd is given as \"0.1905 ± 0.0013 MeV\"; to match the f_Bs entry and the text, the unit should be GeV.","section":"Table 2"},{"comment":"In the expression for α_2^(2), the N_H term contains an extra closing parenthesis after the bracket; please correct the typo.","section":"Eq. (54)"},{"comment":"The sentence \"In Eqs. 107 we set mb to unity\" should read \"In Eqs. (107) we set mb to unity\".","section":"Appendix C"},{"comment":"The abstract describes the result as having \"full dependence\" on the quark masses; since the computation is a semi-analytic expansion in z up to order z^10 (with z^50 for some leading terms), I suggest adding the word \"semi-analytic\" to the abstract to avoid overstatement.","section":"Abstract and Section 4"},{"comment":"The statement that the O(ϵ) terms of the third- and fourth-generation evanescent operators \"can be checked\" to drop out of the physical matching coefficients is asserted but not demonstrated in the text; please provide the check or an explicit reference.","section":"Section 2.3"}],"recommendation":"major_revision","confidential_remarks":"This is a high-quality calculation from a group with a strong track record in the field, and the internal cross-checks are impressive. The main issue I see is the unverified subleading 1/N_c scheme dependence of the physical NNLO matching coefficients, which the paper itself acknowledges only in part. A revision that either proves the full independence from the remaining evanescent constants or provides a numerical estimate of the induced shift in ΔΓ_s would resolve my concern. I also recommend that the abstract be reworded slightly to clarify the semi-analytic nature of the z-expansion and the current-current scope of the NNLO corrections."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Bottom line: this is a genuine three-loop QCD matching computation with an unusually strong verification structure, and the new results are real. The main thing to know: it completes the NNLO current-current contribution to Gamma_12 with full charm-mass dependence, and the phenomenology is honestly hedged.\n\nWhat is new: full z = m_c^2/m_b^2 dependence of the three-loop current-current matching coefficients (semi-analytic to z^10 with z^50 pieces), the new two-loop finite renormalisation constants alpha_1^(2), alpha_2^(2), and a Fierz-symmetric |Delta B| = 2 basis. What is well done: two independent computer setups, analytic master integrals checked with pySecDec, reproduction of NLO and fermionic NNLO benchmarks, IR pole cancellations in the matching, and gauge- and regulator-independence checks for alpha_1,2. Those checks matter; they make the three-loop result credible even though I cannot re-derive the integrals myself. The final Delta_Gamma_s = (0.077 +/- 0.016) ps^-1 has an uncertainty dominated by 1/m_b matrix elements, as the paper states, and the scale-dependence plots show genuine improvement with NNLO.\n\nSoft spots, in proportion. (1) The evanescent scheme issue in the stress test is real: Appendix D demonstrates scheme-independence of the leading O(N_c^2) term but not the subleading 1/N_c terms. For N_c = 3 that is a gap in the argument that the computed Wilson coefficients are directly compatible with four-dimensional lattice matrix elements. It is not, however, a demonstrated numerical failure; no evidence is presented that a different allowed solution shifts Delta_Gamma_s by more than the quoted scale uncertainty. A referee should ask for either a proof or a quantitative check. (2) NNLO penguin contributions are absent, so the \"NNLO\" claim is scoped to the current-current sector. The paper says this explicitly, so it is a boundary, not a defect. (3) The z^6 terms of alpha_1,2 come from private communication and the correction to Ref. [20] is not described; this is a minor transparency issue.\n\nWho this is for: B-meson mixing phenomenologists and anyone doing multiloop matching in EFTs. It deserves a serious referee; I would send it to review. I expect the scheme-independence question can be resolved with extra work, and the paper should be publishable once that is addressed.","headline":"A solid, internally well-checked three-loop matching calculation that completes the NNLO current-current program; the residual 1/N_c evanescent-scheme issue is a real gap but not a demonstrated numerical flaw.","tokens_in":46023,"tokens_out":2870,"would_cite":true,"duration_ms":35701,"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":"The paper reports the three-loop (NNLO) QCD matching coefficients for the current-current operator contribution to the $B_q$-$\\bar B_q$ decay matrix with full dependence on $z=m_c^2/m_b^2$, and uses them to predict…","keywords":["B meson mixing","decay width difference","NNLO QCD","three-loop calculation","Fierz symmetry","evanescent operators","heavy quark expansion","CP asymmetry in flavour-specific decays"],"falsifier":"Take the generic evanescent-operator solution of Appendix D instead of the specific choice in Eqs. (31)-(34), recompute the subleading $1/N_c$ terms of the renormalised physical matrix elements at two loops, and check whether the NNLO matching coefficients change by more than the quoted scale uncertainty; a larger shift would show that the Fierz-symmetry conditions do not uniquely fix the scheme. On the experimental side, a future measurement of $\\Delta\\Gamma_s$ with uncertainty below about $0.005\\,\\text{ps}^{-1}$ that falls outside $0.077\\pm0.016\\,\\text{ps}^{-1}$ would rule out the Standard Model prediction presented here.","tokens_in":44728,"feed_emoji":"⚛️","tokens_out":8135,"duration_ms":80200,"temperature":0.7,"pith_summary":"This paper claims that the three-loop (NNLO) QCD corrections induced by two current-current operators to the off-diagonal decay matrix $\\Gamma_{12}$ in $B_q$-$\\bar B_q$ mixing have been computed, with full dependence on $z=m_c^2/m_b^2$ expressed as a semi-analytic expansion to order $z^{10}$ (and to $z^{50}$ for the leading terms). The computation required constructing the $|\\Delta B|=2$ effective theory so that Fierz symmetry is preserved in $D\\neq 4$ dimensions, which fixes the evanescent-operator definitions that enter the renormalisation. The new corrections stabilise the renormalisation-scale dependence and reduce the perturbative uncertainty of the leading $1/m_b$ term in $\\Delta\\Gamma_s$ to the level of the current experimental error. A sympathetic reader would care because $\\Gamma_{12}$ controls the observable width differences and flavour-specific CP asymmetries, and the updated Standard Model predictions sharpen tests of $B$-mixing that do not rely on a global CKM fit.","feed_headline":"Three-loop QCD corrections tighten B_s width-difference prediction","feed_subtitle":"Full charm-mass dependence at NNLO shrinks the perturbative error of the leading 1/m_b term to the experimental level.","key_machinery":"The load-bearing object is the Fierz-symmetric set of evanescent operators in the $|\\Delta B|=2$ theory: operators with extra Dirac matrices (first to fourth generation) that vanish in four dimensions but are needed to renormalise the physical operators $Q$, $Q_S$, and $\\tilde Q_S$ in $D\\neq 4$. Their $O(\\epsilon)$ and $O(\\epsilon^2)$ coefficients are fixed by three conditions: equality of renormalised physical matrix elements with their Fierz transforms at one- and two-loop order, flavour-number independence of the operator definitions, and correct large-$N_c$ counting. The second-generation coefficients are pinned down only up to a solution space, and the paper verifies that the leading $O(N_c^2)$ term of the renormalised physical matrix elements is independent of the remaining constants. A second load-bearing element is the finite renormalisation of $R_0=\\frac12 Q+Q_S+\\tilde Q_S$, whose matrix element must be $1/m_b$ suppressed; the constants $\\alpha_1,\\alpha_2$ are extracted at two loops by separating UV from IR divergences with a gluon mass. These two elements make the matching coefficients IR-finite order by order in $\\alpha_s$ and compatible with four-dimensional lattice matrix elements.","core_discovery":"The central claim is that the three-loop (NNLO) QCD matching coefficients for the contribution of two current-current operators to $\\Gamma^q_{12}$ are now known with full dependence on $z=m_c^2/m_b^2$, provided as a semi-analytic expansion to order $z^{10}$ (and to $z^{50}$ for the leading terms). The calculation is performed in a $|\\Delta B|=2$ effective theory whose evanescent operators are fixed so that Fierz symmetry holds in $D\\neq 4$ dimensions: the renormalised matrix elements of physical operators equal those of their Fierz transforms at one- and two-loop level, the operator definitions do not depend on the number of quark flavours, and the large-$N_c$ limit is correctly reproduced. With the finite renormalisation of $R_0$ enforcing its $1/m_b$ suppression, the matching yields finite NNLO Wilson coefficients. Combined with lattice bag parameters and $1/m_b$ matrix elements, these coefficients give $\\Delta\\Gamma_s=(0.077\\pm0.016)\\,\\text{ps}^{-1}$, $\\Delta\\Gamma_d=(0.00211\\pm0.00045)\\,\\text{ps}^{-1}$, $a_{\\rm fs}^s=(2.28\\pm0.14)\\times10^{-5}$, and $a_{\\rm fs}^d=-(5.21\\pm0.32)\\times10^{-4}$, with the uncertainty of $\\Delta\\Gamma_s$ dominated by the sub-leading terms of the $1/m_b$ expansion.","pith_inferences":["The same Fierz-symmetry construction could be applied to the uncalculated penguin-operator contributions at NNLO; the paper estimates their numerical effect at the central scale is small, but including them may further stabilise the renormalisation-scale dependence.","The solution-space freedom in the second-generation evanescent operators (Eqs. (31)-(34) versus the generic solution in Appendix D) offers a concrete test: recomputing the matching coefficients with different admissible constants would reveal whether the subleading $1/N_c$ terms are truly scheme-independent, as the leading $O(N_c^2)$ term is claimed to be.","Because the semi-analytic expansion to $z^{10}$ is effectively exact for the physical charm-to-bottom mass ratio, the practical bottleneck for $\\Delta\\Gamma_s$ is the lattice determination of the dimension-7 operator matrix elements; improving those inputs will translate directly into a more precise Standard Model prediction."],"forward_implications":["The perturbative uncertainty of the leading $1/m_b$ term in $\\Delta\\Gamma_s$ drops to the level of the current experimental error, so further progress on $\\Delta\\Gamma_s$ now hinges mostly on the matrix elements of the $1/m_b$-suppressed operators.","Because the ratio $\\Delta\\Gamma_q/\\Delta M_q$ is almost independent of $|V_{cb}|$, the NNLO prediction sharpens a probe of new physics that does not require resolving the $|V_{cb}|$ puzzle.","The NNLO corrections reduce the renormalisation-scale variation of $\\Delta\\Gamma_s/\\Delta M_s$ by roughly a factor of two in both the MS and PS schemes (about 15% and 17% in the $[2.1,8.4]$ GeV interval, versus 33% and 26% at NLO).","The Standard Model values $a_{\\rm fs}^s=(2.28\\pm0.14)\\times10^{-5}$ and $a_{\\rm fs}^d=-(5.21\\pm0.32)\\times10^{-4}$ are now stable under scale variation, making a future measurement of $a_{\\rm fs}^d$ a direct test of the CKM-apex constraint from $B$-mixing observables alone."],"supporting_citations":[{"why":"Supplies the NLO current-current matching coefficients exact in $z$ that are extended here to NNLO.","marker":"[11]"},{"why":"Provides the renormalisation of the $|\\Delta B|=2$ theory and the NLO penguin contributions used in the phenomenological analysis.","marker":"[16]"},{"why":"Provides the fermionic NNLO contributions from closed charm loops that are reproduced and used as an input for the $z$-dependent terms.","marker":"[17]"},{"why":"The previous numerical NNLO result for $\\Delta\\Gamma_s$ that this paper reproduces and improves with full charm-mass dependence.","marker":"[20]"},{"why":"The technical companion supplying the integral families, master-integral expansions, and projectors used for the three-loop calculation.","marker":"[21]"},{"why":"The ancillary computer-readable files containing the NNLO matching coefficients that the final predictions are based on.","marker":"[30]"},{"why":"Supplies lattice matrix elements for the leading-power bag parameters and for the $1/m_b$-suppressed operators in the $B_s$ and $B_d$ systems.","marker":"[56]"},{"why":"Supplies the lattice QCD matrix elements for the width difference beyond leading order, which dominate the remaining uncertainty.","marker":"[79]"}],"fun_headline_variants":["NNLO QCD tightens B_s width-difference prediction","Full charm-mass dependence at three loops sharpens B_s mixing","Three-loop QCD narrows B_s decay-width uncertainty","New NNLO result refines B_s width difference","Precise three-loop QCD for B_s width difference"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The result depends on the assumption that the chosen Fierz-symmetry-preserving renormalisation scheme is the one that matches the four-dimensional lattice operator basis: if another admissible choice of the remaining evanescent-operator constants altered the subleading colour-suppressed terms of the NNLO matching coefficients, the quoted prediction would shift at order $\\alpha_s^2$ by more than the stated uncertainty.","fun_headline_variants_meta":{"raw":{"variants":["NNLO QCD tightens B_s width-difference prediction","Full charm-mass dependence at three loops sharpens B_s mixing","Three-loop QCD narrows B_s decay-width uncertainty","New NNLO result refines B_s width difference","Precise three-loop QCD for B_s width difference"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.00088,"raw_usage":{"total_tokens":3893,"prompt_tokens":1123,"completion_tokens":2770,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":739,"completion_tokens_details":{"reasoning_tokens":2687}},"tokens_in":739,"tokens_out":2770,"duration_ms":20892,"temperature":1.0,"reasoning_tokens":2687,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-07T13:03:39.822549+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Take the generic evanescent-operator solution of Appendix D instead of the specific choice in Eqs. (31)-(34), recompute the subleading $1/N_c$ terms of the renormalised physical matrix elements at two loops, and check whether the NNLO matching coefficients change by more than the quoted scale uncertainty; a larger shift would show that the Fierz-symmetry conditions do not uniquely fix the scheme. On the experimental side, a future measurement of $\\Delta\\Gamma_s$ with uncertainty below about $0.005\\,\\text{ps}^{-1}$ that falls outside $0.077\\pm0.016\\,\\text{ps}^{-1}$ would rule out the Standard Model prediction presented here.","supporting_citations":[],"review_version":1}