skip to main content

Explore Research Products in the NSF-PAR

Title: URVA and Local Mode Analysis of an Iridium Pincer Complex Efficiently Catalyzing the Hydrogenation of Carbon Dioxide
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
Award ID(s):
2102461
NSF-PAR ID:
10430205
Author(s) / Creator(s):
;
Date Published:
Journal Name:
Inorganics
Volume:
10
Issue:
12
ISSN:
2304-6740
Page Range / eLocation ID:
234
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. ; ( , Catalysts)
    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. ; ( , Catalysts)
    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
  3. ; ; ; ; ; ( , Nanoscale)
    null (Ed.)
    Palladium catalyzed cross-coupling reactions represent a significant advancement in contemporary organic synthesis as these reactions are of strategic importance in the area of pharmaceutical drug discovery and development. Supported palladium-based catalysts are highly sought-after in carbon–carbon bond forming catalytic processes to ensure catalyst recovery and reuse while preventing product contamination. This paper reports the development of heterogeneous Pd-based bimetallic catalysts supported on fumed silica that have high activity and selectivity matching those of homogeneous catalysts, eliminating the catalyst's leaching and sintering and allowing efficient recycling of the catalysts. Palladium and base metal (Cu, Ni or Co) contents of less than 1.0 wt% loading are deposited on a mesoporous fumed silica support (surface area SA BET = 350 m 2 g −1 ) using strong electrostatic adsorption (SEA) yielding homogeneously alloyed nanoparticles with an average size of 1.3 nm. All bimetallic catalysts were found to be highly active toward Suzuki cross-coupling (SCC) reactions with superior activity and stability for the CuPd/SiO 2 catalyst. A low CuPd/SiO 2 loading (Pd: 0.3 mol%) completes the conversion of bromobenzene and phenylboronic acid to biphenyl in 30 minutes under ambient conditions in water/ethanol solvent. In contrast, monometallic Pd/SiO 2 (Pd: 0.3 mol%) completes the same reaction in three hours under the same conditions. The combination of Pd with the base metals helps in retaining the Pd 0 status by charge donation from the base metals to Pd, thus lowering the activation energy of the aryl halide oxidative addition step. Along with its exceptional activity, CuPd/SiO 2 exhibits excellent recycling performance with a turnover frequency (TOF) of 280 000 h −1 under microwave reaction conditions at 60 °C. Our study demonstrates that SEA is an excellent synthetic strategy for depositing ultra-small Pd-based bimetallic nanoparticles on porous silica for SCC. This avenue not only provides highly active and sintering-resistant catalysts but also significantly lowers Pd contents in the catalysts without compromising catalytic activity, making the catalysts very practical for large-scale applications. 
    more » « less
  4. ; ; ( , Molecules)
    null (Ed.)
    The formation of amide bonds represents one of the most fundamental processes in organic synthesis. Transition-metal-catalyzed activation of acyclic twisted amides has emerged as an increasingly powerful platform in synthesis. Herein, we report the transamidation of N-activated twisted amides by selective N–C(O) cleavage mediated by air- and moisture-stable half-sandwich Ni(II)–NHC (NHC = N-heterocyclic carbenes) complexes. We demonstrate that the readily available cyclopentadienyl complex, [CpNi(IPr)Cl] (IPr = 1,3-bis(2,6-diisopropylphenyl)imidazol-2-ylidene), promotes highly selective transamidation of the N–C(O) bond in twisted N-Boc amides with non-nucleophilic anilines. The reaction provides access to secondary anilides via the non-conventional amide bond-forming pathway. Furthermore, the amidation of activated phenolic and unactivated methyl esters mediated by [CpNi(IPr)Cl] is reported. This study sets the stage for the broad utilization of well-defined, air- and moisture-stable Ni(II)–NHC complexes in catalytic amide bond-forming protocols by unconventional C(acyl)–N and C(acyl)–O bond cleavage reactions. 
    more » « less
  5. ; ; ; ; ( , Chemical Science)
    Pyridine and quinoline undergo selective C–H activation in the 2-position with Rh and Ir complexes of a boryl/bis(phosphine) PBP pincer ligand, resulting in a 2-pyridyl bridging the transition metal and the boron center. Examination of this reactivity with Rh and Ir complexes carrying different non-pincer ligands on the transition metal led to the realization of the possible isomerism derived from the 2-pyridyl fragment connecting either via B–N/C–M bonds or via B–C/N–M bonds. This M–C/M–N isomerism was systematically examined for four structural types. Each of these types has a defined set of ligands on Rh/Ir besides 2-pyridyl and PBP. A pair of M–C/M–N isomers for each type was computationally examined for Rh and for Ir, totaling 16 compounds. Several of these compounds were isolated or observed in solution by experimental methods, in addition to a few 2-quinolyl variants. The DFT predictions concerning the thermodynamic preference within each M–C/M–N isomeric match the experimental findings very well. In two cases where DFT predicts <2 kcal mol −1 difference in free energy, both isomers were experimentally observed in solution. Analysis of the structural data, of the relevant Wiberg bond indices, and of the ETS-NOCV partitioning of the interaction of the 2-pyridyl fragment with the rest of the molecule points to the strength of the M–C(pyridyl) bond as the dominant parameter determining the relative M–C/M–N isomer favorability. This M–C bond is always stronger for the analogous Ir vs. Rh compounds, but the nature of the ligand trans to it has a significant influence, as well. DFT calculations were used to evaluate the mechanism of isomerization for one of the molecule types. 
    more » « less
  • Citation Formats
  • MLA
  • APA
  • Chicago
  • BibTeX
  • Save / Share this Record
  • Send to Email