LAUSR

In this study, the oxygen reduction reaction (ORR) process of dual-metal active site catalysts (FeMN6-Gra, M = Mn, Ni, Co, or Cu) mediated by p-block elements was investigated using density functional theory calculations. The obtained results demonstrate that, in most cases, the B-doped FeMN6-Gra (M = Mn, Ni, Co, or Cu) catalysts exhibit higher catalytic performance than their undoped counterparts. Among the investigated catalysts, FeNiN6-Gra doping by B modulates the adsorption strength of the metal center on the oxygen-containing intermediates, showing the largest increase in the onset potential (from 0.66 to 0.94 V). Importantly, we found a new law that B-doping affects the total charge of the metal adsorption site and the four surrounding N atoms and that there is a linear relationship between the total charge and the Gibbs free energy. Transition state analysis shows that the energy barrier of the thermodynamic rate-determining step (*OH hydrogenation to H2O) in the FeNiN6B1-Gra-catalyzed ORR process is 0.17 eV, which is smaller than that of the FeNiN6-Gra-catalyzed process (0.28 eV). Overall, the results demonstrate that B-doping can improve the activity of FeMN6-Gra catalysts and provide a new method for the future development of efficient electrocatalysts.