Ions in water: The microscopic structure of a concentrated HCl solution

Fast Structural Dynamics in Concentrated HCl Solutions: From Proton Hopping to the Bulk Viscosity

Concentrated acid solutions, particularly HCl, have been studied extensively to examine the proton hopping and infrared spectral signatures of hydronium ions. Much less attention has been given to the structural dynamics of concentrated HCl solutions. Here, we apply optical heterodyne detected-optical Kerr effect (OHD-OKE) measurements to examine HCl concentration-dependent dynamics from moderate (0.8 m) to very high (15.5 m) concentrations and compare the results to the dynamics of NaCl solutions, as Na+ is similar in size to the hydronium cation. Both HCl and NaCl OHD-OKE signals decay as triexponentials at all concentrations, in contrast to pure water, which decays as a biexponential. Two remarkable features of the HCl dynamics are the following: (1) the bulk viscosity is linearly related to the slowest decay constant, t3, and (2) the concentration-dependent proton hopping times, determined by ab initio MD simulations and 2D IR chemical exchange experiments, both obtained from the literature, fall on the same line as the slowest structural dynamics relaxation time, t3, within experimental error. The structural dynamics of hydronium/chloride/water clusters, with relaxation times t3, are responsible for the concentration dependence of microscopic property of proton hopping and the macroscopic bulk viscosity. The slowest time constant (t3), which does not have a counterpart in pure water, is 3 ps at 0.8 m and increases by a factor of ∼2 by 15.5 m. The two fastest HCl decay constants, t1 and t2, are similar to those of pure water and increase mildly with the concentration.

Anionic Effects on Concentrated Aqueous Lithium Ion Dynamics

The dynamics, orientational anisotropy, diffusivity, viscosity, and density were measured for concentrated lithium salt solutions, including lithium chloride (LiCl), lithium bromide (LiBr), lithium nitrite (LiNO2), and lithium nitrate (LiNO3), with methyl thiocyanate as an infrared vibrational probe molecule, using two-dimensional infrared spectroscopy (2D IR), nuclear magnetic resonance (NMR) spectroscopy, and viscometry. The 2D IR, NMR, and viscosity results show that LiNO2 exhibits longer correlation times, lower diffusivity, and nearly 4 times greater viscosity compared to those of the other lithium salt solutions of the same concentration, suggesting that nitrite anions may strongly facilitate structure formation via strengthening water-ion network interactions, directly impacting bulk solution properties at sufficiently high concentrations. Additionally, the LiNO2 and LiNO3 solutions show significantly weakened chemical interactions between the lithium cations and the methyl thiocyanate when compared with those of the lithium halide salts.

Alkali hydroxide (LiOH, NaOH, KOH) in water: Structural and vibrational properties, including neutron scattering results

Structural and vibrational properties of aqueous solutions of alkali hydroxides (LiOH, NaOH, and KOH) are computed using quantum molecular dynamics simulations for solute concentrations ranging between 1 and 10M. Element-resolved partial radial distribution functions, neutron and x-ray structure factors, and angular distribution functions are computed for the three hydroxide solutions as a function of concentration. The vibrational spectra and frequency-dependent conductivity are computed from the Fourier transforms of velocity autocorrelation and current autocorrelation functions. Our results for the structure are validated with the available neutron data for 17M concentration of NaOH in water [Semrouni et al., Phys. Chem. Chem. Phys. 21, 6828 (2019)]. We found that the larger ionic radius [rLi+ Collapse

Key Words Collapse

MESH Headings Collapse

Grants Collapse

Affiliation(s)

Ruru Ma Collaboratory for Advanced Computing and Simulations, University of Southern California, Los Angeles, California 90007-0242, USA
Nitish Baradwaj Collaboratory for Advanced Computing and Simulations, University of Southern California, Los Angeles, California 90007-0242, USA
Ken-Ichi Nomura Collaboratory for Advanced Computing and Simulations, University of Southern California, Los Angeles, California 90007-0242, USA
Aravind Krishnamoorthy Department of Mechanical Engineering, Texas A&M University, College Station, Texas 77843-3123, USA
Rajiv K Kalia Collaboratory for Advanced Computing and Simulations, University of Southern California, Los Angeles, California 90007-0242, USA
Aiichiro Nakano Collaboratory for Advanced Computing and Simulations, University of Southern California, Los Angeles, California 90007-0242, USA
Priya Vashishta Collaboratory for Advanced Computing and Simulations, University of Southern California, Los Angeles, California 90007-0242, USA

Collapse

Effect of Water on a Hydrophobic Deep Eutectic Solvent

Deep eutectic solvents (DESs) formed by hydrogen bond donors and acceptors are a promising new class of solvents. Both hydrophilic and hydrophobic binary DESs readily absorb water, making them ternary mixtures, and a small water content is always inevitable under ambient conditions. We present a thorough study of a typical hydrophobic DES formed by a 1:2 mole ratio of tetrabutyl ammonium chloride and decanoic acid, focusing on the effects of a low water content caused by absorbed water vapor, using multinuclear NMR techniques, molecular modeling, and several other physicochemical techniques. Already very low water contents cause dynamic nanoscale phase segregation, reduce solvent viscosity and fragility, increase self-diffusion coefficients and conductivity, and enhance local dynamics. Water interferes with the hydrogen-bonding network between the chloride ions and carboxylic acid groups by solvating them, which enhances carboxylic acid self-correlation and ion pair formation between tetrabutyl ammonium and chloride. Simulations show that the component molar ratio can be varied, with an effect on the internal structure. The water-induced changes in the physical properties are beneficial for most prospective applications but water creates an acidic aqueous nanophase with a high halide ion concentration, which may have chemically adverse effects.

Interaction Energy Analysis of Monovalent Inorganic Anions in Bulk Water Versus Air/Water Interface

Soft anions exhibit surface activity at the air/water interface that can be probed using surface-sensitive vibrational spectroscopy, but the structural implications of this surface activity remain a matter of debate. Here, we examine the nature of anion-water interactions at the air/water interface using a combination of molecular dynamics simulations and quantum-mechanical energy decomposition analysis based on symmetry-adapted perturbation theory. Results are presented for a set of monovalent anions, including Cl-, Br-, I-, CN-, OCN-, SCN-, NO2-, NO3-, and ClOn- (n=1,2,3,4), several of which are archetypal examples of surface-active species. In all cases, we find that average anion-water interaction energies are systematically larger in bulk water although the difference (with respect to the same quantity computed in the interfacial environment) is well within the magnitude of the instantaneous fluctuations. Specifically for the surface-active species Br-(aq), I-(aq), ClO4-(aq), and SCN-(aq), and also for ClO-(aq), the charge-transfer (CT) energy is found to be larger at the interface than it is in bulk water, by an amount that is greater than the standard deviation of the fluctuations. The Cl-(aq) ion has a slightly larger CT energy at the interface, but NO3-(aq) does not; these two species are borderline cases where consensus is lacking regarding their surface activity. However, CT stabilization amounts to <20% of the total induction energy for each of the ions considered here, and CT-free polarization energies are systematically larger in bulk water in all cases. As such, the role of these effects in the surface activity of soft anions remains unclear. This analysis complements our recent work suggesting that the short-range solvation structure around these ions is scarcely different at the air/water interface from what it is in bulk water. Together, these observations suggest that changes in first-shell hydration structure around soft anions cannot explain observed surface activities.

Structure and dynamics of ions in dipolar solvents: a coarse-grained simulation study

We employ the coarse-grained molecular dynamics simulation to investigate the fundamental structural and dynamic properties of the ionic solution with and without the application of an external electric field. Our simulations, in which the solvent molecules are treated as Stockmayer fluids and the ions are modeled as spheres, can effectively account for the multi-body correlations between ion-ion, ion-dipole, and dipole-dipole interactions, which are often ignored by the mean-field theories or coarse-grained simulations based on a dielectric continuum. By focusing on the coupling between effects of ion solvation, electrostatic correlations and applied electric field, we highlight some nontrivial microscopic molecular features of the systems, such as the reorganization of the dipolar solvent, clustering of the ions, and diffusions of ions and dipolar solvent molecules. Particularly, our simulation indicates the nonmonotonic dependence of the ionic clustering and ion diffusion rates on the dipolar nature of the solvent molecules, as well as the amplification of these tendencies caused by the electric field application. This work provides insights into the fundamental understanding of physicochemical properties for ion-containing liquids and contributes to the design and development of ion-containing materials.

Neutron total scattering investigation on the dissolution mechanism of trehalose in NaOH/urea aqueous solution

Trehalose is chosen as a model molecule to investigate the dissolution mechanism of cellulose in NaOH/urea aqueous solution. The combination of neutron total scattering and empirical potential structure refinement yields the most probable all-atom positions in the complex fluid and reveals the cooperative dynamic effects of NaOH, urea, and water molecules in the dissolution process. NaOH directly interacts with glucose rings by breaking the inter- and intra-molecular hydrogen bonding. Na+, thus, accumulates around electronegative oxygen atoms in the hydration shell of trehalose. Its local concentration is thereby 2-9 times higher than that in the bulk fluid. Urea molecules are too large to interpenetrate into trehalose and too complex to form hydrogen bonds with trehalose. They can only participate in the formation of the hydration shell around trehalose via Na+ bridging. As the main component in the complex fluid, water molecules have a disturbed tetrahedral structure in the presence of NaOH and urea. The structure of the mixed solvent does not change when it is cooled to -12 °C. This indicates that the dissolution may be a dynamic process, i.e., a competition between hydration shell formation and inter-molecule hydrogen bonding determines its dissolution. We, therefore, predict that alkali with smaller ions, such as LiOH, has better solubility for cellulose.

Ion hydration and association in aqueous potassium tetrahydroxyborate solutions

Potassium tetrahydroxyborate solution is a significant material in the borate solution family, but there is limited knowledge about hydration structures and interactions of K⁺, [B(OH)₄^-], and water. In this study, the X-ray diffraction measurements of potassium tetrahydroxyborate solutions have been made. The experimental structure factors are subjected to Empirical Potential Structure Refinement (EPSR) modeling to reveal the details of ion hydration and association in the aqueous solutions. This study shows that the O(W)-O(W) distance of water in the studied solutions ranges from 2.82 to 2.76 Å with a coordination number that ranges from 4.7 ± 1.4 to 3.1 ± 1.3 when the value of the water-salt molar ratio (WSR) is decreased from 30 to 6. The addition of ions slightly affects the tetrahedral structure of water even when the concentration of ions is high. The first hydration distance of K⁺ remained at ∼2.67 Å, whereas the value of the coordination number (CN) decreased from 5.4 ± 1.3 to 3.9 ± 1.5 when the concentration of the borate solution was increased. The hydration ability of K⁺ was weak and almost did not have a fixed local hydration structure. The pair distribution function (PDF) of g_B-O(W)(r) shows that [B(OH)₄^-] has a broad hydration distance from 2.9 to 5.4 Å because of the complex interactive relationship between K⁺, [B(OH)₄^-] and water. There is a competitive hydration between K⁺ and [B(OH)₄^-]. Both the X-ray diffraction and DFT-based calculations confirm that the main species is monodentate contact ion pairs when WSR = 30, bidentate contact ion pairs when WSR = 14, and triple contact ion pairs when WSR = 6. These results will provide a new understanding about potassium tetrahydroxyborate solution.

Multistate Reactive Molecular Dynamics Simulations of Proton Diffusion in Water Clusters and in the Bulk

The molecular mechanics with proton transfer (MMPT) force field is combined with multistate adiabatic reactive molecular dynamics (MS-ARMD) to describe proton transport in the condensed phase. Parametrization for small protonated water clusters based on electronic structure calculations at the MP2/6-311+G(2d,2p) level of theory and refinement by comparing with infrared spectra for a protonated water tetramer yields a force field which faithfully describes the minimum energy structures of small protonated water clusters. In protonated water clusters up to (H₂O)₁₀₀H⁺, the proton hopping rate is around 100 hops/ns. This rate converges for 21 ≤ n ≤ 31, and no further speedup in bulk water is found. This indicates that bulklike behavior requires the solvation of a Zundel motif by ∼25 water molecules, which corresponds to the second solvation sphere. For smaller cluster sizes, the number of available states (i.e., the number of proton acceptors) is too small and slows down proton-transfer rates. The cluster simulations confirm that the excess proton is typically located on the surface. The free-energy surface as a function of the weights of the two lowest states and a configurational parameter suggests that the "special pair" plays a central role in rapid proton transport. The barriers between this minimum-energy structure and the Zundel and Eigen minima are sufficiently low (∼1 kcal/mol, consistent with recent experiments and commensurate with a hopping rate of ∼100/ns or 1 every 10 ps), leading to a highly dynamic environment. These findings are also consistent with recent experiments which find that Zundel-type hydration geometries are prevalent in bulk water.

Ion Pairing in HCl-Water Clusters: From Electronic Structure Investigations to Multiconfigurational Force-Field Development

In the bulk, condensed-phase HCl exists as a dissociated Cl^- ion and a proton that is delocalized over solvating water molecules. However, in the gas phase, HCl is covalent, and even on the introduction of hydrating water molecules, the HCl covalent state dominates small clusters and is relevant at larger clusters including 21 water molecules. Electronic structure calculations (at the MP2 level) and ab initio metadynamics simulations (at the DFT level) have been carried out on HCl-(H₂O)_n clusters with n = 2-22 to investigate distinct solvation environments in clusters from covalent HCl structure, to contact ion pairs and solvent-separated ion pairs. The data were further used to train and validate a multiconfigurational force-field for HCl-water clusters that incorporates covalent HCl states into the MS-EVB3.2 formalism. Additionally, the many-body interaction of the Cl^- ion with water and the excess proton was modeled by the introduction of two geometric three-body terms that incorporates the dominant many-body interaction in an efficient noniterative manner.

Nucleation Curve of Carbon Dioxide Hydrate from a Linear Cooling Ramp Method

Formation of gas hydrate is a first-order phase transition that starts with nucleation. Understanding of nucleation is of interest to many in chemical and petroleum industries, as nucleation, while beneficial in many chemical processes, is detrimental in flow assurance of oil and natural gas pipelines. A primary difficulty in the investigation of gas hydrate nucleation has been the inability of researchers to compare nucleation rates of gas hydrates across various systems of different scales and complexities, which in turn has been limiting the ability of researchers to study the nucleation process itself. In this study, a first-generation high-pressure automated lag time apparatus (HP-ALTA MkI) was used to determine the nucleation curve of structure I (sI) - forming carbon dioxide hydrate. The instrument subjected a quiescent water sample of well-defined dimensions to a large number of linear cooling ramps under isobaric conditions, and detected and recorded carbon dioxide hydrate formation temperature distributions. A survival curve was constructed from the measured ensemble, and a nucleation curve was derived from the survival curve using the empirical model-independent method we had previously reported. The nucleation rate of carbon dioxide hydrate was found to be significantly greater than that of pure methane hydrate or that of natural gas hydrate over the entire range of subcooling investigated. We provide a new physical interpretation of an experimentally determined nucleation curve and, by doing so, solve one of the outstanding puzzles of the HP-ALTA technology.

Tracking Aqueous Proton Transfer by Two-Dimensional Infrared Spectroscopy and ab Initio Molecular Dynamics Simulations

Proton transfer in water is ubiquitous and a critical elementary event that, via proton hopping between water molecules, enables protons to diffuse much faster than other ions. The problem of the anomalous nature of proton transport in water was first identified by Grotthuss over 200 years ago. In spite of a vast amount of modern research effort, there are still many unanswered questions about proton transport in water. An experimental determination of the proton hopping time has remained elusive due to its ultrafast nature and the lack of direct experimental observables. Here, we use two-dimensional infrared spectroscopy to extract the chemical exchange rates between hydronium and water in acid solutions using a vibrational probe, methyl thiocyanate. Ab initio molecular dynamics (AIMD) simulations demonstrate that the chemical exchange is dominated by proton hopping. The observed experimental and simulated acid concentration dependence then allow us to extrapolate the measured single step proton hopping time to the dilute limit, which, within error, gives the same value as inferred from measurements of the proton mobility and NMR line width analysis. In addition to obtaining the proton hopping time in the dilute limit from direct measurements and AIMD simulations, the results indicate that proton hopping in dilute acid solutions is induced by the concerted multi-water molecule hydrogen bond rearrangement that occurs in pure water. This proposition on the dynamics that drive proton hopping is confirmed by a combination of experimental results from the literature.

Resolving local configurational contributions to X-ray and neutron radial distribution functions within solutions of concentrated electrolytes - a case study of concentrated NaOH

Extreme conditions of complex materials often lead to a manifold of local environments that challenge characterization and require new advances at the intersection of modern experimental and theoretical techniques. In this contribution, highly caustic and viscous aqueous NaOD solutions were characterized with a combination of X-ray and neutron radial distribution function (RDF) analyses, molecular dynamics simulations and sub-ensemble analysis. While this system has been the topic of some study, the current work expands upon the state of knowledge regarding the extent to which water is perturbed within this chemically extreme solution. Further, we introduce analyses that goes beyond merely identifying the different local environments (ion solvation and coordination environments) that are present, but toward understanding their relative contributions to the ensemble solution RDF. This integrated approach yields unique insight into the experimental sensitivity of RDFs to changes in local geometries, the composition of solvation environments about ions, and the challenge of experimentally differentiating the ensemble of all superimposed local environments-a feature of increasing importance within the extreme condition of high ionic strength.

Correlated dynamics in aqueous proton diffusion

The aqueous proton displays an anomalously large diffusion coefficient that is up to 7 times that of similarly sized cations. There is general consensus that the proton achieves its high diffusion through the Grotthuss mechanism, whereby protons hop from one molecule to the next. A main assumption concerning the extraction of the timescale of the Grotthuss mechanism from experimental results has been that, on average, there is an equal probability for the proton to hop to any of its neighboring water molecules. Herein, we present ab initio simulations that show this assumption is not generally valid. Specifically, we observe that there is an increased probability for the proton to revert back to its previous location. These correlations indicate that the interpretation of the experimental results need to be re-examined and suggest that the timescale of the Grotthuss mechanism is significantly shorter than was previously thought.

Anticorrelated Contributions to Pre-edge Features of Aluminate Near-Edge X-ray Absorption Spectroscopy in Concentrated Electrolytes

Ion pairing within complex solutions and electrolytes is a difficult phenomenon to measure and investigate, yet it has significant impact upon macroscopic processes, such as crystal formation. Traditional methods of detecting and characterizing ion pairing are sensitive to contact ion pairs, may require minimum concentrations that limit applicability, and can have difficulty in characterizing solutions with many components. Because of its element specificity and sensitivity to local environment, X-ray absorption near edge structure (XANES) is a promising tool for investigating ion pairing in complex solutions. In concentrated sodium aluminate solutions, a shift in the pre-edge shoulder correlated to sodium concentration is observed, and the physical origins of that shift are investigated using energy specific time-dependent density functional theory of subensembles obtained from ab initio molecular dynamics. Two transitions are found to contribute to the pre-edge feature, yet they are anticorrelated with respect to the sodium···aluminate distance. Unexpectedly, this causes Al XANES to be an effective probe for longer-range ion interactions than the traditional counterparts of NMR or vibrational spectroscopies. Given the nature of the transitions involved, this observation may be extended to other systems where ion-ion interactions dominate; however, a complete understanding of the contributing transitions is necessary for accurate analysis of XANES pre-edge features in concentrated electrolytes.

OH Stretching Dynamics in Hydroxide Aqueous Solutions

The concept of ions being either water "structure makers" or water "breakers" seems to be inconsistent with the existence of a critical number of water molecules per ion dictating the properties of an aqueous solution, independent of the ion identity. To investigate this issue, Raman spectra of hydroxide aqueous solutions in the region of the OH stretching mode have been obtained under ambient conditions and at concentrations ranging from extreme dilution to the solubility limit. Spectra have been analyzed with a relatively model-free approach, in terms of a superposition of contributions due to the vibrations of the OH^- ions, with two contributions due to the solvent. One of these latter contributions falls at wavenumbers very close to that of the OH^- stretching band, sharing with it its concentration dependence of the full width at half maximum (FWHM). The other contribution due to the solvent is very broad, with increasing FWHM with increasing ion concentration. In the light of these observations, an interpretation of the Raman spectra, based on the possibility of distinguishing the self and distinct contributions, is proposed. The present analysis is supported by structural data on the same solutions and puts into evidence relevant structural and dynamical changes occurring when the number of water molecules available per solute is below ∼20, irrespective of the ion identity.

Classical Molecular Dynamics with Mobile Protons

An important limitation of standard classical molecular dynamics simulations is the inability to make or break chemical bonds. This restricts severely our ability to study processes that involve even the simplest of chemical reactions, the transfer of a proton. Existing approaches for allowing proton transfer in the context of classical mechanics are rather cumbersome and have not achieved widespread use and routine status. Here we reconsider the combination of molecular dynamics with periodic stochastic proton hops. To ensure computational efficiency, we propose a non-Boltzmann acceptance criterion that is heuristically adjusted to maintain the correct or desirable thermodynamic equilibria between different protonation states and proton transfer rates. Parameters are proposed for hydronium, Asp, Glu, and His. The algorithm is implemented in the program CHARMM and tested on proton diffusion in bulk water and carbon nanotubes and on proton conductance in the gramicidin A channel. Using hopping parameters determined from proton diffusion in bulk water, the model reproduces the enhanced proton diffusivity in carbon nanotubes and gives a reasonable estimate of the proton conductance in gramicidin A.

Decomposition of the Experimental Raman and Infrared Spectra of Acidic Water into Proton, Special Pair, and Counterion Contributions

Textbooks describe excess protons in liquid water as hydronium (H₃O⁺) ions, although their true structure remains lively debated. To address this question, we have combined Raman and infrared (IR) multivariate curve resolution spectroscopy with ab initio molecular dynamics and anharmonic vibrational spectroscopic calculations. Our results are used to resolve, for the first time, the vibrational spectra of hydrated protons and counterions and reveal that there is little ion-pairing below 2 M. Moreover, we find that isolated excess protons are strongly IR active and nearly Raman inactive (with vibrational frequencies of ∼1500 ± 500 cm^-1), while flanking water OH vibrations are both IR and Raman active (with higher frequencies of ∼2500 ± 500 cm^-1). The emerging picture is consistent with Georg Zundel's seminal work, as well as recent ultrafast dynamics studies, leading to the conclusion that protons in liquid water are primarily hydrated by two flanking water molecules, with a broad range of proton hydrogen bond lengths and asymmetries.

Multistate empirical valence bond study of temperature and confinement effects on proton transfer in water inside hydrophobic nanochannels

Microscopic characteristics of an aqueous excess proton in a wide range of thermodynamic states, from low density amorphous ices (down to 100 K) to high temperature liquids under the critical point (up to 600 K), placed inside hydrophobic graphene slabs at the nanometric scale (with interplate distances between 3.1 and 0.7 nm wide) have been analyzed by means of molecular dynamics simulations. Water-proton and carbon-proton forces were modeled with a multistate empirical valence bond method. Densities between 0.07 and 0.02 Å(-3) have been considered. As a general trend, we observed a competition between effects of confinement and temperature on structure and dynamical properties of the lone proton. Confinement has strong influence on the local structure of the proton, whereas the main effect of temperature on proton properties is observed on its dynamics, with significant variation of proton transfer rates, proton diffusion coefficients, and characteristic frequencies of vibrational motions. Proton transfer is an activated process with energy barriers between 1 and 10 kJ/mol for both proton transfer and diffusion, depending of the temperature range considered and also on the interplate distance. Arrhenius-like behavior of the transfer rates and of proton diffusion are clearly observed for states above 100 K. Spectral densities of proton species indicated that in all states Zundel-like and Eigen-like complexes survive at some extent. © 2016 Wiley Periodicals, Inc.

Role of Presolvation and Anharmonicity in Aqueous Phase Hydrated Proton Solvation and Transport

Results from condensed phase ab initio molecular dynamics (AIMD) simulations suggest a proton transfer reaction is facilitated by "presolvation" in which the hydronium is transiently solvated by four water molecules, similar to the typical solvation structure of water, by accepting a weak hydrogen bond from the fourth water molecule. A new version 3.2 multistate empirical valence bond (MS-EVB 3.2) model for the hydrated excess proton incorporating this presolvation behavior is therefore developed. The classical MS-EVB simulations show similar structural properties as those of the previous model but with significantly improved diffusive behavior. The inclusion of nuclear quantum effects in the MS-EVB also provides an even better description of the proton diffusion rate. To quantify the influence of anharmonicity, a second model (aMS-EVB 3.2) is developed using the anharmonic aSPC/Fw water model, which provides similar structural properties but improved spectroscopic responses at high frequencies.

Propensity of Hydrated Excess Protons and Hydroxide Anions for the Air-Water Interface

Significant effort has been undertaken to better understand the molecular details governing the propensity of ions for the air-water interface. Facilitated by computationally efficient reactive molecular dynamics simulations, new and statistically conclusive molecular-scale results on the affinity of the hydrated excess proton and hydroxide anion for the air-water interface are presented. These simulations capture the dynamic bond breaking and formation processes (charge defect delocalization) that are important for correctly describing the solvation and transport of these complex species. The excess proton is found to be attracted to the interface, which is correlated with a favorable enthalpic contribution and consistent with reducing the disruption in the hydrogen bond network caused by the ion complex. However, a recent refinement of the underlying reactive potential energy function for the hydrated excess proton shows the interfacial attraction to be weaker, albeit nonzero, a result that is consistent with the experimental surface tension measurements. The influence of a weak hydrogen bond donated from water to the protonated oxygen, recently found to play an important role in excess hydrated proton transport in bulk water, is seen to also be important for this study. In contrast, the hydroxide ion is found to be repelled from the air-water interface. This repulsion is characterized by a reduction of the energetically favorable ion-water interactions, which creates an enthalpic penalty as the ion approaches the interface. Finally, we find that the fluctuation in the coordination number around water sheds new light on the observed entropic trends for both ions.

Proton transfer in liquid water confined inside graphene slabs

The microscopic structure and dynamics of an excess proton in water constrained in narrow graphene slabs between 0.7 and 3.1 nm wide has been studied by means of a series of molecular dynamics simulations. Interaction of water and carbon with the proton species was modeled using a multistate empirical valence bond Hamiltonian model. The analysis of the effects of confinement on proton solvation structure and on its dynamical properties has been considered for varying densities. The system is organized in one interfacial and a bulk-like region, both of variable size. In the widest interplate separations, the lone proton shows a marked tendency to place itself in the bulk phase of the system, due to the repulsive interaction with the carbon atoms. However, as the system is compressed and the proton is forced to move to the vicinity of graphene walls it moves closer to the interface, producing a neat enhancement of the local structure. We found a marked slowdown of proton transfer when the separation of the two graphene plates is reduced. In the case of lowest distances between graphene plates (0.7 and 0.9 nm), only one or two water layers persist and the two-dimensional character of water structure becomes evident. By means of spectroscopical analysis, we observed the persistence of Zundel and Eigen structures in all cases, although at low interplate separations a signature frequency band around 2500 cm^{-1} suffers a blue shift and moves to characteristic values of asymmetric hydronium ion vibrations, indicating some unstability of the typical Zundel-Eigen moieties and their eventual conversion to a single hydronium species solvated by water.

Molecular simulation of water and hydration effects in different environments: challenges and developments for DFTB based models

We discuss the description of water and hydration effects that employs an approximate density functional theory, DFTB3, in either a full QM or QM/MM framework. The goal is to explore, with the current formulation of DFTB3, the performance of this method for treating water in different chemical environments, the magnitude and nature of changes required to improve its performance, and factors that dictate its applicability to reactions in the condensed phase in a QM/MM framework. A relatively minor change (on the scale of kBT) in the O-H repulsive potential is observed to substantially improve the structural properties of bulk water under ambient conditions; modest improvements are also seen in dynamic properties of bulk water. This simple change also improves the description of protonated water clusters, a solvated proton, and to a more limited degree, a solvated hydroxide. By comparing results from DFTB3 models that differ in the description of water, we confirm that proton transfer energetics are adequately described by the standard DFTB3/3OB model for meaningful mechanistic analyses. For QM/MM applications, a robust parametrization of QM-MM interactions requires an explicit consideration of condensed phase properties, for which an efficient sampling technique was developed recently and is reviewed here. The discussions help make clear the value and limitations of DFTB3 based simulations, as well as the developments needed to further improve the accuracy and transferability of the methodology.

Dynamical aspects of intermolecular proton transfer in liquid water and low-density amorphous ices

The microscopic dynamics of an excess proton in water and in low-density amorphous ices has been studied by means of a series of molecular dynamics simulations. Interaction of water with the proton species was modelled using a multistate empirical valence bond Hamiltonian model. The analysis of the effects of low temperatures on proton diffusion and transfer rates has been considered for a temperature range between 100 and 298 K at the constant density of 1 g cm(-3). We observed a marked slowdown of proton transfer rates at low temperatures, but some episodes are still seen at 100 K. In a similar fashion, mobility of the lone proton gets significantly reduced when temperature decreases below 273 K. The proton transfer in low-density amorphous ice is an activated process with energy barriers between 1-10 kJ/mol depending of the temperature range considered and eventually showing Arrhenius-like behavior. Spectroscopic data indicated the survival of both Zundel and Eigen structures along the whole temperature range, revealed by significant spectral frequency shifts.

Solvation structures of protons and hydroxide ions in water

X-ray Raman spectroscopy (XRS) combined with small-angle x-ray scattering (SAXS) were used to study aqueous solutions of HCl and NaOH. Hydrated structures of H(+) and OH(-) are not simple mirror images of each other. While both ions have been shown to strengthen local hydrogen bonds in the hydration shell as indicated by XRS, SAXS suggests that H(+) and OH(-) have qualitatively different long-range effects. The SAXS structure factor of HCl (aq) closely resembles that of pure water, while NaOH (aq) behaves similar to NaF (aq). We propose that protons only locally enhance hydrogen bonds while hydroxide ions induce tetrahedrality in the overall hydrogen bond network of water.

Lifetimes of excess protons in water using a dissociative water potential

Molecular dynamics simulations using a dissociative water potential were applied to study transport of excess protons in water and determine the applicability of this potential to describe such behavior. While originally developed for gas-phase molecules and bulk liquid water, the potential is transferrable to nanoconfinement and interface scenarios. Applied here, it shows proton behavior consistent with ab initio calculations and empirical models specifically designed to describe proton transport. Both Eigen and Zundel complexes are observed in the simulations showing the Eigen-Zundel-Eigen-type mechanism. In addition to reproducing the short-time rattling of the excess proton between the two oxygens of Zundel complexes, a picosecond-scale lifetime was also found. These longer-lived H3O(+) ions are caused by the rapid conversion of the local solvation structure around the transferring proton from a Zundel-like form to an Eigen-like form following the transfer, effectively severing the path along which the proton can rattle. The migration of H(+) over long times (>100 ps) deviates from the conventional short-time multiexponentially decaying lifetime autocorrelation model and follows the t(-3/2) power-law behavior. The potential function employed here matches many of the features of proton transport observed in ab initio molecular dynamics simulations as well as the highly developed empirical valence bond models, yet is computationally very efficient, enabling longer time and larger systems to be studied.