Structure and self-diffusivity of mixed-cation electrolytes between neutral and charged graphene sheets.

Graphene-based applications, such as supercapacitors or capacitive deionization, take place in an aqueous environment, and they benefit from molecular-level insights into the behavior of aqueous electrolyte solutions in single-digit graphene nanopores with a size comparable to a few molecular diameters. Under single-digit graphene nanoconfinement (smallest dimension <2 nm), water and ions behave drastically different than in the bulk. Most aqueous electrolytes in the graphene-based applications as well as in nature contain a mix of electrolytes. We study several prototypical aqueous mixed alkali-chloride electrolytes containing an equimolar fraction of Li/Na, Li/K, or Na/K cations confined between neutral and positively or negatively charged parallel graphene sheets. The strong hydration shell of small Li+ vs a larger Na+ or large K+ with weaker or weak hydration shells affects the interplay between the ions's propensity to hydrate or dehydrate under the graphene nanoconfinement and the strength of the ion-graphene interactions mediated by confinement-induced layered water. We perform molecular dynamics simulations of the confined mixed-cation electrolytes using the effectively polarizable force field for electrolyte-graphene systems and focused on a relation between the electrochemical adsorption and structural properties of the water molecules and ions and their diffusion behavior. The simulations show that the one-layer nanoslits have the biggest impact on the ions' adsorption and the water and ions' diffusion. The positively charged one-layer nanoslits only allow for Cl- adsorption and strengthen the intermolecular bonding, which along with the ultrathin confinement substantially reduces the water and Cl- diffusion. In contrast, the negatively charged one-layer nanoslits only allow for adsorption of weakly hydrated Na+ or K+ and substantially break up the non-covalent bond network, which leads to the enhancement of the water and Na+ or K+ diffusion up to or even above the bulk diffusion. In wider nanoslits, cations adsorb closer to the graphene surfaces than Cl-'s with preferential adsorption of a weakly hydrated cation over a strongly hydrated cation. The positive graphene charge has an intuitive effect on the adsorption of weakly hydrated Na+'s or K+'s and Cl-'s and a counterintuitive effect on the adsorption of strongly hydrated Li+'s. On the other hand, the negative surface charge has an intuitive effect on the adsorption of both types of cations and only mild intuitive or counterintuitive effects on the Cl- adsorption. The diffusion of water molecules and ions confined in the wider nanoslits is reduced with respect to the bulk diffusion, more for the positive graphene charge, which strengthened the intermolecular bonding, and less for the negative surface charge, which weakened the non-covalent bond network.

ReaxFF molecular dynamics of graphene oxide/NaCl aqueous solution interfaces.

In this work, the interaction of NaCl aqueous solution with graphene (G), graphene oxide (GO), and graphite oxide (GTO) is studied using the ReaxFF module of Amsterdam Modeling Suite (AMS) software. We consider four models using the NaCl aqueous solution, containing a graphene sheet (G), a single sheet of GO with epoxide and hydroxyl groups on its surface, 4 layers of GO to model GTO, and a bulk NaCl solution as a reference. The structural and dynamical properties of G, GO, and GTO were quantified by analyzing the functional groups, radial distribution functions, density profiles and diffusivities of water and ions. Due to the reactive force field, the systems underwent spontaneous modification of surface functional groups during the first 750 ps after which the structure stabilizes (the energy stabilizes in less than 400 ps). Pristine graphene in contact with the NaCl solution formed hydroxyl groups on the edges, i.e., converted to partially reduced graphene oxide. The epoxy groups (Oe) on the initial GO were rather unstable, leading to a reduction of their number, however, there was an increase in the number of hydroxyl groups (Oh), mainly at the edges. The interactions of NaCl with the carbon-based sheets are rather weak, including GO and GTO which are decorated with numerous functional groups. Diffusion coefficients of water agree with the available data, but discrepancies in Na+ and Cl- diffusivity compared to other references underscore the need for further development in the dynamic parameters of the reactive force field used. In essence, our research provides specific data previously unreported, laying a foundation for advancing water desalination system design. The study's novelty lies in its realistic approach to graphene/graphene oxide modification, comprehensive characterization, and the application of the reactive force field to explore the graphene oxide-NaCl aqueous interface, contributing to the development of a practical membrane system for water desalination.

Structure and self-diffusivity of alkali-halide electrolytes in neutral and charged graphene nanochannels.

Understanding the microscopic behaviour of aqueous electrolyte solutions in graphene-based ultrathin nanochannels is important in nanofluidic applications such as water purification, fuel cells, and molecular sensing. Under extreme confinement (<2 nm), the properties of water and ions differ drastically from those in the bulk phase. We studied the structural and diffusion behaviour of prototypical aqueous solutions of electrolytes (LiCl, NaCl, and KCl) confined in both neutral and positively-, and negatively-charged graphene nanochannels. We performed molecular dynamics simulations of the solutions in the nanochannels with either one, two- or three-layer water structures using the effectively polarisable force field for graphene. We analysed the structure and intermolecular bond network of the confined solutions along with their relation to the self-diffusivity of water and ions. The simulations show that Na and K cations can more easily rearrange their solvation shells under the graphene nanoconfinement and adsorb on the graphene surfaces or dissolve in the confinement-induced layered water than the Li cation. The negative surface charge together with the presence of ions orient water molecules with hydrogens towards the graphene surfaces, which in turn weakens the intermolecular bond network. The one-layer nanochannels have the biggest effect on the water structure and intermolecular bonding as well as on the adsorption of ions with only co-ions entering these nanochannels. The self-diffusivity of confined water is strongly reduced with respect to the bulk water and decreases with diminishing nanochannel heights except for the negatively-charged one-layer nanochannel. The self-diffusivity of ions also decreases with the reducing the nanochannel heights except for the self-diffusivity of cations in the negatively-charged one-layer nanochannel, evidencing cooperative diffusion of confined water and ions. Due to the significant break-up of the intermolecular bond network in the negatively-charged one-layer nanochannel, self-diffusion coefficients of water and cations exceed those for the two- and three-layer nanochannels and become comparable to the bulk values.

Evolution of the electrical double layer with electrolyte concentration probed by second harmonic scattering.

Investigating the electrical double layer (EDL) structure has been a long-standing challenge and has seen the emergence of several sophisticated techniques able to probe selectively the few molecular layers of a solid/water interface. While a qualitative estimation of the thickness of the EDL can be obtained using simple theoretical models, following experimentally its evolution is not straightforward and can be even more complicated in nano- or microscale systems, particularly when changing the ionic concentration by several orders of magnitude. Here, we bring insight into the structure of the EDL of SiO2 nanoparticle suspensions and its evolution with increasing ionic concentration using angle-resolved second harmonic scattering (AR-SHS). Below millimolar salt concentrations, we can successively characterize inner-sphere adsorption, diffuse layer formation, and outer-sphere adsorption. Moreover, we show for the first time that, by appropriately selecting the nanoparticle size, it is possible to retrieve information also in the millimolar range. There, we observe a decrease in the magnitude of the surface potential corresponding to a compression in the EDL thickness, which agrees with the results of several other electroanalytical and optical techniques. Molecular dynamics simulations suggest that the EDL compression mainly results from the diffuse layer compression rather than outer-sphere ions (Stern plane) moving closer to the surface.

The "good," the "bad," and the "hidden" in neutron scattering and molecular dynamics of ionic aqueous solutions.

We characterize a concentrated 7.3 m CaCl2 solution, combining neutron diffraction with chloride isotopic substitution (Cl-NDIS) in null water and molecular dynamics (MD) simulations. We elucidate the solution structure, thermodynamic properties, and extent of ion pairing previously suggested as concentration-dependent and often not observed at lower concentrations. Our Cl-NDIS measurements designate the solvent-shared ion pairing as dominant and the contact ion pairing (CIP) as insignificant even under conditions close to the solubility limit. The MD models parameterized against neutron diffraction with calcium isotopic substitution (Ca-NDIS) overestimate CIP despite successfully reproducing most of the Cl-NDIS signal. This drawback originates from the fact that Ca2+-Cl- interactions were primarily "hidden" in the Ca-NDIS signal due to overlapping with Ca2+-Ow and Ca2+-Hw contributions to the total scattering. Contrary, MD models with moderate CIP and possessing generally good performance at high concentrations fail to reproduce the NDIS measurements accurately. Therefore, the electronic polarization, introduced in most of the recent MD models via scaling ionic charges, resolves some but not all parameterization drawbacks. We conclude that despite improving the quality of MD models "on average," the question "which model is the best" has not been answered but replaced by the question "which model is better for a given research." An overall "good" model can still be inappropriate or, in some instances, "bad" and, unfortunately, produce erroneous results. The accurate interpretation of several NDIS datasets, complemented by MD simulations, can prevent such mistakes and help identify the strengths, weaknesses, and convenient applications for corresponding computational models.

Second Harmonic Scattering Reveals Ion-Specific Effects at the SiO2 and TiO2 Nanoparticle/Aqueous Interface.

Ion-specific effects play a crucial role in controlling the stability of colloidal systems and regulating interfacial processes. Although mechanistic pictures have been developed to explain the electrostatic structure of solid/water colloidal interfaces, ion-specific effects remain poorly understood. Here we quantify the average interfacial water orientation and the electrostatic surface potential around 100 nm SiO2 and TiO2 colloidal particles in the presence of NaCl, RbCl, and CaCl2 using polarimetric angle-resolved second harmonic scattering. We show that these two parameters can be used to establish the ion adsorption mechanism in a low ionic strength regime (<1 mM added salt). The relative differences between salts as a function of the ionic strength demonstrate cation- and surface-specific preferences for inner- vs outer-sphere adsorption. Compared to monovalent Rb+ and Na+, Ca2+ is found to be preferentially adsorbed as outer-sphere on SiO2 surfaces, while a dominant inner-sphere adsorption is observed for Ca2+ on TiO2. Molecular dynamics simulations performed on crystalline SiO2 and TiO2 surfaces support the experimental conclusions. This work contributes to the understanding of the electrostatic environment around colloidal nanoparticles on a molecular level by providing insight into ion-specific effects with micromolar sensitivity.

Surface Potential and Interfacial Water Order at the Amorphous TiO2 Nanoparticle/Aqueous Interface.

Colloidal nanoparticles exhibit unique size-dependent properties differing from their bulk counterpart, which can be particularly relevant for catalytic applications. To optimize surface-mediated chemical reactions, the understanding of the microscopic structure of the nanoparticle-liquid interface is of paramount importance. Here we use polarimetric angle-resolved second harmonic scattering (AR-SHS) to determine surface potential values as well as interfacial water orientation of ∼100 nm diameter amorphous TiO2 nanoparticles dispersed in aqueous solutions, without any initial assumption on the distribution of interfacial charges. We find three regions of different behavior with increasing NaCl concentration. At very low ionic strengths (0-10 µM), the Na+ ions are preferentially adsorbed at the TiO2 surface as inner-sphere complexes. At low ionic strengths (10-100 µM), a distribution of counterions equivalent to a diffuse layer is observed, while at higher ionic strengths (>100 µM), an additional layer of hydrated condensed ions is formed. We find a similar behavior for TiO2 nanoparticles in solutions of different basic pH. Compared to identically sized SiO2 nanoparticles, the TiO2 interface has a higher affinity for Na+ ions, which we further confirm with molecular dynamics simulations. With its ability to monitor ion adsorption at the surface with micromolar sensitivity and changes in the surface potential, AR-SHS is a powerful tool to investigate interfacial properties in a variety of catalytic and photocatalytic applications.

Oxalic Acid Adsorption on Rutile: Molecular Dynamics and ab Initio Calculations.

Detailed analysis of the adsorption of oxalic acid ions, that is, oxalate and hydrogenoxalate, on the rutile (110) surface was carried out using molecular dynamics augmented by free energy calculations and supported by ab initio calculations. The predicted adsorption on perfect nonhydroxylated and hydroxylated surfaces with surface charge density from neutral to +0.208 C/m2 corresponding to pH values of about 6 and 3.7, respectively, agrees with experimental adsorption data and charge-distribution multisite ion complexation model predictions obtained using the most favorable surface complexes identified in our simulations. We found that outer-sphere complexes are the most favorable, owing to strong hydrogen binding of oxalic acid ions with surface hydroxyls and physisorbed water. The monodentate complex, the most stable among inner-sphere complexes, was about 15 kJ/mol higher in energy, but separated by a large energy barrier. Other inner-sphere complexes, including some previously suggested in the literature as likely adsorption structures such as bidentate and chelate complexes, were found to be unstable both by classical and by ab initio modeling. Both the surfaces and (hydrogen)oxalate ions were modeled using charges scaled to 75% of the nominal values in accord with the electronic continuum theory and our earlier parameterization of (hydrogen)oxalate ions, which showed that nominal charges exaggerate ion-water interactions.

Oxalic Acid Adsorption on Rutile: Experiments and Surface Complexation Modeling to 150 °C.

Here, we characterize oxalate adsorption by rutile in NaCl media (0.03 and 0.30 m) and between pH 3 and 10 over a wide temperature range which includes the near hydrothermal regime (10-150 °C). Oxalate adsorption increases with decreasing pH (as is typical for anion binding by metal oxides), but systematic trends with respect to ionic strength or temperature are absent. Surface complexation modeling (SCM) following the CD-MUSIC formalism, and as constrained by molecular modeling simulations and IR spectroscopic results from the literature, is used to interpret the adsorption data. The molecular modeling simulations, which include molecular dynamics simulations supported by free-energy and ab initio calculations, reveal that oxalate binding is outer-sphere, albeit via strong hydrogen bonds. Conversely, previous IR spectroscopic results conclude that various types of inner-sphere complexes often predominate. SCMs constrained by both the molecular modeling results and the IR spectroscopic data were developed, and both fit the adsorption data equally well. We conjecture that the discrepancy between the molecular simulation and IR spectroscopic results is due to the nature of the rutile surfaces investigated, that is, the perfect (110) crystal faces for the molecular simulations and various rutile powders for the IR spectroscopy studies. Although the (110) surface plane is most often dominant for rutile powders, a variety of steps, kinks, and other types of surface defects are also invariably present. Hence, we speculate that surface defect sites may be primarily responsible for inner-sphere oxalate adsorption, although further study is necessary to prove or disprove this hypothesis.

Surface Characterization of Colloidal Silica Nanoparticles by Second Harmonic Scattering: Quantifying the Surface Potential and Interfacial Water Order.

The microscopic description of the interface of colloidal particles in solution is essential to understand and predict the stability of these systems, as well as their chemical and electrochemical reactivity. However, this description often relies on the use of simplified electrostatic mean field models for the structure of the interface, which give only theoretical estimates of surface potential values and do not provide properties related to the local microscopic structure, such as the orientation of interfacial water molecules. Here we apply polarimetric angle-resolved second harmonic scattering (AR-SHS) to 300 nm diameter SiO2 colloidal suspensions to experimentally determine both surface potential and interfacial water orientation as a function of pH and NaCl concentration. The surface potential values and interfacial water orientation change significantly over the studied pH and salt concentration range, whereas zeta-potential (ζ) values remain constant. By comparing the surface and ζ-potentials, we find a layer of hydrated condensed ions present in the high pH case, and for NaCl concentrations ≥1 mM. For milder pH ranges (pH < 11), as well as for salt concentrations <1 mM, no charge condensation layer is observed. These findings are used to compute the surface charge densities using the Gouy-Chapman and Gouy-Chapman-Stern models. Furthermore, by using the AR-SHS data, we are able to determine the preferred water orientation in the layer directly in contact with the silica interface. Molecular dynamics simulations confirm the experimental trends and allow deciphering of the contributions of water layers to the total response.

Modeling of solid-liquid interfaces using scaled charges: rutile (110) surfaces.

Electronic continuum correction (ECC) has been proven to bring significant improvement in the modeling of interactions of ions (especially multivalent) in aqueous solutions. We present a generalization and the first application of this approach to modeling solid-liquid interfaces, which are omnipresent in physical chemistry, geochemistry, and biophysics. Scaling charges of the top layer of surface atoms makes the existing solid models compatible with the ECC models of ions and molecules, allowing the use of modified force fields for a more accurate investigation of interactions of various metal and metal-oxide surfaces with aqueous solutions, including complex biomolecules and multivalent ions. We have reparametrized rutile (110) models with different surface charge densities (from 0 to -0.416 C m-2) and adopted/developed scaled charge force fields for ions, namely Na+, Rb+, Sr2+, and Cl-. A good agreement of the obtained molecular dynamics (MD) data with X-ray experiments and previously reported MD results was observed, but changes in the occupancy of various adsorption sites were observed and discussed in detail.

Force field parametrization of hydrogenoxalate and oxalate anions with scaled charges.

Models of the hydrogenoxalate (bioxalate, charge -1) and oxalate (charge -2) anions were developed for classical molecular dynamics (CMD) simulations and parametrized against ab initio molecular dynamics (AIMD) data from our previous study (Kroutil et al. (2016) J Mol Model 22:210). The interactions of the anions with water were described using charges scaled according to the electronic continuum correction approach with rescaling of nonbonded parameters (ECCR), and those descriptions of anion interactions were found to agree well with relevant AIMD and experimental results. The models with full RESP charges showed excessively strong electrostatic interactions between the solute and water molecules, leading to an overstructured solvation shell around the anions and thus to a diffusion coefficient that was much too low. The effect of charge scaling was more evident for the oxalate dianion than for the hydrogenoxalate anion. Our work provides CMD models for ions of oxalic acid and extends previous studies that showed the importance of ECCR for modeling divalent ions and ions of organic compounds. Graphical abstract The radial distribution function between a water oxygen (Ow) and an oxygen of the oxalate dianion (Ox) significantly improved when scaled charges were applied to the anion.

Molecular Origins of the Zeta Potential.

The zeta potential (ZP) is an oft-reported measure of the macroscopic charge state of solid surfaces and colloidal particles in contact with solvents. However, the origin of this readily measurable parameter has remained divorced from the molecular-level processes governing the underlying electrokinetic phenomena, which limits its usefulness. Here, we connect the macroscopic measure to the microscopic realm through nonequilibrium molecular dynamics simulations of electroosmotic flow between parallel slabs of the hydroxylated (110) rutile (TiO2) surface. These simulations provided streaming mobilities, which were converted to ZP via the commonly used Helmholtz-Smoluchowski equation. A range of rutile surface charge densities (0.1 to -0.4 C/m2), corresponding to pH values between about 2.8 and 9.4, in RbCl, NaCl, and SrCl2 aqueous solutions, were modeled and compared to experimental ZPs for TiO2 particle suspensions. Simulated ZPs qualitatively agree with experiment and show that "anomalous" ZP values and inequalities between the point of zero charge derived from electrokinetic versus pH titration measurements both arise from differing co- and counterion sorption affinities. We show that at the molecular level the ZP arises from the delicate interplay of spatially varying dynamics, structure, and electrostatics in a narrow interfacial region within about 15 Å of the surface, even in dilute salt solutions. This contrasts fundamentally with continuum descriptions of such interfaces, which predict the ZP response region to be inversely related to ionic strength. In reality the properties of this interfacial region are dominated by relatively immobile and structured water. Consequently, viscosity values are substantially greater than in the bulk, and electrostatic potential profiles are oscillatory in nature.

Pt···H Nonclassical Interaction in Water-Dissolved Pt(II) Complexes: Coaction of Electronic Effects with Solvent-Assisted Stabilization.

The structure of the hydration shell of cisplatin, cis-[Pt(NH3)2Cl2], and its aquated derivatives cis-[Pt(NH3)2Cl(H2O)](+), cis-[Pt(NH3)2OH(H2O)](+), and cis-[Pt(NH3)2(H2O)2](2+) were studied by a number of density functional molecular dynamics (DFT-MD) simulations (from 30 to 250 ps) in which Pt(II) complexes were immersed in a periodic box with 72 explicit water molecules. Furthermore, Pt(II) complex-water binding energy curves and full DFT optimizations of clusters derived from the lowest potential energy DFT-MD frames offered a deeper insight into the structure of the first hydration shell and electronic changes connected with the formation of a nonclassical Pt···H-O-H (Pt···Hw) hydrogen bond (inverse hydration). The probability of a Pt···Hw interaction decreases with increasing charge of the platinum complex due to disadvantageous electrostatics. The main stabilization comes from the charge transfer being followed by polarization and dispersion. Ligands form a framework for the network of H-bond interactions between the solvent molecules, which play an important role in the promotion/suppression of the formation of the Pt···Hw interactions. In the +2 charged diaqua complex the Pt···Hw interaction is still attractive but cannot compete with classical H bonds between solvent molecules. Thus, the formation of a Pt···Hw interaction is the result of a suitable solvent H-bonding network and the probability of its incidence decreases with increasing flexibility of the solvent.

Determination of the distance-dependent viscosity of mixtures in parallel slabs using non-equilibrium molecular dynamics.

We generalize a technique for determination of the shear viscosity of mixtures in planar slabs using non-equilibrium computer simulations by applying an external force parallel to the surface generating Poiseuille flow. The distance-dependent viscosity of the mixture, given as a function of the distance from the surface, is determined by analysis of the resulting velocity profiles of all species. We present results for a highly non-ideal water + methanol mixture in the whole concentration range between rutile (TiO(2)) walls. The bulk results are compared to the existing equilibrium molecular dynamics and experimental data while the inhomogeneous viscosity profiles at the interface are interpreted using the structural data and information on hydrogen bonding.

Influence of the environment on the specificity of the Mg(II) binding to uracil.

Interactions of uracil with the Mg(2+) ion were studied theoretically in the gas phase and in solution. The bare Mg(2+) prefers bidentate N-C═O binding sites stabilizing rare keto-enol forms of the base. Hydration and/or phosphate binding of the Mg(2+) ion shield its positive charge, which leads to preference of monodentate binding to the oxygen keto atoms, shifting fully the equilibrium between the tautomers back toward the canonical diketo tautomer. In solution, a direct inner-sphere metal binding to uracil is not clearly advantageous compared to the outer-sphere metal binding. Similar trends were also obtained for the Ca(2+) ion. Results are supported by the natural bond orbital (NBO) and atoms in molecule (AIM) analyses and the combined extended transition-state energy decomposition analysis and natural orbitals for chemical valence (ETS-NOCV).

Comparison of cation adsorption by isostructural rutile and cassiterite.

Macroscopic net proton charging curves for powdered rutile and cassiterite specimens with the (110) crystal face predominant, as a function of pH in RbCl and NaCl solutions, trace SrCl(2) in NaCl, and trace ZnCl(2) in NaCl and Na Triflate solutions, are compared to corresponding molecular-level information obtained from static DFT optimizations and classical MD simulations, as well as synchrotron X-ray methods. The similarities and differences in the macroscopic charging behavior of rutile and cassiterite largely reflect the cation binding modes observed at the molecular level. Cation adsorption is primarily inner-sphere on both isostructural (110) surfaces, despite predictions that outer-sphere binding should predominate on low bulk dielectric constant oxides such as cassiterite (ε(bulk) ≈ 11). Inner-sphere adsorption is also significant for Rb(+) and Na(+) on neutral surfaces, whereas Cl(-) binding is predominately outer-sphere. As negative surface charge increases, relatively more Rb(+), Na(+), and especially Sr(2+) are bound in highly desolvated tetradentate fashion on the rutile (110) surface, largely accounting for enhanced negative charge development relative to cassiterite. Charging curves in the presence of Zn(2+) are very steep but similar for both oxides, reflective of Zn(2+) hydrolysis (and accompanying proton release) during the adsorption process, and the similar binding modes for ZnOH(+) on both surfaces. These results suggest that differences in cation adsorption between high and low bulk dielectric constant oxides are more subtly related to the relative degree of cation desolvation accompanying inner-sphere binding (i.e., more tetradentate binding on rutile), rather than distinct inner- and outer-sphere adsorption modes. Cation desolvation may be favored at the rutile (110) surface in part because inner-sphere water molecules are bound further from and less tightly than on the cassiterite (110) surface. Hence, their removal upon inner-sphere cation binding is relatively more favorable.

Surface protonation at the rutile (110) interface: explicit incorporation of solvation structure within the refined MUSIC model framework.

The detailed solvation structure at the (110) surface of rutile (alpha-TiO2) in contact with bulk liquid water has been obtained primarily from experimentally verified classical molecular dynamics (CMD) simulations of the ab initio-optimized surface in contact with SPC/E water. The results are used to explicitly quantify H-bonding interactions, which are then used within the refined MUSIC model framework to predict surface oxygen protonation constants. Quantum mechanical molecular dynamics (QMD) simulations in the presence of freely dissociable water molecules produced H-bond distributions around deprotonated surface oxygens very similar to those obtained by CMD with nondissociable SPC/E water, thereby confirming that the less computationally intensive CMD simulations provide accurate H-bond information. Utilizing this H-bond information within the refined MUSIC model, along with manually adjusted Ti-O surface bond lengths that are nonetheless within 0.05 A of those obtained from static density functional theory (DFT) calculations and measured in X-ray reflectivity experiments (as well as bulk crystal values), give surface protonation constants that result in a calculated zero net proton charge pH value (pHznpc) at 25 degrees C that agrees quantitatively with the experimentally determined value (5.4+/-0.2) for a specific rutile powder dominated by the (110) crystal face. Moreover, the predicted pHznpc values agree to within 0.1 pH unit with those measured at all temperatures between 10 and 250 degrees C. A slightly smaller manual adjustment of the DFT-derived Ti-O surface bond lengths was sufficient to bring the predicted pHznpcvalue of the rutile (110) surface at 25 degrees C into quantitative agreement with the experimental value (4.8+/-0.3) obtained from a polished and annealed rutile (110) single crystal surface in contact with dilute sodium nitrate solutions using second harmonic generation (SHG) intensity measurements as a function of ionic strength. Additionally, the H-bond interactions between protolyzable surface oxygen groups and water were found to be stronger than those between bulk water molecules at all temperatures investigated in our CMD simulations (25, 150 and 250 degrees C). Comparison with the protonation scheme previously determined for the (110) surface of isostructural cassiterite (alpha-SnO2) reveals that the greater extent of H-bonding on the latter surface, and in particular between water and the terminal hydroxyl group (Sn-OH) results in the predicted protonation constant for that group being lower than for the bridged oxygen (Sn-O-Sn), while the reverse is true for the rutile (110) surface. These results demonstrate the importance of H-bond structure in dictating surface protonation behavior, and that explicit use of this solvation structure within the refined MUSIC model framework results in predicted surface protonation constants that are also consistent with a variety of other experimental and computational data.

Electric double layer at metal oxide surfaces: static properties of the cassiterite-water interface.

The structure of water at the (110) surface of cassiterite (alpha-SnO2) at ambient conditions was studied by means of molecular dynamics simulations and X-ray crystal truncation rod experiments and interpreted with the help of the revised MUSIC model of surface protonation. The interactions of the metal oxide in the simulations were described by a recently developed classical force field based on the SPC/E model of water. Two extreme cases of completely hydroxylated and nonhydroxylated surfaces were considered along with a mixed surface with 50% dissociation. To study the dependence of the surface properties on pH, neutral and negatively charged variants of the surfaces were constructed. Axial and lateral density distributions of water for different types of surfaces were compared to each other and to experimental axial density distributions found by X-ray experiments. Although significant differences were found between the structures of the studied interfaces, the axial distances between Sn and O atoms are very similar and therefore could not be clearly distinguished by the diffraction technique. The explanation of structures observed in the density distributions was provided by a detailed analysis of hydrogen bonding in the interfacial region. It revealed qualitatively different hydrating patterns formed at neutral hydroxylated and nonhydroxylated surfaces and suggested a preference for the dissociative adsorption of water. At negatively charged surfaces, however, the situation can be reversed by the electric field stabilizing a hydrogen bond network similar to that found at the neutral nonhydroxylated surface. Comparison with previously studied rutile (alpha-TiO2) surfaces provided insight into the differences between the hydration of these two metal oxides, and an important role was ascribed to their different lattice parameters. A link to macroscopic properties was provided by the revised MUSIC surface protonation model. Explicit use of the Sn-O bond lengths based on ab initio calculations and H-bond configurations as inputs led to the prediction of a pH of zero net-proton induced surface charge (pHpzc) that agrees very well with those determined experimentally (about 4.4 at 298 K).

From dimer to condensed phases at extreme conditions: accurate predictions of the properties of water by a Gaussian charge polarizable model.

Water exhibits many unusual properties that are essential for the existence of life. Water completely changes its character from ambient to supercritical conditions in a way that makes it possible to sustain life at extreme conditions, leading to conjectures that life may have originated in deep-sea vents. Molecular simulation can be very useful in exploring biological and chemical systems, particularly at extreme conditions for which experiments are either difficult or impossible; however this scenario entails an accurate molecular model for water applicable over a wide range of state conditions. Here, we present a Gaussian charge polarizable model (GCPM) based on the model developed earlier by Chialvo and Cummings [Fluid Phase Equilib. 150, 73 (1998)] which is, to our knowledge, the first that satisfies the water monomer and dimer properties, and simultaneously yields very accurate predictions of dielectric, structural, vapor-liquid equilibria, and transport properties, over the entire fluid range. This model would be appropriate for simulating biological and chemical systems at both ambient and extreme conditions. The particularity of the GCPM model is the use of Gaussian distributions instead of points to represent the partial charges on the water molecules. These charge distributions combined with a dipole polarizability and a Buckingham exp-6 potential are found to play a crucial role for the successful and simultaneous predictions of a variety of water properties. This work not only aims at presenting an accurate model for water, but also at proposing strategies to develop classical accurate models for the predictions of structural, dynamic, and thermodynamic properties.