The catalytic effects of iridium pincer complexes for the hydrogenation of carbon dioxide were investigated with the Unified Reaction Valley Approach (URVA), exploring the reaction mechanism along the reaction path traced out by the reacting species on the potential energy surface. Further details were obtained with the Local Mode Analysis performed at all stationary points, complemented by the Natural Bond Orbital and Bader’s Quantum Atoms in Molecules analyses. Each of the five reaction paths forming the catalytic cycle were calculated at the DFT level complemented with DLPNO-CCSD(T) single point calculations at the stationary points. For comparison, the non-catalytic reaction was also investigated. URVA curvature profiles identified all important chemical events taking place in the non-catalyzed reaction and in the five reactions forming the catalytic cycle, and their contribution to the activation energy was disclosed. The non-catalytic reaction has a large unfavorable activation energy of 76.3 kcal/mol, predominately caused by HH bond cleave in the H2 reactant. As shown by our study, the main function of the iridium pincer catalyst is to split up the one–step non-catalytic reaction into an energy efficient multistep cycle, where HH bond cleavage is replaced by the cleavage of a weaker IrH bond with a small contribution to the activation energy. The dissociation of the final product from the catalyst requires the cleavage of an IrO bond, which is also weak, and contributes only to a minor extent to the activation energy. This, in summary, leads to the substantial lowering of the overall activation barrier by about 50 kcal/mol for the catalyzed reaction. We hope that this study inspires the community to add URVA to their repertoire for the investigation of catalysis reactions.
more »
« less
This content will become publicly available on June 8, 2024
Reaction mechanism – explored with the unified reaction valley approach
One of the ultimate goals of chemistry is to understand and manipulate chemical reactions, which implies the ability to monitor the reaction and its underlying mechanism at an atomic scale. In this article, we introduce the Unified Reaction Valley Approach (URVA) as a tool for elucidating reaction mechanisms, complementing existing computational procedures. URVA combines the concept of the potential energy surface with vibrational spectroscopy and describes a chemical reaction via the reaction path and the surrounding reaction valley traced out by the reacting species on the potential energy surface on their way from the entrance to the exit channel, where the products are located. The key feature of URVA is the focus on the curving of the reaction path. Moving along the reaction path, any electronic structure change of the reacting species is registered by a change in the normal vibrational modes spanning the reaction valley and their coupling with the path, which recovers the curvature of the reaction path. This leads to a unique curvature profile for each chemical reaction, with curvature minima reflecting minimal change and curvature maxima indicating the location of important chemical events such as bond breaking/formation, charge polarization and transfer, rehybridization, etc. A decomposition of the path curvature into internal coordinate components or other coordinates of relevance for the reaction under consideration, provides comprehensive insight into the origin of the chemical changes taking place. After giving an overview of current experimental and computational efforts to gain insight into the mechanism of a chemical reaction and presenting the theoretical background of URVA, we illustrate how URVA works for three diverse processes, (i) [1,3] hydrogen transfer reactions; (ii) α-keto-amino inhibitor for SARS-CoV-2 M pro ; (iii) Rh-catalyzed cyanation. We hope that this article will inspire our computational colleagues to add URVA to their repertoire and will serve as an incubator for new reaction mechanisms to be studied in collaboration with our experimental experts in the field.
more »
« less
- Award ID(s):
- 2102461
- NSF-PAR ID:
- 10430211
- Date Published:
- Journal Name:
- Chemical Communications
- Volume:
- 59
- Issue:
- 47
- ISSN:
- 1359-7345
- Page Range / eLocation ID:
- 7151 to 7165
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
-
In this work, we investigated the catalytic effects of a Sharpless dimeric titanium (IV)–tartrate–diester catalyst on the epoxidation of allylalcohol with methyl–hydroperoxide considering four different orientations of the reacting species coordinated at the titanium atom (reactions R1–R4) as well as a model for the non-catalyzed reaction (reaction R0). As major analysis tools, we applied the URVA (Unified Reaction Valley Approach) and LMA (Local Mode Analysis), both being based on vibrational spectroscopy and complemented by a QTAIM analysis of the electron density calculated at the DFT level of theory. The energetics of each reaction were recalculated at the DLPNO-CCSD(T) level of theory. The URVA curvature profiles identified the important chemical events of all five reactions as peroxide OO bond cleavage taking place before the TS (i.e., accounting for the energy barrier) and epoxide CO bond formation together with rehybridization of the carbon atoms of the targeted CC double bond after the TS. The energy decomposition into reaction phase contribution phases showed that the major effect of the catalyst is the weakening of the OO bond to be broken and replacement of OH bond breakage in the non-catalyzed reaction by an energetically more favorable TiO bond breakage. LMA performed at all stationary points rounded up the investigation (i) quantifying OO bond weakening of the oxidizing peroxide upon coordination at the metal atom, (ii) showing that a more synchronous formation of the new CO epoxide bonds correlates with smaller bond strength differences between these bonds, and (iii) elucidating the different roles of the three TiO bonds formed between catalyst and reactants and their interplay as orchestrated by the Sharpless catalyst. We hope that this article will inspire the computational community to use URVA complemented with LMA in the future as an efficient mechanistic tool for the optimization and fine-tuning of current Sharpless catalysts and for the design new of catalysts for epoxidation reactions.more » « less
-
In this study we investigate the Diels–Alder reaction between methyl acrylate and butadiene, which is catalyzed by BF3 Lewis acid in explicit water solution, using URVA and Local Mode Analysis as major tools complemented with NBO, electron density and ring puckering analyses. We considered four different starting orientations of methyl acrylate and butadiene, which led to 16 DA reactions in total. In order to isolate the catalytic effects of the BF3 catalyst and those of the water environment and exploring how these effects are synchronized, we systematically compared the non-catalyzed reaction in gas phase and aqueous solution with the catalyzed reaction in gas phase and aqueous solution. Gas phase studies were performed at the B3LYP/6-311+G(2d,p) level of theory and studies in aqueous solution were performed utilizing a QM/MM approach at the B3LYP/6-311+G(2d,p)/AMBER level of theory. The URVA results revealed reaction path curvature profiles with an overall similar pattern for all 16 reactions showing the same sequence of CC single bond formation for all of them. In contrast to the parent DA reaction with symmetric substrates causing a synchronous bond formation process, here, first the new CC single bond on the CH2 side of methyl acrylate is formed followed by the CC bond at the ester side. As for the parent DA reaction, both bond formation events occur after the TS, i.e., they do not contribute to the energy barrier. What determines the barrier is the preparation process for CC bond formation, including the approach diene and dienophile, CC bond length changes and, in particular, rehybridization of the carbon atoms involved in the formation of the cyclohexene ring. This process is modified by both the BF3 catalyst and the water environment, where both work in a hand-in-hand fashion leading to the lowest energy barrier of 9.06 kcal/mol found for the catalyzed reaction R1 in aqueous solution compared to the highest energy barrier of 20.68 kcal/mol found for the non-catalyzed reaction R1 in the gas phase. The major effect of the BF3 catalyst is the increased mutual polarization and the increased charge transfer between methyl acrylate and butadiene, facilitating the approach of diene and dienophile and the pyramidalization of the CC atoms involved in the ring formation, which leads to a lowering of the activation energy. The catalytic effect of water solution is threefold. The polar environment leads also to increased polarization and charge transfer between the reacting species, similar as in the case of the BF3 catalyst, although to a smaller extend. More important is the formation of hydrogen bonds with the reaction complex, which are stronger for the TS than for the reactant, thus stabilizing the TS which leads to a further reduction of the activation energy. As shown by the ring puckering analysis, the third effect of water is space confinement of the reacting partners, conserving the boat form of the six-member ring from the entrance to the exit reaction channel. In summary, URVA combined with LMA has led to a clearer picture on how both BF3 catalyst and aqueous environment in a synchronized effort lower the reaction barrier. These new insights will serve to further fine-tune the DA reaction of methyl acrylate and butadiene and DA reactions in general.more » « less
-
We report absolute integral cross section (ICS) measurements using a dual-source merged-fast-beams apparatus to study the titular reactions over the relative translational energy range of E r ∼ 0.01–10 eV. We used photodetachment of C − to produce a pure beam of atomic C in the ground electronic 3 P term, with statistically populated fine-structure levels. The H 2 + and D 2 + were formed in an electron impact ionization source, with well known vibrational and rotational distributions. The experimental work is complemented by a theoretical study of the CH 2 + electronic system in the reactant and product channels, which helps to clarify the possible reaction mechanisms underlying the ICS measurements. Our measurements provide evidence that the reactions are barrierless and exoergic. They also indicate the apparent absence of an intermolecular isotope effect, to within the total experimental uncertainties. Capture models, taking into account either the charge-induced dipole interaction potential or the combined charge-quadrupole and charge-induced dipole interaction potentials, produce reaction cross sections that lie a factor of ∼4 above the experimental results. Based on our theoretical study, we hypothesize that the reaction is most likely to proceed adiabatically through the 1 4 A′ and 1 4 A′′ states of CH 2 + via the reaction C( 3 P) + H 2 + ( 2 Σ+g) → CH + ( 3 Π) + H( 2 S). We also hypothesize that at low collision energies only H 2 + ( v ≤ 2) and D 2 + ( v ≤ 3) contribute to the titular reactions, due to the onset of dissociative charge transfer for higher vibrational v levels. Incorporating these assumptions into the capture models brings them into better agreement with the experimental results. Still, for energies ≲0.1 eV where capture models are most relevant, the modified charge-induced dipole model yields reaction cross sections with an incorrect energy dependence and lying ∼10% below the experimental results. The capture cross section obtained from the combined charge-quadrupole and charge-induced dipole model better matches the measured energy dependence but lies ∼30–50% above the experimental results. These findings provide important guidance for future quasiclassical trajectory and quantum mechanical treatments of this reaction.more » « less
-
Abstract In this work, we investigated atmospheric pressure plasma jet (APPJ)-assisted methane oxidation over a Ni-SiO 2 /Al 2 O 3 catalyst. We evaluated possible reaction mechanisms by analyzing the correlation of gas phase, surface and plasma-produced species. Plasma feed gas compositions, plasma powers, and catalyst temperatures were varied to expand the experimental parameters. Real-time Fourier-transform infrared spectroscopy was applied to quantify gas phase species from the reactions. The reactive incident fluxes generated by plasma were measured by molecular beam mass spectroscopy using an identical APPJ operating at the same conditions. A strong correlation of the quantified fluxes of plasma-produced atomic oxygen with that of CH 4 consumption, and CO and CO 2 formation implies that O atoms play an essential role in CH 4 oxidation for the investigated conditions. With the integration of APPJ, the apparent activation energy was lowered and a synergistic effect of 30% was observed. We also performed in-situ diffuse reflectance infrared Fourier-transform spectroscopy to analyze the catalyst surface. The surface analysis showed that surface CO abundance mirrored the surface coverage of CH n at 25 °C. This suggests that CH n adsorbed on the catalyst surface as an intermediate species that was subsequently transformed into surface CO. We observed very little surface CH n absorbance at 500 °C, while a ten-fold increase of surface CO and stronger CO 2 absorption were seen. This indicates that for a nickel catalyst at 500 °C, the dissociation of CH 4 to CH n may be the rate-determining step in the plasma-assisted CH 4 oxidation for our conditions. We also found the CO vibrational frequency changes from 2143 cm −1 for gas phase CO to 2196 cm −1 for CO on a 25 °C catalyst surface, whereas the frequency of CO on a 500 °C catalyst was 2188 cm −1 . The change in CO vibrational frequency may be related to the oxidation of the catalyst.more » « less
