Brownian dynamics simulation of a model simple electrolyte in solvents of low dielectric constant

Brownian dynamics simulations are performed to investigate the ionic transport of model simple electrolytes, in which ions are interacting with each other through the repulsive core and Coulombic interactions. The equivalent conductivity and self-diffusion coefficient show minima as the function of the number density of ions when the dielectric constant of the solvent is low. Although the minimum of the former is in harmony with various experiments, no experiment has ever been reported on that of the latter. The analysis of time-dependent transport coefficients reveals that the presence of the minima is ascribed to the slow dynamics, rather than to static association models. The inclusion of a model function that resembles the short-range part of the potential of mean force induced by solvent affects the transport coefficients qualitatively, which suggests the importance of solvent-induced potential of mean force in the conduction mechanism of electrolytes in solvents of low dielectric constant.     

Electrostatics of B-DNA in NaCl and CaCl2 solutions: ion size, interionic correlation, and solvent dielectric saturation effects

The epsilon-modified Poisson-Boltzmann (-MPB) equations ( J. Phys. Chem. B, 2007, 111, 5264) have been solved on a three-dimensional grid for an all-atom geometry model of B-DNA. The approach is based on the implicit solvent model including finite sizes of hydrated ions and a dielectric approximation of the ion hydration shell. Results were obtained for the detailed geometry model of B-DNA in dilute and moderately concentrated solutions of NaCl and CaCl(2). All -MPB parameters of ions and dielectric medium were extracted from published results of all-atom molecular dynamics simulations. The study allows evaluations of the ion size, interionic correlation, and the solvent dielectric saturation effects on the ion distributions around DNA. It unambiguously suggests that the difference between the -MPB and Poisson-Boltzmann distributions of ions is low for Na(+) counterions. Such a difference in the case of divalent counterions Ca(2+) is dramatic: the dielectric saturation of the ion hydration shell leads to point-like adsorption of Ca(2+) on the phosphate groups of DNA. The -MPB equations were also applied to calculate the energy of interaction between two B-DNA molecules. Results agree with previously published simulations and experimental data. Some aspects of ion specificity of polyelectrolyte properties are discussed.     

On the interactions of ions with the air/water interface

In the vicinity of a charged interface, the Poisson-Boltzmann approach considers that the ions obey Boltzmann distributions in a mean electrical field that satisfies the Poisson equation. However, the boundary between two dielectrics generates additional interactions between ions and the interface. The traditional models of ion hydration interactions, that assume that water is a homogeneous dielectric, predict that these interactions are repulsive for all kinds of ions, since all ions should prefer the medium with a larger dielectric constant, where they are better hydrated. In reality, the interactions between the ions and the neighboring water molecules can generate additional short-range ion-hydration interactions, which are either repulsive (for structure-making ions) or attractive (for structure-breaking ions). In the present paper, various models for the ion-hydration forces are examined and compared with the results of molecular dynamics simulations. At large ionic strengths, the latter results could be reproduced qualitatively only when short-ranged attractions between the structure-breaking ions and the interface were taken into account.     

Selective coalescence of bubbles in simple electrolytes

Simple ions in electrolytes exhibit different degrees of affinity for the approach to the free surface of water. This results in strong ion-specific effects that are particularly dramatic in the selective inhibition of bubble coalescence. I present here the calculation of electrostatic interaction between free surfaces of electrolytes caused by the ion accumulation or depletion near a surface. When both anion and cation are attracted to the surface (like H+ and Cl- in HCl solutions), van der Waals attraction facilitates approach of the surfaces and the coalescence of air bubbles. When only an anion or cation is attracted to the surface (like Cl- in NaCl solutions), an electric double layer forms, resulting in repulsive interaction between free surfaces. I applied the method of effective potentials (evaluated from published ion density profiles obtained in simulations) to calculate the ionic contribution to the surface-surface interaction in NaCl and HCl solutions. In NaCl, but not in HCl, the double-layer interaction creates a repulsive barrier to the approach of bubbles, in agreement with the experiments. Moreover, the concentration where ionic repulsion in NaCl becomes comparable in magnitude to the short-range hydrophobic attraction corresponds to the experimentally found transition region toward the inhibition of coalescence.     

计算离子液体溶液汽液相平衡的分子热力学模型

用平均球近似理论、微扰理论和UNIFAC基团贡献方法分别考虑离子之间的长程静电作用、离子与溶剂之间的中程静电作用以及所有粒子之间的短程作用,本文提出了一种新的分子热力学模型,可用于离子液体溶液中溶剂活度系数的计算.通过对含烷基咪唑磷酸酯类离子液体与水、甲醇或乙醇组成的9个二元体系的饱和蒸汽压数据进行关联,获得了相关的模型参数,即溶剂的分子直径和基团之间的交互作用能参数.溶剂活度系数及饱和蒸汽压的计算结果与实验值的平均偏差为1.40%,符合良好,因此本模型可望用于含离子液体体系汽液相平衡的预测.     

Development and validation of an integrated computational approach for the study of ionic species in solution by means of effective two-body potentials. The case of Zn2+, Ni2+, and Co2+ in aqueous solutions

In this paper we have developed an effective computational procedure for the structural and dynamical investigation of ions in aqueous solutions. Quantum mechanical potential energy surfaces for the interaction of a transition metal ion with a water molecule have been calculated taking into account the effect of bulk solvent by the polarizable continuum model (PCM). The effective ion-water interactions have been fitted by suitable analytical potentials, and have been utilized in molecular dynamics (MD) simulations to obtain structural and dynamical properties of the ionic aqueous solutions. This procedure has been successfully applied to the Co2+-H2O open-shell system and, for the first time, Co-oxygen and Co-hydrogen pair potential functions have been determined and employed in MD simulations. The reliability of the whole procedure has been assessed by applying it also to the Zn2+ and Ni2+ aqueous solutions, and the structural and dynamical properties of the three systems have been calculated by means of MD simulations and have been found to be in very good agreement with experimental results. The structural parameters of the first solvation shells issuing from the MD simulations provide an effective complement to extended X-ray absorption fine structure (EXAFS) experiments.     

Insights into Tris‐(2‐Hydroxylethyl)methylammonium Methylsulfate Aqueous Solutions

We report herein a combined experimental–computational study on tris‐(2‐hydroxylethyl)methylammonium methylsulfate in water solutions, as a representative ionic liquid of the aqueous‐solution behavior of hydroxylammonium‐based ionic liquids. Relevant thermophysical properties were measured as a function of mixture composition and temperature. Classical molecular dynamics simulations were performed to infer microscopic structural features. The reported results for ionic liquid in water‐rich solutions show that it behaves as isolated non‐interacting ions solvated by water molecules, through well‐defined solvation shells, exerting a disrupting effect on the water hydrogen bonding network. Nevertheless, as ionic liquid concentration increase, interionic association increases, even for diluted water solutions, evolving from the typical behavior of strong electrolytes in solution toward large interacting structures. For ionic‐liquid‐rich mixtures, water exerts a minor disrupting effect on the fluid’s structuring because it occupies regions around each ion (developing water–ion hydrogen bonds) but without significantly weakening anion–cation interactions.     

Transport coefficients of aqueous dodecyltrimethylammonium bromide solutions: comparison between experiments, analytical calculations and numerical simulations

We study dynamical properties of ionic species in aqueous solutions of dodecyltrimethylammonium bromide, for several concentrations below and above the critical micellar concentration (cmc). New experimental determinations of the electrical conductivity are given which are compared to results obtained from an analytical transport theory; transport coefficients of ions in these solutions above the cmc are also computed from Brownian dynamics simulations. Analytical calculations as well as the simulation treat the solution within the framework of the continuous solvent model. Above the cmc, three ionic species are considered: the monomer surfactant, the micelle and the counterion. The analytical transport theory describes the structural properties of the electrolyte solution within the mean spherical approximation and assumes that the dominant forces which determine the deviations of transport processes from the ideal behavior (i.e., without any interactions between ions) are hydrodynamic interactions and electrostatic relaxation forces. In the simulations, both direct interactions and hydrodynamic interactions between solutes are taken into account. The interaction potential is modeled by pairwise repulsive 1/r(12) interactions and Coulomb interactions. The input parameters of the simulation (radii and self-diffusion coefficients of ions at infinite dilution) are partially obtained from the analytical transport theory which fits the experimental determinations of the electrical conductivity. Both the electrical conductivity of the solution and the self-diffusion coefficients of each species computed from Brownian dynamics are compared to available experimental data. In every case, the influence of hydrodynamic interactions (HIs) on the transport coefficients is investigated. It is shown that HIs are crucial to obtain agreement with experiments. In particular, the self-diffusion coefficient of the micelle, which is the largest and most charged species in the present system, is enhanced when HIs are included whereas the diffusion coefficients of the monomer and the counterion are roughly not influenced by HIs.     

Effects of dielectric saturation and ionic screening on the proton self-diffusion coefficients in perfluorosulfonic acid membranes

Proton transport in perfluorosulfonic acid (PFSA) membranes is investigated through a statistical mechanical model that includes the effects of the interaction of the tethered sulfonate groups with both the water and solvated protons. We first derive a potential that describes the electrostatic field due to the dissociated sulfonic acid groups by extending the work of Gronbech-Jensen et al. [ Mol. Phys. 92, 941 (1997)] to a finite array of point charges. A highly convergent series is obtained which includes the effects of screening due to the protons. We then investigate the effects of both dielectric saturation and two distinct formulations of ionic screening on the proton self-diffusion coefficient in Nafion membranes over a range of water contents. Our computations show that the two phenomena (i.e., dielectric saturation and ionic screening) under constant temperature conditions result in canceling affects. Our calculations provide a radial dependence of the proton mobility suggesting that the dominant self-diffusion occurs in the central region of the pores, well separated from the sulfonate groups. Through comparison of our calculated diffusion coefficients with the experimental values we derived a slightly smaller average separation distance of the hydronium ion from the sulfonate ions than suggested by either electronic structure calculations or multistate empirical valence bond molecular-dynamics simulations.     

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.     

Theoretical calculations of the ionic strength dependence of the ionic product of water based on a mean spherical approximation

A modified version of the restricted primitive model for electrolyte solutions based on the mean spherical approximation (MSA) is applied to estimate the ionic strength dependence of the ionic product of water in aqueous solutions containing different salts, which are commonly used as background electrolytes (NaCl, KCl, KNO₃, and NaC10₄). The modification involves the use of permittivity of the solvent as concentration-dependent parameter and a single average effective diameter. This is a way of including effects originated from the solvent which do not exist in the primitive model. In the case of potassium nitrate and sodium perchlorate, a complete methodology to calculate the effective diameter and density dependence of the dielectric constant has been proposed and developed. Fits between calculated and experimental pK_w values are possible over wide concentration ranges using a single adjustable parameter, namely, the average hard core diameter of water.     

Nucleophilic substitution reactions at liquid/liquid interfaces: molecular dynamics simulation of a model S(N)1 dissociation reaction at the water/carbon tetrachloride interface

The ionic dissociation step of the nucleophilic substitution reaction t-BuCl --> t-Bu(+) + Cl(-) is studied at the water/carbon tetrachloride interface using molecular dynamics computer simulations. The empirical valence bond approach is used to couple two diabatic states, covalent and ionic, in the electronically adiabatic limit. The umbrella sampling technique is used to calculate the potential of mean force along the reaction coordinate (defined as the t-Bu to Cl distance) at several interface regions of varying distances from the Gibbs dividing surface. We find a significant increase of the ionic dissociation barrier height and of the reaction free energy at the interface relative to bulk water. This is shown to be due to the reduced polarity of the interface which causes a destabilization of the pure ionic state. However, deformation to the neat interface structure in the form of water protrusions into the organic phase may provide partial stabilization of the ionic species. The importance of these structural effects is examined by repeating the calculations with an artificially smooth interface. The destabilization of the ionic state at the interface also manifests itself with a rapid (picosecond time scale) recombination dynamics of the ions to form the parent molecule followed by a slow vibrational relaxation.     

Interactions between colloidal particles in amphiphilic mixtures: a density functional theory study

We present a density functional theory study of interactions between spherical colloidal particles in amphiphile solutions. Theory is found to be in good agreement with previously published molecular dynamics simulations. It is used to analyze the effect of the amphiphile solution bulk density, the chain length, and the solvent mole fraction on the potential of mean force between the particles. The general features of the potential of mean force are rationalized in terms of formation of layers and bilayers of amphiphilic molecules in the intercolloidal gap. Theory yields the same general trends as observed in simulations and in experiments. In particular, the computed mean force changes its character from repulsive to attractive and back to repulsive as the solvent mole fraction is gradually increased.     

Incorporating the excluded solvent volume and surface charges for computing solvation free energy

Gauss's law or Poisson's equation is conventionally used to calculate solvation free energy. However, the near‐solute dielectric polarization from Gauss's law or Poisson's equation differs from that obtained from molecular dynamics (MD) simulations. To mimic the near‐solute dielectric polarization from MD simulations, the first‐shell water was treated as two layers of surface charges, the densities of which are proportional to the electric field at the solvent molecule that is modeled as a hard sphere. The intermediate water was treated as a bulk solvent. An equation describing the solvation free energy of ions using this solvent scheme was derived using the TIP3P water model. © 2013 Wiley Periodicals, Inc.     

Structural and dynamic properties of concentrated alkali halide solutions: a molecular dynamics simulation study

The physicochemical properties of alkali halide solutions have long been attributed to the collective interactions between ions and water molecules in the solution, yet the structure of water in these systems and its effect on the equilibrium and dynamic properties of these systems are not clearly understood. Here, we present a systematic view of water structure in concentrated alkali halide solutions using molecular dynamics simulations. The results of the simulations show that the size of univalent ions in the solution has a significant effect on the dynamics of ions and other transport properties such as the viscosity that are correlated with the structural properties of water in aqueous ionic solution. Small cations (e.g., Li+) form electrostatically stabilized hydrophilic hydration shells that are different from the hydration shells of large ions (e.g., Cs+) which behave more like neutral hydrophobic particles, encapsulated by hydrogen-bonded hydration cages. The properties of solutions with different types of ion solvation change in different ways as the ion concentration increases. Examples of this are the diffusion coefficients of the ions and the viscosities of solutions. In this paper we use molecular dynamics (MD) simulations to study the changes in the equilibrium and transport properties of LiCl, RbCl, and CsI solutions at concentrations from 0.22 to 3.97 M.     

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.     

Electrolyte solutions between uncharged walls

The anisotropic hypernetted chain approximation is applied to calculate the ion concentration profiles, the pair distribution functions and the surface-surface interaction between two uncharged slabs immersed in a primitive model electrolyte solution. The dielectric constant of the slabs is much smaller than that of the solvent. The full complexity of this problem is considered without further approximations. At short surface separations the ions are largely expelled from the space between the slabs due to repulsive image charge interactions. The pressure results clearly demonstrate the interdependence of the image charge interactions and the dispersion interactions between the slabs.     

On the solvation of ions in small water droplets

The solvations of positively and negatively charged model ions in water droplets have been studied using Monte Carlo simulations performed with a polarizable intermolecular potential function model. Special focus has been placed on the position of the ion in the water droplet. It was found that the sign of the ionic charge is of minor importance but an increased ionic charge localizes the ion to the central regions of the droplet, whereas a large polarizability and a large ionic radius favor locations close to the surface of the water droplet.     

Quantitative characterization of ion pairing and cluster formation in strong 1:1 electrolytes

Aqueous solutions of 1:1 strong electrolytes are considered to be the prototype for complete ionic dissociation. Nonetheless, clustering of strong 1:1 electrolytes has been widely reported in all atom molecular dynamics simulations, and their presence is indirectly implicated in a diverse range of experimental results. Is there a physical basis for nonidealities such as ion pairing and cluster formation in aqueous solutions of strong 1:1 electrolytes? We attempt to answer this question by direct comparison of results from detailed molecular dynamics simulations to experimentally observed properties of 1:1 electrolytes. We report the analysis of a series of lengthy molecular dynamics simulations of alkali-halide solutions carried out over a wide range of physiologically relevant concentrations using explicit representations of water molecules. We find evidence for pronounced nonideal behavior of ions at all concentrations in the form of ion pairs and clusters which are in rapid equilibrium with dissociated ions. The phenomenology for ion pairing seen in these simulations is congruent with the multistep scheme proposed by Eigen and Tamm based on data from ultrasonic absorption experiments. For a given electrolyte, we show that the dependence of cluster populations on concentration can be described through a single set of equilibrium constants. We assess the accuracy of calculated ion pairing constants by favorable comparison to estimates obtained by Fuoss and co-workers and based on conductometric experiments. Ion pairs and clusters form on length scales where the size of individual water molecules is as important as the hard core radius of ions. Ion pairing results as a balance between the favorable Coulomb interactions and the unfavorable partial desolvation of ions needed to form a pair.