Lujain
Alrais
a,
Indranil
Dutta
a,
Amol
Hengne
b,
Priyanka
Chakraborty
a,
Edy
Abou-Hamad
c,
Shibo
Xi
b,
Mohammad Misbahur
Rahman
a,
Jia
Zhang
d,
Benjamin W. J.
Chen
*d,
Jean-Marie
Basset
*ae and
Kuo-Wei
Huang
*ab
aCenter for Renewable Energy and Storage Technologies (CREST) and Division of Physical Science and Engineering, King Abdullah University of Science and Technology, Thuwal 23955-6900, Saudi Arabia. E-mail: jeanmarie.basset@kaust.edu.sa; hkw@kaust.edu.sa
bInstitute of Sustainability for Chemicals, Energy and Environment (ISCE2), Agency for Science, Technology and Research (A*STAR), Singapore 138634, Singapore
cCore Labs, King Abdullah University of Science and Technology, Thuwal 23955-6900, Saudi Arabia
dInstitute of High-Performance Computing (IHPC), Agency for Science, Technology, and Research (A*STAR), Fusionopolis Way, #16–16 Connexis, Singapore 138632, Singapore. E-mail: benjamin_chen@ihpc.a-star.edu.sg
eEcole National Supérieure de Chimie de Paris (ENSCP), 75231 Cedex 05 Paris, France
First published on 14th February 2025
The dehydrogenation of formic acid can provide an efficient pathway for hydrogen generation in the presence of a suitable catalyst. Homogeneous catalysts have been extensively studied and utilized for highly active and selective processes compared to conventional heterogeneous catalysis, which often shows lower reactivity and selectivity. However, the latter is preferred for practical applications, considering its easy separation and recyclability. By incorporating a homogeneous organometallic complex on an appropriate support, the unique features of both catalysts can be combined and utilized effectively. Herein, we investigate the immobilization of an iridium picolinamide complex (1) supported on 3D fibrous modified silica that demonstrates high accessibility. The support involves a tetracoordinate aluminum hydride site featuring a strong Lewis acidic nature. A study of the interaction and coordination sites around the surface fragment was conducted via various techniques, including elemental analysis, FT-IR, solid-state NMR, XAS, and first-principles calculations, which provided informative data. We explored the use of solid additives in a solvent-free reaction medium and avoided utilizing volatile bases to achieve process feasibility with a high TOF of 40000 h−1.
Surface organometallic chemistry (SOMC) can provide a practical strategy for surface preparation and processing.36–39 SOMC is concerned with the investigation into the structure, reactivity, and properties of the surface fragments. This approach allows us to understand the reactions of organometallic complexes with the surface of a support to prepare single-site catalysts on well-defined surfaces with several advantages such as improved stability and resistance to aggregation.
In this regard, the utilization of tetracoordinate aluminum hydride [(Si–O–Si
)(
Si–O–)2Al–H] [Al-H]40 as an anchor site for organometallic complexes has shown interesting results, and its Lewis acidic nature can alter complex performance.41,42
Nowadays, researchers are advancing toward greener and more cost-effective methods for FA dehydrogenation by focusing on additive-free and solvent-free procedures43,44 because using solvents often leads to waste generation and can pose environmental hazards. Conducting the reaction under neat conditions and eliminating the solvents from the reaction medium would lead to a straightforward process with fewer purification steps and further reduce the operational cost. However, ensuring effective catalyst-to-substrate interaction without solvents, especially when using heterogeneous catalysts, can be challenging.
More importantly, reports on the neat FA dehydrogenation by immobilized catalysis are scarce. Our previous investigation has provided insight into the positive influence of the [Al-H] surface species, in which the immobilized iridium PN3P-pincer complex has shown superior activity and stability compared to its homogeneous analog.45 In the quest for an alternative catalytic system, we evaluated the amine linkage (NH) effect by employing a picolinamide iridium complex 1 supported on [Al-H]. We investigated the impact of changes at the coordination sites around the surface fragment and the effect on reactivity. The FT-IR, solid-state NMR, XAS, and first-principles density functional theory (DFT) calculations provide informative data to characterize the obtained immobilized complex 2. The immobilized catalyst has been tested under solvent-free/neat conditions, showing unprecedented reactivity in FA dehydrogenation with a high TOF of 40000 h−1.
Immobilized catalyst 2 was synthesized (details in the ESI†) and the obtained pale-yellow material was stored in a glovebox for further characterization (Scheme 1). Throughout the gas analysis, it was found that H2 was released during the immobilization reaction (Fig. S7†). Elemental analysis of 2 suggested an iridium loading of 4.1 wt% corresponding to 0.2 mmol g−1. The carbon and nitrogen values of C/Ir and N/Ir were found to be 22.4 and 3.9, respectively, in close agreement with the theoretical values of 20 and 4 (Table 1).
Element | wt% | mmol g−1 | M/Ir |
---|---|---|---|
a Carbon value after subtracting the residual support carbon. | |||
Ir | 4.1 | 0.2 | |
Al | 3.7 | 1.3 | 6.5 |
C | 5.3 | 4.5 | 22.5 |
N | 1.1 | 0.78 | 3.9 |
The FT-IR spectra of the immobilized catalyst 2 and the support [Al-H] were recorded (Fig. 2). The FT-IR spectrum of the support (a) showed a small peak at 3743 cm−1, which was attributed to the remaining silanol groups Si–OH. Signals of Si–H and Si–H2 were observed at 2252–2170 cm−1, which remained intact after the immobilization. The Al-hydride signal was detected at 1943 cm−1, confirming the successful preparation of the support. The spectra of complex 2 (b) indicated a reduction in the Al-H signal, and the emergence of all characteristic peaks corresponded to the successful immobilization of the organometallic complex on the support. Furthermore, the set of bands in the range of 3060–2800 cm−1 and 1490–1360 cm−1 were ascribed to the stretching and bending vibrations of the C–H bond of the alkyl and aryl groups. The broad peak at 3430 cm−1 was attributed to the vibration of the N–H fragment.
![]() | ||
Fig. 2 The FT-IR spectra of (a) the [Al-H] support. (b) The immobilized catalyst 2. (c) Subtraction of b − a. |
The peak at 1690 cm−1 was assigned to the vibration of the CO bond, whereas the signals between 1600 and 1520 cm−1 resulted from the vibration of the C
C and C
N bonds within the homogeneous complex.
Further characterization was conducted by implementing multinuclear solid state (SS) NMR. The 1H-MAS SS-NMR spectrum of 2 exhibited significant broad signals at δ 1.0 and δ 1.7 ppm arising from the protons of the alkyl fragments δH(CH3) and δH(CH2) within the structure of the complex. The peak at 4.4 ppm was assigned to δH(N–H, Al-H, and Si–Hx). The aromatic protons’ δH(HCC) signal appeared at δ 7.82 ppm (Fig. 3a). The 13C CP/MAS solid-state NMR spectra (Fig. 3b) showed the characteristic peaks of the immobilized catalyst at δ 7.0 and δ 22.8 ppm, which correspond to various alkyl groups (C–H3 and C–H2) bonded to the ligand. These peaks correlate with the proton signal at δ 1.7 ppm observed in the 1H–13C HETCOR solid-state NMR spectra (Fig. S8†). The signal at δ 38.2 ppm was assigned to the methyl groups linked to the amine ((CH3)2–N) fragment. The peaks at δ 85.5 and δ 111.0 ppm were attributed to the aromatic carbons (C
C) and (N–C
C) within the pyridine and methylcyclopentadienyl groups. Signals overlapping at δ 146.5 ppm were assigned to (C
N) and (C–N) and the peak at δ 154.4 ppm was attributed to the (N–C
O) moiety. These results were consistent with the homogeneous complex 1, confirming the successful immobilization of 1 on the [Al-H] support.
Additionally, the 1H–1H DQ solid-state NMR spectra (Fig. S9†) showed an autocorrelation on the diagonal of proton resonances from the different alkyl groups and aromatic protons. In other words, at least two protons of these fragments are in close proximity to each other. Outside the diagonal, a strong correlation at δ 5.4 ppm was attributed to the protons near the CH2–N–H group, [5.4 ppm in F1: δH(NH) + δH(CH) = 4.4 + 1.0]. The 1H–H1 TQ solid-state NMR spectra of 2 (Fig. S9†) were also recorded; an autocorrelation at 5 ppm in F1 corresponded to the methyl groups. These observations offer valuable insights, and to the best of our understanding, the proposed structure of 2 in Scheme 1 is the most suitable option based on available data.
To further investigate the coordination environment between complex 1 and the support, 27Al solid-state NMR was conducted to study the coordination state of the Al site on the surface. Fig. 4 displays a broad signal at δ 40.16 ppm, corresponding to a distorted tetra-coordination arrangement. A small shoulder at δ 6.87 ppm attributed to the penta-coordinated Al site was also observed, suggesting an acid–base interaction between the surface species and complex 1.
A notable broadband in the region of δ 150–200 ppm was associated with Al-Alkyl, remaining on the support surface, arising from the incomplete transformation of Al-isobutyl to Al-hydride during the support preparation procedure.40,46,47 The high-angle annular dark field imaging (HAADF) TEM analysis (Fig. 5) confirmed the preservation of the fibrous 3D-mesophase of the support and demonstrated a well and uniform distribution of the elements on the support surface, as evidenced by EDS elemental mapping. Additionally, the analysis revealed the presence of iridium at 3.27 wt%, which closely matches the results obtained through elemental analysis (Fig. S10†). To elucidate the electronic structure and coordination environment at the atomic level, X-ray absorption spectroscopy measurement was conducted. The X-ray absorption near-edge structure (XANES) spectra of 2 and the reference compound at the Ir L-edge were recorded (Fig. S11†). IrO2 was used as a reference for the +4 oxidation state, while IrCl3 served as a reference for the +3 oxidation state. The EXANS energies for 2 and the reference IrCl3 were found to be identical and equal to 11218.17 eV, suggesting an average oxidation state of +3 for the supported system. Subsequently, analysis of the extended X-ray fine structure (EXAFS) Fourier transforms at the iridium L-edge for 2 (Fig. S12†) showed spectral features corresponding to those of the homogeneous Ir complex, thus suggesting the successful formation of the desired Ir complex on the support surface.
![]() | ||
Fig. 5 HAADF-STEM image of 2: (a) images and elemental mapping of aluminum (b), iridium (c), nitrogen (d), silicon (e), and superimposed elements (f). |
Entry | Catalyst | FA![]() ![]() |
TOF | Yield |
---|---|---|---|---|
a General conditions: FA (1.2 ml, 30 mmol), catalyst 2 (10 mg, 2 μmol), temperature (90 °C). b FA (1.2 ml, 30 mmol), [Al-H] (20 mg). c FA (0.6 ml, 15 mmol), catalyst 2 (10 mg, 2 μmol), temperature (90 °C). d Catalyst 1 (1.2 mg, 2 μmol), FA (1.2 ml, 30 mmol), temperature (90 °C). The maximum TOF was estimated based on the H2 production in the first 10 min. | ||||
1b | [Al-H] | 1![]() ![]() |
No reaction | 2% |
2 | 2 | 1![]() ![]() |
No reaction | 7% |
3 | 2 | 2![]() ![]() |
37![]() |
90% |
4 | 2 | 1![]() ![]() |
40![]() |
88% |
5c | 2 | 1![]() ![]() |
17![]() |
81% |
6d | 1 | 1![]() ![]() |
45![]() |
93% |
The optimized catalytic performance test of catalyst 2 was conducted under solvent-free/neat conditions, and a TON of 25000 was recorded. Despite our several attempts to optimize different parameters to increase the TON, it did not further improve and indicated limited stability of the immobilized catalyst under these harsh conditions. We examined the batch-wise process; however, catalyst 2 decomposes after only a few cycles of dehydrogenation upon accumulation of FA and an increase in the acidity of the medium. Likewise, an optimized catalytic performance test was conducted using the homogeneous complex. Catalyst 1 demonstrated limited stability and a TON of 16
000 was achieved.
To gain a deeper understanding, a control experiment was carried out at the end of the stability test. The reaction mixture was filtered, and the used solid catalyst was subjected to elemental analysis which revealed a significant decrease in iridium loading from 0.2 to 0.005 mmol g−1. On the other hand, the aluminum loading was found to be 1.2 mmol g−1, a negligible decrease compared to the fresh catalyst. The liquid phase, which contained the leached catalyst, was dried first to remove the residual FA and water. Afterward, it was dissolved in d-chloroform and analyzed using liquid-state NMR analysis. The spectrum of the leached catalyst displays a significant change compared to the homogeneous analog 1, and the ligand appears rather fragmented. As far as we can tell, this is potentially caused by the hydrolysis of the ligand backbone, leading to the rupture of the catalyst core (Fig. S17 and S18†).
![]() | ||
Fig. 7 First-principles modeling of FA decomposition over catalysts 1 and 2. (a) Atomistic models of the activated catalysts. The H* state is illustrated for simplicity. Ligands of the catalysts are depicted as sticks and O and H atoms of the adsorbates and spectator H2O molecules are shown in green and light blue, respectively, for clarity. (b) Potential energy diagrams for FA decomposition. Energies correspond to Gibbs free energies relative to the initial states. Insets show the illustrations of catalytic intermediates, with the adsorbed species highlighted in blue. Atomistic models of the intermediates are provided in Fig. S20 and S21† for catalysts 1 and 2, respectively. |
We model the FA dehydrogenation reactivity of catalysts 1 and 2 in their activated state, where Cl− is exchanged with HCOO−. This is consistent with earlier studies showing that similar catalysts undergo an activation process involving this ligand exchange.50,51 Once activated, the catalysts can undergo multiple turnovers with FA in the catalytic cycle (Fig. S19†). As water is always present in commercially available neat FA, we additionally considered 2 explicit spectator H2O molecules in the system. Our previous works showed that water can play key roles by significantly lowering the barriers of several elementary steps in the reaction via hydrogen bonding and proton shuffling.41,45 It should be noted that as CO is not produced in the reaction, we have only considered the formate (HCOO*) pathway and not the carboxyl (COOH*) pathway, as the decomposition of COOH* is known to easily generate CO*.52
From our calculated potential energy diagrams (PEDs), the reaction mechanism for FA dehydrogenation has been found to be similar for both catalysts (Fig. 7b). After activation, HCOO* is adsorbed on the catalysts. From this initial state, HCOO* then flips to form HCOOm*, where the C–H bond points downwards facing the Ir center. Deprotonation of HCOOm* releases CO2 and leaves behind a hydridic H* on the Ir center. Next, an HCOOH molecule binds to the catalyst via hydrogen bonding of its CO group with the –NH2+ group on the catalyst ligand. The acidic –OH group of HCOOH then reacts with the hydridic H* on Ir via proton shuffling mediated by a water molecule. This acid–base reaction releases H2 in a spontaneous process for 2 and via an intermediate 2H* state for 1. Finally, the initial state is regenerated as the HCOO− formed binds to the Ir vacancy created by the loss of H2.
The energetic spans of the computed reaction mechanisms over catalysts 1 and 2 are 0.90 and 0.94 eV, respectively (Fig. 7b). This is consistent with the slightly lower TOF over the immobilized catalyst 2. The rate-determining intermediate (RDI) is the H* species, which is likely bound with an HCOO− counterion. This agrees with the NMR analysis of dried catalyst 2 after the reaction with FA. The higher energetic span for 2 is attributed to the slightly higher barrier for H2 production from the H* species. However, the PEDs are in general similar, which is consistent with our finding that the calculated Bader charges of Ir on 1 and 2 are the same (+0.69 and +0.67|e−|, respectively). Overall, our simulations reveal that the reactivity of catalyst 1 is not significantly affected by grafting to the Al center as the atoms attached to Al are relatively far away from the active Ir center, allowing 2 to maintain a high TOF similar to 1.
DFT calculations revealed similar mechanisms of FA decomposition over catalysts 1 and 2. The results showed that the high TOF of 2 can likely be attributed to the grafting site being located away from the active Ir center, which prevents a decrease in the reactivity of 1. However, catalyst 2 demonstrated limited stability. As outlined in the proposed mechanism, which is supported by DFT calculations, the reaction pathway involved the protonation of the N–H side arm, consequently weakening the bond between complex 1 and the support and facilitating leaching from the surface.
Moreover, analysis of the deactivated catalyst after the optimized catalytic performance test indicated a significant reduction in the iridium loading. We hypothesized that this is potentially caused by the hydrolysis of the ligand backbone, leading to the rupture of the catalyst core. As the grafting anchoring point is positioned away from the Ir active site, it results in minimal changes to the coordination sites. The data presented here exemplify the importance of understanding the structure–reactivity relationship. Consequently, significant advancement can be made in the design of catalysts to tailor catalysts for specific applications.
Footnote |
† Electronic supplementary information (ESI) available. CCDC 2412317. For ESI and crystallographic data in CIF or other electronic format see DOI: https://doi.org/10.1039/d4dt03521f |
This journal is © The Royal Society of Chemistry 2025 |