By employing the perturbative QCD (PQCD) factorization approach, we made a systematic investigation for the CP-averaged branching ratios and the CP-violating asymmetries of the thirteenB 0 s → P P decays ( here P = (π, K, η, η ′ )) with the inclusion of all currently known next-to-leading order (NLO) contributions, and compared our results with the measured values or the theoretical predictions from other different approaches. We focused on the examination of the effects of the NLO contributions and found the following points:decays, the NLO contributions can provide a (40 − 150)% enhancement or a 20% reduction to the corresponding leading order (LO) PQCD predictions for their decay rates and result in a much better agreement between the PQCD predictions and the experimental measurements; (b) for the pure annihilation decaȳ B 0 s → π + π − , the PQCD prediction still remain consistent with the data after the inclusion of the small NLO reduction; (c) the PQCD predictions for the ratiobecome agree very well with the measured ones after the inclusion of a 40% NLO reduction; (d) for B 0 s → K + K − andB 0 s → π − K + decays, the NLO PQCD predictions for their CP-violating asymmetries do agree very well with the measured values in both the sign and the magnitude; and (e) for allB 0 s → P P decays, we also compared our results with those obtained in the QCD factorization approach and soft-collinear effective theory and discussed their similarities and differences.