Elucidating Heavy Atom Tunneling Kinetics in the Cope Rearrangement of Semibullvalene.

In this work, we characterize the temperature dependence of kinetic properties in heavy atom tunneling reactions by means of molecular dynamics simulations, including nuclear quantum effects (NQEs) via Path Integral theory. To this end, we consider the prototypical Cope rearrangement of semibullvalene. The reaction was studied in the 25-300K temperature range observing that the inclusion of NQEs modifies the temperature behavior of both free energy barriers and dynamical recrossing factors with respect to classical dynamics. Notably, while in classical simulations the activation free energy shows a very weak temperature dependence, it becomes strongly dependent on temperature when NQEs are included. This temperature behavior shows a transition from a regime where the quantum effects are limited and can mainly be traced back to zero point energy, to a low temperature regime where tunneling plays a dominant role. In this regime, the free energy curve tunnels below the potential energy barrier along the reaction coordinate,  allowing much faster reaction rates. Finally, the temperature dependence of the rate constants obtained from molecular dynamics simulations was compared with available experimental data and with semi-classical transition state theory calculations, showing comparable behaviors and similar transition temperatures from  thermal to (deep) tunneling regime.

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.

Boron Tunneling in the "Weak" Bond-Stretch Isomerization of N-B Lewis Adducts.

Some nitrile-boron halide adducts exhibit a double-well potential energy surface with two distinct minima: a "long bond" geometry (LB, a van der Waals interaction mostly based on electrostatics, but including a residual charge transfer component) and a "short bond" structure (SB, a covalent dative bond). This behavior can be considered as a "weak" form of bond stretch isomerism. Our computations reveal that complexes RCN-BX3 (R=CH3 , FCH2 , BrCH2 , and X=Cl, Br) exhibit a fast interconversion from LB to SB geometries even close to the absolute zero thanks to a boron atom tunneling mechanism. The computed half-lives of the meta-stable LB compounds vary between minutes to nanoseconds at cryogenic conditions. Accordingly, we predict that the long bond structures are practically impossible to isolate or characterize, which agrees with previous matrix-isolation experiments.

Heavy-Atom Tunneling Through Crossing Potential Energy Surfaces: Cyclization of a Triplet 2-Formylarylnitrene to a Singlet 2,1-Benzisoxazole.

Not long ago, the occurrence of quantum mechanical tunneling (QMT) chemistry involving atoms heavier than hydrogen was considered unreasonable. Contributing to the shift of this paradigm, we present here the discovery of a new and distinct heavy-atom QMT reaction. Triplet syn-2-formyl-3-fluorophenylnitrene, generated in argon matrices by UV-irradiation of an azide precursor, was found to spontaneously cyclize to singlet 4-fluoro-2,1-benzisoxazole. Monitoring the transformation by IR spectroscopy, temperature-independent rate constants (k≈1.4×10-3  s-1 ; half-life of ≈8 min) were measured from 10 to 20 K. Computational estimated rate constants are in fair agreement with experimental values, providing evidence for a mechanism involving heavy-atom QMT through crossing triplet to singlet potential energy surfaces. Moreover, the heavy-atom QMT takes place with considerable displacement of the oxygen atom, which establishes a new limit for the heavier atom involved in a QMT reaction in cryogenic matrices.

Heavy-Atom Tunneling in Semibullvalenes: How Driving Force, Substituents, and Environment Influence the Tunneling Rates.

The Cope rearrangement of selectively deuterated isotopomers of 1,5-dimethylsemibullvalene 2 a and 3,7-dicyano-1,5-dimethylsemibullvalene 2 b were studied in cryogenic matrices. In both semibullvalenes the Cope rearrangement is governed by heavy-atom tunneling. The driving force for the rearrangements is the small difference in the zero-point vibrational energies of the isotopomers. To evaluate the effect of the driving force on the tunneling probability in 2 a and 2 b, two different pairs of isotopomers were studied for each of the semibullvalenes. The reaction rates for the rearrangement of 2 b in cryogenic matrices were found to be smaller than the ones of 2 a under similar conditions, whereas differences in the driving force do not influence the rates. Small curvature tunneling (SCT) calculations suggest that the reduced tunneling rate of 2 b compared to that of 2 a results from a change in the shape of the potential energy barrier. The tunneling probability of the semibullvalenes strongly depends on the matrix environment; however, for 2 a in a qualitatively different way than for 2 b.

Fast Heavy-Atom Tunneling in Trifluoroacetyl Nitrene.

Chemical reactions involving quantum mechanical tunneling (QMT) increasingly attract the attention of scientists. In contrast to the hydrogen-tunneling as frequently observed in chemistry and biology, tunneling solely by heavy atoms is rare. Herein, we report heavy-atom tunneling in trifluoroacetyl nitrene, CF3 C(O)N. The carbonyl nitrene CF3 C(O)N in the triplet ground state was generated in cryogenic matrices by laser (193 or 266 nm) photolysis of CF3 C(O)N3 and characterized by IR and EPR spectroscopy. In contrast to the theoretically predicted activation barriers (>10 kcal mol-1 ), CF3 C(O)N undergoes rapid rearrangement into CF3 NCO with half-life times of less than 10 min and unprecedentedly large 14 N/15 N kinetic isotope effects (1.18-1.33) in solid Ar, Ne, and N2 matrices even at 2.8 K. The tunneling disappearance of CF3 C(O)N becomes much slower in the chemically active toluene and in 2-methyltetrahydrofuran at 5 K.

Heavy-Atom Tunneling Calculations in Thirteen Organic Reactions: Tunneling Contributions are Substantial, and Bell's Formula Closely Approximates Multidimensional Tunneling at ≥250†K.

Multidimensional tunneling calculations are carried out for 13 reactions, to test the scope of heavy-atom tunneling in organic chemistry, and to check the accuracy of one-dimensional tunneling models. The reactions include pericyclic, cycloaromatization, radical cyclization and ring opening, and SN 2. When compared at the temperatures that give the same effective rate constant of 3×10-5  s-1 , tunneling accounts for 25-95 % of the rate in 8 of the 13 reactions. Values of transmission coefficients predicted by Bell's formula, κBell , agree well with multidimensional tunneling (canonical variational transition state theory with small curvature tunneling), κSCT . Mean unsigned deviations of κBell vs. κSCT are 0.08, 0.04, 0.02 at 250, 300 and 400 K. This suggests that κBell is a useful first choice for predicting transmission coefficients in heavy-atom tunnelling.