Bouncing off walls - widths of exit channels from shallow minima can dominate selectivity control

Revisiting a classic carbocation - DFT, coupled-cluster, and ab initio molecular dynamics computations on barbaralyl cation formation and rearrangements

Density functional theory computations were used to model the formation and rearrangement of the barbaralyl cation (C9H+ 9). Two highly delocalized minima were located for C9H+ 9, one of C s symmetry and the other of D 3h symmetry, with the former having lower energy. Quantum chemistry-based NMR predictions affirm that the lower energy structure is the best match with experimental spectra. Partial scrambling was found to proceed through a C 2 symmetric transition structure associated with a barrier of only 2.3 kcal mol-1. The full scrambling was found to involve a C 2v symmetric transition structure associated with a 5.0 kcal mol-1 barrier. Ab initio molecular dynamics simulations initiated from the D 3h C9H+ 9 structure revealed its connection to six minima, due to the six-fold symmetry of the potential energy surface. The effects of tunneling and boron substitution on this complex reaction network were also examined.

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.

A 21st Century View of Allowed and Forbidden Electrocyclic Reactions

In 1965, Woodward and Hoffmann proposed a theory to predict the stereochemistry of electrocyclic reactions, which, after expansion and generalization, became known as the Woodward-Hoffmann Rules. Subsequently, Longuet-Higgins and Abrahamson used correlation diagrams to propose that the stereoselectivity of electrocyclizations could be explained by the correlation of reactant and product orbitals with the same symmetry. Immediately thereafter, Hoffmann and Woodward applied correlation diagrams to explain the mechanism of cycloadditions. We describe these discoveries and their evolution. We now report an investigation of various electrocyclic reactions using DFT and CASSCF. We track the frontier molecular orbitals along the intrinsic reaction coordinate and modeled trajectories and examine the correlation between HOMO and LUMO for thermally forbidden systems. We also investigate the electrocyclizations of several highly polarized systems for which the Houk group had predicted that donor-acceptor substitution can induce zwitterionic character, thereby providing low-energy pathways for formally forbidden reactions. We conclude with perspectives on the field of pericyclic reactions, including a refinement as the meaning of Woodward and Hoffmann's "Violations. There are none!" Lastly, we comment on the burgeoning influence of computations on all fields of chemistry.

Metal-Free Arylation of Benzothiophenes at C4 by Activation as their Benzothiophene S-Oxides

Benzothiophenes, activated by oxidation to the corresponding S-oxides, undergo C-H/C-H-type coupling with phenols to give C4 arylation products. While an electron-withdrawing group at C3 of the benzothiophene is important, the process operates without a directing group and a metal catalyst, thus rendering it compatible with sensitive functionalities-e.g. halides and formyl groups. Quantum chemical calculations suggest a formal stepwise mechanism involving heterolytic cleavage of an aryloxysulfur species to give a π-complex of the corresponding benzothiophene and a phenoxonium cation. Subsequent addition of the phenoxonium cation to the C4 position of the benzothiophene is favored over the addition to C3; Fukui functions predict that the major regioisomer is formed at the more electron-rich position between C3 and C4. Varied selective manipulation of the benzothiophene products showcase the synthetic utility of the metal-free arylation process.

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.

Understanding chemistry: from "heuristic (soft) explanations and reasoning by analogy" to "quantum chemistry"

"Soft theories," i.e., "heuristic models based on reasoning by analogy" largely drove chemistry understanding for 150 years or more. But soft theories have their limitations and with the expansion of chemistry in the mid-20th century, more and more inexplicable (by soft theory) experimental results were being obtained. In the past 50 years, quantum chemistry, most often in the guise of applied theoretical chemistry including computational chemistry, has provided (a) the underlying "hard evidence" for many soft theories and (b) the explanations for chemical phenomena that were unavailable by soft theories. In this publication, we define "hard theories" as "theories derived from quantum chemistry." Both soft and hard theories can be qualitative and quantitative, and the "Houk quadrant" is proposed as a helpful categorization tool. Furthermore, the language of soft theories is often used appropriately to describe quantum chemical results. A valid and useful way of doing science is the appropriate use and application of both soft and hard theories along with the best nomenclature available for successful communication of results and ideas.

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.

Portable Models for Entropy Effects on Kinetic Selectivity

Differences in entropies of competing transition states can direct kinetic selectivity. Understanding and modeling such entropy differences at the molecular level is complicated by the fact that entropy is statistical in nature; i.e., it depends on multiple vibrational states of transition structures, the existence of multiple dynamically accessible pathways past these transition structures, and contributions from multiple transition structures differing in conformation/configuration. The difficulties associated with modeling each of these contributors are discussed here, along with possible solutions, all with an eye toward the development of portable qualitative models of use to experimentalists aiming to design reactions that make use of entropy to control kinetic selectivity.

Entropic Path Sampling: Computational Protocol to Evaluate Entropic Profile along a Reaction Path

Fleeting intermediates constitute dynamically stepwise mechanisms. They have been characterized in molecular dynamics trajectories, but whether these intermediates form a free energy minimum to become entropic intermediates remains elusively defined. We developed a computational protocol known as entropic path sampling to evaluate the entropic variation of reacting species along a reaction path based on an ensemble of trajectories. Using cyclopentadiene dimerization as a model reaction, we observed an entropy maximum along the reaction path which originates from an enhanced conformational flexibility as the reacting species enter into a flat energy region. As the reacting species further approach product formation, unfavorable entropic restriction fails to offset the potential energy drop, resulting in no free energy minimum along the post-TS pathway. Our results show that cyclopentadiene dimerization involves an entropy maximum that leads to dynamic intermediates with elongated lifetimes, but the reaction does not involve entropic intermediates.

Dynamic Effects on Migratory Aptitudes in Carbocation Reactions

Carbocation rearrangement reactions are of great significance to synthetic and biosynthetic chemistry. In pursuit of a scale of inherent migratory aptitude that takes into account dynamic effects, both uphill and downhill ab initio molecular dynamics (AIMD) simulations were used to examine competing migration events in a model system designed to remove steric and electronic biases. The results of these simulations were combined with detailed investigations of potential energy surface topography and variational transition state theory calculations to reveal the importance of nonstatistical dynamic effects on migratory aptitude.