Density functional and coupled cluster calculations show that facile thermal configuration change in [12]annulene occurs via a twist-coupled bond-shifting mechanism. The transition state for this process is highly aromatic with Möbius topology. At the CCSD(T)/cc-pVDZ//BH&HLYP/6-311+G** level, the isomerization of tri-trans-[12]annulene 1a (CTCTCT) to its di-trans isomer 2 (CCCTCT) via such a mechanism has a barrier of 18.0 kcal/mol, in good agreement with earlier experiments. Two other aromatic Möbius bond-shifting transition states were located that result in configuration change for other [12]annulene conformers. This mechanism contrasts sharply with diradical configuration change for acyclic polyenes and with planar bond-shifting mechanisms generally assumed for annulenes. This constitutes evidence that neutral Möbius aromatic annulenes play a role in the dynamic processes of neutral [4n]annulenes.

We report density-functional and coupled-cluster calculations on conformation change and degenerate bond shifting in [10]annulene isomers 1-5. At the CCSD(T)/cc-pVDZ//CCSD/6-31G level, conversion of the twist (1) to the heart (2) has a barrier of 10.1 kcal/mol, compared to Ea = 16.2 kcal/mol for degenerate "two-twist" bond shifting in 1. Pseudorotation in the all-cis boat isomer (3) proceeds with a negligible barrier. The naphthalene-like isomer 4 has a 3.9 kcal/mol barrier to degenerate bond shifting. The azulene-like isomer 5 is the only species for which the nature of the bond-equalized form (5-eq) depends on the method. At the CCSD(T)/cc-pVDZ//CCSD/6-31G level, 5-eq is 1.2 kcal/mol more stable than the bond-alternating form 5-alt. Conversion of 5-eq to 4 has a barrier of 12.6 kcal/mol. Despite being significantly nonplanar, both 5-eq and the transition state for bond shifting in 4 are highly aromatic based on magnetic susceptibility exaltations. On the basis of a detailed consideration of these mechanisms and barriers, we can now, with greater confidence, rule out 4 and 5 as candidates to explain the NMR spectra observed by Masamune. Our results support Masamune's original assignments for both isolated isomers.

Dynamics calculations are described for carbocation rearrangements involving product-forming pathways with post-transition-state bifurcations. We show that noncovalent interactions with associated benzene rings (a simple model of aromatic amino acid side chains) can switch inherent dynamical tendencies for competing modes of disrotation, establishing that meaningful changes in dynamically controlled product selectivity can be achieved with few weak noncovalent interactions.

The lifetimes of secondary carbocations proposed to occur along reaction coordinates for terpene-forming carbocation rearrangements were estimated from direct dynamics simulations using density functional theory. Classical secondary carbocations supported by bicyclo[2.2.1] and bicyclo[2.2.2] frameworks have distinct characteristic lifetimes, ca. 40 and 90 fs, respectively. The fusion of additional rings to these frameworks was found to have little effect on these lifetimes, despite altering the potential energy surfaces for rearrangement. Inherent dynamical tendencies of secondary carbocations are shown to be manipulable by alkylation and enforced intermolecular interactions.

No abstract

Although the aromatic aza-Claisen rearrangement is a general strategy for accessing substituted aromatic amines, there are no highly enantioselective examples of this process. We report the first Brønsted acid catalyzed enantioselective indole aza-Claisen rearrangement for the synthesis of chiral 3-amino-2-substituted indoles. We present evidence for an arene CH-O interaction as a source of activation and stereoinduction, which is an unprecedented phenomenon in enantioselective Brønsted acid catalysis. The products of this reaction can be transformed into 3-aminooxindoles, which are prevalent in many biologically active small molecules.

Density functional and ab initio methods have been used to study the mechanisms for key dynamic processes of the experimentally known S4-symmetric [16]annulene (1a). Using BH&HLYP/6-311+G** and B3LYP/6-311+G**, we located two viable stepwise pathways with computed energy barriers (Ea = 8-10 kcal/mol) for conformational automerization of 1a, in agreement with experimental data. The transition states connecting these conformational minima have Möbius topology and serve as starting points for non-degenerate pi-bond shifting (configuration change) via Möbius aromatic transition states. The key transition state, TS1-2, that connects the two isomers of [16]annulene (CTCTCTCT, 1 --> CTCTTCTT, 2) has an energy, relative to the S4 isomer, that ranged from 6.9 kcal/mol (B3LYP/6-311+G**) to 16.7 kcal/mol (BH&HLYP/6-311+G**), bracketing the experimental barrier. At our best level of theory, CCSD(T)/cc-pVDZ(est), this barrier is 13.7 kcal/mol. Several other Möbius bond-shifting transition states, as well as Möbius topology conformational minima, were found with BH&HLYP energies within 22 kcal/mol of 1a, indicating that many possibilities exist for facile thermal configuration change in [16]annulene. This bond-shifting mechanism and the corresponding low barriers contrast sharply with those observed for cis/trans isomerization in acyclic polyenes, which occurs via singlet diradical transition states. All Möbius bond-shifting transition states located in [16]- and [12]annulene were found to have RHF --> UHF instabilities with the BH&HLYP method but not with B3LYP. This result appears to be an artifact of the BH&HLYP method. These findings support the idea that facile thermal configuration change in [4n]annulenes can be accounted for by mechanisms involving twist-coupled bond shifting.

An introduction to the application of quantum chemical dynamics calculations to mechanistic problems in the field of terpene biosynthesis is provided. A bare bones introduction to the fundamentals of chemical dynamics is followed by a brief account of previous applications to terpene-forming carbocation reactions, a discussion of questions in this field that dynamics calculations may help answer, and a description of current problems to which dynamics calculations are being applied.