  By employing the QCD factorization approach, we calculate the new physics
contributions to the branching radios of the two-body charmless $ B \to PV$ and
$B \to VV$ decays in the framework of the minimal supergravity (mSUGRA) model.
we choose three typical sets of the mSUGRA input parameters in which the Wilson
coefficient $C_{7\gamma}(m_b)$ can be either SM-like (the case A and C) or has
a flipped-sign (the case B). We found numerically that (a) the SUSY
contributions are always very small for both case A and C; (b) for those
tree-dominated decays, the SUSY contributions in case B are also very small;
(c) for those QCD penguin-dominated decay modes, the SUSY contributions in case
B can be significant, and can provide an enhancement about $30% \sim 260%$ to
the branching ratios of $B \to K^*(\pi,\phi,\rho)$ and $K \phi$ decays, but a
reduction about $30% \sim 80%$ to $ B\to K(\rho, \omega)$ decays; and (d) the
large SUSY contributions in the case B may be masked by the large theoretical
errors dominated by the uncertainty from our ignorance of calculating the
annihilation contributions in the QCD factorization approach.
