Contribution of the Darwin operator to non-leptonic decays of heavy quarks

We compute the Darwin operator contribution ($$ 1/{m}_b^3 $$ 1 / m b 3 correction) to the width of the inclusive non-leptonic decay of a B meson (B+, Bd or Bs), stemming from the quark flavour-changing transition b → $$ {q}_1{\overline{q}}_2{q}_3 $$ q 1 q ¯ 2 q 3 , where q1, q2 = u, c and q3 = d, s. The key ideas of the computation are the local expansion of the quark propagator in the external gluon field including terms with a covariant derivative of the gluon field strength tensor and the standard technique of the Heavy Quark Expansion (HQE). We confirm the previously known expressions of the $$ 1/{m}_b^3 $$ 1 / m b 3 contributions to the semi-leptonic decay b → $$ {q}_1\mathrm{\ell}{\overline{\nu}}_{\mathrm{\ell}} $$ q 1 ℓ ν ¯ ℓ , with ℓ = e, μ, τ and of the $$ 1/{m}_b^2 $$ 1 / m b 2 contributions to the non-leptonic modes. We find that this new term can give a sizeable correction of about −4 % to the non-leptonic decay width of a B meson. For Bd and Bs mesons this turns out to be the dominant correction to the free b-quark decay, while for the B+ meson the Darwin term gives the second most important correction — roughly 1/2 to 1/3 of the phase space enhanced Pauli interference contribution. Due to the tiny experimental uncertainties in lifetime measurements the incorporation of the Darwin term contribution is crucial for precision tests of the Standard Model.