Vibrational Analysis of the H5O2+ Infrared Spectrum Using Molecular and Driven Molecular Dynamics

Graph Theoretic Molecular Fragmentation for Multidimensional Potential Energy Surfaces Yield an Adaptive and General Transfer Machine Learning Protocol

Over a series of publications we have introduced a graph-theoretic description for molecular fragmentation. Here, a system is divided into a set of nodes, or vertices, that are then connected through edges, faces, and higher-order simplexes to represent a collection of spatially overlapping and locally interacting subsystems. Each such subsystem is treated at two levels of electronic structure theory, and the result is used to construct many-body expansions that are then embedded within an ONIOM-scheme. These expansions converge rapidly with many-body order (or graphical rank) of subsystems and have been previously used for ab initio molecular dynamics (AIMD) calculations and for computing multidimensional potential energy surfaces. Specifically, in all these cases we have shown that CCSD and MP2 level AIMD trajectories and potential surfaces may be obtained at density functional theory cost. The approach has been demonstrated for gas-phase studies, for condensed phase electronic structure, and also for basis set extrapolation-based AIMD. Recently, this approach has also been used to derive new quantum-computing algorithms that enormously reduce the quantum circuit depth in a circuit-based computation of correlated electronic structure. In this publication, we introduce (a) a family of neural networks that act in parallel to represent, efficiently, the post-Hartree-Fock electronic structure energy contributions for all simplexes (fragments), and (b) a new k-means-based tessellation strategy to glean training data for high-dimensional molecular spaces and minimize the extent of training needed to construct this family of neural networks. The approach is particularly useful when coupled cluster accuracy is desired and when fragment sizes grow in order to capture nonlocal interactions accurately. The unique multidimensional k-means tessellation/clustering algorithm used to determine our training data for all fragments is shown to be extremely efficient and reduces the needed training to only 10% of data for all fragments to obtain accurate neural networks for each fragment. These fully connected dense neural networks are then used to extrapolate the potential energy surface for all molecular fragments, and these are then combined as per our graph-theoretic procedure to transfer the learning process to a full system energy for the entire AIMD trajectory at less than one-tenth the cost as compared to a regular fragmentation-based AIMD calculation.

Permutationally invariant polynomial representation of polarizability tensor surfaces for linear regression analysis

A linearly parameterized functional form for a Cartesian representation of molecular dipole polarizability tensor surfaces (PTS) is described. The proposed expression for the PTS is a linearization of the recently reported power series ansatz of the original Applequist model, which by construction is non-linear in parameter space. This new approach possesses (i) a unique solution to the least-squares fitting problem; (ii) a low level of the computational complexity of the resulting linear regression procedure, comparable to those of the potential energy and dipole moment surfaces; and (iii) a competitive level of accuracy compared to the non-linear PTS model. Calculations of CH₄ PTS, with polarizabilities fitted to 9000 training set points with the energies up to 14,000 cm^-1 show an impressive level of accuracy of the linear PTS model obtained with ~1600 parameters: ~1% versus 0.3% RMSE for the non-linear vs. linear model on a test set of 1000 configurations.

Intramolecular Proton Transfer in the Hydrogen Oxalate Anion and the Cooperativity Effects of the Low-Frequency Vibrations: A Driven Molecular Dynamics Study

We report first-principles molecular dynamics (MD) and dipole-driven molecular dynamics (μ-DMD) simulations of the hydrogen oxalate anion at the MP2/aug-cc-pVDZ level of theory. We examine the role of vibrational coupling between the OH stretching bands, that is, the fundamental and a few combination bands spanning the 2900-3100 cm^-1 range, and several of the low-frequency bending and stretching fundamental modes. The low-frequency modes between 300 and 825 cm^-1 play a crucial role in the proton-transfer motion. Strong involvement of CO₂ and CCO bending and the CC stretching vibrations indicate that these large amplitude motions cause the shortening of the O···O distance and thus promote H⁺ transfer to the other oxygen by bringing it over the 3.4 kcal/mol barrier. Analysis of resonant μ-DMD trajectories shows that the complex spectral feature near 825 cm^-1, closely corresponding to both an overtone of two quanta of 425 cm^-1 and a combination band of low-frequency CO₂ rocking (300 cm^-1) and CCO bending (575 cm^-1) modes, is involved in the proton transfer. μ-DMD shows that exciting the system at these mode combinations leads to faster barrier activation than exciting at the OH fundamental mode.

Coupled Local-Mode Approach for the Calculation of Vibrational Spectra: Application to Protonated Water Clusters

Spectroscopic studies of protonated water clusters (PWCs) have yielded enormous insights into the fundamental nature of the hydrated proton. Here, we introduce a new coupled local-mode (CLM) approach to calculate PWC OH stretch vibrational spectra. The CLM method combines a sampling of representative configurations from density functional theory (DFT)-based ab initio molecular dynamics (AIMD) simulations with DFT calculations of local-mode vibrational frequencies and couplings. Calculations of inhomogeneous OH stretch vibrational spectra for H⁺(H₂O)₄ and H⁺(H₂O)₂₁ agree well with experiment and higher-level calculations, and decompositions of the calculated spectra in terms of the coupled modes aids in the interpretation of the spectra. This observation is consistent with the idea that capturing anharmonicity and coupling is as important to accuracy as the underlying level of electronic structure theory. The CLM calculations can easily discern the configuration that dominates the experimental measurement for H⁺(H₂O)₅, which can adopt several low-energy conformations.

Demystifying the Diffuse Vibrational Spectrum of Aqueous Protons Through Cold Cluster Spectroscopy

The ease with which the pH is routinely determined for aqueous solutions masks the fact that the cationic product of Arrhenius acid dissolution, the hydrated proton, or H⁺(aq), is a remarkably complex species. Here, we review how results obtained over the past 30 years in the study of H⁺⋅(H₂O)_n cluster ions isolated in the gas phase shed light on the chemical nature of H⁺(aq). This effort has also revealed molecular-level aspects of the Grotthuss relay mechanism for positive-charge translocation in water. Recently developed methods involving cryogenic cooling in radiofrequency ion traps and the application of two-color, infrared-infrared (IR-IR) double-resonance spectroscopy have established a clear picture of how local hydrogen-bond topology drives the diverse spectral signatures of the excess proton. This information now enables a new generation of cluster studies designed to unravel the microscopic mechanics underlying the ultrafast relaxation dynamics displayed by H⁺(aq).

Analysis of the Proton Transfer Bands in the Infrared Spectra of Linear N2H+···OC and N2D+···OC Complexes Using Electric Field-Driven Classical Trajectories

In this work, we describe ab initio calculations and assignment of infrared (IR) spectra of hydrogen-bonded ion-molecular complexes that involve a fluxional proton: the linear N₂H⁺···OC and N₂D⁺···OC complexes. Given the challenges of describing fluxional proton dynamics and especially its IR activity, we use electric field-driven classical trajectories, i.e., the driven molecular dynamics (DMD) method that was developed by us in recent years and for similar applications, in conjunction with high-level electronic structure theory. Namely, we present a modified and a numerically efficient implementation of DMD specifically for direct (or "on the fly") calculations, which we carry out at the MP2-F12/AVDZ level of theory for the potential energy surface (PES) and MP2/AVDZ for the dipole moment surfaces (DMSs). Detailed analysis of the PES, DMS, and the time-dependence of the first derivative of the DMS, referred to as the driving force, for the highly fluxional vibrations involving H⁺/D⁺ revealed that the strongly non-harmonic PES and non-linear DMS yield remarkably complex vibrational spectra. Interestingly, the classical trajectories reveal a doublet in the proton transfer part of the spectrum with the two peaks at 1800 and 1980 cm^-1. We find that their shared intensity is due to a Fermi-like resonance interaction, within the classical limit, of the H⁺ parallel stretch fundamental and an H⁺ perpendicular bending overtone. This doublet is also observed in the deuterated species at 1360 and 1460 cm^-1.

Reduced rovibrational coupling Cartesian dynamics for semiclassical calculations: Application to the spectrum of the Zundel cation

We study the vibrational spectrum of the protonated water dimer, by means of a divide-and-conquer semiclassical initial value representation of the quantum propagator, as a first step in the study of larger protonated water clusters. We use the potential energy surface from the work of Huang et al. [J. Chem. Phys. 122, 044308 (2005)]. To tackle such an anharmonic and floppy molecule, we employ fully Cartesian dynamics and carefully reduce the coupling to global rotations in the definition of normal modes. We apply the time-averaging filter and obtain clean power spectra relative to suitable reference states that highlight the spectral peaks corresponding to the fundamental excitations of the system. Our trajectory-based approach allows for the physical interpretation of the very challenging proton transfer modes. We find that it is important, for such a floppy molecule, to selectively avoid initially exciting lower energy modes, in order to obtain cleaner spectra. The estimated vibrational energies display a mean absolute error (MAE) of ∼29 cm^-1 with respect to available multiconfiguration time-dependent Hartree calculations and MAE ∼ 14 cm^-1 when compared to the optically active experimental excitations of the Ne-tagged Zundel cation. The reasonable scaling in the number of trajectories for Monte Carlo convergence is promising for applications to higher dimensional protonated cluster systems.

Assignment of Infrared-Active Combination Bands in the Vibrational Spectra of Protonated Molecular Clusters Using Driven Classical Trajectories: Application to N4H+ and N4D. J Phys Chem A 2019;123:5613-5620. [PMID: 31244119 DOI: 10.1021/acs.jpca.9b04480]

We investigate the utility of the driven molecular dynamics (DMD) approach to complex molecular vibrations by applying it to linear clusters with several degenerate vibrational modes and infrared (IR) intense combination bands. Here, the prominent features in N₄H⁺ and N₄D⁺ IR spectra, reported and described by others previously, have been characterized for the first time by DMD using recently published high-level potential and dipole moment surfaces. Namely, the calculations closely correlate the parallel proton stretch vibration in N₄H⁺, at 750 cm^-1, with the one observed experimentally at 743 cm^-1. Second, the intense IR-active combination bands found in experimental spectra within 900-1100 cm^-1 have been properly recovered by DMD at 950 cm^-1 as strongly IR-active and confirmed as consisting of H⁺ asymmetric stretch and N₂···N₂ intermolecular symmetric stretch modes. Furthermore, we show that certain combination bands involving overtone transitions may be recovered by DMD using a hard-driving regime, such as the 1409 cm^-1 band measured in N₄H⁺, revealed by DMD at 1375 cm^-1, and assigned to a progressive combination of the parallel H⁺ stretch and two quanta of N₂···N₂ stretch, in agreement with quantum mechanical studies reported previously by others.

Introductory lecture: advances in ion spectroscopy: from astrophysics to biology

This introduction provides a historical context for the development of ion spectroscopy over the past half century by following the evolution of experimental methods to the present state-of-the-art. Rather than attempt a comprehensive review, we focus on how early work on small ions, carried out with fluorescence, direct absorption, and photoelectron spectroscopy, evolved into powerful technologies that can now address complex chemical problems ranging from catalysis to biophysics. One of these developments is the incorporation of cooling and temperature control to enable the general application of "messenger tagging" vibrational spectroscopy, first carried out using ionized supersonic jets and then with buffer gas cooling in radiofrequency ion traps. Some key advances in the application of time-resolved pump-probe techniques to follow ultrafast dynamics are also discussed, as are significant benchmarks in the refinement of ion mobility to allow spectroscopic investigation of large biopolymers with well-defined shapes. We close with a few remarks on challenges and opportunities to explore molecular level mechanics that drive macroscopic behavior.

Decoding the spectroscopic features and time scales of aqueous proton defects

Acid solutions exhibit a variety of complex structural and dynamical features arising from the presence of multiple interacting reactive proton defects and counterions. However, disentangling the transient structural motifs of proton defects in the water hydrogen bond network and the mechanisms for their interconversion remains a formidable challenge. Here, we use simulations treating the quantum nature of both the electrons and nuclei to show how the experimentally observed spectroscopic features and relaxation time scales can be elucidated using a physically transparent coordinate that encodes the overall asymmetry of the solvation environment of the proton defect. We demonstrate that this coordinate can be used both to discriminate the extremities of the features observed in the linear vibrational spectrum and to explain the molecular motions that give rise to the interconversion time scales observed in recent nonlinear experiments. This analysis provides a unified condensed-phase picture of the proton structure and dynamics that, at its extrema, encompasses proton sharing and spectroscopic features resembling the limiting Eigen [H₃O(H₂O)₃]⁺ and Zundel [H(H₂O)₂]⁺ gas-phase structures, while also describing the rich variety of interconverting environments in the liquid phase.

IR spectral assignments for the hydrated excess proton in liquid water

The local environmental sensitivity of infrared (IR) spectroscopy to a hydrogen-bonding structure makes it a powerful tool for investigating the structure and dynamics of excess protons in water. Although of significant interest, the line broadening that results from the ultrafast evolution of different solvated proton-water structures makes the assignment of liquid-phase IR spectra a challenging task. In this work, we apply a normal mode analysis using density functional theory of thousands of proton-water clusters taken from reactive molecular dynamics trajectories of the latest generation multistate empirical valence bond proton model (MS-EVB 3.2). These calculations are used to obtain a vibrational density of states and IR spectral density, which are decomposed on the basis of solvated proton structure and the frequency dependent mode character. Decompositions are presented on the basis of the proton sharing parameter δ, often used to distinguish Eigen and Zundel species, the stretch and bend character of the modes, the mode delocalization, and the vibrational mode symmetry. We find there is a wide distribution of vibrational frequencies spanning 1200-3000 cm^-1 for every local proton configuration, with the region 2000-2600 cm^-1 being mostly governed by the distorted Eigen-like configuration. We find a continuous red shift of the special-pair O⋯H⁺⋯O stretching frequency, and an increase in the flanking water bending intensity with decreasing δ. Also, we find that the flanking water stretch mode of the Zundel-like species is strongly mixed with the flanking water bend, and the special pair proton oscillation band is strongly coupled with the bend modes of the central H₅O2+moiety.

Deconstructing Prominent Bands in the Terahertz Spectra of H7O3+ and H9O4+: Intermolecular Modes in Eigen Clusters

We report experimental vibrational action spectra (210-2200 cm^-1) and calculated IR spectra, using recent ab initio potential energy and dipole moment surfaces, of H₇O₃⁺ and H₉O₄⁺. We focus on prominent far-IR bands, which postharmonic analyses show, arise from two types of intermolecular motions, i.e., hydrogen bond stretching and terminal water wagging modes, that are similar in both clusters. The good agreement between experiment and theory further validates the accuracy of the potential and dipole moment surfaces, which was used in a recent theoretical study that concluded that infrared photodissociation spectra of the cold clusters correspond to the Eigen isomer. The comparison between theory and experiment adds further confirmation of the need of postharmonic analysis for these clusters.

Communication: Trapping a proton in argon: Spectroscopy and theory of the proton-bound argon dimer and its solvation

Ion-molecule complexes of the form H⁺Ar_n are produced in pulsed-discharge supersonic expansions containing hydrogen and argon. These ions are analyzed and mass-selected in a reflectron spectrometer and studied with infrared laser photodissociation spectroscopy. Infrared spectra for the n = 3-7 complexes are characterized by a series of strong bands in the 900-2200 cm^-1 region. Computational studies at the MP2/aug-cc-pVTZ level examine the structures, binding energies, and infrared spectra for these systems. The core ion responsible for the infrared bands is the proton-bound argon dimer, Ar-H⁺-Ar, which is progressively solvated by the excess argon. Anharmonic vibrational theory is able to reproduce the vibrational structure, identifying it as arising from the asymmetric proton stretch in combination with multiple quanta of the symmetric argon stretch. Successive addition of argon shifts the proton vibration to lower frequencies, as the charge is delocalized over more ligands. The Ar-H⁺-Ar core ion has a first solvation sphere of five argons.

MP4 Study of the Anharmonic Coupling of the Shared Proton Stretching Vibration of the Protonated Water Dimer in Equilibrium and Transition States

The structure and harmonic and anharmonic IR spectra of the protonated water dimer (PWD) were calculated in C₁, C₂, and C_s symmetry at the MP4/acc-pVTZ level of theory. We found that structure and IR spectra are practically identical in C₂ and C₁ symmetry, demonstrating that an equilibrium C₁ configuration of the PWD is not realized. Anharmonic coupling of the shared proton stretching vibration with all other modes in the PWD in C₂ and C_s symmetry was the focus of this investigation. For this purpose, 28 two-dimensional potential energy surfaces (2D PES) were built at the MP4/acc-pVTZ level of theory and the corresponding vibrational Schrödinger equations were solved using the DVR method. Differences in the coupling of the investigated mode with other modes in the C₂ and C_s configurations, along with some factors that determine the red- or blue-shift of the stretching vibration frequency, were analyzed. We obtained a rather reasonable value of the stretching frequency of the bridging proton (1058.4 cm^-1) unperturbed by Fermi resonance. The Fermi resonance between the fundamental vibration ν₇ and the combined vibration ν₂ + ν₆ of the same symmetry was analyzed through anharmonic second-order perturbation theory calculations, as well as by 3D PES constructed using Q₂, Q₆, and Q₇ as normal coordinates. A significant (up to 50%) transfer of intensity from the fundamental vibration to the combined one was found. We have estimated the frequency of the bridging proton stretching vibration in the C_s configuration of the PWD based on calculations of the intrinsic anharmonicity and anharmonic double modes interactions at the MP4/acc-pVTZ level of theory (1261 cm^-1).

Experimental evidence of a 3-centre, 2-electron covalent bond character of the central O–H–O fragment on the Zundel cation in crystals of Zundel nitranilate tetrahydrate

Charge density of the Zundel cation in the solid state reveals a covalent nature of its central O–H–O fragment.

Ab Initio Studies of Structural and Vibrational Properties of Protonated Water Cluster H7O3(+) and Its Deuterium Isotopologues: An Application of Driven Molecular Dynamics

In this work, we present infrared (IR) spectra of H7O3(+) and its deuterium isotopomers calculated by direct molecular dynamics (MD) simulations at the B3LYP/6-31+G** computational level. The calculated spectra obtained at 100, 300, and 500 K were compared to available experimental observations, and spectral features were assigned using normal-mode analysis (NMA) and driven molecular dynamics (DMD). Spectral peaks at 2410 and 2540 cm(-1) were assigned to asymmetric and symmetric stretches of the bridging hydrogen (BH) using NMA. The weak spectral features at 2166 and 2275 cm(-1) were assigned to a combination band of BH asymmetric stretch, H2O in phase wagging, OO stretch, and H3O(+) rocking vibrations by DMD simulations. Our observation of BH stretch vibrations as low as 2166 cm(-1) is in good agreement with the assignment of the low-resolution spectrum obtained by Schwarz at 2200-2300 cm(-1) [Schwarz, H. A. J. Chem. Phys. 1977, 67, 5525-5534] and vibrational predissociation spectrum by Lee et al. ∼2300 cm(-1) [Okumura, M.; Yeh, L. I.; Myers, J. D.; Lee, Y. T. J. Chem. Phys. 1990, 94, 3416-3427].

Driven molecular dynamics studies of the shared proton motion in the H5O2+·Ar cluster: the effect of argon tagging and deuteration on vibrational spectra

We report IR spectra of H5O2(+) and H5O2(+)·Ar and their deuterium isotopologues using ab initio molecular dynamics. The trajectories were propagated as microcanonical (NVE) ensembles at energies corresponding to temperatures 50 and 100 K. The potential energy surface is calculated on-the-fly at the MP2/aug-cc-pVDZ level of theory. The calculations show that adding an argon atom to H5O2(+) introduces symmetry breaking in the Zundel core ion, causes blueshift in the shared proton vibration by about 200 cm(-1), and leads to the splitting of the OH stretch vibrations into four bands. Driven molecular dynamics (DMD) method is used to assign the spectrum by coupling the dipole moment to an external electric field oscillating at frequency ω. The broad feature at 1100 cm(-1) in the H5O2(+)·Ar spectrum is ascribed to the large amplitude shared proton vibration coupled with torsion and wag modes. MD MP2 simulations predict the H/D redshift in the shared proton vibration and water bending vibration to be about 280 and 460 cm(-1), respectively, in good agreement with experimental observations.

Fundamental frequency from classical molecular dynamics

We give a theoretical validation for calculating fundamental frequencies of a molecule from classical molecular dynamics (MD) when its anharmonicity is small enough to be treated by perturbation theory. We specifically give concrete answers to the following questions: (1) What is the appropriate initial condition of classical MD to calculate the fundamental frequency? (2) From that condition, how accurately can we extract fundamental frequencies of a molecule? (3) What is the benefit of using ab initio MD for frequency calculations? Our analytical approaches to those questions are classical and quantum normal form theories. As numerical examples we perform two types of MD to calculate fundamental frequencies of H2O with MP2/aug-cc-pVTZ: one is based on the quartic force field and the other one is direct ab initio MD, where the potential energies and the gradients are calculated on the fly. From those calculations, we show comparisons of the frequencies from MD with the post vibrational self-consistent field calculations, second- and fourth-order perturbation theories, and experiments. We also apply direct ab initio MD to frequency calculations of C-H vibrational modes of tetracene and naphthalene. We conclude that MD can give the same accuracy in fundamental frequency calculation as second-order perturbation theory but the computational cost is lower for large molecules.

A Molecular Dynamic Simulation of Hydrated Proton Transfer in Perfluorosulfonate Ionomer Membranes (Nafion 117)

A molecular dynamic model based on Lennard-Jones Potential, the interaction force between two particles, molecular diffusion, and radial distribution function (RDF) is presented. The diffusion of the hydrated ion, triggered by both Grotthuss and vehicle mechanisms, is used to study the proton transfer in Nafion 117. The hydrated ion transfer mechanisms and the effects of the temperature, the water content in the membrane, and the electric field on the diffusion of the hydrated ion are analyzed. The molecular dynamic simulation results are in good agreement with those reported in the literature. The modeling results show that when the water content in Nafion 117 is low, H3O+is the main transfer ion among the different hydrated ions. However, at higher water content, the hydrated ion in the form of H+(H2O)2is the main transfer ion. It is also found that the negatively charged sulfonic acid group as the fortified point facilitates the proton transfer in Nafion 117 better than the free water molecule. The diffusion of the hydrated ion can be improved by increasing the cell temperature, the water content in Nafion, and the electric field intensity.

Vibrational mode assignment of finite temperature infrared spectra using the AMOEBA polarizable force field

The Driven Molecular Dynamics approach has been adapted and associated with the AMOEBA polarizable force field to assign and visualize vibrational modes in infrared spectra obtained by molecular dynamics simulations.

Nitranilic acid hexahydrate, a novel benchmark system of the Zundel cation in an intrinsically asymmetric environment: spectroscopic features and hydrogen bond dynamics characterised by experimental and theoretical methods

Nitranilic acid (2,5-dihydroxy-3,6-dinitro-2,5-cyclohexadiene-1,4-dione) as a strong dibasic acid in acidic aqueous media creates the Zundel cation, H5O2(+). The structural unit in a crystal comprises (H5O2)2(+) (2,5-dihydroxy-3,6-dinitro-1,4-benzoquinonate)(2-) dihydrate where the Zundel cation reveals no symmetry, being an ideal case for studying proton dynamics and its stability. The Zundel cation and proton transfer dynamics are studied by variable-temperature X-ray diffraction, IR and solid-state NMR spectroscopy, and various quantum chemical methods, including periodic DFT calculations, ab initio molecular dynamics simulation, and quantization of nuclear motion along three fully coupled internal coordinates. The Zundel cation features a short H-bond with the O···O distance of 2.433(2) Å with an asymmetric placement of hydrogen. The proton potential is of a single well type and, due to the non-symmetric surroundings, of asymmetric shape. The formation of the Zundel cation is facilitated by the electronegative NO2 groups. The employed spectroscopic techniques supported by calculations confirm the presence of a short H-bond with a complex proton dynamics.

Probing Small Protonated Water Clusters in Acetonitrile Solutions by 1H NMR

In a previous publication by Kalish et al. (J. Phys. Chem. A 115 (2011) 4063) the existence of well defined small protonated water clusters in acetonitrile has been established by IR spectroscopy. Here we report on a 1H NMR study of triflic acid, CF3SO3H, in acetonitrile-water solutions. Using NMR we are able to corroborate the general solvation scheme we have proposed for the hydrated proton in acetonitrile as a function of the molar ratio between the strong mineral acid and water, n = [H2O]/[acid]. According to this scheme, backed now by both IR absorption spectroscopy and NMR measurements, the very strong triflic acid completely dissociates in acetonitrile/water solutions to yield the aqueous proton and the triflate anion when n > 1. Furthermore, increasing n results in the proton solvated in increasingly larger water clusters formed within the acetonitrile solution.

Clearly distinguishable by NMR are the smallest protonated water clusters, the protonated water monomer, H3 +O, and the protonated water dimer, H5 +O2, which dominate the solution for n = 1,2,3. For larger n the NMR study indicates the gradual increase of the average protonated water cluster size as a function of n while the proton inner solvation core more closely retaining the characteristics of a deformed protonated water dimer, (H2O-H+⋯OH2) s than that of the protonated water monomer (H3 +O) s .

Water-assisted proton transfer in ferredoxin I

The role of water molecules in assisting proton transfer (PT) is investigated for the proton-pumping protein ferredoxin I (FdI) from Azotobacter vinelandii. It was shown previously that individual water molecules can stabilize between Asp(15) and the buried [3Fe-4S](0) cluster and thus can potentially act as a proton relay in transferring H(+) from the protein to the μ(2) sulfur atom. Here, we generalize molecular mechanics with proton transfer to studying proton transfer reactions in the condensed phase. Both umbrella sampling simulations and electronic structure calculations suggest that the PT Asp(15)-COOH + H(2)O + [3Fe-4S](0) → Asp(15)-COO(-) + H(2)O + [3Fe-4S](0) H(+) is concerted, and no stable intermediate hydronium ion (H(3)O(+)) is expected. The free energy difference of 11.7 kcal/mol for the forward reaction is in good agreement with the experimental value (13.3 kcal/mol). For the reverse reaction (Asp(15)-COO(-) + H(2)O + [3Fe-4S](0)H(+) → Asp(15)-COOH + H(2)O + [3Fe-4S](0)), a larger barrier than for the forward reaction is correctly predicted, but it is quantitatively overestimated (23.1 kcal/mol from simulations versus 14.1 from experiment). Possible reasons for this discrepancy are discussed. Compared with the water-assisted process (ΔE ≈ 10 kcal/mol), water-unassisted proton transfer yields a considerably higher barrier of ΔE ≈ 35 kcal/mol.

Apparent stoichiometry of water in proton hydration and proton dehydration reactions in CH3CN/H2O solutions

Gradual solvation of protons by water is observed in liquids by mixing strong mineral acids with various amounts of water in acetonitrile solutions, a process which promotes rapid dissociation of the acids in these solutions. The stoichiometry of the reaction XH(+) + n(H(2)O) = X + (H(2)O)(n)H(+) was studied for strong mineral acids (negatively charged X, X = ClO(4)¯, Cl¯, Br¯, I¯, CF(3)SO(3)¯) and for strong cationic acids (uncharged X, X = R*NH(2), H(2)O). We have found by direct quantitative analysis preference of n = 2 over n = 1 for both groups of proton transfer reactions at relatively low water concentrations in acetonitrile. At high water concentrations, we have found that larger water solvates must also be involved in the solvation of the proton while the spectral features already observed for n = 2, H(+)(H(2)O)(2), remain almost unchanged at large n values up to at least 10 M of water.

Comparison of cationic, anionic and neutral hydrogen bonded dimers

Short Strong Hydrogen Bonds (SSHBs) play an important role in many fields of physics, chemistry and biology. Since it is known that SSHBs exist in many biological systems, the role of hydrogen bonding motifs has been particularly interesting in enzyme catalysis, bio-metabolism, protein folding and proton transport phenomena. To explore the characteristic features of neutral, anionic and cationic hydrogen bonds, we have carried out theoretical studies of diverse homogeneous and heterogeneous hydrogen bonded dimers including water, peroxides, alcohols, ethers, aldehydes, ketones, carboxylic acids, anhydrides, and nitriles. Geometry optimization and harmonic frequency calculations are performed at the levels of Density Functional Theory (DFT) and Møller-Plesset second order perturbation (MP2) theory. First principles Car-Parrinello molecular dynamics (CPMD) simulations are performed to obtain IR spectra derived from velocity- and dipole-autocorrelation functions. We find that the hydrogen bond energy is roughly inversely proportional to the fourth power of the r(O/N-H) distance. Namely, the polarization of the proton accepting O/N atom by the proton-donating H atom reflects most of the binding energy in these diverse cation/anion/neutral hydrogen bonds. The present study gives deeper insight into the nature of hydrogen-bonded dimers including SSHBs.