skip to main content
US FlagAn official website of the United States government
dot gov icon
Official websites use .gov
A .gov website belongs to an official government organization in the United States.
https lock icon
Secure .gov websites use HTTPS
A lock ( lock ) or https:// means you've safely connected to the .gov website. Share sensitive information only on official, secure websites.


Title: Ensemble Effects on Allylic Oxidation within Explicit Solvation Environments
Umbrella-sampling density functional theory molecular dynamics (DFT-MD) has been employed to study the full catalytic cycle of the allylic oxidation of cyclohexene using a Cu(II) 7-amino-6-((2- hydroxybenzylidene)amino)quinoxalin-2-ol complex in acetonitrile to create cyclohexenone and H2O as products. After the initial H-atom abstraction step, two different reaction pathways have been identified that are distinguished by the participation of alkyl hydroperoxide (referred to as the “open” cycle) versus the methanol side-product (referred to as the “closed” cycle) within the catalyst recovery process. Importantly, both pathways involve dehydrogenation and re-hydrogenation of the -NH2 group bound to the Cu-site - a feature that is revealed from the ensemble sampling of configurations of the reactive species that are stabilized within the explicit solvent environment of the simulation. Estimation of the energy span from the experimental turnover frequency yields an approximate value of 22.7 kcal/mol at 350 K. Whereas the closed cycle value is predicted to be 26.2 kcal/mol, the open cycle value at 16.5 kcal/mol. Both pathways are further consistent with the equilibrium between Cu(II) and Cu(III) than what has previously been observed. In comparison to prior static DFT calculations, the ensemble of of both solute and solvent configurations has helped to reveal a breadth of processes that underpin the full catalytic cycle yielding a more comprehensive understanding of the importance of radical reactions and catalysis recovery.  more » « less
Award ID(s):
1934725
PAR ID:
10244157
Author(s) / Creator(s):
Date Published:
Journal Name:
Dalton transactions
ISSN:
1477-9226
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. A recent advance in the synthesis of alkenylated arenes was the demonstration that the Pd(OAc)2 catalyst precursor gives >95% selectivity toward styrene from ethylene and benzene under optimized conditions using excess Cu(II) carboxylate as the in situ oxidant [ Organometallics 2019, 38(19), 3532−3541]. To understand the mechanism underlying this catalysis, we applied density functional theory (DFT) calculations in combination with experimental studies. From DFT calculations, we determined the lowest-energy multimetallic Pd and Pd–Cu mixed metal species as possible catalyst precursors. From the various structures, we determined the cyclic heterotrinuclear complex PdCu2(μ-OAc)6 to be the global minimum in Gibbs free energy under conditions of excess Cu(II). For cyclic PdCu2(μ-OAc)6 and the parent [Pd(μ-OAc)2]3, we evaluated the barriers for benzene C–H activation through concerted metalation deprotonation (CMD). The PdCu2(μ-OAc)6 cyclic trimer leads to a CMD barrier of 33.5 kcal/mol, while the [Pd(μ-OAc)2]3 species leads to a larger CMD barrier at >35 kcal/mol. This decrease in the CMD barrier arises from the insertion of Cu(II) into the trimetallic species. Because cyclic PdCu2(μ-OAc)6 is likely the predominant species under experimental conditions (the Cu to Pd ratio is 480:1 at the start of catalysis) with a predicted CMD barrier within the range of the experimentally determined activation barrier, we propose that cyclic PdCu2(μ-OAc)6 is the Pd species responsible for catalysis and report a full reaction mechanism based on DFT calculations. For catalytic conversion of benzene and ethylene to styrene at 120 °C using Pd(OAc)2 as the catalyst precursor and Cu(OPiv)2 (OPiv = pivalate) as the oxidant, an induction period of ∼1 h was observed, followed by catalysis with a turnover frequency of ∼2.3 × 10–3 s–1. In situ1H NMR spectroscopy experiments indicate that during the induction period, Pd(OAc)2 is likely converted to cyclic PdCu2(η2-C2H4)3(μ-OPiv)6, which is consistent with the calculations and consistent with the proposal that the active catalyst is the ethylene-coordinated heterotrinuclear complex cyclic PdCu2(η2-C2H4)3(μ-OPiv)6. 
    more » « less
  2. A Cu or Fe-catalyzed radical ring expansion/coupling cascade was developed. It provides a new access to medium-sized and macrolactones with CN, N3, SCN, and X groups. DFT calculations suggest a Cu(i)/Cu(ii)/Cu(iii) catalytic cycle for the cyanation. 
    more » « less
  3. 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
  4. Aseries of novel 1,4-disubstituted 1,2,3-triazoles were synthesized from an (R)-carvone terminal alkyne derivative via a Cu (I)-catalyzed azide–alkyne cycloaddition reaction using CuSO4,5H2O as the copper (II) source and sodium ascorbate as a reducing agent which reduces Cu (II) into Cu (I). All the newly synthesized 1,2,3-triazoles 9a–h were fully identified on the basis of their HRMS and NMR spectral data and then evaluated for their cell growth inhibition potential by MTS assay against HT-1080 fibrosarcoma, A-549 lung carcinoma, and two breast adenocarcinoma (MCF-7 and MDA-MB-231) cell lines. Compound 9d showed notable cytotoxic effects against the HT-1080 and MCF-7 cells with IC50 values of 25.77 and 27.89 µM, respectively, while compound 9c displayed significant activity against MCF-7 cells with an IC50 value of 25.03 µM. Density functional calculations at the B3LYP/6-31G* level of theory were used to confirm the high reactivity of the terminal alkyne as a dipolarophile. Quantum calculations were also used to investigate the mechanism of both the uncatalyzed and copper (I)-catalyzed azide–alkyne cycloaddition reaction (CuAAC). The catalyzed reaction gives complete regioselectivity via a stepwise mechanism streamlining experimental observations. The calculated free-energy barriers 4.33 kcal/mol and 29.35 kcal/mol for the 1,4- and 1,5-regioisomers, respectively, explain the marked regioselectivity of the CuAAC reaction. 
    more » « less
  5. Fluorocarbons have been shown experimentally by Baker and coworkers to combine with the cyclopentadienylcobalt (CpCo) moiety to form fluoroolefin and fluorocarbene complexes as well as fluorinated cobaltacyclic rings. In this connection density functional theory (DFT) studies on the cyclopentadienylcobalt fluorocarbon complexes CpCo(L)(C n F 2n ) (L = CO, PMe 3 ; n = 3 and 4) indicate structures with perfluoroolefin ligands to be the lowest energy structures followed by perfluorometallacycle structures and finally by structures with perfluorocarbene ligands. Thus, for the CpCo(L)(C 3 F 6 ) (L = CO, PMe 3 ) complexes, the perfluoropropene structure has the lowest energy, followed by the perfluorocobaltacyclobutane structure and the perfluoroisopropylidene structure less stable by 8 to 11 kcal mol −1 , and the highest energy perfluoropropylidene structure less stable by more than 12 kcal mol −1 . For the two metal carbene structures Cp(L)CoC(CF 3 ) 2 and Cp(L)CoCF(C 2 F 5 ), the former is more stable than the latter, even though the latter has Fischer carbene character. For the CpCo(L)(C 4 F 8 ) (L = CO, PMe 3 ) complexes, the perfluoroolefin complex structures have the lowest energies, followed by the perfluorometallacycle structures at 10 to 20 kcal mol −1 , and the structures with perfluorocarbene ligands at yet higher energies more than 20 kcal mol −1 above the lowest energy structure. This is consistent with the experimentally observed isomerization of the perfluorinated cobaltacyclobutane complexes CpCo(PPh 2 Me)(–CFR–CF 2 –CF 2 –) (R = F, CF 3 ) to the perfluoroolefin complexes CpCo(PPh 2 Me)(RCFCF 2 ) in the presence of catalytic quantities of HN(SO 2 CF 3 ) 2 . Further refinement of the relative energies by the state-of-the-art DLPNO-CCSD(T) method gives results essentially consistent with the DFT results summarized above. 
    more » « less