Active role of hydrogen bond and ambient water in Meyer–Schuster rearrangements in high‐temperature water

Meyer–Schuster rearrangements of 2‐phenyl‐3‐butyn‐2‐ol with H₃O⁺ and (H₂O)₆ model in high‐temperature water (HTW) have been investigated by the use of density functional theory calculations. In the substrate 2‐phenyl‐3‐butyn‐2‐ol catalyzed by H₃O⁺ and (H₂O)₆, the Meyer–Schuster rearrangements were predicted by the frontier molecular orbital theory. The results show that the rearrangement does not involve the carbonium ion intermediates, but the first transition state is carboniumion like. Dehydration and hydration may occur via the intermolecular proton relay along the hydrogen‐bond chains and the second step of reaction path is a total acid–base catalytic process. Based on the results, a model considered both HTW ambient and water molecules are proposed to represent mechanisms of other reactions in HTW. © 2011 Wiley Periodicals, Inc. Int J Quantum Chem, 2011     

Insights into the dynamics of evaporation and proton migration in protonated water clusters from Large‐scale Born–Oppenheimer direct dynamics

Large‐scale on‐the‐fly Born–Oppenheimer molecular dynamics simulations using recent advances in linear scaling electronic structure theory and trajectory integration techniques have been performed for protonated water clusters around the magic number (H₂O)_nH⁺, for n = 20 and 21. Besides demonstrating the feasibility and efficiency of the computational approach, the calculations reveal interesting dynamical details. Elimination of water molecules is found to be fast for both cluster sizes but rather insensitive to the initial geometry. The water molecules released acquire velocities compatible with thermal energies. The proton solvation shell changes between the well‐known Eigen and Zundel motifs and is characterized by specific low‐frequency vibrational modes, which have been quantified. The proton transfer mechanism largely resembles that of bulk water but one interesting variation was observed. © 2012 Wiley Periodicals, Inc.     

Proton dissociation and transfer in hydrated phosphoric acid clusters

The dynamics and mechanisms of proton dissociation and transfer in hydrated phosphoric acid (H₃PO₄) clusters under excess proton conditions were studied based on the concept of presolvation using the H₃PO₄–H₃O⁺–nH₂O complexes (n = 1–3) as the model systems and ab initio calculations and Born–Oppenheimer molecular dynamics (BOMD) simulations at the RIMP2/TZVP level as model calculations. The static results showed that the smallest, most stable intermediate complex for proton dissociation (n = 1) is formed in a low local‐dielectric constant environment (e.g., ε = 1), whereas proton transfer from the first to the second hydration shell is driven by fluctuations in the number of water molecules in a high local‐dielectric constant environment (e.g., ε = 78) through the Zundel complex in a linear H‐bond chain (n = 3). The two‐dimensional potential energy surfaces (2D‐PES) of the intermediate complex (n = 1) suggested three characteristic vibrational and ¹H NMR frequencies associated with a proton moving on the oscillatory shuttling and structural diffusion paths, which can be used to monitor the dynamics of proton dissociation in the H‐bond clusters. The BOMD simulations over the temperature range of 298–430 K validated the proposed proton dissociation and transfer mechanisms by showing that good agreement between the theoretical and experimental data can be achieved with the proposed rate‐determining processes. The theoretical results suggest the roles played by the polar solvent and iterate that insights into the dynamics and mechanisms of proton transfer in the protonated H‐bond clusters can be obtained from intermediate complexes provided that an appropriate presolvation model is selected and that all of the important rate‐determining processes are included in the model calculations. © 2015 Wiley Periodicals, Inc.     

A Gas‐Phase CanMn4−nO4+ Cluster Model for the Oxygen‐Evolving Complex of Photosystem II

One of the fundamental processes in nature, the oxidation of water, is catalyzed by a small CaMn₃O₄?MnO cluster located in photosystem II (PS II). Now, the first successful preparation of a series of isolated ligand‐free tetrameric Ca_nMn_4?nO₄⁺ (n=0–4) cluster ions is reported, which are employed as structural models for the catalytically active site of PS II. Gas‐phase reactivity experiments with D₂O and H₂¹⁸O in an ion trap reveal the facile deprotonation of multiple water molecules via hydroxylation of the cluster oxo bridges for all investigated clusters. However, only the mono‐calcium cluster CaMn₃O₄⁺ is observed to oxidize water via elimination of hydrogen peroxide. First‐principles density functional theory (DFT) calculations elucidate mechanistic details of the deprotonation and oxidation reactions mediated by CaMn₃O₄⁺ as well as the role of calcium.     

The Hydrated Excess Proton in the Zundel Cation H5O2+: The Role of Ultrafast Solvent Fluctuations

The nature of the excess proton in liquid water has remained elusive after decades of extensive research. In view of ultrafast structural fluctuations of bulk water scrambling the structural motifs of excess protons in water, we selectively probe prototypical protonated water solvates in acetonitrile on the femtosecond time scale. Focusing on the Zundel cation H₅O₂⁺ prepared in room‐temperature acetonitrile, we unravel the distinct character of its vibrational absorption continuum and separate it from OH stretching and bending excitations in transient pump‐probe spectra. The infrared absorption continuum originates from a strong ultrafast frequency modulation of the H⁺ transfer vibration and its combination and overtones. Vibrational lifetimes of H₅O₂⁺ are found to be in the sub‐100 fs range, much shorter than those of unprotonated water. Theoretical results support a picture of proton hydration where fluctuating electrical interactions with the solvent and stochastic thermal excitations of low‐frequency modes continuously modify the proton binding site while affecting its motions.     

Proton transfer at the carboxylic sites of amino acids: A single water molecule catalyzed process

Ab initio calculations at MP2 level of theory were used to study the proton transfer at the carboxylic sites of amino acids, in the isolated, mono‐ and di‐hydrated forms. In the case of water dimer, two interaction modes with glycine neutral structures (see Fig. 3 ) were explored, corresponding to the concerted and stepwise reaction pathways. Their transition states can be described as (H₂O? H? OH₂)⁺ [Fig. 4 (a)] and (H₂O‐‐‐H? OH₂)⁺ [Fig. 4 (b)], respectively. The energy analysis indicated that the concerted pathway is preferred. In the isolated, mono‐ and di‐hydrated glycine complexes, the activation barriers of the proton transfer at the carboxylic sites were calculated to be 34.49, 16.59, and 13.36 kcal mol^?1, respectively. It was thus shown that the proton transfer is significantly assisted and catalyzed by water monomer so that it can take place at room temperature. Instead, the further addition of water molecules plays solvent effects rather than catalytic effects to this proton transfer process. The above results obtained with discrete water molecules were supported by the solvent continuum calculated data. It was also observed that the heavy dependence of the solvent continuum models on dipole moments may produce misleading results. © 2008 Wiley Periodicals, Inc. Int J Quantum Chem, 2009     

A DFT study of proton transfers for the reaction of phenol and hydroxyl radical leading to dihydroxybenzene and H2O in the water cluster

Reactions of phenol and hydroxyl radical were studied under the aqueous environment to investigate the antioxidant characters of phenolic compounds. M06‐2X/6‐311 + G(d,p) calculations were carried out, where proton transfers via water molecules were examined carefully. Stepwise paths from phenol + OH^• + (H₂O)_n (n = 3, 7, and 12) to the phenoxyl radical (Ph O^•) via dihydroxycyclohexadienyl radicals (ipso, ortho, meta, and para OH‐adducts) were obtained. In those paths, the water dimer was computed to participate in the bond interchange along hydrogen bonds. The concerted path corresponding to the hydrogen atom transfer (HAT, apparently Ph OH + OH^• → Ph O^• + H₂O) was found. In the path, the hydroxyl radical located on the ipso carbon undergoes the charge transfer to prompt the proton (not hydrogen) transfer. While the present new mechanism is similar to the sequential proton loss electron transfer (SPLET) one, the former is of the concerted character. Tautomerization reactions of ortho or para (OH)C₆H₅=O + (H₂O)_n → C₆H₄(OH)₂(H₂O)_n were traced with n = 2, 3, 4, and 14. The n = 3 (and n = 14) model of ortho and para was calculated to be most likely by the strain‐less hydrogen‐bond circuit.     

Thermodynamics of hydrogen-isotope-exchange reactions. Proton solvation inn-hexanol and water

This paper reports our results for the direct experimental determination of the equilibrium constant for the hydrogen-isotope-exchange reaction, 1/2D₂(g)+HCl(hexOH)=1/2H₂(g)+DCl(hexOD), where hexOH isn-hexanol and hexOD isn-hexanol with deuterium substitution in the alcohol function. The reaction was studied in electrochemical double cells without liquid junction for which the net cell reaction is the above isotope-exchange reaction. The experimentally determined value of ε° (296.0°K) for this cell is 4.03±0.95 mV (strong electrolyte standard states, mole-fraction composition scale); the value of the equilibrium constant for the reaction is 1.17±0.05. The contributions of isotope-exchange and transfer effects to the magnitude of the standard Gibbs energy change for the above reaction and for the analogous reaction 1/2D₂(g)+HCl(aq)=DCl(daq)+1/2H₂(g) are considered. Our results support the conclusion of Heinzinger and Weston that the formulation of the solvated proton in water as H₃O⁺, as opposed to H₉O₄ ⁺, is sufficient for the interpretation of the thermodynamics of hydrogen-isotope-exchange reactions in water. We also find that the formulation of the solvated proton inn-hexanol as ROH ₂ ⁺ is sufficient for the interpretation of our results on the thermodynamics of hydrogen-isotope-exchange inn-hexanol.     

AM1 studies on the potential energy surface for the proton transfer in protonated water clusters,H+(H2O)n

The potential energy surfaces for the proton transfer processes in H⁺(H₂O)_n with n=2 ～ 11 have been studied using the semiempirical AM1 method. Two model systems were adopted: branched and linear systems. The branched system showed a tendency to form a bulk cluster, while the linear system showed a tendency toward a constant barrier height with increasing number of water molecules in the model system. The potential energy surfaces were discussed using Marcus theory. In the case of H⁺ (H₂O)_n with n=10 and 11, the intrinsic barrier to the proton transfer was found to be around 1.0 kcal/mol.     

An application of the reaction field theory to hydrated metal cations in the framework of the MNDO,AM1, and PM3 methods

A reaction field theory, combined with the MNDO, AM1, and PM3 molecular orbital methods, was applied to hydration phenomena of metal cationic species. The first hydration shell was treated explicitly by using a supermolecular model, [M(H₂O)_n]^m+, and its surrounding medium was described with a continuum dielectric. Hydration free energies were evaluated as a sum of the contributions from the electrostatic interaction with the bulk medium, the hydrated cluster formation, the cavity formation, and the vaporization of water molecules forming the cluster. As a whole, calculated hydration energies were in good agreement with the corresponding experimental data over various kinds of metal cationic species. © 1995 by John Wiley & Sons, Inc.     

A Hybrid statistical mechanics—quantum chemical model for proton transfer in 5‐azauracil and 6‐azauracil in water solution

A hybrid statistical physics—quantum‐chemical methodology was implemented to study the water‐assisted intramolecular proton‐transfer processes in 5‐ and 6‐azauracils in aqueous solutions. The solvent effects were included in the model by explicit inclusion of two pairs of water molecules, which model the relevant part of the first hydration shell around the solute. The position of these water molecules was initially estimated by carrying out a classical Metropolis of dilute water solutions of the title compounds and subsequently analyzing solute–solvent intermolecular interactions in the Monte Carlo‐generated configurations. Sequentially to the statistical physics simulation, ab initio quantum mechanical (QM) level of theory was implemented. The effects of the water as solvent (at ab initio QM level) were introduced at two different levels—using solute–solvent clusters (four‐water molecules) and using the same clusters embedded in an external continuum. Full geometry optimizations of these complexes were carried out at MP2/6–31 + G(d, p) and conductor‐polarizable continuum model (C‐PCM)/MP2/6–31 + G(d, p). Single point calculations were performed at CCSD(T)/6–31 + G(d, p)//MP2/6–31 + G(d, p) computational level to obtain more accurate energies. According to our calculations hydrated azauracils should exist in three forms: mainly dioxo form and two hydroxy forms. The calculated proton transfer activation energies for tautomeric reactions of 5‐azauracil and 6‐azauracil show different pictures for these two compounds. According to C‐PCM/MP2/6–31 +G(d, p) data, water‐assisted proton transfer in 5‐azauracil realizes through two parallel reactions: 1,3,5‐triazine‐2,4(1H,3H)‐dione → 6‐hydroxy‐1,3,5‐triazin‐2(1H)‐one and 1,3,5‐triazine‐2,4(1H,3H)‐dione → 4‐hydroxy‐1,3,5‐triazin‐2(1H)‐one. Tautomeric equilibrium in 6‐azauracil in water could occur by two contiguous reactions: 1,2,4‐triazine‐3,5(2H,4H)‐dione → 5‐hydroxy‐1,2,4‐triazin‐3(2H)‐one and 5‐hydroxy‐1,2,4‐triazin‐3(2H)‐one → 3‐hydroxy‐1,2,4‐triazin‐5(2H)‐one. The proton transfer investigated reactions in 5‐ and 6‐azauracils involve concerted atomic movement. © 2015 Wiley Periodicals, Inc.     

Double‐Oxygen‐Atom Transfer in Reactions of CemO2m+ (m=2–6) with C2H2

Cerium oxide cluster cations (Ce_mO_n⁺, m=2–16; n=2m, 2m±1 and 2m±2) are prepared by laser ablation and reacted with acetylene (C₂H₂) in a fast‐flow reactor. A time‐of‐flight mass spectrometer is used to detect the cluster distribution before and after the reactions. Reactions of stoichiometric Ce_mO_2m⁺ (m=2–6) with C₂H₂ produce Ce_mO_2m?2⁺ clusters, which indicates a “double‐oxygen‐atom transfer” reaction Ce_mO_2m⁺+C₂H₂→Ce_mO_2m?2⁺+(CHO)₂ (ethanedial). A single‐oxygen‐atom transfer reaction channel is also identified as Ce_mO_2m⁺+C₂H₂→Ce_mO_2m?1⁺+C₂H₂O (at least for m=2 and 3). Density functional theory calculations are performed to study reaction mechanisms of Ce₂O₄⁺+C₂H₂, and the calculated results confirm that both the single‐ and double‐oxygen‐atom transfer channels are thermodynamically and kinetically favourable.     

Effects of zero point vibration on the reaction dynamics of water dimer cations following ionization

Reactions of water dimer cation following ionization have been investigated by means of a direct ab initio molecular dynamics method. In particular, the effects of zero point vibration and zero point energy (ZPE) on the reaction mechanism were considered in this work. Trajectories were run on two electronic potential energy surfaces (PESs) of : ground state (²A″‐like state) and the first excited state (²A′ ‐ like state). All trajectories on the ground‐state PES lead to the proton‐transferred product: H₂O⁺(Wd)‐H₂O(Wa) → OH(Wd)‐H₃O⁺(Wa), where Wd and Wa refer to the proton donor and acceptor water molecules, respectively. Time of proton transfer (PT) varied widely from 15 to 40 fs (average time of PT = 30.9 fs). The trajectories on the excited‐state PES gave two products: an intermediate complex with a face‐to‐face structure (H₂O‐OH₂)⁺ and a PT product. However, the proton was transferred to the opposite direction, and the reverse PT was found on the excited‐state PES: H₂O(Wd)‐H₂O⁺ (Wa) → H₃O⁺(Wd)‐OH(Wa). This difference occurred because the ionizing water molecule in the dimer switched between the ground and excited states. The reaction mechanism of and the effects of ZPE are discussed on the basis of the results. © 2017 Wiley Periodicals, Inc.     

The Energy Barrier to Recombination in Hydrated Plasma

An abnormally high potential barrier that separates the H₃O⁺ and Cl^? ions in a cluster of water molecules was revealed. The profile of the barrier was calculated by computer simulation. The calculation was based on a detailed model of intermolecular interactions developed on the basis of experimental data on the free energy and entropy for the addition reactions of water molecules in the vapor phase to the hydration shell of ions in conjunction with the results of quantum-chemical calculations.     

Structure and energy of the positively ionized water clusters

Geometry optimization of small (H₂O)_n⁺ clusters (n ≤ 4) at the UHF/4–31 + + G^** level indicates that the cations consist of two fragments: the OH radical and the H_2n−1 O⁺_n−1 ion. The latter can be considered as a thermodynamically stable combination of a distorted H₃O⁺ ion and (n−2) H₂O molecules. The H bond between the fragments becomes weaker with increasing cluster size. Extrapolation of the adiabatic ionization potentials calculated for the (H₂O)_n oligomers (n ≤ 4) at the MP2 level to an infinite cluster size provides the value of approximately 8.7 eV, which can be presumably necessary for the ionization of liquid water in a vacuum. © 1997 John Wiley & Sons, Inc.     

On the Gas‐Phase Reactivity of Complexed OH+ with Halogenated Alkanes

OH⁺ is an extraordinarily strong oxidant. Complexed forms (L? OH⁺), such as H₂OOH⁺, H₃NOH⁺, or iron–porphyrin‐OH⁺ are the anticipated oxidants in many chemical reactions. While these molecules are typically not stable in solution, their isolation can be achieved in the gas phase. We report a systematic survey of the influence on L on the reactivity of L? OH⁺ towards alkanes and halogenated alkanes, showing the tremendous influence of L on the reactivity of L? OH⁺. With the help of with quantum chemical calculations, detailed mechanistic insights on these very general reactions are gained. The gas‐phase pseudo‐first‐order reaction rates of H₂OOH⁺, H₃NOH⁺, and protonated 4‐picoline‐N‐oxide towards isobutane and different halogenated alkanes C_nH_2n+1Cl (n=1–4), HCF₃, CF₄, and CF₂Cl₂ have been determined by means of Fourier transform ion cyclotron resonance meaurements. Reaction rates for H₂OOH⁺ are generally fast (7.2×10^?10–3.0×10^?9 cm³ mol^?1 s^?1) and only in the cases HCF₃ and CF₄ no reactivity is observed. In contrast to this H₃NOH⁺ only reacts with tC₄H₉Cl (k_obs=9.2×10^?10), while 4‐CH₃‐C₅H₄N‐OH⁺ is completely unreactive. While H₂OOH⁺ oxidizes alkanes by an initial hydride abstraction upon formation of a carbocation, it reacts with halogenated alkanes at the chlorine atom. Two mechanistic scenarios, namely oxidation at the halogen atom or proton transfer are found. Accurate proton affinities for HOOH, NH₂OH, a series of alkanes C_nH_2n+2 (n=1–4), and halogenated alkanes C_nH_2n+1Cl (n=1–4), HCF₃, CF₄, and CF₂Cl₂, were calculated by using the G3 method and are in excellent agreement with experimental values, where available. The G3 enthalpies of reaction are also consistent with the observed products. The tendency for oxidation of alkanes by hydride abstraction is expressed in terms of G3 hydride affinities of the corresponding cationic products C_nH_2n+1⁺ (n=1–4) and C_nH_2nCl⁺ (n=1–4). The hypersurface for the reaction of H₂OOH⁺ with CH₃Cl and C₂H₅Cl was calculated at the B3 LYP, MP2, and G3_m* level, underlining the three mechanistic scenarios in which the reaction is either induced by oxidation at the hydrogen or the halogen atom, or by proton transfer.     

Tautomerism and proton transfer in photoionized acetaldehyde and acetaldehyde–water clusters

Understanding the gas‐phase chemistry of acetaldehyde can be challenging because the molecule can assume several tautomeric forms, namely keto, enol and carbene. The two last forms are the most stable ionic forms. Here, insight into the gas‐phase cluster ion chemistry of homogeneous acetaldehyde and mixed water–acetaldehyde clusters is provided by mass spectrometry/vacuum ultraviolet photoionization combined with density functional theory calculations. (AA)_nH⁺ clusters (AA = acetaldehyde) and mixed (AA)_nH₃O⁺ clusters were detected using tunable vacuum ultraviolet photoionization. Barrierless proton transfers were observed during the geometry optimization of the most stable dimer structures and helped to explain the cluster ion chemistry induced by photoionization, namely the formation of deprotonated tautomers and protonated keto tautomers. Water was found to catalyze the keto–enol and keto–carbene isomerizations and facilitate the proton transfer from the ionized enol or carbene part of the cluster to the neutral keto part, resulting in protonated keto structures. The production of protonated keto structures was identified to be the main fragmentation channel following ionization of the homogeneous acetaldehyde cluster and a channel for ionized mixed clusters as well. These findings are significant for a broad range of fields, including current atmospheric models, because acetaldehyde is one of the most prominent organic species in the troposphere and ions play a crucial role in aerosol formation. Copyright © 2014 John Wiley & Sons, Ltd.     

Theoretical study on the role of cooperative solvent molecules in the neutral hydrolysis of ketene

The results of a theoretical study of the reaction mechanism for the neutral hydration of ketene, H₂C=C=O + (n + 1) H₂O → CH₃COOH + nH₂O (n = 0–4), in solution are presented. All structures were optimized and characterized at the MP2(fc)/6-31 + G* level of theory, and then re-optimized by MP2(fc)/6-311 ++G**, and the effect of the bulk solvent is taken into account according to the conductor-like polarized continuum model (CPCM) using the gas MP2(fc)/6-311 ++G** geometries. Energies were refined for five-water hydration at higher level of theory, QCISD(T)(fc)/6-311 ++G**//MP2(fc)/6-311 ++G**. In the combined supermolecular/continuum model, one water molecule directly attacks the central C-atom, and the other four explicit water molecules are divided into two groups, one acting as catalyst(s) by participating in the proton transfer to reduce the tension of proton transfer ring, and the other being placed near the non-reactive oxygen or carbon atom in order to catalyze the hydration by engaging in hydrogen-bonding to the substrate (the so-called cooperative effect). Between the two possible nucleophilic addition reactions of water molecule, across the C=O bond or the C=C bond, the former one is preferred. Our calculations suggest that the favorable hydrolysis mechanism of ketene involves a sort of eight-membered ring transition structure formed by a three-water proton transfer loop, and a cooperative dimeric water near the non-reactive carbon-atom. The best-estimated in the present paper for the rate-determining barrier in solution, $ \Updelta G_{\text{sol}}^{ \ne } $ (298 K), is about 58 kJ/mol, reasonably close to the available experimental result.     

Oxygen reduction on a platinum cluster

The reaction Pt₅O₂(ads)+H⁺(aq)+e⁻→Pt₅–OOH is analyzed on the basis of density functional theory calculations. It is found that the electron transfer process takes place gradually as the hydronium ion gets close to the adsorbed oxygen. At a certain small O_ads–H⋯O_water distance, the barrier for proton and electron transfer becomes negligible. The effect of an electric field on the reaction is studied by charging the metal/adsorbate complex, and that of a solvation shell on the proton transfer process is explored by using the H₃O⁺(H₂O)₂ ion cluster to model the hydrated proton.     

Challenges in predicting ΔrxnG in solution: Hydronium,hydroxide, and water autoionization

Standard non‐semiempirical continuum‐dielectric orbital‐based methods horribly overpredict, by 26‐50 kcal mol⁻¹, the Gibbs energy for the water autoionization reaction 2 H₂O_(l) → H₃O⁺_(aq) + OH^–_(aq). Here, we demonstrate these errors, fully investigate the reasons for these errors, and show that the use of 4 explicit solvent within the continuum (the “semicontinuum,” “cluster‐continuum,” or “hybrid” technique) can reduce the error of a standard continuum model from 50 to 2 kcal mol⁻¹. Results from pure cluster, pure continuum (several versions including semiempirical ones), and semicontinuum modeling are each presented and discussed. We recommend use of 3 waters around hydronium and 4 waters around hydroxide with standard continua whenever these ions are involved in reaction. To the possible surprise of some, time‐consuming molecular‐dynamics simulations are not needed to reproduce this problematic energy.