▎ 摘　　要

Heteroatom doped graphene as a single-atom catalyst for oxygen reduction reaction (ORR) has received extensive attention in recent years. In this paper, the ORR activity of defective graphene anchoring single heteroatom (IIIA, IVA, VA, VIA and VIIA) was systematically investigated using a dispersion-corrected density functional theory method. For all of the 34 catalysts, 14 of which were further analyzed, and the Gibbs free energy of each elementary reaction was calculated. According to the scaling relationship between G(OOH*) and G(OH*), we further analyzed the rate-determining step of the remaining 20 catalysts. The results show that when the ORR reaction proceeds in the path O-2 OOH O OH H2O, the reaction energy barriers are lower than 0.8 eV for Te-SV, Sb-DV, Pb-SV, Pb-DV, As-SV, As-DV, B-SV, Sn-SV and N-SV. Our result provides a theoretical basis for further exploration of carbon-based single-atom catalysts for ORR.