Unusual KIE and dynamics effects in the Fe-catalyzed hetero-Diels-Alder reaction of unactivated aldehydes and dienes

A Synergistic Bimetallic Ti/Co-Catalyzed Isomerization of Epoxides to Allylic Alcohols Enabled by Two-State Reactivity

Isomerization of epoxides into versatile allylic alcohols is an atom-economical synthetic method to afford vicinal bifunctional groups. Comprehensive density functional theory (DFT) calculations were carried out to elucidate the complex mechanism of a bimetallic Ti/Co-catalyzed selective isomerization of epoxides to allyl alcohols by examining several possible pathways. Our results suggest a possible mechanism involving (1) radical-type epoxide ring opening catalyzed by Cp2Ti(III)Cl leading to a Ti(IV)-bound β-alkyl radical, (2) hydrogen-atom transfer (HAT) catalyzed by the Co(II) catalyst to form the Ti(IV)-enolate and Co(III)-H intermediate, (3) protonation to give the alcohols, and (4) proton abstraction to form the Co(I) species followed by electron transfer to regenerate the active Co(II) and Ti(III) species. Moreover, bimetallic catalysis and two-state reactivity enable the key rate-determining HAT step. Furthermore, a subtle balance between dispersion-driven bimetallic processes and entropy-driven monometallic processes determines the most favorable pathway, among which the monometallic process is energetically more favorable in all steps except the vital hydrogen-atom transfer step. Our study should provide an in-depth mechanistic understanding of bimetallic catalysis.

Running Wild through Dirhodium Tetracarboxylate-Catalyzed Combined CH(C)-Functionalization/Cope Rearrangement Landscapes: Does Post-Transition-State Dynamic Mismatching Influence Product Distributions?

A special type of C-H functionalization can be achieved through C-H insertion combined with Cope rearrangement (CHCR) in the presence of dirhodium catalysts. This type of reaction was studied using density functional theory and ab initio molecular dynamics simulations, the results of which pointed to the dynamic origins of low yields observed in some experiments. These studies not only reveal intimate details of the complex reaction network underpinning CHCR reactions but also further cement the generality of the importance of nonstatistical dynamic effects in controlling Rh2L4-promoted reactions.

Molecular Dynamics of Iron Porphyrin-Catalyzed C-H Hydroxylation of Ethylbenzene

Quasi-classical molecular dynamics (MD) simulations were carried out to study the mechanism of iron porphyrin-catalyzed hydroxylation of ethylbenzene. The hydrogen atom abstraction from ethylbenzene by iron-oxo species is the rate-determining step, which generates the radical pair of iron-hydroxo species and the benzylic radical. In the subsequent radical rebound step, the iron-hydroxo species and benzylic radical recombine to form the hydroxylated product, which is barrierless on the doublet energy surface. In the gas-phase quasi-classical MD study on the doublet energy surface, 45% of the reactive trajectories lead directly to the hydroxylated product, and this increases to 56% in implicit solvent model simulations. The percentage of reactive trajectories leading to the separated radical pair is 98-100% on high-spin (quartet/sextet) energy surfaces. The low-spin state reactivity dominates in the hydroxylation of ethylbenzene, which is dynamically both concerted and stepwise, since the time gap between C-H bond cleavage and C-O bond formation ranges from 41 to 619 fs. By contrast, the high-spin state catalysis is an energetically stepwise process, which has a negligible contribution to the formation of hydroxylation products.

Investigation of Iron(III) Tetraphenylporphyrin as a Redox Flow Battery Anolyte: Unexpected Side Reactivity with the Electrolyte

Redox flow batteries (RFBs) present an opportunity to bridge the gap between the intermittent availability of green energy sources and the need for on-demand grid level energy storage. While aqueous vanadium-based redox flow batteries have been commercialized, they are limited by the constraints of using water as an electrochemical solvent. Nonaqueous redox flow battery systems can be used to produce high voltage batteries due to the larger electrochemical window in nonaqueous solvents and the ability to tune the redox properties of active materials through functionalization. Iron porphyrins, a class of organometallic macrocycles, have been the subject of many studies for their photocatalytic and electrocatalytic properties in nonaqueous solvents. Often, iron porphyrins can undergo multiple redox events making them interesting candidates for use as anolytes in asymmetrical redox flow batteries or as both catholyte and anolyte in symmetrical redox flow battery systems. Here the electrochemical properties of Fe(III)TPP species relevant to redox flow battery electrolytes are investigated including solubility, electrochemical properties, and charge/discharge cycling. Commonly used support electrolyte salts can have reactivities that are often overlooked beyond their conductivity properties in nonaqueous solvents. Parasitic reactions with the cations of common support electrolytes are highlighted herein, which underscore the careful balance required to fully assess the potential of novel RFB electrolytes.

Electron-rich benzofulvenes as effective dipolarophiles in copper(i)-catalyzed asymmetric 1,3-dipolar cycloaddition of azomethine ylides

A series of benzofulvenes without any electron-withdrawing substituents were employed as 2π-type dipolarophiles for the first time to participate in Cu(i)-catalyzed asymmetric 1,3-dipolar cycloaddition (1,3-DC) reactions of azomethine ylides. An intrinsic non-benzenoid aromatic characteristic from benzofulvenes serves as a key driving force for activation of the electron-rich benzofulvenes. Utilizing the current methodology, a wide range of multi-substituted chiral spiro-pyrrolidine derivatives containing two contiguous all-carbon quaternary centers were formed in good yield with exclusive chemo-/regioselectivity and high to excellent stereoselectivity. Computational mechanistic studies elucidate the origin of the stereochemical outcome and the chemoselectivity, in which the thermostability of these cycloaddition products is the major factor.

Origin of the Second-Order Proton Catalysis of Ferriin Reduction in Belousov-Zhabotinsky Reactions: Density Functional Studies of Ferroin and Ferriin Aggregates with Outer Sphere Ligands Sulfate, Bisulfate, and Sulfuric Acid

The detailed mechanisms of Belousov-Zhabotinsky oscillating reactions continue to present grand challenges, even after half a century of study. The origin of the pH dependence of the oscillation pattern had never been rigorously identified. In our recent kinetic study of one of the key Belousov-Zhabotinsky reactions, the iron-catalyzed bromate oxidation of malonic acid, compelling agreement between experiments and kinetic simulations was achieved only with the inclusion of second-order proton catalysis of the reduction of the [Fe(phen)₃]³⁺ species. After exhausting all other avenues in search of an explanation of this proton catalysis, we considered the possibility that the parent iron-phenanthroline complexes could aggregate with neutral and anionic outer sphere ligands (OSLs) in the highly concentrated sulfuric acid solution, and we hypothesized that OSL protonation would increase the capacity of the aggregated complex to oxidize the organic fuel. We performed potential energy surface analyses at the SMD(APFD/6-311G*) level of complexes of the types [Fe(phen)₃(SO₄^2-)_m(HSO₄^-)_n(H₂SO₄)_o]^(c-2m-n)+ for ferriin (c = 3) and ferroin (c = 2) aggregated with m sulfate, n bisulfate, and o sulfuric acid OSLs. We present structures of the OSL aggregates, develop a nomenclature for their description, and characterize their electronic structure. The structural chemistry provides the foundation to discuss the ferroin/ferriin redox couple with emphasis on the relationship between the vertical electron affinities of ferriin aggregates and their OSL protonation states. For proton catalysis to manifest itself, double-protonation paths that are slightly endergonic should be present, and proton affinities of aggregated OSLs allow the identification of such double-protonation chains. As a first test of our mechanistic proposal for the second-order proton catalysis of the Belousov-Zhabotinsky reaction, the results presented here provide compelling evidence in support of the importance of outer sphere ligation of the iron catalyst.

Enantioselective oxa-Diels-Alder Sequences of Dendralenes

Diene-transmissive hetero-Diels-Alder sequences involving carbonyl dienophiles are reported for the first time. High enantioselectivities are achieved in the reaction of phenylglyoxal with a broad range of dendralene structures, through the optimization of a Pd²⁺ catalyst system. The initial catalyst-controlled enantioselective oxa-Diels-Alder (ODA) cycloaddition to a [3]dendralene generates a dihydropyran carrying a semicyclic diene. This participates in a subsequent catalyst or substrate-controlled Diels-Alder reaction to generate sp³ -rich fused polycyclic systems containing both heterocycles and carbocycles. Computational investigations reveal a concerted asynchronous mechanism. π-Complexation of a diene C=C bond to Pd²⁺ occurs in both the pre-transition state (TS) complex and in cycloaddition TSs, controlling stereoselectivity. A formal enantioselective [4+2]cycloaddition of a CO₂ dienophile is demonstrated.

New Insights and Predictions into Complex Homogeneous Reactions Enabled by Computational Chemistry in Synergy with Experiments: Isotopes and Mechanisms

Homogeneous catalysis and biocatalysis have been widely applied in synthetic, medicinal, and energy chemistry as well as synthetic biology. Driven by developments of new computational chemistry methods and better computer hardware, computational chemistry has become an essentially indispensable mechanistic "instrument" to help understand structures and decipher reaction mechanisms in catalysis. In addition, synergy between computational and experimental chemistry deepens our mechanistic understanding, which further promotes the rational design of new catalysts. In this Account, we summarize new or deeper mechanistic insights (including isotope, dispersion, and dynamical effects) into several complex homogeneous reactions from our systematic computational studies along with subsequent experimental studies by different groups. Apart from uncovering new mechanisms in some reactions, a few computational predictions (such as excited-state heavy-atom tunneling, steric-controlled enantioswitching, and a new geminal addition mechanism) based on our mechanistic insights were further verified by ensuing experiments.The Zimmerman group developed a photoinduced triplet di-π-methane rearrangement to form cyclopropane derivatives. Recently, our computational study predicted the first excited-state heavy-atom (carbon) quantum tunneling in one triplet di-π-methane rearrangement, in which the reaction rates and ¹²C/¹³C kinetic isotope effects (KIEs) can be enhanced by quantum tunneling at low temperatures. This unprecedented excited-state heavy-atom tunneling in a photoinduced reaction has recently been verified by an experimental ¹²C/¹³C KIE study by the Singleton group. Such combined computational and experimental studies should open up opportunities to discover more rare excited-state heavy-atom tunneling in other photoinduced reactions. In addition, we found unexpectedly large secondary KIE values in the five-coordinate Fe(III)-catalyzed hetero-Diels-Alder pathway, even with substantial C-C bond formation, due to the non-negligible equilibrium isotope effect (EIE) derived from altered metal coordination. Therefore, these KIE values cannot reliably reflect transition-state structures for the five-coordinate metal pathway. Furthermore, our density functional theory (DFT) quasi-classical molecular dynamics (MD) simulations demonstrated that the coordination mode and/or spin state of the iron metal as well as an electric field can affect the dynamics of this reaction (e.g., the dynamically stepwise process, the entrance/exit reaction channels).Moreover, we unveiled a new reaction mechanism to account for the uncommon Ru(II)-catalyzed geminal-addition semihydrogenation and hydroboration of silyl alkynes. Our proposed key gem-Ru(II)-carbene intermediates derived from double migrations on the same alkyne carbon were verified by crossover experiments. Additionally, our DFT MD simulations suggested that the first hydrogen migration transition-state structures may directly and quickly form the key gem-Ru-carbene structures, thereby "bypassing" the second migration step. Furthermore, our extensive study revealed the origin of the enantioselectivity of the Cu(I)-catalyzed 1,3-dipolar cycloaddition of azomethine ylides with β-substituted alkenyl bicyclic heteroarenes enabled by dual coordination of both substrates. Such mechanistic insights promoted our computational predictions of the enantioselectivity reversal for the corresponding monocyclic heteroarene substrates and the regiospecific addition to the less reactive internal C═C bond of one diene substrate. These predictions were proven by our experimental collaborators. Finally, our mechanistic insights into a few other reactions are also presented. Overall, we hope that these interactive computational and experimental studies enrich our mechanistic understanding and aid in reaction development.

Post-spin crossing dynamics determine the regioselectivity in open-shell singlet biradical recombination

Comprehensive computational studies reveal unique dynamic effects in a multi-spin-state reaction that determine the regioselectivity of a biradical recombination process.

Quasiclassical Direct Dynamics Trajectory Simulations of Organometallic Reactions

Homogeneous metal-mediated organometallic reactions represent a very large and diverse reaction class. Density functional theory calculations are now routinely carried out and reported for analyzing organometallic mechanisms and reaction pathways. While density functional theory calculations are extremely powerful to understand the energy and structure of organometallic reactions, there are several assumptions in their use and interpretation to define reaction mechanisms and to analyze reaction selectivity. Almost always it is assumed that potential energy structures calculated with density functional theory adequately describe mechanisms and selectivity within the framework of statistical theories, for example, transition state theory and RRKM theory. However, these static structures and corresponding energy landscapes do not provide atomic motion information during reactions that could reveal nonstatistical intermediates without complete intramolecular vibrational redistribution and nonintrinsic reaction coordinate (non-IRC) pathways. While nonstatistical intermediates and non-IRC reaction pathways are now relatively well established for organic reactions, these dynamic effects have heretofore been highly underexplored in organometallic reactions. Through a series of quasiclassical density functional theory direct dynamics trajectory studies, my group has recently demonstrated that dynamic effects occur in a variety of fundamental organometallic reactions, especially bond activation reactions. For example, in the C-H activation reaction between methane and [Cp*(PMe₃)Ir^III(CH₃)]⁺, while the density functional theory energy landscape showed a two-step oxidative cleavage and reductive coupling mechanism, trajectories revealed a mixture of this two-step mechanism and a dynamic one-step mechanism that skipped the [Cp*(PMe₃)Ir^V(H)(CH₃)₂]⁺ intermediate. This study also showed that despite a methane σ-complex being located on the density functional theory surface before oxidative cleavage and after reductive coupling, this intermediate is always skipped and should not be considered an intermediate during reactive trajectories. For non-IRC reaction pathways, quasiclassical direct dynamics trajectories showed that for the isomerization of [Tp(NO)(PMe₃)W(η²-benzene)] to [Tp(NO)(PMe₃)W(H)(Ph)], there are many dynamic reaction pathway connections due to a relatively flat energy landscape and π coordination is not necessary for C-H bond activation through oxidative cleavage. Trajectories also showed that dynamic effects are important in selectivity for ethylene C-H activation versus π coordination in reaction with Cp(PMe₃)₂Re, and trajectories provide a more quantitative model of selectivity than transition state theory. Quasiclassical trajectories examining Au-catalyzed monoallylic diol cyclizations showed dynamic coupling of several reaction steps that include alkoxylation π bond addition, proton shuttling, and water elimination reaction steps. Overall, these studies highlight the need to use direct dynamics trajectory simulations to consider atomic motion during reactions to understand organometallic reaction mechanisms and selectivity.

On the Question of Stepwise [4+2] Cycloaddition Reactions and Their Stereochemical Aspects

Even at the end of the twentieth century, the view of the one-step [4+2] cycloaddition (Diels-Alder) reaction mechanism was widely accepted as the only possible one, regardless of the nature of the reaction components. Much has changed in the way these reactions are perceived since then. In particular, multi-step mechanisms with zwitterionic or diradical intermediates have been proposed for a number of processes. This review provided a critical analysis of such cases.

Selective oxidation of alcohol-d 1 to aldehyde-d 1 using MnO2

The selective oxidation of alcohol-d₁ to prepare aldehyde-d₁ was newly developed by means of NaBD₄ reduction/activated MnO₂ oxidation. Various aldehyde-d₁ derivatives including aromatic and unsaturated aldehyde-d₁ can be prepared with a high deuterium incorporation ratio (up to 98% D). Halogens (chloride, bromide, and iodide), alkene, alkyne, ester, nitro, and cyano groups in the substrates are tolerated under the mild conditions.

A facile method for deutrium incorporation into aldehydes by mild reduction of NaBD₄ of aldehydes and MnO₂ oxidation (98% D) is disclosed.

Ligand-Directed Conformational Control over Porphyrinic Zirconium Metal-Organic Frameworks for Size-Selective Catalysis

Zirconium-based metal-organic frameworks (Zr-MOFs) have aroused enormous interest owing to their superior stability, flexible structures, and intriguing functions. Precise control over their crystalline structures, including topological structures, porosity, composition, and conformation, constitutes an important challenge to realize the tailor-made functionalization. In this work, we developed a new Zr-MOF (PCN-625) with a csq topological net, which is similar to that of the well-known PCN-222 and NU-1000. However, the significant difference lies in the conformation of porphyrin rings, which are vertical to the pore surfaces rather than in parallel. The resulting PCN-625 exhibits two types of one-dimensional channels with concrete diameters of 2.03 and 0.43 nm. Furthermore, the vertical porphyrins together with shrunken pore sizes could limit the accessibility of substrates to active centers in the framework. On the basis of the structural characteristics, PCN-625(Fe) can be utilized as an efficient heterogeneous catalyst for the size-selective [4 + 2] hetero-Diels-Alder cycloaddition reaction. Due to its high chemical stability, this catalyst can be repeatedly used over six times. This work demonstrates that Zr-MOFs can serve as tailor-made scaffolds with enhanced flexibility for target-oriented functions.

Carbon dioxide-promoted palladium-catalyzed dehydration of primary allylic alcohols: access to substituted 1,3-dienes

A carbon dioxide promoted dehydration reaction of primary allylic alcohols gives synthetically important substituted 1,3-dienes in good yields under milder conditions compared with the reaction using a heterogeneous catalyst or carbon monoxide.

Theoretical investigation of the effect of hydrogen bonding on the stereoselectivity of the Diels–Alder reaction

Here, we report the computational prediction of high exo selectivities in a series of Diels–Alder reactions with H-bonding interaction.