Heavy-Atom Tunneling Processes during Denitrogenation of 2,3-Diazabicyclo[2.2.1]hept-2-ene and Ring Closure of Cyclopentane-1,3-diyl Diradical. Stereoselectivity in Tunneling and Matrix Effect

Heavy-atom tunnelling in benzene isomers: how many tricyclic species are truly stable?

The variety of possible benzene isomers may provide a fundamental basis for understanding structural and reactivity patterns in organic chemistry. However, the vast majority of these isomers remain unsynthesized, while most of the experimentally known species are only moderately stable. Consequently, there is a high probability that the theoretically proposed isomers would also be barely metastable, a factor that must be taken into account if their creation in the laboratory is sought. In this work, we studied the kinetic stability of all 73 hypothetical tricyclic benzene isomers, especially focusing on their nuclear quantum effects. With this in mind, we evaluated which species are theoretically possible to synthesize, detect, and isolate. Our computations predict that 26% of the previously deemed stable molecules are completely unsynthesizable due to their intrinsic quantum tunnelling instability pushing for their unimolecular decomposition even close to the absolute zero. Five more systems would be detectable, but they will slowly and inevitably degrade, while seven more supposedly stable systems will break apart in barrierless mechanisms.

A Perspective on the Investigation of Spectroscopy and Kinetics of Complex Molecular Systems with Semiclassical Approaches

In this Perspective we show that semiclassical methods provide a rigorous hierarchical way to study the vibrational spectroscopy and kinetics of complex molecular systems. The time averaged approach to spectroscopy and the semiclassical transition state theory for kinetics, which have been first adopted and then further developed in our group, provide accurate quantum results on rigorous physical grounds and can be applied even when dealing with a large number of degrees of freedom. In spectroscopy, the multiple coherent, divide-and-conquer, and adiabatically switched semiclassical approaches have practically permitted overcoming issues related to the convergence of results. In this Perspective we demonstrate the possibility of studying the semiclassical vibrational spectroscopy of a molecule adsorbed on an anatase (101) surface, a system made of 51 atoms. In kinetics, the semiclassical transition state theory is able to account for anharmonicity and the coupling between the reactive and bound modes. Our group has developed this technique for practical applications involving the study of phenomena like kinetic isotope effect, heavy atom tunneling, and elusive conformer lifetimes. Here, we show that our multidimensional anharmonic quantum approach is able to tackle on-the-fly the thermal kinetic rate constant of a 135 degree-of-freedom system. Overall, semiclassical methods open up the possibility to describe at the quantum mechanical level systems characterized by hundreds of degrees of freedom leading to the accurate spectroscopic and kinetic description of biomolecules and complex molecular systems.

Differential Tunneling-Driven and Vibrationally-Induced Reactivity in Isomeric Benzazirines

Quantum mechanical tunneling of heavy-atoms and vibrational excitation chemistry are unconventional and scarcely explored types of reactivity. Once fully understood, they might bring new avenues to conduct chemical transformations, providing access to a new world of molecules or ways of exquisite reaction control. In this context, we present here the discovery of two isomeric benzazirines exhibiting differential tunneling-driven and vibrationally-induced reactivity, which constitute exceptional results for probing into the nature of these phenomena. The isomeric 6-fluoro- and 2-fluoro-4-hydroxy-2H-benzazirines (3-a and 3'-s) were generated in cryogenic krypton matrices by visible-light irradiation of the corresponding triplet nitrene 3 2-a, which was produced by UV-light irradiation of its azide precursor. The 3'-s was found to be stable under matrix dark conditions, whereas 3-a spontaneously rearranges (τ1/2 ∼64 h at 10 and 20 K) by heavy-atom tunneling to 3 2-a. Near-IR-light irradiation at the first OH stretching overtone frequencies (remote vibrational antenna) of the benzazirines induces the 3'-s ring-expansion reaction to a seven-member cyclic ketenimine, but the 3-a undergoes 2H-azirine ring-opening reaction to triplet nitrene 3 2-a. Computations demonstrate that 3-a and 3'-s have distinct reaction energy profiles, which explain the different experimental results. The spectroscopic direct measurement of the tunneling of 3-a to 3 2-a constitutes a unique example of an observation of a species reacting only by nitrogen tunneling. Moreover, the vibrationally-induced sole activation of the most favorable bond-breaking/bond-forming pathway available for 3-a and 3'-s provides pioneer results regarding the selective nature of such processes.

Transformations of Strained Three-Membered Rings a Common, Yet Overlooked, Motif in Heavy-Atom Tunneling Reactions

Quantum mechanical tunneling has long been recognized as an important phenomenon when considering transformations dominated by a lightweight hydrogen atom. Tunneling of heavier atoms like carbon, initially dismissed as negligible, has seen a quickly increasing number of computationally predicted and/or experimentally confirmed examples over the last decade, thus highlighting its importance for a wide variety of reactions. However, no common structural motif has been pointed out within these seemingly unconnected examples, strongly limiting the predictability of the impact of heavy-atom tunneling on a given reaction. This Concept article will provide this perspective and showcase how the recognition of the formation and cleavage of three-membered rings as common motif can inform the prediction of and research into heavy-atom tunneling reactions.

Heavy Atom Tunneling in Organic Reactions at Coupled Cluster Potential Accuracy with a Parallel Implementation of Anharmonic Constant Calculations and Semiclassical Transition State Theory

We describe and test on some organic reactions a parallel implementation strategy to compute anharmonic constants, which are employed in semiclassical transition state theory reaction rate calculations. Our software can interface with any quantum chemistry code capable of a single point energy estimate, and it is suitable for both minimum and transition state geometry calculations. After testing the accuracy and comparing the efficiency of our implementation against other software, we use it to estimate the semiclassical transition state theory (SCTST) rate constant of three reactions of increasing dimensionality, known as examples of heavy atom tunneling. We show how our method is improved in efficiency with respect to other existing implementations. In conclusion, our approach allows SCTST rates and heavy atom tunneling at a high level of electronic structure theory (up to CCSD(T)) to be evaluated. This work shows how crucial the possibility to perform high level ab initio rate evaluations can be.

Stereoelectronic and dynamical effects dictate nitrogen inversion during valence isomerism in benzene imine

Non-classical processes such as heavy-atom tunneling and post transition-state dynamics govern stereoselectivity for benzene imine ⇌ 1H-azepine.

Direct Detection of Singlet Cyclopentane-1,3-diyl Diradicals By Infrared and Ultraviolet-Visible Spectroscopy at Cryogenic Temperature and Their Photoreactivity

Photolysis of a 7,7-difluoro-1,4-diphenyl-2,3-diazabicyclo[2.2.1]hept-2-ene derivative (AZ1) using a 365 nm light-emitting diode in an Ar matrix at 4 K resulted in the formation of a planar singlet 2,2-difluoro-1,3-diphenylcyclopentane-1,3-diyl diradical derivative, S-DR1-pl (λ_max = 520 nm). A singlet cyclopentane-1,3-diyl diradical system (S-DR1-pl) was directly detected by steady-state infrared (IR) spectroscopy. Due to the photolability of S-DR1-pl, initial photolysis of AZ1 also yielded the ring-closed product ret-CP1 and migration products trans-MG1 and/or cis-MG1, which were observed using IR spectra. Monitoring of prolonged photolysis using IR and ultraviolet-visible (UV-vis) spectra demonstrated the formation of the allylic cation CT1 (λ_max = 470 nm). On the other hand, photolysis of a 7,7-dimethoxy-1,4-diphenyl-2,3-diazabicyclo[2.2.1]hept-2-ene derivative (AZ2) yielded a puckered conformer (instead of planar) of the corresponding diradical S-DR2-puc, which was detected by IR and UV-vis spectroscopy in an Ar matrix at 4 K. This spectroscopic characterization opens a new strategy to obtain more detailed information about the structure and reactivity of singlet cyclopentane-1,3-diyl diradicals.

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.

Switch chemistry at cryogenic conditions: quantum tunnelling under electric fields

While the influence of intramolecular electric fields is a known feature in enzymes, the use of oriented external electric fields (EEF) to enhance or inhibit molecular reactivity is a promising topic still in its infancy. Herein we will explore computationally the effects that EEF can provoke in simple molecules close to the absolute zero, where quantum tunnelling (QT) is the sole mechanistic option. We studied three exemplary systems, each one with different reactivity features and known QT kinetics: π bond-shifting in pentalene, Cope rearrangement in semibullvalene, and cycloreversion of diazabicyclohexadiene. The kinetics of these cases depend both on the field strength and its direction, usually giving subtle but remarkable changes. However, for the cycloreversion, which suffers large changes on the dipole through the reaction, we also observed striking results. Between the effects caused by the EEF on the QT we observed an inversion of the Arrhenius equation, deactivation of the molecular fluxionality, and stabilization or instantaneous decomposition of the system. All these effects may well be achieved, literally, at the flick of a switch.

Impact of the macrocyclic structure and dynamic solvent effect on the reactivity of a localised singlet diradicaloid with π-single bonding character

Localised singlet diradicals are key intermediates in bond homolysis processes. Generally, these highly reactive species undergo radical–radical coupling reaction immediately after their generation. Therefore, their short-lived character hampers experimental investigations of their nature. In this study, we implemented the new concept of “stretch effect” to access a kinetically stabilised singlet diradicaloid. To this end, a macrocyclic structure was computationally designed to enable the experimental examination of a singlet diradicaloid with π-single bonding character. The kinetically stabilised diradicaloid exhibited a low carbon–carbon coupling reaction rate of 6.4 × 10³ s⁻¹ (155.9 μs), approximately 11 and 1000 times slower than those of the first generation of macrocyclic system (7.0 × 10⁴ s⁻¹, 14.2 μs) and the parent system lacking the macrocycle (5 × 10⁶ s⁻¹, 200 ns) at 293 K in benzene, respectively. In addition, a significant dynamic solvent effect was observed for the first time in intramolecular radical–radical coupling reactions in viscous solvents such as glycerin triacetate. This theoretical and experimental study demonstrates that the stretch effect and solvent viscosity play important roles in retarding the σ-bond formation process, thus enabling a thorough examination of the nature of the singlet diradicaloid and paving the way toward a deeper understanding of reactive intermediates.

An extremely long-lived localised singlet diradical with π-single bonding character is found in a macrocyclic structure that retards the radical–radical coupling reaction by the “stretch and solvent-dynamic effects”.