Solvation Effects on Quantum Tunneling Reactions

H-Tunneling Rotamerization in Glycine Imine

Quantum mechanical tunneling governs the chemical reactivity at cryogenic temperatures. Here we present the near-infrared (NIR) light-induced generation of a higher energy conformer of glycine imine and its H-tunneling CO bond rotamerization in solid argon (Ar), para-hydrogen (p-H2), and dinitrogen (N2) at cryogenic temperatures. The tunneling half-life for the CO bond rotamerization highly depends on the host matrix and is approximately 5 h in Ar and 18 h in both p-H2 and N2. Surprisingly, experiments in p-H2 revealed a much longer half-life than in Ar, possibly due to the formation of a dinitrogen complex. Deuteration of the carboxylic acid group completely inhibits the CO bond D-tunneling rotamerization. We performed Wentzel - Kramer - Brillouin (WKB) and canonical variational transition state theory (CVT) calculations, incorporating multidimensional small curvature tunneling corrections (SCT), using the B3LYP/cc-pVTZ level of theory. The gas-phase tunneling half-lives are 11 h according to the one-dimensional WKB model and 1 h according to the multidimensional CVT/SCT model, both aligning well with our experimental results.

Hydrogen-Bonded Complexes of HPN⋅ and HNP⋅ Radicals with Carbon Monoxide

Phosphorus mononitride (PN) is a carrier of phosphorus in the interstellar medium. As the simplest derivatives of PN, the radical species HPN⋅ and HNP⋅ have remained elusive. Herein, we report the generation, characterization, and photochemistry of HPN⋅ and HNP⋅ in N2-matrix at 3 K. Specifically, HPN⋅ was formed as a weakly bonded complex with CO in the matrix by 254 nm photolysis of the novel phosphinyl radical HPNCO⋅. The ⋅NPH-CO complex is extremely unstable, as it undergoes spontaneous isomerization to the lower-energy isomer ⋅PNH-CO through fast quantum mechanical tunneling (QMT) with a half-life of 6.1 min at 3 K. Upon further irradiation at 254 nm, the reverse conversion of ⋅PNH-CO to ⋅NPH-CO along with dehydrogenation to yield PN was observed. The characterization ⋅NPH-CO and ⋅PNH-CO with matrix-isolation IR spectroscopy is supported by D, 15N, and 13C isotope labeling and quantum chemical calculations at the XYGJ-OS/AVTZ level of theory, and the mechanism by hydrogen atom tunneling is consistent with multidimensional instanton theory calculations.

Quantum Tunneling: History and Mystery of Large Amplitude Motions over a Century

Large amplitude motions (LAMs), most notably represented by proton tunneling, mark a significant departure from small amplitude vibrations where protons merely oscillate around their equilibrium positions. These substantial displacements require tunneling through potential energy barriers, leading to splittings in, e.g., rotational spectra. Since Hund's pioneering work in 1927, proton tunneling has offered a unique glimpse into the internal dynamics of gas-phase molecules, with microwave spectroscopy being the key technique for such investigations. The ubiquous LAM type is methyl internal rotation, characterized by 3-fold potentials arising from the interaction between methyl rotors and their molecular frame, with the barrier hindering methyl torsion and the orientation of the torsional axis being defining features. Investigating methyl internal rotations plays a key role in fields ranging from molecular physics, where the methyl rotor serves as a sensitive probe for molecular structures, to atmospheric chemistry and astrophysics, where methyl-containing species have been detected in the Earth's atmosphere and interstellar environments and even discussed as potential probes for effects beyond the standard model of physics. Despite nearly a century of study, modeling methyl internal rotations with appropriate model Hamiltonians and fully understanding the origins of these motions, particularly the factors that influence torsional barriers, remain partially unresolved, reflecting the enduring mystery of quantum tunneling. This Perspective reviews the history of LAMs, highlights advances in decoding their complex spectra, and explores future research directions aimed at uncovering the remaining mysteries of these fascinating motions.

Exceptions, Paradoxes, and Their Resolutions in Chemical Reactivity

Progress in chemistry has primarily been framed through inductive processes, leading to the frequent emergence of exceptions and unexpected reactivities. These anomalies─ranging from surprising reactivity trends and paradoxical understandings to unanticipated parameter influences and unexpected successes or failures in synthetic methods─offer valuable insights that can drive scientific discovery. While it is commonly accepted that such exceptions can drive progress, many have been passively accepted without being explored for opportunities. Although numerous chemists have addressed exceptions and refined chemical models and theories, employing a systematic framework for actively exploring and understanding the underlying causes of exceptions could resolve paradoxes in broader contexts and create a greater impact than treating anomalies as isolated occurrences. This perspective demonstrates a proactive epistemic approach to uncovering the opportunities presented by exceptions and promotes deliberate, thoughtful engagement with paradoxes and anomalies. While the examples primarily focus on physical organic chemistry, the concepts are broadly applicable across various fields in chemical science. The thinking framework presented here is neither exhaustive nor prescriptive, but it outlines one of many potentially possible ways to inspire further development in how these anomalies could be harnessed for advancement.

Electrical monitoring of single-event protonation dynamics at the solid-liquid interface and its regulation by external mechanical forces

Detecting chemical reaction dynamics at solid-liquid interfaces is important for understanding heterogeneous reactions. However, there is a lack of exploration of interface reaction dynamics from the single-molecule perspective, which can reveal the intrinsic reaction mechanism underlying ensemble experiments. Here, single-event protonation reaction dynamics at a solid-liquid interface are studied in-situ using single-molecule junctions. Molecules with amino terminal groups are used to construct single-molecule junctions. An interfacial cationic state present after protonation is discovered. Real-time electrical measurements are used to monitor the reversible reaction between protonated and deprotonated states, thereby revealing the interfacial reaction mechanism through dynamic analysis. The protonation reaction rate constant has a linear positive correlation with proton concentration, whereas the deprotonation reaction rate constant has a linear negative correlation. In addition, external mechanical forces can effectively regulate the protonation reaction process. This work provides a single-molecule perspective for exploring interface science, which will contribute to the development of heterogeneous catalysis and electrochemistry.

Heavy-Atom Tunneling in Ring-Closure Reactions of Beryllium Ozonide Complexes

Quantum mechanical tunneling (QMT) has long been recognized as crucial for understanding chemical reaction mechanisms, particularly in reactions involving light atoms like hydrogen. However, recent findings have expanded this understanding to include heavy-atom tunneling reactions. In this report, we present the observation of two heavy-atom tunneling reactions involving the spontaneous conversions from end-on bonded beryllium ozonide complexes, OBeOOO (A) and BeOBeOOO (C), to their corresponding side-on bonded ozonide isomers, OBe(η2-O3) (B) and BeOBe(η2-O3) (D), respectively, in a cryogenic neon matrix. This discovery is supported by the weak temperature dependence of the rate constants and unusually large 16O/18O kinetic isotope effects. Quantum chemistry calculations reveal extremely low barriers (<1 kcal/mol) for both ring-closure reactions. Additionally, instanton theory calculations on both reactions unveil that the tunneling processes involve the concerted motion of all four oxygen atoms.

Hydrogen Tunneling Exhibiting Unexpectedly Small Primary Kinetic Isotope Effects

Probing quantum mechanical tunneling (QMT) in chemical reactions is crucial to understanding and developing new transformations. Primary H/D kinetic isotopic effects (KIEs) beyond the semiclassical maximum values of 7-10 (room temperature) are commonly used to assess substantial QMT contributions in one-step hydrogen transfer reactions, because of the much greater QMT probability of protium vs. deuterium. Nevertheless, we report here the discovery of a reaction model occurring exclusively by H-atom QMT with residual primary H/D KIEs. 2-Hydroxyphenylnitrene, generated in N2 matrix, was found to isomerize to an imino-ketone via sequential (domino) QMT involving anti to syn OH-rotamerization (rate determining step) and [1,4]-H shift reactions. These sequential QMT transformations were also observed in the OD-deuterated sample, and unexpected primary H/D KIEs between 3 and 4 were measured at 3 to 20 K. Analogous residual primary H/D KIEs were found in the anti to syn OH-rotamerization QMT of 2-cyanophenol in a N2 matrix. Evidence strongly indicates that these intriguing isotope-insensitive QMT reactivities arise due to the solvation effects of the N2 matrix medium, putatively through coupling with the moving H/D tunneling particle. Should a similar scenario be extrapolated to conventional solution conditions, then QMT may have been overlooked in many chemical reactions.

Unraveling Water Solvation Effects with Quantum Mechanics/Molecular Mechanics Semiclassical Vibrational Spectroscopy: The Case of Thymidine

We introduce a quantum mechanics/molecular mechanics semiclassical method for studying the solvation process of molecules in water at the nuclear quantum mechanical level with atomistic detail. We employ it in vibrational spectroscopy calculations because this is a tool that is very sensitive to the molecular environment. Specifically, we look at the vibrational spectroscopy of thymidine in liquid water. We find that the C═O frequency red shift and the C═C frequency blue shift, experienced by thymidyne upon solvation, are mainly due to reciprocal polarization effects, that the molecule and the water solvent exert on each other, and nuclear zero-point energy effects. In general, this work provides an accurate and practical tool to study quantum vibrational spectroscopy in solution and condensed phase, incorporating high-level and computationally affordable descriptions of both electronic and nuclear problems.

Nuclear Quantum Effects in Proton or Hydrogen Transfer

Proton or hydrogen transfers, basic chemical reactions, proceed either by thermally activated barrier crossing or via tunneling. Studies of molecules undergoing single or double proton or hydrogen transfer in the ground or excited electronic state reveal that tunneling can dominate under conditions usually considered to favor the thermal process. Moreover, the tunneling probability strongly varies for excitation of certain vibrational modes, which changes the effective barrier and/or proton transfer distance. When the reaction is fast compared to vibrational relaxation, the mode selectivity can still be maintained for molecules in solutions at 293 K. These observations point to dangers of relating the calculated minimum energy paths and the associated barriers to the experimentally obtained activation energies. The multidimensional character of the reaction coordinate is obvious; it can dramatically change for slowly and rapidly relaxing environments. We postulate that the hydrogen bond definition should be extended by specifically including the role of molecular vibrations.

Phosphinidenes: Fundamental Properties and Reactivity

Phosphinidenes are heavy congeners of nitrenes that have been broadly used as in situ reagents in synthetic phosphorus chemistry and also serve as versatile ligands in coordination with transition metals. However, the detection of free phosphinidenes is largely challenged by their high reactivity and also the lack of suitable synthetic methods, rendering the knowledge about the fundamental properties of this class of low-valent phosphorus compounds limited. Recently, an increasing number of free phosphinidenes bearing prototype structural and bonding properties have been prepared for the first time, thus enabling the exploration of their distinct reactivity from the nitrene analogues. This Concept article will discuss the experimental approaches for the generation of the highly unstable phosphinidenes and highlight their distinct reactivity from the nitrogen analogues so as to stimuate future studies about their potential applications in phosphorus chemistry.

The Intrinsic Barrier Width and Its Role in Chemical Reactivity

Chemical reactions are in virtually all cases understood and explained on the basis of depicting the molecular potential energy landscape, i.e., the change in atomic positions vs the free-energy change. With such landscapes, the features of the reaction barriers solely determine chemical reactivities. The Marcus dissection of the barrier height (activation energy) on such a potential into the thermodynamically independent (intrinsic) and the thermodynamically dependent (Bell-Evans-Polanyi) contributions successfully models the interplay of reaction rate and driving force. This has led to the well-known and ubiquitously used reactivity paradigm of "kinetic versus thermodynamic control". However, an analogous dissection concept regarding the barrier width is absent. Here we define and outline the concept of intrinsic barrier width and the driving force effect on the barrier width and report experimental as well as theoretical studies to demonstrate their distinct roles. We present the idea of changing the barrier widths of conformational isomerizations of some simple aromatic carboxylic acids as models and use quantum mechanical tunneling (QMT) half-lives as a read-out for these changes because QMT is particularly sensitive to barrier widths. We demonstrate the distinct roles of the intrinsic and the thermodynamic contributions of the barrier width on QMT half-lives. This sheds light on resolving conflicting trends in chemical reactivities where barrier widths are relevant and allows us to draw some important conclusions about the general relevance of barrier widths, their qualitative definition, and the consequences for more complete descriptions of chemical reactions.

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.

First-Principles Study of the Dispersion Effects in the Structures and Keto-Enol Tautomerization of Curcumin

Noncovalent interactions, such as dispersion, play a significant role in the stability of flexible molecules, such as curcumin. This study revealed the importance of dispersion correction in the structure and keto-enol tautomerization of curcumin, which has rarely been addressed in computational studies. We rigorously constructed all possible unique curcumin conformers in the enol and keto forms within the first-principles framework. Regardless of the different environments, we carefully explained the agreement between the computational geometry (in the gas phase) and the experimental measurement (in the polymorph) by using dispersion correction. The calculation results for the aqueous solution of conformational abundance, thermochemistry, and reaction kinetics support the experimental observations after considering the dispersion correction. The study also suggests a water-catalyzed mechanism for keto-enol tautomerization, where dispersion correction plays a role in decreasing the energy barrier and making the keto form thermochemically and kinetically favorable. Our results could be helpful in future computational studies to find a method for increasing the aqueous solubility of curcumin; hence, the potential of curcumin as a multifunctional medicine can be fully achieved.

The Marcus dimension: identifying the nuclear coordinate for electron transfer from ab initio calculations

The Marcus model forms the foundation for all modern discussion of electron transfer (ET). In this model, ET results in a change in diabatic potential energy surfaces, separated along an ET nuclear coordinate. This coordinate accounts for all nuclear motion that promotes electron transfer. It is usually assumed to be dominated by a collective asymmetric vibrational motion of the redox sites involved in the ET. However, this coordinate is rarely quantitatively specified. Instead, it remains a nebulous concept, rather than a tool for gaining true insight into the ET pathway. Herein, we describe an ab initio approach for quantifying the ET coordinate and demonstrate it for a series of dinitroradical anions. Using sampling methods at finite temperature combined with density functional theory calculations, we find that the electron transfer can be followed using the energy separation between potential energy surfaces and the extent of electron localization. The precise nuclear motion that leads to electron transfer is then obtained as a linear combination of normal modes. Once the coordinate is identified, we find that evolution along it results in a change in diabatic state and optical excitation energy, as predicted by the Marcus model. Thus, we conclude that a single dimension of the electron transfer described in Marcus-Hush theory can be described as a well-defined nuclear motion. Importantly, our approach allows the separation of the intrinsic electron transfer coordinate from other structural relaxations and environmental influences. Furthermore, the barrier separating the adiabatic minima was found to be sufficiently thin to enable heavy-atom tunneling in the ET process.

Valve turning towards on-cycle in cobalt-catalyzed Negishi-type cross-coupling

Ligands and additives are often utilized to stabilize low-valent catalytic metal species experimentally, while their role in suppressing metal deposition has been less studied. Herein, an on-cycle mechanism is reported for CoCl2bpy2 catalyzed Negishi-type cross-coupling. A full catalytic cycle of this kind of reaction was elucidated by multiple spectroscopic studies. The solvent and ligand were found to be essential for the generation of catalytic active Co(I) species, among which acetonitrile and bipyridine ligand are resistant to the disproportionation events of Co(I). Investigations, based on Quick-X-Ray Absorption Fine Structure (Q-XAFS) spectroscopy, Electron Paramagnetic Resonance (EPR), IR allied with DFT calculations, allow comprehensive mechanistic insights that establish the structural information of the catalytic active cobalt species along with the whole catalytic Co(I)/Co(III) cycle. Moreover, the acetonitrile and bipyridine system can be further extended to the acylation, allylation, and benzylation of aryl zinc reagents, which present a broad substrate scope with a catalytic amount of Co salt. Overall, this work provides a basic mechanistic perspective for designing cobalt-catalyzed cross-coupling reactions.

Hydrogen Bonding Shuts Down Tunneling in Hydroxycarbenes: A Gas-Phase Study by Tandem-Mass Spectrometry, Infrared Ion Spectroscopy, and Theory

Hydroxycarbenes can be generated and structurally characterized in the gas phase by collision-induced decarboxylation of α-keto carboxylic acids, followed by infrared ion spectroscopy. Using this approach, we have shown earlier that quantum-mechanical hydrogen tunneling (QMHT) accounts for the isomerization of a charge-tagged phenylhydroxycarbene to the corresponding aldehyde in the gas phase and above room temperature. Herein, we report the results of our current study on aliphatic trialkylammonio-tagged systems. Quite unexpectedly, the flexible 3-(trimethylammonio)propylhydroxycarbene turned out to be stable─no H-shift to either aldehyde or enol occurred. As supported by density functional theory calculations, this novel QMHT inhibition is due to intramolecular H-bonding of a mildly acidic α-ammonio C-H bonds to the hydroxyl carbene's C-atom (C:···H-C). To further support this hypothesis, (4-quinuclidinyl)hydroxycarbenes were synthesized, whose rigid structure prevents this intramolecular H-bonding. The latter hydroxycarbenes underwent "regular" QMHT to the aldehyde at rates comparable to, e.g., methylhydroxycarbene studied by Schreiner et al. While QMHT has been shown for a number of biological H-shift processes, its inhibition by H-bonding disclosed here may serve for the stabilization of highly reactive intermediates such as carbenes, even as a mechanism for biasing intrinsic selectivity patterns.

A highly accurate full-dimensional ab initio potential surface for the rearrangement of methylhydroxycarbene (H3C-C-OH)

We report here a full-dimensional machine learning global potential surface (PES) for the rearrangement of methylhydroxycarbene (H₃C-C-OH, 1t). The PES is trained with the fundamental invariant neural network (FI-NN) method on 91 564 ab initio energies calculated at the UCCSD(T)-F12a/cc-pVTZ level of theory, covering three possible product channels. FI-NN PES has the correct symmetry properties with respect to permutation of four identical hydrogen atoms and is suitable for dynamics studies of the 1t rearrangement. The averaged root mean square error (RMSE) is 11.4 meV. Six important reaction pathways, as well as the energies and vibrational frequencies at the stationary geometries on these pathways are accurately preproduced by our FI-NN PES. To demonstrate the capacity of the PES, we calculated the rate coefficient of hydrogen migration in -CH₃ (path A) and hydrogen migration of -OH (path B) with instanton theory on this PES. Our calculations predicted the half-life of 1t to be 95 min, which is excellent in agreement with experimental observations.

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.

Infrared Spectra and Phototransformations of meta-Fluorophenol Isolated in Argon and Nitrogen Matrices

Monomers of meta-fluorophenol (mFP) were trapped from the gas phase into cryogenic argon and nitrogen matrices. The estimated relative energies of the two conformers are very close, and in the gas phase they have nearly equal populations. Due to the similarity of their structures (they only differ in the orientation of the OH group), the two conformers have also similar predicted vibrational signatures, which makes the vibrational characterization of the individual rotamers challenging. In the present work, it has been established that in an argon matrix only the most stable trans conformer of mFP exists (the OH group pointing away from the fluorine atom). On the other hand, the IR spectrum of mFP in a nitrogen matrix testifies to the simultaneous presence in this matrix of both the trans conformer and of the higher-energy cis conformer (the OH group pointing toward the fluorine atom), which is stabilized by interaction with the matrix gas host. We found that the exposition of the cryogenic N₂ matrix to the Globar source of the infrared spectrometer affects the conformational populations. By collecting experimental spectra, either in the full mid-infrared range or only in the range below 2200 cm^-1, we were able to reliably distinguish two sets of experimental bands originating from individual conformers. A comparison of the two sets of experimental bands with computed infrared spectra of the conformers allowed, for the first time, the unequivocal vibrational identification of each of them. The joint implementation of computational vibrational spectroscopy and matrix-isolation infrared spectroscopy proved to be a very accurate method of structural analysis. Some mechanistic insights into conformational isomerism (the quantum tunneling of hydrogen atom and vibrationally-induced conformational transformations) have been addressed. Finally, we also subjected matrix-isolated mFP to irradiations with UV light, and the phototransformations observed in these experiments are also described.

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.