Excess proton solvation and delocalization in a hydrophilic pocket of the proton conducting polymer membrane nafion

Solvation properties of the hydrated excess proton are studied in a hydrophilic pocket of Nafion 117 through a series of molecular dynamics simulations. The multistate empirical valence bond (MS-EVB) methodology, which enables the delocalization of the excess proton through the Grotthuss hopping mechanism, was employed for one of the excess protons in the simulation cell. Simulations were performed such that "classical" nondissociable hydronium cations and a single excess proton treated with the MS-EVB methodology were at a concentration ratio of 39:1. Two degrees of hydration of the Nafion polymer electrolyte membrane were simulated, each displaying the same marked difference between the solvation structures of the classical versus MS-EVB treated (Grotthuss shuttling) excess proton species. These differences are attributed to the solvent dynamics needed to transfer the cation between the solvent separated and contact pair positions about the sulfonic acid counterion. The results demonstrate that it is generally impossible to describe the low pH conditions in the hydrophilic domains of Nafion without the explicit treatment of Grotthuss delocalization in the underlying molecular dynamics model for the excess protons.     

A molecular dynamics study of a nafion polyelectrolyte membrane and the aqueous phase structure for proton transport

A molecular dynamics simulation study of hydrated Nafion at water contents ranging from 5 to 20 wt % was performed to examine the structure and dynamics of the hydrated polyelectrolyte system. The simulations show that the system forms segregated hydrophobic regions consisting primarily of the polymer backbone and hydrophilic regions with an inhomogeneous water distribution. We find that the water clustering strongly depends on the water content. At low water content, only isolated small water clusters are formed. As the water content increases, it becomes increasingly possible that a predominant majority of water molecules form a single cluster, suggesting that the hydrophilic regions become connected. We characterize the atomic structures formed within the system by various atomic pair correlation functions. The water structure factor shows a peak at q values corresponding to an intercluster distance about 2.5 nm and greater. With increasing water content, the distance moves to larger values, consistent with findings from scattering experiments. We find that the degree of solvation of hydronium ions by water molecules is a strong function of water content. At 5 wt %, a majority of the hydronium ions are hydrated by no more than two water molecules, prohibiting structural diffusion. As water content increases, the hydronium ions continue to become increasingly hydrated, resulting in structures capable of forming eigen ions, a necessary step in structural diffusion. Addressing the experimentally observed fact that conductivity in these membranes abruptly drops near 5 wt %, we find that both the local structure of the poorly hydrated hydronium ions and the disconnected nature of the global morphology of the water nanonetwork at low water content should contribute to poor conductivity.     

From a localized H3O radical to a delocalized H3O+···e- solvent-separated pair by sequential hydration

The impact of microhydration on the electronic structure and reactivity of the H(3)O moiety is investigated by ab initio calculations. In the gas phase, H(3)O is a radical with spin density localized on its hydrogen end, which is only kinetically stable and readily decomposes into a water molecule and a hydrogen atom. When solvated by a single water molecule, H(3)O preserves to a large extent its radical character, however, two water molecules are already capable to shift most of the spin density to the solvent. With three solvating water molecules this shift is practically completed and the system is best described as a solvent-separated pair of a hydronium cation and a hydrated electron. The electronic structure of this system and its proton transfer reactivity leading to formation of a hydrogen atom already resemble those of a proton-electron pair in bulk water.     

Hydration of Y3+ ion: a Car-Parrinello molecular dynamics study

The solvation shell structure of Y3+ and the dynamics of the hydrated ion in an aqueous solution of 0.8 M YCl3 are studied in two conditions with and without an excess proton by using first principles molecular dynamics method. We find that the first solvation shell around Y3+ contains eight water molecules forming a square antiprism as expected from x-ray absorption near edge structure in both the conditions we examined. A detailed analysis relying upon localized orbitals reveals that the complexation of water molecules with yttrium cation leads to a substantial amount of charge redistribution particularly on the oxygen atoms, giving rise to the chemical shifts of approximately -20 ppm in 17O nuclear magnetic resonance relative to the computed nuclear shieldings of the bulk water.     

Adsorption, ionization, and migration of hydrogen chloride on ice films at temperatures between 100 and 140 K

Adsorption of hydrogen chloride (HCl) on water ice films is studied in the temperature range of 100-140 K by using Cs+ reactive ion scattering (Cs+ RIS), low energy sputtering (LES), and temperature-programmed-desorption mass spectrometry (TPDMS). At 100 K, HCl on ice partially dissociates to hydronium and chloride ions and the undissociated HCl exists in two distinct molecular states (alpha- and beta-states). Upon heating of the ice films, HCl molecules in the alpha-state desorb at 135-150 K, whereas those in the beta-state first become ionized and then desorb via recombinative reaction of ions at 170 K. An adsorption kinetics study reveals that HCl adsorption into the ionized state is slightly favored over adsorption into the molecular states at 100 K, leading to earlier saturation of the ionized state. Between the two molecular states, the beta-state is formed first, and the alpha-state appears only at high HCl coverage. At 140 K, ionic dissociation of HCl is completed. The resulting hydronium ion can migrate into the underlying sublayer to a depth <4 bilayers, suggesting that the migration is assisted by self-diffusion of water molecules near the surface. When HCl is covered by a water overlayer at 100 K, its ionization efficiency is enhanced, but a substantial portion of HCl remains undissociated as molecules or contact ion pairs. The observation suggests that three-dimensional surrounding by water molecules does not guarantee ionic dissociation of HCl. Complete ionization of HCl requires additional thermal energy to separate the hydronium and chloride ions.     

Molecular modeling of proton transport in the short-side-chain perfluorosulfonic acid ionomer

An explanation for the superior proton conductivity of low equivalent weight (EW) short-side-chain (SSC) perfluorosulfonic acid membranes is pursued through the determination of hydrated morphology and hydronium ion diffusion coefficients using classical molecular dynamics (MD) simulations. A unique force field set for the SSC ionomer was derived from torsion profiles determined from ab initio electronic structure calculations of an oligomeric fragment consisting of two side chains. MD simulations were performed on a system consisting of a single macromolecule of the polymer (EW of 580) with the general formula F3C-[CF(OCF2CF2SO3H)-(CF2)7]40-CF3 at hydration levels corresponding to 3, 6, and 13 water molecules per sulfonic acid group. Examination of the hydrated morphology indicates the formation of hydrogen bond "bridges" between distant sulfonate groups without significant bending of the polytetrafluoroethylene backbone. Pair correlation functions of the system identify the presence of ion cages consisting of hydronium ions hydrogen-bonded to three sulfonate groups at the lowest water content. Such structures exhibit very low S-OH3+ separations, well below 4 A and severely inhibit vehicular diffusion of the protons. The number of sulfonate groups in the first solvation shell of a given hydronium ion correlates well with the differences between Nafion and the SSC polymer (Hyflon). The calculated hydronium ion diffusion coefficients of 2.84 x 10-7, 1.36 x 10-6, and 3.47 x 10-6 cm2/s for water contents of 3, 6, and 13, respectively, show only good agreement to experimentally measured values at the lowest water content, underscoring the increasing contribution of proton shuttling or hopping at the higher hydration levels. At the highest water content, the vehicular diffusion accounts for only about 1/5 of the total proton transport similar to that observed in Nafion.     

Special pair dance and partner selection: elementary steps in proton transport in liquid water

Conditional and time-dependent radial distribution functions reveal the details of the water structure surrounding the hydronium during the proton mobility process. Using this methodology for classical multistate empirical valence bond (MS-EVB) and ab initio molecular dynamics trajectories, as well as quantal MS-EVB trajectories, we supply statistical proof that proton hops in liquid water occur by a transition from the H3O+[3H2O] Eigen-complex, via the H5O2+ Zundel-complex, to a H3O+[3H2O] centered on a neighboring water molecule. In the "resting period" before a transition, there is a distorted hydronium with one of its water ligands at a shorter distance and another at a longer distance than average. The identity of this "special partner" interchanges rapidly within the three first-shell water ligands. This is coupled to cleavage of an acceptor-type hydrogen bond. Just before the transition, a partner is selected by an additional translation of the H3O+ moiety in its direction, possibly enabled by loosening of donor-type hydrogen bonds on the opposite side. We monitor the transition in real time, showing how the average structure is converted to a distorted H5O2+ cation constituting the transitional complex for proton hopping between water molecules.     

Protons in supercritical water: a multistate empirical valence bond study

Molecular dynamics simulations have been performed to analyze microscopic details related to aqueous solvation of excess protons along the supercritical T = 673 K isotherm, spanning a density interval from a typical liquid down to vapor environments. The simulation methodology relies on a multistate empirical valence bond Hamiltonian model that includes a proton translocation mechanism. Our results predict a gradual stabilization of the solvated Eigen cation [H(3)O.(H(2)O)(3)](+) at lower densities, in detriment of the symmetric Zundel dimer [H.(H(2)O)(2)](+). At all densities, the average solvation structure in the close vicinity of the hydronium is characterized by three hydrogen bond acceptor water molecules and presents minor changes in the solute water distances. Characteristic times for the proton translocation jumps have been computed using population relaxation time correlation functions. Compared to room temperature results, the rates at high densities are 4 times faster and become progressively slower in steamlike environments. Diffusion coefficients for the excess proton have also been computed. In agreement with conductometric data, our results show that contributions from the Grotthus mechanism to the overall proton transport diminish at lower densities and predict that in steamlike environments, the proton diffusion is almost 1 order of magnitude slower than that for pure water. Spectroscopic information for the solvated proton is accordant to the gradual prevalence of proton localization in Eigen-like structures at lower densities.     

Effects of electric fields on proton transport through water chains

Molecular dynamics simulations on quantum energy surfaces are carried out to study the effects of perturbing electric fields on proton transport (PT) in protonated water chains. As an idealized model of a hydrophobic cavity in the interior of a protein the water molecules are confined into a carbon nanotube (CNT). The water chain connects a hydrated hydronium ion (H3O+) at one end of the CNT and an imidazole molecule at the other end. Without perturbing electric fields PT from the hydronium proton donor to the imidazole acceptor occurs on a picosecond time scale. External perturbations to PT are created by electric fields of varying intensities, normal to the CNT axis, generated by a neutral pair of charges on the nanotube wall. For fields above approximately 0.5 VA, the hydronium ion is effectively trapped at the CNT center, and PT blocked. Fields of comparable strength are generated inside proteins by nearby polar/charged amino acids. At lower fields the system displays a rich dynamic behavior, where the excess charge shuttles back and forth along the water chain before reaching the acceptor group on the picosecond time scale. The effects of the perturbing field on the proton movement are analyzed in terms of structural and dynamic properties of the water chain. The implications of these observations on PT in biomolecular systems and its control by external perturbing fields are discussed.     

Interfacial structures of acidic and basic aqueous solutions

Phase-sensitive sum-frequency vibrational spectroscopy was used to study water/vapor interfaces of HCl, HI, and NaOH solutions. The measured imaginary part of the surface spectral responses provided direct characterization of OH stretch vibrations and information about net polar orientations of water species contributing to different regions of the spectrum. We found clear evidence that hydronium ions prefer to emerge at interfaces. Their OH stretches contribute to the "ice-like" band in the spectrum. Their charges create a positive surface field that tends to reorient water molecules more loosely bonded to the topmost water layer with oxygen toward the interface, and thus enhances significantly the "liquid-like" band in the spectrum. Iodine ions in solution also like to appear at the interface and alter the positive surface field by forming a narrow double-charge layer with hydronium ions. In NaOH solution, the observed weak change of the "liquid-like" band and disappearance of the "ice-like" band in the spectrum indicates that OH(-) ions must also have excess at the interface. How they are incorporated in the interfacial water structure is, however, not clear.     

Influence of Hydronium Ions in Zeolites on Sorption

In the presence of sufficient concentrations of water, stable, hydrated hydronium ions are formed in the pores and at the surface of solid acids such as zeolites. For a medium‐pore zeolite, such as zeolite MFI, hydrated hydronium ions consist of eight water molecules and have an effective volume of 0.24 nm³. In their presence, larger organic molecules can only adsorb in the portions of the pore that are not occupied by hydronium ions. As a consequence, the available pore volume decreases proportionally to the concentration of the hydronium ions. The higher charge density (the increasing ionic strength) that accompanies an increasing concentration of hydronium ions leads to an increase in the activity coefficients of the adsorbed substrates, thus, weakening the interactions between the organic part of the molecules and the zeolite and favoring the interactions with polar groups. The quantitative understanding of these interactions makes it possible to link a collective property such as hydrophilicity and hydrophobicity of zeolites to specific interactions on molecular level.     

A density functional theory for studying ionization processes in water clusters

A generalized Kohn-Sham (GKS) approach to density functional theory (DFT), based on the Baer-Neuhauser-Livshits range-separated hybrid, combined with ab initio motivated range-parameter tuning is used to study properties of water dimer and pentamer cations. The water dimer is first used as a benchmark system to check the approach. The present brand of DFT localizes the positive charge (hole), stabilizing the proton transferred geometry in agreement with recent coupled-cluster calculations. Relative energies of various conformers of the water dimer cation compare well with previously published coupled cluster results. The GKS orbital energies are good approximations to the experimental ionization potentials of the system. Low-lying excitation energies calculated from time-dependent DFT based on the present method compare well with recently published high-level "equation of motion-coupled-cluster" calculations. The harmonic frequencies of the water dimer cation are in good agreement with experimental and wave function calculations where available. The method is applied to study the water pentamer cation. Three conformers are identified: two are Eigen type and one is a Zundel type. The structure and harmonic vibrational structure are analyzed. The ionization dynamics of a pentamer water cluster at 0 K shows a fast <50 fs transient for transferring a proton from one of the water molecules, releasing a hydroxyl radical and creating a protonated tetramer carrying the excess hole.     

Thermal decomposition of hydronium jarosite (H3O)Fe3(SO4)2(OH)6

Thermogravimetry combined with mass spectrometry has been used to study the thermal decomposition of a synthetic hydronium jarosite. Five mass loss steps are observed at 262, 294, 385, 557 and 619°C. The mass loss step at 557°C is sharp and marks a sharp loss of sulphate as SO₃ from the hydronium jarosite. Mass spectrometry through evolved gases confirms the first three mass loss steps to dehydroxylation, the fourth to a mass loss of the hydrated proton and a sulphate and the final step to the loss of the remaining sulphate. Changes in the molecular structure of the hydronium jarosite were followed by infrared emission spectroscopy. This technique allows the infrared spectrum at the elevated temperatures to be obtained. Infrared emission spectroscopy confirms the dehydroxylation has taken place by 400 and the sulphate loss by 650°C. Jarosites are a group of minerals formed in evaporite deposits and form a component of the efflorescence. The minerals can function as cation and heavy metal collectors. Hydronium jarosite has the potential to act as a cation collector by the replacement of the proton with a heavy metal cation.     

Hydrogen detachment of the hydrated hydrohalogen acids upon attaching an excess electron

High level ab initio calculations are employed to investigate the excess electron attachment to the hydrated hydrohalogen acids. The excess electron leads to the dissociation of hydrogen halide acids, which results in the release of a hydrogen radical. Neutral HCl, HBr, and HI are dissociated by tetrahydration. Upon binding an excess electron, these hydrated hydrohalogen acids show that (i) the H-X bond strength weakens with redshifted H-X stretching frequencies, (ii) HX can have a bound-electron state, a dissociated structure, or a zwitter-ionic structure, and (iii) HClHBr is dissociated by tri/mono-hydration, while HI is dissociated even without hydration. This dissociation is in contrast to the case of electron attachment to hydrated hydrogen fluoric acids for which HF is not dissociated by more than ten water molecules.     

Generating Excess Protons in Microsolvated Acid Clusters under Ambient Conditions: An Issue of Configurational Entropy versus Internal Energy

Acid dissociation, and thus liberation of excess protons in small water droplets, impacts on diverse fields such as interstellar, atmospheric or environmental chemistry. At cryogenic temperatures below 1 K, it is now well established that as few as four water molecules suffice to dissociate the generic strong acid HCl, yet temperature-driven recombination sets in simply upon heating that cluster. Here, the fundamental question is posed of how many more water molecules are required to stabilize a hydrated excess proton at room temperature. Ab initio path integral simulations disclose that not five, but six water molecules are needed at 300 K to allow for HCl dissociation independently from nuclear quantum effects. In order to provide the molecular underpinnings of these observations, the classical and quantum free energy profiles were decomposed along the dissociation coordinate in terms of the corresponding internal energy and entropy profiles. What decides in the end about acid dissociation, and thus ion pair formation, in a specific microsolvated water cluster at room temperature is found to be a fierce competition between classical configurational entropy and internal energy, where the former stabilizes the undissociated state whereas the latter favors dissociation. It is expected that these are generic findings with broad implications on acid–base chemistry depending on temperature in small water assemblies.     

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.     

The electronic structure of the hydrated proton: a comparative X-ray absorption study of aqueous HCl and NaCl solutions

The oxygen K edge X-ray absorption spectra of aqueous HCl and NaCl solutions reveal distinct perturbations of the local water molecules by the respective solutes. While the addition of NaCl leads to large spectral changes, the effect of HCl on the observed X-ray absorption spectrum is surprisingly small. Density functional theory calculations suggest that this difference primarily reflects a strong blue shift of the hydrated proton (in either the Eigen (H9O4+) or Zundel (H2O5+) forms) spectrum relative to that of H2O, indicating the tighter binding of electrons in H3O+. This spectral shift counteracts the spectral changes that arise from direct electrostatic perturbation of water molecules in the first solvation shell of Cl-. Consequently, the observed spectral changes effected by HCl addition are minimal compared to those engendered by NaCl. Additionally, these results indicate that the effect of monovalent cations on the nature of the unoccupied orbitals of water molecules in the first solvation shell is negligible, in contrast to the large effects of monovalent anions.     

Structure of electrolyte solutions sorbed in carbon nanospaces, studied by the replica RISM theory

The replica RISM theory is used to investigate the structure of electrolyte solutions confined in carbonized polyvinylidene chloride (PVDC) nanoporous material, compared to bulk electrolyte solution. Comparisons are made between the models of electrolyte solution sorbed in the carbonized PVDC material and a single carbon nanosphere in bulk electrolyte solution. Particular attention is paid to the chemical potential balance between the species of the sorbed electrolyte solution and the bulk solution in contact with the nanoporous material. As a result of the strong hydrophobicity of the carbonized PVDC material in the absence of activating chemical groups, the densities of water and ions sorbed in the material are remarkably low compared to those in the ambient bulk solution. The interaction between water molecules and cations becomes strong in nanospaces. It turns out that, in carbon nanopores, a cation adsorbed at the carbon surface is fully surrounded by the hydration shell of water molecules which separates the cation and the surface. Distinctively, an anion is adsorbed in direct contact with the carbon surface, which squeezes a part of its hydration shell out. The tendency increases toward smaller cations, which are characterized as "positive hydration" ions. In the bulk, cations are not hydrated so strongly and behave similarly to anions. The results suggest that the specific capacitance of an electric double-layer supercapacitor with nanoporous electrodes is intimately related to the solvation structure of electrolyte solution sorbed in nanopores, which is affected by the microscopic structure of the nanoporous electrode.     

Hydrogen bonding mediated ion pairs of some aprotic ionic liquids and their structural transition in aqueous solution

Ion pair speciation of ionic liquids(ILs) has an important effect on the physical and chemical properties of ILs and recognition of the structure of ion pairs in solution is essential. It has been reported that ion pairs of some ILs can be formed by hydrogen bonding interactions between cations and anions of them. Considering the fact that far-IR(FIR) spectroscopy is a powerful tool in indicating the intermolecular and intramolecular hydrogen bonding, in this work, this spectroscopic technique has been combined with molecular dynamic(MD) simulation and nuclear magnetic resonance hydrogen spectroscopy(~1H NMR) to investigate ion pairs of aprotic ILs [Bmim][NO_3], [BuPy][NO_3], [Pyr_(14)][NO_3], [PP_(14)][NO_3] and [Bu-choline][NO_3] in aqueous IL mixtures. The FIR spectra have been assigned with the aid of density functional theory(DFT) calculations, and the results are used to understand the effect of cationic nature on the structure of ion pairs. It is found that contact ion pairs formed in the neat aprotic ILs by hydrogen bonding interactions between cation and anion, were still maintained in aqueous solutions up to high water mole fraction(say 0.80 for [BuPy][NO3]). When water content was increased to a critical mole fraction of water(say 0.83 for [BuPy][NO3]), the contact ion pairs could be transformed into solvent-separated ion pairs due to the formation of the hydrogen bonding between ions and water. With the further dilution of the aqueous ILs solution, the solvent-separated ion pairs was finally turned into free cations and free anions(fully hydrated cations or anions). The concentrations of the ILs at which the contact ion pairs were transformed into solvent-separated ion pairs and solvent-separated ion pairs were transformed into free ions(fully hydrated ion) were dependent on the cationic structures. These information provides direct spectral evidence for ion pair structures of the aprotic ILs in aqueous solution. MD simulation and ~1H NMR results support the conclusion drawn from FIR spectra investigations.