Hofmeister effects: interplay of hydration, nonelectrostatic potentials, and ion size

The classical Derjaguin-Landau-Verwey-Overbeek (DLVO) theory of colloids, and corresponding theories of electrolytes, are unable to explain ion specific forces between colloidal particles quantitatively. The same is true generally, for surfactant aggregates, lipids, proteins, for zeta and membrane potentials and in adsorption phenomena. Even with fitting parameters the theory is not predictive. The classical theories of interactions begin with continuum solvent electrostatic (double layer) forces. Extensions to include surface hydration are taken care of with concepts like inner and outer Helmholtz planes, and "dressed" ion sizes. The opposing quantum mechanical attractive forces (variously termed van der Waals, Hamaker, Lifshitz, dispersion, nonelectrostatic forces) are treated separately from electrostatic forces. The ansatz that separates electrostatic and quantum forces can be shown to be thermodynamically inconsistent. Hofmeister or specific ion effects usually show up above ≈10(-2) molar salt. Parameters to accommodate these in terms of hydration and ion size had to be invoked, specific to each case. Ionic dispersion forces, between ions and solvent, for ion-ion and ion-surface interactions are not explicit in classical theories that use "effective" potentials. It can be shown that the missing ionic quantum fluctuation forces have a large role to play in specific ion effects, and in hydration. In a consistent predictive theory they have to be included at the same level as the nonlinear electrostatic forces that form the skeletal framework of standard theory. This poses a challenge. The challenges go further than academic theory and have implications for the interpretation and meaning of concepts like pH, buffers and membrane potentials, and for their experimental interpretation. In this article we overview recent quantitative developments in our evolving understanding of the theoretical origins of specific ion, or Hofmeister effects. These are demonstrated through an analysis that incorporates nonelectrostatic ion-surface and ion-ion dispersion interactions. This is based on ab initio ionic polarisabilities, and finite ion sizes quantified through recent ab initio work. We underline the central role of ionic polarisabilities and of ion size in the nonelectrostatic interactions that involve ions, solvent molecules and interfaces. Examples of mechanisms through which they operate are discussed in detail. An ab initio hydration model that accounts for polarisabilities of the tightly held hydration shell of "cosmotropic" ions is introduced. It is shown how Hofmeister effects depend on an interplay between specific surface chemistry, surface charge density, pH, buffer, and counterion with polarisabilities and ion size. We also discuss how the most recent theories on surface hydration combined with hydrated nonelectrostatic potentials may predict experimental zeta potentials and hydration forces.     

Effective interaction potentials for alkali and alkaline earth metal ions in SPC/E water and polarization model of hydrated ions

In the first paper (J. Phys. Chem. B, 2006, 110, 10878), effective ion-ion potentials in SPC/E water were obtained for Me-Me, Me-Cl-, and Cl(-)-Cl- pairs, where Me is Li+, Na+, K+, Mg2+, Ca2+, Sr2+, and Ba2+ cations. In this second part of the study of effective interionic potentials, ion-ion distribution functions obtained from implicit-water Monte Carlo simulations of electrolyte solution with these potentials have been explored. This analysis verifies the range of applicability of the primitive model of electrolyte. It is shown that this approximation can be applied to monovalent electrolyte solutions in a wide range of concentrations, whereas the nature of ion-ion interactions is notably different for 2:1 electrolytes. An improved model of ions is discussed. The model includes approximations of the ion hydration shell polarization and specific short-range ion-ion interaction. It allows approximation of the potential of mean force acting on ions in strong electric fields of highly charged macromolecules and bilayers.     

离子质电比和相差异因子对离子半径的综合标度

根据万有引力势与电势的关系式和系统的质电比(单位电量的质量)Sr的物理意义, 研究了离子半径r与离子的Sr和相差异因子的关系. 对于阳离子, r与lgSr和相差异因子呈线性关系; 对于稳定构型阴离子, r与Sr和相差异因子也存在定量关系. 采用回归分析方法, 给出稳定构型和非稳定构型阳离子半径计算公式, 以及稳定构型阴离子半径计算公式. 从相关系数R和回归方程的显著性检验(F)都可说明r与Sr和相差异因子密切相关, 其中拟合的96种元素的138种阳离子半径数据与具有代表性参考值相比, 平均绝对误差为2.0 pm, 相对误差为2.5%. 并预测出较为合理的稀有气体等30种离子半径数据. 同时给出一条获取离子半径(包括复杂离子)数据的新途径.     

Coarse-graining intermolecular interactions in dispersions of highly charged colloids

Effective pair potentials between charged colloids, obtained from Monte Carlo simulations of two single colloids in a closed cell at the primitive model level, are shown to reproduce accurately the structure of aqueous salt-free colloidal dispersions, as determined from full primitive model simulations by Linse et al. (Linse, P.; Lobaskin, V. Electrostatic Attraction and Phase Separation in Solutions of Like-Charged Colloidal Particles. Phys. Rev. Lett.1999, 83, 4208). Excellent agreement is obtained even when ion-ion correlations are important and is in principle not limited to spherical particles, providing a potential route to coarse-grained colloidal interactions in more complex systems.     

Effective solvent mediated potentials of Na+ and Cl- ions in aqueous solution: temperature dependence

The effective solvent-mediated potentials for Na(+) and Cl(-) ions in aqueous solution were calculated in a wide range of temperatures from 0 to 100 °C. The potentials have been determined using the inverse Monte Carlo approach, from the ion-ion radial distribution functions computed in 50 ns molecular dynamics simulations of ions and explicit water molecules. We further separated the effective potentials into a short-range part and an electrostatic long-range part represented by a coulombic potential with some dielectric permittivity. We adjusted the value of the dielectric permittivity to provide the fastest possible decay of the short-range potentials at larger distances. The obtained temperature dependence of the dielectric permittivity follows well the experimental data. We show also that the largest part of the temperature dependence of the effective potentials can be attributed to the temperature-dependent dielectric permittivity.     

Effective interaction potentials for alkali and alkaline earth metal ions in SPC/E water and prediction of mean ion activity coefficients

The potential of mean force (PMF) acting between two simple ions surrounded by SPC/E water have been determined by molecular dynamics (MD) simulations using a spherical cavity approach. Such effective ion-ion potentials were obtained for Me-Me, Me-Cl-, and Cl(-)-Cl- pairs, where Me is a Li+, Na+, K+, Mg2+, Ca2+, Sr2+, and Ba2+ cation. The ionic sizes estimated from the effective potentials are not pairwise additive, a feature in the frequently used primitive model for electrolytes. The effective potentials were used in Monte Carlo (MC) simulations with implicit water to calculate mean ion activity coefficients of LiCl, NaCl, KCl, MgCl2, CaCl2, SrCl2, and BaCl2. Predicted activities were compared with experimental ones in the electrolyte concentration range 0.1-1 M. A qualitative agreement for LiCl and a satisfactory agreement for NaCl were found, whereas the predictions for KCl by two K+ models were less coherent. In the case of alkaline earth metal ions, all experimental activities were successfully reproduced at c = 0.1 M. However, at higher concentrations, similar deviations occurred for all divalent cations, suggesting that the dependence of the permittivity on the salt concentration and the polarization deficiency arising from the ordering of water molecules in the ion hydration shells are important in such systems.     

Kinetic lattice grand canonical Monte Carlo simulation for ion current calculations in a model ion channel system

An algorithm in which kinetic lattice grand canonical Monte Carlo simulations are combined with mean field theory (KLGCMC/MF) is presented to calculate ion currents in a model ion channel system. In this simulation, the relevant region of the system is treated by KLGCMC simulations, while the rest of the system is described by modified Poisson-Boltzmann mean field theory. Calculation of reaction field due to induced charges on the channel/water and membrane/water boundaries is carried out using a basis-set expansion method [Im and Roux, J. Chem. Phys. 115, 4850 (2001)]. Calculation of ion currents, electrostatic potentials, and ion concentrations, as obtained from the KLGCMC/MF simulations, shows good agreement with Poisson-Nernst-Planck (PNP) theory predictions when the channel and membrane have the same dielectric constant as water. If the channel and membrane have a lower dielectric constant than water, however, there is a considerable difference between the KLGCMC/MF and PNP predictions. This difference is attributed to the reaction field, which is missing in PNP theory. It is demonstrated that the reaction field as well as fixed charges in the channel play key roles in selective ion transport. Limitations and further development of the current KLGCMC/MF approach are also discussed.     

胶体粒子表面有效电荷的实验测量

胶体溶液中带电胶粒的有效电荷是计算粒子间相互作用势的一个重要参数. 在本文中用电导率-粒子数密度关系法和电导滴定法分别研究了七种粒径及表面带电情况均不相同的聚苯乙烯粒子, 结果显示两种测量方法得到的有效电荷数值具有较好的一致性, 误差在7%以内. 同时发现, 经验公式计算的有效电荷差不多是实验值的2倍, 表明文献中的经验公式对于本文所研究的胶体粒子体系不适用.     

A finite expansion method for the calculation and interpretation of molecular electrostatic potentials

Because it is useful to have the molecular electrostatic potential as an element in a complex scheme to assess the toxicity of large molecules, efficient and reliable methods are needed for the calculation and characterization of these potentials. A multicenter multipole expansion of the molecular electron charge density calculated with a limited Gaussian basis set is shown here to have only a finite number of nonzero terms from which the molecular electrostatic potential can be calculated. The discrete contributions to the electrostatic potentials from the terms of this expansion provide a physically meaningful decomposition of the potential and a means for its characterization. With pyrrole as an example, the electrostatic potential calculated from this finite expansion of the electron density is compared to that obtained from exact calculations from the same wave function. Good agreement is obtained at distances greater than 1.5 A from any atom in the molecule. In contrast, rearrangement of the terms into an expansion corresponding only to Mulliken atomic charges and dipoles yields a decomposition that produces electrostatic potentials which agree less well with the exact potential. This discrepancy is attributable to the neglect of terms due to higher moments.     

Solvent-averaged potentials for alkali-, earth alkali-, and alkylammonium halide aqueous solutions

We derive effective, solvent-free ion-ion potentials for alkali-, earth alkali-, and alkylammonium halide aqueous solutions. The implicit solvent potentials are parametrized to reproduce experimental osmotic coefficients. The modeling approach minimizes the amount of input required from atomistic (force field) models, which usually predict large variations in the effective ion-ion potentials at short distances. For the smaller ion species, the reported potentials are composed of a Coulomb and a Weeks-Chandler-Andersen term. For larger ions, we find that an additional, attractive potential is required at the contact minimum, which is related to solvent degrees of freedom that are usually not accounted for in standard electrostatics models. The reported potentials provide a simple and accurate force field for use in molecular dynamics and Monte Carlo simulations of (poly-)electrolyte systems.     

A density functional theory approach to noncovalent interactions via interacting monomer densities

A recently proposed "DFT + dispersion" treatment (Rajchel et al., Phys. Rev. Lett., 2010, 104, 163001) is described in detail and illustrated by more examples. The formalism derives the dispersion-free density functional theory (DFT) interaction energy and combines it with the dispersion energy from separate DFT calculations. It consists of the self-consistent polarization of DFT monomers restrained by the exclusion principle via the Pauli blockade technique. Within the monomers a complete exchange-correlation potential should be used, but between them only the exact exchange operates. The application to a wide range of molecular complexes from rare-gas dimers to hydrogen-bonds to π-electron interactions shows good agreement with benchmark values.     

BROMOC-D: Brownian Dynamics/Monte-Carlo Program Suite to Study Ion and DNA Permeation in Nanopores

A theoretical framework is presented to model ion and DNA translocation across a nanopore confinement under an applied electric field. A combined Grand Canonical Monte Carlo Brownian Dynamics (GCMC/BD) algorithm offers a general approach to study ion permeation through wide molecular pores with a direct account of ion-ion and ion-DNA correlations. This work extends previously developed theory by incorporating the recently developed coarse-grain polymer model of DNA by de Pablo and colleagues [Knotts, T. A.; Rathore, N.; Schwartz, D. C.; de Pablo, J. J. J. Chem. Phys. 2007, 126] with explicit ions for simulations of polymer dynamics. Atomistic MD simulations were used to guide model developments. The power of the developed scheme is illustrated with studies of single-stranded DNA (ss-DNA) oligomer translocation in two model cases: a cylindrical pore with a varying radius and a well-studied experimental system, the staphylococcal α-hemolysin channel. The developed model shows good agreement with experimental data for model studies of two homopolymers: ss-poly(dA)(n) and ss-poly(dC)(n). The developed protocol allows for direct evaluation of different factors (charge distribution and pore shape and size) controlling DNA translocation in a variety of nanopores.     

Explanation of Ionic Sequences in Various Phenomena. II. Reversal of Colloid Charge and Ion Binding

Abstract

Previous proposed models for the structure of hydrated ions and the calculated values of the effective dielectric constant of such hydrated ions were used to explain the reversal of colloid charge and ion binding phenomena. In contrast to the conclusions made by Bungenberg de Jong, it is shown that the more soluble the counter -cation or counter-anion of a colloid charge, the greater is the ability of the counterion to reverse the electrical charge of the colloid. The reversal of charge phenomenon is therefore associated with the counterions's solubility, not its insolubility. The solubility sequence is determined by whether or not the carboxylate, sulfate, or phosphate ion is positively (A regions) or negatively hydrated. The phosphate group of DNA or RNA must be associated with a base by means of ion-ion bonds in order to produce the observed reversal of charge sequence. Just as in the reversal of charge phenomenon, the ion-binding phenomenon involves the electrostatic attraction of a counterion with the polyelectrolyte rather than a binding or insolubilization of the counterion. The reverse ion-binding sequence can be obtained if one dialyzes extensively in the presence of sufficient salt before physical measurements are made. This is because the solubility of a counterion determines the true electrostatic charge of the polymer. In other words, different concentrations of salt arise in the dialysis bag when different counter-ions are added because the activity coefficient of the counterion is determined by the solubility of the ion-ion complex between the counterion and the colloid′s charged group.     

Osmotic coefficients of atomistic NaCl (aq) force fields

Solvated ions are becoming increasingly important for (bio)molecular simulations. But there are not much suitable data to validate the intermediate-range solution structure that ion-water force fields produce. We compare six selected combinations of four biomolecular Na-Cl force fields and four popular water models by means of effective ion-ion potentials. First we derive an effective potential at high dilution from simulations of two ions in explicit water. At higher ionic concentration multibody effects will become important. We propose to capture those by employing a concentration dependent dielectric permittivity. With the so obtained effective potentials we then perform implicit solvent simulations. We demonstrate that our effective potentials accurately reproduce ion-ion coordination numbers and the local structure. They allow us furthermore to calculate osmotic coefficients that can be directly compared with experimental data. We show that the osmotic coefficient is a sensitive and accurate measure for the effective ion-ion interactions and the intermediate-range structure of the solution. It is therefore a suitable and useful quantity for validating and parametrizing atomistic ion-water force fields.     

计算脂肪族胺、醇和醚的气相碱性的经验公式

本文给出计算脂肪族胺、醇和醚气相碱性的经验公式。由公式得出的PA值与ICR实验值接近。     

Modeling the electrophoresis of peptides and proteins: improvements in the "bead method" to include ion relaxation and "finite size effects"

A bead model methodology developed in our lab (Xin et al. J. Phys. Chem. B 2006, 110, 1038) and applicable to modeling the free solution electrophoretic mobility of peptides and proteins is generalized in two significant ways. First, an approximate account is taken of the relaxation effect, which makes the methodology applicable to more highly charged peptides and proteins than was previously possible. Second, a more accurate account is taken of the finite size of the beads making up the model structure. This improvement makes the method applicable at higher salt concentrations and/or to models consisting of larger sized subunits. The relaxation effect is accounted for by correcting "unrelaxed" mobilities on the basis of model size and average electrostatic surface, or zeta potential. Correction factors are estimated using those of spheres with the same hydrodynamic radius and zeta potential as the model structure. The correction factors of spheres are readily determined. The more general methodology is first applied to two sets of peptides (74 different peptides total) varying in size from 2 to 42 amino acids. The sets also cover a wide range of net charges. It is shown that accounting for finite bead size results in a small change in model mobilities under the conditions of the experiments (35 mM monovalent salt). The correction for ion relaxation, however, can be significant for highly charged peptides and improves agreement between model and experimental mobilities. Our correction procedure is also tested by examining the electrophoretic mobility of a particular protein "charge ladder" (Carbeck et al. J. Am. Chem. Soc. 1999, 121, 10,671), where the protein charge is varied over a wide range yet the conformation remains essentially constant. In summary, the effects of ion relaxation can be significant if the absolute electrophoretic mobility of a peptide exceeds approximately 0.20 cm2/(kV s).     

An Exact Solution to the Electrostatic Interaction between an Ion-Penetrable Sphere and an Ion-Penetrable Rod

No exact solution for the free energy of electrostatic interaction for a charged sphere and rod geometry in an electrolyte solution has yet been proposed. This geometry is interesting because it can be applied to describe macromolecules interacting with a random fiber-matrix for modeling of hindered transport in diffusional systems. Here we present an analytical approach that yields an exact solution to the problem for ion-penetrable-also called "soft"-sphere and infinitely long rod. This solution is compared to a published finite-element analysis of the same system with nonpenetrable-also called "hard"-sphere and infinitely long rod maintaining a constant surface charge density restriction. For any ionic strength or ratio of rod radius to sphere radius the ion-penetrable method yields an electrostatic free energy of interaction which is lower than that given by the analysis for hard bodies. This free energy is significantly lower for most parameter value combinations and therefore suggests that one should carefully examine the system being modeled to determine if it is approximated better by a hard body or ion-penetrable body approach. Copyright 2000 Academic Press.     

The solvation of Li and Na in acetonitrile from ab initio-derived many-body ion–solvent potentials

Several Li⁺- and Na⁺-acetonitrile models were derived from ab initio calculations at the counterpoise-corrected MP2/TZV++(d,p) level for distorted ion-(MeCN)_n clusters with n=1, 4 and 6. Two different many-body ion-acetonitrile models were constructed: an effective three-body potential for use with the six-site effective pair model of Böhm et al., and an effective polarizable many-body model. The polarizable acetonitrile model used in the latter model is a new empirical model which was also derived in the present paper. Mainly for comparative purposes, two ion-acetonitrile pair potentials were also constructed from the ab initio cluster calculations: one pure pair potential and one effective pair potential. Using all these potential models, MD simulations in the NPT ensemble were performed for the pure acetonitrile liquid and for Li⁺(MeCN) and Na⁺(MeCN) solutions with 1 ion in 512 solvent molecules and with a simulation time of at least 120 ps per system. Thermodynamic properties, solvation-shell structure and the self-diffusion coefficient of the ions and of the solvent molecules were calculated and compared between the different models and with experimental data, where available. The Li⁺ ion is found to be four-coordinated when the new many-body potentials are used, in contrast to the six-coordinated structure obtained for the pure pair and effective pair potentials. The coordination number of Na⁺ is close to six for all the models derived here, although the coordination number becomes slightly smaller with the many-body potentials. For both ions, the solvent molecules in the first shell point their nitrogen ends towards the cation, while in the second shell the opposite orientation is the most common.     

Extraction of site–site bridge functions and effective pair potentials from simulations of polar molecular liquids

We develop an efficient method to extract site–site bridge functions from molecular simulations. The method is based on the inverse solution of the reference site interaction model. Using the exact long‐range asymptotics of site–site direct correlation functions defined by the site–site Ornstein–Zernike equations, we regularize the ill‐posed inverse problem, and then calculate site–site bridge functions and effective pair potentials for ambient water, methanol, and ethanol. We have tested the proposed algorithm and checked its performance. Our study has revealed various peculiarities of the site–site bridge functions, such as long‐range behavior, strong dependence on the electrostatic interactions. Using the obtained data, we have calculated thermodynamic properties of the solvents, namely, isothermal compressibility, internal energy, and Kirkwood‐Buff integrals. The obtained values are in excellent agreement not only with molecular simulations but also with available experimental data. Further extensions of the method are discussed. © 2014 Wiley Periodicals, Inc.