QCD corrections to the Wilson coefficients C9 and C10 in two-Higgs   doublet models
  In this letter we present the analytic results for the two-loop corrections to the Wilson coefficients C_9(mu_W) and C_10(mu_W) in type-I and type-II two-Higgs-doublet models at the matching scale mu_W. These corrections are important ingredients for next-to-next-to-leading logarithmic predictions of various observables related to the decays B -> X_s l^+ l^- in these models. In scenarios with moderate values of tan(beta) neutral Higgs boson contributions can be safely neglected for e,mu. Therefore we concentrate on the contributions mediated by charged Higgs bosons. 
S. Schilling, C. Greub, N. Salzmann and B. Toedtli