  We study the effect of scalar leptoquarks on some rare decays of $B_s$ mesons
involving the quark level transition $b\rightarrow s$ $l^{+} l^{-}.$ In
particular we consider the decays $B_{s} \rightarrow \mu^{+} \mu^{-}$, $\bar
B_d^0 \to X_s \mu^+ \mu^-$ and $B_{s} \rightarrow \phi$ $\mu^{+} \mu^{-}$. The
leptoquark parameter space is constrained using the recently measured branching
ratio of the $B_{s} \rightarrow \mu^{+} \mu^{-}$ process at LHCb and CMS
experiments. Using such parameters we obtain the branching ratio, forward
backward asymmetry and the CP asymmetry parameters in the angular distribution
of $B_s \rightarrow \phi$ $\mu^+ \mu^-$ process.
