Probing new physics through $B_s^{*}\to {\mu }^{+}{\mu }^{-}$ decay

We perform a model independent analysis of new physics in $B_s^{*}\to {\mu }^{+}{\mu }^{-}$ decay. We intend to identify new physics operator(s) which can provide large enhancement in the branching ratio of $B_s^{*}\to {\mu }^{+}{\mu }^{-}$ above its standard model prediction. For this, we consider new physics in the form of vector, axial-vector, scalar and pseudoscalar operators. We find that scalar and pseudoscalar operators do not contribute to the branching ratio of $B_s^{*}\to {\mu }^{+}{\mu }^{-}$. We perform a global fit to all relevant $b\to s{\mu }^{+}{\mu }^{-}$ data for different new physics scenarios. For each of these scenarios, we predict $\mathrm{Br}(B_s^{*}\to {\mu }^{+}{\mu }^{-})$. We find that a significant enhancement in $\mathrm{Br}(B_s^{*}\to {\mu }^{+}{\mu }^{-})$ is not allowed by any of these new physics operators. In fact, for all new physics scenarios providing a good fit to the data, the branching ratio of $B_s^{*}\to {\mu }^{+}{\mu }^{-}$ is suppressed as compared to the standard model (SM) value. Hence the present $b\to s{\mu }^{+}{\mu }^{-}$ data indicates that the future measurement of $\mathrm{Br}(B_s^{*}\to {\mu }^{+}{\mu }^{-})$ is expected to be suppressed in comparison to the standard model prediction.