Active Role of Phosphorus in the Hydrogen Evolving Activity of Nickel Phosphide (0001) Surfaces

Optimizing catalysts for the hydrogen evolution reaction (HER) is a critical step toward the efficient production of H 2 (g) fuel from water. It has been demonstrated experimentally that transition-metal phosphides, specifically nickel phosphides Ni 2 P and Ni 5 P 4, efficiently catalyze the HER at a small fraction of the cost of archetypal Pt-based electrocatalysts. However, the HER mechanism on nickel phosphides remains unclear. We explore, through density functional theory with thermodynamics, the aqueous reconstructions of Ni 2 P(0001) and Ni 5 P 4 (0001)/(0001̅), and we find that the surface P content on Ni 2 P(0001) depends on the applied potential, which has not been considered previously. At −0.21 V ≥ U ≥ −0.36 V versus the standard hydrogen electrode and pH = 0, a PH x -enriched Ni 3 P 2 termination of Ni 2 P(0001) is found to be most stable, consistent with its P-rich ultrahigh-vacuum reconstructions. Above and below this potential range, the stoichiometric Ni 3 P 2 surface is instead passivated by H at the Ni 3 -hollow sites. On the other hand, Ni 5 P 4 (0001̅) does not favor additional P. Instead, the Ni 4 P 3 bulk termination of Ni 5 P 4 (0001̅) is passivated by H at both the Ni 3 and P 3 -hollow sites. We also found that the most HER-active surfaces are Ni 3 P 2 +P+(7/3)H of Ni 2 P(0001) and Ni 4 P 3 +4H of Ni 5 P 4 (0001̅) due to weak H adsorption at P catalytic sites, in contrast with other computational investigations that propose Ni as or part of the active site. By looking at viable catalytic cycles for HER on the stable reconstructed surfaces, and calculating the reaction free energies of the associated elementary steps, we calculate that the overpotential on the Ni 4 P 3 +4H surface of Ni 5 P 4 (0001̅) (−0.16 V) is lower than that of the Ni 3 P 2 +P+(7/3)H surface of Ni 2 P(0001) (−0.21 V). This is due to the abundance of P 3 -hollow sites on Ni 5 P 4 and the limited surface stability of the P-enriched Ni 2 P(0001) surface phase. The trend in the calculated catalytic overpotentials, and the potential-dependent bulk and surface stabilities explain why the nickel phosphides studied here perform almost as well as Pt, and why Ni 5 P 4 is more active than Ni 2 P toward HER, as is found in the experimental literature. This study emphasizes the importance of considering aqueous surface stability in predicting the HER-active sites, mechanism, and overpotential, and highlights the primary role of P in HER catalysis on transition-metal phosphides.