Starting from the most general effective |∆B| = 1 Hamiltonian describing non-leptonic b-quark decays \( b\to {q}_1{\overline{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 (BSM) effects in non-leptonic, tree-level, b-quark decays relevant for B meson lifetimes and lifetime ratios such as \( \tau \left({B}_s^0\right)/\tau \left({B}_d^0\right) \) . Such BSM contributions naturally arise in generic 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 \( {\alpha}_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.
