Relationships between Product Ratios in Ambimodal Pericyclic Reactions and Bond Lengths in Transition Structures

Bispericyclic Ambimodal Dimerization of Pentafulvene: The Origin of Asynchronicity and Kinetic Selectivity of the Endo Transition State

The propensity of fulvenes to undergo dimerization has long been known, although the in-depth mechanism and electronic behavior during dimerization are still elusive. Herein, we made an attempt to gain insights into the reactivity of pentafulvene for Diels-Alder (DA) and [6 + 4]-cycloadditions via conventional and ambimodal routes. The result emphasizes that pentafulvene dimerization preferentially proceeds through a unique bifurcation mechanism where two DA pathways merge together to produce two degenerate [4 + 2]-cycloadducts from a single TS. Despite the [6 + 4]-cycloadduct being thermodynamically preferred, [4 + 2]-cycloaddition reactions are kinetically driven. Singlet biradicaloid is involved in through-space 6e- delocalization as a secondary orbital interaction that originates asynchronicity and stabilizes the bispericyclic transition state (TS). The transformation of various actively participating intrinsic bonding orbitals (IBOs) unambiguously forecasts the formation of multiple products from a single TS and rationalizes the mechanism of ambimodal reactions that are rather difficult to probe with other analyses. The changes in active IBOs clearly distinguish the conventional reactions from bifurcation reactions and can be employed to characterize and confirm the ambimodal mechanism. This report gains a crucial theoretical insight into the mechanism of bifurcation, the origin of asynchronicity, and electronic behavior in ambimodal TS, which will certainly be of enormous value for future studies.

Periselectivity and ambimodal transition states in cycloadditions of tetrachloro-o-benzoquinone with 6,6-dimethylfulvene

The reaction mechanism of cycloadditions of tetrachloro-o-benzoquinone with 6,6-dimethylfulvene were systematically investigated with density functional theory calculations. It was found that conditional primary interactions stabilize the ambimodal transition states in the endo pathways. Ambimodal transition states lead to [6 + 4]/[4 + 2] adducts or [4 + 2]/[2 + 4] adducts, which interconvert through 3,3-sigmatropic shift reactions. The substituent effects on periselectivity were also investigated.

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.

Dynamically or Kinetically Controlled? Computational Study of the Mechanisms of Electrophilic Aminoalkenylation of Heteroaromatics with Keteniminium Ions

Quantum chemical calculations and molecular dynamics simulations were applied to study the electrophilic aminoalkenylation of heteroaromatics with keniminium ions. Post-transition state bifurcation (PTSB) was found in the electrophilic addition step for the aminoalkenylation of pyrroles and indoles, and the selectivity for these reactions was dynamically controlled. However, the aminoalkenylation of furan was kinetically controlled because no apparent PTSB was found in the electrophilic addition step. The substituents on the keteniminium ions can also affect the dynamic results for the aminoalkenylations to pyrroles: the C2-aminoalkenylated product is much more favored over the C3-aminoalkenylated product for keteniminium ions with electron-donating substituents, while the product ratio (C2 product/C3 product) decreased when stronger electron-withdrawing substituents were applied.

Binomial Uncertainty in Molecular Dynamics-Based Reactions Analysis

Molecular Dynamics-based reaction analysis is an indispensable tool for studying processes defying the transition-state theory (TST), where the product ratios do not follow energies of transition states. The main class of such processes is ambimodal reactions, which have a post-transition-state bifurcation, so that several products form via a single transition state. Multiple runs of molecular dynamics allow one to sample the space of possibilities and ultimately predict the product ratio without relying on TST; however, no techniques for estimating the reliability of the prediction were proposed so far. Here we show that dynamics runs follow the same rules as die rolls, which paves a simple way for estimating their uncertainty and, accordingly, the number of runs necessary to achieve the required accuracy. Remarkably, we find that the majority of such studies carried out in the last 5 years use far too few runs, so that the product ratios predicted in them can be off by >50% in more than 50% of cases.

Post-Transition-State Dynamic Effects in the Transmetalation of Pd(II)-F to Pd(II)-CF3

The observation of post-transition-state dynamic effects in the context of metal-based transformation is rare. To date, there has been no reported case of a dynamic effect for the widely employed class of palladium-mediated coupling reactions. We performed an experimental and computational study of the trifluoromethylation of Pd(II)F, which is a key step in the Pd(0)/Pd(II)-catalyzed trifluoromethylation of aryl halides or acid fluorides. Our experiments show that the cis/trans speciation of the formed Pd(II)CF3 is highly solvent- and transmetalation reagent-dependent. We employed GFN2-xTB- and B3LYP-D3-based molecular dynamics trajectory calculations (with and without explicit solvation) along with high-level QM calculations and found that depending on the medium, different transmetalation mechanisms appear to be operative. A statistically representative number of Born-Oppenheimer molecular dynamics (MD) simulations suggest that in benzene, a difluorocarbene is generated in the transmetalation with R3SiCF3, which subsequently recombines with the Pd via two distinct pathways, leading to either the cis- or trans-Pd(II)CF3. Conversely, GFN2-xTB simulations in MeCN suggest that in polar/coordinating solvents an ion-pair mechanism is dominant. A CF3 anion is initially liberated and then rebinds with the Pd(II) cation to give a cis- or trans-Pd(II). In both scenarios, a single transmetalation transition state gives rise to both cis- and trans-species directly, owing to bifurcation after the transition state. The potential subsequent cis- to trans isomerization of the Pd(II)CF3 was also studied and found to be strongly inhibited by free phosphine, which in turn was experimentally identified to be liberated through displacement by a polar/coordinating solvent from the cis-Pd(II)CF3 complex. The simulations also revealed how the variation of the Pd-coordination sphere results in divergent product selectivities.

Reaction dynamics as the missing puzzle piece: the origin of selectivity in oxazaborolidinium ion-catalysed reactions

The selectivity in a group of oxazaborolidinium ion-catalysed reactions between aldehyde and diazo compounds cannot be explained using transition state theory. VRAI-selectivity, developed to predict the outcome of dynamically controlled reactions, can account for both the chemo- and the stereo-selectivity in these reactions, which are controlled by reaction dynamics. Subtle modifications to the substrate or catalyst substituents alter the potential energy surface, leading to changes in predominant reaction pathways and altering the barriers to the major product when reaction dynamics are considered. In addition, this study suggests an explanation for the mysterious inversion of enantioselectivity resulting from the inclusion of an orthoiPrO group in the catalyst.

Origins of Periselectivity and Regioselectivity in Ambimodal Tripericyclic [8+6]/[6+4]/[4+2] Intramolecular Cycloadditions of a Heptafulvenyl-Fulvene

Quantum mechanical calculations and molecular dynamics simulations have elucidated the reaction mechanism for intramolecular cycloadditions of a heptafulvenyl-fulvene tethered by a trimethylene chain. Prior experiments by Liu and Houk reported the formation of only an endo-[8+6] cycloadduct at 185 °C. Liu et al. later reported an exo-[4+2] Diels-Alder cycloadduct as the major product at 140 °C (Tetrahedron, 1999, 55, 9171). Cycloadditions involve Diels-Alder and an ambimodal intramolecular tripericyclic [8+6]/[6+4]/[4+2] cycloaddition. The mechanistic details explain the experimental reports of temperature dependence on the periselectivity of intramolecular cycloadditions. Additional calculations with multireference-based methods CASSCF and NEVPT2 highlight the artifacts of DFT methods and single-reference wavefunction-based CCSD(T) in the description of complete potential energy surface involving various cycloadditions of the heptafulvenyl-fulvene.

Four-Dimensional-Spacetime Atomistic Artificial Intelligence Models

We demonstrate that AI can learn atomistic systems in the four-dimensional (4D) spacetime. For this, we introduce the 4D-spacetime GICnet model, which for the given initial conditions (nuclear positions and velocities at time zero) can predict nuclear positions and velocities as a continuous function of time up to the distant future. Such models of molecules can be unrolled in the time dimension to yield long-time high-resolution molecular dynamics trajectories with high efficiency and accuracy. 4D-spacetime models can make predictions for different times in any order and do not need a stepwise evaluation of forces and integration of the equations of motions at discretized time steps, which is a major advance over traditional, cost-inefficient molecular dynamics. These models can be used to speed up dynamics, simulate vibrational spectra, and obtain deeper insight into nuclear motions, as we demonstrate for a series of organic molecules.

Barrier Heights for Diels-Alder Transition States Leading to Pentacyclic Adducts: A Benchmark Study of Crowded, Strained Transition States of Large Molecules

Theoretical characterization of reactions of complex molecules depends on providing consistent accuracy for the relative energies of intermediates and transition states. Here we employ the DLPNO-CCSD(T) method with core-valence correlation, large basis sets, and extrapolation to the CBS limit to provide benchmark values for Diels-Alder transition states leading to competitive strained pentacyclic adducts. We then used those benchmarks to test a diverse set of wave function and density functional methods for the absolute and relative barrier heights of these transition states. Our results show that only a few of the tested density functionals can predict the absolute barrier heights satisfactorily, although relative barrier heights are more accurate. The most accurate functionals tested are ωB97M-V, M11plus, ωB97X-V, PBE-D3(0), M11, and MN15 with MUDs from best estimates less than 3.0 kcal. These findings can guide selection of density functionals for future studies of crowded, strained transition states of large molecules.

Dynamics study of the post-transition-state-bifurcation process of the (HCOOH)H+ → CO + H3O+/HCO+ + H2O dissociation: application of machine-learning techniques

The process of protonated formic acid dissociating from the transition state was studied using ring-polymer molecular dynamics (RPMD), classical MD, and quasi-classical trajectory (QCT) simulations. Temperature had a strong influence on the branching fractions for the HCO⁺ + H₂O and CO + H₃O⁺ dissociation channels. The RPMD and classical MD simulations showed similar behavior, but the QCT dynamics were significantly different owing to the excess energies in the quasi-classical trajectories. Machine-learning analysis identified several key features in the phase information of the vibrational motions at the transition state. We found that the initial configuration and momentum of a hydrogen atom connected to a carbon atom and the shrinking coordinate of the CO bond at the transition state play a role in the dynamics of HCO⁺ + H₂O production.

Accelerated Entropic Path Sampling with a Bidirectional Generative Adversarial Network

The role of entropy in mediating the dynamic outcomes of chemical reactions remains largely unknown. To evaluate the change of entropy along post-transition state paths, we have previously developed entropic path sampling that computes configurational entropy from an ensemble of reaction trajectories. However, one major caveat of this approach lies in its high computational demand: about 2000 trajectories are needed to converge the computation of an entropic profile. Here, by leveraging a deep generative model, we developed an accelerated entropic path sampling approach that evaluates entropic profiles using merely a few hundred reaction dynamic trajectories. The new method, called bidirectional generative adversarial network-entropic path sampling, can enhance the estimation of probability density functions of molecular configurations by generating pseudo-molecular configurations that are statistically indistinguishable from the true data. The method was established using cyclopentadiene dimerization, in which we reproduced the reference entropic profiles (derived from 2480 trajectories) using merely 124 trajectories. The method was further benchmarked using three reactions with symmetric post-transition-state bifurcation, including endo-butadiene dimerization, 5-fluoro-1,3-cyclopentadiene dimerization, and 5-methyl-1,3-cyclopentadiene dimerization. The results indicate the existence of a "hidden entropic intermediate", which is a dynamic species that binds to a local entropic maximum where no free energy minimum is formed.

Design of C1-symmetric tridentate ligands for enantioselective dearomative [3 + 2] annulation of indoles with aminocyclopropanes

Chiral polycyclic indolines are widely present in natural products and have become the focus of extensive synthetic efforts. Here, we show the catalytic asymmetric dearomative [3 + 2] annulation of indoles with donor-acceptor aminocyclopropanes to construct tricyclic indolines. Key to the success of the reaction is the rational design of C₁-symmetric bifunctional tridentate imidazoline-pyrroloimidazolone pyridine ligand. Under 5 mol% of Ni(OTf)₂-ligand complex, diverse tricyclic indolines containing cyclopentamine moieties are obtained in good chemoselectivities, high diastereoselectivities, and excellent enantioselectivities. An unusual cis-configuration ligand is superior to the trans-configuration ligand and the corresponding C₂-symmetric tridentate nitrogen ligands in the annulation reaction. Mechanistic studies by control experiments and density functional theory calculations reveal a dual activation manner, where Ni(II) complex activates the aminocyclopropane via coordination with the geminal diester, and imidazolidine NH forms a H-bond with the succinimide moiety.

Molecular Understanding and Practical In Silico Catalyst Design in Computational Organocatalysis and Phase Transfer Catalysis-Challenges and Opportunities

Through the lens of organocatalysis and phase transfer catalysis, we will examine the key components to calculate or predict catalysis-performance metrics, such as turnover frequency and measurement of stereoselectivity, via computational chemistry. The state-of-the-art tools available to calculate potential energy and, consequently, free energy, together with their caveats, will be discussed via examples from the literature. Through various examples from organocatalysis and phase transfer catalysis, we will highlight the challenges related to the mechanism, transition state theory, and solvation involved in translating calculated barriers to the turnover frequency or a metric of stereoselectivity. Examples in the literature that validated their theoretical models will be showcased. Lastly, the relevance and opportunity afforded by machine learning will be discussed.

Computational Design of a Tetrapericyclic Cycloaddition and the Nature of Potential Energy Surfaces with Multiple Bifurcations

An ambimodal transition state (TS) that leads to formation of four different pericyclic reaction products ([4 + 6]-, [2 + 8]-, [8 + 2]-, and [6 + 4]-cycloadducts) without any intervening minima has been designed and explored with DFT computations and quasiclassical molecular dynamics. Direct dynamics simulations propagated from the ambimodal TS show the evolution of trajectories to give the four cycloadducts. The topography of the PES is a key factor in product selectivity. A good correlation is observed between geometrical resemblance of the products to the ambimodal TS (measured by the RMSD) and the ratio of products formed in the dynamics simulations.

Symmetric Post-Transition State Bifurcation Reactions with Berry Pseudomagnetic Fields

We investigate how the Berry force (i.e., the pseudomagnetic force operating on nuclei as induced by electronic degeneracy and spin-orbit coupling (SOC)) might modify a post-transition state bifurcation (PTSB) reaction path and affect product selectivity for situations when multiple products share the same transition state. To estimate the magnitude of this effect, Langevin dynamics are performed on a model system with a valley-ridge inflection (VRI) point in the presence of a magnetic field (that mimics the Berry curvature). We also develop an analytic model for such selectivity that depends on key parameters such as the surface topology, the magnitude of the Berry force, and the nuclear friction. Within this dynamical model, static electronic structure calculations (at the level of generalized Hartree-Fock with spin-orbit coupling (GHF+SOC) theory) suggest that electronic spin induced Berry force effects may indeed lead to noticeable changes in methoxy radical isomerization.

Prediction of Kinetic Product Ratios: Investigation of a Dynamically Controlled Case

Of the various factors influencing kinetically controlled product ratios, the role of nonstatistical dynamics is arguably the least well understood. In this paper, reactions were chosen in which dynamics played a dominant role in product selection, by design. Specifically, the reactions studied were the ring openings of cyclopropylidene to allene and tetramethylcyclopropylidene to tetramethylallene (2,4-dimethylpenta-2,3-diene). Both reactions have intrinsic reaction coordinates that bifurcate symmetrically, leading to products that are enantiomeric once the atoms are uniquely labeled. The question addressed in the study was whether the outcomes─that is, which product well on the potential energy surface was selected─could be predicted from their initial conditions for individual trajectories in quasiclassical dynamics simulations. Hybrid potentials were developed based on cooperative interaction between molecular mechanics and artificial neural networks, trained against data from electronic structure calculations. These potentials allowed simulations of both gas-phase and condensed-phase reactions. The outcome was that, for both reactions, prediction of initial selection of product wells could be made with >95% success from initial conditions of the trajectories in the gas phase. However, when trajectories were run for longer, looking for "final" products for each trajectory, the predictability dropped off dramatically. In the gas-phase simulations, this drop off was caused by trajectories hopping between product wells on the potential energy surface. That behavior could be suppressed in condensed phases, but then new uncertainty was introduced because the intermolecular interactions between solute and bath, necessary to permit intermolecular energy transfer and cooling of the hot initial products, often led to perturbations of the initial directions of trajectories on the potential energy surface. It would consequently appear that a general ability to predict outcomes for reactions in which nonstatistical dynamics dominate remains a challenge even in the age of sophisticated machine-learning capabilities.

Nickel-Catalyzed Homologation of Vinylidene Difluoride (CH2═CF2): Selective β-F vs β-H Elimination

Hydrofluoroolefins (HFOs) constitute the newest generation of fluorocarbon refrigerants and foam-blowing agents due to their reduced global warming potential vs their saturated analogues. To identify new synthetic routes to HFOs, we show that reactions of bulky Ni(0) phosphine and -NHC complexes with vinylidene difluoride (VF2) afford μ-fluoro-1,1,3-trifluorobut-3-enyl Ni complexes. Moreover, addition of triisopropylsilane allows for reductive elimination of the reduced product─2,4,4-trifluoro-1-butene─demonstrating the Ni-catalyzed hydrodefluorodimerization of VF2. Accompanying DFT calculations identify the T-shaped nickelacyclopentane intermediate that spontaneously undergoes selective intramolecular β-F (vs β-H) elimination.

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

Rh-catalyzed C-H insertion reactions to form β-lactones suffer from post-transition state bifurcations, with the same transition states leading to ketones and ketenes via fragmentation in addition to β-lactones. In such a circumstance, traditional transition state theory cannot predict product selectivity, so we employed ab initio molecular dynamics simulations to do so and provide a framework for rationalizing the origins of said selectivity. Weak interactions between the catalyst and substrate were studied using energy decomposition and noncovalent interaction analyses, which unmasked an important role of the 2-bromophenyl substituent that has been used in multiple β-lactone-forming C-H insertion reactions. Small and large catalysts were shown to behave differently, with the latter providing a means of overcoming dynamically preferred fragmentation by lowering the barrier for the recombination of the product fragments in the grip of the large catalyst active site cavity.

Collective Synthesis of Illudalane Sesquiterpenes via Cascade Inverse Electron Demand (4 + 2) Cycloadditions of Thiophene S,S-Dioxides

Thiophene S,S-dioxides are underutilized tools for the de novo construction of benzene rings in organic synthesis. We report a collective synthesis of nine illudalane sesquiterpenes using bicyclic thiophene S,S-dioxides as generalized precursors to the indane core of the natural products. Exploiting furans as unusual dienophiles in this inverse electron demand Diels-Alder cascade, this concise and convergent approach enables the synthesis of these targets in as little as five steps. Theoretical studies rationalize the reactivity of thiophene S,S-dioxides with both electron-poor and electron-rich dienophiles and reveal reaction pathways involving either nonpolar pericyclic or bifurcating ambimodal cycloadditions. Overall, this work demonstrates the wider potential of thiophene S,S-dioxides as convenient and flexible precursors to polysubstituted arenes.

Kinetic Analysis of a Reaction Path Network Including Ambimodal Transition States: A Case Study of an Intramolecular Diels-Alder Reaction

This study proposes a methodology for the kinetic analysis of a reaction path network including ambimodal transition states (TSs), through which an ensemble of trajectories bifurcates to multiple minima in a phenomenon called dynamical bifurcation. The proposed methodology consists of three techniques: an automated reaction path search to construct a reaction path network including ambimodal TSs, an ab initio molecular dynamics simulation to evaluate the branching ratio, and the definition of rate constants incorporating this ratio. Applying the procedure to a Diels-Alder reaction, it was found that the inclusion of dynamical bifurcations is necessary to explain the experimental reaction yield of a byproduct. In addition, it was verified that the products take 10¹³ s to reach thermal equilibrium and that the experimental selectivity is determined by the dynamical bifurcations.

Computational investigation of cycloadditions between cyclopentadiene and tropone-3,4-dimethylester

Thermally promoted cycloaddition reactions of tropone-3,4-dimethylester and cyclopentadiene have been investigated using density functional theory calculations at the M06-2X level and the CBS-QB3 method.

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.

Influence of Water and Enzyme on the Post-Transition State Bifurcation of NgnD-Catalyzed Ambimodal [6+4]/[4+2] Cycloaddition

The enzyme NgnD catalyzes an ambimodal cycloaddition that bifurcates to [6+4]- and [4+2]-adducts. Both products have been isolated in experiments, but it remains unknown how enzyme and water influence the bifurcation selectivity at the femtosecond time scale. Here, we study the impact of water and enzyme on the post-transition state bifurcation of NgnD-catalyzed [6+4]/[4+2] cycloaddition by integrating quantum mechanics/molecular mechanics quasiclassical dynamics simulations and biochemical assays. The ratio of [6+4]/[4+2] products significantly differs in the gas phase, water, and enzyme. Biochemical assays were employed to validate computational predictions. The study informs how water and enzyme affect the bifurcation selectivity through perturbation of the reaction dynamics in the femtosecond time scale, revealing the fundamental roles of condensed media in dynamically controlling the chemical selectivity for biosynthetic reactions.

BH9, a New Comprehensive Benchmark Data Set for Barrier Heights and Reaction Energies: Assessment of Density Functional Approximations and Basis Set Incompleteness Potentials

The calculation of accurate reaction energies and barrier heights is essential in computational studies of reaction mechanisms and thermochemistry. To assess methods regarding their ability to predict these two properties, high-quality benchmark sets are required that comprise a reasonably large and diverse set of organic reactions. Due to the time-consuming nature of both locating transition states and computing accurate reference energies for reactions involving large molecules, previous benchmark sets have been limited in scope, the number of reactions considered, and the size of the reactant and product molecules. Recent advances in coupled-cluster theory, in particular local correlation methods like DLPNO-CCSD(T), now allow the calculation of reaction energies and barrier heights for relatively large systems. In this work, we present a comprehensive and diverse benchmark set of barrier heights and reaction energies based on DLPNO-CCSD(T)/CBS called BH9. BH9 comprises 449 chemical reactions belonging to nine types common in organic chemistry and biochemistry. We examine the accuracy of DLPNO-CCSD(T) vis-a-vis canonical CCSD(T) for a subset of BH9 and conclude that, although there is a penalty in using the DLPNO approximation, the reference data are accurate enough to serve as a benchmark for density functional theory (DFT) methods. We then present two applications of the BH9 set. First, we examine the performance of several density functional approximations commonly used in thermochemical and mechanistic studies. Second, we assess our basis set incompleteness potentials regarding their ability to mitigate basis set incompleteness errors. The number of data points, the diversity of the reactions considered, and the relatively large size of the reactant molecules make BH9 the most comprehensive thermochemical benchmark set to date and a useful tool for the development and assessment of computational methods.

Bifurcating reactions: distribution of products from energy distribution in a shared reactive mode

Bifurcating reactions yield two different products emerging from one single transition state and are therefore archetypal examples of reactions that cannot be described within the framework of the traditional Eyring's transition state theory (TST). With the growing number and importance of these reactions in organic and biosynthetic chemistry, there is also an increasing demand for a theoretical tool that would allow for the accurate quantification of reaction outcome at low cost. Here, we introduce such an approach that fulfils these criteria, by evaluating bifurcation selectivity through the energy distribution within the reactive mode of the key transition state. The presented method yields an excellent agreement with experimentally reported product ratios and predicts the correct selectivity for 89% of nearly 50 various cases, covering pericyclic reactions, rearrangements, fragmentations and metal-catalyzed processes as well as a series of trifurcating reactions. With 71% of product ratios determined within the error of less than 20%, we also found that the methodology outperforms three other tested protocols introduced recently in the literature. Given its predictive power, the procedure makes reaction design feasible even in the presence of complex non-TST chemical steps.

Reactive Mode Composition Factor (RMCF) analysis is a powerful tool to forecast the product distribution of bifurcating reactions through analysis of the kinetic energy distribution within the first transition state traversed by the reacting system.

Mechanisms and Dynamics of Synthetic and Biosynthetic Formation of Delitschiapyrones: Solvent Control of Ambimodal Periselectivity

The mechanism and dynamics for the formation of the delitschiapyrone family of natural products are studied by density functional theory (DFT) calculations and quasiclassical molecular dynamics simulations with DFT and xTB. In the uncatalyzed reaction, delitschiapyrones A and B are formed by Diels-Alder reactions through a single transition state and a post-transition state bifurcation that favors formation of delitschiapyrone B. In water and most likely in the enzyme, the acidic hydroxyquinone ionizes, and the resulting conjugate base undergoes cycloaddition preferentially to delitschiapyrone A. We demonstrate a new type of biosynthetic transformation and variable selectivity from a (4 + 2)/(4 + 3) ambimodal transition state.

Pathway Bifurcations in the Activation of Allylic Halides by Palladium and Their Influence on the Dynamics of η1 and η3 Allyl Intermediates

Transition-metal-catalyzed allylic substitution often exhibits complex product selectivity patterns, which have been primarily attributed to π ↔ σ ↔ π isomerization of the η¹ and η³ allyl intermediates. Product selectivity may be even further complicated if η¹- and η³-allyls share a single transition state (TS), leading to their formation resulting in a post-transition-state bifurcation (PTSB). In this work, density functional theory calculations using ab initio molecular dynamics (AIMD) have been carried out that support the presence of a PTSB in Pd-catalyzed allylic halide activation directly influencing product selectivity. The AIMD results initiated from the TS predict the η¹-allyl to be favored in the gas phase and a low dielectric (ε < 2.5) for trialkylphosphines, while the selectivity shifts toward the η³-allyl in higher dielectrics. The minimum energy path is also predicted to shift in product preference, consistent with the dynamics predictions. The bifurcation in allylic chloride activation is predicted to largely favor the η³-allyl at any solvent polarity. A PTSB was also discovered to be present in Ni and Pt allylic activation but with less bifurcation. These results offer a unique view into the mechanism of metal-catalyzed allylic substitution.

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.

Computational Exploration of the Mechanism of Critical Steps in the Biomimetic Synthesis of Preuisolactone A, and Discovery of New Ambimodal (5 + 2)/(4 + 2) Cycloadditions

Computational studies with ωB97X-D density functional theory of the mechanisms of the steps in Trauner's biomimetic synthesis of preuisolactone A have elaborated and refined mechanisms of several unique processes. An ambimodal transition state has been identified for the cycloaddition between an o-quinone and a hydroxy-o-quinone; this leads to both (5 + 2) (with H shift) and (4 + 2) cycloaddition products, which can in principle interconvert via α-ketol rearrangements. The origins of periselectivity of this ambimodal cycloaddition have been investigated computationally with molecular dynamics simulations and tested further by an experimental study. In the presence of bicarbonate ions, the deprotonated hydroxy-o-quinone leads to only the (5 + 2) cycloaddition adduct. A new mechanism for a benzilic acid rearrangement resulting in ring contraction is proposed.

Cycloadditions of Cyclopentadiene and Cycloheptatriene with Tropones: All Endo-[6+4] Cycloadditions Are Ambimodal

The cycloadditions of cyclopentadiene and cycloheptatriene with tropone are some of the earliest published examples of [6+4] cycloaddition reactions. We report quantum mechanical studies (ωB97X-D and DLPNO-CCSD(T)) of transition structures and products of these reactions, as well as quasi-classical molecular dynamics simulations of reaction trajectories. The study reveals that these cycloadditions involve ambimodal transition states resulting in a web of products by pericyclic interconversion pathways. Combined with these studies, calculations of simple parent systems and a thorough meta-analysis of literature examples reveal the general concept that all endo-[6+4] cycloadditions are ambimodal.

A Future Perspective on Phototriggered Isomerizations of Transition Metal Sulfoxides and Related Complexes

Photochromic molecules are examples of light-activated bistable molecules. We highlight the design criteria for a class of ruthenium and osmium sulfoxide complexes that undergo phototriggered isomerization of the bound sulfoxide. The mode of action in these complexes is an excited-state isomerization of the sulfoxide from S-bonded to O-bonded. We discuss the basic mechanism for this transformation and highlight specific examples that demonstrate the effectiveness and efficiency of the isomerization. We subsequently discuss future research directions within the field of phototriggered sulfoxide isomerizations on transition metal polypyridine complexes. These efforts involve new synthetic directions, including the choice of metal as well as new ambidentate ligands for isomerization.

Natural polyenic macrolactams and polycyclic derivatives generated by transannular pericyclic reactions: optimized biogenesis challenging chemical synthesis

Covering from 1992 to the end of 2020-11-20.Genetically-encoded polyenic macrolactams, which are constructed by Nature using hybrid polyketide synthase/nonribosomal peptide synthase (PKSs/NRPSs) assembly lines, are part of the large collection of natural products isolated from bacteria. Activation of cryptic (i.e., silent) gene clusters in these microorganisms has more recently allowed to generate and eventually isolate additional members of the family. Having two unsaturated fragments separated by short saturated chains, the primary macrolactam is posited to undergo transannular reactions and further rearrangements thus leading to the generation of a structurally diverse collection of polycyclic (natural) products and oxidized derivatives. The review will cover the challenges that scientists face on the isolation of these unstable compounds from the cultures of the producing microorganisms, their structural characterization, biological activities, optimized biogenetic routes, as well as the skeletal rearrangements of the primary structures of the natural macrolactams derived from pericyclic reactions of the polyenic fragments. The efforts of the synthetic chemists to emulate Nature on the successful generation and structural confirmation of these natural products will also be reported.

Huisgen's 1,3-Dipolar Cycloadditions to Fulvenes Proceed via Ambimodal [6+4]/[4+2] Transition States

Huisgen's 1960 announcement of the concept of 1,3-dipolar cycloadditions was published the year before Alder's study of the reaction of diazomethane and dimethylfulvene. The diazomethane reaction was studied again in 1970 by Houk et al. and shown to give a [6+4] adduct. Padwa's nitrile ylide cycloaddition to dimethylfulvene (1978) gave [6+4] and [4+2] adducts. We performed computational studies of these reactions with density functional theory (DFT) and show that they involve ambimodal [6+4]/[4+2] transition states that can lead to either type of cycloadduct from one transition state. We dedicate this paper to the extraordinary life and humanity of Rolf Huisgen, and to the undying influence of his discoveries on chemistry.

Dynamical Trajectory Study of the Transannular [6+4] and Ambimodal Cycloaddition in the Biosynthesis of Heronamides

We report the dynamical trajectory study of the transannular [6+4] and ambimodal cycloaddition proposed in the biosynthesis of heronamide natural products. The originally proposed bifurcation of the potential energy surface is found to strongly favor the formation of the [6+4] product, both in the gas phase and in an explicit water environment, as evidenced by our trajectory simulations. The detailed information on how the bonds are formed and how water influences the bonding dynamics was analyzed at the femtosecond time scale.

An excited state dynamics driven reaction: wavelength-dependent photoisomerization quantum yields in [Ru(bpy)2(dmso)2]2

We report the excited-state behavior of a structurally simple bis-sulfoxide complex, cis-S,S-[Ru(bpy)2(dmso)2]2+, as investigated by femtosecond pump-probe spectroscopy. The results reveal that a single photon prompts phototriggered isomerization of one or both dmso ligands to yield a mixture of cis-S,O-[Ru(bpy)2(dmso)2]2+ and cis-O,O-[Ru(bpy)2(dmso)2]2+. The quantum yields of isomerization of each product and relative product distribution are dependent upon the excitation wavelength, with longer wavelengths favoring the double isomerization product, cis-O,O-[Ru(bpy)2(dmso)2]2+. Transient absorption measurements on cis-O,O-[Ru(bpy)2(dmso)2]2+ do not reveal an excited-state isomerization pathway to produce either the S,O or S,S isomers. Femtosecond pulse shaping experiments reveal no change in the product distribution. Pump-repump-probe transient absorption spectroscopy of cis-S,S-[Ru(bpy)2(dmso)2]2+ shows that a pump-repump time delay of 3 ps dramatically alters the S,O : O,O product ratio; pump-repump-probe transient absorption spectroscopy of cis-O,O-[Ru(bpy)2(dmso)2]2+ with a time delay of 3 ps uncovers an excited-state isomerization pathway to produce the S,O isomer. In conjunction with low-temperature steady-state emission spectroscopy, these results are interpreted in the context of an excited-state bifurcating pathway, in which the isomerization product distribution is determined not by thermodynamics, but rather as a dynamics driven reaction.

AFIR explorations of transition states of extended unsaturated systems: automatic location of ambimodal transition states

Based on the artificial force induced reaction (AFIR) method, we proposed a procedure to systematically explore ambimodal transition states (TSs) that cause the dynamical bifurcation.

The expanding world of biosynthetic pericyclases: cooperation of experiment and theory for discovery

Covering: 2000 to 2018 Pericyclic reactions are a distinct class of reactions that have wide synthetic utility. Before the recent discoveries described in this review, enzyme-catalyzed pericyclic reactions were not widely known to be involved in biosynthesis. This situation is changing rapidly. We define the scope of pericyclic reactions, give a historical account of their discoveries as biosynthetic reactions, and provide evidence that there are many enzymes in nature that catalyze pericyclic reactions. These enzymes, the "pericyclases," are the subject of this review.