Hanna H.
Cramer‡
a,
Shubhajit
Das‡
b,
Matthew D.
Wodrich
bc,
Clémence
Corminboeuf
*bcd,
Christophe
Werlé
*ae and
Walter
Leitner
*af
aMax Planck Institute for Chemical Energy Conversion, Stiftstr. 34–36, 45470 Mülheim an der Ruhr, Germany. E-mail: christophe.werle@cec.mpg.de; walter.leitner@cec.mpg.de
bLaboratory for Computational Molecular Design Institute of Chemical Sciences and Engineering, École Polytechnique Fédérale de Lausanne (EPFL), 1015 Lausanne, Switzerland. E-mail: clemence.corminboeuf@epfl.ch
cNational Centre for Competence in Research – Catalysis (NCCR-Catalysis), École Polytechnique Fédérale de Lausanne (EPFL), 1015 Lausanne, Switzerland
dNational Centre for Computational Design and Discovery of Novel Materials (MARVEL), Ecole Polytechnique Fédérale de Lausanne (EPFL), 1015 Lausanne, Switzerland
eRuhr University Bochum, Universitätsstr. 150, 44801 Bochum, Germany
fInstitut für Technische und Makromolekulare Chemie (ITMC), RWTH Aachen University, Worringer Weg 2, 52074 Aachen, Germany
First published on 9th February 2023
The stepwise catalytic reduction of carbon dioxide (CO2) to formic acid, formaldehyde, and methanol opens non-fossil pathways to important platform chemicals. The present article aims at identifying molecular control parameters to steer the selectivity to the three distinct reduction levels using organometallic catalysts of earth-abundant first-row metals. A linear scaling relationship was developed to map the intrinsic reactivity of 3d transition metal pincer complexes to their activity and selectivity in CO2 hydrosilylation. The hydride affinity of the catalysts was used as a descriptor to predict activity/selectivity trends in a composite volcano picture, and the outstanding properties of cobalt complexes bearing bis(phosphino)triazine PNP-type pincer ligands to reach the three reduction levels selectively under different reaction conditions could thus be rationalized. The implications of the composite volcano picture were successfully experimentally validated with selected catalysts, and the challenging intermediate level of formaldehyde could be accessed in over 80% yield with the cobalt complex 6. The results underpin the potential of tandem computational-experimental approaches to propel catalyst design for CO2-based chemical transformations.
Fig. 1 Stepwise reduction of CO2 to the formic acid, formaldehyde, and methanol level via hydrosilylation or hydroboration. |
The reduction sequence can be described as a three-step cascade reaction (Fig. 1) that involves the first two-electron reduction to formate derivatives (cycle I), followed by the second and third reduction to acetals (cycle II) and methoxides (cycle III). While many catalytic systems facilitate the reduction to the formate level, those that overcome the kinetic barriers of further reduction to the formaldehyde and methanol levels are scarcer.4 Furthermore, the overreduction of formaldehyde to methanol derivatives is difficult to suppress, and obtaining high selectivity for this important product level remains particularly challenging.5
Transition metal pincer complexes have emerged as promising candidates to overcome these obstacles since they offer many possibilities for precise steric and electronic tuning to fulfill the requirements of the desired catalytic pathway (Fig. 2).6 Early examples are the nickel and cobalt pincer complexes 1 and 2, for which the product distribution was primarily determined by the chosen reductant (Fig. 2A and B). Turculet and co-workers studied a series of group 10 PSiP pincer hydride complexes in the hydroboration of CO2 with HBPin (Fig. 2C).4j While the platinum complexes were largely inactive, palladium analogs were found to generate mainly the formate derivatives, and nickel derivatives lead to the formaldehyde level. A variety of PCP and PSiP pincer hydride complexes catalyzed the hydroboration of CO2 with HBPin mainly to the formate level, as reported by Hazari and co-workers (Fig. 2D, left).4q They showed that (tBuPCP)NiH 3 (Fig. 2D, right) even selectively facilitated the two-, four- or six-electron reduction of CO2, whereby the choice of the reducing agent largely determined the preferred product level. Our group reported a cobalt PNP triazine pincer complex 4 that selectively accessed the formate-, formaldehyde-, or methanol reduction level with the same reducing agent, depending on the appropriate adjustment of the reaction conditions (Fig. 2E).4r,t
Fig. 2 Selected examples of pincer complexes for the catalytic hydroboration and hydrosilylation of CO2. |
While these studies provide valuable mechanistic insight for a few individual catalytic systems, general conclusions about the relationship between catalyst design, activity, and product distribution remain limited. The systematic evaluation of molecular control factors that influence product distribution is of high interest, since it can reveal how the catalyst architecture dictates the preferred reduction level. These insights are very valuable for developing efficient catalysts that (a) specifically target one particular product or (b) allow access to multiple product levels in an adaptive manner. The present study, therefore, goes beyond mechanistic studies on individual catalysts and aims to unravel general structure–reactivity relationships. To achieve this, it combines computational screening using a volcano-type scaling relationship with experimental validation for selected catalysts to contribute to this fundamental challenge in catalyst design for CO2 conversion.
Generally, the relative ratios of the energy spans of cycles I–III determine the observed product distribution in addition to external parameters such as temperature. These energy span ratios will differ significantly for different metals and ligand frameworks, however. Screening potential catalyst candidates would therefore require calculating the free energy profile for all catalyst candidates individually at a high expense of time and computational resources. These computational efforts can be significantly reduced by estimating the free energies of the intermediates and transition states with molecular volcano plots, which relate the catalytic activity to a unique descriptor for each catalyst candidate.7 Calculating the descriptor estimates the entire free energy profile by linear free energy scaling relationships (LFESRs) used to construct the volcano plot. This way, the descriptor value of an individual catalyst allows for a rapid assessment of the relative activity by its position on the plot. Expanding this analysis to multiple catalytic cycles and comparing the corresponding activities for an individual catalyst provides insight into its expected selectivity.
We examined a series of isoelectronic d8 pincer complexes with Fe(0), Co(I), and Ni(II) metal centers selected due to their natural abundance and low costs (computational details can be found in the ESI†).8 The combination of eleven pincer ligands L1–L11 (see Fig. 3B) and three metals furnishes 33 complexes with the general formula [L–M–H]n as I1. The charge n of the complexes results from the combination of neutral or anionic ligands and the formal oxidation state of the metal center. The ligand library includes PNP (L1–L7), NNN (L8), PCP (L9), OCO (L10), and NCN pincer ligands (L11), featuring various substitution patterns which provide catalysts with diverse steric and electronic environments around the metal center. Fig. 4 depicts the Gibbs free energy profile for a representative pincer catalyst, L6-Co-H. The respective oxidative addition products, I3 and I5, undergo an increase of oxidation state by +2 upon their formation.
Fig. 4 Gibbs free energy profile of L6-Co-H (SMD(benzene)/M06/def2-TZVP//M06/def2-SVPD; computational details are available in the ESI†). |
The hydride affinity (ΔGH-) was chosen as a descriptor to represent the intrinsic reactivity of the complexes to map out the activity of the pincer complexes towards the individual cycles and hence also selectivity. The correlation of ΔGH- with the catalytic performance of transition metal hydrides, including hydride transfer and hydrogen splitting barriers in the hydrogenation of CO2 to formate, has been demonstrated in previous studies.9 This work establishes an expanded theoretical framework that features the hydride affinity as a descriptor for multiple product levels beyond formate, while it provides quantitative information not only on the activity but also on the selectivity. ΔGH- can be easily estimated from the free energy change of the hydride exchange equilibrium concerning the trityl cation as the reference hydride acceptor (eqn (1)). According to the sign convention used here, a more negative ΔGH- value indicates a weaker metal–hydride bond.
I1 + [Ph3C]+ ⇌ [I1−H]+ + Ph3CH |
ΔGH− = (ΔG[I1−H]+ + ΔGPh3CH) − (ΔGI1 + ΔG[Ph3C]+) | (1) |
Providing a measure of the metal–hydride interaction strength, ΔGH- is connected with the relative Gibbs free energies for the intermediates and transition states of all three cycles by LFESRs (see Fig. S1†), as exemplified for I3 in Fig. 5A. Note that it is a solely catalyst-based property and, thus, independent of a specific chemical transformation which not only streamlines the screening process but also showcases the transferability of our approach. For the present case, ΔGH- can also act as a descriptor in the reduction of formaldehyde and methanol derivatives, including elementary steps that do not directly involve a hydride, such as the reductive eliminations involving TS3, TS6, or TS7. As shown in Fig. 5A, the descriptor ΔGH- gradually decreases for catalysts with a more negative charge, as expected from a weaker metal hydride bond. Accordingly, anionic NCN, OCO, and PCP ligands exhibit ΔGH- significantly more negative relative to the ligand L1 for each metal (Fig. 5B). While the differences are less pronounced for ligands with the same charge, the influence of the electronic properties of the coordinating atoms is apparent. For example, replacing triazine (L1) with the stronger donor pyridine (L4) leads to a more negative ΔGH- in agreement with a higher electron density at the metal and a weakened M–H bond. Replacing phosphine (L3) with nitrogen arms (L8) leads to a substantial decrease of ΔGH-, reflecting the higher basicity of amines compared to phosphines. In contrast, electronegative oxygen atoms, as in L2 and L7, or para-CF3 substitution, as in L5, induce a higher ΔGH-.
Three individual volcano plots (one for each cycle) are derived from the LFESRs (Fig. 6) to analyze the general trends in hydrosilylation activity and selectivity. The step-by-step procedure for constructing molecular volcano plots can be found in a previously published extensive and general protocol.10 While LFESRs and volcano plots were recently employed in homogeneous catalysis by some of us,7j,k,11 Norskov and co-workers have previously applied them in the fields of solid-state heterogeneous catalysis and electrocatalysis.12 Here, ΔGH- is plotted along the X-axes of these volcanoes while the Y-axes correspond to the negative of the respective catalytic cycle's energy span (δG), as defined by eqn (2).13 δG is calculated from the relative energies TTDTS and ITDI of the turnover-determining transition state (TDTS) and of the turnover-determining intermediate (TDI), respectively, which are the pair of transition state and intermediate that constitutes the largest barrier in the catalytic cycle. The total Gibbs free energy ΔGR is added if the TDI appears after the TDTS in the catalytic cycle. Consequently, the catalysts having a lower energy span are more active than the others and are positioned higher in the volcano plots of Fig. 6.
(2) |
Each volcano plot in Fig. 6 can be divided into different regions depending on the combinations of TDI/TDTS that determine the overall barrier under the framework of the energy span model in the cycle. The oxidative addition products I3 and I5 on the left side of the volcano (more negative ΔGH-) are too strongly stabilized and thus control the activity together with the high-lying transition states TS3 and TS6. These features align with the expected behavior for catalysts with a high electron density at the metal center since the stabilization of the oxidative addition products renders the forward step of reductive elimination energetically unfavorable. In contrast, complexes with a low electron density at the metal center (more positive or positive ΔGH- values) have stronger metal–hydride bonds, and thus, I1 dominates the right sides of the volcano plots. The volcano tops for cycles I (Fig. 6A) and II (Fig. 6B) exhibit peaks at ΔGH- = −48 kcal mol−1 and –73 kcal mol−1, respectively. This suggests that electron-rich complexes are more active in catalyzing the second reduction step in comparison to the first one. Cycle III (Fig. 6C) features a plateau (ΔGH- = −99 kcal mol−1 to −50 kcal mol−1) that is governed by I5 and TS7 involved in the formation of formaldehyde and diphenylsiloxane. This reaction step is associated with bond-forming and breaking events that do not occur in direct proximity to the metal center, and therefore it is less affected by the electronic properties of the metal.
Since the overall product yield of all three C1 products is mainly determined by the activity of the first reduction step of CO2 to formate, candidates that appear close to the volcano peak of the cycle I would be expected to form the most productive catalysts. This area includes the neutral complexes of cobalt and nickel L1-Co-H, L4-Co-H, and L11-Ni-H. The complexes 2, 3, and 4 that were previously reported as active CO2 hydrosilylation and hydroboration catalysts (Fig. 2) structurally resemble these complexes in accordance with the predictions from the volcano plots. Notably, the catalyst L9-Ni-H appears at a larger distance to the volcano top and is thus expected to be less active. This aligns with the observation that the related complex 1 was reported inactive in the hydrosilylation of CO2 with phenylsilane at room temperature.4d
To obtain insights into the selectivity of pincer catalysts in CO2 hydrosilylation, we superimposed the three graphs to produce a single composite volcano plot (Fig. 6D). The non-overlaying positions of the three volcano tops emphasize that the ideal hydride affinity value is distinct for each product level. Based on the relative magnitude of the energy spans of the three cycles δGI, δGII, and δGIII, the composite volcano can be divided into two distinct regions. In the green region (i.e., –45 kcal mol−1 < ΔGH- < 33 kcal mol−1), the energy span of the cycles gradually increases in the order I < II < III, suggesting that stepwise access to the three product levels should be feasible by control of external parameters. The sequential ascending energy barriers in the three cycles imply that the catalysts located in this region will be responsive to changes in reaction temperature: At elevated temperature, the selectivity could gradually shift from the formate to the formaldehyde and the methanol level. As the hydride affinity goes down, the decrease in δGII and δGIII is steeper (indicated by the higher slope of the volcano lines) compared to δGI, gradually favoring the formation of bis(silyl)acetal and formaldehyde. Below ΔGH- < −52 kcal mol−1 (blue region), the relative span ordering changes to I > II ≈ III, suggesting that formate production is more challenging than the subsequent formation of further reduced products. Thus, any generated silyl formate is expected to be rapidly further reduced to bis(silyl)acetal and formaldehyde, preventing isolation of the silyl formate independently of the temperature.
Fig. 7 Top region of three volcano plots with indicated positions of selected catalyst candidates (A) and practical application of selected catalysts in CO2 hydrosilylation (B). |
In the case of the nickel NCN pincer complex 5, the overall yield for C1 products drastically increases with the temperature from 3% (r.t.) to 58% (80 °C). This suggests that forming the active hydride complex from precursor 5 may require higher temperatures or longer reaction times. In accordance with the prediction, 5 turns out to be more suitable than the related PCP complex 1 which was reported inactive in the hydrosilylation of CO2 with phenylsilane.4d In contrast to the use of the cobalt catalysts 6, 4, and 7, only traces of silyl formate are observed with complex 5 but mixtures of acetal and methoxide are obtained directly. The relatively highest activity for reduction beyond the formate level is in line with the lowest hydride affinity of 5 among the four examined catalysts. To increase the catalyst selectivity towards the methanol level, consequently, an ideal catalyst would employ an even lower hydride affinity and appear close to the volcano tops of all cycles. This would apply to the start of the plateau region of cycle III at −50 kcal mol−1. As shown previously for complex 4,4r the methanol product yield might be increased additionally by employing a higher reaction temperature and/or longer reaction time.
The cobalt complexes 6, 4, and 7 reach overall yields between 83% and 94% of C1 reduction products already at room temperature as predicted from their position on the volcano plot. The selectivity is strongly temperature-dependent and the quantities of bis(silyl)acetals and methoxysilanes increase at the expense of the silyl formate yield at higher temperatures in all cases, which is in line with the relative energy spans of the three catalytic cycles I–III. Notably, the reduction capability beyond the formate level drastically increases in the series of 7, 4, and 6, as reflected in room temperature yields for bis(silyl)acetals increasing substantially from 6% (7), 24% (4), to 47% (6). In particular, complex 6 reaches 80% and 81% yield towards the formaldehyde level as the most challenging product at 60 °C and 80 °C, respectively. This significant improvement over the previously reported complex 4 would have been difficult to predict using conventional ligand effects due to the complexity of the reaction network. Notably, only a few catalytic systems are reported to facilitate the reduction of CO2 to formaldehyde and they often require precious metals, elevated pressure, (Lewis) acidic additives, or long reaction times.4f,g,4n,14a–d In contrast, the base-metal catalyst 6 operates at ambient pressure, in the absence of additives, and reaches 81% yield at 80 °C in only 2 h.
Footnotes |
† Electronic supplementary information (ESI) available: General considerations, experimental methods, synthetic details, and copies of NMR spectra. See DOI: https://doi.org/10.1039/d2sc06793e |
‡ H. H. C. and S. D. contributed equally to this work. |
This journal is © The Royal Society of Chemistry 2023 |