Theoretical Insights into Hydrogen Production from Formic Acid Catalyzed by Pt-Group Single-Atom Catalysts

TL;DR

This study uses theoretical calculations to explore how single-atom catalysts made of platinum-group metals can efficiently produce hydrogen from formic acid.

Contribution

The paper introduces a new electronic descriptor for designing selective single-atom catalysts based on d-band center positions and metal-support interactions.

Findings

Pd and Au SACs show significantly lower energy barriers for dehydrogenation compared to dehydration.

Rh SACs have nearly equal energy barriers for both reaction pathways, showing poor selectivity.

Metal-support interactions and charge transfer strongly influence catalytic performance and stability.

Abstract

The rational development of single-atom catalysts (SACs) for selective formic acid dehydrogenation (FAD) requires an atomic-scale understanding of metal–support interactions and electronic modulation. In this study, spin-polarized density functional theory (DFT) calculations were performed to systematically examine platinum-group SACs anchored on graphitic carbon nitride (g-C3N4). The findings reveal that Pd and Au SACs exhibit superior selectivity toward the dehydrogenation pathway, lowering the free energy barrier by 1.42 eV and 1.39 eV, respectively, compared to the competing dehydration route. Conversely, Rh SACs demonstrate limited selectivity due to nearly equivalent energy barriers for both reaction pathways. Stability assessments indicate robust metal–support interactions driven by d–p orbital hybridization, while a linear correlation is established between the d-band center position relative to the Fermi level and catalytic selectivity. Additionally, charge transfer (ranging from 0.029 to 0.467 e) substantially modulates the electronic structure of the active sites. These insights define a key electronic descriptor for SAC design and offer a mechanistic framework for optimizing selective hydrogen production.

Figures

Full text

1. Introduction

Energy is the cornerstone of modern civilization, driving both economic development and technological advancement. However, the extensive reliance on fossil fuels has led to an energy crisis and severe environmental degradation, necessitating the development of sustainable alternatives [1]. Hydrogen, with its high gravimetric energy density, abundance, and environmentally benign combustion byproduct (water), has emerged as a promising energy carrier [2]. Despite these advantages, practical limitations in hydrogen storage and transportation hinder its widespread implementation [3].

Conventional storage methods, including high-pressure and cryogenic systems, pose significant safety risks [4]. In contrast, chemical hydrogen storage via liquid organic hydrogen carriers (LOHCs) offers a safer and more efficient alternative [5], enabling reversible hydrogen release under mild catalytic conditions. Among LOHCs, formic acid is particularly attractive due to its high hydrogen content (4.4 wt%), ambient stability, low toxicity, and non-flammability [6]. Furthermore, formic acid can be synthesized via photocatalytic or electrocatalytic CO_2_ reduction [7], contributing to a sustainable carbon-neutral cycle.

Formic acid decomposition follows two competing pathways [8]: dehydrogenation (HCOOH → H_2_ + CO_2_, ΔG = −32.8 kJ/mol) and dehydration (HCOOH → H_2_O + CO, ΔG = −14.9 kJ/mol). Selectively promoting the dehydrogenation route is critical for high-purity hydrogen production, as CO from the dehydration pathway poisons Pt-based formic acid dehydrogenation (FAD) catalysts [9]. This challenge has driven extensive research efforts to develop catalysts that achieve near-complete dehydrogenation under mild conditions.

Since Coffey et al. first reported the catalytic dehydrogenation of formic acid in the 1960s [10], a multitude of numerous homogeneous catalysts, particularly those based on Ir and Ru, have demonstrated remarkable activity in this process. However, these systems suffer from drawbacks such as high ligand costs, poor stability, and complex recovery processes, limiting their large-scale application [11]. In contrast, heterogeneous catalysts, particularly those based on noble metals such as Au, Pt, and Pd, offer advantages in recyclability and scalability [12,13]. However, each metal presents challenges: Pd is susceptible to CO poisoning, Au exhibits low intrinsic activity at low temperatures, and non-precious metals, while stable, generally lack selectivity. To enhance catalytic performance, strategies such as alloying (e.g., Ag, Cu, Cr, Mn, Co) [14] and support engineering (e.g., activated carbon, carbon nanotubes, metal–organic frameworks, TiO_2_) [15] have been explored. These approaches aim to optimize noble metal utilization while balancing activity and stability. However, challenges such as structural degradation and mass transport limitations persist [16]. The scarcity and cost of noble metals have further spurred interest in non-precious metal catalysts, including Co, Ni, and Cu-based systems [17]. While these catalysts demonstrate long-term stability, they often exhibit poor selectivity, sluggish kinetics, and a tendency toward dehydration side reactions [18]. Consequently, the development of efficient formic acid dehydrogenation (FAD) catalysts remains an ongoing challenge.

Recent advancements in single-atom catalysts (SACs) have shown significant promise by maximizing atomic efficiency and enhancing selectivity. Bulusev et al. [19] synthesized single Pd atoms supported on N-doped carbon, achieving excellent catalytic stability and activity at 60 °C. Density functional theory (DFT) calculations indicated that Pd atoms preferentially bind to pyridinic N-edge sites, enhancing catalytic performance. Han et al. [20] developed Pt-modified Te nanowires that achieved 100% selective hydrogen production, while Yang et al. [20] demonstrated that single-atom Pt/TiN catalysts effectively modulate the activity and selectivity of small organic molecules, favoring the direct dehydrogenation pathway in formic acid decomposition. Additionally, Shi et al. [20] reported a non-precious metal Co-based SAC with a hydrogen production rate of 1403.8 mL g⁻^1^ h⁻^1^, which is 15 times higher than that of commercial Pd/C catalysts.

Despite these advances, the high surface energy of single atoms and atomic clusters often leads to aggregation, resulting in catalyst deactivation. Thus, the selection of appropriate support material is crucial for stabilizing isolated metal atoms and preventing agglomeration. Among potential supports, graphitic carbon nitride (g-C_3_N_4_) has emerged as an excellent candidate due to its structural stability, unique electronic properties, and high porosity, which facilitate strong metal–support interactions [21]. SACs based on Pt-group metals (Pt, Pd, Ru) [22] have demonstrated remarkable catalytic activity and H_2_ selectivity, attributed to their electronic structure and well-defined metal–support interactions.

In this study, we investigate Pt-group SACs (M = Ru, Rh, Pd, Ag, Ir, Pt, Au) anchored on g-C_3_N_4_ using spin-polarized DFT calculations. By correlating d-band center positions, charge redistribution, and reaction energetics, we identify Pd and Au SACs as promising catalysts that favor dehydrogenation over dehydration, thereby suppressing CO formation. This work not only reveals a linear correlation relationship between the d-band center and selectivity but also provides valuable insights for the rational design of efficient SACs for hydrogen production.

2. Methods

Density functional theory (DFT) calculations were carried out using the hybrid Gaussian and plane wave (GPW) method [23] within the QUICKSTEP module of the CP2K 2023.1 package [24]. We utilized the double-zeta valence plus polarization (DZVP-MOLOPT-SR-GTH) basis set along with the Goedecker–Teter–Hutter (GTH) pseudopotential [25,26,27] for calculation. Electron density was expanded in an auxiliary plane wave basis set with a cutoff energy of 500 Ry. Exchange–correlation interactions were described using the spin-polarized Perdew–Burke–Ernzerhof (PBE) functional [28] under the generalized gradient approximation (GGA), with dispersion corrections applied via the DFT-D3 method of Grimme et al. [29].

The Kohn–Sham equations were solved using conventional matrix diagonalization, with Fermi–Dirac smearing (electronic temperature: 300 K) applied to enhance convergence. To improve computational efficiency, diagonalization routines from the ELPA library were employed [29]. Geometry optimizations were performed with energy and force convergence thresholds of 1.0 × 10^−6^ eV and 3.0 × 10⁻^4^ eV Å^−1^, respectively, while static calculations employed a tighter convergence criterion of 1.0 × 10^−8^ eV.

Brillouin zone sampling was achieved using a 3 × 3 × 1 Monkhorst-Pack k-point mesh for both g-C_3_N_4_ and M@g-C_3_N_4_ SACs. A 20 Å vacuum layer was introduced to eliminate spurious periodic interactions. Thermodynamic corrections and density of states (DOS) analyses were conducted using the Shermo v.2.6 [30] and Multiwfn v.3.8 [31] software packages.

3. Results and Discussion

3.1. Structure Stability and Metal–Support Interactions

Pt-group metals were chosen for their well-documented catalytic performance in formic acid dehydrogenation (FAD). The g-C_3_N_4_ substrate, consisting of 24 carbon and 32 nitrogen atoms, serves as an ideal optimal framework for single-atom anchoring. Previous studies [32] have identified the hollow site as the most thermodynamically favorable position for metal incorporation. Upon adsorption, a stable coordination environment is established between the metal atom and surrounding nitrogen atoms, effectively stabilizing the M@g-C_3_N_4_ SACs.

Figure 1 illustrates the atomic structures of pristine g-C_3_N_4_ and various M@g-C_3_N_4_ SACs. While the undoped substrate retains a planar structure, metal doping induces localized distortions. Notably, Ag, Au, Pd, and Rh occupy nearly central positions within the g-C_3_N_4_ cavity, indicating symmetric coordination with the six surrounding nitrogen atoms. In contrast, Ir, Pt, and Ru exhibit off-center displacements, with Ir and Ru preferentially coordinating with two nitrogen atoms, whereas Pt binds to three nitrogen atoms, reflecting distinct metal–support interactions. The most significant factor contributing to this trend in M@g-C_3_N_4_ structures is the size of the cationic metal atoms [33]. The ionic radius of early-transition metal atoms (such as Ru in the fifth period, Ir and Pt in the sixth period) is larger than that of late-transition metal atoms (e.g., Rh, Pd, and Ag in the fifth period; Au in the sixth period). To accommodate these larger early-transition metal atoms ions, flexible g-C_3_N_4_ undergoes structural distortions.

The stability of SACs is crucial for their catalytic performance. To evaluate their thermodynamic stability, the binding energy (E_b_) of each SAC was computed as follows [34]:

$$[eqn]$$

where E_M@g-C3N4_, E_M_, and E_g-C3N4_ represent the total energies of the SACs system, isolated metal atom, and pristine g-C_3_N_4_ substrate, respectively. A negative binding energy confirms the thermodynamic stability of M@g-C_3_N_4_ SACs, with more negative values, with the magnitude of these energies reflecting the strength of the metal–support interactions, indicating greater structural stability.

As shown in Figure 2, all SACs exhibit negative binding energies, reinforcing their thermodynamic feasibility. The stability trend is as follows: Ag > Ru > Rh > Pt > Ir > Pd > Au. As illustrated in Figure S1, the d orbitals of metal atoms in M@g-C_3_N_4_ SACs exhibit strong electronic interactions with the p orbitals of C and N in the g-C_3_N_4_ framework [35]. These interactions contribute to the remarkable stability of M@g-C_3_N_4_ SACs and serve as a fundamental basis for understanding their catalytic performance in FAD reactions.

3.2. Adsorption Behavior and Pathways for FAD

The adsorption behavior of formic acid and its key intermediates (HCOO, COOH, and CO) on the M@g-C_3_N_4_ SACs surfaces was evaluated to understand the reaction mechanism. The adsorption energy is defined as [36]

$$[eqn]$$

where E_ad_, E_M@g-C3N4+x_, E_M@g-C3N4_, and E_x_ correspond to the adsorption energy, total energy of the SACs with an adsorbate, total energy of pristine SAC, and the energy of gas-phase adsorbates, respectively.

We first computed the adsorption energies (Figure 3a and Table S1) of HCOOH in both cis (HCOOH_cis_) and trans (HCOOH_trans_) configurations, as well as its dehydrogenation and dehydration intermediates. These include the bidentate (HCOO_bi_) and monodentate (HCOO_mo_) formats of formate, which have been identified through in situ infrared spectroscopy [37], along with COOH (cis and trans) [38], and the final CO product. Our calculations reveal that HCOOH adsorption is exothermic across all SACs, with E_ad_ values ranging from −0.01 to −1.01 eV, suggesting its favorable interaction with the catalyst surface. The bidentate formate (HCOO_bi_) intermediate preferentially adsorbs on Pt and Ru SACs, whereas its monodentate counterpart (HCOO_mo_) exhibits weaker binding across all SACs, indicating that monodentate HCOO may be a crucial reaction intermediate in formic acid decomposition. Similarly, COOH adsorption is exothermic (−0.20 to −3.23 eV), while CO strongly adsorbs to all SACs, particularly on Pt (−3.08 eV) and Rh (−3.09 eV), where E_ad_ values exceed −3.00 eV. To further refine our thermodynamic analysis, we incorporated zero-point energy and vibrational contributions to compute the adsorption free energy (ΔG), as shown in Figure 3b and Table S2. Notably, while ΔG for CO adsorption remained negative across all SACs, the inclusion of water molecules led to positive ΔG values on Ir and Pt SACs, indicating potential CO desorption in aqueous environments.

To elucidate the catalytic mechanism, we investigated two competing pathways for formic acid decomposition: dehydrogenation (Path 1, via the HCOO intermediate) and dehydration (Path 2, via the COOH intermediate), as illustrated in Scheme 1.

Path 1 (dehydrogenation): Formic acid initially adsorbs onto the SAC surface in a trans configuration, undergoes O–H bond cleavage to form bidentate HCOO, and then transitions to a monodentate structure before C–H bond cleavage releases H_2_ and CO_2_. Path 2 (dehydration): The reaction proceeds through cis-configured formic acid, followed by C–H bond cleavage yielding trans-configured COOH, which subsequently transitions to a cis configuration before undergoing O–O bond cleavage to generate H_2_O and CO.

3.3. Catalytic Activity and Selectivity

The free energy profiles (Figure 4) highlight significant differences between these pathways on different SACs. Among the various SACs, the maximum elementary step free energy barriers for dehydrogenation were 1.81 eV (Ag), 0.48 eV (Au), 0.49 eV (Ir), 0.51 eV (Pd), 0.28 eV (Pt), 0.45 eV (Rh), and 1.01 eV (Ru), while the free energy barriers for dehydration were 1.18 eV (Ag), 1.90 eV (Au), 1.24 eV (Ir),1.90 eV (Pd), 1.21 eV (Pt), 0.49 eV (Rh), and 1.98 eV (Ru).

As shown in Figure 4, the rate-determining step (RDS) in the dehydrogenation pathway of formic acid involves O–H bond cleavage in HCOOH for Ag, Au, and Ir SACs, whereas for Pd, Pt, Rh, and Ru SACs, the RDS corresponds to a configurational transformation of the adsorbed HCOO intermediate. In the competing dehydration pathway, the RDS is assigned to C–H bond cleavage in formic acid for Ag and Ru SACs, to an adsorption configuration transition of COOH for Au and Pd SACs, and to O–C bond cleavage in COOH for Ir, Pt, and Rh SACs.

The Au and Pd SACs display notably lower energy barriers for dehydrogenation compared to dehydration, with reductions of approximately 1.42 eV for Au and 1.39 eV for Pd. This energetic preference strongly favors the production of H_2_ over the formation of CO, which is a critical advantage considering CO’s detrimental effect on catalyst longevity, making them promising candidates for selective H_2_ generation, as shown in Table S3. In contrast, Rh SACs exhibit similar energy barriers for both pathways, leading to non-selective behavior that may compromise the purity of the hydrogen produced.

To understand the origin of Pd SAC selectivity, we analyzed the density of states (DOS) and d-band center of M@g-C_3_N_4_ SACs (Figure 5). The sharp, well-defined d-orbital peaks suggest near-free-electron behavior [39].

As shown in Figure 6a, the d-band center exhibits a volcano-type relationship with atomic number for period IV metal SACs, whereas it shifts monotonically positive for period V metals. More importantly, as shown in Figure 6b, d-band center correlates linearly with formic acid decomposition selectivity, which indicates the SACs with a d-band center closer to the Fermi level exhibit higher selectivity toward dehydrogenation. This trend suggests that the electronic structure can be used as a predictive descriptor for catalyst performance. The proximity of the d-band center to the Fermi level directly influences the strength of the intermediate adsorption reaction on the catalyst [40]. Consequently, tuning the electronic structure can effectively alter the reactivity and selectivity of the catalytic material.

Further, as shown in Figure 7, the Mulliken charge analysis indicates that charge transfer occurs from the metal atom to the g-C_3_N_4_ support, with values of 0.467 (Ag), −0.024 (Au), 0.164 (Ir), 0.222 (Pd), 0.073 (Pt), 0.244 (Rh), and 0.029 e (Ru). This redistribution of charge modifies the electronic environment at the active sites, influencing both the adsorption strengths of intermediates and the overall reaction kinetics. The modulation of electronic properties through metal–support interactions provides a clear pathway for tailoring catalyst performance, emphasizing the potential of electronic descriptors in guiding the rational design of next-generation SACs.

Overall, our results demonstrate that the careful selection of metal species and support material can effectively steer the reaction pathway toward selective dehydrogenation of formic acid. The combined insights from structural, energetic, and electronic analyses pave the way for designing catalysts that optimize hydrogen production while minimizing undesirable side reactions. Future work will focus on experimental validation and the exploration of bimetallic systems to further enhance catalytic performance.

4. Conclusions

This study elucidates the structure–activity relationships governing Pt-group SACs supported on g-C_3_N_4_ for formic acid dehydrogenation. Pd and Au SACs emerge as highly selective catalysts, achieving dehydrogenation barriers below 0.51 eV while suppressing CO formation through unfavorable dehydration pathways. Rh SACs, despite high activity, suffer from poor selectivity due to nearly equivalent barriers for both reaction routes. Moreover, Pd and Au emerge as promising catalytic candidates for hydrogen generation from formic acid, owing to their relatively lower cost compared to Rh. The thermodynamic stability of SACs is attributed to strong metal–support interactions, evidenced by negative binding energies and d–p orbital hybridization. A linear correlation between the d-band center position and selectivity highlights the critical role of electronic structure in catalytic performance. Charge redistribution (0.029–0.467 e) further fine-tunes the metal centers’ reactivity. These results underscore the potential of g-C_3_N_4_-supported SACs for efficient hydrogen production and provide a roadmap for optimizing catalyst design through electronic modulation and support engineering. Future work should focus on experimental validation and exploring synergistic effects of bimetallic SACs.