$B$-meson decay width up to $1/m_b^3$ corrections within and beyond the Standard Model

Starting from the most general effective Hamiltonian describing non-leptonic $b$-quark decays $b\to q_1 \bar q_2 q_3$, we compute analytic expressions for all matching coefficients of the two-quark operator contributions in the heavy quark expansion~(HQE) of a $B$ meson, up to mass-dimension-six. In addition, we calculate the weak-annihilation contributions, which enter the matching of four-quark operators in the HQE at dimension-six and were previously missing. Our results complete the calculation of beyond Standard Model effects in non-leptonic, tree-level, $b$-quark decays relevant for $B$ meson lifetimes and lifetime ratios such as $τ(B^0_s)/τ(B^0_d)$. Such effects naturally arise in extensions of the Standard Model (SM) that aim to address the observed tensions between experimental measurements and theoretical predictions based on QCD factorisation in several colour-allowed non-leptonic $B$-meson decays. As a by-product of our calculation, we also determine the matching coefficients in the HQE induced by the QCD-penguin operators within the SM, including both the interference between current-current and penguin operators and the contributions quadratic in the penguin operators. Owing to the suppression of the QCD-penguin Wilson coefficients within the SM, these effects are typically regarded as corrections of order $α_s$ and $α_s^2$ in the strong coupling, respectively. Our results reproduce the known expressions at dimension-three and provide new results for the coefficients of the chromomagnetic operator at dimension-five and of the Darwin operator at dimension-six.

💡 Research Summary

The paper presents a comprehensive calculation of the heavy‑quark expansion (HQE) for non‑leptonic b‑quark decays of the form b→q₁ q̄₂ q₃, covering all operator contributions up to mass dimension six. Starting from the most general effective Hamiltonian that includes the standard current‑current operators (Q₁, Q₂) and the QCD penguin operators (Q₃–Q₆), the authors extend the framework to accommodate generic beyond‑Standard‑Model (BSM) Dirac and colour structures.

Within the HQE, the total decay width Γ(B) is expanded as Γ₃ + Γ₅⟨O₅⟩/m_b² + Γ₆⟨O₆⟩/m_b³ + …, where Γ₃ corresponds to the free‑quark decay, Γ₅ to the chromomagnetic operator (dimension‑5), and Γ₆ to the Darwin operator together with four‑quark operators (dimension‑6). While the dimension‑3 coefficients are well known, the dimension‑5 and dimension‑6 coefficients have only been partially available in the literature, especially for the QCD‑penguin sector and for weak‑annihilation (WA) contributions.

The authors compute, at leading order in α_s, analytic expressions for all matching coefficients of the two‑quark operators (the “two‑loop” part) and for the four‑quark operators (the “one‑loop” part) up to dimension six. They provide new results for the chromomagnetic operator coefficient at dimension‑5 and for the Darwin operator coefficient at dimension‑6. The Darwin operator mixes with four‑quark operators; to handle the infrared (IR) divergences arising from soft‑gluon emission off light quarks, the calculation is performed using two independent IR regulators: a small finite light‑quark mass and dimensional regularisation. The agreement between the two methods validates the subtraction of IR singularities.

A major novelty is the inclusion of weak‑annihilation diagrams that contribute to the matching of dimension‑6 four‑quark operators. These diagrams were missing in previous analyses and are essential for a complete BSM treatment because many BSM scenarios (e.g., left‑handed W′ bosons, di‑quark exchanges, two‑Higgs‑doublet models) generate operators with non‑standard colour or chirality structures that feed into WA. The paper supplies the full set of WA matching coefficients for generic Dirac structures.

As a by‑product, the authors also compute the SM penguin‑induced matching coefficients up to dimension‑6, including both the interference between current‑current and penguin operators and the quadratic penguin contributions. Although penguin Wilson coefficients are numerically small in the SM, the authors’ results give the complete α_s and α_s² corrections, reproducing known dimension‑3 results and extending them to higher dimensions.

Phenomenologically, the new coefficients allow a refined prediction of B‑meson lifetimes and, crucially, of lifetime ratios such as τ(B_s⁰)/τ(B_d⁰). The authors show that SU(3)_F‑breaking effects from dimension‑6 operators can shift this ratio by 0.1–0.3 %, a size comparable to current experimental uncertainties (permille level). Moreover, generic BSM contributions entering at dimension‑6 can modify the ratio by up to about 1 % within present bounds, providing a powerful indirect probe of new physics.

In summary, the work completes the HQE matching at order 1/m_b³ for non‑leptonic b‑decays, delivering the missing chromomagnetic, Darwin, and weak‑annihilation pieces, and furnishing the full SM penguin sector up to dimension‑6. These results constitute a solid theoretical foundation for high‑precision lifetime studies and for constraining BSM scenarios through inclusive B‑meson observables. Future extensions could include NLO‑QCD corrections to the new coefficients, a systematic treatment of subleading colour structures, and improved non‑perturbative matrix elements from lattice QCD or HQET sum rules.