Energetics of proton transfer in liquid water. I. Ab initio study for origin of many-body interaction and potential energy surfaces

The energetics of proton transfer in liquid water investigated by using ab initio calculation. The molecular electronic interaction of hydrated proton clusters in classified into many-body interaction elements by a new energy decomposition method. It is found that up to three-body molecular interaction is essential to describe the potential energy surface. The three-body effect mainly arises from the (non-classical) charge transfer and strongly depends on their configuration. Higher than three-body effects are small enough to be neglected. To simulate the liquid state reactions, two cluster models including all water molecules up to the second shell in the proton transfer reactions are employed. It is shown that these proton transfer reactions only involve small potential energy barriers of a few kcal/mol or less when structural rearrangement of the solvent is induced along the proton movement.     

Hydrated excess proton at water-hydrophobic interfaces

The behavior of the hydrated excess proton is investigated at the water-vapor, water-hydrophobic wall, and water-carbon tetrachloride interfaces through molecular dynamics simulations using the third-generation multistate empirical valence bond model (MS-EVB3). The MS-EVB3 simulations show a surface preference of the excess proton at the water-vapor interface, consistent with the discovery of this effect using an earlier version of the MS-EVB model (Petersen et al. J. Phys. Chem. B 2004, 108, 14804) and with the experimental results. The preference of the hydrated excess proton for other water-hydrophobic interfaces is also analyzed for the first time. The hydrated proton structures and charge defect delocalization effects at these interfaces are discussed in detail. By decomposing the free energy profiles into the internal energy and entropic contributions, the thermodynamic (free energy) driving forces for the surface preference of the excess proton are also elaborated. These results indicate that the "rigid" hydrated proton structures at all the interfaces are energetically (as opposed to entropically) stabilized due to the "amphiphilic" nature of the hydrated excess proton, resulting in its overall interfacial concentration enhancement. The effects of acid concentration and nuclear quantization are also explored.     

Orientational dynamics of water in the nafion polymer electrolyte membrane and its relationship to proton transport

Orientational anisotropies are calculated from molecular dynamics simulations of bulk water and the Na(+) and H(+) forms of hydrated Nafion and then compared with corresponding experimental values. The extended jump model of Laage and Hynes is applied to water reorientations for each system, and the anisotropies are explored as a product of hydrogen bond restricted "wobble-in-a-cone" reorientations and that due to the discrete jumps of hydrogen bond reorganization. Additionally, the timescales of hydrogen bond switching and proton transport are presented for bulk water and the H(+) form of hydrated Nafion. The short time scale of proton hopping is found to be independent of Nafion water loading, suggesting the short time dynamics of proton hopping are relatively insensitive to the level of hydration. Furthermore, the long time decay for the forward rate of hydrogen bond switching is shown to be identical to the long time decay in the forward rate of proton hopping, for bulk water and all water loadings of Nafion investigated, suggesting a unified process.     

Ab initio integrated multi-center molecular orbitals method for large cluster systems: total energy and normal vibration

A new computational scheme integrating multi-center ab initio molecular orbitals for determining total energy and normal vibration of large cluster systems is presented. This method can be used to treat large cluster systems such as solvents by quantum mechanics. The geometry parameters, the total energies, the relative energies, and the normal vibrations for four models of water cluster, the hydrated hydronium ion complex, and the transition state of proton transfer are calculated by the present method and are compared with those obtained by the full ab initio MO method. The results agree very well. The scheme proposed in this article is also intended to be used in modeling computer cluster systems using parallel algorithms.     

Unusual "amphiphilic" association of hydrated protons in strong acid solution

The solvation structure of the hydrated excess proton in concentrated aqueous HCl solution is studied using the self-consistent iterative multi-state empirical valence bond method. At 0.43-0.85 M concentrations, hydronium cations are found to form unusual cation pairs. This behavior is consistent with our earlier finding that hydronium cations can have an "amphiphilic" character due in part to the asymmetric nature of their hydrogen bonding to nearby water molecules. The existence of these hydronium amphiphilic pairs is further supported by a Car-Parrinello ab initio molecular dynamics simulation at 1.0 M HCl concentration. It is also found that the hydronium cation pairs are stabilized by a delocalization of the hydrated excess proton charge defects involving additional water molecules. At the higher concentrations of 1.68 and 3.26 M, the abundance of such hydronium pairs decreases, and the analysis of the radial distribution functions indicates the possible formation of an aggregate structure with longer-ranged order.     

Microhydration effects on the electronic spectra of protonated polycyclic aromatic hydrocarbons: [naphthalene-(H2O)(n = 1,2)]H+

Vibrational and electronic spectra of protonated naphthalene (NaphH(+)) microsolvated by one and two water molecules were obtained by photofragmentation spectroscopy. The IR spectrum of the monohydrated species is consistent with a structure with the proton located on the aromatic molecule, NaphH(+)-H(2)O. Similar to isolated NaphH(+), the first electronic transition of NaphH(+)-H(2)O (S(1)) occurs in the visible range near 500 nm. The doubly hydrated species lacks any absorption in the visible range (420-600 nm) but absorbs in the UV range, similar to neutral Naph. This observation is consistent with a structure, in which the proton is located on the water moiety, Naph-(H(2)O)(2)H(+). Ab initio calculations for [Naph-(H(2)O)(n)]H(+) confirm that the excess proton transfers from Naph to the solvent cluster upon attachment of the second water molecule.     

sodium Hypochlorite Reaction. II Degradation of N,N-Dimethylhydrazides

The Hofmann degradation² involves treatment of carboxamides with an alkaline hypohalite solution, thereby engendering primary amines with one less carbon. It is well known that formation of N-haloamides constitutes the first step of the degradation. Removal of the remaining proton on the nitrogen atom triggers rearrangement with concomitant expulsion of the halide ion to afford isocyanates which are immediately hydrated and then decarboxylate in situ. Consequently, only unsubstituted amides are capable of undergoing the Hofmann reaction.     

On the nature and signatures of the solvated electron in water

The hydrated electron is one of the simplest chemical transients and has been the subject of intense investigation and speculation since its discovery in 1962 by Hart and Boag. Although extensive kinetic and spectroscopic research on this species has been reported for many decades, its structure, i.e., the dominant electron-water binding motif, and its binding energy remained uncertain. A recent milestone in the research on the hydrated electron was the determination of its binding energy by liquid-jet photoelectron spectroscopy. It turned out that the assumption of a single electron binding motif in liquid water is an oversimplification. In addition to different isomers in cluster spectroscopy and different transient species of unknown structure in ultrafast experiments, long-lived hydrated electrons near the surface of liquid water have recently been discovered. The present article gives an account of recent work on the topic "solvated electrons" from the perspectives of cluster spectroscopy, condensed-phase spectroscopy, as well as theory. It highlights and critically discusses recent findings and their implications for our understanding of electron solvation in aqueous environments.     

Quantum simulation of a hydrated noradrenaline analog with the torsional path integral method

An extended version of the torsional path integral Monte Carlo (TPIMC) method is presented and shown to be useful for studying the conformation of flexible molecules in solvated clusters. The new technique is applied to the hydrated clusters of the 2-amino-1-phenyl-ethanol (APE) molecule. APE + nH2O clusters with n = 0-4 are studied at 100 and 300 K using both classical and quantum simulations. Only at the lower temperature is the hydration number n found to impact the conformational distribution of the APE molecule. This is shown to be a result of the temperature-dependent balance between the internal energy and entropy contributions to the relative conformer free energies. Furthermore, at 100 K, large quantum effects are observed in the calculated conformer populations. A particularly large quantum shift of 30% of the total population is calculated for the APE + 2H2O cluster, which is explained in terms of the relative zero point energy of the lowest-energy hydrated structures for this cluster. Finally, qualitative agreement is found between the reported calculations and recent spectroscopy experiments on the hydrated clusters of APE, including an entropically driven preference for the formation of AG-type hydrated structures and the formation of a water "droplet" in the APE + 4H2O cluster.     

Computer Simulation of Water Clusterization on Chlorine Ions: 2. Microstructure

The molecular structure of Cl^–(H₂O) _n clusters, n = 1–60, in equilibrium with vapor, and the cluster with n = 500 was studied by the Monte Carlo method. The first hydrated layer of a cluster is formed in unsaturated water vapors. The second hydrated layer begins to be formed in saturated vapor. The position of hydrated layers is not changed with an increase in cluster size and coincides with the position of the hydrated layers of ions in aqueous solutions of weak electrolytes. Orientational order in a cluster also has the layered structure. The orientation of molecules between the layers is random. The stability of the first layer is ensured only due to direct interactions with ions, whereas the stability of subsequent layers is due to cooperative interactions between molecules and between molecules and ions. As temperature decreases, the effect of ion displacement to the cluster surface becomes stronger.     

Oxygen reduction on a platinum cluster

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

Proton storage site in bacteriorhodopsin: new insights from quantum mechanics/molecular mechanics simulations of microscopic pK(a) and infrared spectra

Identifying the group that acts as the proton storage/loading site is a challenging but important problem for understanding the mechanism of proton pumping in biomolecular proton pumps, such as bacteriorhodopsin (bR) and cytochrome c oxidase. Recent experimental studies of bR propelled the idea that the proton storage/release group (PRG) in bR is not an amino acid but a water cluster embedded in the protein. We argue that this idea is at odds with our knowledge of protein electrostatics, since invoking the water cluster as the PRG would require the protein to raise the pK(a) of a hydronium by almost 11 pK(a) units, which is difficult considering known cases of pK(a) shifts in proteins. Our recent quantum mechanics/molecular mechanics (QM/MM) simulations suggested an alternative "intermolecular proton bond" model in which the stored proton is shared between two conserved Glu residues (194 and 204). Here we show that this model leads to microscopic pK(a) values consistent with available experimental data and the functional requirement of a PRG. Extensive QM/MM simulations also show that, independent of a number of technical issues, such as the influence of QM region size, starting X-ray structure, and nuclear quantum effects, the "intermolecular proton bond" model is qualitatively consistent with available spectroscopic data. Potential of mean force calculations show explicitly that the stored proton strongly prefers the pair of Glu residues over the water cluster. The results and analyses help highlight the importance of considering protein electrostatics and provide arguments for why the "intermolecular proton bond" model is likely applicable to the PRG in biomolecular proton pumps in general.     

Proton relay in a one-dimensional hydrogen-bonded chain composed of water molecules and a squaric acid derivative.

The crystal structure of a hydrated crystal of bis(squaryl)biphenyl (BSQB*4H2O), in which two squaric acid moieties are connected with a 4,4'-biphenyl unit, was characterized by the presence of a one-dimensional hydrogen-bonded chain composed of BSQB and water molecules. X-ray crystallographic analysis showed that BSQB exists in a dianion form and that, on average, two of the four water molecules are protonated. The enhanced temperature dependence of the thermal parameters of the oxygen atoms of the water molecules suggested dynamic disorder of the water molecules. The solid-state magic angle spinning deuterium NMR spectrum of BSQB*4D2O revealed that deuterons are exchanged between heavy water molecules and oxonium ions with an exchange rate of ca. 700 Hz around 250 K and that deuterons start to migrate in a hydrogen-bonded cluster of water molecules. Ac dielectric measurements were also used to examine the dynamic process in the hydrated crystal. The dielectric permittivity of the crystal dramatically increased above 250 K with a distinct frequency dependence (epsilon' = 4.7 x 10(4) at 340 K and 1 kHz). The frequency dependence of tan delta at 290 K exhibited a maximum at 3.0 kHz, and this maximum shifted to lower frequencies when the temperature of the crystal decreased. These experimental results suggested that in the one-dimensional hydrogen-bonded chain of BSQB*4H2O a proton relay between oxonium ions and water molecules occurred within a cluster of four water molecules and that the relay was transmitted to the adjacent cluster mediated by the modulation of the negative charge distribution of the BSQB dianion. These phenomena were interpreted as the solitonic migration of the charged domain boundaries along the one-dimensional hydrogen-bonded chain.     

Gold Nanostructures Formed in Ionic Clusters of Perfluorinated Ionomer

Summary: Gold nanoparticles and nanorods have been synthesized in the ionic cluster network of Nafion using a liquid‐phase chemical impregnation/reduction process seemingly visualizing the phase‐separated hydrated cluster network of Nafion. The spherical shape of the hydrated clusters is clearly seen by the replicated gold nanoparticles with diameters of 5–6 nm. The gold nanoparticles grow through the Nafion cluster network to form gold nanorods with diameters of 5–6 nm and a characteristic length of 12–15 nm, which compare well to the previously reported small‐angle X‐ray scattering results. Tiny channels between adjacent spherical clusters are also observed, which expand to form a rod structure by a reorganization of ion exchange sites and an elastic deformation of Nafion polymer chains. The gold nanorods become interconnected in a cascadic feature, and a tripod‐shaped nanorod structure is one of the most commonly observed structures of the replicated gold. Although further study should be performed, the synthesized gold may be used to visualize the hydrated cluster and the network structure of Nafion, which could be used as important information in identifying the morphology and ion transport phenomena of Nafion.

Gold nanoparticles and nanorods growing through a hydrated Nafion cluster network.     

The structure of the hydrated electron. Part 2. A mixed quantum/classical molecular dynamics embedded cluster density functional theory: single-excitation configuration interaction study

Adiabatic mixed quantum/classical (MQC) molecular dynamics (MD) simulations were used to generate snapshots of the hydrated electron in liquid water at 300 K. Water cluster anions that include two complete solvation shells centered on the hydrated electron were extracted from the MQC MD simulations and embedded in a roughly 18 Ax18 Ax18 A matrix of fractional point charges designed to represent the rest of the solvent. Density functional theory (DFT) with the Becke-Lee-Yang-Parr functional and single-excitation configuration interaction (CIS) methods were then applied to these embedded clusters. The salient feature of these hybrid DFT(CIS)/MQC MD calculations is significant transfer (approximately 18%) of the excess electron's charge density into the 2p orbitals of oxygen atoms in OH groups forming the solvation cavity. We used the results of these calculations to examine the structure of the singly occupied and the lower unoccupied molecular orbitals, the density of states, the absorption spectra in the visible and ultraviolet, the hyperfine coupling (hfcc) tensors, and the infrared (IR) and Raman spectra of these embedded water cluster anions. The calculated hfcc tensors were used to compute electron paramagnetic resonance (EPR) and electron spin echo envelope modulation (ESEEM) spectra for the hydrated electron that compared favorably to the experimental spectra of trapped electrons in alkaline ice. The calculated vibrational spectra of the hydrated electron are consistent with the red-shifted bending and stretching frequencies observed in resonance Raman experiments. In addition to reproducing the visible/near IR absorption spectrum, the hybrid DFT model also accounts for the hydrated electron's 190-nm absorption band in the ultraviolet. Thus, our study suggests that to explain several important experimentally observed properties of the hydrated electron, many-electron effects must be accounted for: one-electron models that do not allow for mixing of the excess electron density with the frontier orbitals of the first-shell solvent molecules cannot explain the observed magnetic, vibrational, and electronic properties of this species. Despite the need for multielectron effects to explain these important properties, the ensemble-averaged radial wavefunctions and energetics of the highest occupied and three lowest unoccupied orbitals of the hydrated electrons in our hybrid model are close to the s- and p-like states obtained in one-electron models. Thus, one-electron models can provide a remarkably good approximation to the multielectron picture of the hydrated electron for many applications; indeed, the two approaches appear to be complementary.     

Microhydration shell structure in Cl2*-.nH2O clusters: A theoretical study

We present the results of a detailed study on structure and electronic properties of hydrated cluster Cl2*-.nH2O (n = 1-7) based on a nonlocal density functional, namely, Becke's [J. Chem. Phys. 98, 1372 (1993)] half and half hybrid exchange-correlation functional with a split valence 6-311++G(d,p) basis function. Geometry optimizations for all the clusters are carried out with various possible initial guess structures without any symmetry restriction. Several minimum energy structures (conformers) are predicted with a small difference in total energy. There is a competition between the binding of solvent H2O units with Cl2*- dimer radical anion directly through ion-molecule interaction and forming interwater hydrogen-bonding network in Cl2*-.nH2O (n > or = 2) hydrated cluster. Structure having interwater H-bonded network is more stable over the structure where H2O units are connected to the solute dimer radical anion Cl2*- rather independently either by single or double H bonding in a particular size (n) of hydrated cluster Cl2*-.nH2O. At the maximum four solvent H2O units reside in interwater H-bonding network present in these hydrated clusters. It is observed that up to six H2O units are independently linked to the anion having four double H bondings and two single H bondings suggesting the primary hydration number of Cl2*- to be 6. In all these clusters, the odd electron is found to be mostly localized over the two Cl atoms and these two atoms are bound by a three-electron hemibond. Calculated interaction (between solute and different water clusters) and vertical detachment energy profiles show saturation at n = 6 in the hydrated cluster Cl2*-.nH2O (n = 1-7). However, calculated solvation energy increases with the increase in number of solvent H2O molecules in the cluster. Interaction energy varies linearly with vertical detachment energy for the hydrated clusters Cl2*-.nH2O (n < or = 6). Calculation of the vibration frequencies show that the formation of Cl2*(-)-water clusters induces significant shifts from the normal stretching modes of isolated water. A clear difference in the pattern of IR spectra is observed in the O-H stretching region of water from hexa- to heptahydrated cluster.     

Influence of chemical composition and sequence length on the transport properties of proton exchange membranes

One of the integral parts of the fuel cell is the proton exchange membrane. Our research group has been engaged in the past few years in the synthesis of several sulfonated poly(arylene ether) random copolymers. The copolymers were varied in both the bisphenol structure as well as in the functional groups in the backbone such as sulfone and ketones. To compare the effect of sequence length, multiblock copolymers based on poly(arylene ether sulfone)s were synthesized. This paper aims to describe our investigation of the effect of chemical composition, morphology, and ion exchange capacity (IEC) on the transport properties of proton conducting membranes. The key properties examined were proton conductivity, methanol permeability, and water self diffusion coefficient in the membranes. It was observed that under fully hydrated conditions, proton conductivity for both random and block copolymers was a function of IEC and water uptake. However, under partially hydrated conditions, the block copolymers showed improved proton conductivity over the random copolymers. The proton conductivity for the block copolymer series was found to increase with increasing block lengths under partially hydrated conditions. © 2006 Wiley Periodicals, Inc. J Polym Sci Part B: Polym Phys 44: 2226–2239, 2006     

Density functional study of intradimer proton transfers in hydrated adenine dimer ions, A2(+)(H2O)n (n = 0-2)

Structures of mono- and dihydrated adenine dimers and their cations were calculated using B3LYP density functional theory with the 6-31+G(d,p) basis set, in order to help understand photofragmentation experiments of hydrated adenine dimers from the energetics point of view. Several important pathways leading to the major fragmentation product, protonated adenine ion (AH(+)), thermodynamically at minimum costs were investigated at the ground-state electronic potential surface of hydrated adenine dimer cations. Our calculations suggest that the proton transfer from one adenine moiety to the other in hydrated dimer ions readily occurs with negligible barriers in normal hydration conditions. In asymmetrically hydrated ions, however, the proton transfer to more hydrated adenine moieties is kinetically hindered due to heightened transition-state barriers, while the other way is still barrierless. Such directional preference in proton transfer may be characterized as a unique dimer ion property, stemming from the difference in basicity of the two nitrogen atoms involved in the double hydrogen bond that would be equivalent without hydration. We also found that dimer cleavage requires about 4 times larger energy than evaporation of individual water molecules, so it is likely that most solvent molecules evaporate before the eventual dimer cleavage when available internal energy is limited.