Next-to-leading logarithmic QCD contribution of the electromagnetic dipole operator to B¯→Xsγγ with a massive strange quark

We calculate the O(α_s) corrections to the double differential decay width dΓ₇₇/(ds₁ds₂) for the process B¯→X_sγγ, originating from diagrams involving the electromagnetic dipole operator O₇. The kinematical variables s₁ and s₂ are defined as s_i=(p_b-q_i)²/mb_{2, where pb, q1, q2 are the momenta of the b quark and two photons. We introduce a nonzero mass ms for the strange quark to regulate configurations where the gluon or one of the photons become collinear with the strange quark and retain terms which are logarithmic in ms, while discarding terms which go to zero in the limit ms→0. When combining virtual and bremsstrahlung corrections, the infrared and collinear singularities induced by soft and/or collinear gluons drop out. By our cuts the photons do not become soft, but one of them can become collinear with the strange quark. This implies that in the final result a single logarithm of ms survives. In principle, the configurations with collinear photon emission could be treated using fragmentation functions. In a related work we find that similar results can be obtained when simply interpreting ms appearing in the final result as a constituent mass. We do so in the present paper and vary ms between 400 and 600 MeV in the numerics. This work extends a previous paper by us, where only the leading power terms with respect to the (normalized) hadronic mass s3}=(p_b-q₁-q₂)²/mb2 were taken into account in the underlying triple differential decay width dΓ₇₇/(ds₁ds₂ds₃).