Reaction modelling of hydrogen evolution on nickel phosphide catalysts: density functional investigation.

Nickel phosphides (NiP), particularly NiP, are promising catalysts for the acidic hydrogen evolution reaction (HER). Using density functional theory (DFT), we model HER at the potential of zero charge (PZC), incorporating solvation effects an explicit water cluster and implicit surrounding solvent. Comparing the Volmer, Tafel, and Heyrovsky steps under saturated hydrogen coverage on NiP(0001) terminations, we find that the NiP (pristine) surface termination prefers the Volmer-Volmer-Tafel (VVT) pathway with activation energy () of 0.57 eV. Conversely, the NiP + 4P (reconstructed) surface favors the Volmer-Heyrovsky (VH) pathway with = 0.60 eV. For the pristine surface termination, the differential gas-phase hydrogen adsorption free energies (Δ) correlate with the Volmer and Tafel step reaction energies, and a linear Bell-Evans-Polanyi relationship for the calculated activation and reaction energies validates the usefulness of the Δ descriptor for the Volmer step under PZC conditions. Nickel atoms play a crucial role in H production on both pristine and reconstructed surfaces, suggesting that modifications of the Ni sites can be used for catalyst design. Our findings highlight the importance of considering surface reconstruction and solvation effects on the HER catalytic performance.