Complete Analysis on the Short Distance Contribution of B_sto ell^+ell^-γ in Standard Model

Using the $B_s$ meson wave function extracted from non-leptonic $B_s$ decays, we evaluate the short distance contribution of rare decays $B_s\to \ell^+\ell^-\gamma(\ell=e,\mu)$ in the standard model, including all the possible diagrams. We focus on the contribution from four-quark operators which are not taken into account properly in previous researches. We found that the contribution is large, leading to the branching ratio of $B_s\to \ell^+\ell^-\gamma$ being nearly enhanced by a factor 3 and up to $1.7\times 10^{-8}$. The predictions for such processes can be tested in the LHC-b and B factories in near future.