Efficient and durable catalysts are crucial for the oxygen evolution reaction (OER). The discovery of the high OER catalytic activity in Ni12P5 has attracted a great deal of attention recently. Herein, the microscopic mechanism of OER on the surface of Ni12P5 is studied using density functional theory calculations (DFT) and ab initio molecular dynamics simulation (AIMD). Our results demonstrate that the H2O molecule is preferentially adsorbed on the P atom instead of on the Ni atom, indicating that the nonmetallic P atom is the active site of the OER reaction. AIMD simulations show that the dissociation of H from the H2O molecule takes place in steps; the hydrogen bond changes from Oa-H⋯Ob to Oa⋯H-Ob, then the hydrogen bond breaks and an H+ is dissociated. In the OER reaction on nickel phosphides, the rate-determining step is the formation of the OOH group and the overpotential of Ni12P5 is the lowest, thus showing enhanced catalytic activity over other nickel phosphides. Moreover, we found that the charge of Ni and P sites has a linear relationship with the adsorption energy of OH and O, which can be utilized to optimize the OER catalyst.
Keywords: Ni12P5; ab initio molecular dynamics; adsorption energy; density functional theory; nonmetallic active site; oxygen evolution reaction.