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

We compute the Darwin operator contribution (1/mb3 correction) to the width of the inclusive non-leptonic decay of a B meson (B+, B-d or B-s), stemming from the quark flavour-changing transition b -> q1q2q3, where q(1), q(2) = u, c and q(3) = 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/mb3 contributions to the semi-leptonic decay b -> q1lnu l, with l = e, mu, tau and of the 1/mb2 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 B-d and B-s 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.