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.

Electro/Ni Dual-Catalyzed Decarboxylative C(sp3)-C(sp2) Cross-Coupling Reactions of Carboxylates and Aryl Bromide

Paired redox-neutral electrolysis offers an attractive green platform for organic synthesis by avoiding sacrificial oxidants and reductants. Carboxylates are non-toxic, stable, inexpensive, and widely available, making them ideal nucleophiles for C-C cross-coupling reactions. Here, we report the electro/Ni dual-catalyzed redox-neutral decarboxylative C(sp3)-C(sp2) cross-coupling reactions of pristine carboxylates with aryl bromides. At a cathode, a NiII(Ar)(Br) intermediate is formed through the activation of Ar-Br bond by a NiI-bipyridine catalyst and subsequent reduction. At an anode, the carboxylates, including amino acid, benzyl carboxylic acid, and 2-phenoxy propionic acid, undergo oxidative decarboxylation to form carbon-based free radicals. The combination of NiII(Ar)(Br) intermediate and carbon radical results in the formation of C(sp3)-C(sp2) cross-coupling products. The adaptation of this electrosynthesis method to flow synthesis and valuable molecule synthesis was demonstrated. The reaction mechanism was systematically studied through electrochemical voltammetry and density functional theory (DFT) computational studies. The relationships between the electrochemical properties of carboxylates and the reaction selectivity were revealed. The electro/Ni dual-catalyzed cross-coupling reactions described herein expand the chemical space of paired electrochemical C(sp3)-C(sp2) cross-coupling and represent a promising method for the construction of the C(sp3)-C(sp2) bonds because of the ubiquitous carboxylate nucleophiles and the innate scalability and flexibility of electrochemical flow-synthesis technology.

Intramolecular Heteroatom and Styryl Diels-Alder Reactions, Asymmetric Cycloadditions of Chiral 3-Phenylallyl Maleic Esters

Polycyclic aryl naphthalene and tetralin dihydro arylnaphthalene lactone lignans possess anticancer and antibiotic activity. Related furo[3,4-c]pyranones, typified by the sequester-terpenoid isobolivianine, show similar antiproliferative bioactivity. Efficient syntheses of compounds featuring these polycyclic cores have proven challenging due to low yields and poor stereoselectivity. We report the synthesis of chiral cinnamyl but-2-enanoates and 3,3-diphenylallyl-but-2-enoates 1 as new Diels-Alder substrates. These compounds undergo [4 + 2]-cycloadditions to give furo[3,4-c]pyranones 2 in good yield (70%) and diastereoselectivity (7:1), together with naphthyl 3 and dihydronaphthyl tetralins 4 as minor products. Molecular structures and stereochemistries of the major products were verified using X-ray diffraction. Density functional theory calculations revealed that the cycloaddition process involves a bispericyclic/ambimodal process where there is a single transition state that leads to both intramolecular styryl Diels-Alder (ISDA) 3, 4 and intramolecular hetero Diels-Alder (IHDA) cycloadducts 2. With the elevated temperature conditions after cycloaddition, the resulting ISDA cycloadduct either undergoes [3,3]-sigmatropic rearrangement to the more stable major IHDA product or aromatization leading to the phenyltetralin.

Coordination-induced O-H/N-H bond weakening by a redox non-innocent, aluminum-containing radical

Several renewable energy schemes aim to use the chemical bonds in abundant molecules like water and ammonia as energy reservoirs. Because the O-H and N-H bonds are quite strong (>100 kcal/mol), it is necessary to identify substances that dramatically weaken these bonds to facilitate proton-coupled electron transfer processes required for energy conversion. Usually this is accomplished through coordination-induced bond weakening by redox-active metals. However, coordination-induced bond weakening is difficult with earth's most abundant metal, aluminum, because of its redox inertness under mild conditions. Here, we report a system that uses aluminum with a redox non-innocent ligand to achieve significant levels of coordination-induced bond weakening of O-H and N-H bonds. The multisite proton-coupled electron transfer manifold described here points to redox non-innocent ligands as a design element to open coordination-induced bond weakening chemistry to more elements in the periodic table.

Machine Learning Models for Predicting Zirconocene Properties and Barriers

Zr metallocenes have significant potential to be highly tunable polyethylene catalysts through modification of the aromatic ligand framework. Here we report the development of multiple machine learning models using a large library (>700 systems) of DFT-calculated zirconocene properties and barriers for ethylene polymerization. We show that very accurate machine learning models are possible for HOMO-LUMO gaps of precatalysts but the performance significantly depends on the machine learning algorithm and type of featurization, such as fingerprints, Coulomb matrices, smooth overlap of atomic positions, or persistence images. Surprisingly, the description of the bonding hapticity, the number of direct connections between Zr and the ligand aromatic carbons, only has a moderate influence on the performance of most models. Despite robust models for HOMO-LUMO gaps, these types of machine learning models based on structure connectivity type features perform poorly in predicting ethylene migratory insertion barrier heights. Therefore, we developed several relatively robust and accurate machine learning models for barrier heights that are based on quantum-chemical descriptors (QCDs). The quantitative accuracy of these models depends on which potential energy surface structure QCDs were harvested from. This revealed a Hammett-type principle to naturally emerge showing that QCDs from the π-coordination complexes provide much better descriptions of the transition states than other potential-energy structures. Feature importance analysis of the QCDs provides several fundamental principles that influence zirconocene catalyst reactivity.

Dynamical Origin of Rebound versus Dissociation Selectivity during Fe-Oxo-Mediated C-H Functionalization Reactions

The mechanism of catalytic C-H functionalization of alkanes by Fe-oxo complexes is often suggested to involve a hydrogen atom transfer (HAT) step with the formation of a radical-pair intermediate followed by diverging pathways for radical rebound, dissociation, or desaturation. Recently, we showed that in some Fe-oxo reactions, the radical pair is a nonstatistical-type intermediate and dynamic effects control rebound versus dissociation pathway selectivity. However, the effect of the solvent cage on the stability and lifetime of the radical-pair intermediate has never been analyzed. Moreover, because of the extreme complexity of motion that occurs during dynamics trajectories, the underlying physical origin of pathway selectivity has not yet been determined. For the reaction between [(TQA_Cl)FeIVO]+ and cyclohexane, here, we report explicit solvent trajectories and machine learning analysis on transition-state sampled features (e.g., vibrational, velocity, and geometric) that identified the transferring hydrogen atom kinetic energy as the most important factor controlling rebound versus nonrebound dynamics trajectories, which provides an explanation for our previously proposed dynamic matching effect in fast rebound trajectories that bypass the radical-pair intermediate. Manual control of the reaction trajectories confirmed the importance of this feature and provides a mechanism to enhance or diminish selectivity for the rebound pathway. This led to a general catalyst design principle and proof-of-principle catalyst design that showcases how to control rebound versus dissociation reaction pathway selectivity.

ReaLigands: A Ligand Library Cultivated from Experiment and Intended for Molecular Computational Catalyst Design

Computational catalyst design requires identification of a metal and ligand that together result in the desired reaction reactivity and/or selectivity. A major impediment to translating computational designs to experiments is evaluating ligands that are likely to be synthesized. Here, we provide a solution to this impediment with our ReaLigands library that contains >30,000 monodentate, bidentate (didentate), tridentate, and larger ligands cultivated by dismantling experimentally reported crystal structures. Individual ligands from mononuclear crystal structures were identified using a modified depth-first search algorithm and charge was assigned using a machine learning model based on quantum-chemical calculated features. In the library, ligands are sorted based on direct ligand-to-metal atomic connections and on denticity. Representative principal component analysis (PCA) and uniform manifold approximation and projection (UMAP) analyses were used to analyze several tridentate ligand categories, which revealed both the diversity of ligands and connections between ligand categories. We also demonstrated the utility of this library by implementing it with our building and optimization tools, which resulted in the very rapid generation of barriers for 750 bidentate ligands for Rh-hydride ethylene migratory insertion.

Mechanistic studies of Ni-catalyzed electrochemical homo-coupling reactions of aryl halides

Ni-catalyzed electrochemical arylation is an attractive, emerging approach for molecular construction as it uses air-stable Ni catalysts and efficiently proceeds at room temperature. However, the homo-coupling of aryl halide substrates is one of the major side reactions. Herein, extensive experimental and computational studies were conducted to examine the mechanism of Ni-catalyzed electrochemical homo-coupling of aryl halides. The results indicate that an unstable NiII(Ar)Br intermediate formed through oxidative addition of the cathodically generated NiI species with aryl bromide and a consecutive chemical reduction step. For electron-rich aryl halides, homo-coupling reaction efficiency is limited by the oxidative addition step, which can be improved by negatively shifting the redox potential of the Ni-catalyst. DFT computational studies suggest a NiIII(Ar)Br2/NiII(Ar)Br ligand exchange pathway for the formation of a high-valent NiIII(Ar)2Br intermediate for reductive elimination and production of the biaryl product. This work reveals the reaction mechanism of Ni-catalyzed electrochemical homo-coupling of aryl halides, which may provide valuable information for developing cross-coupling reactions with high selectivity.

Dynamic-dependent selectivity in a bisphosphine iron spin crossover C-H insertion/π-coordination reaction

Reaction pathway selectivity is generally controlled by competitive transition states. Organometallic reactions are complicated by the possibility that electronic spin state changes rather than transition states can control the relative rates of pathways, which can be modeled as minimum energy crossing points (MECPs). Here we show that in the reaction between bisphosphine Fe and ethylene involving spin state crossover (singlet and triplet spin states) that neither transition states nor MECPs model pathway selectivity consistent with experiment. Instead, single spin state and mixed spin state quasiclassical trajectories demonstrate nonstatistical intermediates and that C-H insertion versus π-coordination pathway selectivity is determined by the dynamic motion during reactive collisions. This example of dynamic-dependent product outcome provides a new selectivity model for organometallic reactions with spin crossover.

Directing-Group-Free Arene C(sp2)-H Amination Using Bulky Aminium Radicals and DFT Analysis of Regioselectivity

A hydroxylamine-derived electrophilic aminating reagent produces a transient and bulky aminium radical intermediate upon in situ activation by either TMSOTf or TFA and a subsequent electron transfer from an iron(II) catalyst. Density functional theory calculations were used to examine the regioselectivity of arene C-H amination reactions on diversely substituted arenes. The calculations suggest a simple charge-controlled regioselectivity model that enables prediction of the major C(sp2)-H amination product.

Exploring Cuneanes as Potential Benzene Isosteres and Energetic Materials: Scope and Mechanistic Investigations into Regioselective Rearrangements from Cubanes

Cuneane is a strained hydrocarbon that can be accessed via metal-catalyzed isomerization of cubane. The carbon atoms of cuneane define a polyhedron of the C2v point group with six faces─two triangular, two quadrilateral, and two pentagonal. The rigidity, strain, and unique exit vectors of the cuneane skeleton make it a potential scaffold of interest for the synthesis of functional small molecules and materials. However, the limited previous synthetic efforts toward cuneanes have focused on monosubstituted or redundantly substituted systems such as permethylated, perfluorinated, and bis(hydroxymethylated) cuneanes. Such compounds, particularly rotationally symmetric redundantly substituted cuneanes, have limited potential as building blocks for the synthesis of complex molecules. Reliable, predictable, and selective syntheses of polysubstituted cuneanes bearing more complex substitution patterns would facilitate the study of this ring system in myriad applications. Herein, we report the regioselective, AgI-catalyzed isomerization of asymmetrically 1,4-disubstituted cubanes to cuneanes. In-depth DFT calculations provide a charge-controlled regioselectivity model, and direct dynamics simulations indicate that the nonclassical carbocation invoked is short-lived and dynamic effects augment the charge model.

Understanding Formation and Roles of NiII Aryl Amido and NiIII Aryl Amido Intermediates in Ni-Catalyzed Electrochemical Aryl Amination Reactions

Ni-catalyzed electrochemical aryl amination (e-amination) is an attractive, emerging approach to building C-N bonds. Here, we report in-depth experimental and computational studies that examined the mechanism of Ni-catalyzed e-amination reactions. Key NiII-amine dibromide and NiII aryl amido intermediates were chemically synthesized and characterized. The combination of experiments and DFT calculations suggest (1) there is coordination of an amine to the NiII catalyst before the cathodic reduction and oxidative addition steps, (2) a stable NiII aryl amido intermediate is produced from the cathodic half-reaction, a critical step in controlling the selectivity between cross-coupling and undesired homo-coupling reaction pathways, (3) the diazabicycloundecene additive shifts the aryl halide oxidative addition mechanism from a NiI-based pathway to a Ni0-based pathway, and (4) redox-active bromide in the supporting electrolyte functions as a redox mediator to promote the oxidation of the stable NiII aryl amido intermediate to a NiIII aryl amido intermediate. Subsequently, the NiIII aryl amido intermediate undergoes facile reductive elimination to provide a C-N cross-coupling product at room temperature. Overall, our results provide new fundamental understandings about this e-amination reaction and guidance for further development of other Ni-catalyzed electrosynthetic reactions such as C-C and C-O cross-couplings.

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.

Density functional theory and CCSD(T) evaluation of ionization potentials, redox potentials, and bond energies related to zirconocene polymerization catalysts

Quantum-mechanical-based computational design of molecular catalysts requires accurate and fast electronic structure calculations to determine and predict properties of transition-metal complexes. For Zr-based molecular complexes related to polyethylene catalysis, previous evaluation of density functional theory (DFT) and wavefunction methods only examined oxides and halides or select reaction barrier heights. In this work, we evaluate the performance of DFT against experimental redox potentials and bond dissociation enthalpies (BDEs) for zirconocene complexes directly relevant to ethylene polymerization catalysis. We also examined the ability of DFT to compute the fourth atomic ionization potential of zirconium and the effect the basis set selection has on the ionization potential computed with CCSD(T). Generally, the atomic ionization potential and redox potentials are very well reproduced by DFT, but we discovered relatively large deviations of DFT-calculated BDEs compared to experiment. However, evaluation of BDEs with CCSD(T) suggests that experimental values should be revisited, and our CCSD(T) values should be taken as most accurate.

Microwave- and Thermally Promoted Iminyl Radical Cyclizations: A Versatile Method for the Synthesis of Functionalized Pyrrolines

A detailed study of iminyl radical cyclizations of O-aryloximes tethered to alkenes is reported. The reactions can be triggered by either microwave irradiation or conventional heating in an oil bath. A variety of radical traps can be employed, enabling C-C, C-N, C-O, C-S, or C-X bond formation and producing a diverse array of functionalized pyrrolines. Substrates containing an allylic sulfide furnish terminal alkenes by a tandem cyclization-thiyl radical β-elimination pathway. Cyclizations of hydroxylated substrates exhibit moderate diastereoselectivity that in some cases can partially be attributed to intramolecular hydrogen bonding. Computational studies suggested a possible role for thermodynamics in controlling the stereochemistry of cyclizations. The reaction temperature can be lowered from 120 to 100 °C by employing O-(p-tert-butylphenyl)oximes instead of O-phenyloximes as substrates, and these second-generation iminyl radical precursors can be used in a one-pot oxime ether formation-cyclization that is promoted by conventional heating. The functionalized pyrrolines obtained from these reactions can be conveniently transformed in several different ways.

Experimental Demonstration and Density Functional Theory Mechanistic Analysis of Arene C-H Bond Oxidation and Product Protection by Osmium Tetroxide in a Strongly Basic/Nucleophilic Solvent

Reactions that result in the oxy-functionalization of sp² C-H bonds to give phenols are relatively rare. Here we report experiments and density functional theory (DFT) calculations that demonstrate selective C-H bond hydroxylation of nitroarenes to their corresponding mono-phenoxide as the exclusive product using OsO₄ in a highly basic solvent mixture of water, hydroxide, and pyridine. DFT calculations using a mixed explicit/continuum solvent approach indicate that there is likely a mixture of OsO₄-hydroxide/pyridine ground-state structures that have competitive reactivity and that the mechanism involves the nucleophilic addition of an anionic metal-oxo species to the arene followed by a hydride transfer process that is different from the standard [3 + 2] mechanism often invoked for the OsO₄ oxidation of σ and π bonds. This work demonstrates the utility of using a strongly basic solvent for C-H bond oxidation reactions as this effectively converts any reactive phenolic product into the corresponding phenoxide, which is protected and essentially inert to further oxidation by the nucleophilic metal-oxo species.

Demonstration of High-Throughput Building Block and Composition Analysis of Metal-Organic Frameworks

Metal-organic frameworks (MOFs) are composed of inorganic metal-containing nodes and organic linker groups and are promising porous materials for a wide range of applications. More than 90 000 different MOFs have been synthesized with different inorganic nodes, organic linkers, and node-linker connectivity patterns. While databases have been created to catalog this enormous number of structures, they generally do not provide functionality to easily search, sort, and understand MOFs based on composition and building blocks. Because structure-property relationships are critical to identify, here we outline our new program MOFseek and demonstrate that it can perform high-throughput structure and composition analyses of MOF structures. This program enables the fast analysis of tens of thousands of MOFs in databases based on the local chemical environment. We demonstrate the unique capabilities of MOFseek by analyzing the CoRE MOF database of structures.

Cu(II) carboxylate arene C─H functionalization: Tuning for nonradical pathways

We report carbon-hydrogen acetoxylation of nondirected arenes benzene and toluene, as well as related functionalization with pivalate and 2-ethylhexanoate ester groups, using simple copper(II) [Cu(II)] salts with over 80% yield. By changing the ratio of benzene and Cu(II) salts, 2.4% conversion of benzene can be reached. Combined experimental and computational studies results indicate that the arene carbon-hydrogen functionalization likely occurs by a nonradical Cu(II)-mediated organometallic pathway. The Cu(II) salts used in the reaction can be isolated, recycled, and reused with little change in reactivity. In addition, the Cu(II) salts can be regenerated in situ using oxygen and, after the removal of the generated water, the arene carbon-hydrogen acetoxylation and related esterification reactions can be continued, which leads to a process that enables recycling of Cu(II).

Aza-Quasi-Favorskii Reaction: Construction of Highly Substituted Aziridines through a Concerted Multibond Rearrangement Process

A new molecular rearrangement, the aza-Quasi-Favorskii rearrangement, has been developed for the construction of highly substituted aziridines. Electron-deficient O-sulfonyl oximes react readily with α,α-disubstituted acetophenone-derived enolates to furnish highly substituted aziridines via this unprecedented domino process. In-depth computational studies reveal an asynchronous yet concerted nitrenoid-type rearrangement pathway.

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.

Theory and Experiment Demonstrate that Sb(V)-Promoted Methane C-H Activation and Functionalization Outcompete Superacid Protonolysis in Sulfuric Acid

Sb(V) in strong Brønsted acid solvents is traditionally assumed to react with light alkanes through superacid protonolysis, which results in carbocation intermediates, H₂, and carbon oligomerization. In contrast to this general assumption, our density functional theory (DFT) calculations revealed an accessible barrier for C-H activation between methane and Sb(V) in sulfuric acid that could potentially outcompete superacid protonolysis. This prompted us to experimentally examine this reaction in sulfuric acid with oleum, which has never been reported because of presumed superacid reactivity. Reaction of methane at 180 °C for 3 h resulted in very high yields of methyl bisulfate without significant overoxidation. Our DFT calculations show that a C-H activation and Sb-Me bond functionalization mechanism to give methyl bisulfate outcompetes methane protonolysis and many other possible reaction mechanisms, such as electron transfer, proton-coupled electron transfer, and hydride abstraction. Our DFT calculations also explain experimental hydrogen-deuterium exchange studies and the absence of methane carbo-functionalization/oligomerization products. Overall, this work demonstrates that in very strong Brønsted acid solvent, Sb(V) can induce innersphere reaction mechanisms akin to transition metals and outcompete superacid reactivity.

Direct Dynamics Trajectories Reveal Nonstatistical Coordination Intermediates and Demonstrate that σ and π-Coordination Are Not Required for Rhenium(I)-Mediated Ethylene C-H Activation

The C-H activation reaction between Cp(PMe₃)₂Re and ethylene results in kinetic selectivity for the Re-vinyl hydride I over the thermodynamically more stable Cp(PMe₃)₂Re(η²-ethylene) π-complex II. While transition-state and variational transition-state structures were located for individual pathways leading to I and II, DFT and CCSD(T) energies predict a large kinetic selectivity of 10²-10⁴, which is incompatible with the experimental 10:1 ratio. DFT direct quasiclassical trajectories revealed that the transition states do not provide a qualitatively correct reaction mechanism or a quantitatively correct selectivity due to a nonstatistical σ-CH coordination intermediate that precedes the transition states for C-H activation and π coordination. Using metadynamics and quasiclassical direct dynamics, we show that trajectories for the reaction between Cp(PMe₃)₂Re and ethylene result in direct formation of either the Re-vinyl hydride I or the π-complex II. Trajectories leading to the Re-vinyl hydride skip σ-coordination and do not require π-coordination. Consistent with experiments, trajectory selectivity provides a relatively small kinetic selectivity for the Re-vinyl hydride.

Machine learning classification of disrotatory IRC and conrotatory non-IRC trajectory motion for cyclopropyl radical ring opening

Quasiclassical trajectory analysis is now a standard tool to analyze non-minimum energy pathway motion of organic reactions. However, due to the large amount of information associated with trajectories, quantitative analysis of the dynamic origin of reaction selectivity is complex. For the electrocyclic ring opening of cyclopropyl radical, more than 4000 trajectories were run showing that allyl radicals are formed through a mixture of disrotatory intrinsic reaction coordinate (IRC) motion as well as conrotatory non-IRC motion. Geometric, vibrational mode, and atomic velocity transition-state features from these trajectories were used for supervised machine learning analysis with classification algorithms. Accuracy >80% with a random forest model enabled quantitative and qualitative assessment of transition-state trajectory features controlling disrotatory IRC versus conrotatory non-IRC motion. This analysis revealed that there are two key vibrational modes where their directional combination provides prediction of IRC versus non-IRC motion.

Experiments and Direct Dynamics Simulations That Probe η2-Arene/Aryl Hydride Equilibria of Tungsten Benzene Complexes

Key steps in the functionalization of an unactivated arene often involve its dihaptocoordination by a transition metal followed by insertion into the C-H bond. However, rarely are the η²-arene and aryl hydride species in measurable equilibrium. In this study, the benzene/phenyl hydride equilibrium is explored for the {WTp(NO)(PBu₃)} (Bu = n-butyl; Tp = trispyrazoylborate) system as a function of temperature, solvent, ancillary ligand, and arene substituent. Both face-flip and ring-walk isomerizations are identified through spin-saturation exchange measurements, which both appear to operate through scission of a C-H bond. The effect of either an electron-donating or electron-withdrawing substituent is to increase the stability of both arene and aryl hydride isomers. Crystal structures, electrochemical measurements, and extensive NMR data further support these findings. Static density functional theory calculations of the benzene-to-phenyl hydride landscape suggest a single linear sequence for this transformation involving a sigma complex and oxidative cleavage transition state. Static DFT calculations also identified an η²-coordinated benzene complex in which the arene is held more loosely than in the ground state, primarily through dispersion forces. Although a single reaction pathway was identified by static calculations, quasiclassical direct dynamics simulations identified a network of several reaction pathways connecting the η²-benzene and phenyl hydride isomers, due to the relatively flat energy landscape.

Quantum-mechanical transition-state model combined with machine learning provides catalyst design features for selective Cr olefin oligomerization

The use of data science tools to provide the emergence of non-trivial chemical features for catalyst design is an important goal in catalysis science. Additionally, there is currently no general strategy for computational homogeneous, molecular catalyst design. Here, we report the unique combination of an experimentally verified DFT-transition-state model with a random forest machine learning model in a campaign to design new molecular Cr phosphine imine (Cr(P,N)) catalysts for selective ethylene oligomerization, specifically to increase 1-octene selectivity. This involved the calculation of 1-hexene : 1-octene transition-state selectivity for 105 (P,N) ligands and the harvesting of 14 descriptors, which were then used to build a random forest regression model. This model showed the emergence of several key design features, such as Cr-N distance, Cr-α distance, and Cr distance out of pocket, which were then used to rapidly design a new generation of Cr(P,N) catalyst ligands that are predicted to give >95% selectivity for 1-octene.

H/D exchange under mild conditions in arenes and unactivated alkanes with C6D6 and D2O using rigid, electron-rich iridium PCP pincer complexes

The synthesis and characterization of an iridium polyhydride complex (Ir-H4) supported by an electron-rich PCP framework is described. This complex readily loses molecular hydrogen allowing for rapid room temperature hydrogen isotope exchange (HIE) at the hydridic positions and the α-C-H site of the ligand with deuterated solvents such as benzene-d6, toluene-d8 and THF-d8. The removal of 1-2 equivalents of molecular H2 forms unsaturated iridium carbene trihydride (Ir-H3) or monohydride (Ir-H) compounds that are able to create further unsaturation by reversibly transferring a hydride to the ligand carbene carbon. These species are highly active hydrogen isotope exchange (HIE) catalysts using C6D6 or D2O as deuterium sources for the deuteration of a variety of substrates. By modifying conditions to influence the Ir-Hn speciation, deuteration levels can range from near exhaustive to selective only for sterically accessible sites. Preparative level deuterations of select substrates were performed allowing for procurement of >95% deuterated compounds in excellent isolated yields; the catalyst can be regenerated by treatment of residues with H2 and is still active for further reactions.

Machine Learning Analysis of Direct Dynamics Trajectory Outcomes for Thermal Deazetization of 2,3-Diazabicyclo[2.2.1]hept-2-ene

Experimentally, the thermal gas-phase deazetization of 2,3-diazabicyclo[2.2.1]hept-2-ene (1) results in the loss of N₂ and the formation of bicyclo products 3 (exo) and 4 (endo) in a nonstatistical ratio, with preference for the exo product. Here, we report unrestricted M06-2X quasiclassical trajectories initialized from the concerted N₂ ejection transition state that were able to replicate the experimental preference to form 3. We found that the 3:4 ratio results from the relative amounts of very fast (ballistic) exotype trajectories versus trajectories that lead to the 1,3-diradical intermediate 2. These quasiclassical trajectories provided a set of transition-state vibrational, velocity, momenta, and geometric features for the machine learning analysis. A selection of popular supervised classification algorithms (e.g., random forest) provided poor prediction of trajectory outcomes based on only transition-state vibrational quanta and energy features. However, these machine learning models provided more accurate predictions using atomic velocities and atomic positions, attaining ∼70% accuracy using initial conditions and between 85 and 95% accuracy at later reaction time steps. This increased accuracy allowed the feature importance analysis to reveal that, at the later-time analysis, the methylene bridge out-of-plane bending is correlated with trajectory outcomes for the formation of either the exo product or toward the diradical intermediate. Possible reasons for the struggle of machine learning algorithms to classify trajectories based on transition-state features is the heavily overlapping feature values, the finite but very large possible vibrational mode combinations, and the possibility of chaos as trajectories propagate. We examined this chaos by comparing a set of nearly identical trajectories that differed by only a very small scaling of the kinetic energies resulting from the transition-state reaction coordinate.

Synthetic investigations of low-coordinate (N-phosphino-amidinate) nickel chemistry: agostic alkyl complexes and benzene insertion into Ni-H

Treatment of (PN)NiX (X = NHdipp or OtBu; PN = N-phosphinoamidinate ligand) with Me2PhSiH in benzene solvent afforded the crystallographically characterized, antifacial-coordinated, dinuclear species 1, the formation of which corresponds to the hitherto unknown net Ni-H addition of two equivalents of the putative (PN)NiH intermediate across C[double bond, length as m-dash]C units within a single benzene molecule. Computational analysis supports the view of 1 as being comprised of two cationic (PN)NiII fragments ligated by a substituted butadiene dianion μ2-η3:η3-C6H82- bridging group. Also described is the formation and characterization of three-coordinate (PN)Ni(alkyl) complexes stabilized by β-agostic (alkyl = Et, 2; n-Bu, 3; n-hexyl, 4) or γ-agostic (alkyl = neopentyl, 5) interactions, and our efforts to employ 2 and 3 as synthons for the generation of (PN)NiHvia β-hydride elimination. Notably, compound 5 represents both the first crystallographically characterized three-coordinate Ni-alkyl complex featuring a heterobidentate ligation, and the first neutral γ-agostic NiII-alkyl complex.

Semiempirical Peak Fitting Guided by ab Initio Calculations of X-ray Photoelectron Spectroscopy Narrow Scans of Chemisorbed, Fluorinated Silanes

Here, we address the issue of finding correct CF₂/CF₃ area ratios from X-ray photoelectron spectroscopy (XPS) C 1s narrow scans of materials containing -CH₂CH₂(CF₂)_nCF₃ (n = 0, 1, 2, ...) moieties. For this work, we modified silicon wafers with four different fluorosilanes. The smallest had a trifluoropropyl (n = 0) moiety, followed by nonafluorohexyl (n = 3), tridecafluoro (n = 5), and finally, heptadecafluoro (n = 7) moieties. Monolayer deposition of the fluorosilanes was confirmed by spectroscopic ellipsometry, wetting, and XPS. Analysis of the trifluoropropyl (n = 0) surface and a sample of polytetrafluoroethylene provided pure-component XPS spectra for -CF₃ and -(CF₂)_n- moieties, respectively. Initial XPS C 1s peak fitting, which follows the literature precedent, was not entirely adequate. To address this issue, six different fitting approaches with increasing complexity and/or input from the Hartree-Fock theory (HF) were considered. Ultimately, we show that by combining HF results with empirical analyses, we obtain more accurate CF₂/CF₃ area ratios while maintaining high-quality fits.

Mechanistic molecular motion of transition-metal mediated β-hydrogen transfer: quasiclassical trajectories reveal dynamically ballistic, dynamically unrelaxed, two step, and concerted mechanisms

Quasiclassical direct dynamics reveal new dynamical mechanisms for metal-alkyl to ethylene β-hydrogen transfer.

Supermetal: SbF5-mediated methane oxidation occurs by C-H activation and isobutane oxidation occurs by hydride transfer

Sb^VF₅ is generally assumed to oxidize methane through a methanium-to-methyl cation mechanism. However, experimentally no H₂ is observed, and the mechanism of methane oxidation has remained unsolved for several decades. To solve this problem, density functional theory calculations with multiple chemical models (mononuclear and dinuclear) were used to examine methane oxidation by Sb^VF₅ in the presence of CO leading to the methyl acylium cation ([CH₃CO]⁺). While there is a low barrier for methane protonation by [Sb^VF₆]^-[H]⁺ (the combination of Sb^VF₅ and HF) to give the [Sb^VF₅]^-[CH₅]⁺ ion pair, H₂ dissociation is a relatively high energy process, even with CO assistance, and so this protonation pathway is reversible. While Sb-mediated hydride transfer has a reasonable barrier, the C-H activation/σ-bond metathesis mechanism with the formation of an Sb^V-Me intermediate is lower in energy. This pathway leads to the acylium cation by functionalization of the Sb^V-Me intermediate with CO and is consistent with no observation of H₂. Because this C-H activation/metal-alkyl functionalization pathway is higher in energy than methane protonation, it is also consistent with the experimentally observed methane hydrogen-to-deuterium exchange. This is the first time that evidence is presented demonstrating that Sb^VF₅ acts beyond a Bronsted superacid and involves C-H activation with an organometallic intermediate. In contrast to methane, due to the much lower carbocation hydride affinity, isobutane significantly favors hydride transfer to give the tert-butyl carbocation with concomitant Sb^V to Sb^III reduction. In this mechanism, the resulting highly acidic Sb^V-H intermediate provides a route to H₂ through protonation of isobutane, which is consistent with experiments and resolves the longstanding enigma of different experimental results for methane versus isobutane.

Synthesis of Structurally Diverse 3-, 4-, 5-, and 6-Membered Heterocycles from Diisopropyl Iminomalonates and Soft C-Nucleophiles

Herein, we present a general synthetic strategy for the preparation of 3-, 4-, 5-, and 6-membered heterocyclic unnatural amino acid derivatives by exploiting facile Mannich-type reactions between readily available N-alkyl- and N-aryl-substituted diisopropyl iminomalonates and a wide range of soft anionic C-nucleophiles without using any catalyst or additive. Fully substituted aziridines were obtained in a single step when enolates of α-bromo esters were employed as nucleophiles. Enantiomerically enriched azetidines, γ-lactones, and tetrahydroquinolines were obtained via a two-step catalytic asymmetric reduction and cyclization sequence from ketone enolate-derived adducts. Finally, highly substituted γ-lactams were prepared in one pot from adducts obtained using acetonitrile-derived carbanions. Overall, this work clearly demonstrates the utility of iminomalonates as highly versatile building blocks for the practical and scalable synthesis of structurally diverse heterocycles.