$\bar{B}^0_s \to (π^0 η^{(*)}, η^{(*)}η^{(*)})$ decays and the effects of next-to-leading order contributions in the perturbative QCD approach
In this paper, we calculate the branching ratios and CP violating asymmetries of the five $\bar{B}^0_s \to (π^0η^{(*)},η^{(*)}η^{(*)})$ decays, by employing the perturbative QCD (pQCD) factorization approach and with the inclusion of all currently known next-to-leading order (NLO) contributions. We find that (a) the NLO contributions can provide about 100% enhancements to the LO pQCD predictions for the decay rates of $\bar{B}_s^0 \to ηη^\prime$ and $η^\prime η^\prime$ decays, but result in small changes to $Br(\bar{B}_s \to π^0 η^{(*)})$ and $Br(\bar{B}_s \to ηη)$; (b) the newly known NLO twist-2 and twist-3 contributions to the relevant form factors can provide about 10% enhancements to the decay rates of the considered decays; (c) for $\bar{B}_s \to π^0 η^{(*)}$ decays, their direct CP-violating asymmetries $\cala_f^{dir}$ could be enhanced significantly by the inclusion of the NLO contributions; and (d) the pQCD predictions for $Br(\bar{B}_s \to ηη^{(*)})$ and $Br(\bar{B}_s \to η^\prime η^\prime)$ can be as large as $4\times 10^{-5}$, which may be measurable at LHCb or the forthcoming super-B experiments.