  We study the effects of the MSSM contribution on B -> phi K decays using the
perturbative QCD approach. In this approach, strong phases can be calculated,
so that we can predict the values of CP asymmetries with the MSSM contribution.
We predict a large relative strong phase between the penguin amplitude and the
chromomagnetic penguin amplitude. If there is a new CP violating phase in the
chromomagnetic penguin amplitude, then the CP asymmetries may change
significantly from the SM prediction. We parametrize the new physics
contributions that appear in the Wilson coefficients. We maximize the new
physics parameters up to the point where it is limited by experimental
constraints. In the case of the LR insertion, we find that the direct CP
asymmetries can reach about 85% and the indirect CP asymmetry can reach about
-30%.
