▎ 摘　　要

In this study, we investigate the oxygen reduction reaction (ORR) reactivity of nitrogen-doped graphene by using density functional theory (DFT), a computational quantum mechanical technique. Four doping configurations and five models were comprehensively studied: quaternary nitrogen (NQ), pyrrolic nitrogen (N5), two forms of pyridinic nitrogen (N6, N6nH) and three-pyridinic nitrogen (3N6). Models for possible sites during each step of ORR were set up and visualized to provide a platform to calculate the free energy of the reaction pathway to determine the suitability of each doping scenario. Associative mechanisms were displayed by all models except N5, which showed dissociative mechanism. The calculated free energy pathways demonstrate that the ranking of the reactivity for ORR by different nitrogen configurations from most reactive to least reactive is N6, NQ, N6nH, 3N6, and N5. Spin density and charge density aid in describing levels of reactivity.