Computational Prediction of Excited-State Carbon Tunneling in the Two Steps of Triplet Zimmerman Di-π-Methane Rearrangement

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.

Heavy-Atom Tunneling in Bicyclo[4.1.0]hepta-2,4,6-trienes

Heavy-atom tunneling limits the lifetime and observability of bicyclo[4.1.0]hepta-2,4,6-triene, a key intermediate in the rearrangement of phenylcarbene. Bicyclo[4.1.0]hepta-2,4,6-triene had been proposed as the primary intermediate of the rearrangement of phenylcarbene, but despite many efforts evaded its characterization even in cryogenic matrices. By introducing fluorine substituents into the ortho-positions of the phenyl ring of phenylcarbene, the highly strained cyclopropene 1,5-difluorobicyclo[4.1.0]hepta-2,4,6-triene becomes stable enough to be characterized in argon matrices. However, even at 3 K this cyclopropene is only metastable and rearranges via heavy-atom tunneling to the corresponding cycloheptatetraene. Calculations suggest that fluorination is necessary to slow down the tunneling rearrangement of the bicycloheptatriene. The parent bicycloheptatriene rapidly rearranges via heavy-atom tunneling and therefore cannot be detected under matrix isolation conditions.

Quantum Tunneling in Reactions Modulated by External Electric Fields: Reactivity and Selectivity

Quantum tunneling and external electric fields (EEFs) can promote some reactions. However, the synergetic effect of an EEF on a tunneling-involving reaction and its temperature-dependence is not very clear. In this study, we extensively investigated how EEFs affect three reactions that involve hydrogen- or (ground- and excited-state) carbon-tunneling using reliable DFT, DLPNO-CCSD(T1), and variational transition-state theory methods. Our study revealed that oriented EEFs can significantly reduce the barrier and corresponding barrier width (and vice versa) through more electrostatic stabilization in transition states. These EEF effects enhance the nontunneling and tunneling-involving rates. Such EEF effects also decrease the crossover temperatures and quantum tunneling contribution, albeit with lower and thinner barriers. Moreover, EEFs can modulate and switch on/off the tunneling-driven 1,2-H migration of hydroxycarbenes under cryogenic conditions. Furthermore, our study predicts for the first time that EEF/tunneling synergy can control the chemo- or site-selectivity of one molecule bearing two similar/same reactive sites.

New Insights and Predictions into Complex Homogeneous Reactions Enabled by Computational Chemistry in Synergy with Experiments: Isotopes and Mechanisms

Homogeneous catalysis and biocatalysis have been widely applied in synthetic, medicinal, and energy chemistry as well as synthetic biology. Driven by developments of new computational chemistry methods and better computer hardware, computational chemistry has become an essentially indispensable mechanistic "instrument" to help understand structures and decipher reaction mechanisms in catalysis. In addition, synergy between computational and experimental chemistry deepens our mechanistic understanding, which further promotes the rational design of new catalysts. In this Account, we summarize new or deeper mechanistic insights (including isotope, dispersion, and dynamical effects) into several complex homogeneous reactions from our systematic computational studies along with subsequent experimental studies by different groups. Apart from uncovering new mechanisms in some reactions, a few computational predictions (such as excited-state heavy-atom tunneling, steric-controlled enantioswitching, and a new geminal addition mechanism) based on our mechanistic insights were further verified by ensuing experiments.The Zimmerman group developed a photoinduced triplet di-π-methane rearrangement to form cyclopropane derivatives. Recently, our computational study predicted the first excited-state heavy-atom (carbon) quantum tunneling in one triplet di-π-methane rearrangement, in which the reaction rates and ¹²C/¹³C kinetic isotope effects (KIEs) can be enhanced by quantum tunneling at low temperatures. This unprecedented excited-state heavy-atom tunneling in a photoinduced reaction has recently been verified by an experimental ¹²C/¹³C KIE study by the Singleton group. Such combined computational and experimental studies should open up opportunities to discover more rare excited-state heavy-atom tunneling in other photoinduced reactions. In addition, we found unexpectedly large secondary KIE values in the five-coordinate Fe(III)-catalyzed hetero-Diels-Alder pathway, even with substantial C-C bond formation, due to the non-negligible equilibrium isotope effect (EIE) derived from altered metal coordination. Therefore, these KIE values cannot reliably reflect transition-state structures for the five-coordinate metal pathway. Furthermore, our density functional theory (DFT) quasi-classical molecular dynamics (MD) simulations demonstrated that the coordination mode and/or spin state of the iron metal as well as an electric field can affect the dynamics of this reaction (e.g., the dynamically stepwise process, the entrance/exit reaction channels).Moreover, we unveiled a new reaction mechanism to account for the uncommon Ru(II)-catalyzed geminal-addition semihydrogenation and hydroboration of silyl alkynes. Our proposed key gem-Ru(II)-carbene intermediates derived from double migrations on the same alkyne carbon were verified by crossover experiments. Additionally, our DFT MD simulations suggested that the first hydrogen migration transition-state structures may directly and quickly form the key gem-Ru-carbene structures, thereby "bypassing" the second migration step. Furthermore, our extensive study revealed the origin of the enantioselectivity of the Cu(I)-catalyzed 1,3-dipolar cycloaddition of azomethine ylides with β-substituted alkenyl bicyclic heteroarenes enabled by dual coordination of both substrates. Such mechanistic insights promoted our computational predictions of the enantioselectivity reversal for the corresponding monocyclic heteroarene substrates and the regiospecific addition to the less reactive internal C═C bond of one diene substrate. These predictions were proven by our experimental collaborators. Finally, our mechanistic insights into a few other reactions are also presented. Overall, we hope that these interactive computational and experimental studies enrich our mechanistic understanding and aid in reaction development.

Quantum mechanical tunnelling: the missing term to achieve sub-kJ mol-1 barrier heights

To predict barrier heights at low temperatures, it is not enough to employ highly accurate electronic structure methods. We discuss the influence of quantum tunnelling on the comparison of experimental and theoretical activation parameters (E_a, ΔH^‡, ΔG^‡, or ΔS^‡), since the slope-based experimental techniques to obtain them completely neglect the tunnelling component. The intramolecular degenerate rearrangement of four fluxional molecules (bullvalene, barbaralane, semibullvalene, and norbornadienylidene) were considered, systems that cover the range between fast deep tunneling and small but significant shallow tunnelling correction. The barriers were computed with the composite W3lite-F12 method at the CCSDT(Q)/CBS level, and the tunnelling contribution with small curvature tunnelling. While at room temperature the effect is small (∼1 kJ mol^-1), at low temperatures it can be considerable (in the order of tens of kJ mol^-1 at ∼80 K).

Nunes C, Fausto R. Spin-forbidden heavy-atom tunneling in the ring-closure of triplet cyclopentane-1,3-diyl

The putative spin-forbidden heavy-atom tunneling process for the ring closure of cyclopentane-1,3-diyl at cryogenic temperatures is confirmed with calculations employing the weak-coupling formulation of nonadiabatic transition state theory.

Vibrationally Hot and Cold Triplets. Sensitizer-Dependent Dynamics and Localized Vibrational Promotion of a Di-π-methane Rearrangement

Large intramolecular 13C kinetic isotope effects (KIEs) for the di-π-methane rearrangement of benzobarrelene fit with statistical expectations from heavy-atom tunneling when a low-energy sensitizer is employed, but much lower KIEs are observed with higher-energy sensitizers. These results in combination with trajectory studies suggest that the excess vibrational energy available from triplet energy transfer leads to hot and nonstatistical dynamics in the rearrangement.

Degradation Mechanism of Benzo[a]pyrene Initiated by the OH Radical and 1O2: An Insight from Density Functional Theory Calculations

The degradation mechanism of benzo[a]pyrene (BaP) initiated by ^•OH and ¹O₂ in aqueous solution is investigated by density functional theory calculations. The main degradation products are BaP-1,6-quinone, BaP-3,6-quinone, BaP-4,6-quinone, and BaP-6,12-quinone. ^•OH and HO₂ are the main intermediate radical species. At a low initial concentration of ^•OH, ¹O₂ could be a primary driver for BaP degradation. The degradation mechanism includes six consecutive elementary reactions: (1) ¹O₂ initiation forming BaP-6-OO. (2) 1,3 H-shift (H atom shifts to the OO group) that is promoted by H₂O, forming BaP-6-OOH. (3) BaP-6-OOH decomposes into the ^•OH radical and BaP-6-O. (4) ^•OH addition to BaP-6-O forming BaP-6-O-1(3,4,12)-OH. (5) Extracting the H atom from the carbon with the OH group by ¹O₂. (6) Extracting the H atom from the OH group by HO₂. At a high initial concentration of ^•OH, the ^•OH-initiated and ¹O₂-initiated degradation reactions of BaP are both feasible. The degradation mechanism includes six consecutive elementary reactions: (1) ^•OH initiation forming BaP-6-OH or ¹O₂ initiation forming BaP-6-OO. (2) ¹O₂ addition to BaP-6-OH forming BaP-6-OH-12(1,3,4)-OO or ^•OH addition to BaP-6-OO forming BaP-6-OO-12(1,3,4)-OH. (3) Extracting the H atom from the carbon with the OH group by ¹O₂, forming HO₂. (4) 1,3 H-shift (H-shift from the carbon to the OO group), promoted by H₂O. (5) The loss of the OH radical. (6) Abstracting the H atom from the OH group by HO₂. In this paper, the formation of BaP-4,6-quinone via the BaP degradation is first reported. Water participates in the elementary reaction in which the H atom attached on the aromatic ring shifts to the OO group, serving as a bridge that stabilizes the transition state and transports the proton. A comprehensive investigation explains the degradation mechanism of BaP initiated by ^•OH and ¹O₂ in aqueous solution.

Distannabarrelenes with Three Coordinated SnII Atoms

Crystalline 1,4-distannabarrelene compounds [(ADCAr )3 Sn2 ]SnCl3 (3-Ar) (ADCAr ={ArC(NDipp)2 CC}; Dipp=2,6-iPr2 C6 H3 , Ar=Ph or DMP; DMP=4-Me2 NC6 H4 ) derived from anionic dicarbenes Li(ADCAr ) (2-Ar) (Ar=Ph or DMP) have been reported. The cationic moiety of 3-Ar features a barrelene framework with three coordinated SnII atoms at the 1,4-positions, whereas the anionic unit SnCl3 is formally derived from SnCl2 and chloride ion. The all carbon substituted bis-stannylenes 3-Ar have been characterized by NMR spectroscopy and X-ray diffraction. DFT calculations reveal that the HOMO of 3-Ph (ϵ=-6.40 eV) is mainly the lone-pair orbital at the SnII atoms of the barrelene unit. 3-Ar readily react with sulfur and selenium to afford the mixed-valence SnII /SnIV compounds [(ADCAr )3 SnSn(E)](SnCl6 )0.5 (E=S 4-Ar, Ar=Ph or DMP; E=Se 5-Ph).

Calculations Predict That Heavy-Atom Tunneling Dominates Möbius Bond Shifting in [12]- and [16]Annulene

The contribution of heavy-atom tunneling to reactions of [12]- and [16]annulene was probed using small-curvature tunneling rate calculations. At the CCSD(T)/cc-pVDZ//M06-2X/cc-pVDZ level, tunneling is predicted to account for more than 50% of the rate for Möbius bond shifting and ca. 35% of the rate for electrocyclization in [12]annulene, and over 80% of the rate for Möbius bond shifting in [16]annulene, at temperatures at which these reactions have been observed experimentally.

Tunneling by 16 Carbons: Planar Bond Shifting in [16]Annulene

Midsized annulenes are known to undergo rapid π-bond shifting. Given that heavy-atom tunneling plays a role in planar bond shifting of cyclobutadiene, we computationally explored the contribution of heavy-atom tunneling to planar π-bond shifting in the major (CTCTCTCT, 5a) and minor (CTCTTCTT, 6a) known isomers of [16]annulene. UM06-2X/cc-pVDZ calculations yield bond-shifting barriers of ca. 10 kcal/mol. The results also reveal extremely narrow barrier widths, suggesting a high probability of tunneling for these bond-shifting reactions. Rate constants were calculated using canonical variational transition state theory (CVT) as well as with small curvature tunneling (SCT) contributions, via direct dynamics. For the major isomer 5a, the computed SCT rate constant for bond shifting at 80 K is 0.16 s^-1, corresponding to a half-life of 4.3 s, and indicating that bond shifting is rapid at cryogenic temperatures despite a 10 kcal/mol barrier. This contrasts with the CVT rate constant of 8.0 × 10^-15 s^-1 at 80 K. The minor isomer 6a is predicted to undergo rapid bond shifting via tunneling even at 10 K. For both isomers, bond shifting is predicted to be much faster than competing conformation change despite lower barriers for the latter process. The preference for bond shifting represents cases of tunneling control in which the preferred reaction is dominated by heavy-atom motions. At all temperatures below -50 °C, tunneling is predicted to dominate the bond shifting process for both 5a and 6a. Thus, [16]annulene is predicted to be an example of tunneling by 16 carbons. Bond shifting in both isomers is predicted to be rapid at temperatures accessible by solution-phase NMR spectroscopy, and an experiment is proposed to verify these predictions.

Nonconjugated Hydrocarbons as Rigid-Linear Motifs: Isosteres for Material Sciences and Bioorganic and Medicinal Chemistry

Nonconjugated hydrocarbons, like bicyclo[1.1.1]pentane, bicyclo[2.2.2]octane, triptycene, and cubane are a unique class of rigid linkers. Due to their similarity in size and shape they are useful mimics of classic benzene moieties in drugs, so-called bioisosteres. Moreover, they also fulfill an important role in material sciences as linear linkers, in order to arrange various functionalities in a defined spatial manner. In this Review article, recent developments and usages of these special, rectilinear systems are discussed. Furthermore, we focus on covalently linked, nonconjugated linear arrangements and discuss the physical and chemical properties and differences of individual linkers, as well as their application in material and medicinal sciences.

Triplet-sensitised di-π-methane rearrangement of N-substituted 2-azabarrelenones

When irradiated at λ = 366 nm or at λ = 420 nm in the presence of an appropriate sensitiser the title compounds underwent a di-π-methane rearrangement which led to the formation of tricyclic azasemibullvalenones (2a,2a1,2b,4a-tetrahydroazacyclopropa[cd]pentalenones) in yields of 63-87%.

Redox-Neutral Photocatalytic Cyclopropanation via Radical/Polar Crossover

A benchtop stable, bifunctional reagent for the redox-neutral cyclopropanation of olefins has been developed. Triethylammonium bis(catecholato)iodomethylsilicate can be readily prepared on multigram scale. Using this reagent in combination with an organic photocatalyst and visible light, cyclopropanation of an array of olefins, including trifluoromethyl- and pinacolatoboryl-substituted alkenes, can be accomplished in a matter of hours. The reaction is highly tolerant of traditionally reactive functional groups (carboxylic acids, basic heterocycles, alkyl halides, etc.) and permits the chemoselective cyclopropanation of polyolefinated compounds. Mechanistic interrogation revealed that the reaction proceeds via a rapid anionic 3- exo- tet ring closure, a pathway consistent with experimental and computational data.

Reaction Electronic Flux Perspective on the Mechanism of the Zimmerman Di-π-methane Rearrangement

The reaction electronic flux (REF) offers a powerful tool in the analysis of reaction mechanisms. Noteworthy, the relationship between aromaticity and REF can eventually reveal subtle electronic events associated with reactivity in aromatic systems. In this work, this relationship was studied for the triplet Zimmerman di-π-methane rearrangement. The aromaticity loss and gain taking place during the reaction is well acquainted by the REF, thus shedding light on the electronic nature of reactions involving dibenzobarrelenes.

Does 8-Nitroguanine Form 8-Oxoguanine? An Insight from Its Reaction with •OH Radical

8-Nitroguanine (8-nitroG) formed due to nitration of guanine base of DNA plays an important role in mutagenesis and carcinogenesis. In the present contribution, state-of-the-art quantum chemical calculations using M06-2X density functional and domain-based local pair natural orbital-coupled cluster theory with single, double, and perturbative triple excitations (DLPNO-CCSD(T)) methods have been carried out to investigate the mechanism of reaction of ^•OH radical with 8-nitroG leading to the formation of 8-oxoguanine (8-oxoG) (one of the most mutagenic and carcinogenic derivatives of guanine) in gas phase and aqueous media. Calculations of barrier energies and rate constants involved in the addition reactions of ^•OH radical at different sites of 8-nitroguanine show that C8 and C2 sites are the most and least reactive sites, respectively. Relative stability and Boltzmann populations of adducts show that the adduct formed at the C8 site occurs predominantly in equilibrium. Our calculations reveal that 8-nitroG is very reactive toward ^•OH radical and is converted readily into 8-oxoG when attacked by ^•OH radicals, in agreement with available experimental observations.