Tunneling Control: Competition between 6π-Electrocyclization and [1,5]H-Sigmatropic Shift Reactions in Tetrahydro-1H-cyclobuta[e]indene Derivatives

Perturbing π-clouds with Substituents to Study the Effects on Reaction Dynamics of gauche-1,3-Butadiene to Bicyclobutane Electrocyclization

The conical intersection (CI) governs the ultra-fast relaxation of excited states in a radiationless manner and are observed mainly in photochemical processes. In the current work, we investigated the effects of substituents on the reaction dynamics for the conversion of gauche-1,3-butadiene to bicyclobutane via photochemical electrocyclization. We incorporated both electron withdrawing (-F) and donating (-CH₃ ) groups in the conjugated system. In our study, we optimized the minimum energy conical intersection (MECI) geometries using the multi-configurational state-averaged CASSCF approach, whereas, to study the ground state reaction pathways for the substituted derivatives, dispersion corrected, B3LYP-D3 functional was used. The non-adiabatic surface hopping molecular dynamics simulations were performed to observe the behaviour of electronic states involved throughout the photoconversion process. The results obtained from the multi-reference second-order perturbation correction of energy at the XMS-CASPT2 level of theory, topography analysis, and non-adiabatic dynamics suggest that the -CH₃ substituted derivatives can undergo faster thermal conversion to the product in the ground state with a smaller activation energy barrier compared to -F substituted derivative. Our study also reveals that the GBUT to BIBUT conversion follows both conrotatory and disrotatory pathways, whereas, on substitution with -F or -CH₃ , the conversion proceeds via the conrotatory pathway.

Simultaneous Tunneling Control in Conformer-Specific Reactions

We present here a new example of chemical reactivity governed by quantum tunneling, which also highlights the limitations of the classical theories. The syn and anti conformers of a triplet 2-formylphenylnitrene, generated in a nitrogen matrix, were found to spontaneously rearrange to the corresponding 2,1-benzisoxazole and imino-ketene, respectively. The kinetics of both transformations were measured at 10 and 20 K and found to be temperature-independent, providing clear evidence of concomitant tunneling reactions (heavy-atom and H-atom). Computations confirm the existence of these tunneling reaction pathways. Although the energy barrier between the nitrene conformers is lower than any of the observed reactions, no conformational interconversion was observed. These results demonstrate an unprecedented case of simultaneous tunneling control in conformer-specific reactions of the same chemical species. The product outcome is impossible to be rationalized by the conventional kinetic or thermodynamic control.

Controlling Tunneling Oscillation and Quantum Localization in an Asymmetric Double-Well Potential: A Bohmian Perspective

The roles of spatial symmetry and strength of external time-dependent perturbation on the dynamics of a quantum particle, initially localized in one of the wells of an asymmetric double-well potential are studied using the recently developed techniques incorporating quantum theory of motion and time-dependent Fourier grid Hamiltonian methods. The model used here includes a mimic of the related experimental situations which is considered as a perturbation to the static double-well potential. Analysis of localized and delocalized phase space structures and corresponding time-profile of tunneling probability reveal the recipe toward controlling the tunneling oscillations by modulating the parameters of applied perturbation. A study on a stochastic pulsating potential also reveals the root to the quantum localization, even in moderate field strength.

Computational Evidence for Tunneling and a Hidden Intermediate in the Biosynthesis of Tetrahydrocannabinol

Quantum tunneling is computed for a reaction sequence that models the conversion of the ortho-quinone methide of cannabigerolic acid 1 to the decarboxylated product (-)-trans-Δ⁹-tetrahydrocannabinol (THC, 3). This calculation is the first to evaluate multidimensional tunneling in this sequence. Computations were carried out with POLYRATE and GAUSSRATE using B3LYP/6-31G(d,p) to examine the mechanism of THC 3 formation. The pentyl chain on THC 3 and its precursors were replaced with a methyl group to compute tunneling contributions to the rates of four separate steps: (i) initial Diels-Alder reaction of the quinone methide with the trisubstituted alkene end-group of the geranyl 1Z-CH₃ to give 2Z-CH₃, (ii) acid-catalyzed keto-enol tautomerization, which converts 2rZ-CH₃ to 4rZ-CH₃, (iii) carboxyl rotamerization converting 4rZ-CH₃ to 4E-CH₃, and (iv) decarboxylation that converts 4E-CH₃ to 3-CH₃. Tunneling contributions to the rate constants of steps (i)-(iv) are between 19 and 76% at 293 K. In step (ii), nonuniform changes in the zero-point vibrational energy along the reaction path created a shallow minimum in the 0 K free energy. It is a hidden intermediate because it is not a minimum on the potential energy surface and is detectable only when zero-point energy is taken into account along the reaction path. Predicted kinetic isotope effects would be experimentally observable at temperatures that are convenient to use. This is particularly relevant in the decarboxylation stage of the reaction sequence and has important implications because of its role in THC 3 formation.

Designing Force Probes Based on Reversible 6π-Electrocyclizations in Polyenes Using Quantum Chemical Calculations

The conjugated π-system in polyenes can be interrupted by electrocyclic ring-closure reactions. In this work, this 6π-electrocylization is shown by means of density functional calculations to be reversible by the application of an external mechanical pulling force at the terminal ends of the interrupted polyene chain. The test systems were constrained in a fused ring system, thus locking the orientation of three π-bonds and generally promoting 6π-electrocyclic ring-closure reactions. For several systems, the forward reaction is exergonic and the corresponding reaction barrier is comparable to those reported in the literature. The reverse reaction is triggered by an external pulling force of 2 nN (nano-Newton) or less and also becomes exergonic in all investigated polyenes under these force conditions. Moreover, it proceeds via a low reaction barrier when a pulling force of 2 nN is active, indicating that the mechanical force is an efficient stimulus for triggering ring-opening reactions. Analysis of the strain energy induced by this mechanical force confirms an optimal activation of the corresponding C-C σ-bond that breaks upon ring opening when the pulling positions are located on the polyene chain.

Hydrogen tunnelling in the rearrangements of carbenes: the role of dynamical calculations

A tunnelling controlled reaction is studied with semiclassical transition state theory, rationalising the results of experiment.

The unimolecular dissociation of magnesium chloride squarate (ClMgC4O4−) and reductive cyclooligomerisation of CO on magnesium

Cyclooligomerisation to squarate is initiated by the coupling of two CO units in MgCl⁻.

Structural, Electronic, and Spectral Properties of Metal Dimethylglyoximato [M(DMG)2; M = Ni2+, Cu2+] Complexes: A Comparative Theoretical Study

Dimethylglyoxime (DMG) usually forms thermodynamically stable chelating complexes with selective divalent transition-metal ions. Electronic and spectral properties of metal-DMG complexes are highly dependent on the nature of metal ions. Using range-separated hybrid functional augmented with dispersion corrections within density functional theory (DFT) and time-dependent DFT, we present a detailed and comprehensive study on structural, electronic, and spectral (both IR and UV-vis) properties of M(DMG)₂ [M = Ni²⁺, Cu²⁺] complexes. Ni(DMG)₂ results are thoroughly compared with Cu(DMG)₂ and also against available experimental data. Stronger H-bonding leads to greater stability of Ni(DMG)₂ with respect to isolated ions (M²⁺ and DMG^-) compared to Cu(DMG)₂. In contrast, a relatively larger reaction enthalpy for Cu(DMG)₂ formation from chemically relevant species is found than that of Ni(DMG)₂ because of the greater binding enthalpy of [Ni(H₂O)₆]²⁺ than that of [Cu(H₂O)₆]²⁺. In dimers, Ni(DMG)₂ is found to be 6 kcal mol^-1 more stable than Cu(DMG)₂ due to a greater extent of dispersive interactions. Interestingly, a modest ferromagnetic coupling (588 cm^-1) is predicted between two spin-¹/₂ Cu²⁺ ions present in the Cu(DMG)₂ dimer. Additionally, the potential energy curves calculated along the O-H bond coordinate for both complexes suggest asymmetry and symmetry in the H-bonding interactions between the H-bond donor and acceptor O centers in the solid-state and in solution, respectively, well corroborating with early experimental findings. Interestingly, a lower proton transfer barrier is obtained for the Ni(DMG)₂ compared to its Cu-analogue due to stronger H-bonding in the former complex. In fact, relatively weaker H-bonding in Cu(DMG)₂ results in blue-shifted O-H stretching modes compared to that in Ni(DMG)₂. On the other hand, qualitatively similar optical absorption spectra are obtained for both complexes with red-shifted peaks found for the Cu(DMG)₂. Finally, computational models for axial mono- and diligand (aqua and ammonia) coordinated M(DMG)₂ complexes are predicted to be energetically feasible and stable with relatively greater binding stability obtained for the ammonia-coordination.

Thermal decomposition of syn- and anti-dihydropyrenes; functional group-dependent decomposition pathway

Syn and anti dihydropyrene (DHP) are excellent thermochromes, and therefore extensively studied for their thermochromic and photochromic properties, respectively. However, they suffer from thermal decomposition due to thermal instability. In this study, we thoroughly investigated pathways for the thermal decomposition of anti- and syn- dihydropyrenes through computational methods. The decomposition pathways include sigmatropic shift and hemolytic and heterolytic (cationic and anionic) cleavages. The decomposition pathway is influenced not only by the dihydropyrene (syn- or anti-) but also by the functional groups present. For anti-dihydropyrenes, sigmatropic shift is the most plausible pathways for CN and CHO internal groups. The cascade of sigmatropic shifts is followed by elimination to deliver substituted pyrenes. For CH₃- and H- dihydropyrenes, hemolytic cleavage of the internal groups is the most plausible pathway for decomposition to pyrenes. The pathway is changed to heterolytic cleavage when the internal groups on the dihydropyrenes are Cl^-, Br^-, and SMe^-. Comparison of the activation barriers for syn (30.18 kcal mol^-1) and anti (32.10 kcal mol^-1) dimethyldihydropyrenes for radical pathway reveal that decomposition of syn- DHP is more facile over anti-, which is consistent with the experimental observation. The decomposition pathway for syn-dihydropyrene is also hemolytic in cleavage when the internal groups are methyl and hydrogen. Syn-dihydropyrenes (symmetrical or unsymmetrical) bearing CN group do not follow sigmatropic shift, quite contrary to the anti-dihydropyrene. The lack of tendency of the syn-dihydropyrene for sigmatropic shift is rationalized on the planarity of the scaffold. The results of the theoretical study are consistent with the experimental observations. The results here help in understanding the behavior of substituents on the dihydropyrene scaffold, which will be useful in designing new molecules with improved thermal stabilities. Graphical abstract Functional group dependent decomposition pathways of dihydropyrenes.

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.

Computational elucidation of the reaction mechanism for synthesis of pyrrolidinedione derivatives via Nef-type rearrangement – cyclization reaction

This paper reports a quantum chemical study of all stages of a one-pot synthesis of pyrrolidinedione derivatives from nitromethane and coumarin, which includes Michael addition, migration of an oxygen atom (Nef-type rearrangement), and cyclization to a pyrrolidine ring. The energy barrier of deprotonated nitromethane addition to coumarin is 21.7 kJ mol⁻¹, while the barrier of proton transfer from the methylene to the nitro group in the nitromethyl group is notably higher, 197.8 kJ mol⁻¹. The second stage of the reaction, migration of an oxygen atom within the nitromethyl group, occurs with lowest energy barrier, 142.4 kJ mol⁻¹, when it is assisted by an additional water molecule. The last stage – cyclization, passes with a very low energy barrier of 11.9 kJ mol⁻¹ but the tautomerization of the nitrosohydroxymethyl group to the hydroxy-N-hydroxyiminomethyl, necessary for the process, has an energy barrier of 178.4 kJ mol⁻¹. Analogous calculations for the same process with the ethyl ester of 3-coumarin-carboxylic acid as substrate show that the relative energies of the intermediates and transition states are by at most 10–16 kJ mol⁻¹ more stable than the corresponding structures with coumarin.

A quantum chemical study of all stages of pyrrolidinedione derivative synthesis from nitromethane and coumarin, which includes Michael addition, migration of an oxygen atom (Nef-type rearrangement), and cyclization to a pyrrolidine ring is reported.

Isotope-Controlled Selectivity by Quantum Tunneling: Hydrogen Migration versus Ring Expansion in Cyclopropylmethylcarbenes

Using the tunneling-controlled reactivity of cyclopropylmethylcarbene, we demonstrate the viability of isotope-controlled selectivity (ICS), a novel control element of chemical reactivity where a molecular system with two conceivable products of tunneling exclusively produces one or the other, depending only on isotopic composition. Our multidimensional small-curvature tunneling (SCT) computations indicate that, under cryogenic conditions, 1-methoxycyclopropylmethylcarbene shows rapid H-migration to 1-methoxy-1-vinylcyclopropane, whereas deuterium-substituted 1-methoxycyclopropyl-d₃-methylcarbene undergoes ring expansion to 1-d₃-methylcyclobutene. This predicted change in reactivity constitutes the first example of a kinetic isotope effect that discriminates between the formation of two products.

Experimental and computational kinetic investigations for the reactions of Cl atoms with unsaturated ketones in the gas phase

Cl atom initiated photo oxidation of unsaturated ketones.

Kinetic investigations of Cl atom initiated photo-oxidation reactions of cyclic unsaturated hydrocarbons in the gas phase: an experimental and theoretical study

Cl atom initiated photo oxidation kinetics of cyclohexene and cycloheptene.