A new puzzle in non-leptonic B decays

We build a set of new observables using closely related non-leptonic penguin-mediated B d and B s decays: B_d,s→K^∗ 0K^∗ 0,B_d,s→K^0K^0,B_d,s→K^0K^∗ 0 and B_d,s→K^0K^∗ 0 together with their CP conjugate partners. These optimised observables are designed to reduce hadronic uncertainties, mainly coming from form factors and power-suppressed contributions, and thus maximize their sensitivity to New Physics (NP). The deviations observed with respect to the Standard Model (SM) in the ratios of branching ratios of B_d,s→K^∗ 0K^∗ 0(2.6σ) and B_d,s→K^0K^0(2.4σ) can be explained by simple NP scenarios involving the Wilson coefficients 𝒞 4 and 𝒞 6 associated with QCD penguin operators and the coefficient 𝒞 8 g of the chromomagnetic operator. The optimised observables for B_d,s→K^0K^∗ 0 and B_d,s→K^0K^∗ 0 show distinctive patterns of deviations with respect to their SM predictions under these NP scenarios. The pattern of deviations of individual branching ratios, though affected by significant hadronic uncertainties, suggests that NP is needed both in b → d and b → s transitions. We provide the regions for the Wilson coefficients consistent with both optimised observables and individual branching ratios. The NP scenarios considered to explain the deviations of B_d,s→K^∗ 0K^∗ 0 and B_d,s→K^0K^0 can yield deviations up to an order of magnitude among the observables that we introduced for B_d,s→K^0K^∗ 0 and B_d,s→K^0K^∗ 0 . Probing these new observables experimentally by improving the measurements of individual branching ratios of penguin-mediated decays would confirm the consistency of the deviations already observed and provide a highly valuable hint of NP in the non-leptonic sector.