C-H Insertion in Dirhodium Tetracarboxylate-Catalyzed Reactions despite Dynamical Tendencies toward Fragmentation: Implications for Reaction Efficiency and Catalyst Design

Site- and enantioselective B-H functionalization of carboranes

Functionalization of carboranes, icosahedral boron-carbon molecular clusters, is of great interest as they have wide applications in medicinal and materials chemistry. Thus, site- and enantioselective synthesis of carboranes requires complete control of the reaction. Herein, we describe the asymmetric Rh(II)-catalyzed insertion reactions of carbenes into cage B-H bond of carboranes. This reaction thereby generates carboranes possessing a carbon-stereocenter adjacent to cage boron of the carborane, in excellent site- and enantioselectivity under mild reaction conditions. The fully computed transition structures of Rh(II)-catalyzed carbene insertion process through density functional theory are reported. These B-H insertion transition structures, in conjunction with topographical proximity surfaces analyses, visually reveal the region between the carborane and the phthalimide ligands responsible for the selectivities of this reaction.

Hexafluoroisopropanol Solvent Effects on Enantioselectivity of Dirhodium Tetracarboxylate-Catalyzed Cyclopropanation

In recent years, additives that modulate both reactivity and selectivity in rhodium-catalyzed reactions of aryldiazoacetates have become increasingly prominent. 1,1,1,3,3,3-Hexafluoroisopropanol (HFIP) has been shown to have a profound effect on rhodium carbene reactivity and selectivity, especially on enabling carbene cyclopropanation in the presence of various nucleophilic poisons. HFIP also has a variable influence on the enantioselectivity of the reactions catalyzed by chiral dirhodium tetracarboxylates, and this study examines the fundamental properties of the rhodium carbene/HFIP system through experimentation, density functional theory (DFT), and molecular dynamics (MD) simulations. These studies revealed that the C4-symmetric bowl-shaped catalysts, which have been previously considered to be relatively rigid, experience far greater flexibility in this hydrogen bonding media, resulting in distortion of the bowl-shaped catalysts. These studies explain why even though a majority of the catalysts have a drop in enantioselectivity in HFIP, some catalysts, such as Rh2(TCPTAD)4, lead to a switch in enantioselectivity, whereas others, such as Rh2(NTTL)4, lead to a considerably enhanced enantioselectivity.

Dirhodium Complexes Heterochiral-at-the-Metal Centers: An Alternative Type of Paddlewheel Catalyst for Asymmetric Synthesis

The prototype of an entirely new class of dirhodium paddlewheel complexes is disclosed, which is chiral at both inequivalent metal centers although it carries just three different achiral μ-bridging equatorial ligands. One of them is a carboxamidate, which ensures selective carbene formation at the Rh[O3N] face of the catalyst; moreover, the -NH group engages the ester carbonyl of the emerging rhodium carbene intermediate derived from α-stannylated α-diazoacetate in interligand hydrogen bonding, which is critically important for controlling the stereochemical course of the ensuing [2 + 1] cycloaddition. The new heterochiral catalyst furnished stannylated cyclopropane derivatives with excellent diastereo- and enantioselectivity; moreover, a cyclopropane with a quaternary trifluoromethylated stereocenter was obtained with an unprecedented level of asymmetric induction.

Synthetically Relevant Post-Transition State Bifurcation Leading to Diradical and Zwitterionic Intermediates: Controlling Nonstatistical Kinetic Selectivity through Solvent Effects

A post-transition state surface intersection (PTSSI) between radical and zwitterionic states that causes a bifurcation in the reaction pathway was discovered through density functional theory calculations on potential energy surfaces and ab initio molecular dynamics simulations of cycloadditions between a bicyclobutane and a triazolinedione (BCB-TAD). It was predicted that changes to the solvent polarity would enable control over the dynamic selectivity in this system; indeed, experimental evidence supported this prediction. This work not only provides new insights into an unusual type of post-transition state bifurcation, but also demonstrates how the nonstatistical dynamic effects that control selectivity for such reactions can be manipulated rationally to increase the yields of synthetically useful reactions.

Molecular Dynamics of the Davies Ambimodal C-H Functionalization/Cope Rearrangement Reaction

The mechanism of the dirhodium-catalyzed combined C-H functionalization/Cope rearrangement (CH/Cope) reaction discovered by the Davies group has been investigated with density functional theory (DFT) calculations and quasi-classical molecular dynamics (MD) simulations. Computations from the Davies group previously showed that there is a post-transition state bifurcation leading to a direct CH reaction and also to the CH/Cope product. While this work was in preparation, the Tantillo group and the Ess group independently reported quantum mechanical and molecular dynamics studies on the dirhodium-tetracarboxylate-catalyzed diazoester CH/Cope and CH insertion reactions with 1,3-cyclohexadiene and 1,4-cyclohexadiene, respectively. The Tantillo group cited "dynamic mismatching" to explain the origins of the low yield of CH/Cope products in some experiments; the Ess group explained the origins of product selectivity from the perspective of TS vibrational modes and their synchronization that occurs at the entropic intermediates. We report quasi-classical trajectories for the reaction of the carbene with 1-methylcyclohexene that afford both the CH/Cope and C-H insertion products. After passing through the transition state that involves mostly hydrogen transfer, momentum drives the reaction trajectories toward the CH/Cope products.

Quantum Chemical Interrogation of Reactions Promoted by Dirhodium Tetracarboxylate Catalysts─Mechanism, Selectivity, and Nonstatistical Dynamic Effects

ConspectusRh2L4 catalysts have risen in popularity in the world of organic synthesis, being used to accomplish a variety of reactions, including C-H insertion and cyclopropanation, and often doing so with high levels of stereocontrol. While the mechanisms and origins of selectivity for such reactions have been examined with computational quantum chemistry for decades, only recently have detailed pictures of the dynamic behavior of reacting Rh2L4-complexed molecules become accessible. Our computational studies on Rh2L4 catalyzed reactions are described here, with a focus on C-H insertion reactions of Rh2L4-carbenes. Several issues complicate the modeling of these reactions, each providing an opportunity for greater understanding and each revealing issues that should be incorporated into future rational design efforts. First, the fundamental mechanism of C-H insertion is discussed. While early quantum chemical studies pointed to transition structures with 3-center [C-H-C] substructures and asynchronous hydride transfer/C-C bond formation, recent examples of reactions with particularly flat potential energy surfaces and even discrete zwitterionic intermediates have been found. These reactions are associated with systems bearing π-donating groups at the site of hydride transfer, allowing for an intermediate with a carbocation substructure at that site to be selectively stabilized. Second, the possible importance of solvent coordination at the Rh atom distal to the carbene is discussed. While effects on reactivity and selectivity were found to be small, they turn out not to be negligible in some cases. Third, it is shown that, in contrast to many other transition metal promoted reactions, many Rh2L4 catalyzed reactions likely involve dissociation of the Rh2L4 catalyst before key chemical steps leading to products. When to expect dissociation is associated with specific features of substrates and the product-forming reactions in question. Often, dissociation precedes transition structures for pericyclic reactions that involve electrons that would otherwise bind to Rh2L4. Finally, the importance of nonstatistical dynamic effects, characterized through ab initio molecular dynamics studies, in some Rh2L4 catalyzed reactions is discussed. These are reactions where transition structures are shown to be followed by flat regions, very shallow minima, and/or pathways that bifurcate, all allowing for trajectories from a single transition state to form multiple different products. The likelihood of encountering such a situation is shown to be associated again with the likelihood of formation of zwitterionic structures along reaction paths, but ones for which pathways to multiple products are expected to be associated with very low or no barriers. The connection between these features and reduced yields of desired products are highlighted, as are the means by which some Rh2L4 catalysts modulate dynamic behavior to produce particular products in high yield.

Trajectory-Based Time-Resolved Mechanism for Benzene Reductive Elimination from Cyclopentadienyl Mo/W Phenyl Hydride Complexes

Calculated potential energy structures and landscapes are very often used to define the sequence of reaction steps in an organometallic reaction mechanism and interpret kinetic isotope effect (KIE) measurements. Underlying most of this structure-to-mechanism translation is the use of statistical rate theories without consideration of atomic/molecular motion. Here we report direct dynamics simulations for an organometallic benzene reductive elimination reaction, where nonstatistical intermediates and dynamic-controlled pathways were identified. Specifically, we report single spin state as well as mixed spin state quasiclassical direct dynamics trajectories in the gas phase and explicit solvent for benzene reductive elimination from Mo and W bridged cyclopentadienyl phenyl hydride complexes ([Me2Si(C5Me4)2]M(H)(Ph), M = Mo and W). Different from the energy landscape mechanistic sequence, the dynamics trajectories revealed that after the benzene C-H bond forming transition state (often called reductive coupling), σ-coordination and π-coordination intermediates are either skipped or circumvented and that there is a direct pathway to forming a spin flipped solvent caged intermediate, which occurs in just a few hundred femtoseconds. Classical molecular dynamics simulations were then used to estimate the lifetime of the caged intermediate, which is between 200 and 400 picoseconds. This indicates that when the η2-π-coordination intermediate is formed, it occurs only after the first formation of the solvent-caged intermediate. This dynamic mechanism intriguingly suggests the possibility that the solvent-caged intermediate rather than a coordination intermediate is responsible (or partially responsible) for the inverse KIE value experimentally measured for W. Additionally, this dynamic mechanism prompted us to calculate the kH/kD KIE value for the C-H bonding forming transition states of Mo and W. Surprisingly, Mo gave a normal value, while W gave an inverse value, albeit small, due to a much later transition state position.

Vibrational synchronization and its reaction pathway influence from an entropic intermediate in a dirhodium catalyzed allylic C-H activation/Cope rearrangement reaction

In reactions with consecutive transition states without an intermediate, and an energy surface bifurcation, atomic motion generally determines product selectivity. Understanding this dynamic-based selectivity can be straightforward if there is extremely fast descent from the first transition state to a product. However, in cases where a nonstatistical roaming/entropic intermediate occurs prior to product formation the motion that influences selectivity can be difficult to identify. Here we report quasiclassical direct dynamics trajectories for the dirhodium catalyzed reaction between styryldiazoacetate and 1,4-cyclohexadiene and prior experiments by Davies showed competitive allylic C-H insertion and Cope products. Trajectories confirmed the proposed energy surface bifurcation and revealed that dirhodium vinylcarbenoid when reacting with 1,4-cyclohexadiene can induce either a dynamically concerted pathway or a dynamically stepwise pathway with a nonstatistical entropic tight ion-pair intermediate. In the dynamically stepwise reaction pathway C-H insertion versus Cope selectivity is highly influenced by whether or not vibrational synchronization occurs in the nonstatistical entropic intermediate. This vibrational synchronization highlights the possible need for an entropic intermediate to have organized transition state-like motion to proceed to a product.

Analogies between photochemical reactions and ground-state post-transition-state bifurcations shed light on dynamical origins of selectivity

Revealing the origins of kinetic selectivity is one of the premier tasks of applied theoretical organic chemistry, and for many reactions, doing so involves comparing competing transition states. For some reactions, however, a single transition state leads directly to multiple products, in which case non-statistical dynamic effects influence selectivity control. The selectivity of photochemical reactions-where crossing between excited-state and ground-state surfaces occurs near ground-state transition structures that interconvert competing products-also should be controlled by the momentum of the reacting molecules as they return to the ground state in addition to the shape of the potential energy surfaces involved. Now, using machine-learning-assisted non-adiabatic molecular dynamics and multiconfiguration pair-density functional theory, these factors are examined for a classic photochemical reaction-the deazetization of 2,3-diazabicyclo[2.2.2]oct-2-ene-for which we demonstrate that momentum dominates the selectivity for hexadiene versus [2.2.2] bicyclohexane products.

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.

Exploring Downhill Bifurcations in [3,3]-Sigmatropic Rearrangement by Finding Transitions from an Uphill Bifurcation to a Downhill Bifurcation

Downhill bifurcation is a phenomenon in which an ensemble of trajectories passing through a transition state (TS), called an ambimodal TS, bifurcates into multiple products. Finding downhill bifurcations for unreported pairs of chemical transformations is essential, because they affect reaction selectivity. Marx et al. reported that perturbations such as applying mechanical stress or changing a substituent cause a transition from an uphill bifurcation to a downhill bifurcation in the ring-opening reaction of cyclopropane derivatives (ChemPhysChem, 2018, 19, 837-847). Investigating the occurrence of this phenomenon in other reactions, especially in pericyclic reactions, is interesting for understanding and controlling the reaction selectivity considering downhill bifurcations. In this study, we proposed a method for finding perturbation-induced downhill bifurcations and applied it to three pericyclic reactions. The transition from an uphill bifurcation to a downhill bifurcation occurred in two of the three pericyclic reactions, one of which was previously unreported. Interestingly, the occurrence of a downhill bifurcation by a perturbation depended on the directions of the intrinsic reaction coordinate paths of the two TSs when they emerged from the reactant minimum. Our method can be applied in mechanistic studies to avoid the risk of overlooking downhill bifurcations.

Modulating Escape Channels of Cycloheptatrienyl Rhodium Carbenes To Form Semibullvalene

We describe the various escape channels available to dirhodium carbene intermediates from cycloheptatrienyl diazo compounds located with density functional theory. An intramolecular cyclopropanation would, in principle, provide a new route to semibullvalenes (SBVs). A detailed exploration of the potential energy surface reveals that methylating carbon-7 suppresses a competing β-hydride migration pathway to heptafulvene products, giving SBV formation a reasonable chance. During our explorations, we additionally discovered unusual spirononatriene, spironorcaradiene, and metal-stabilized 9-barbaralyl cation structures as local minima.

Kinetic, Thermodynamic, and Dynamic Control in Normal vs. Cross [2 + 2] Cycloadditions of Ene-Keteniminium Ions: Computational Understanding, Prediction, and Experimental Verification

Almost all reported intramolecular [2 + 2] reactions of ene-keteniminium ions gave normal [2 + 2] products with a fused bicycle framework, but not cross [2 + 2] products with a bicyclo[3.1.1]heptane skeleton, a highly pursued bioisostere in pharmaceutical chemistry. How to rationalize this and design new cross [2 + 2] reactions? Theoretical studies using density functional theory, high-level ab initio single-point energy calculations, and molecular dynamics showed that this [2 + 2] reaction has all three patterns of regiochemical control: the reaction is controlled either kinetically, thermodynamically, or dynamically. A carbocation model of forming endo and exo carbocations has been proposed to rationalize the reaction outcomes, revealing that the tethers (between alkenes and keteniminium ions), substituents (on the alkenes), and alkene configurations in ene-keteniminium ions play critical roles. These understandings were further used to predict that introducing a substituent in the terminal position of alkene with a trans configuration in ene-keteniminium ions can realize the cross [2 + 2] reaction, which is dynamically controlled for alkyl substituents or kinetically controlled for aryl substituents. These and more other predictions were realized experimentally, and many cross [2 + 2] products with a bicyclo[3.1.1]heptane skeleton can be achieved. Both molecular dynamics and new experiments have also been applied to correct a key but misassigned [2 + 2] product reported in the literature, further supporting the insightful mechanisms reported here.

Direct Dynamics Trajectories Demonstrate Dynamic Matching and Nonstatistical Radical Pair Intermediates during Fe-Oxo-Mediated C-H Functionalization Reactions

The generally proposed mechanism for the reaction between non-heme Fe-oxo complexes and alkane C-H bonds involves a hydrogen atom transfer (HAT) reaction step with a radical pair intermediate that then has competitive radical rebound, dissociation, or desaturation pathways. Here, we report density functional theory-based quasiclassical direct dynamics trajectories that examine post-HAT reaction dynamics. Trajectories revealed that the radical pair intermediate can be a nonstatistical type intermediate without complete internal vibrational redistribution and post-HAT selectivity is generally determined by dynamic effects. Fast rebound trajectories occur through dynamic matching between the rotational motion of the newly formed Fe-OH bond and collision with the alkane radical, and all of this occurs through a nonsynchronous dynamically concerted process that circumvents the radical pair intermediate structure. For radical pair dissociation, trajectories proceeded to the radical pair intermediate for a very brief time, followed by complete dissociation. These trajectories provide a new viewpoint and model to understand the inherent reaction pathway selectivity for non-heme Fe-oxo-mediated C-H functionalization reactions.