skip to main content


Title: BF3–Catalyzed Diels–Alder Reaction between Butadiene and Methyl Acrylate in Aqueous Solution—An URVA and Local Vibrational Mode Study
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
Award ID(s):
2102461
NSF-PAR ID:
10342539
Author(s) / Creator(s):
;
Date Published:
Journal Name:
Catalysts
Volume:
12
Issue:
4
ISSN:
2073-4344
Page Range / eLocation ID:
415
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. 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
  2. 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
  3. We report and interpret recently recorded high-resolution infrared spectra for the fundamentals of the CH 2 scissors and CH stretches of gas phase cyclopentane at −26.1 and −50 °C, respectively. We extend previous theoretical studies of this molecule, which is known to undergo barrierless pseudorotation due to ring puckering, by constructing local mode Hamiltonians of the stretching and scissor vibrations for which the frequencies, couplings, and linear dipoles are calculated as functions of the pseudorotation angle using B3LYP/6-311++(d,p) and MP2/cc-pVTZ levels of theory. Symmetrization ( D 5 h ) of the vibrational basis sets leads to simple vibration/pseudorotation Hamiltonians whose solutions lead to good agreement with the experiment at medium resolution, but which miss interesting line fractionation when compared to the high-resolution spectra. In contrast to the scissor motion, pseudorotation leads to significant state mixing of the CH stretches, which themselves are Fermi coupled to the scissor overtones. 
    more » « less
  4. The Diels–Alder (DA) reaction, a classic cycloaddition reaction involving a diene and a dienophile to form a cyclohexene, is among the most versatile organic reactions. Theories have predicted thermodynamically unfavorable DA reactions on pristine graphene owing to its low chemical reactivity. We hypothesized that metals like Ni could enhance the reactivity of graphene towards DA reactions through charge transfer. The results indeed showed that metal substrates enhanced the reactivity of graphene in the DA reactions with a diene, 2,3-dimethoxy butadiene (DMBD), and a dienophile, maleic anhydride (MAH), with the activity enhancement in the order of Ni > Cu, and both are more reactive than graphene supported on silicon wafer. The rate constants were estimated to be two times higher for graphene supported on Ni than on silicon wafer. The computational results support the experimentally obtained rate trend of Ni > Cu, both predicted to be greater than unsupported graphene, which is explained by the enhanced graphene–substrate interaction reflected in charge transfer effects with the strongly interacting Ni. This study opens up a new avenue for enhancing the chemical reactivity of pristine graphene through substrate selection. 
    more » « less
  5. Per- and polyfluoroalkyl substances (PFAS) are a group of synthetic organofluorine surfactants that are resistant to typical methods of degradation. Thermal techniques along with other novel, less energy-intensive techniques are currently being investigated for the treatment of PFAS-contaminated matrices. Non-equilibrium plasma is one technique that has shown promise for the treatment of PFAS-contaminated water. To better tailor non-equilibrium plasma systems for this application, knowledge of the energy required for mineralization, and in turn the roles that plasma reactive species and heat can play in this process, would be useful. In this study, fundamental thermodynamic equations were used to estimate the enthalpies of reaction (480 kJ/mol) and formation (−4640 kJ/mol) of perfluorooctanoic acid (PFOA, a long-chain legacy PFAS) in water. This enthalpy of reaction estimate indicates that plasma reactive species alone cannot catalyze the reaction; because the reaction is endothermic, energy input (e.g., heat) is required. The estimated enthalpies were used with HSC Chemistry software to produce a model of PFOA defluorination in a 100 mg/L aqueous solution as a function of enthalpy. The model indicated that as enthalpy of the reaction system increased, higher PFOA defluorination, and thus a higher extent of mineralization, was achieved. The model results were validated using experimental results from the gliding arc plasmatron (GAP) treatment of PFOA or PFOS-contaminated water using argon and air, separately, as the plasma gas. It was demonstrated that PFOA and PFOS mineralization in both types of plasma required more energy than predicted by thermodynamics, which was anticipated as the model did not take kinetics into account. However, the observed trends were similar to that of the model, especially when argon was used as the plasma gas. Overall, it was demonstrated that while energy input (e.g., heat) was required for the non-equilibrium plasma degradation of PFOA in water, a lower energy barrier was present with plasma treatment compared to conventional thermal treatments, and therefore mineralization was improved. Plasma reactive species, such as hydroxyl radicals (⋅OH) and/or hydrated electrons (e−(aq)), though unable to accelerate an endothermic reaction alone, likely served as catalysts for PFOA mineralization, helping to lower the energy barrier. In this study, the activation energies (Ea) for these species to react with the alpha C–F bond in PFOA were estimated to be roughly 1 eV for hydroxyl radicals and 2 eV for hydrated electrons.

     
    more » « less