Computational Study Explains How Intermetallic Compounds Catalyze Hydrogen Evolution Reactions

Density functional theory calculations analyzing 145 ordered L12-type intermetallic compounds have revealed that transition metal dopants enhance hydrogen evolution reaction activity primarily by electronically tuning surrounding noble metal host atoms rather than acting as direct adsorption centers. In a study published in Catal, researchers constructed periodic superlattice models composed of noble matrix metals—specifically rhodium, palladium, iridium, and platinum, alongside copper—paired with 29 distinct transition metals across the 3d, 4d, and 5d series to evaluate intrinsic catalytic mechanisms. The computational investigations demonstrated that 70.9 percent of the thermodynamically optimal active sites on these intermetallic surfaces are dominated by matrix noble metal atoms. By examining d-orbital hybridization and metal-metal bonding, the theoretical models explain the underlying origin of the volcano-shaped relationship between hydrogen adsorption free energy and transition metal placement across the periodic table, providing concrete physical principles for designing efficient electrochemical catalysts with lower noble metal contents.

The research addresses pressing challenges in sustainable hydrogen energy production, where water electrolysis serves as a scalable pathway for clean fuel generation but remains constrained by the extreme scarcity and high cost of platinum-group metals. While pure platinum benchmarks commercial hydrogen evolution electrocatalysis, ordered L12-type intermetallic compounds offer a structural architecture to reduce precious metal consumption while maintaining high reaction rates. These intermetallics adopt a face-centered cubic superlattice possessing the stoichiometric formula M3TM, wherein matrix metal atoms occupy face-centered positions while transition metal dopants reside at corner vertex sites. Unlike lower-symmetry intermetallic phases such as B2 or L10 that suffer from coordination distortion, the perfect octahedral coordination around vertex transition metals in L12 lattices induces crystal field splitting into t2g and eg states. Optimal spatial overlap allows t2g orbitals to form strong covalent hybridization with neighboring atoms, imparting thermodynamic stability and enabling fine-grained electronic modulation across adjacent surface sites.

To evaluate catalyst stability and structural trends across the 145 modeled compounds, the authors systematically evaluated formation energies alongside equilibrium lattice parameters. Formation energy calculations identified a pronounced periodic volcano trend correlated with the position of the transition metal dopant within the periodic table, where compounds containing transition metals in the middle of each period displayed substantially higher formation energies than early or late combinations. Conversely, copper-based intermetallics exhibited elevated formation energies across all series because copper’s completely filled 3d electronic shell restricts interfacial charge-transfer capacity. Lattice constant analyses demonstrated a regular structural evolution that first decreased and subsequently increased across each transition metal period, tracking variations in dopant atomic radii. When examining host matrix metals longitudinally, unit cell lattice constants followed the descending hierarchy of platinum, palladium, iridium, rhodium, and copper, matching precisely the measured atomic radii of these elements, ranging from 1.30 angstroms for platinum to 1.17 angstroms for copper.

Atomic size discrepancies between constituent species inevitably induce lattice distortion and surface strain that directly alter the electronic structure of catalytic active sites. The researchers resolved surface distortions into two complementary strain components, designated epsilon(H0) and epsilon(F0), corresponding to distinct adsorption regions. Because geometric constraints enforce continuity across the face-centered cubic lattice, local expansion around one site is accompanied by contraction in the adjacent sublattice, establishing a strictly reciprocal relationship. Using a defined strain formula based on distorted vacant atom spacing relative to the bulk lattice parameter, calculations established that rhodium-, palladium-, iridium-, and platinum-based intermetallics display compressive strain in the H0 surface region, whereas copper-based systems show tensile strain. This divergence occurs because the subsurface H0 layer is dominated by transition metal atoms while copper has the smallest matrix atomic radius at 1.17 angstroms. When dopants feature unusually large atomic radii—such as yttrium, cadmium, mercury, and gold—they induce universal tensile strain and structural mismatches exceeding 30 to 48 percent.

Surface catalytic activity was evaluated across eight potential active sites designated T0, T1, B0, B1, H0, H1, F0, and F1, representing top, bridge, hexagonal close-packed, and face-centered cubic hollow configurations containing varying dopant counts. Initial proton placements followed by structural relaxation revealed that top and hollow sites maintain high thermodynamic stability, retaining adsorbed hydrogen atoms in place. In sharp contrast, bridge sites such as B0 and B1 proved largely metastable, with adsorbed hydrogen atoms migrating spontaneously toward nearby hollow sites where stronger atomic bonds and lower total energies are achieved. Introducing transition metals with nearly empty or completely filled d orbitals substantially elevated the occurrence of metastable adsorption sites across the intermetallic surfaces, linked to lower dopant electronegativity and localized electron transfer toward host metal atoms. Notably, the top T0 site on most copper-based intermetallics exhibited metastability, causing hydrogen atoms to migrate to surrounding F0 hollow positions with energy reductions of 0.89 to 1.16 electron volts.

Projected partial density of states calculations detailed how orbital coupling between matrix and dopant atoms systematically drives electronic modulation. Across each transition metal period, d orbitals of matrix atoms in rhodium, palladium, iridium, and platinum systems initially undergo an energy decrease, or redshift, accompanied by peak broadening reflecting orbital delocalization, followed subsequently by a blueshift. This progression generates a characteristic volcano-shaped trajectory in the d-band center spanning approximately 1 electron volt. The investigation delineated three distinct electronic mechanisms: transition metals with low d-electron counts like scandium donate charge, upshifting the d-band center; fully occupied dopants like zinc induce localized orbital states near the Fermi level that similarly upshift the center; and intermediate transition metals like cobalt generate covalent orbital hybridization that downshifts the d-band center toward favorable catalytic levels. Meanwhile, copper-based systems exhibited an alternating redshift-blueshift-redshift profile attributable to copper’s fully filled 3d10 configuration, while dopant d-band centers decreased monotonically across each period.

Proton adsorption free energy calculations across all 145 intermetallics revealed a periodic volcano response for noble metal matrices and an inverted volcano profile for copper-based systems, directly mirroring the Sabatier principle for optimal catalysis. By applying the criterion where the most active site minimizes absolute hydrogen adsorption free energy, the authors constructed volcano plots relating free energy to reaction overpotentials. Adopting pure platinum as a benchmark with a reference free energy of -0.09 electron volts, the research team established an activity window spanning positive and negative 0.09 electron volts around zero. This performance boundary screened 74 high-performing intermetallic compounds. Matrix metal d-band centers fluctuated within a narrow range of 0.5 to 1.5 electron volts, whereas transition metal d-band centers spanned 10 electron volts. Consequently, 70.9 percent of the optimal sites were situated on matrix metal sublattices, confirming that transition metal additives act primarily as electronic regulators rather than direct catalytic reaction centers.

To identify practical catalyst formulations combining superior reaction kinetics with physical durability, the authors integrated hydrogen adsorption free energy values with negative formation energies, which characterize thermal stability. Applying these dual criteria narrowed the field from 74 active candidates down to 14 promising intermetallic compounds predicted to outperform conventional platinum catalysts in hydrogen evolution efficiency. Among the most prominent candidates identified in the screening were Pt3Ti, which exhibited an adsorption free energy of -4.58 × 10⁻4 electron volts, Pd3Zr with an adsorption energy of -0.001 electron volts, and Rh3Cr with an adsorption value of -0.008 electron volts. Structural analysis of these leading compounds demonstrated that the single-atom top site, T0, serves as the dominant catalytic center across numerous intermetallic lattices. By tuning the electronic environment of the matrix noble metal through strategic transition metal pairing, these compounds achieve balanced proton adsorption and desorption kinetics while lowering precious metal requirements.

The study also evaluated the applicability of classical d-band center theory, which predicts that metal-adsorbate binding strengthens as the d-band center approaches the Fermi level. Evaluating linear correlations confirmed that monometallic T0 sites and multiatomic hollow configurations like H0, F0, H1, and F1 conform closely to d-band center predictions after d-orbital weighting normalization. However, a significant breakdown occurred at mixed T1 adsorption sites when paired with transition metals possessing less than half-filled d shells, such as scandium, titanium, and vanadium with d1 through d3 configurations. In these cases, the d-band center resides above the Fermi level, causing severe deviations in adsorption free energy and inferior catalytic activity. This discrepancy arises because unoccupied transition metal d orbitals fail to adequately populate antibonding states with hydrogen 1s orbitals, weakening chemical adsorption. Conversely, when the d-band center lies below the Fermi level—as seen with d4 through d10 dopants—the linear correlation recovers completely.

All first-principles calculations were performed using the plane-wave pseudopotential density functional theory framework implemented within the CASTEP software package, utilizing the Perdew-Burke-Ernzerhof generalized gradient approximation for electron exchange and correlation. Structural models were constructed by cleaving four-layer (111) surface slabs from 4 × 4 × 1 supercells, incorporating a 20-angstrom vacuum buffer along the perpendicular z-axis to eliminate unphysical interactions between periodic images. Optimization protocols employed ultrasoft pseudopotentials with a 440 electron volt kinetic energy cutoff and a 3 × 2 × 1 Monkhorst-Pack k-point grid, allowing the top two atomic layers to relax fully while holding bottom layers fixed until Hellmann-Feynman forces dropped below 0.01 electron volts per angstrom and total energy converged within 5 × 10⁻⁶ electron volts per atom. Adsorption free energies incorporated zero-point vibrational energy and entropy corrections. Overall, the resulting framework elucidates how atomic ordering, surface lattice strain, and orbital hybridization regulate hydrogen evolution, guiding experimental synthesis of durable, noble-metal-efficient electrocatalysts.

Subject of Research: Theoretical mechanisms of hydrogen evolution catalysis in transition metal-doped intermetallic compounds

Article Title: Theoretical explorations of hydrogen evolution reaction mechanisms in transition metal-doped intermetallic compounds

Article References: Kou, J., Sun, M., Yu, H., Wu, X., Xing, Z., Liang, S., & Huang, B. (2025). Theoretical explorations of hydrogen evolution reaction mechanisms in transition metal-doped intermetallic compounds. Catal, 1(1), Article 3. https://doi.org/10.1007/s44422-025-00003-9

Image Credits: AI Generated

DOI: 10.1007/s44422-025-00003-9

Keywords: hydrogen evolution reaction, intermetallic compounds, density functional theory, electrocatalysis, d-band center theory, transition metals, noble metal catalysts, surface lattice strain, Theoretical, explorations, hydrogen, evolution