CP asymmetries in B->φKS and B->η'KS in MSSM

We study the B->φKS and B->η'KS decays in MSSM by calculating hadronic matrix elements of operators with QCD factorization approach and including neutral Higgs boson (NHB) contributions. We calculate the Wilson coefficients of operators including the new operators which are induced by NHB penguins at LO using the MIA with double insertions. We calculate the αs order hadronic matrix elements of the new operators for B->φKS and B->η'KS. It is shown that the recent experimental results on the time-dependent CP asymmetries in B->φKS and B->η'KS, SφKS is negative and Sη'KS is positive, which cannot be explained in SM, can be explained in MSSM if there are flavor non-diagonal squark mass matrix elements of 2nd and 3rd generations whose size satisfies all relevant constraints from known experiments (B->Xsγ, Bs->μ+μ-, B->Xsμ+μ-, B->Xsg, ΔMs, etc.). In particular, we find that one can explain the experimental results with a flavor non-diagonal mass insertion of chirality LL or LR when αs corrections of hadronic matrix elements of operators are included, in contrast with the claim in the literature. At the same time, the branching ratios for the two decays can also be in agreement with experimental measurements