  By employing the QCD factorization approach, we calculate the new physics
contributions to the branching radios of the two-body charmless $ B \to PP$
decays in the framework of the minimal supergravity (mSUGRA) model. Within the
considered parameter space, we find that (a) the supersymmetric (SUSY)
corrections to the Wilson coefficients $C_k$ ($k=3-6$) are very small and can
be neglected safely, but the leading order SUSY contributions to
$C_{7\gamma}(M_W)$ and $C_{8g}(M_W)$ can be rather large and even change the
sign of the corresponding coefficients in the standard model; (b) the possible
SUSY contributions to those penguin-dominated decays in mSUGRA model can be as
large as $30-50%$; (c) for the well measured $B \to K \pi$ decays, the
significant SUSY contributions play an important rule to improve the
consistency of the theoretical predictions with the data; (d) for $B \to K
\eta'$ decays, the theoretical predictions of the corresponding branching
ratios become consistent with the data within one standard deviation after the
inclusion of the large SUSY contributions in the mSUGRA model.
