Nonfactorizable charming-loop contribution to FCNC Bs→γ l+l- decay

We present the first theoretical calculation of nonfactorizable charm-quark loop contributions to the B_s→γ l⁺l^- amplitude. We calculate the relevant form factors, H_A^{,V NF}(k^'2,k²), and provide convenient parametrizations of our results in the form of fit functions of two variables, k^'2 and k², applicable in the region below hadron resonances, k^'2<M_J^{/ψ 2} and k²<M_ϕ². We report that factorizable and nonfactorizable charm contributions to the B_s→γ l⁺l^- amplitude have opposite signs. To compare the charm and the top contributions, it is convenient to express nonfactorizable charming loop contribution as a nonuniversal (i.e., dependent on the reaction) q²-dependent correction Δ^NFC₇(q²) to the Wilson coefficient C₇. For the B_s→γ l⁺l^- amplitude, the correction is found to be positive, Δ^NFC₇(q²)/C₇>0 .