Tunneling Control of Chemical Reactions: The Third Reactivity Paradigm

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.

Hydrogen tunneling with an atypically small KIE measured in the mediated decomposition of the Co(CH3COOH)+ complex

Quantum mechanical tunneling (QMT) is a well-documented phenomenon in the C-H bond activation mechanism and is commonly identified by large KIE values. Herein we present surprising findings in the kinetic study of hydrogen tunneling in the Co+ mediated decomposition of acetic acid and its perdeuterated isotopologue, conducted with the energy resolved single photon initiated dissociative rearrangement reaction (SPIDRR) technique. Following laser activation, the reaction proceeds along parallel product channels Co(CH4O)+ + CO and Co(C2H2O)+ + H2O. An energetic threshold is observed in the energy dependence of the unimolecular microcanonical rate constants, k(E). This is interpreted as the reacting population surmounting a rate-limiting Eyring barrier in the reaction's potential energy surface. Measurements of the heavier isotopologue's reaction kinetics supports this interpretation. Kinetic signatures measured at energies below the Eyring barrier are attributed to H/D QMT. The below-the-barrier tunneling kinetics presents an unusually linear energy dependence and a staggeringly small tunneling KIE of ∼1.4 over a wide energy range. We explain this surprising observation in terms of a narrow tunneling barrier, wherein the electronic structure of the Co+ metal plays a pivotal role in enhanced reactivity by promoting efficient tunneling. These results suggest that hydrogen tunneling could play important functions in transition metal chemistry, such as that found in enzymatic mechanisms, even if small KIE values are measured.

Quantum tunnelling effect in the cis-trans isomerization of uranyl tetrahydroxide

The role of quantum tunnelling (QT) in the proton transfer kinetics of the uranyl tetrahydroxide (UTH, [UO2(OH)4]2-) cis to trans isomerization was computationally studied under three possible reaction pathways. The first pathway involved a direct proton transfer from the hydroxide ligand to the oxo atom. In the other two pathways, one or two water molecules were added to the second sphere. The first H2O, bound by hydrogen bonds to the ligands, acts as a bridge enabling a proton shuttling, a concerted hopping of a proton from the hydroxide to the oxo atom similar to the Grotthuss mechanism. In the third pathway, the second water molecule does not participate in the H-transfer chain, but works as an anchor for the first water molecule, limiting its movement and therefore enhancing the QT. Since experimentally the reaction occurs in water, the first two pathways (no water or one H2O) serve only as models of the gas phase behaviour, while the third pathway will always be thermodynamically and kinetically preferred. The effects were investigated in the gas phase as well as in a continuum aqueous model, including the H/D Kinetic Isotope Effect (KIE). The results indicate that at very low temperatures, QT is the only mechanism that permits the reaction kinetics, consistent with the large computed KIE. At higher temperatures, thermally activated tunnelling competes with the classical crossing over the potential barrier.

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.

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.

Phenyl Radical Activates Molecular Hydrogen Through Protium and Deuterium Tunneling

Activating dihydrogen, H2, is a challenging endeavor typically achieved using transition metal centers. Pure main-group compounds capable of this are rare and have emerged in recent decades. These systems rely on synergistic donor-acceptor interactions with H2's antibonding σ* and bonding σ orbital. An alternative (hydrocarbon) radical-mediated activation is problematic because the H-H bond is stronger (104.2 kcal mol-1) than most C-H bonds. Here, we explore using the phenyl radical to tackle this, forming benzene with a C-H bond energy (112.9 kcal mol-1) that provides a thermodynamic driving force. We mainly observe a benzene-HI complex upon photolysis of iodobenzene in an H2-doped neon matrix at 4.4 K despite a barrier of 7.6 kcal mol-1, while phenyl radical forms in case of the heavier D2 isotopologue. When D2 molecules are allowed to diffuse, mono-deuterated benzene accumulates within hours. Computations using path integral-based instanton theory highlight that primarily the transferred hydrogen atom is moving during the reaction which greatly increases the tunneling probability. In excellent agreement with the experimental results, we predict significant tunneling rate constants for both isotopologues, H2 and D2, featuring a strong kinetic isotope effect of up to four orders of magnitude at the lowest temperatures.

Nuclear Quantum Effects on the Nature of Hydroboration Selectivity: Experimental Effects of First-Collision Tunneling

The understanding of selectivity in reactions exhibiting nonstatistical dynamics is impeded by the limitations of trajectory studies with regard to nuclear quantum effects, especially tunneling. We described here the use of ring-polymer molecular dynamics (RPMD) to account for an unusual regiochemical isotope effect on the regioselectivity of hydroborations of alkenes with BH3/BD3. RPMD is able to account for the experimental observation, while statistical approaches and classical trajectories fail. The combination of experiment and RPMD trajectories suggests that tunneling in the initial collision of reactants is the major source of the nonstatistical selectivity.

Intramolecular hydrogen-atom tunneling in matrix-isolated heterocyclic compounds: 2-thiouracil and its analogues

Hydrogen-atom tunneling leading to spontaneous tautomeric conversion in monomeric heterocyclic molecules without intramolecular hydrogen bonds has been experimentally detected for the first time. For monomers of 2-thiouracil, 6-aza-2-thiouracil and 1-methyl-2-thiouracil isolated in low-temperature matrices, higher-energy thiol forms were generated upon UV (λ = 305 nm) excitation of the most stable thione tautomers. When the matrices were subsequently kept in the dark and at low temperature, hydrogen-atom tunneling occurred, leading to the thiol → thione conversion. During this process, the photoproduced thiol form spontaneously converted into the lowest energy thione tautomer.

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.

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.

Norcaradiene-Cycloheptatriene Equilibrium: A Heavy-Atom Quantum Tunneling Case

The equilibrium between norcaradiene and cycloheptatriene, which has captivated chemists for more than half a century, is revisited by state-of-the-art quantum chemical calculations. Our theoretical data significantly deviate from the experimental results (J. Am. Chem. Soc., 1981, 26, 7791-7792), especially at low temperatures, where isomerization is dominated by heavy-atom tunneling. This effect results in an extremely short half-life for norcaradiene, rendering it undetectable. This work sheds light on this equilibrium, updating the kinetic and thermodynamic data while also expanding the repertoire of organic reactions controlled by this exotic quantum effect.

Quantum Tunneling Instability in Pericyclic Reactions: The Cheletropic, Coarctate, and Ene Cases

Some retro-pericyclic reactions, as a result of their high exothermicity and short trajectories, are the perfect ground for heavy atom tunneling molecular decompositions, also known as "quantum tunneling instability" (QTI). Considering this effect, in our first installment [Frenklach, A.; Amlani, H.; Kozuch, S. Quantum Tunneling Instability in Pericyclic Reactions. J. Am. Chem. Soc. 2024, 146 (17), 11823-11834, DOI: 10.1021/jacs.4c00608], we computed several retro-Diels-Alder reactions, predicting that many studied reactants cannot be isolated. Herein, we will explore the QTI of retro-cheletropic, coarctate, and ene exemplars, where again we hypothesize the impossibility to detect their reactants.

Quantum Tunneling Instability in Pericyclic Reactions

Several cycloreversion reactions of the retro-Diels-Alder type were computationally assessed to understand their quantum tunneling (QT) reactivity. N2, CO, and other leaving groups were considered based on their strong exothermicity, as it reduces their thermodynamic and kinetic stabilities. Our results indicate that for many of these reactions, it is essential to take into account their QT decomposition rate, which can massively weaken their molecular stability and shorten their half-lives even at deep cryogenic temperatures. In practical terms, this indicates that many supposedly stable molecules will actually be unsynthesizable or unisolable, and therefore trying to prepare or detect them would be a futile attempt. In addition, we discuss the importance of tunneling to correctly understand the enthalpy of activation and the collective atomic effect on the tunneling kinetic isotope effects to test if third-row atoms can tunnel in a chemical reaction. This project raises the question of the importance of in silico chemistry to guide in vitro chemistry, especially in cases where the latter cannot solve its own uncertainties.

Optimization of Quantum Nuclei Positions with the Adaptive Nuclear-Electronic Orbital Approach

The use of multicomponent methods has become increasingly popular over the last years. Under this framework, nuclei (commonly protons) are treated quantum mechanically on the same footing as the electronic structure problem. Under the use of atomic-centered orbitals, this can lead to some complications as the ideal location of the nuclear basis centers must be optimized. In this contribution, we propose a straightforward approach to determine the position of such centers within the self-consistent cycle of a multicomponent calculation, making use of individual proton charge centroids. We test the method on model systems including the water dimer, a protonated water tetramer, and a porphine system. Comparing to numerical gradient calculations, the adaptive nuclear-electronic orbital (NEO) procedure is able to converge the basis centers to within a few cents of an Ångström and with less than 0.1 kcal/mol differences in absolute energies. This is achieved in one single calculation and with a small added computational effort of up to 80% compared to a regular NEO- self-consistent field run. An example application for the human transketolase proton wire is also provided.

Proton tunneling in the dissociation of H2+ and its asymmetric isotopologues driven by circularly polarized THz laser pulses

Light-induced deprotonation of molecules is an important process in photochemical reactions. Here, we theoretically investigate the tunneling deprotonation of H2+ and its asymmetric isotopologues driven by circularly polarized THz laser pulses. The quasi-static picture shows that the field-dressed potential barrier is significantly lowered for the deprotonation channel when the mass asymmetry of the diatomic molecule increases. Our numerical simulations demonstrate that when the mass symmetry breaks, the tunneling deprotonation is significantly enhanced and the proton tunneling becomes the dominant dissociation channel in the THz driving fields. In addition, the simulated nuclear momentum distributions show that the emission of the proton is directed by the effective vector potential for the deprotonation channel and, meanwhile, the angular distribution of the emitting proton is affected by the alignment and rotation of the molecule induced by the rotating field.

Infrared spectroscopy of the α-hydroxyethyl radical isolated in cryogenic solid media

The α-hydroxyethyl radical (CH3·CHOH, 2A) is a key intermediate in ethanol biochemistry, combustion, atmospheric chemistry, radiation chemistry, and astrochemistry. Experimental data on the vibrational spectrum of this radical are crucially important for reliable detection and understanding of the chemical dynamics of this species. This study represents the first detailed experimental report on the infrared absorption bands of the α-hydroxyethyl radical complemented by ab initio computations. The radical was generated in solid para-H2 and Xe matrices via the reactions of hydrogen atoms with matrix-isolated ethanol molecules and radiolysis of isolated ethanol molecules with x rays. The absorption bands with maxima at 3654.6, 3052.1, 1425.7, 1247.9, 1195.6 (1177.4), and 1048.4 cm-1, observed in para-H2 matrices appearing upon the H· atom reaction, were attributed to the OHstr, α-CHstr, CCstr, COstr + CCObend, COstr, and CCstr + CCObend vibrational modes of the CH3·CHOH radical, respectively. The absorption bands with the positions slightly red-shifted from those observed in para-H2 were detected in both the irradiated and post-irradiation annealed Xe matrices containing C2H5OH. The results of the experiments with the isotopically substituted ethanol molecules (CH3CD2OH and CD3CD2OH) and the quantum-chemical computations at the UCCSD(T)/L2a_3 level support the assignment. The photolysis with ultraviolet light (240-300 nm) results in the decay of the α-hydroxyethyl radical, yielding acetaldehyde and its isomer, vinyl alcohol. A comparison of the experimental and theoretical results suggests that the radical adopts the thermodynamically more stable anti-conformation in both matrices.

Heme-copper and Heme O2-derived synthetic (bioinorganic) chemistry toward an understanding of cytochrome c oxidase dioxygen chemistry

Cytochrome c oxidase (CcO), also widely known as mitochondrial electron-transport-chain complex IV, is a multi-subunit transmembrane protein responsible for catalyzing the last step of the electron transport chain, dioxygen reduction to water, which is essential to the establishment and maintenance of the membrane proton gradient that drives ATP synthesis. Although many intermediates in the CcO catalytic cycle have been spectroscopically and/or computationally authenticated, the specifics regarding the IP intermediate, hypothesized to be a heme-Cu (hydro)peroxo species whose O-O bond homolysis is supported by a hydrogen-bonding network of water molecules, are largely obscured by the fast kinetics of the A (FeIII-O2•-/CuI/Tyr) → PM (FeIV=O/CuII-OH/Tyr•) step. In this review, we have focused on the recent advancements in the design, development, and characterization of synthetic heme-peroxo‑copper model complexes, which can circumvent the abovementioned limitation, for the investigation of the formation of IP and its O-O cleavage chemistry. Novel findings regarding (a) proton and electron transfer (PT/ET) processes, together with their contributions to exogenous phenol induced O-O cleavage, (b) the stereo-electronic tunability of the secondary coordination sphere (especially hydrogen-bonding) on the geometric and spin state alteration of the heme-peroxo‑copper unit, and (c) a plausible mechanism for the Tyr-His cofactor biogenesis, are discussed in great detail. Additionally, since the ferric-superoxide and the ferryl-oxo (Compound II) species are critically involved in the CcO catalytic cycle, this review also highlights a few fundamental aspects of these heme-only (i.e., without copper) species, including the structural and reactivity influences of electron-donating trans-axial ligands and Lewis acid-promoted H-bonding.

Rapid Thalidomide Racemization Is Related to Proton Tunneling Reactions via Water Bridges

Quantum mechanical tunneling (QMT) can play an important role in light element-related chemical reactions; however, its influence on racemization is not fully understood. Herein, we demonstrate that the role of QMT is decisive for rapid racemization of the well-known thalidomide molecule in aqueous environments, increasing the reaction rate constants of the most likely racemization pathways by 87-149 times at approximately body temperature and achieving good agreement between theoretical calculations and experimental observations. In addition, the kinetic isotope effect values fit well with those of previous experiments. These results are attributed to enhanced tunneling probability due to the alteration of potential barriers for proton transfer reactions via water bridges. This work highlights the significance of the QMT effect in racemization and its potential impact on drug safety, providing a fundamental perspective for understanding chirality-related issues in biological systems.

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.

The enol of propionic acid

We demonstrate the gas-phase synthesis of prop-1-ene-1,1-diol, the hitherto unreported higher energy tautomer of propionic acid. The enol was trapped in an argon matrix and characterized by IR and UV/Vis spectroscopy in combination with density functional theory computations. Upon photolysis, the enol rearranges to propionic acid and methylketene.

Hydrogen-Atom Tunneling in a Homochiral Environment

The role-exchanging concerted torsional motion of two hydrogen atoms in the homochiral dimer of trans-1,2-cyclohexanediol was characterized through a combination of broadband rotational spectroscopy and theoretical modeling. The results reveal that the concerted tunneling motion of the hydrogen atoms leads to the inversion of the sign of the dipole moment components along the a and b principal axes, due to the interchange motion that cooperatively breaks and reforms one intermolecular hydrogen bond. This motion is also coupled with two acceptor switching motions. The energy difference between the two ground vibrational states arising from this tunneling motion was determined to be 29.003(2) MHz. The corresponding wavefunctions suggest that the two hydrogen atoms are evenly delocalized on two equivalent potential wells, which differs from the heterochiral case where the hydrogen atoms are confined in separate wells, as the permutation-inversion symmetry breaks down. This intriguing contrast in hydrogen-atom behavior between homochiral and heterochiral environments could further illuminate our understanding of the role of chirality in intermolecular interactions and dynamics.

Modulating Escape Channels of Cycloheptatrienyl Rhodium Carbenes To Form Semibullvalene

We describe the various escape channels available to dirhodium carbene intermediates from cycloheptatrienyl diazo compounds located with density functional theory. An intramolecular cyclopropanation would, in principle, provide a new route to semibullvalenes (SBVs). A detailed exploration of the potential energy surface reveals that methylating carbon-7 suppresses a competing β-hydride migration pathway to heptafulvene products, giving SBV formation a reasonable chance. During our explorations, we additionally discovered unusual spirononatriene, spironorcaradiene, and metal-stabilized 9-barbaralyl cation structures as local minima.

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.

Accelerated Storage for Shelf-Life Prediction of Lyophiles: Temperature Dependence of Degradation of Amorphous Small Molecular Weight Drugs and Proteins

Prediction of lyophilized product shelf-life using accelerated stability data requires understanding the temperature dependence of the degradation rate. Despite the abundance of published studies on stability of freeze-dried formulations and other amorphous materials, there are no definitive conclusions on the type of pattern one can expect for the temperature dependence of degradation. This lack of consensus represents a significant gap which may impact development and regulatory acceptance of freeze-dried pharmaceuticals and biopharmaceuticals. Review of the literature demonstrates that the temperature dependence of degradation rate constants in lyophiles can be represented by the Arrhenius equation in most cases. In some instances there is a break in the Arrhenius plot around the glass transition temperature or a related characteristic temperature. The majority of the activation energies (E_a), which are reported for various degradation pathways in lyophiles, falls in the range of 8 to 25 kcal/mol. The degradation E_a values for lyophiles are compared with the E_a for relaxation processes and diffusion in glasses, as wells as solution chemical reactions. Collectively, analysis of the literature demonstrates that the Arrhenius equation represents a reasonable empirical tool for analysis, presentation, and extrapolation of stability data for lyophiles, provided that specific conditions are met.

Molecular Understanding and Practical In Silico Catalyst Design in Computational Organocatalysis and Phase Transfer Catalysis-Challenges and Opportunities

Through the lens of organocatalysis and phase transfer catalysis, we will examine the key components to calculate or predict catalysis-performance metrics, such as turnover frequency and measurement of stereoselectivity, via computational chemistry. The state-of-the-art tools available to calculate potential energy and, consequently, free energy, together with their caveats, will be discussed via examples from the literature. Through various examples from organocatalysis and phase transfer catalysis, we will highlight the challenges related to the mechanism, transition state theory, and solvation involved in translating calculated barriers to the turnover frequency or a metric of stereoselectivity. Examples in the literature that validated their theoretical models will be showcased. Lastly, the relevance and opportunity afforded by machine learning will be discussed.

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.

Prediction of Kinetic Product Ratios: Investigation of a Dynamically Controlled Case

Of the various factors influencing kinetically controlled product ratios, the role of nonstatistical dynamics is arguably the least well understood. In this paper, reactions were chosen in which dynamics played a dominant role in product selection, by design. Specifically, the reactions studied were the ring openings of cyclopropylidene to allene and tetramethylcyclopropylidene to tetramethylallene (2,4-dimethylpenta-2,3-diene). Both reactions have intrinsic reaction coordinates that bifurcate symmetrically, leading to products that are enantiomeric once the atoms are uniquely labeled. The question addressed in the study was whether the outcomes─that is, which product well on the potential energy surface was selected─could be predicted from their initial conditions for individual trajectories in quasiclassical dynamics simulations. Hybrid potentials were developed based on cooperative interaction between molecular mechanics and artificial neural networks, trained against data from electronic structure calculations. These potentials allowed simulations of both gas-phase and condensed-phase reactions. The outcome was that, for both reactions, prediction of initial selection of product wells could be made with >95% success from initial conditions of the trajectories in the gas phase. However, when trajectories were run for longer, looking for "final" products for each trajectory, the predictability dropped off dramatically. In the gas-phase simulations, this drop off was caused by trajectories hopping between product wells on the potential energy surface. That behavior could be suppressed in condensed phases, but then new uncertainty was introduced because the intermolecular interactions between solute and bath, necessary to permit intermolecular energy transfer and cooling of the hot initial products, often led to perturbations of the initial directions of trajectories on the potential energy surface. It would consequently appear that a general ability to predict outcomes for reactions in which nonstatistical dynamics dominate remains a challenge even in the age of sophisticated machine-learning capabilities.

Equilibrating parent aminomercaptocarbene and CO2 with 2-amino-2-thioxoacetic acid via heavy-atom quantum tunneling

The search for methods to bind CO2 and use it synthetically as a C1-building block under mild conditions is an ongoing endeavor of great urgency. The formation of heterocyclic carbene-carbon dioxide adducts occurs rapidly when the carbene is generated in solution in the presence of CO2. Here we demonstrate the reversible formation of a complex of the hitherto unreported aminomercaptocarbene (H2N-C̈-SH) with CO2 isolated in solid argon by photolysis of 2-amino-2-thioxoacetic acid. Remarkably, the complex disappears in the dark as deduced by time-dependent matrix infrared measurements, and equilibrates back to the covalently bound starting material. This kinetically excluded process below ca. 8 K is made possible through heavy-atom quantum mechanical tunneling, as also evident from density functional theory and ab initio computations at the CCSD(T)/cc-pVTZ level of theory. Our results provide insight into CO2 activation using a carbene and emphasize the role of quantum mechanical tunneling in organic processes, even involving heavy atoms.

Mactanamide and lariciresinol as radical scavengers and Fe2+ ion chelators - A DFT study

A DFT based kinetic study of OOH radical scavenging potency of mactanamide (MA) and lariciresinol (LA), two natural polyphenols, indicates their nearly equal potential via the proton coupled electron transfer (PCET) mechanism in lipid media. Contribution of C-H bond breaking to this potency is negligible compared to O-H bond breaking, contrary to recent claims. The predicted potency of both compounds is not sufficient to protect biological molecules from oxidative damage in lipid media. In aqueous media, the scavenging potency of MA and LA phenoxide anions via the single electron transfer (SET) mechanism is much higher and may contribute to the protection of lipids, proteins, and DNA from OOH radical damage. Also, MA and LA have the potential to chelate catalytic Fe2+ ions, thus suppressing the formation of dangerous OH radicals via Fenton-type reactions. The monoanionic species of MA and LA show stronger monodentate chelating ability with Fe2+ ion compared to its neutral form. The dianionic specie LA2- exhibited the highest chelation ability with Fe2+ ion via bidentate 1:2 coordination. However, direct radical scavenging and metal chelation could be rarely operative in vivo because MA and LA presumably achieve very low concentrations in systemic circulation.

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.

Unraveling sulfur chemistry in interstellar carbon oxide ices

Formyl radical (HCO•) and hydroxycarbonyl radical (HOCO•) are versatile building blocks in the formation of biorelevant complex organic molecules (COMs) in interstellar medium. Understanding the chemical pathways for the formation of HCO• and HOCO• starting with primordial substances (e.g., CO and CO2) is of vital importance in building the complex network of prebiotic chemistry. Here, we report the efficient formation of HCO• and HOCO• in the photochemistry of hydroxidooxidosulfur radical (HOSO•)–a key intermediate in SO2 photochemistry–in interstellar analogous ices of CO and CO2 at 16 K through hydrogen atom transfer (HAT) reactions. Specifically, 266 nm laser photolysis of HOSO• embedded in solid CO ice yields the elusive hydrogen‑bonded complexes HCO•···SO2 and HOCO•···SO, and the latter undergoes subsequent HAT to furnish CO2···HOS• under the irradiation conditions. Similar photo-induced HAT of HOSO• in solid CO2 ice leads to the formation of HOCO•···SO2. The HAT reactions of HOSO• in astronomical CO and CO2 ices by forming reactive acyl radicals may contribute to understanding the interplay between the sulfur and carbon ice-grain chemistry in cold molecular clouds and also in the planetary atmospheric chemistry.

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.

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.

Theoretical Investigation of Dissociation versus Intramolecular Rearrangements in Aminohydroxymethylene

Aminohydroxymethylene (H₂N-C̈-OH) is the simplest aminooxycarbene which is a heteroatom stabilized carbene. This highly reactive molecule was prepared in an Ar matrix in a recent experimental work. Unimolecular reactivity of this astrochemically important molecule was investigated and only fragmentations were identified contrary to the observations of both fragmentations and intramolecular rearrangements in other hydroxycarbenes. These rearrangement reactions form the corresponding imine and carbonyl compounds. In the present work, direct chemical dynamics simulations of unimolecular chemistry of aminohydroxymethylene were performed in the gas phase to study atomic level dissociation mechanisms. Classical trajectories were generated on-the-fly using potentials and gradients computed at the density functional B3LYP/6-31+G* level of electronic structure theory. Simulation results showed that intramolecular rearrangements accompany fragmentations during the unimolecular decay process of aminohydroxymethylene. However, the average lifetime of the intermediate isomers were found to be only few picoseconds which might not have been long enough for detection in the experiments.

Halogen-atom and group transfer reactivity enabled by hydrogen tunneling

The generation of carbon radicals by halogen-atom and group transfer reactions is generally achieved using tin and silicon reagents that maximize the interplay of enthalpic (thermodynamic) and polar (kinetic) effects. In this work, we demonstrate a distinct reactivity mode enabled by quantum mechanical tunneling that uses the cyclohexadiene derivative γ-terpinene as the abstractor under mild photochemical conditions. This protocol activates alkyl and aryl halides as well as several alcohol and thiol derivatives. Experimental and computational studies unveiled a noncanonical pathway whereby a cyclohexadienyl radical undergoes concerted aromatization and halogen-atom or group abstraction through the reactivity of an effective H atom. This activation mechanism is seemingly thermodynamically and kinetically unfavorable but is rendered feasible through quantum tunneling.

Antioxidant activity of caffeic acid: thermodynamic and kinetic aspects on the oxidative degradation pathway

Caffeic acid is a phenolic secondary metabolite from plants, which is known for its antioxidant properties. The effective mitigation of methanol-induced oxidative stress by caffeic acid depends on the direct radical scavenging as well as the formation of new metabolites via oxidative degradation. Herein, thermodynamic and kinetic aspects of the oxidative degradation pathway of caffeic acid in the presence of radical CH₃O^• and its isomer, ^•CH₂OH are discussed for the first time, employing density functional theory (DFT). The direct radical scavenging activity of caffeic acid against these radicals is verified via hydrogen atom transfer (HAT) and radical adduct formation (RAF) mechanisms. HAT is predicted to be more feasible than RAF mechanism as per the computed data. Additionally, energetic details of the proposed oxidative degradation pathway of radical adduct intermediates toward the formation of a cyclic metabolite is analyzed. Kinetic studies indicated a significant tunneling contribution to the H abstraction pathways having high activation barriers. Further, our results imply that the newly formed metabolites exhibit comparable antioxidant activity with that of caffeic acid.

Upgrading the peroxi-coagulation treatment of complex water matrices using a magnetically assembled mZVI/DSA anode: Insights into the importance of ClO radical

The electrochemical technologies for water treatment have flourished over the last decades. However, it is still challenging to treat the actual complex water effluents by a single electrochemical process, often requiring coupling of technologies. In this study, an upgraded peroxi-coagulation (PC) process with a magnetically assembled mZVI/DSA anode has been devised for the first time. COD, NH₃-N and total phosphorous were simultaneously and effectively removed from livestock wastewater. The advantages, influence of key parameters and evolution of electrogenerated species were systematically investigated to fully understand this novel PC process. The fluorescent substances in livestock wastewater could also be almost removed under optimal conditions (300 mA, 0.2 g ZVI particles and pH 6.8). The interaction between OH and active chlorine yielded ClO with a high steady-state concentration of 6.85 × 10^-13 M, which did not cause COD removal but accelerated the oxidation of NH₃-N. The Mulliken population suggested that OH and NH₃-N had similar electron-donor behavior, whereas ClO acted as an electron-withdrawing species. Besides, although the energy barrier for the reaction between OH and NH₃-N (17.0 kcal/mol) was lower than that with ClO (18.8 kcal/mol), considering the tunneling in the H abstraction reaction, the Skodje-Truhlar method adopted for calculations evidenced a 17-fold faster NH₃-N oxidation rate with ClO. In summary, this work describes an advantageous single electrochemical process for the effective treatment of a complex water matrix.

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.

Heavy-Atom Quantum Tunnelling in Spin Crossovers of Nitrenes

We simulate two recent matrix-isolation experiments at cryogenic temperatures, in which a nitrene undergoes spin crossover from its triplet state to a singlet state via quantum tunnelling. We detail the failure of the commonly applied weak-coupling method (based on a linear approximation of the potentials) in describing these deep-tunnelling reactions. The more rigorous approach of semiclassical golden-rule instanton theory in conjunction with double-hybrid density-functional theory and multireference perturbation theory does, however, provide rate constants and kinetic isotope effects in good agreement with experiment. In addition, these calculations locate the optimal tunnelling pathways, which provide a molecular picture of the reaction mechanism. The reactions involve substantial heavy-atom quantum tunnelling of carbon, nitrogen and oxygen atoms, which unexpectedly even continues to play a role at room temperature.

Photolytic insertion of carbon monoxide into nitrosyl chloride: formation of nitrosoformyl chloride

Nitrosocarbonyls are exotic intermediates that remain scarcely characterized. By UV photolysis (365 nm) of nitrosyl chloride (ClNO) embedded in solid CO ice at 20 K, the elusive nitrosoformyl chloride (ClC(O)NO) has been synthesized via CO-insertion into the Cl-N bond in ClNO. The characterization of ClC(O)NO with matrix-isolation IR spectroscopy is supported by ¹³C and ¹⁵N isotope labeling and quantum chemical calculations at the B3LYP/6-311+G(3df) level of theory. Upon subsequent laser irradiation at 266 nm, CO-elimination in ClC(O)NO occurs by reformation of ClNO. In line with the calculated potential energy surface for ClC(O)NO at the CCSD(T)-F12a/aug-cc-pVTZ//B3LYP/6-311+G(3df) level, the observed IR frequencies and the corresponding isotopic shifts coincide with the calculated values for the lowest-energy planar conformer, in which the CO and NO moities adopt trans configuration with respect to the C-N bond. Furthermore, the CO-insertion in ClNO involves a stepwise pathway by first homolytic cleavage of the Cl-N bond in ClNO (→ Cl˙ + ˙NO), followed by successive CO-trapping (CO + Cl˙ → ClCO˙) and radical combination (ClCO˙ + ˙NO → ClC(O)NO) inside the solid CO-matrix cages.

Iterative Synthesis of 2-Deoxyoligosaccharides Enabled by Stereoselective Visible-Light-Promoted Glycosylation

The photoinitiated intramolecular hydroetherification of alkenols has been used to form C-O bonds, but the intermolecular hydroetherification of alkenes with alcohols remains an unsolved challenge. We herein report the visible-light-promoted 2-deoxyglycosylation of alcohols with glycals. The glycosylation reaction was completed within 2 min in a high quantum yield (ϕ=28.6). This method was suitable for a wide array of substrates and displayed good reaction yields and excellent stereoselectivity. The value of this protocol was further demonstrated by the iterative synthesis of 2-deoxyglycans with α-2-deoxyglycosidic linkages up to a 20-mer in length and digoxin with β-2-deoxyglycosidic linkages. Mechanistic studies indicated that this reaction involved a glycosyl radical cation intermediate and a photoinitiated chain process.

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.

New Insights and Predictions into Complex Homogeneous Reactions Enabled by Computational Chemistry in Synergy with Experiments: Isotopes and Mechanisms

Homogeneous catalysis and biocatalysis have been widely applied in synthetic, medicinal, and energy chemistry as well as synthetic biology. Driven by developments of new computational chemistry methods and better computer hardware, computational chemistry has become an essentially indispensable mechanistic "instrument" to help understand structures and decipher reaction mechanisms in catalysis. In addition, synergy between computational and experimental chemistry deepens our mechanistic understanding, which further promotes the rational design of new catalysts. In this Account, we summarize new or deeper mechanistic insights (including isotope, dispersion, and dynamical effects) into several complex homogeneous reactions from our systematic computational studies along with subsequent experimental studies by different groups. Apart from uncovering new mechanisms in some reactions, a few computational predictions (such as excited-state heavy-atom tunneling, steric-controlled enantioswitching, and a new geminal addition mechanism) based on our mechanistic insights were further verified by ensuing experiments.The Zimmerman group developed a photoinduced triplet di-π-methane rearrangement to form cyclopropane derivatives. Recently, our computational study predicted the first excited-state heavy-atom (carbon) quantum tunneling in one triplet di-π-methane rearrangement, in which the reaction rates and ¹²C/¹³C kinetic isotope effects (KIEs) can be enhanced by quantum tunneling at low temperatures. This unprecedented excited-state heavy-atom tunneling in a photoinduced reaction has recently been verified by an experimental ¹²C/¹³C KIE study by the Singleton group. Such combined computational and experimental studies should open up opportunities to discover more rare excited-state heavy-atom tunneling in other photoinduced reactions. In addition, we found unexpectedly large secondary KIE values in the five-coordinate Fe(III)-catalyzed hetero-Diels-Alder pathway, even with substantial C-C bond formation, due to the non-negligible equilibrium isotope effect (EIE) derived from altered metal coordination. Therefore, these KIE values cannot reliably reflect transition-state structures for the five-coordinate metal pathway. Furthermore, our density functional theory (DFT) quasi-classical molecular dynamics (MD) simulations demonstrated that the coordination mode and/or spin state of the iron metal as well as an electric field can affect the dynamics of this reaction (e.g., the dynamically stepwise process, the entrance/exit reaction channels).Moreover, we unveiled a new reaction mechanism to account for the uncommon Ru(II)-catalyzed geminal-addition semihydrogenation and hydroboration of silyl alkynes. Our proposed key gem-Ru(II)-carbene intermediates derived from double migrations on the same alkyne carbon were verified by crossover experiments. Additionally, our DFT MD simulations suggested that the first hydrogen migration transition-state structures may directly and quickly form the key gem-Ru-carbene structures, thereby "bypassing" the second migration step. Furthermore, our extensive study revealed the origin of the enantioselectivity of the Cu(I)-catalyzed 1,3-dipolar cycloaddition of azomethine ylides with β-substituted alkenyl bicyclic heteroarenes enabled by dual coordination of both substrates. Such mechanistic insights promoted our computational predictions of the enantioselectivity reversal for the corresponding monocyclic heteroarene substrates and the regiospecific addition to the less reactive internal C═C bond of one diene substrate. These predictions were proven by our experimental collaborators. Finally, our mechanistic insights into a few other reactions are also presented. Overall, we hope that these interactive computational and experimental studies enrich our mechanistic understanding and aid in reaction development.

Hydroxy Mercapto Methylene: The Missing H2CSO Isomer

We report the isolation of hydroxy mercapto methylene (HO-C̈-SH) under cryogenic conditions via pyrolysis of 2-ethoxy-2-thioxo-acetic acid. The two most stable carbene rotamers form via extrusion of ethylene and CO₂ from this precursor. This donor-stabilized carbene represents a hitherto uncharacterized CH₂SO species and the first spectroscopically characterized free mercapto carbene. CCSD(T)/cc-pVTZ computations support our findings.

IR-induced and tunneling reactions in cryogenic matrices: the (incomplete) story of a successful endeavor

In this article, IR-induced and tunneling-driven reactions observed in cryogenic matrices are described in a historical perspective, the entangling of the two types of processes being highlighted. The story of this still ongoing fascinating scientific endeavor is presented here following closely our own involvement in the field for more than 30 years, and thus focuses mostly on our work. It is, because of this reason, also an incomplete story. Nevertheless, it considers a large range of examples, from very selective IR-induced conformational isomerizations to IR-induced bond-breaking/bond-forming reactions and successful observations of rare heavy atom tunneling processes. As a whole, this article provides a rather general overview of the major progress achieved in the field.