Next-to-leading logarithmic QCD contribution of the electromagnetic dipole operator to B¯→Xsγγ

We calculate the set of 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)²/m_b², where p_b, q₁, and q₂ are the momenta of the b quark and two photons. While the (renormalized) virtual corrections are worked exactly for a certain range of s₁ and s₂, we retain in the gluon bremsstrahlung process only the leading power with respect to the (normalized) hadronic mass s₃=(p_b-q₁-q₂)²/m_b² in the underlying triple differential decay width dΓ₇₇/(ds₁ds₂ds₃). The double differential decay width, based on this approximation, is free of infrared and collinear singularities when combining virtual and bremsstrahlung corrections. The corresponding results are obtained analytically. When retaining all powers in s₃, the sum of virtual and bremsstrahlung corrections contains uncanceled 1/ɛ singularities (which are due to collinear photon emission from the s quark), and other concepts, which go beyond perturbation theory, such as parton fragmentation functions of a quark or a gluon into a photon, are needed which is beyond the scope of our paper.