$B\to X_s\tau^+\tau^-$ in a two Higgs doublet model
  The inclusive rate and forward-backward asymmetry of dilepton angular distribution for a B-meson to decay to strange hadronic final states and a $\tau^+\tau^-$ pair in a two Higgs doublet model are computed. In particular, contributions of neutral Higgs bosons to the decay are included. QCD corrections to the effective Hamiltonian for $B\rightarrow X_s\tau^+\tau^-$ are calculated using the leading logarithmic approximation. 
Yuan-Ben Dai, Chao-Shang Huang, Han-Wen Huang