Direct Proton Transfer Research Articles

Biocatalytic Carbenoid N-H Insertions Via Engineered Myoglobins: A Systematic Reaction Mechanism Study

Recently engineered myoglobins have demonstrated excellent biocatalytic activities for a broad range of carbenoid N-H insertion reactions with up to >99% yield. Such biocatalysts offer various tuning points beyond protein residues, such as metal center, and axial ligand, to regulate different chemical reactivities. However, there is no computational mechanistic study to reveal its basic reaction mechanism, which could be utilized to study effects of different catalyst component for rational biocatalyst design. The reported computational mechanistic papers of other heme protein catalyzed N-H insertion have limited information without a careful comparison of all possible pathways and just focused on certain steps. Building on our group’s recent mechanistic work on a number of heme carbene transfer reactions, we performed a quantum chemical study to systematically study all possible reaction pathways for myoglobin catalyzed N-H insertions, which starts from the reactants to final released products and thus constitute a complete reaction pathway study. We examined ylide (proton transfer), hydrogen atom transfer, and hydride transfer pathways. In the favored ylide pathways, we calculated the direct proton transfer from amine to carbene and indirect paths with both concerted and stepwise proton transfer and dissociation components. For the indirect pathways in which the proton resides on the carbonyl oxygen of the original carbene substituent, subsequent water-assisted proton transfer from this carbonyl oxygen to carbene’s carbon to form the final product was also investigated. The most favorable pathway from this systematic reaction mechanism study was further supported by new experimental work from our collaborators, the original developer of myoglobin-based biocatalysts. These novel mechanistic results provide important mechanistic information for future biocatalyst design with enhanced reactivities.

Proton Transfer Kinetics at a Defect-Free Liquid-Liquid Interface

Proton transfer at electrochemical interfaces is fundamentally important for many areas of science and technology, yet kinetic measurements of this elementary step are often convoluted by inhomogeneous electrode surface structures. We show that facilitated proton transfer at the interface between two immiscible electrolyte solutions (ITIES) can serve as a model system to study proton transfer kinetics in the absence of defects found at solid|electrolyte interfaces. Diffusion-controlled micropipette voltammetry revealed that 2,6-diphenylpyridine (DPP) facilitated direct proton transfer across the liquid|liquid interface and voltammetry at nanopipette-supported interfaces yielded activation-controlled ion transfer currents. The exchange current density for proton transfer at the ITIES is comparable to exchange current densities recently reported for H-UPD on Pt (111) and provides the first direct means of comparing ion transfer kinetics between solid|liquid and liquid|liquid interfaces, supporting the development of a general theory of ion transfer kinetics.

QM/MM Free Energy Calculations of Long-Range Biological Protonation Dynamics by Adaptive and Focused Sampling.

Water-mediated proton transfer reactions are central for catalytic processes in a wide range of biochemical systems, ranging from biological energy conversion to chemical transformations in the metabolism. Yet, the accurate computational treatment of such complex biochemical reactions is highly challenging and requires the application of multiscale methods, in particular hybrid quantum/classical (QM/MM) approaches combined with free energy simulations. Here, we combine the unique exploration power of new advanced sampling methods with density functional theory (DFT)-based QM/MM free energy methods for multiscale simulations of long-range protonation dynamics in biological systems. In this regard, we show that combining multiple walkers/well-tempered metadynamics with an extended system adaptive biasing force method (MWE) provides a powerful approach for exploration of water-mediated proton transfer reactions in complex biochemical systems. We compare and combine the MWE method also with QM/MM umbrella sampling and explore the sampling of the free energy landscape with both geometric (linear combination of proton transfer distances) and physical (center of excess charge) reaction coordinates and show how these affect the convergence of the potential of mean force (PMF) and the activation free energy. We find that the QM/MM-MWE method can efficiently explore both direct and water-mediated proton transfer pathways together with forward and reverse hole transfer mechanisms in the highly complex proton channel of respiratory Complex I, while the QM/MM-US approach shows a systematic convergence of selected long-range proton transfer pathways. In this regard, we show that the PMF along multiple proton transfer pathways is recovered by combining the strengths of both approaches in a QM/MM-MWE/focused US (FUS) scheme and reveals new mechanistic insight into the proton transfer principles of Complex I. Our findings provide a promising basis for the quantitative multiscale simulations of long-range proton transfer reactions in biological systems.

Impact of an Ionic Liquid Solution on Horseradish Peroxidase Activity.

Horseradish peroxidase (HRP) is an enzyme that oxidizes pollutants from wastewater. A previous report indicated that peroxidases can have an enhancement in initial enzymatic activity in an aqueous solution of 0.26 M 1-ethyl-3-methylimidazolium ethyl sulfate ([EMIm][EtSO4]) at neutral pH. However, the atomistic details remain elusive. In the enzymatic landscape of HRP, compound II (Cpd II) plays a key role and involves a histidine (H42) residue. Cpd II exists as oxoferryl (2a) or hydroxoferryl (2b(FeIV)) forms, where 2a is the predominantly observed form in experimental studies. Intriguingly, the ferric 2b(FeIII) form seen in synthetic complexes has not been observed in HRP. Here, we have investigated the structure and dynamics of HRP in pure water and aqueous [EMIm][EtSO4] (0.26 M), as well as the reaction mechanism of 2a to 2b conversion using polarizable molecular dynamics (MD) simulations and quantum mechanics/molecular mechanics (QM/MM) calculations. When HRP is solvated in aq [EMIm][EtSO4], the catalytic water displaces, and H42 directly orients over the ferryl moiety, allowing a direct proton transfer (PT) with a significant energy barrier reduction. Conversely, in neat water, the reaction of 2a to 2b follows the previously reported mechanism. We further investigated the deprotonated form of H42. Analysis of the electric fields at the active site indicates that the aq [EMIm][EtSO4] medium facilitates the reaction by providing a more favorable environment compared with the system solvated in neat water. Overall, the atomic level supports the previous experimental observations and underscores the importance of favorable electric fields in the active site to promote catalysis.

Using the phospha-Michael reaction for making phosphonium phenolate zwitterions.

The reactions of 2,4-di-tert-butyl-6-(diphenylphosphino)phenol and various Michael acceptors (acrylonitrile, acrylamide, methyl vinyl ketone, several acrylates, methyl vinyl sulfone) yield the respective phosphonium phenolate zwitterions at room temperature. Nine different zwitterions were synthesized and fully characterized. Zwitterions with the poor Michael acceptors methyl methacrylate and methyl crotonate formed, but could not be isolated in pure form. The solid-state structures of two phosphonium phenolate molecules were determined by single-crystal X-ray crystallography. The bonding situation in the solid state together with NMR data suggests an important contribution of an ylidic resonance structure in these molecules. The phosphonium phenolates are characterized by UV-vis absorptions peaking around 360 nm and exhibit a negative solvatochromism. An analysis of the kinetics of the zwitterion formation was performed for three Michael acceptors (acrylonitrile, methyl acrylate, and acrylamide) in two different solvents (chloroform and methanol). The results revealed the proton transfer step necessary to stabilize the initially formed carbanion as the rate-determining step. A preorganization of the carbonyl bearing Michael acceptors allowed for reasonable fast direct proton transfer from the phenol in aprotic solvents. In contrast, acrylonitrile, not capable of forming a similar preorganization, is hardly reactive in chloroform solution, while in methanol the corresponding phosphonium phenolate is formed.

Computational study of the dimerization of glyphosate: mechanism and effect of solvent.

A computational study on the structure and stability of different series of glyphosate (Glyph) dimers comprising nonionized (N) and zwitterionic structures (Z) for neutral monomers, followed by an analysis of energetics of Glyph dimerization process have been performed by means of quantum chemical calculations in different media. Optimized geometries for energy minima, as well as relative potential and free energies of the possible various conformers of each series of Glyph dimer were computed as a function of the medium at B3LYP-D3/6-311++G(2d,2p) level. The solvation model based on density (SMD) is employed for all solution phase computations. Non-ionized dimers (DN), anion-cation (AC) and either zwitterion-zwitterion (DZP and DZC) or non-ionized-phosphonate zwitterion (NZP) ionized neutral forms of Glyph dimer are predicted to exist in the gas phase and in solution in large contrast to Glyph monomers. The DZC dimer form exhibiting a centrosymmetric arrangement of two carboxylate zwitterion units was found to be the most stable dimer structure in all media. In aqueous solution, the DZP and AC dimer type structures are significantly stabilized by hydration. The tautomerisms between DZC, DZP and AC dimer type structures have been investigated in the gas phase and in solution. The DZC type structures are more prone to experience proton transfer in water than in the gas phase and in cyclohexane. The mechanism for the tautomerization process in neutral ionized Glyph dimers proceeds via two direct proton transfer paths: DZP ⇋ AC ⇋ DZC. Results show that solvents play a key role in modulating the energetics of the dimerization process of Glyph. Solvation in cyclohexane, favors the dimerization process however, hydration opposed it. In aqueous solution, the mechanism of the dimerization of Glyph from its phosphonate zwitterionic monomer form (ZP) could be described by a set of equilibria including direct proton transfer paths as follows: 2ZP ⇋ DZP ⇋ AC ⇋ DZC. According to our results, in aqueous solution, DZC Glyph dimers and their corresponding DZP and AC tautomers should be present in higher concentration than phosphonate zwitterionic Glyph monomers for high Glyph concentration, a fact that seems controversial in the literature.

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.

Initiation Step in the Condensed Phase Decomposition Process of Ammonium Perchlorate

AbstractAmmonium perchlorate (AP) is a commonly used oxidizer in solid propellant compositions. As a result, numerous experimental and theoretical studies have been carried out in order to better understand the behavior of its decomposition in both the liquid and gas phases. In this work, the first step of Ammonium perchlorate (AP) decomposition in the condensed phase has been investigated using quantum mechanics‐based calculations. The process of proton transfer from ammonium ion (NH4+) to perchlorate ion (ClO4−) was investigated in‐depth as it is thought to be the preferred reaction for the initiation of AP decomposition in the condensed phase. The calculations were carried out using density functional theory (DFT) at the B3LYP level in conjunction with the 6‐311G++(d,p) basis set. The current study revealed that proton transfer from ammonium cation to perchlorate anion takes place in two steps. The first step involves the decomposition of ClO4− to ClO2− and O2 followed by the second step which is H+ transfer from NH4+ to ClO2− giving NH3 and HClO2. In order to better understand protonation and deprotonation, it is essential to know the proton affinity (PA) of the involved species. In this study, the CBS‐QB3 method is used to determine the PA of various energetic ionic compounds. These compounds include guanidinium azotetrazolate, triaminoguanidinium azotetrazolate, guanidinium nitrate, ammonium nitrate, and ammonium perchlorate. Based on the results of the calculations, it is concluded that direct proton transfer process does not occur in the condensed phase because the proton affinities of all cations are significantly higher than those of the corresponding anions. In conjunction with experimental studies, this modelling study can help in the development of a better understanding of AP decomposition in the condensed phase.

Comparison of H2O Adsorption and Dissociation Behaviors on Rutile (110) and Anatase (101) Surfaces Based on ReaxFF Molecular Dynamics Simulation.

The relationship between structure and reactivity plays a dominant role in water dissociation on the various TiO2 crystallines. To observe the adsorption and dissociation behavior of H2O, the reaction force field (ReaxFF) is used to investigate the dynamic behavior of H2O on rutile (110) and anatase (101) surfaces in an aqueous environment. Simulation results show that there is a direct proton transfer between the adsorbed H2O (H2Oad) and the bridging oxygen (Obr) on the rutile (110) surface. Compared with that on the rutile (110) surface, an indirect proton transfer occurs on the anatase (101) surface along the H-bond network from the second layer of water. This different mechanism of water dissociation is determined by the distance between the 5-fold coordinated Ti (Ti5c) and Obr of the rutile and anatase TiO2 surfaces, resulting in the direct or indirect proton transfer. Additionally, the hydrogen bond (H-bond) network plays a crucial role in the adsorption and dissociation of H2O on the TiO2 surface. To describe interfacial water structures between TiO2 and bulk water, the double-layer model is proposed. The first layer is the dissociated H2O on the rutile (110) and anatase (101) surfaces. The second layer forms an ordered water structure adsorbed to the surface Obr or terminal OH group through strong hydrogen bonding (H-bonding). Affected by the H-bond network, the H2O dissociation on the rutile (110) surface is inhibited but that on the anatase (101) surface is promoted.

The 46Ar(3He, d)47K direct reaction as a probe of the 46Ar proton wavefunction

The discrepancy between shell-model calculations and intermediate-energy Coulomb excitation measurements in 46Ar still stands as an unsolved puzzle in understanding the N = 28 shell evolution. This phenomenon has significant relevance considering the remarkable achievements of the shell model and the SDPF-U interaction in the region which is able to predict the fading of the N = 28 shell gap in neutron-rich 44S. Recent measurements narrowed down this discrepancy to an overestimation of the proton amplitude to the quadrupole transition matrix element. The current work aims to propose a different perspective on the puzzle, by studying a direct proton-transfer reaction on 46Ar as a means to directly probe the proton wavefunction of the ground state this isotope. By measuring the amount of l = 0 transfer to the ground state (1/2+) of 47K with respect to the l = 2 to the first excited state (3/2+), we aim to gain insight into the ground state proton wavefunction of 46Ar. We will present a brief description of the experiment performed at the SPIRAL1 facility in GANIL (France). The experimental apparatus allowed a full reconstruction of the two-body reaction thanks to the combination of AGATA, VAMOS, MUGAST, CATS2, and HECTOR.

Room Temperature Wearable Gas Sensors for Fabrication and Applications

AbstractThe recent surge in demand for human–machine interaction (HMI), Internet of Things (IoTs), and artificial intelligence (AI) has created both opportunities and challenges for room‐temperature wearable gas sensors. These sensors serve as a source of perceptual information and can be easily integrated into wearable electronic devices due to their portability and miniaturization. In recent years, various types of wearable room temperature gas sensors have been developed for fields like environmental monitoring, healthcare, smart home, industrial security, food safety monitoring, and public security. These sensors not only adjust to the movements of human effortlessly but also have reduced power consumption. Therefore, room temperature wearable gas sensors hold great promise for the development of integrated intelligent gas sensing system worn on the human body. These sensors can be fabricated using various sensing materials to detect diverse target gases. This review provides a comprehensive summary of the preparation of sensing materials with extraordinary sensing capabilities at room temperature. Additionally, the article includes a brief discussion of the sensing mechanism, employing four models to explain it: oxygen adsorption, direct electron transfer, proton transfer, and ions conduction. Finally, this article discusses the various applications and future perspectives of room‐temperature wearable gas sensors.

Mechanistic roles of metal- and ligand-protonated species in hydrogen evolution with [Cp*Rh] complexes

Protonation reactions involving organometallic complexes are ubiquitous in redox chemistry and often result in the generation of reactive metal hydrides. However, some organometallic species supported by η5-pentamethylcyclopentadienyl (Cp*) ligands have recently been shown to undergo ligand-centered protonation by direct proton transfer from acids or tautomerization of metal hydrides, resulting in the generation of complexes bearing the uncommon η4-pentamethylcyclopentadiene (Cp*H) ligand. Here, time-resolved pulse radiolysis (PR) and stopped-flow spectroscopic studies have been applied to examine the kinetics and atomistic details involved in the elementary electron- and proton-transfer steps leading to complexes ligated by Cp*H, using Cp*Rh(bpy) as a molecular model (where bpy is 2,2'-bipyridyl). Stopped-flow measurements coupled with infrared and UV-visible detection reveal that the sole product of initial protonation of Cp*Rh(bpy) is [Cp*Rh(H)(bpy)]+, an elusive hydride complex that has been spectroscopically and kinetically characterized here. Tautomerization of the hydride leads to the clean formation of [(Cp*H)Rh(bpy)]+. Variable-temperature and isotopic labeling experiments further confirm this assignment, providing experimental activation parameters and mechanistic insight into metal-mediated hydride-to-proton tautomerism. Spectroscopic monitoring of the second proton transfer event reveals that both the hydride and related Cp*H complex can be involved in further reactivity, showing that [(Cp*H)Rh] is not necessarily an off-cycle intermediate, but, instead, depending on the strength of the acid used to drive catalysis, an active participant in hydrogen evolution. Identification of the mechanistic roles of the protonated intermediates in the catalysis studied here could inform design of optimized catalytic systems supported by noninnocent cyclopentadienyl-type ligands.

The Catalytic Mechanism of Acetoacetate Decarboxylase: A Detailed Study of Schiff Base Formation, Protonation States, and Their Impact on Catalysis.

The enzyme acetoacetate decarboxylase (AAD) has a crucial function in the process of decarboxylating the substrate acetoacetate (AA). It has been extensively studied over the years, but its exact catalytic mechanism has remained partly unsolved due to the difficulty in assessing reaction intermediates. In this study, we combine molecular dynamics and electronic structure calculations to rediscover its catalytic mechanism. Our results show that the presence of the substrate, the acetoacetate, significantly influences the electrostatic potential of the active site. Furthermore, our simulations show that the decarboxylation reaction can take place by means of a direct proton transfer instead of via an enamine intermediate, which is thought to be strictly necessary. This work provides new insights into the role of the electrostatic interactions on the catalytic activity of AAD and for the first time connects it to the catalytic mechanism of other decarboxylases.

Site-Averaged Ab Initio Kinetics: Importance Learning for Multistep Reactions on Amorphous Supports

Single-atom centers on amorphous supports include catalysts for polymerization, partial oxidation, metathesis, hydrogenolysis, and more. The disordered environment makes each site different, and the kinetics exponentially magnifies these differences to make ab initio site-averaged kinetics calculations extremely difficult. This work extends the importance learning algorithm for efficient and precise site-averaged kinetics estimates to ab initio calculations and multistep reaction mechanisms. Specifically, we calculate site-averaged proton transfer relaxation rates on an ensemble of cluster models representing Brønsted acid sites on silica-alumina. We include direct and water-assisted proton transfer pathways and simultaneously estimate the water adsorption and activation enthalpies for forward and backward proton transfers. We use density functional theory (DFT) to obtain a site-averaged rate, somewhat like a turnover frequency, for the proton transfer relaxation rate. Finally, we show that importance learning can provide orders-of-magnitude acceleration over standard sampling methods for site-averaged rate calculations in cases where the rate is dominated by a few highly active sites.

Cyanide Binding to [FeFe]-Hydrogenase Stabilizes the Alternative Configuration of the Proton Transfer Pathway.

Hydrogenases are H2 converting enzymes that harbor catalytic cofactors in which iron (Fe) ions are coordinated by biologically unusual carbon monoxide (CO) and cyanide (CN- ) ligands. Extrinsic CO and CN- , however, inhibit hydrogenases. The mechanism by which CN- binds to [FeFe]-hydrogenases is not known. Here, we obtained crystal structures of the CN- -treated [FeFe]-hydrogenase CpI from Clostridium pasteurianum. The high resolution of 1.39 Å allowed us to distinguish intrinsic CN- and CO ligands and to show that extrinsic CN- binds to the open coordination site of the cofactor where CO is known to bind. In contrast to other inhibitors, CN- treated crystals show conformational changes of conserved residues within the proton transfer pathway which could allow a direct proton transfer between E279 and S319. This configuration has been proposed to be vital for efficient proton transfer, but has never been observed structurally.

Hemiacetal-based dynamic systems: a new mechanistic insight.

The formation of hemiacetals from pyrazine trifluoromethylketone as a model receptor and four simple alcohols was studied by using quantum chemical calculations and NMR spectroscopy. Free energy profiles for four types of mechanistic pathways were calculated and discussed with respect to kinetic and thermodynamic measurements. We show that hemiacetal formation is facilitated by an assisted proton transfer process via a pseudo eight-membered transition state which brings the theory and experiment into close agreement. Also, a newly proposed mechanistic pathway for hemiacetal formation via a five-membered transition state leading to zwitterionic intermediates is discussed. Direct proton transfer in a pseudo four-membered transition state can be ruled out due to the high energy of transition states with respect to other mechanistic pathways. We also show that in the case of hemiacetals, water and alcohol molecules cannot account sufficiently for the H-transfer process via six-membered transition states.

Contribution of the Collective Excitations to the Coupled Proton and Energy Transport along Mitochondrial Cristae Membrane in Oxidative Phosphorylation System.

The results of many experimental and theoretical works indicate that after transport of protons across the mitochondrial inner membrane (MIM) in the oxidative phosphorylation (OXPHOS) system, they are retained on the membrane-water interface in nonequilibrium state with free energy excess due to low proton surface-to-bulk release. This well-established phenomenon suggests that proton trapping on the membrane interface ensures vectorial lateral transport of protons from proton pumps to ATP synthases (proton acceptors). Despite the key role of the proton transport in bioenergetics, the molecular mechanism of proton transfer in the OXPHOS system is not yet completely established. Here, we developed a dynamics model of long-range transport of energized protons along the MIM accompanied by collective excitation of localized waves propagating on the membrane surface. Our model is based on the new data on the macromolecular organization of the OXPHOS system showing the well-ordered structure of respirasomes and ATP synthases on the cristae membrane folds. We developed a two-component dynamics model of the proton transport considering two coupled subsystems: the ordered hydrogen bond (HB) chain of water molecules and lipid headgroups of MIM. We analytically obtained a two-component soliton solution in this model, which describes the motion of the proton kink, corresponding to successive proton hops in the HB chain, and coherent motion of a compression soliton in the chain of lipid headgroups. The local deformation in a soliton range facilitates proton jumps due to water molecules approaching each other in the HB chain. We suggested that the proton-conducting structures formed along the cristae membrane surface promote direct lateral proton transfer in the OXPHOS system. Collective excitations at the water-membrane interface in a form of two-component soliton ensure the coupled non-dissipative transport of charge carriers and elastic energy of MIM deformation to ATP synthases that may be utilized in ATP synthesis providing maximal efficiency in mitochondrial bioenergetics.

Theoretical Insights into Enantioselective [3 + 2] Cycloaddition between Cinnamaldehyde and Cyclic N-Sulfonyl Trifluoromethylated Ketimine Catalyzed by N-Heterocyclic Carbene.

The density functional theory (DFT) method was used to investigate the mechanism and the origin of stereoselectivity of N-heterocyclic carbene (NHC)-catalyzed [3 + 2] cycloaddition between enals and cyclic imine N-sulfonyl trifluoromethyl ketimines at the M06-2X/SMD/6-311+G(d,p)//M06-2X/SMD/6-31G (d,p) level. The results show that the favorable reaction path consists of five steps: nucleophilic attack, proton transfer, the formation of the C-C bond, the tautomerism of the enol intermediate, the formation of the five-membered ring, and the regeneration of the catalyst. For the process of proton transfer, the base-assisted reaction can reduce the activation free energy and make the reaction easier to occur compared with the direct proton transfer process. The formation of the C-C bond is the crucial step of stereoselectivity, in which two chiral centers and four configurations of intermediates (RR/RS/SR/SS) were generated. The free energy barriers obtained and the noncovalent interaction analysis confirm that the dominant configuration is SS, becoming the final trans-type product observed in experiment. Furthermore, through the analyses of the conceptual DFT and natural atomic charges, it is revealed that NHC acts as a double catalyst, which can not only increase the nucleophilicity of reactants by Lewis base but also activate the C-H bond and promote the proton transfer process. The understanding of the mechanism obtained in this study should be helpful to the other organic catalytic reactions with high stereoselectivity.

A density functional theory study on the mechanism of simultaneous trifluoromethylation and oximation of aryl-substituted ethylenes

The effects of different substituents, located at the para position of the aromatic ring and at the β-carbon atom of styrenes, on difunctionalizations involving trifluoromethylation and oxime formation are investigated, showing that the difunctionalization reaction has a good adaptability to such reactants containing a range of substituents. This is important in the actual production process. It was found that proton transfer in the final tautomerism step involving transformation of a nitroso intermediate into an oxime is the rate-limiting step. The solvent effect did not influence the rate-limiting step significantly. Compared with direct proton transfer in a vacuum, the energy barrier of the final tautomerism step decreased from 57.80 kcal mol−1 in vacuum to 12.98 kcal mol−1 in water occurring via mediated proton transfer, which declines by 77.5%. When water participates in the rate-limiting steps in organic solvents, the energy barrier also decreases significantly, which indicates that a small amount of water in the organic solvent is conducive to the reaction.

Solvated Nuclear-Electronic Orbital Structure and Dynamics.

Nonadiabatic dynamical processes such as proton-coupled electron transfer and excited state intramolecular proton transfer have been the subject of much research. One of the promising theoretical methods to describe these processes is the nuclear-electronic orbital (NEO) approach. This approach inherently accounts for nuclear quantum effects within quantum chemistry calculations, and it has recently been extended to directly simulate nonadiabatic processes with the development of real-time NEO methods. These processes can also be significantly dependent on the surrounding chemical environment, however, and capturing the effects of the environment is often necessary for analyzing experimentally relevant systems. This work couples the NEO density functional theory and real-time time-dependent density functional theory approaches with solvation through the polarizable continuum model. The effects of this coupling are investigated for ground state properties, solvent-dependent vibrational frequencies, and direct excited state intramolecular proton transfer dynamics.