HCl hydrates as model systems for protonated water

Ab initio molecular dynamics simulations are presented of vibrational dynamics and spectra of crystal HCl hydrates. Depending on the composition, the hydrates include distinct protonated water forms, which in their equilibrium structures approximate either the Eigen ion H3O+(H2O)3 (in the hexahydrate) or the Zundel H2O...H+...OH2 ion (in the di- and trihydrate). Thus, the hydrates offer the opportunity to study spectra and dynamics of distinct species of protonated water trapped in a semirigid solvating environment. The experimentally measured spectra are reproduced quite well by BLYP/DZVP-level calculations employing Fourier transform of the system dipole. The large overall width (800-1000 cm-1) of structured proton bands reflects a broad range of solvating environments generated by crystal vibrations. The aqueous HCl solution was also examined in search of an objective criterion for separating the contributions of "Zundel-like" and "Eigen-like" protonated forms. It is suggested that no such criterion exists since distributions of proton-related structural properties appear continuous and unimodal. Dipole derivatives with respect to OH and O...H+ stretches in water and protonated water were also investigated to advance the understanding of the corresponding IR intensities. The effects of H bonding and solvation on the intensities were analyzed with the help of the Wannier centers' representation of electron density.     

Eigen and Zundel forms of small protonated water clusters: structures and infrared spectra

The spectral properties of protonated water clusters, especially the difference between Eigen (H3O+) and Zundel (H5O2+) conformers and the difference between their unhydrated and dominant hydrated forms are investigated with the first principles molecular dynamics simulations as well as with the high level ab initio calculations. The vibrational modes of the excess proton in H3O+ are sensitive to the hydration, while those in H5O2+ are sensitive to the messenger atom such as Ar (which was assumed to be weakly bound to the water cluster during acquisitions of experimental spectra). The spectral feature around approximately 2700 cm-1 (experimental value: 2665 cm-1) for the Eigen moiety appears when H3O+ is hydrated. This feature corresponds to the hydrating water interacting with H3O+, so it cannot appear in the Eigen core. Thus, H3O+ alone would be somewhat different from the Eigen forms in water. For the Zundel form (in particular, H5O2+), there have been some differences in spectral features among different experiments as well as between experiments and theory. When an Ar messenger atom is introduced at a specific temperature corresponding to the experimental condition, the calculated vibrational spectra for H5O2+.Ar are in good agreement with the experimental infrared spectra showing the characteristic Zundel frequency at approximately 1770 cm-1. Thus, the effect of hydration, messenger atom Ar, and temperature are crucial to elucidating the nature of vibrational spectra of Eigen and Zundel forms and to assigning the vibrational modes of small protonated water clusters.     

HCl solvation at the surface and within methanol clusters/nanoparticles II: Evidence for molecular wires

Condensed-phase solvation of HCl on and within methanol nanoparticles was investigated by Fourier transform infrared (FTIR) spectroscopy, on-the-fly molecular dynamics as implemented in the density functional code Quickstep (which is part of the CP2K package), and ab initio calculations. Adsorption and solvation stages are identified and assigned with the help of calculated infrared spectra obtained from the simulations. The results have been further checked with MP2-level ab initio calculations. The range of acid solvation states extends from the single-coordinated slightly stretched HCl to proton-sharing with Zundel-like methanol O...H+...X- states, and finally to MeOH2+...Cl- units with full proton transfer. Furthermore, once the proton moves to methanol, it is mobilized along methanol molecular chains. Since the proton dynamics reflects the evolving local structures, the "proton" spectra display broad bands usually with underlying continua.     

Ab initio molecular dynamics simulation of the energy-relaxation process of the protonated water dimer

Relaxation processes of the energy-rich protonated water dimer H+(H2O)2 were investigated by the ab initio molecular dynamics (AIMD) method. At first, the energy-rich H+(H2O)2 was reproduced by simulating a collision reaction between the protonated water monomer H3O+ and H2O. Next it was collided with N2 in order to observe the effects of intramolecular vibration redistribution and intermolecular energy transfer. Forty-eight AIMD simulations of the collision of H+(H2O)2 with N2 were performed by changing the initial orientation and the time interval between two collisions. It was revealed that the amount of energy transferred from H+(H2O)2 to N2 decreased the longer the time interval. The relationship between the intermolecular energy transfer and the vibrational states was examined with the use of an energy-transfer spectrogram (ETS), which is an analysis technique combining energy density analysis and short-time Fourier transform. The ETS demonstrates a characteristic vibrational mode for the energy transfer, which corresponds to the stretching of the hydrogen bond between H+(H2O)2 and N2 in an active complex.     

Infrared photodissociation spectroscopy of protonated acetylene and its clusters

The protonated acetylene cation, C2H3+, (also known as the vinyl cation) and the proton-bound acetylene dimer cation (C4H5+) are produced by a pulsed supersonic nozzle/pulsed electrical discharge cluster source. The parent ions are also generated with weakly attached argon "tag" atoms, e.g., C2H3+Ar and C4H5+Ar. These ions are mass selected in a specially designed reflectron time-of-flight mass spectrometer and studied with infrared laser photodissociation spectroscopy in the 800-3600 cm-1 region. Vibrational resonances are detected for both ions in the C-H stretching region. C2H3+ has a strong vibrational resonance near 2200 cm-1 assigned to the bridged proton stretch of the nonclassical ion, while C4H5+ has no such free-proton vibration. Instead, C4H5+ has resonances near 1300 cm-1, consistent with a symmetrically shared proton in a di-bridged structure. Although the shared proton structure is not the lowest energy isomer of C4H5+, this species is apparently stabilized under the supersonic beam conditions. Larger clusters containing additional acetylene units are also investigated via the elimination of acetylene. These species have new IR bands indicating that rearrangement reactions have taken place to produce core C4H5+ ions with the methyl cyclopropane cation structure and/or the protonated cyclobutadiene isomer. Ab initio (MP2) calculations provide structures and predicted spectra consistent with all of these experiments.     

Ab initio molecular dynamics of protonated dialanine and comparison to infrared multiphoton dissociation experiments

Finite temperature Car-Parrinello molecular dynamics simulations are performed for the protonated dialanine peptide in vacuo, in relation to infrared multiphoton dissociation experiments. The simulations emphasize the flexibility of the different torsional angles at room temperature and the dynamical exchange between different conformers which were previously identified as stable at 0 K. A proton transfer occurring spontaneously at the N-terminal side is also observed and characterized. The theoretical infrared absorption spectrum is computed from the dipole time correlation function, and, in contrast to traditional static electronic structure calculations, it accounts directly for anharmonic and finite temperature effects. The comparison to the experimental infrared multiphoton dissociation spectrum turns out very good in terms of both band positions and band shapes. It does help the identification of a predominant conformer and the attribution of the different bands. The synergy shown between the experimental and theoretical approaches opens the door to the study of the vibrational properties of complex and floppy biomolecules in the gas phase at finite temperature.     

Compatibility between methanol and water in the three-dimensional cage formation of large-sized protonated methanol-water mixed clusters

Infrared spectra of large-sized protonated methanol-water mixed clusters, H(+)(MeOH)(m)(H(2)O)(n) (m=1-4, n=4-22), were measured in the OH stretch region. The free OH stretch bands of the water moiety converged to a single peak due to the three-coordinated sites at the sizes of m+n=21, which is the magic number of the protonated water cluster. This is a spectroscopic signature for the formation of the three-dimensional cage structure in the mixed cluster, and it demonstrates the compatibility of a small number of methanol molecules with water in the hydrogen-bonded cage formation. Density functional theory calculations were carried out to examine the relative stability and structures of selected isomers of the mixed clusters. The calculation results supported the microscopic compatibility of methanol and water in the hydrogen-bonded cage development. The authors also found that in the magic number clusters, the surface protonated sites are energetically favored over their internal counterparts and the excess proton prefers to take the form of H(3)O(+) despite the fact that the proton affinity of methanol is greater than that of water.     

Vibrational analysis of the H5O2+ infrared spectrum using molecular and driven molecular dynamics

Standard molecular and driven molecular dynamics are used to analyze prominent spectral features in the H5O2+ infrared spectrum. In the driven method, the molecular Hamiltonian is augmented with a time-dependent term, mu x epsilon(0) sin(omegat), where mu is the dipole moment of H5O2+, epsilon0 is the electric field, and omega is the frequency. The magnitude of the electric field determines whether the driving is mild (the harmonic limit) or strong (anharmonic motion and mode coupling). We analyze the spectrum in the wavenumber range from 600 to 1900 cm(-1), where recent experimental measurements are available for H5O2+. On the basis of the simulations, we have assigned the broad feature around 1000 cm(-1) to the proton transfer coupled with the torsion motion. Intense absorption near 1780 cm(-1) is assigned to the H2O monomer bend coupled with proton transfer.     

How Does a Membrane Protein Achieve a Vectorial Proton Transfer Via Water Molecules?

We present a detailed mechanism for the proton transfer from a protein‐bound protonated water cluster to the bulk water directed by protein side chains in the membrane protein bacteriorhodopsin. We use a combined approach of time‐resolved Fourier transform infrared spectroscopy, molecular dynamics simulations, and X‐ray structure analysis to elucidate the functional role of a hydrogen bond between Ser193 and Glu204. These two residues seal the internal protonated water cluster from the bulk water and the protein surface. During the photocycle of bacteriorhodopsin, a transient protonation of Glu204 leads to a breaking of this hydrogen bond. This breaking opens the gate to the extracellular bulk water, leading to a subsequent proton release from the protonated water cluster. We show in detail how the protein achieves vectorial proton transfer via protonated water clusters in contrast to random proton transfer in liquid water.     

Comprehensive analysis of the hydrogen bond network morphology and OH stretching vibrations in protonated methanol-water mixed clusters, H(+)(MeOH)1(H2O)n (n = 1-8)

Density functional theory (DFT) calculations of protonated methanol-water mixed clusters, H (+)(MeOH) 1(H 2O) n ( n = 1-8), were extensively carried out to analyze the hydrogen bond structures of the clusters. Various structural isomers were energy optimized, and their relative energies with zero point energy corrections and temperature dependence of the free energies were examined. Coexistence of different morphological isomers was suggested. Infrared spectra were simulated on the basis of the optimized structures. The infrared spectra were also experimentally measured for n = 3-9 in the OH stretching vibrational region. The observed broad bands in the hydrogen-bonded OH stretch region were assigned in comparison with the simulations. From the DFT calculations, the preferential proton location was also investigated. Clear correlations between the excess proton location and the cluster morphology were found.     

Gas-phase reactions of protonated tryptophan

The gas phase reactions of protonated tryptophan have been examined in a quadrupole ion trap using a combination of collision induced dissociation, hydrogen-deuterium exchange, regiospecific deuterium labeling and molecular orbital calculations (at the B3LYP/6-31G* level of theory). The loss of ammonia from protonated tryptophan is observed as the primary fragmentation pathway, with concomitant formation of a [M + H - NH(3)](+) ion by nucleophilic attack from the C3 position of the indole side chain. Hydrogen-deuterium exchange and regiospecific deuterium labeling reveals that scrambling of protons in the C2 and C4 positions of the indole ring, via intramolecular proton transfer from the thermodynamically preferred site of protonation at the amino nitrogen, precedes ammonia loss. Molecular orbital calculations have been employed to demonstrate that the activation barriers to intramolecular proton transfer are lower than that for NH(3) loss.     

Reactive force fields for proton transfer dynamics

A force field-inspired method based on fitted, high-quality multidimensional potential energy surfaces to follow proton transfer (PT) reactions in molecular dynamics simulations is presented. In molecular mechanics with proton transfer (MMPT) a system is partitioned into a region where proton transfer takes place and the remaining degrees of freedom which are treated with a conventional force field. The implementation of the method and applications to specific chemically and biologically relevant scenarios are presented. MMPT is developed in view of two primary areas in mind: to follow the molecular dynamics of proton transfer in the condensed phase on realistic time scales and to adapt the shape (morphing) of the potential energy surface for specific applications. MMPT is applied to PT in protonated ammonia dimer, double proton transfer in 2-pyridone-2-hydroxypyridine, and the first step of PT from a protein side-chain towards a buried [3Fe4S] cluster in ferredoxin I. Specific findings of the work include the fundamental role of the N-N vibration as the gating mode for PT in NH4+...NH3 and the qualitative understanding of PT from the protein to a metastable active-site water molecule in Ferredoxin I.     

Proton mobility and main fragmentation pathways of protonated lysylglycine

Theoretical model calculations were performed to validate the 'mobile proton' model for protonated lysylglycine (KG). Detailed scans carried out at various quantum chemical levels of the potential energy surface (PES) of protonated KG resulted in a large number of minima belonging to various protonation sites and conformers. Transition structures corresponding to proton transfer reactions between different protonation sites were determined, to obtain some energetic and structural insight into the atomic details of these processes. The rate coefficients of the proton transfer reactions between the isomers were calculated using the Rice-Ramsperger-Kassel-Marcus (RRKM) method in order to obtain a quantitative measure of the time-scale of these processes. Our results clearly indicate that the added proton is less mobile for protonated KG than for peptides lacking a basic amino acid residue. However, the energy needed to reach the energetically less favorable but-from the point of view of backbone fragmentation-critical amide nitrogen protonation sites is available in tandem mass spectrometers operated under low-energy collision conditions. Using the results of our scan of the PES of protonated KG, the dissociation pathways corresponding to the main fragmentation channels for protonated KG were also determined. Such pathways include loss of ammonia and formation of a protonated alpha-amino-epsilon-caprolactam. The results of our theoretical modeling, which revealed all the atomic details of these processes, are in agreement with the available experimental results.     

Mechanisms of proton transfer in Nafion: elementary reactions at the sulfonic acid groups

Proton transfer reactions at the sulfonic acid groups in Nafion were theoretically studied, using complexes formed from triflic acid (CF3SO3H), H3O+ and H2O, as model systems. The investigations began with searching for potential precursors and transition states at low hydration levels, using the test-particle model (T-model), density functional theory (DFT) and ab initio calculations. They were employed as starting configurations in Born-Oppenheimer molecular dynamics (BOMD) simulations at 298 K, from which elementary reactions were analyzed and categorized. For the H3O+-H2O complexes, BOMD simulations suggested that a quasi-dynamic equilibrium could be established between the Eigen and Zundel complexes, and that was considered to be one of the most important elementary reactions in the proton transfer process. The average lifetime of H3O+ obtained from BOMD simulations is close to the lowest limit, estimated from low-frequency vibrational spectroscopy. It was demonstrated that proton transfer reactions at -SO3H are not concerted, due to the thermal energy fluctuation and the existence of various quasi-dynamic equilibria, and -SO3H could directly and indirectly mediate proton transfer reactions through the formation of proton defects, as well as the -SO3- and -SO3H2+ transition states.     

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.     

Gas-phase nazarov cyclization of protonated 2-methoxy and 2-hydroxychalcone: An example of intramolecular proton-transport catalysis

Upon CA, ESI generated [M + H]⁺ ions of chalcone (benzalacetophenone) and 3-phenyl-indanone both undergo losses of H₂O, CO, and the elements of benzene. CA of the [M + H]⁺ ions of 2-methoxy and 2-hydroxychalcone, however, prompts instead a dominant loss of ketene. In addition, CA of the [M + H]⁺ ions of 2-methoxy-β-methylchalcone produces an analogous loss of methylketene instead. Furthermore, the [M + D]⁺ ion of 2-methoxychalcone upon CA eliminates only unlabeled ketene, and the resultant product, the [M + D − ketene]⁺ ion, yields only the benzyl-d ₁ cation upon CA. We propose that the 2-methoxy and 2-hydroxy (ortho) substituents facilitate a Nazarov cyclization to the corresponding protonated 3-aryl-indanones by mediating a critical proton transfer. The resultant protonated indanones then undergo a second proton transport catalysis facilitated by the same ortho substituents producing intermediates that eliminate ketene to yield 2-methoxy- or 2-hydroxyphenyl-phenyl-methylcarbocations, respectively. The basicity of the ortho substituent is important; for example, replacement of the ortho function with a chloro substituent does not provide an efficient catalyst for the proton transports. The Nazarov cyclization must compete with an alternate cyclization, driven by the protonated carbonyl group of the chalcone that results in losses of H₂O and CO. The assisted proton transfer mediated by the ortho substituent shifts the competition in favor of the Nazarov cyclization. The proposed mechanisms for cyclization and fragmentation are supported by high-mass resolving power data, tandem mass spectra, deuterium labeling, and molecular orbital calculations.     

FTIR adsorption studies of H2O and CH3OH in the isostructural H-SSZ-13 and H-SAPO-34: formation of H-bonded adducts and protonated clusters

The acidity of the isostructural H-SSZ-13 and H-SAPO-34 has been investigated by transmission FTIR spectroscopy using H2O and CH3OH as molecular probes. Interactions between the zeolitic samples and the probe molecules led to perturbations and proton transfers directly related to the acidity of the materials. The entire set of acidic sites in H-SSZ-13 interacts with H2O and CH3OH to give H-bonded adducts or protonated species. H3O+ is not formed in appreciable amounts upon H2O adsorption on H-SSZ-13, but at high coverages H2O generates clusters that have a proton affinity sufficiently high to abstract protons from the zeolite framework. Parallel experiments carried out for H-SAPO-34 showed that the H2O clusters abstract protons from Br?nsted sites only to a minor extent. Moving to CH3OH, even if it has a higher proton affinity than H2O and should expectingly experience an easier protonation, proton transfer is totally absent in H-SAPO-34 under our set of conditions. The clear evidence of methanol protonation in H-SSZ-13 definitely states the strong acidic character of this material. When irreversibly adsorbed CH3OH is present in H-SSZ-13, an appreciable amount of (CH3)2O is formed upon heating to 573 K. Compared to its SAPO analogue, the present set of data indisputably points to H-SSZ-13 as the strongest Br?nsted acidic material.     

Infrared spectra of homogeneous and heterogeneous proton-bound dimers in the gas phase

Infrared multiphoton dissociation spectra of three homogeneous and two heterogeneous proton-bound dimers were recorded in the gas phase. Comparison of the experimental infrared spectra recorded in the fingerprint region of the proton-bound dimers with spectra predicted by electronic structure calculations shows that all modes which are observed contain motion of the proton oscillating between the two monomers. The O-H-O asymmetric stretch for the homogeneous dimers is shown to occur at around 800 cm-1. As expected, the O-H-O asymmetric stretching modes for the heterogeneous proton-bound dimers are observed to shift to significantly higher energy with respect to those for the homogeneous proton-bound dimers due to the asymmetry of the O-H-O moeity. This shift is shown to be predictable from the difference in proton affinities between the two monomers. Density functional predictions of the infrared spectra based on the harmonic oscillator model are demonstrated to predict the observed spectra of the homogeneous proton-bound dimers with reasonable accuracy. Calculations of the structure and infrared spectrum of protonated diglyme at the B3LYP/6-31+G** level and basis also agree well with an infrared spectrum recorded previously. For both heterogeneous proton-bound dimers, however, the predicted spectra are blue-shifted with respect to experiment.     

Intramolecular condensation reactions in protonated dipeptides: Carbon monoxide,water, and ammonia losses in competition

The elimination of carbon monoxide and water from a series of protonated dipeptides, [XxxYyy + H](+), is investigated by tandem mass spectrometry experiments and density functional theory. The combined results show that CO loss occurs on the a(1)-y(1) pathway, which begins by rearrangement of the added proton to the amide N-atom and creates the proton-bound dimer of an amino acid (Yyy) and an imine (that from Xxx residue). The loss of H(2)O is initiated from a tautomer in which the added proton has migrated to the hydroxyl group of the C-terminus, thereby promoting the formation of an ion with protonated oxazolone structure (a nominal b(2) ion). The highest yields of [XxxYyy + H - CO](+) and [XxxYyy + H - H(2)O](+) are observed at threshold energies. As the internal energy of the protonated dipeptides increases, these primary products are depleted by consecutive dissociations yielding mostly backbone fragments. Specifically, [XxxYyy + H - CO](+) decomposes to y(1) (protonated Yyy) and a(1) (immonium ion of Xxx residue), while [XxxYyy + H - H(2)O](+) produces a(2) and the immonium ions of residues Xxx (a(1)) and Yyy ("internal" immonium ion). Water loss takes place more efficiently when the more basic residue is at the C-terminal position. Increasing the basicity of the N-terminal residue enhances the extent of CO versus H(2)O loss and introduces the competitive elimination of NH(3). The dissociations leading to eliminations of small neutrals (CO, H(2)O, etc.) generally proceed over transition states that lie higher in energy than the corresponding dissociation products. The excess energy is disposed of either in translational or rovibrational modes of the products, depending on the stability of the incipient noncovalent assemblies emerging during the cleavage of the small neutrals.