Unravelling the solvent response to neutral and charged solutes

Multiscale Modeling of Aqueous Electric Double Layers

From the stability of colloidal suspensions to the charging of electrodes, electric double layers play a pivotal role in aqueous systems. The interactions between interfaces, water molecules, ions and other solutes making up the electrical double layer span length scales from Ångströms to micrometers and are notoriously complex. Therefore, explaining experimental observations in terms of the double layer's molecular structure has been a long-standing challenge in physical chemistry, yet recent advances in simulations techniques and computational power have led to tremendous progress. In particular, the past decades have seen the development of a multiscale theoretical framework based on the combination of quantum density functional theory, force-field based simulations and continuum theory. In this Review, we discuss these theoretical developments and make quantitative comparisons to experimental results from, among other techniques, sum-frequency generation, atomic-force microscopy, and electrokinetics. Starting from the vapor/water interface, we treat a range of qualitatively different types of surfaces, varying from soft to solid, from hydrophilic to hydrophobic, and from charged to uncharged.

Orientational Ordering in Nano-confined Polar Liquids

Water and other polar liquids exhibit nanoscale structuring near charged interfaces. When a polar liquid is confined between two charged surfaces, the interfacial solvent layers begin to overlap, resulting in solvation forces. Here, we perform molecular dynamics simulations of polar liquids with different dielectric constants and molecular shapes and sizes confined between charged surfaces, demonstrating strong orientational ordering in the nanoconfined liquids. To rationalize the observed structures, we apply a coarse-grained continuum theory that captures the orientational ordering and solvation forces of those liquids. Our findings reveal the subtle behavior of different nanoconfined polar liquids and establish a simple law for the decay length of the interfacial orientations of the solvents, which depends on their molecular size and polarity. These insights shed light on the nature of solvation forces, which are important in colloid and membrane science, scanning probe microscopy, and nano-electrochemistry.

Polar liquids at charged interfaces: A dipolar shell theory

The structure of polar liquids and electrolytic solutions, such as water and aqueous electrolytes, at interfaces underlies numerous phenomena in physics, chemistry, biology, and engineering. In this work, we develop a continuum theory that captures the essential features of dielectric screening by polar liquids at charged interfaces, including decaying spatial oscillations in charge and mass, starting from the molecular properties of the solvent. The theory predicts an anisotropic dielectric tensor of interfacial polar liquids previously studied in molecular dynamics simulations. We explore the effect of the interfacial polar liquid properties on the capacitance of the electrode/electrolyte interface and on hydration forces between two plane-parallel polarized surfaces. In the linear response approximation, we obtain simple formulas for the characteristic decay lengths of molecular and ionic profiles at the interface.

Effects of electrical image potentials and solvation energy on salt ions near a metallic or dielectric wall

Electrical image potentials near a metallic or a dielectric wall of higher dielectric constant than that of the solution are attractive, and therefore, could concentrate salt ions near the wall. In fact, ions in room temperature ionic liquids have been observed to precipitate near a metallic surface (but not near a nonmetallic surface). It will be argued that a likely reason for why precipitation of ions in salt water due to electrical image forces has not as yet been observed is that the solvation of the ions is reduced near the wall. This results in an energy barrier. This reduction occurs because of the large decrease near the wall of the dielectric constant of water normal to the wall. The conditions under which ions are able to get past the resulting energy barrier and concentrate at a solid wall, either as a result of a reduction in this barrier due to screening at high ion concentration or as a result of thermal activation over the barrier will be explored.

Dipolar Poisson models in a dual view

In this work, we study the continuum theories of dipolar-Poisson models. Both the standard dipolar-Poisson model and the dipolar-Poisson-Langevin model, which keeps the dipolar density fixed, are non-convex functionals of the scalar electrostatic potential ϕ. Applying the Legendre transform approach introduced by Maggs [Europhys. Lett. 98, 16012 (2012)], the dual functionals of these models are derived and are given by convex vector-field functionals of the dielectric displacement D and the polarization field P. We compare the convex functionals in P-space to the non-convex functionals in electric field E-space and apply them to the classic problem of the solvation of point-like ions. Since the dipolar-Poisson model does not properly describe polarization saturation, we argue that only the dipolar-Poisson-Langevin functional can be used to provide a nonlinear generalization of the harmonic polarization functional used in the theory of Marcus for the electron transfer rate to nonlinear regimes. We show that the model can be quantitatively parameterized by molecular dynamics simulations.

Ultra-fast charging in aluminum-ion batteries: electric double layers on active anode

With the rapid iteration of portable electronics and electric vehicles, developing high-capacity batteries with ultra-fast charging capability has become a holy grail. Here we report rechargeable aluminum-ion batteries capable of reaching a high specific capacity of 200 mAh g-1. When liquid metal is further used to lower the energy barrier from the anode, fastest charging rate of 104 C (duration of 0.35 s to reach a full capacity) and 500% more specific capacity under high-rate conditions are achieved. Phase boundaries from the active anode are believed to encourage a high-flux charge transfer through the electric double layers. As a result, cationic layers inside the electric double layers responded with a swift change in molecular conformation, but anionic layers adopted a polymer-like configuration to facilitate the change in composition.

Effects of the dielectric properties of the ceramic-solvent interface on the binding of proteins to oxide ceramics: a non-local electrostatic approach

The rapid development of nanoscience and nanotechnology has raised many fundamental questions that significantly impede progress in these fields. In particular, understanding the physicochemical processes at the interface in aqueous solvents requires the development and application of efficient and accurate methods. In the present work we evaluate the electrostatic contribution to the energy of model protein-ceramic complex formation in an aqueous solvent. We apply a non-local (NL) electrostatic approach that accounts for the effects of the short-range structure of the solvent on the electrostatic interactions of the interfacial systems. In this approach the aqueous solvent is considered as a non-ionic liquid, with the rigid and strongly correlated dipoles of the water molecules. We have found that an ordered interfacial aqueous solvent layer at the protein- and ceramic-solvent interfaces reduces the charging energy of both the ceramic and the protein in the solvent, and significantly increases the electrostatic contribution to their association into a complex. This contribution in the presented NL approach was found to be significantly shifted with respect to the classical model at any dielectric constant value of the ceramics. This implies a significant increase of the adsorption energy in the protein-ceramic complex formation for any ceramic material. We show that for several biocompatible ceramics (for example HfO2, ZrO2, and Ta2O5) the above effect predicts electrostatically induced protein-ceramic complex formation. However, in the framework of the classical continuum electrostatic model (the aqueous solvent as a uniform dielectric medium with a high dielectric constant ∼80) the above ceramics cannot be considered as suitable for electrostatically induced complex formation. Our results also show that the protein-ceramic electrostatic interactions can be strong enough to compensate for the unfavorable desolvation effect in the process of protein-ceramic complex formation.

Origin of Asymmetric Solvation Effects for Ions in Water and Organic Solvents Investigated Using Molecular Dynamics Simulations: The Swain Acity-Basity Scale Revisited

The asymmetric solvation of ions can be defined as the tendency of a solvent to preferentially solvate anions over cations or cations over anions, at identical ionic charge magnitudes and effective sizes. Taking water as a reference, these effects are quantified experimentally for many solvents by the relative acity (A) and basity (B) parameters of the Swain scale. The goal of the present study is to investigate the asymmetric solvation of ions using molecular dynamics simulations, and to connect the results to this empirical scale. To this purpose, the charging free energies of alkali and halide ions, and of their hypothetical oppositely charged counterparts, are calculated in a variety of solvents. In a first set of calculations, artificial solvent models are considered that present either a charge or a shape asymmetry at the molecular level. The solvation asymmetry, probed by the difference in charging free energy between the two oppositely charged ions, is found to encompass a term quadratic in the ion charge, related to the different solvation structures around the anion and cation, and a term linear in the ion charge, related to the solvation structure around the uncharged ion-sized cavity. For these simple solvent models, the two terms are systematically counteracting each other, and it is argued that only the quadratic term should be retained when comparing the results of simulations involving physical solvents to experimental data. In a second set of calculations, 16 physical solvents are considered. The theoretical estimates for the acity A are found to correlate very well with the Swain parameters, whereas the correlation for B is very poor. Based on this observation, the Swain scale is reformulated into a new scale involving an asymmetry parameter Σ, positive for acitic solvents and negative for basitic ones, and a polarity parameter Π. This revised scale has the same predictive power as the original scale, but it characterizes asymmetry in an absolute sense, the atomistic simulations playing the role of an extra-thermodynamic assumption, and is optimally compatible with the simulation results. Considering the 55 solvents in the Swain set, it is observed that a moderate basity (Σ between -0.9 and -0.3, related to electronic polarization) represents the baseline for most solvents, while a highly variable acity (Σ between 0.0 and 3.0, related to hydrogen-bond donor capacity modulated by inductive effects) represents a landmark of protic solvents.

Hydration Free Energies of Molecular Ions from Theory and Simulation

We present a theoretical/computational framework for accurate calculation of hydration free energies of ionized molecular species. The method is based on a molecular theory, 3D-RISM, combined with a recently developed pressure correction (PC+). The 3D-RISM/PC+ model can provide ∼3 kcal/mol hydration free energy accuracy for a large variety of ionic compounds, provided that the Galvani potential of water is taken into account. The results are compared with direct atomistic simulations. Several methodological aspects of hydration free energy calculations for charged species are discussed.

Molecular origin of high free energy barriers for alkali metal ion transfer through ionic liquid–graphene electrode interfaces

We study mechanisms of solvent-mediated ion interactions with charged surfaces in ionic liquids by molecular dynamics simulations, in an attempt to reveal the main trends that determine ion–electrode interactions in ionic liquids.

Dielectric response of the water hydration layer around spherical solutes

We calculate the local dielectric function ɛ(r) inside the hydration layer around a spherical solute (i) from molecular dynamics simulations including explicit solutes and (ii) theoretically using the nonlocal dielectric function of bulk water which includes the radial electric field, but not the explicit solute. The observed agreement between the two concepts shows that while ɛ(r) is strongly different from bulk, this difference is not due to restructuring of the hydrogen bond network but is mostly a consequence of the field geometry. The dielectric response differs for anions and cations, yielding a natural explanation for the well-known charge asymmetry of ionic solvation in agreement with experimental data.

Communication: modeling charge-sign asymmetric solvation free energies with nonlinear boundary conditions

We show that charge-sign-dependent asymmetric hydration can be modeled accurately using linear Poisson theory after replacing the standard electric-displacement boundary condition with a simple nonlinear boundary condition. Using a single multiplicative scaling factor to determine atomic radii from molecular dynamics Lennard-Jones parameters, the new model accurately reproduces MD free-energy calculations of hydration asymmetries for: (i) monatomic ions, (ii) titratable amino acids in both their protonated and unprotonated states, and (iii) the Mobley "bracelet" and "rod" test problems [D. L. Mobley, A. E. Barber II, C. J. Fennell, and K. A. Dill, "Charge asymmetries in hydration of polar solutes," J. Phys. Chem. B 112, 2405-2414 (2008)]. Remarkably, the model also justifies the use of linear response expressions for charging free energies. Our boundary-element method implementation demonstrates the ease with which other continuum-electrostatic solvers can be extended to include asymmetry.

Fast Computation of Solvation Free Energies with Molecular Density Functional Theory: Thermodynamic-Ensemble Partial Molar Volume Corrections

Molecular density functional theory (MDFT) offers an efficient implicit-solvent method to estimate molecule solvation free-energies, whereas conserving a fully molecular representation of the solvent. Even within a second-order approximation for the free-energy functional, the so-called homogeneous reference fluid approximation, we show that the hydration free-energies computed for a data set of 500 organic compounds are of similar quality as those obtained from molecular dynamics free-energy perturbation simulations, with a computer cost reduced by 2-3 orders of magnitude. This requires to introduce the proper partial volume correction to transform the results from the grand canonical to the isobaric-isotherm ensemble that is pertinent to experiments. We show that this correction can be extended to 3D-RISM calculations, giving a sound theoretical justification to empirical partial molar volume corrections that have been proposed recently.

Simple electrolyte solutions: comparison of DRISM and molecular dynamics results for alkali halide solutions

Using the dielectrically consistent reference interaction site model (DRISM) of molecular solvation, we have calculated structural and thermodynamic information of alkali-halide salts in aqueous solution, as a function of salt concentration. The impact of varying the closure relation used with DRISM is investigated using the partial series expansion of order-n (PSE-n) family of closures, which includes the commonly used hypernetted-chain equation (HNC) and Kovalenko-Hirata closures. Results are compared to explicit molecular dynamics (MD) simulations, using the same force fields, and to experiment. The mean activity coefficients of ions predicted by DRISM agree well with experimental values at concentrations below 0.5 m, especially when using the HNC closure. As individual ion activities (and the corresponding solvation free energies) are not known from experiment, only DRISM and MD results are directly compared and found to have reasonably good agreement. The activity of water directly estimated from DRISM is nearly consistent with values derived from the DRISM ion activities and the Gibbs-Duhem equation, but the changes in the computed pressure as a function of salt concentration dominate these comparisons. Good agreement with experiment is obtained if these pressure changes are ignored. Radial distribution functions of NaCl solution at three concentrations were compared between DRISM and MD simulations. DRISM shows comparable water distribution around the cation, but water structures around the anion deviate from the MD results; this may also be related to the high pressure of the system. Despite some problems, DRISM-PSE-n is an effective tool for investigating thermodynamic properties of simple electrolytes.

Affine-response model of molecular solvation of ions: Accurate predictions of asymmetric charging free energies

Two mechanisms have been proposed to drive asymmetric solvent response to a solute charge: a static potential contribution similar to the liquid-vapor potential, and a steric contribution associated with a water molecule's structure and charge distribution. In this work, we use free-energy perturbation molecular-dynamics calculations in explicit water to show that these mechanisms act in complementary regimes; the large static potential (∼44 kJ/mol/e) dominates asymmetric response for deeply buried charges, and the steric contribution dominates for charges near the solute-solvent interface. Therefore, both mechanisms must be included in order to fully account for asymmetric solvation in general. Our calculations suggest that the steric contribution leads to a remarkable deviation from the popular "linear response" model in which the reaction potential changes linearly as a function of charge. In fact, the potential varies in a piecewise-linear fashion, i.e., with different proportionality constants depending on the sign of the charge. This discrepancy is significant even when the charge is completely buried, and holds for solutes larger than single atoms. Together, these mechanisms suggest that implicit-solvent models can be improved using a combination of affine response (an offset due to the static potential) and piecewise-linear response (due to the steric contribution).

Profile of the static permittivity tensor of water at interfaces: consequences for capacitance, hydration interaction and ion adsorption

We derive the theoretical framework to calculate the dielectric response tensor and determine its components for water adjacent to hydrophilic and hydrophobic surfaces using molecular dynamics simulations. For the nonpolarizable water model used, linear response theory is found to be applicable up to an external perpendicular field strength of ∼2 V/nm, which is well beyond the experimental dielectric breakdown threshold. The dipole contribution dominates the dielectric response parallel to the interface, whereas for the perpendicular component it is essential to keep the quadrupole and octupole terms. Including the space-dependent dielectric function in a mean-field description of the ion distribution at a single charged interface, we reproduce experimental values of the interfacial capacitance. At the same time, the dielectric function decreases the electrostatic part of the disjoining pressure between two charged surfaces, unlike previously thought. The difference in interfacial polarizability between hydrophilic and hydrophobic surfaces can be quantized in terms of the dielectric dividing surface. Using the dielectric dividing surface and the Gibbs dividing surface positions to estimate the free energy of a single ion close to an interface, ion-specific adsorption effects are found to be more pronounced at hydrophobic surfaces than at hydrophilic surfaces, in agreement with experimental trends.

Integral Equation Theory of Biomolecules and Electrolytes

The so-called three-dimensional version (3D-RISM) can be used to describe the interactions of solvent components (here we treat water and ions) with a chemical or biomolecular solute of arbitrary size and shape. Here we give an overview of the current status of such models, describing some aspects of “pure” electrolytes (water plus simple ions) and of ionophores, proteins and nucleic acids in the presence of water and salts. Here we focus primarily on interactions with water and dissolved salts; as a practical matter, the discussion is mostly limited to monovalent ions, since studies of divalent ions present many difficult problems that have not yet been addressed. This is not a comprehensive review, but covers a few recent examples that illustrate current issues.

Differential geometry based solvation model II: Lagrangian formulation

Solvation is an elementary process in nature and is of paramount importance to more sophisticated chemical, biological and biomolecular processes. The understanding of solvation is an essential prerequisite for the quantitative description and analysis of biomolecular systems. This work presents a Lagrangian formulation of our differential geometry based solvation models. The Lagrangian representation of biomolecular surfaces has a few utilities/advantages. First, it provides an essential basis for biomolecular visualization, surface electrostatic potential map and visual perception of biomolecules. Additionally, it is consistent with the conventional setting of implicit solvent theories and thus, many existing theoretical algorithms and computational software packages can be directly employed. Finally, the Lagrangian representation does not need to resort to artificially enlarged van der Waals radii as often required by the Eulerian representation in solvation analysis. The main goal of the present work is to analyze the connection, similarity and difference between the Eulerian and Lagrangian formalisms of the solvation model. Such analysis is important to the understanding of the differential geometry based solvation model. The present model extends the scaled particle theory of nonpolar solvation model with a solvent-solute interaction potential. The nonpolar solvation model is completed with a Poisson-Boltzmann (PB) theory based polar solvation model. The differential geometry theory of surfaces is employed to provide a natural description of solvent-solute interfaces. The optimization of the total free energy functional, which encompasses the polar and nonpolar contributions, leads to coupled potential driven geometric flow and PB equations. Due to the development of singularities and nonsmooth manifolds in the Lagrangian representation, the resulting potential-driven geometric flow equation is embedded into the Eulerian representation for the purpose of computation, thanks to the equivalence of the Laplace-Beltrami operator in the two representations. The coupled partial differential equations (PDEs) are solved with an iterative procedure to reach a steady state, which delivers desired solvent-solute interface and electrostatic potential for problems of interest. These quantities are utilized to evaluate the solvation free energies and protein-protein binding affinities. A number of computational methods and algorithms are described for the interconversion of Lagrangian and Eulerian representations, and for the solution of the coupled PDE system. The proposed approaches have been extensively validated. We also verify that the mean curvature flow indeed gives rise to the minimal molecular surface and the proposed variational procedure indeed offers minimal total free energy. Solvation analysis and applications are considered for a set of 17 small compounds and a set of 23 proteins. The salt effect on protein-protein binding affinity is investigated with two protein complexes by using the present model. Numerical results are compared to the experimental measurements and to those obtained by using other theoretical methods in the literature.