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

Nickel phosphides (NixPy), particularly Ni2P, 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 via an explicit water cluster and implicit surrounding solvent. Comparing the Volmer, Tafel, and Heyrovsky steps under saturated hydrogen coverage on Ni2P(0001) terminations, we find that the Ni3P2 (pristine) surface termination prefers the Volmer-Volmer-Tafel (VVT) pathway with activation energy (Ea) of 0.57 eV. Conversely, the Ni3P2 + 4P (reconstructed) surface favors the Volmer-Heyrovsky (VH) pathway with Ea = 0.60 eV. For the pristine surface termination, the differential gas-phase hydrogen adsorption free energies (Delta Gdiff) 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 Delta Gdiff descriptor for the Volmer step under PZC conditions. Nickel atoms play a crucial role in H2 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. Volmer-Tafel (left) and Volmer-Heyrovsky (right) reaction cycles of hydrogen evolution reaction on Ni3P2 and it's reconstructed surface.