Breakdown of linear dielectric theory for the interaction between hydrated ions and graphene

Many vital processes taking place in electrolytes, such as nanoparticle self-assembly, water purification, and the operation of aqueous supercapacitors, rely on the precise many-body interactions between surfaces and ions in water. Here we study the interaction between a hydrated ion and a charge-neutral graphene layer using atomistic molecular dynamics simulations. For small separations, the ion–graphene repulsion is of nonelectrostatic nature, and for intermediate separations, van der Waals attraction becomes important. Contrary to prevailing theory, we show that nonlinear and tensorial dielectric effects become non-negligible close to surfaces, even for monovalent ions. This breakdown of standard isotropic linear dielectric theory has important consequences for the understanding and modeling of charged objects at surfaces

Electrokinetic, electrochemical and electronic surface potentials of the pristine water liquid-vapor interface

Although conceptually simple, the air-water interface displays rich behavior. Different definitions of the electrostatic potential, each relevant for distinct experimental scenarios, lead to widely varying surface potential magnitudes and even different signs. Based on quantum-chemical density-functional molecular dynamics simulations, all relevant surface potentials are evaluated and compared. The spatially averaged surface potential, accessible to electron holography, is dominated by the trace of the water molecular quadrupole moment and amounts to more than + 4 V inside the water phase, very different from results obtained with force-field water models. The surface potential inside a cavity is much smaller, less than 200 mV in magnitude, and depends specifically on the cavity radius. This is the electrochemical surface potential relevant for ion transfer reactions and ion surface adsorption. Charge transfer between water molecules leads to pronounced surface potentials as well. However, when probing electrophoresis by explicitly applying a lateral electric field, the zeta potential turns out to be zero. Thus, charge transfer between water molecules does not translate to a non-zero electrophoretic mobility at the pristine vapor-liquid water interface

Electrokinetic, electrochemical, and electrostatic surface potentials of the pristine water liquid–vapor interface

Although conceptually simple, the air–water interface displays rich behavior and is subject to intense experimental and theoretical investigations. Different definitions of the electrostatic surface potential as well as different calculation methods, each relevant for distinct experimental scenarios, lead to widely varying potential magnitudes and sometimes even different signs. Based on quantum-chemical density-functional-theory molecular dynamics (DFT-MD) simulations, different surface potentials are evaluated and compared to force-field (FF) MD simulations. As well explained in the literature, the laterally averaged electrostatic surface potential, accessible to electron holography, is dominated by the trace of the water molecular quadrupole moment, and using DFT-MD amounts to +4.35 V inside the water phase, very different from results obtained with FF water models which yield negative values of the order of −0.4 to −0.6 V. Thus, when predicting potentials within water molecules, as relevant for photoelectron spectroscopy and non-linear interface-specific spectroscopy, DFT simulations should be used. The electrochemical surface potential, relevant for ion transfer reactions and ion surface adsorption, is much smaller, less than 200 mV in magnitude, and depends specifically on the ion radius. Charge transfer between interfacial water molecules leads to a sizable surface potential as well. However, when probing electrokinetics by explicitly applying a lateral electric field in DFT-MD simulations, the electrokinetic ζ-potential turns out to be negligible, in agreement with predictions using continuous hydrodynamic models. Thus, interfacial polarization charges from intermolecular charge transfer do not lead to significant electrokinetic mobility at the pristine vapor–liquid water interface, even assuming these transfer charges are mobile in an external electric field

Molecular dynamics simulations of the evaporation of hydrated ions from aqueous solution

Although important for atmospheric processes and gas-phase catalysis, very little is known about the hydration state of ions in the vapor phase. Here we study the evaporation energetics and kinetics of a chloride ion from liquid water by molecular dynamics simulations. As chloride permeates the interface, a water finger forms and breaks at a chloride separation of ≈ 2.8 nm from the Gibbs dividing surface. For larger separations from the interface, about 7 water molecules are estimated to stay bound to chloride in saturated water vapor, as corroborated by continuum dielectrics and statistical mechanics models. This ion hydration significantly reduces the free-energy barrier for evaporation. The effective chloride diffusivity in the transition state is found to be about 6 times higher than in bulk, which reflects the highly mobile hydration dynamics as the water finger breaks. Both effects significantly increase the chloride evaporation flux from the quiescent interface of an electrolyte solution, which is predicted from reaction kinetic theory

Influence of polar co-solutes and salt on the hydration of lipid membranes

The influence of the co-solutes TMAO, urea, and NaCl on the hydration repulsion between lipid membranes is investigated in a combined experimental/simulation approach. Pressure–hydration curves obtained via sorption experiments reveal that the repulsion significantly increases when the membranes are loaded with co-solutes, most strongly for TMAO. As a result, the co-solutes retain additional water molecules and therefore provide membranes with a fluid and more physiological environment. The experimental data are quantitatively reproduced in complementary solvent-explicit atomistic molecular dynamics simulations, which yield the chemical potential of water. Simulation analysis reveals that the additional repulsion arises from the osmotic pressure generated by the co-solutes, an effect which is maximal for TMAO, due to its unfavorable interactions with the lipid headgroup layer and its extraordinarily high osmotic coefficient

Universal and Nonuniversal Aspects of Electrostatics in Aqueous Nanoconfinement

Dielectric water properties, which significantly change in confinement, determine electrostatic interactions and thereby influence all molecular forces and chemical reactions. We present comparative simulations of water between graphene sheets, decanol monolayers, and phospholipid and glycolipid bilayers. Generally, dielectric profiles strongly differ in perpendicular and parallel surface directions and for large surface separation decay to the bulk value 1-2 nm away from the surface. Polar surface groups enhance the local interfacial dielectric response and for phospholipid bilayers induce a giant parallel contribution. A mapping on a box model with asymptotically determined effective water layer widths demonstrates that the perpendicular effective dielectric constant for all systems decreases for confinement below a nanometer, while the parallel one stays rather constant. The confinement-dependent perpendicular effective dielectric constant for graphene is in agreement with experimental data only if the effective water layer width is suitably adjusted. The interactions between two charges at small separation depend on the product of parallel and perpendicular effective water dielectric components; for large separation the interactions depend on the confining medium. For metallic confining media the interactions at large separation decay exponentially with a decay length that depends on the ratio of the effective parallel and perpendicular water dielectric components

Consistent description of ion-specificity in bulk and at interfaces by solvent implicit simulations and mean-field theory

Solvent-implicit Monte Carlo (MC) simulations and mean-field theory are used to predict activity coefficients and excess interfacial tensions for NaF, NaCl, NaI, KF, KCl, and KI solutions in good agreement with experimental data over the entire experimentally available concentration range. The effective ionic diameters of the solvent-implicit simulation model are obtained by fits to experimental activity coefficient data. The experimental activity coefficients at high salt concentrations are only reproduced if the ion-specific concentration-dependent decrement of the dielectric constant is included. The dielectric-constant dependent contribution of the single-ion solvation free energy to the activity coefficient is significant and is included. To account for the ion-specific excess interfacial tension of salt solutions, in addition to nonideal solution effects and the salt-concentration-dependent dielectric decrement, an ion-specific ion–interface interaction must be included. This ion–interface interaction, which acts in addition to the dielectric image-charge repulsion, is modeled as a box potential, is considerably more long-ranged than the ion radius, and is repulsive for all ions considered except iodide, in agreement with previous findings and arguments. By comparing different models that include or exclude bulk non-ideal solution behavior, dielectric decrement effects, and ion– interface interaction potentials, we demonstrate how bulk and interfacial ion-specific effects couple and partially compensate each other. Our MC simulations, which correctly include ionic correlations and interfacial dielectric image-charge repulsion, are used to determine effective ion–surface interaction potentials that can be used in a modified Poisson–Boltzmann theory

Consistent description of ion-specificity in bulk and at interfaces by solvent implicit simulations and mean-field theory

Solvent-implicit Monte Carlo (MC) simulations and mean-field theory are used to predict activity coefficients and excess interfacial tensions for NaF, NaCl, NaI, KF, KCl, and KI solutions in good agreement with experimental data over the entire experimentally available concentration range. The effective ionic diameters of the solvent-implicit simulation model are obtained by fits to experimental activity coefficient data. The experimental activity coefficients at high salt concentrations are only reproduced if the ion-specific concentration-dependent decrement of the dielectric constant is included. The dielectric-constant dependent contribution of the single-ion solvation free energy to the activity coefficient is significant and is included. To account for the ion-specific excess interfacial tension of salt solutions, in addition to non-ideal solution effects and the salt-concentration-dependent dielectric decrement, an ion-specific ion–interface interaction must be included. This ion–interface interaction, which acts in addition to the dielectric image-charge repulsion, is modeled as a box potential, is considerably more long-ranged than the ion radius, and is repulsive for all ions considered except iodide, in agreement with previous findings and arguments. By comparing different models that include or exclude bulk non-ideal solution behavior, dielectric decrement effects, and ion–interface interaction potentials, we demonstrate how bulk and interfacial ion-specific effects couple and partially compensate each other. Our MC simulations, which correctly include ionic correlations and interfacial dielectric image-charge repulsion, are used to determine effective ion–surface interaction potentials that can be used in a modified Poisson–Boltzmann theory