Revisiting $B_s \to \mu^+\mu^-$ in the two-Higgs doublet models with $Z_2$ symmetry

Abstract

We revisit the rare leptonic decay $B_s \to \mu^+ \mu^-$ in the two-Higgs doublet models with a softly broken $Z_2$ symmetry, namely type-I, type-II, type-X and type-Y 2HDMs. We have derived the relevant full one-loop Wilson coefficients of the four 2HDMs from the recent calculation in the aligned two-Higgs doublet model by Li, Lu and Pich, which could be mapped to all the four 2HDMs for both large and small $\tan\beta$. It is found that, a new term associated with the soft $Z_2$ symmetry breaking parameter $M$ can be enhanced by $\tan^2\beta$ in the type-II 2HDM, which has not been considered in the literature. Imposing both theoretical and experimental constraints, we have renewed the bounds on the parameter spaces of the four 2HDMs. Different from our previous paper, however, we find that all the four 2HDMs give sizable and similar contributions to $\overline{\mathcal B}(B_s \to \mu^+ \mu^-)$ within the stringently restricted parameter spaces, but very tiny as regards the mass-eigenstate rate asymmetry $\mathcal A_{\Delta\Gamma}$; this makes it unfeasible to discriminate the four types of 2HDM with the correlations between the observables in $B_s \to \mu^+ \mu^-$ decay.