Theoretical framework for lattice QCD computations of B→Kℓ+ℓ- and Bs→ℓ+ℓ-γ decays rates, including contributions from charming penguin diagrams

We develop a strategy for computing the B - Kl+l- and BS - gamma l+l- decay amplitudes using lattice QCD (where lt are charged leptons). We focus on those terms which contain complex contributions to the amplitude, due to on-shell intermediate states propagating between the weak operator and electromagnetic current(s). Such terms, which are generally estimated using model calculations and represent significant uncertainties in the phenomenological predictions for these decays, cannot be computed using standard lattice QCD techniques. It has recently been shown that such contributions can be computed using spectral density methods and our proposed strategy, which we discuss in detail, is built on this approach. The complex contributions include the "charming penguins" [matrix elements of the current-current operators 2 defined in Eq. (6) below], in which the charm-quark loop can propagate long distances, particularly close to the region of charmonium resonances. They also include the contributions from the chromomagnetic operator [O8 in standard notation, defined in Eq. (8) below]. We discuss the renormalization of the ultraviolet divergences, and in particular those which arise due to "contact" terms, and explain how those which appear as inverse powers of the lattice spacing can be subtracted nonperturbatively. We apply the spectral density methods in an instructive exploratory computation of charm-quark loop [the diagram in Fig. 1(a) below] and discuss the prospects and strategies for the reliable determination of the amplitudes in future dedicated computations; computations which are, however, beyond the scope of the present paper.