Open Access Article
Sebastian Seidel
a,
Keanu V. A. Birkelbach
a,
J. Mädickea,
M. Pilaski
b,
Peter J. C. Hausoul
a and
Regina Palkovits
*a
aInstitut für Technische und Makromolekulare Chemie, RWTH Aachen University, Worringerweg 2, 52074 Aachen, Germany. E-mail: Palkovits@itmc.rwth-aachen.de
bZBT - Zentrum für Brennstoffzellentechnik GmbH, Carl-Benz-Straße 201, 47057 Duisburg, Germany
First published on 19th February 2026
A synthetic methodology for the immobilization of solid molecular catalysts (SMCs) in cylindrically shaped porous glass substrates was developed to facilitate their use in a continuous plug-flow reactor. The glass supports were coated with tailor-made alkoxysilanes to enable covalent attachment of, ‘in pore’ polymerized SMCs. Inorganic–organic hybrid catalysts based on BINAP-, terpyridine- and triazine were successfully prepared and tested in formic acid dehydrogenation after impregnation with 0.1 wt% Ru or Ir. In a plug-flow reactor, TOFs of up to 11
500 h−1 and selectivities of 75 ppm CO could be observed at 160 °C, while at 100 °C the activity remained at 6200 h−1 with a CO concentration of 79 ppm.
A strategy unifying these concepts is the utilization of solid molecular catalysts (SMCs) for the synthesis of formic acid from CO2 and H2 as well as formic acid dehydrogenation (FAD) to release H2 on demand. SMCs combine the high activity of homogeneous catalysts and the straightforward handling and recycling of heterogeneous catalysts and thus offer unique advantages for technical applications.14,15 To date, several different immobilization strategies, including covalent attachment, ionic interactions and adsorption, have been developed (see Table S2 for a more detailed comparison).14–22 For FAD, immobilized catalysts, based on surface-grafted silica (Fig. 1) were reported as early as 2009. Gan et al. reported a Ru(TPPTS)L4@PPh2(CH2)n/SiO2 catalyst with turnover frequencies (TOFs) of up to 2780 h−1 at 110 °C in 13 wt% aqueous FA. Another notable silica immobilized catalyst is the IrH3@PN3P/SiO2 pincer catalyst reported by Alrais et al. While utilizing a mixture of FA and caesium formate at 90 °C, a TOF of up to 13
290 h−1 and a turnover number (TON) of up to 540
000 were obtained.12 In addition, several noble-metal-free immobilized catalysts have been reported in recent years, decreasing the economic burden of the catalytically active material at the cost of more sophisticated ligand synthesis. One prominent example is the Fe(BF4)2@polyRPhphos/SiO2 catalyst reported by Stathi et al.13 It exhibits TOFs of up to 7600 h−1 and TONs of up to 8500 in the decomposition of 40 wt% FA in propylene carbonate at 90 °C. The technical application of catalysts supported on SiO2 is usually limited by the stability of the support towards solvents, acidic or basic conditions and subsequent leaching of the ligand and metal species in continuous operation.23
![]() | ||
Fig. 1 Examples of molecular catalysts immobilized on SiO2 (left: Ru(TPPTS)@PPh2(CH2)n/SiO2,11 middle: IrH3@PN3P/SiO2,12 and right: Fe(BF4)2@polyRPhphos/SiO2 13). | ||
Several highly active SMCs were reported based on fully organic covalent triazine frameworks (CTFs), covalent organic frameworks as well as polyaromatic frameworks, circumventing the degradation of the catalysts under acidic conditions (Fig. 2).24–32 Organic macromolecular systems often exhibit a pronounced swelling behaviour. Furthermore, shaping of such catalyst materials for industrial application remains largely unstudied.33 Depending on the material, several strategies for shaping are viable. Materials can be shaped into monoliths during or after synthesis by using pressure or a binder and subsequently cutting or grinding them to the required shape.34
![]() | ||
| Fig. 2 Ru(FA)2@pBINAP (left), IrCl3@PyrTerpy and Ir(Cp*)Cl@CTF by Gascon et al.24 | ||
Another option is coating the surface of a rigid support such as silica, glass or zeolite with the active catalyst phase, to yield inorganic–organic hybrid catalysts.28 The latter has previously been investigated by Bavykina et al. for CTF systems supported on cordierite monoliths for FAD of 3 M aqueous FA with TOFs of up to 207
000 h−1 at 80 °C in a continuously stirred tank reactor (CSTR).28–35 Several further continuous reaction systems for FAD have previously been reported in literature. Ruiz-López et al. describe a CSTR for liquid-phase FAD with bimetallic PdRu particles (5 wt% RuPd/C3N4) supported on graphitic carbon nitride, exhibiting a TOF of 30 h−1 and a TON of 60 for the decomposition of 1 M aqueous FA at 60 °C.36 Sponholz et al. developed a homogeneous reactor system based on RuCl2-bis(1,2-bis(diphenylphosphino)ethane), which, in a CSTR with connected gas purification and the addition of N,N-dimethyloctylamine as an activating base, exhibited a TOF of 1000 h−1 over a reaction time of 45 days (TON = 1
000
000) at 25 °C.37 By increasing the amount of FA added and the temperature to 60 °C, this system achieved a TOF of 16
000 h−1 at 6.5 bar. Sawahara et al. previously reported a Cp*Ir@PEI flow reaction system capable of FAD of 20 M aqueous FA at 80 °C using a syringe pump with a TON of 332
000 over 2100 h and CO concentrations of below 0.1 ppm.38
Here we present the preparation of shaped hybrid inorganic–organic catalysts via surface modification and in-pore Suzuki coupling of catalytically active phase in porous borosilicate cylinders. Three different catalysts based on Ru/BINAP, Ir/terpyridine and Ir/CTF are prepared and compared in batch formic acid dehydrogenation reactions. To demonstrate the applicability of the approach as well as long term stability of the catalyst, the best performing system is investigated using a continuous plug flow reactor.
![]() | ||
| Fig. 3 Coating procedure of SiO2 with alkoxysilanes (PhSi(OMe)3 and 4-BrPhSi(OEt)3) and subsequent Suzuki coupling of surface grafted BrPh@SiO2 with 4-fluoro phenyl pinacol. | ||
Silica support (Saint Gobain) was ground to a powder and subsequently coated with the silanes trimethylsilyl chloride (TMSCl) (TMS@SIO2), trimethoxyphenylsilane (PhSi(OMe)3) or triethoxy-(4-bromophenyl)silane (BrPhSi(OEt)3).42,43 The specific surface area (SSA) of the utilized SiO2 remained comparable with 228 m2 g−1 for the unmodified material and 223 m2 g−1 for BrPh@SiO2. After the reaction, the DRIFTS spectrum (Fig. 4, left) of TMS@SiO2 exhibits the expected signals at 2970 cm−1, corresponding to stretching vibration modes of C–H bonds in surface –CH3 groups formed during the silylation process.44,45
![]() | ||
| Fig. 4 DRIFTS (left, inset: zoom on CH-region of the three modified materials), solid-state 29SI (middle) and 13C (right) NMR spectra of pristine SiO2, TMS@SiO2, Ph@SiO2 and BrPh@SiO2. | ||
Coating with PhSi(OMe)3 leads to additional signals corresponding to Carom.–Carom. at 1433 cm−1 and 1595 cm−1, as well as Carom.–H vibration bands at 3060 cm−1 and 3077 cm−1 thus validating the coating procedure.46 Usage of the target silane BrPhSi(OEt)3 leads to a vibration band at 1580 cm−1, attributed to a C–C stretching vibration mode (Fig. S1 and S2).47 Additionally, a band emerges at 1486 cm−1, also attributed to BrPhSi(OEt)3. Subsequently, solid state magic angle spinning (MAS) 29Si NMR (Fig. 4, middle) was employed to further elucidate the structure of BrPh@SiO2. In addition to the Qn-type Si(OSiO3)n(R)4−n signals (Q2 = −92 ppm, Q3 = −100 ppm and Q4 = −110 ppm) exhibited by the unmodified silica material, a Tn-type silicon signal at T3 = −67 ppm Si(OSiO3)n(R)3−n corresponding to the BrPhSi(OEt)3 coating could be observed.41,48 13C NMR reveals the signals expected for BrPhSi(OEt)3 at 125 ppm, 130 ppm and 134 ppm (Fig. 4, right). Additionally, signals can be observed at 16 ppm and 57 ppm resulting from residual ethoxy groups, suggesting incomplete coupling of BrPhSi(OEt)3 to the surface. XRD reveals no significant change due to the coating (Fig. S3).
A subsequent Suzuki-cross coupling reaction was explored as a pathway to post-functionalize the surface-bound fragments. Due to the lability of the model support SiO2 towards the basic conditions during the reaction, protection of the surface was achieved by treating the BrPhSi(OEt)3 modified support (BrPh@SiO2) with TMSCl and washing with EtOH (TMS + BrPh@SiO2). This led to an additional signal in the 29Si-SS-NMR (Fig. 4, middle) at 15 ppm with a concurrent decrease in the Q2 signal at −92 ppm and a relative increase in the Q4 signal at −110 ppm. This is in line with the behaviour expected for a subsequent protection of exposed Si atoms with TMSCl.49 13C NMR (Fig. 4, bottom) reveals a relative increase in the residual ethoxy signals at 16 ppm and 57 ppm, attributed to a acidic ethoxylation of free surface silanol groups. Subsequently, the silane-functionalized material was coupled with 2-(4-fluorophenyl) pinacol boron ester at 85 °C in DMF. The structure of the resulting material was investigated using DRIFTS spectroscopy, 13C, 19F as well as 29Si NMR. In 29Si NMR a decrease of the TMS signal can be observed which is attributed to base induced hydrolytic cleavage of the support. This is underlined by the decrease in the signal of the base labile surface bound ethoxy groups observed in 13C NMR. Nonetheless, the signal attributed to a T3-type aromatic silicon atom can be observed at −67 ppm. 13C NMR reveals additional signals at 165 ppm, 140 ppm and 115 ppm attributed to the 4-fluorobiphenyl moieties. In 19F NMR, a broad signal at −120 ppm can be attributed to the surface bound 4-fluorobiphenyl, whereas a sharp peak at −115 ppm represents unreacted substrate adsorbed to the surface (Fig. S5). This confirms that Suzuki cross-coupling can be successfully applied to a surface functionalized silica, albeit with some loss of the functional group density.
:
36% H2O2 7
:
3) and Piranha base (36% H2O2
:
30% NH4OH
:
H2O 1
:
1
:
5) using the procedure described for the silica materials.39
DRIFTS spectra of the coated materials show no significant change in comparison to the spectra of unmodified material (Fig. S6). Analysis of the resulting glass using 29Si NMR was performed and revealed one broad signal at −100 ppm attributed to the bulk material (Fig. S7).50,51
Additionally, a weak shoulder between −70 ppm and −85 ppm can be observed, which is attributed to a surface bound T-type organosilane. No exact structure determination was possible due to the overall low fraction of reactive surface Si atoms in comparison to the inaccessible silicon in the bulk material due to the low surface area of the shaped glass supports.
In addition, commercially available microscope slides were coated accordingly to allow characterization via ATR-IR (Fig. 5). The strong band near 910 cm−1 is attributed to the asymmetric stretching vibration of BØ4-tetrahedra (Ø denotes a bridging oxygen atom).52–54 In the region below 900 cm−1, deformation modes are found and attributed to the various constituents of the network such as B–O–B and Si–O–Si bridges or mixed Si–O–B. The spectrum for BrPh@glass differs from pristine glass by very weak additional v(C–C) bands at 1579 cm−1 as well as 1483 cm−1 which are attributed to the coating with BrPhSi(OEt)3 and were previously observed for BrPh@SiO2.
![]() | ||
| Fig. 5 ATR-IR spectrum of cleaned microscope slides before coating (top) and glass surface-coated with BrPhSi(OEt)3. Insert: Magnified C–C vibrational area of ATR-IR spectrum. | ||
To further verify the successful silane coating on glassy supports, the contact angle with H2O was measured for the different silane coatings as well as uncoated glass (Fig. S8). Unmodified glass exhibits a low contact angle of 29.8°, in line with a hydrophilic surface featuring Si–OH groups.40 After coating with TMSCl, an angle of 99.8° was observed, demonstrating a decreased hydrophilicity. After coating the surface with PhSi(OMe)3 a contact angle of 89.8° was observed, while the BrPhSi(OEt)3 coating active in Suzuki-coupling reactions leads to a contact angle of 100.0°. Thus, it could be shown that the surface coating method proceeds successfully for all the studied glasses.
Three supported porous polymeric catalyst systems were obtained by polymerization of three different monomers and subsequent metal loading (Fig. 6): 1.) 5,5′-dibromo-2,2′-bis(diphenylphosphino)-1,1′-binaphthyl (BINAPO), yielding a ruthenium/polyphosphine catalyst (Ru@pBINAP/G),56 2.) 6,6″-dibromo-4′-(1H-pyrrol-2-yl)-2,2′-6′,2″-terpyridine, yielding a iridium/terpyridine-based solid molecular catalyst (Ir@PyrTerpy/G),56 3.) 2,4,6-tris(4-(pinacol boron ester)phenyl)-1,3,5-triazine yielding a CTF-type supported iridium catalyst (Ir@CTF/G).57
![]() | ||
| Fig. 6 Cross coupling-based synthesis methods for BINAP-, terpyridine- and CTF-based hybrid catalyst materials. | ||
Due to insufficient solubility of the monomers in the solvent DMF, the porous support was pre-loaded with the respective polymerisation solution in DMF by incipient wetness impregnation. Thereby, a sufficient amount of [Pd(PPh3)4] and Cs2CO3 in addition to the monomers could be introduced into the cylindrical or pelletized porous borosilicate supports under argon flow. After 72 h at 85 °C the resulting crude hybrid supports were obtained. Following an oxidative cleaning procedure utilizing H2O2/HCl to remove Pd(0) leftover in the organic fraction, further processing was performed depending on the organic coating. In the case of Ru@pBINAP/G, the BINAP oxide containing material, this involved the reduction of the phosphine oxides to phosphines by transfer (de)oxygenation using PPh3/HSiCl3.58 After washing with toluene to remove leftover PPh3 and subsequent impregnation with Ru(methyl-allyl)2COD (0.1 wt%), Ru@pBINAP/G could be obtained as a yellow-beige shaped hybrid catalyst. The metal uptake for all materials was investigated using ICP of the metal solution after filtration a loading of 0.065 wt% was found for Ru@pBINAP/G (Table 1). Analogously, Ir@PyrTerpy/G was obtained after the polymerisation and oxidative treatment by wet impregnation with IrCl3·H2O in MeOH (0.087 wt% Ir uptake) without prior reduction as a yellow to orange solid. Ir@CTF/G was obtained after polymerisation by wet impregnation with Ir(COD)(acac) in MeOH followed by reduction in H2 at 400 °C, obtaining the hybrid catalyst (0.04 wt% Ir uptake).29
| Entry | SSA [m2 g−1] | Measured metal loading [wt%] |
|---|---|---|
| Glass | 1.0 | — |
| Ru@pBINAP/G | 4.9 | 0.065 ± 0.002 |
| Ir@PyrTerpy/G | 3.3 | 0.087 ± 0.006 |
| Ir@CTF/G | 4.2 | 0.040 ± 0.008 |
Non-destructive analysis of the hybrid catalysts as-synthesized using IR (Fig. S13) or XRD (Fig. S14) was inconclusive due to the small influence of the surface modification on the bulk spectra. For this reason, physicochemical and spectroscopic characterizations were carried out after ball-milling the catalysts into a powder. Nitrogen physisorption reveals a slight increase in specific surface area from 1.0 m2 g−1 of the unmodified glass to 4.9 m2 g−1 for Ru@pBINAP/G, 3.3 m2 g−1 for Ir@PyrTerpy/G and 4.2 m2 g−1 for Ir@CTF/G. This change exemplifies the formation of additional, high surface area structures of the target polymer networks attached to the glass shape precursor. Thermogravimetric analysis (TGA) in N2 reveals a mass loss of less than 1% up to 850 °C (Fig. S15).
As previous studies on comparable polymers revealed thermal stabilities up to 400 °C, this observation is best explained by the small mass fraction of the polymer compared to the even more temperature stable glass precursor.55,56,59 At the same time, this underlines how little polymer is required to both significantly increase the SSA of the catalyst material and to successfully uptake metal. The morphology of the resulting hybrid materials was investigated using scanning electron microscopy/energy dispersive X-ray spectroscopy (SEM/EDX) (Fig. 7 and 8 and S10–S12). The SEM images of the hybrid catalysts reveal a distribution of the organic material over the surface of the glassy support. The polymers primarily fill the free pores of the inorganic carrier and form an interpenetrating network. The close interaction of the polymers with glass the surface also suggests, as was verified using the model studies above, that the polymers are also chemically bound to the surface. Differences in the polymer morphology between the three catalysts can be observed in SEM. While Ru@pBINAP/G and Ir@CTF/G present as net-like structures that cover the entire glass surface, the terpyridine-based system Ir@PyrTerpy/G forms a rather rigid, block-like structure with large cracks in the surface. EDX analysis of Ru@pBINAP/G reveals the presence of P as well as Ru atoms inside the organic domain with a concurrent absence in the glass-domain, in line with metal loaded BINAP moieties. X-ray diffractograms of the hybrid catalysts are dominated by the spectrum of the amorphous glass support (Fig. S12). Analysis of Ru@pBINAP/G using 31P NMR spectroscopy after ball-milling reveals signals at −16 ppm, assigned to reduced phosphines as well as at 24 ppm, assigned to phosphine oxides (Fig. 9).
![]() | ||
| Fig. 7 Scanning electron microscopy images of the hybrid catalysts Ru@pBINAP/G (left), IrCl3@PyrTerpy/G (middle), Ir(acac)@CTF/G (right). | ||
200 h−1 (63% of SMC) and a CO-concentration of 330 ppm. A study on the recycling of Ru@pBINAP/G revealed significantly improved handling of the catalyst material, albeit a decrease in catalytic activity was observed after each cycle (Fig. 10, bottom). The TOF decreased to 19
100 h−1 with a concurrent decrease of CO-concentration to 450 ppm over four catalyst cycles. The lower relative decrease in TOF is attributed to decreased loss of active catalyst material due to the simplified recycling of the hybrid catalysts. The increase in CO-concentration is likely caused by oxidation of the phosphines. In contrast, the terpyridine-based material exhibited a more pronounced decrease in activity from the free SMC of 40
250 h−1 and a CO-concentration of 208 ppm down to 11
000 h−1 with an improved CO-concentration of 99 ppm for Ir@PyrTerpy/G at 160 °C. The CTF-based material Ir@CTF/G exhibited the lowest TOF of 8700 h−1 and highest CO-concentration of 1830 ppm as SMC, and accordingly also as hybrid catalyst with a TOF of 6300 h−1 and a CO-concentration of 3200 ppm.
τ = contact time (min)
X(FA)HPLC = conversion of FA as determined by HPLC analysis (—)
ρ(FAaq.) = density of aqueous formic acid solution (g mL−1)
(FA) = volumetric flow of FA set using HPLC pump (ml min−1)
S(Dehy.) = selectivity towards FA dehydrogenation.
As selectivity in all cases was determined to be above 99.9%,
S(Dehy.) is approximated by 1.
The pressure of the reactor was set to 5 bar during FAD reactions, while the flowrates of the aqueous FA and the reaction temperature were screened. The CO-concentration was determined by offline gas chromatography and the conversion of the reaction determined from the liquid phase by offline HPLC and metal leaching by ICP-MS. Ru@pBINAP/G and Ir@PyrTerpy/G were investigated in continuous operation with 10.0 wt% FA. To this end, after an induction phase at 160 °C, the temperature and contact time of the system was screened. Therefore, the flow rates of aqueous FA were varied between 1.0 and 2.5 mL min−1 at 100 °C, 130 °C and 160 °C respectively and the conversion and CO-concentration were determined (Fig. 12). Initially, Ru@pBINAP/G was investigated. To investigate if deactivation occurred during the reaction, the initial conversion at 160 °C was compared to the conversion at the end of the catalyst run. No significant deactivation could be observed until the TON of 11
000 was reached after 300 min. After the catalyst conditioning, the temperature was lowered to 100 °C and the contact time set to 1.4 min (1.0 mL min−1 FA). Even at relatively low temperatures, a H2 productivity (
(H2)) of 44 mL min−1 with a CO-concentration of 380 ppm was achieved. With increasing flow rate and thus decreasing contact time at a constant temperature, the conversion of the reaction decreased, while CO-concentration remained constant throughout the experiment. ICP-MS reveals a Ru metal loss less than 0.33 ± 0.19 mol%. This extremely low leaching, contrasted by an observable catalyst deactivation at higher TON, highlights the fact that even the products of the Ru deactivation pathway preferably remain on the catalyst surface, resulting in exceptional metal retention, which is one of the key drivers of cost in industrial application.
![]() | ||
| Fig. 12 Exemplary reaction scheme for continuous FAD in the PFR with the hybrid catalyst Ru@pBINAP/G. | ||
Evaluation of the total amount of hydrogen formed under the various conditions reveals a temperature of 130 °C with a flow rate of 2.5 mL min−1 as the optimum operating point. At this point, the highest hydrogen yield of 93 mL min−1 with a CO-concentration of 476 ppm was achieved (Fig. 13). The reaction at 160 °C exhibits a CO-concentration of 456 ppm with a H2 productivity of 90 mL min−1. Due to this, operation at 130 °C is advantageous for technical applications due to the improved energy efficiency resulting from the lower reaction temperature.
![]() | ||
| Fig. 13 Evaluation of the real H2 yields at different temperatures and FA flow rates during continuous FAD experiment with Ru@pBINAP/G (left) and Ir@PyrTerpy/G (right) in the PFR. | ||
Similarly, the properties of the terpyridine-based hybrid catalyst Ir@PyrTerpy/G were also investigated in continuous operation. Up until a TON of 40
000 was reached, no significant deactivation was observable, again highlighting the exceptional stability of the system (Fig. S17). After the induction phase, the reaction temperature and FA flow rate were varied akin to Ru@pBINAP/G. At 100 °C and 1.0 mL min−1 FA flowrate, Ir@PyrTerpy/G already exhibits an H2 productivity of 53 mL min−1 with a much better CO-concentration of just 58 ppm. Hydrogen productivity further increased up to 115 mL min−1 at 75 ppm CO at 160 °C (Fig. 13). This corresponds to a TOF of 11
000 h−1, analogous to the batch experiments. ICP-MS reveals an Ir leaching of 0.38 ± 0.27%, revealing that both SMCs enable exceptional metal retention.
Analysing and comparing the CO-concentration at different H2 productivities and temperatures for Ru@pBINAP/G and Ir@PyrTerpy/G (Fig. 14), it can be noted that the target CO-concentration is not achieved by the BINAP-based system Ru@pBINAP/G under any of the reaction conditions investigated. Ir@PyrTerpy/G, on the other hand, has a large temperature and productivity window at which H2 can be released with sufficiently low CO-concentration. The necessary purity of less than 100 ppm CO in the produced H2 for the utilization of more CO-resistant proton exchange membrane (PEM) anodes such as Pt–WO3/C is achieved throughout all of the reaction parameters screened for Ir@PyrTerpy/G.60 Further optimization to the catalyst preparation and the reaction parameters to decrease the amount of CO below 10 ppm is required to allow the use of product gas in integrated energy storage facilities with platinum-based PEM fuel cells.5,61
![]() | ||
| Fig. 14 Concentration of side product CO at different H2 productivities and reaction temperatures of Ru@pBINAP/G (left) and Ir@PyrTerpy/G (right). | ||
200 h−1 and selectivities of up to 330 ppm for a BINAP containing hybrid catalyst, and TOFs up to 11
000 h−1 and CO contents below 100 ppm for the terpyridine system, proving superior over the CTF-analogue. While the BINAP-system showed particularly promising activities in batch experiments, its limited stability proved challenging in subsequent, scaled-up continuous reactions employing a plug-flow reactor. In contrast, the highest H2 productivities of up to 115 mL min−1, no deactivation and good selectivity at 75 ppm CO could be observed utilizing the exceptionally stable terpyridine-based catalyst Ir@PyrTerpy/G at 160 °C. Only very low leaching (0.38%) of the metal precursors could be observed under continuous operation for both systems, highlighting the exceptional metal retention of such catalysts, with important implications regarding costs en route to practical application.
Supplementary information: additional IR-spectra, XRD results, synthetic methodology, catalytic procedure. See DOI: https://doi.org/10.1039/d5cy01382h.
| This journal is © The Royal Society of Chemistry 2026 |