Concerted Electronic and Nuclear Fluxes During Coherent Tunnelling in Asymmetric Double-Well Potentials

Restriction on molecular fluxionality by substitution: A case study for the 1,10-dicyanobullvalene

We show herein that 1,10-dicyano substitution restricts the paragon fluxionality of bullvalene to just 14 isomers which isomerize along a single cycle. The restricted fluxionality of 1,10-dicyanobullvalene (DCB) is investigated by means of: (i) Bonding analyses of the isomer structures using the adaptive natural density partitioning (AdNDP). (ii) Quantum dynamical simulations of the isomerizations along the cyclic intrinsic reaction coordinate of the potential energy surface (PES). The PES possesses 14 equivalent potential wells supporting 14 isomers which are separated by 14 equivalent potential barriers supporting 14 transition states. Accordingly, at low temperatures, DCB appears as a hindered molecular rotor, without any delocalization of the wavefunction in the 14 potential wells, without any nuclear spin isomers, and with completely negligible tunneling. These results are compared and found to differ from those for molecular boron rotors. (iii) Born-Oppenheimer molecular dynamics (BOMD) simulations of thermally activated isomerizations. (iv) Calculations of the rate constants in the frame of transition state theory (TST) with reasonable agreement achieved with the BOMD results. (v) Simulations of the equilibration dynamics using rate equations for the isomerizations with TST rate coefficients. Accordingly, in the long-time limit, isomerizations of the 14 isomers, each with Cs symmetry, approach the "14 Cs → C7v" thermally averaged structure. This is a superposition of the 14 equally populated isomer structures with an overall C7v symmetry. By extrapolation, the results for DCB yield working hypotheses for so far un-explored properties e.g. for the equilibration dynamics of C10H10.

Solvation Effects on Quantum Tunneling Reactions

A decisive factor for obtaining high yields and selectivities in organic synthesis is the choice of the proper solvent. Solvent selection is often guided by the intuitive understanding of transition state-solvent interactions. However, quantum-mechanical tunneling can significantly contribute to chemical reactions, circumventing the transition state and thus depriving chemists of their intuitive handle on the reaction kinetics. In this Account, we aim to provide rationales for the effects of solvation on tunneling reactions derived from experiments performed in cryogenic matrices.The tunneling reactions analyzed here cover a broad range of prototypical organic transformations that are subject to strong solvation effects. Examples are the hydrogen tunneling probability for the cis-trans isomerization of formic acid which is strongly reduced upon formation of hydrogen-bonded complexes and the [1,2]H-shift in methylhydroxycarbene where a change in product selectivity is predicted upon interaction with hydrogen bond acceptors.Not only hydrogen but also heavy atom tunneling can exhibit strong solvent effects. The direction of the nearly degenerate valence tautomerization between benzene oxide and oxepin was found to reverse upon formation of a halogen or hydrogen bond with ICF₃ or H₂O. But even in the absence of strong noncovalent interactions such as hydrogen or halogen bonding, solvation can have a decisive effect on tunneling as evidenced by the Cope rearrangement of semibullvalenes via heavy-atom tunneling. Can quantum tunneling be catalyzed? The acceleration of the ring expansion of 1H-bicyclo[3.1.0.]-hexa-3,5-dien-2-one by complexation with Lewis acids provides a proof-of-concept for tunneling catalysis.Two concepts are central for the explanation and prediction of solvation effects on tunneling phenomena: a simple approach expands the Born-Oppenheimer approximation by separating nuclear degrees of freedom into intra- and intermolecular degrees. Intermolecular movements represent the slowest motions within molecular aggregates, thus effectively freezing the position of the solvent in relation to the reactant during the tunneling process. Another useful approach is to treat reactants and products by separate single-well potentials, where the intersection represents the transition state. Thus, stabilization of the reactants via solvation should result in an increase in barrier heights and widths which in turn lowers tunneling probabilities. These simple models can predict trends in tunneling kinetics and provide a rational basis for controlling tunneling reactions via solvation.

Hydrogen Delocalization in an Asymmetric Biomolecule: The Curious Case of Alpha-Fenchol

Rotational microwave jet spectroscopy studies of the monoterpenol α-fenchol have so far failed to identify its second most stable torsional conformer, despite computational predictions that it is only very slightly higher in energy than the global minimum. Vibrational FTIR and Raman jet spectroscopy investigations reveal unusually complex OH and OD stretching spectra compared to other alcohols. Via modeling of the torsional states, observed spectral splittings are explained by delocalization of the hydroxy hydrogen atom through quantum tunneling between the two non-equivalent but accidentally near-degenerate conformers separated by a low and narrow barrier. The energy differences between the torsional states are determined to be only 16(1) and 7(1) cm-1hc for the protiated and deuterated alcohol, respectively, which further shrink to 9(1) and 3(1) cm-1hc upon OH or OD stretch excitation. Comparisons are made with the more strongly asymmetric monoterpenols borneol and isopinocampheol as well as with the symmetric, rapidly tunneling propargyl alcohol. In addition, the third-in contrast localized-torsional conformer and the most stable dimer are assigned for α-fenchol, as well as the two most stable dimers for propargyl alcohol.

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.

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.

Open system dynamics using Gaussian-based multiconfigurational time-dependent Hartree wavefunctions: Application to environment-modulated tunneling

A variational approach for the quantum dynamics of statistical mixtures is developed, which is based upon the representation of the natural states of the mixture in terms of hybrid Gaussian-based Multiconfiguration Time-Dependent Hartree (G-MCTDH) wavefunctions. The method, termed ρG-MCTDH, is combined with a treatment of dissipation and decoherence based on the nonstochastic open-system Schrödinger equations. The performance and the convergence properties of the approach are illustrated for a two-dimensional tunneling system, where the primary tunneling coordinate, represented by flexible single-particle functions, is resonantly coupled to a second harmonic mode, represented by Gaussian wave packets. The harmonic coordinate is coupled to the environment and two different processes are studied: (i) vibrational relaxation at zero temperature described by a master equation in the Lindblad form and (ii) thermalization induced by the Caldeira-Leggett master equation. In the second case, the evolution from a quantum tunneling regime to a quasistationary classical-limit distribution, driven by the heat bath, is visualized using a flux analysis.

Curly arrows, electron flow, and reaction mechanisms from the perspective of the bonding evolution theory

Despite the usefulness of curly arrows in chemistry, their relationship with real electron density flows is still imprecise, and even their direct connection to quantum chemistry is still controversial. The paradigmatic description - from first principles - of the mechanistic aspects of a given chemical process is based mainly on the relative energies and geometrical changes at the stationary points of the potential energy surface along the reaction pathway; however, it is not sufficient to describe chemical systems in terms of bonding aspects. Probing the electron density distribution during a chemical reaction can provide important insights, enabling us to understand and control chemical reactions. This aim has required an extension of the relationships between the concepts of traditional chemistry and those of quantum mechanics. Bonding evolution theory (BET), which combines the topological analysis of the electron localization function (ELF) and Thom's catastrophe theory (CT), provides a powerful method that offers insight into the molecular mechanism of chemical rearrangements. In agreement with the laws of physical and aspects of quantum theory, BET can be considered an appropriate tool to tackle chemical reactivity with a wide range of possible applications. In this work, BET is applied to address a long-standing problem: the ability to monitor the flow of electron density. BET analysis shows a connection between quantum mechanics and bond making/forming processes. Likewise, the present approach retrieves the classical curly arrows used to describe the rearrangements of chemical bonds and provides detailed physical grounds for this type of representation. We demonstrate this procedure using the test set of prototypical examples of thermal ring apertures, and the degenerated Cope rearrangement of semibullvalene.