  Taking account of the constraint from radiative two body decays
$l_i\rightarrow l_j\gamma$, we investigate the Lepton Flavor Violation decays
$B^0_q\rightarrow\bar{l_l}l_k$ in the framework of the minimal extension of the
Standard Model with one neutral singlet scalar. The couplings $C_{e\mu}$,
$C_{e\tau}$ and $C_{\mu\tau}$ between the different generation leptons and
scalar $S^0$ are constrained by the current bounds of $l_i\rightarrow
l_j\gamma$. The numerical results show that the theoretical prediction of
$B^0_q\rightarrow\bar{l_l}l_k$ strongly depend on the couplings $C_{qb}$
($q=d\;or\; s$) between down type quarks and new scalar. The contributions from
couplings $C_{uc}$, $C_{ut}$ and $C_{ct}$ between up type quark and new scalar
are less dominant.
