Deciphering the infrared spectrum of the protonated water pentamer and the hybrid Eigen-Zundel cation

Intra-cluster Charge Migration upon Hydration of Protonated Formic Acid Revealed by Anharmonic Analysis of Cold Ion Vibrational Spectra

The rates of many chemical reactions are accelerated when carried out in micron-sized droplets, but the molecular origin of the rate acceleration remains unclear. One example is the condensation reaction of 1,2-diaminobenzene with formic acid to yield benzimidazole. The observed rate enhancements have been rationalized by invoking enhanced acidity at the surface of methanol solvent droplets with low water content to enable protonation of formic acid to generate a cationic species (protonated formic acid or PFA) formed by attachment of a proton to the neutral acid. Because PFA is the key feature in this reaction mechanism, vibrational spectra of cryogenically cooled, microhydrated PFA·(H2O)n=1-6 were acquired to determine how the extent of charge localization depends on the degree of hydration. Analysis of these highly anharmonic spectra with path integral ab initio molecular dynamics simulations reveals the gradual displacement of the excess proton onto the water network in the microhydration regime at low temperatures with n = 3 as the tipping point for intra-cluster proton transfer.

Magic Numbers and Stabilities of Photoionized Water Clusters: Computational and Experimental Characterization of the Nanosolvated Hydronium Ion

The stability and distributions of small water clusters generated in a supersonic beam expansion are interrogated by tunable vacuum ultraviolet (VUV) radiation generated at a synchrotron. Time-of-flight mass spectrometry reveals enhanced population of various protonated water clusters (H⁺(H₂O)_n) based upon ionization energy and photoionization distance from source, suggesting there are "magic" numbers below the traditional n = 21 that predominates in the literature. These intensity distributions suggest that VUV threshold photoionization (11.0-11.5 eV) of neutral water clusters close to the nozzle exit leads to a different nonequilibrium state compared to a skimmed molecular beam. This results in the appearance of a new magic number at 14. Metadynamics conformer searches coupled with modern density functional calculations are used to identify the global minimum energy structures of protonated water clusters between n = 2 and 21, as well as the manifold of low-lying metastable minima. New lowest energy structures are reported for the cases of n = 5, 6, 11, 12, 16, and 18, and special stability is identified by several measures. These theoretical results are in agreement with the experiments performed in this work in that n = 14 is shown to exhibit additional stability, based on the computed second-order stabilization energy relative to most cluster sizes, though not to the extent of the well-known n = 21 cluster. Other cluster sizes that show some additional energetic stability are n = 7, 9, 12, 17, and 19. To gain insight into the balance between ion-water and water-water interactions as a function of the cluster size, an analysis of the effective two-body interactions (which sum exactly to the total interaction energy) was performed. This analysis reveals a crossover as a function of cluster size between a water-hydronium-dominated regime for small clusters and a water-water-dominated regime for larger clusters around n = 17.

Spectral signatures of excess-proton waiting and transfer-path dynamics in aqueous hydrochloric acid solutions

The theoretical basis for linking spectral signatures of hydrated excess protons with microscopic proton-transfer mechanisms has so far relied on normal-mode analysis. We introduce trajectory-decomposition techniques to analyze the excess-proton dynamics in ab initio molecular-dynamics simulations of aqueous hydrochloric-acid solutions beyond the normal-mode scenario. We show that the actual proton transfer between two water molecules involves for relatively large water-water separations crossing of a free-energy barrier and thus is not a normal mode, rather it is characterized by two non-vibrational time scales: Firstly, the broadly distributed waiting time for transfer to occur with a mean value of 200–300 fs, which leads to a broad and weak shoulder in the absorption spectrum around 100 cm⁻¹, consistent with our experimental THz spectra. Secondly, the mean duration of a transfer event of about 14 fs, which produces a rather well-defined spectral contribution around 1200 cm⁻¹ and agrees in location and width with previous experimental mid-infrared spectra.

The spectroscopic signatures of excess protons in HCl solutions are studied by ab initio simulations and THz experiments. Two contributions beyond the normal-mode scenario are identified that reflect proton-waiting and proton-transfer processes.

Towards complete assignment of the infrared spectrum of the protonated water cluster H+(H2O)21

The spectroscopic features of protonated water species in dilute acid solutions have been long sought after for understanding the microscopic behavior of the proton in water with gas-phase water clusters H⁺(H₂O)_n extensively studied as bottom-up model systems. We present a new protocol for the calculation of the infrared (IR) spectra of complex systems, which combines the fragment-based Coupled Cluster method and anharmonic vibrational quasi-degenerate perturbation theory, and demonstrate its accuracy towards the complete and accurate assignment of the IR spectrum of the H⁺(H₂O)₂₁ cluster. The site-specific IR spectral signatures reveal two distinct structures for the internal and surface four-coordinated water molecules, which are ice-like and liquid-like, respectively. The effect of inter-molecular interaction between water molecules is addressed, and the vibrational resonance is found between the O-H stretching fundamental and the bending overtone of the nearest neighboring water molecule. The revelation of the spectral signature of the excess proton offers deeper insight into the nature of charge accommodation in the extended hydrogen-bonding network underpinning this aqueous cluster.

Protonated water species have been the subject of numerous experimental and computational studies. Here the authors provide a nearly complete assignment of the experimental IR spectrum of the H⁺(H₂O)₂₁ water cluster based on high-level wavefunction theory and anharmonic vibrational quasi-degenerate perturbation theory.

Nature of hydrated proton vibrations revealed by nonlinear spectroscopy of acid water nanodroplets

We use polarization-resolved femtosecond pump-probe spectroscopy to investigate the vibrations of hydrated protons in anionic (AOT) and cationic (CTAB/hexanol) reverse micelles in the frequency range 2000-3500 cm-1. For small AOT micelles the dominant proton hydration structure consists of H3O+ with two OH groups donating hydrogen bonds to water molecules, and one OH group donating a weaker hydrogen bond to sulfonate. For cationic reverse micelles, we find that the absorption at frequencies >2500 cm-1 is dominated by asymmetric proton-hydration structures in which one of the OH groups of H3O+ is more weakly hydrogen-bonded to water than the other two OH groups.

Understanding and Tracking the Excess Proton in Ab Initio Simulations; Insights from IR Spectra

Proton transport in aqueous media is ubiquitously important in chemical and biological processes. Although ab initio molecular dynamics (AIMD) simulations have made great progress in characterizing proton transport, there has been a long-standing challenge in defining and tracking the excess proton, or more properly, the center of excess charge (CEC) created when a hydrogen nucleus distorts the electron distributions of water molecules in a delocalized and highly dynamic nature. Yet, defining (and biasing) such a CEC is essential when combining AIMD with enhanced sampling methods to calculate the relevant macroscopic properties via free-energy landscapes, which is the standard practice for most processes of interest. Several CEC formulas have been proposed and used, but none have yet been systematically tested or rigorously derived. In this paper, we show that the CEC can be used as a computational tool to disentangle IR features of the solvated excess proton from its surrounding solvent, and in turn, how correlating the features in the excess charge spectrum with the behavior of CEC in simulations enables a systematic evaluation of various CEC definitions. We present a new definition of CEC and show how it overcomes the limitations of those currently available both from a spectroscopic point of view and from a practical perspective of performance in enhanced sampling simulations.

Global and local minima of protonated acetonitrile clusters

Potential energy surfaces of protonated acetonitrile clusters have been explored to locate global and local minima energy structures. The structures are stabilized by strong hydrogen bonds, anti-parallel dimers, dipole–dipole and CH⋯N interactions.

Disentangling the Complex Vibrational Mechanics of the Protonated Water Trimer by Rational Control of Its Hydrogen Bonds

The vibrational spectrum of the protonated water trimer, H⁺(H₂O)₃, is surprisingly complex, with many strong features in the expected region of the fundamentals associated with two H-bonded OH groups on the H₃O⁺ core ion. Here we follow how the bands in this region of the spectrum evolve when the energies of the fundamentals in the H-bonded OH stretches are systematically increased by the attachment of increasingly strongly bound "tag" molecules (He, Ar, D₂, N₂, CO, and H₂O) to the free OH position on the hydronium core ion of H⁺(H₂O)₃, as well as by replacement of the hydrogen atom in the nonbonded OH group on hydronium with methyl and ethyl groups. This allows for the incremental transformation of the complex band pattern observed in H⁺(H₂O)₃ into that of the "Eigen" structure of the protonated water tetramer. Differences among the trajectories of the various bands provide an empirical way to disentangle features primarily due to the displacements of the OH stretches bound to the hydronium core from those arising from anharmonic coupling to states involving one or more quanta in lower frequency modes. The latter are found to be dramatically enhanced when the nominal frequencies of the intermolecular OH stretching modes approach those of the intramolecular bends of the H₃O⁺ and H₂O constituents in both H and D isotopologues.

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.

Molecular modeling and assignment of IR spectra of the hydrated excess proton in isotopically dilute water

Infrared (IR) spectroscopy of the water O-H stretch has been widely used to probe both the local hydrogen-bonding structure and dynamics of aqueous systems. Although of significant interest, the IR spectroscopy of excess protons in water remains difficult to assign as a result of extensive and strong intermolecular interactions in hydrated proton complexes. As an alternate approach, we develop a mixed quantum-classical model for the vibrational spectroscopy of the excess proton in isotopically dilute water that draws on frozen proton-water clusters taken from reactive molecular dynamics trajectories of the latest generation multi-state empirical valence bond proton model (MS-EVB 3.2). A semi-empirical single oscillator spectroscopic map for the instantaneous transition frequency and transition dipole moment is constructed using potential energy surfaces for the O-H stretch coordinate of the excess proton using electronic structure calculations. Calculated spectra are compared with experimental spectra of dilute H⁺ in D₂O obtained from double-difference FTIR to demonstrate the validity of the map. The model is also used to decompose IR spectra into contributions from different aqueous proton configurations. We find that the O-H transition frequency continuously decreases as the oxygen-oxygen length for a special pair proton decreases, shifting from Eigen- to Zundel-like configurations. The same shift is accompanied by a shift of the flanking water stretches of the Zundel complex to higher frequency than the hydronium O-H vibrations.

DFT Study of Hydrogen-Bonding Interaction, Solvation Effect, and Electric-Field Effect on Raman Spectra of Hydrated Proton

Strong hydrogen-bonding interaction and Raman spectra of hydrated proton have been investigated using hybrid density functional theory method B3LYP. The solvation model of density (SMD) approach is employed in the present calculation to simulate hydrated protons in aqueous solution. Focusing on the different hydrogen-bonded Eigen-water and Zundel-water interactions, we present a better assignment of Raman signals of hydrated proton on the basis of vibrational analysis in different environments. Our results showed that B3LYP calculations could give a good prediction for characteristic vibrational frequencies of Eigen and Zundel isomers in liquid phase. The O-H stretching vibrational frequencies from Eigen and Zundel units are very sensitive to hydrogen-bonding interaction with solvent water molecules. Moreover, the solvation effect and the external electric-field effect lead to the proton deviating from the central position of Zundel structure and finally resulting in a transition to Eigen one in aqueous solution. Furthermore, by combining theoretical prediction and Raman scattering theory, we calculate absolute Raman intensities of characteristic signals based on the polarizability tensor derivatives of hydrated proton clusters. This is very helpful to infer the microstructure of hydrated protons in aqueous solution by using Raman measurements.

Protons and Hydroxide Ions in Aqueous Systems

Understanding the structure and dynamics of water's constituent ions, proton and hydroxide, has been a subject of numerous experimental and theoretical studies over the last century. Besides their obvious importance in acid-base chemistry, these ions play an important role in numerous applications ranging from enzyme catalysis to environmental chemistry. Despite a long history of research, many fundamental issues regarding their properties continue to be an active area of research. Here, we provide a review of the experimental and theoretical advances made in the last several decades in understanding the structure, dynamics, and transport of the proton and hydroxide ions in different aqueous environments, ranging from water clusters to the bulk liquid and its interfaces with hydrophobic surfaces. The propensity of these ions to accumulate at hydrophobic surfaces has been a subject of intense debate, and we highlight the open issues and challenges in this area. Biological applications reviewed include proton transport along the hydration layer of various membranes and through channel proteins, problems that are at the core of cellular bioenergetics.

Gas phase vibrational spectroscopy of the protonated water pentamer: the role of isomers and nuclear quantum effects

We use cryogenic ion trap vibrational spectroscopy to study the structure of the protonated water pentamer, H⁺(H₂O)₅, and its fully deuterated isotopologue, D⁺(D₂O)₅.

Structures and spectroscopy of protonated ammonia clusters at different temperatures

Protonated ammonia clusters are all Eigen structures and the first solvation shell of the related ammonium ion core is saturated by four ammonia molecules.

Disentangling the Contribution of Multiple Isomers to the Infrared Spectrum of the Protonated Water Heptamer

We use infrared/infrared double-resonance population labeling (IR(2)MS(2)) spectroscopy in the spectral region of the free and hydrogen-bonded OH stretching fundamentals (2880-3850 cm(-1)) to identify the number and to isolate the vibrational signatures of individual isomers contributing to the gas-phase IR spectra of the cryogenically cooled protonated water clusters H(+)(H2O)n·H2/D2 with n = 7-10. For n = 7, four isomers are identified and assigned. Surprisingly, the IR(2)MS(2) spectra of the protonated water octa-, nona-, and decamer show no evidence for multiple isomers. The present spectra support the prediction that the quasi-2D to 3D structural transition occurs in between n = 8 and 9 in the cold cluster regime. However, the same models have difficulty explaining the remarkable size dependence of the isomer population reported here.