Structure and vibrational spectra of H[sup +](H[sub 2]O)[sub 8]: Is the excess proton in a symmetrical hydrogen bond?

Electron binding energies and Dyson orbitals of OnH2n+1 +,0,- clusters: Double Rydberg anions, Rydberg radicals, and micro-solvated hydronium cations

Ab initio electron propagator methods are employed to predict the vertical electron attachment energies (VEAEs) of OH₃ ⁺(H₂O)_n clusters. The VEAEs decrease with increasing n, and the corresponding Dyson orbitals are diffused over exterior, non-hydrogen bonded protons. Clusters formed from OH₃ ^- double Rydberg anions (DRAs) and stabilized by hydrogen bonding or electrostatic interactions between ions and polar molecules are studied through calculations on OH₃ ^-(H₂O)_n complexes and are compared with more stable H^-(H₂O)_n+1 isomers. Remarkable changes in the geometry of the anionic hydronium-water clusters with respect to their cationic counterparts occur. Rydberg electrons in the uncharged and anionic clusters are held near the exterior protons of the water network. For all values of n, the anion-water complex H^-(H₂O)_n+1 is always the most stable, with large vertical electron detachment energies (VEDEs). OH₃ ^-(H₂O)_n DRA isomers have well separated VEDEs and may be visible in anion photoelectron spectra. Corresponding Dyson orbitals occupy regions beyond the peripheral O-H bonds and differ significantly from those obtained for the VEAEs of the cations.

Ab Initio Study of Ionized Water Radical Cation (H2O)8+ in Combination with the Particle Swarm Optimization Method

The structures of cationic water clusters (H₂O)₈⁺ have been globally explored by the particle swarm optimization method in combination with quantum chemical calculations. Geometry optimization and vibrational analysis for the 15 most interesting clusters were computed at the MP2/aug-cc-pVDZ level and infrared spectrum calculation at MPW1K/6-311++G** level. Special attention was paid to the relationships between their configurations and energies. Both MP2 and B3LYP-D3 calculations revealed that the cage-like structure is the most stable, which is different from a five-membered ring lowest energy structure but agrees well with a cage-like structure in the literature. Furthermore, our obtained cage-like structure is more stable by 0.87 and 1.23 kcal/mol than the previously reported structures at MP2 and B3LYP-D3 levels, respectively. Interestingly, on the basis of their relative Gibbs free energies and the temperature dependence of populations, the cage-like structure predominates only at very low temperatures, and the most dominating species transforms into a newfound four-membered ring structure from 100 to 400 K, which can contribute greatly to the experimental infrared spectrum. By topological analysis and reduced density gradient analysis, we also investigated the structural characteristics and bonding strengths of these water cluster radical cations.

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(2)O)(n)H(+) , 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.

Role of charge transfer in the structure and dynamics of the hydrated proton

Although it has long been recognized that multiple water molecules strongly associate with an extra proton in bulk water, some models and conceptual frameworks continue to utilize the classical hydronium ion (H(3)O(+)) as a fundamental building block. In this work, the nature of the hydronium ion in aqueous systems is examined using an ab initio energy decomposition analysis (EDA) that evaluates both the magnitude of and energetic stabilization due to charge transfer among H(3)O(+) and the surrounding water molecules. The EDA is performed on structures extracted from dynamical bulk-phase simulations and used to determine how frequently the pure hydronium ion, where the excess charge is primarily localized on H(3)O(+), occurs under dynamic conditions. The answer is essentially never. The energetic stabilization of H(3)O(+) due to charge delocalization to neighboring water molecules is found to be much larger (16-49 kcal/mol) than for other ions (even Li(+)) and to constitute a substantial portion (20-52%) of the complex's binding energy. The charge defect is also shown to have intrinsic dynamical asymmetry and to display dynamical signatures that can be related to features appearing in IR spectra.

Quantum Chemistry Study of H+(H2O)8: A Global Search for Its Isomers by the Scaled Hypersphere Search Method, and Its Thermal Behavior

The structures of the protonated water cluster H+(H2O)8 have been globally explored by the scaled hypersphere search method. On the Hartree-Fock potential energy surface 174 isomers were found, among which 168 were computed to be minima at the B3LYP/6-31+G** level, and their energies were further refined at the level of MP2/6-311++G(3df,2p). The global minimum on the potential energy surface computed at the B3LYP/6-31+G** level shows a cagelike structure with the "Eigen" motif, while the lowest-free-energy isomer has a five-membered-ring structure at 170 K and a chain form at 273 K. The present results are well in line with previous experimental findings. In addition, the ADMP (atom-centered density matrix propagation) simulation indicates that the extra proton in the lowest-free-energy isomer (170 K), which has a five-membered ring and the "Zundel" feature, is often in an asymmetrical hydrogen bond.

Exploration on Regulating Factors for Proton Transfer along Hydrogen-Bonded Water Chains

Proton transfer along a single-file hydrogen-bonded water chain is elucidated with a special emphasis on the investigation of chain length, side water, and solvent effects, as well as the temperature and pressure dependences. The number of water molecules in the chain varies from one to nine. The proton can be transported to the acceptor fragment through the single-file hydrogen-bonded water wire which contains at most five water molecules. If the number of water molecule is more than five, the proton is trapped by the chain in the hydroxyl-centered H(7)O(3) (+) state. The farthest water molecule involved in the formation of H(7)O(3) (+) is the fifth one away from the donor fragment. These phenomena reappear in the molecular dynamics simulations. The energy of the system is reduced along with the proton conduction. The proton transfer mechanism can be altered by excess proton. The augmentation of the solvent dielectric constant weakens the stability of the system, but favors the proton transfer. NMR spin-spin coupling constants can be used as a criterion in judging whether the proton is transferred or not. The enhancement of temperature increases the thermal motion of the molecule, augments the internal energy of the system, and favors the proton transfer. The lengthening of the water wire increases the entropy of the system, concomitantly, the temperature dependence of the Gibbs free energy increases. The most favorable condition for the proton transfer along the H-bonded water wire is the four-water contained chain with side water attached near to the acceptor fragment in polar solvent under higher temperature.

Computational studies of aqueous-phase photochemistry and the hydrated electron in finite-size clusters

A survey of recent ab initio calculations on excited electronic states of water clusters and various chromophore-water clusters is given. Electron and proton transfer processes in these systems have been characterized by the determination of electronic wave functions, minimum-energy reaction paths and potential-energy profiles. It is pointed out that the transfer of a neutral hydrogen atom (leading to biradicals) rather than the transfer of a proton (leading to ion pairs) is the generic excited-state reaction mechanism in these systems. The hydrated hydronium radical, (H3O)(aq), plays a central role in this scenario. The electronic and vibrational spectra of H3O(H2O)(n) clusters and the decay mechanism of these metastable species have been investigated in some detail. The results suggest that (H3O)(aq) could be the carrier of the characteristic spectroscopic properties of the hydrated electron in liquid water.

Cluster Model for the Ionic Product of Water: Accuracy and Limitations of Common Density Functional Methods

In the present study, the performance of six popular density functionals (B3LYP, PBE0, BLYP, BP86, PBE, and SVWN) for the description of the autoionization process in the water octamer was studied. As a benchmark, MP2 energies with complete basis sets limit extrapolation and CCSD(T) correction were used. At this level, the autoionized structure lies 28.5 kcal.mol(-1) above the neutral water octamer. Accounting for zero-point energy lowers this value by 3.0 kcal.mol(-1). The transition state of the proton transfer reaction, lying only 0.7 kcal.mol(-1) above the energy of the ionized system, was identified at the MP2/aug-cc-pVDZ level of theory. Different density functionals describe the reactant and product with varying accuracy, while they all fail to characterize the transition state. We find improved results with hybrid functionals compared to the gradient-corrected ones. In particular, B3LYP describes the reaction energetics within 2.5 kcal.mol(-1) of the benchmark value. Therefore, this functional is suggested to be preferably used both for Carr-Parinello molecular dynamics and for quantum mechanics/molecular mechanics (QM/MM) simulations of autoionization of water.

Hydration and Dissociation of Hydrogen Fluoric Acid (HF)

The hydration and dissociation phenomena of HF(H(2)O)(n)() (n < or = 10) clusters have been studied by using both the density functional theory with the 6-311++G[sp] basis set and the Møller-Plesset second-order perturbation theory with the aug-cc-pVDZ+(2s2p/2s) basis set. The structures for n > or = 8 are first reported here. The dissociated form of the hydrogen-fluoric acid in HF(H(2)O)(n) clusters is found to be less stable at 0 K than the undissociated form until n = 10. HF may not be dissociated at 0 K solely by water molecules because the HF H bond is stronger than the OH H bond, against the expectation that the dissociated HF(H(2)O)(n) would be more stable than the undissociated one in the presence of a number of water molecules. The dissociation would be possible for only a fraction of a number of hydrated HF clusters by the Boltzmann distribution at finite temperatures. This is in sharp contrast to other hydrogen halide acids (HCl, HBr, HI) showing the dissociation phenomena at 0 K for n > or = 4. The IR spectra of dissociated and undissociated structures of HF(H(2)O)(n) are compared. The structures and binding energies of HF(H(2)O)(n) are found to be similar to those of (H(2)O)(n+1). It is interesting that HF(H(2)O)(n=5,6,10) are slightly less stable compared with other sizes of clusters, just like the fact that (H(2)O)(n=6,7,11) are slightly less stable. The present study would be useful for the experimental/spectroscopic investigation of not only the dissociation phenomena of HF but also the similarity of the HF-water clusters to the water clusters.

Alternative Low-Energy Mechanisms for Isotopic Exchange in Gas-Phase D2O-H+(H2O)nReactions

Molecular-dynamics (MD) trajectories and high-level ab initio methods have been used to study the low-energy mechanism for D(2)O-H(+)(H(2)O)(n) reactions. At low collisional energies, MD simulations show that the collisional complexes are long-lived and undergo fast monomolecular isomerization, converting between different isomers within 50-500 ps. Such processes, primarily involving water-molecule shifts along a water chain, require the surmounting of very-low-energy barriers and present sizable non- Rice-Ramsperger-Kassel-Marcus (RRKM) effects, which are interpreted as a lack of randomization of the internal kinetic energy. Interestingly, the rate of water shifts was found to increase upon increasing the size of the cluster. Based on these findings, we propose to incorporate the following steps into the mechanism for low-energy isotopic scrambling these D(2)O-H(+)(H(2)O)(n) reactions: a) formation of the collisional complex [H(+)(H(2)O)(n)D(2)O]* in a vibro-rotational excited state; b) incorporation of the heavy-water molecule in the cluster core as HD(2)O(+) by means of isomerization involving molecular shifts; c) displacement of solvation molecules from the first shell of HD(2)O(+) inducing de-deuteration (shift of a D(+) to a neighbor water molecule); d) reorganization of the clusters and/or expulsion of one of the isotopic variants of water (H(2)O, HDO or D(2)O) from the periphery of the complex.

The properties of ion-water clusters. I. The protonated 21-water cluster

The ab initio atom-centered density-matrix propagation approach and the multistate empirical valence bond method have been employed to study the structure, dynamics, and rovibrational spectrum of a hydrated proton in the "magic" 21 water cluster. In addition to the conclusion that the hydrated proton tends to reside on the surface of the cluster, with the lone pair on the protonated oxygen pointing "outwards," it is also found that dynamical effects play an important role in determining the vibrational properties of such clusters. This result is used to analyze and complement recent experimental and theoretical studies.

Spectral Signatures of Hydrated Proton Vibrations in Water Clusters

The ease with which the pH of water is measured obscures the fact that there is presently no clear molecular description for the hydrated proton. The mid-infrared spectrum of bulk aqueous acid, for example, is too diffuse to establish the roles of the putative Eigen (H3O+) and Zundel (H5O2+) ion cores. To expose the local environment of the excess charge, we report how the vibrational spectrum of protonated water clusters evolves in the size range from 2 to 11 water molecules. Signature bands indicating embedded Eigen or Zundel limiting forms are observed in all of the spectra with the exception of the three- and five-membered clusters. These unique species display bands appearing at intermediate energies, reflecting asymmetric solvation of the core ion. Taken together, the data reveal the pronounced spectral impact of subtle changes in the hydration environment.

Structure of protonated water clusters: Low-energy structures and finite temperature behavior

The structure of protonated water clusters H+(H2O)n (n=5-22) are examined by two Monte Carlo methods in conjunction with the OSS2 potential [L. Ojamae, I. Shavitt, and S. J. Singer J. Chem. Phys. 109, 5547 (1998)]. The basin-hopping method is employed to explore the OSS2 potential energy surface and to locate low-energy structures. The topology of the "global minimum," the most stable low-energy structure, changes from single ring to multiple ring to polyhedral cage as the cluster size grows. The temperature dependence of the cluster geometry is examined by carrying out parallel tempering Monte Carlo simulations. Over the temperature range we studied (25-330 K), all water clusters undergo significant structural changes. The trends are treelike structures dominating at high temperature and single-ring structures appearing in slightly lower temperatures. For n> or =7, an additional transition from single ring to multiple rings appears as the temperature decreases. Only for n> or =16 do polyhedral structures dominate the lowest temperature range. Our results indicate very dynamic structural changes at temperature range relevant to atmospheric chemistry and current experiments. The structures and properties of medium-sized protonated clusters in this temperature range are far from their global minimum cousins. The relevance of these findings to recent experiments and theoretical simulations is also discussed.

Dissociation chemistry of hydrogen halides in water

To understand the mechanism of aqueous acid dissociation, which plays a fundamental role in aqueous chemistry, the ionic dissociation of HX acids (X=F, Cl, Br, and I) in water clusters up to hexamer is examined using density functional theory and Møller-Plesset second-order perturbation methods (MP2). Further accurate analysis based on the coupled clusters theory with singles and doubles excitations agrees with the MP2 results. The equilibrium structures, binding energies, electronic properties, stretching frequencies, and rotational constants of HX(H(2)O)(n) and X(-)(H(3)O)(+)(H(2)O)(n-1) are calculated. The dissociated structures of HF and HCl can be formed for n>/=4, while those of HBr and HI can be formed for n>/=3. Among these, the dissociated structures of HX (X=Cl, Br, and I) are more stable than the undissociated ones for n>/=4, while such cases for HF would require much more than six water molecules, in agreement with previous reports. The IR spectra of stable clusters including anharmonic frequencies are predicted to facilitate IR experimental studies. Undissociated systems have X-H stretching modes which are highly redshifted by hydration. Dissociated hydrogen halides show three characteristic OH stretching modes of hydronium moiety, which are redshifted from the OH stretching modes of water molecules.