We calculate the branching ratios and CP-violating asymmetries for B 0 → ωη, ωη ′ , φη and φη ′ decays in the perturbative QCD (pQCD) factorization approach. The pQCD predictions for the CP-averaged branching ratios are Br(B 0 → ωη) = 2.7 +1.1 −1.0 × 10 −7 , Br(B 0 → ωη ′ ) = 0.75 +0.37 −0.33 × 10 −7 , and Br(B 0 → φη) = 6.3 +3.3 −1.9 × 10 −9 , Br(B 0 → φη ′ ) = 7.3 +3.5 −2.6 × 10 −9 which are consistent with currently available experimental upper limits. The inclusion of the gluonic contribution can change the branching ratios of B → ω(φ)η ′ decays by about 10%. The direct CP-violating asymmetries for B 0 → ωη and ωη ′ decays are generally large in size.

We calculate the CP averaged branching ratios and CP-violating asymmetries for B 0 → ηη, ηη ′ and η ′ η ′ decays by employing the perturbative QCD (pQCD) factorization approach. The pQCD predictions for the CP-averaged branching ratios are Br(B 0 → ηη) ≈ 0.67 × 10 −7 , Br(B 0 → ηη ′ ) ≈ 0.18 × 10 −7 , and Br(B 0 → η ′ η ′ ) ≈ 0.11 × 10 −7 , which are consistent with currently available experimental upper limits. We also predict large CP-violating asymmetries for the considered three decay modes, which can be tested by the future B meson experiments.