Simulation of electric double layers with multivalent counterions: Ion size effect

Effect of Mixed Morphology (Simple Cubic, Face-Centered Cubic, and Body-Centered Cubic)-Based Electrodes on the Electric Double Layer Capacitance of Supercapacitors

Supercapacitors store energy due to the formation of an electric double layer (EDL) at the interface of the electrodes and electrolyte. The present article deals with the finite element study of equilibrium electric double layer capacitance (EDLC) in the mixed morphology electrodes comprising all three fundamental crystal structures, simple cubic (SC), body-centered cubic (BCC), and face-centered cubic morphologies (FCC). Mesoporous-activated carbon forms the electrode in the supercapacitor with (C2H5)4NBF4/propylene carbonate organic electrolyte. Electrochemical interference is clearly demonstrated in the supercapacitors with the formation of the potential bands, as in the case of interference theory due to the increasing packing factor. The effects of electrode thickness varying from a wide range of 50 nm to 0.04 mm on specific EDLC have been discussed in detail. The interfacial geometry of the unit cell in contact with the electrolyte is the most important parameter determining the properties of the EDL. The critical thickness of the electrodes is 1.71 μm in all the morphologies. Polarization increases the interfacial potential and leads to EDL formation. The Stern layer specific capacitance is 167.6 μF cm-2 in all the morphologies. The maximum capacitance is in the decreasing order of interfacial geometry, as FCC > BCC > SC, dependent on the packing factor. The minimum transmittance in all the morphologies is 98.35%, with the constant figure of merit at higher electrode thickness having applications in the chip interconnects. The transient analysis shows that the interfacial current decreases with increasing polarization in the EDL. The capacitance also decreases with the increase of the scan rate.

Competition between excluded-volume and electrostatic interactions for nanogel swelling: effects of the counterion valence and nanogel charge

The equilibrium distribution of monovalent and trivalent ions within a thermo-responsive charged nanogel is investigated using Monte Carlo simulations and Ornstein–Zernike equations.

Specific ion effects on the electrokinetic properties of iron oxide nanoparticles: experiments and simulations

We report experimental and simulation studies on ion specificity in aqueous colloidal suspensions of positively charged, bare magnetite nanoparticles. Magnetite has the largest saturation magnetization among iron oxides and relatively low toxicity, which explain why it has been used in multiple biomedical applications. Bare magnetite is hydrophilic and the sign of the surface charge can be changed by adjusting the pH, its isoelectric point being in the vicinity of pH = 7. Electrophoretic mobility of our nanoparticles in the presence of increasing concentrations of different anions showed that anions regarded as kosmotropic are more efficient in decreasing, and even reversing, the mobility of the particles. If the anions were ordered according to the extent to which they reduced the particle mobility, a classical Hofmeister series was obtained with the exception of thiocyanate, whose position was altered. Monte Carlo simulations were used to predict the diffuse potential of magnetite in the presence of the same anions. The simulations took into account the ion volume, and the electrostatic and dispersion forces among the ions and between the ions and the solid surface. Even though no fitting parameters were introduced and all input data were estimated using Lifshitz theory of van der Waals forces or obtained from the literature, the predicted diffusion potentials of different anions followed the same order as the mobility curves. The results suggest that ionic polarizabilities and ion sizes are to a great extent responsible for the specific ion effects on the electrokinetic potential of iron oxide particles.

Forces between different latex particles in aqueous electrolyte solutions measured with the colloidal probe technique

In this article, a compilation of results on direct force measurements between colloidal particles in monovalent salts carried out with the colloidal probe technique based on Atomic Force Microscopy was presented. The interaction forces between similar and dissimilar particles was studied and it was concluded that, in general, these force profiles may be satisfactorily quantified by the DLVO theory down to distances of few nanometers. However, in the specific case where the charge of one of the involved particle is close to neutral, it was found that the surface potential of this particle may change its sign depending on the sign of charge of the opposite particle. In this respect, the assumption that the surface potential of a particle is a property only related to the particle surface features and the bulk properties is called into question. Microsc. Res. Tech. 80:144-152, 2017. © 2016 Wiley Periodicals, Inc.

Charge Relaxation Dynamics of an Electrolytic Nanocapacitor

Understanding ion relaxation dynamics in overlapping electric double layers (EDLs) is critical for the development of efficient nanotechnology-based electrochemical energy storage, electrochemomechanical energy conversion, and bioelectrochemical sensing devices as well as the controlled synthesis of nanostructured materials. Here, a lattice Boltzmann (LB) method is employed to simulate an electrolytic nanocapacitor subjected to a step potential at t = 0 for various degrees of EDL overlap, solvent viscosities, ratios of cation-to-anion diffusivity, and electrode separations. The use of a novel continuously varying and Galilean-invariant molecular-speed-dependent relaxation time (MSDRT) with the LB equation recovers a correct microscopic description of the molecular-collision phenomena and enhances the stability of the LB algorithm. Results for large EDL overlaps indicated oscillatory behavior for the ionic current density, in contrast to monotonic relaxation to equilibrium for low EDL overlaps. Further, at low solvent viscosities and large EDL overlaps, anomalous plasmalike spatial oscillations of the electric field were observed that appeared to be purely an effect of nanoscale confinement. Employing MSDRT in our simulations enabled modeling of the fundamental physics of the transient charge relaxation dynamics in electrochemical systems operating away from equilibrium wherein Nernst-Einstein relation is known to be violated.

Ionic Size Effects: Generalized Boltzmann Distributions, Counterion Stratification, and Modified Debye Length

Near a charged surface, counterions of different valences and sizes cluster; and their concentration profiles stratify. At a distance from such a surface larger than the Debye length, the electric field is screened by counterions. Recent studies by a variational mean-field approach that includes ionic size effects and by Monte Carlo simulations both suggest that the counterion stratification is determined by the ionic valence-to-volume ratios. Central in the mean-field approach is a free-energy functional of ionic concentrations in which the ionic size effects are included through the entropic effect of solvent molecules. The corresponding equilibrium conditions define the generalized Boltzmann distributions relating the ionic concentrations to the electrostatic potential. This paper presents a detailed analysis and numerical calculations of such a free-energy functional to understand the dependence of the ionic charge density on the electrostatic potential through the generalized Boltzmann distributions, the role of ionic valence-to-volume ratios in the counterion stratification, and the modification of Debye length due to the effect of ionic sizes.

Competitive adsorption and ordered packing of counterions near highly charged surfaces: From mean-field theory to Monte Carlo simulations

Competitive adsorption of counterions of multiple species to charged surfaces is studied by a size-effect-included mean-field theory and Monte Carlo (MC) simulations. The mean-field electrostatic free-energy functional of ionic concentrations, constrained by Poisson's equation, is numerically minimized by an augmented Lagrangian multiplier method. Unrestricted primitive models and canonical ensemble MC simulations with the Metropolis criterion are used to predict the ionic distributions around a charged surface. It is found that, for a low surface charge density, the adsorption of ions with a higher valence is preferable, agreeing with existing studies. For a highly charged surface, both the mean-field theory and the MC simulations demonstrate that the counterions bind tightly around the charged surface, resulting in a stratification of counterions of different species. The competition between mixed entropy and electrostatic energetics leads to a compromise that the ionic species with a higher valence-to-volume ratio has a larger probability to form the first layer of stratification. In particular, the MC simulations confirm the crucial role of ionic valence-to-volume ratios in the competitive adsorption to charged surfaces that had been previously predicted by the mean-field theory. The charge inversion for ionic systems with salt is predicted by the MC simulations but not by the mean-field theory. This work provides a better understanding of competitive adsorption of counterions to charged surfaces and calls for further studies on the ionic size effect with application to large-scale biomolecular modeling.

Mean-field description of ionic size effects with nonuniform ionic sizes: a numerical approach

Ionic size effects are significant in many biological systems. Mean-field descriptions of such effects can be efficient but also challenging. When ionic sizes are different, explicit formulas in such descriptions are not available for the dependence of the ionic concentrations on the electrostatic potential, that is, there is no explicit Boltzmann-type distributions. This work begins with a variational formulation of the continuum electrostatics of an ionic solution with such nonuniform ionic sizes as well as multiple ionic valences. An augmented Lagrange multiplier method is then developed and implemented to numerically solve the underlying constrained optimization problem. The method is shown to be accurate and efficient, and is applied to ionic systems with nonuniform ionic sizes such as the sodium chloride solution. Extensive numerical tests demonstrate that the mean-field model and numerical method capture qualitatively some significant ionic size effects, particularly those for multivalent ionic solutions, such as the stratification of multivalent counterions near a charged surface. The ionic valence-to-volume ratio is found to be the key physical parameter in the stratification of concentrations. All these are not well described by the classical Poisson-Boltzmann theory, or the generalized Poisson-Boltzmann theory that treats uniform ionic sizes. Finally, various issues such as the close packing, limitation of the continuum model, and generalization of this work to molecular solvation are discussed.

The weighted correlation approach for density functional theory: a study on the structure of the electric double layer

Within the framework of density functional theory, a weighted correlation approach is developed in order to obtain the density distributions of an inhomogeneous fluid. It results in a formally exact expression, by means of the concept of a weighted pair correlation function, used to evaluate the change of the single-particle direct correlation function of the system relative to that of a reference state. When applying the approach for practical use, however, an approximation of the pair correlation function has to be made, along with an appropriate definition of a weight function. Noticeably, combining this approach with fundamental measure theory gives rise to a new method, which we call the FMT/WCA-k(2) approach, for studying the structural and thermodynamic properties of a charged hard-sphere fluid subjected to a spatially varying external potential. Application of the FMT/WCA-k(2) approach in a range of electrolyte concentrations and surface charge densities, against the Monte Carlo simulations, shows that it is superior to the typical approaches of density functional theory in predicting the ionic density profiles of both counter-ions and co-ions near a highly charged surface. It is capable of capturing the fine features of the structural properties of the electric double layers, to well reproduce the layering effect and the charge inversion phenomenon, also in strongly coupled cases where divalent counter-ions are involved. In addition, it is found that the FMT/WCA-k(2) approach even has an advantage over the anisotropic, hyper-netted chain approaches in giving better agreement with the Monte Carlo results.

Effects of counterion size on the attraction between similarly charged surfaces

Interaction between similarly charged surfaces can be attractive at high electrostatic coupling constants Ξ = l(B)Z(2)/μ(GC), where l(B) is the Bjerrum length, μ(GC) the Gouy-Chapman length, and Z the valency of counterions. While this effect has been studied previously in detail, as a function of surface charge density and valency of the pointlike counterions, much less is known about the effect of counterion size. We apply the Wang-Landau sampling Monte Carlo (MC) simulation method to compute the free energy F as a function of the scaled distance between the plates D̃=D/μ(GC) for a range of Ξ and scaled counterion radii R̃=R/μ(GC). We find that for large Ξ and small ion radius, there is a global equilibrium distance D̃=D̃(eq)=2(1+R̃), correctly giving the expected value at the point counterion limit. With increasing R̃ the global minimum in F(D̃) changes to a metastable state and finally this minimum vanishes when R̃ reaches a critical value, which depends on Ξ. We present a state diagram indicating approximate boundaries between these three regimes. The Wang-Landau MC method, as it is applied here, offers a possibility to study a wide spectrum of extended problems, which cannot be treated by the use of contact value theorem.

Additional considerations about the role of ion size in charge reversal

The effect of the ion size on the charge reversal process is studied via canonical Monte Carlo simulation. To this end, a primitive model of electrolyte is used to analyze the electric double layer formed by an asymmetric electrolyte in the presence of a charged planar wall. Different values of ion diameters and surface charge densities are used so as to determine the conditions at which the charge reversal first occurs. For each case, the apparent surface charge density is calculated as a function of the distance from the charged wall for the different electrolyte concentrations in order to establish the minimal salt concentration required for the charge reversal. We will refer to this electrolyte concentration as the reversal concentration and will show how it depends on the surface charge density and on the ion size. From the apparent surface charge density profiles, the distance from the wall at which the charge reversal arises as well as its intensity can be also inferred.

Effect of surface charge on colloidal charge reversal

The objective of this research work is to understand the effect of the surface charge density on the charge reversal phenomenon. To this end, we use experimental results and computer simulations. In particular, we measure the electrophoretic mobility of latex particles (macroions) in the presence of a multivalent electrolyte. We have focused on the electrolyte concentration range at which a reversal in the electrophoretic mobility is expected to happen. In particular, the role of the surface charge on the charge reversal process is looked into from several latexes with the same functional group but different surface charge densities. Although the mechanism responsible for the colloidal charge reversal is still a controversial issue, it is proved that ionic correlations are behind the appearance of such phenomenon (especially near the macroion surface). This conclusion can be inferred from a great variety of theoretical models. According to them, one of the factors that determine the charge reversal is the surface charge density of the macroions. However, this feature has been rarely analyzed in experiments. Our results appear therefore as a demanded survey to test the validity of the theoretical predictions. Moreover, we have also performed Monte Carlo simulations that take the ion size into account. The correlation found between experiments and simulations is fairly good. The combination of these techniques provides new insight into the colloidal charge reversal phenomena showing the effect of surface charge.

Effect of ion dispersion forces on the electric double layer of colloids: a Monte Carlo simulation study

In this work, the effect of ionic dispersion forces on the electric double layer of colloids is evaluated through Monte Carlo simulations. Particularly, the influence of these forces on the zeta-potential (as a representative electrokinetic property) is assessed. Ion polarizability is included in the primitive model with the help of the Lifshitz theory. In this way, ion specificity is not considered by means of phenomenological (and unknown a priori) parameters. Our results reveal that the ionic van der Waals forces are responsible (to some extent) for the specificity of the zeta-potential. In any case, the specific ion effects due to ion polarizability are strongly influenced by ion size. Furthermore, a preliminary study on the effect of ionic dehydration shows how this phenomenon improves the qualitative agreement between experimental data and simulations achieved in considering ionic dispersion forces.

Generalized Poisson-Fermi formalism for investigating size correlation effects with multiple ions

We establish a generalized Poisson-Fermi formalism to compute the electrostatic potential next to charged surfaces in the presence of multiple ion species with different sizes. A generalized Fermi-like ion distribution is deduced from the excess free energy, after expansion of the functional entropy of free space in which the ions have all the same size. The ion distribution is expressed in terms of the bulk volume fractions of each ion species rather than their bulk concentrations so as to account for the excluded volumes. We present size correlations effects such as underscreening and ion stratification, which have not been investigated before with such a simple theory. The change of dielectric properties across the space, arising from the finite spatial occupancy of ions, can be solved self-consistently through the Bruggeman model. The generalized Poisson-Fermi formalism is anticipated to be useful for interpreting electrophoretic mobility measurements and for computing the electrostatic potential over the surface of biomolecules in ionic solutions.

Surface charge density and electrokinetic potential of highly charged minerals: experiments and Monte Carlo simulations on calcium silicate hydrate

In this paper, we are concerned with the charging and electrokinetic behavior of colloidal particles exhibiting a high surface charge in the alkaline pH range. For such particles, a theoretical approach has been developed in the framework of the primitive model. The charging and electrokinetic behavior of the particles are determined by the use of a Monte Carlo simulation in a grand canonical ensemble and compared with those obtained through the mean field theory. One of the most common colloidal particles has been chosen to test our theoretical approach. That is calcium silicate hydrate (C-S-H) which is the main component of hydrated cement and is known for being responsible for cement cohesion partly due to its unusually high surface charge density. Various experimental techniques have been used to determine its surface charge and electrokinetic potential. The experimental and simulated results are in excellent agreement over a wide range of electrostatic coupling, from a weakly charged surface in contact with a reservoir containing monovalent ions to a highly charged one in contact with a reservoir with divalent ions. The electrophoretic measurements show a charge reversal of the C-S-H particles at high pH and/or high calcium concentration in excellent agreement with simulation predictions. Finally, both simulation and experimental results clearly demonstrate that the mean field theory fails not only quantitatively but also qualitatively to describe a C-S-H dispersion under realistic conditions.

Effect of the surface charge discretization on electric double layers: A Monte Carlo simulation study

The structure of the electric double layer in contact with discrete and continuously charged planar surfaces is studied within the framework of the primitive model through Monte Carlo simulations. Three different discretization models are considered together with the case of uniform distribution. The effect of discreteness is analyzed in terms of charge density profiles. For point surface groups, a complete equivalence with the situation of uniformly distributed charge is found if profiles are exclusively analyzed as a function of the distance to the charged surface. However, some differences are observed moving parallel to the surface. Significant discrepancies with approaches that do not account for discreteness are reported if charge sites of finite size placed on the surface are considered.

Monte Carlo simulation of electrolytes in the constant voltage ensemble

The authors studied the structural, electrostatic, and electromechanical properties of the terlamellar structure composed of the anode, the cathode, and the electrolyte layer separating them. They used the Monte Carlo simulation technique in the constant voltage ensemble, where the electrical potential difference between the anode and the cathode is introduced as an external field. For ions, they used the primitive models of different sizes and valences in order to investigate how they affect the physical properties when an electrical field is applied between the electrodes. For electrodes, they used impermeable and permeable models, which mimic planar and porous electrodes, respectively. The asymmetry between the anions and the cations in size or valence was found to be responsible for the asymmetry in the concentration profile, the potential drop, and the stress distribution, in comparing the anode and the cathode sides. The charging/discharging process in the planar and porous electrodes is discussed at molecular level.

Testing one component plasma models on colloidal overcharging phenomena

In this paper, the mechanisms of overcharging of a colloidal macroion in the presence of multivalent counterions are investigated by means of Monte Carlo simulations. This computational technique appears as a powerful tool for probing the validity of semianalytical models developed for this issue. In particular, the simulations performed are compared with the predictions of two different models based on the one component plasma (OCP) theory. Therein, the multivalent ionic atmosphere confined at the macroion surface is approximated by a two-dimensional Wigner crystal. These kinds of models are largely used in the literature since (in some cases) they present quite simple equations to describe the electric double layer (EDL) of macroions with different geometries in the presence of much smaller (but still multivalent) ions. In this sense, charge inversion phenomena of membranes, polyelectrolytes, DNA molecules, etc., are straightforwardly predicted in terms of these expressions. Unfortunately, comparisons between these predictions and experimental results are scarce, mostly due to the difficulty to reproduce the experimental conditions in the laboratory. Accordingly, the goal of the present paper is to simulate EDLs under real conditions (in which overcharging phenomena are expected to happen) and use the results obtained in this way for comparing with those obtained from OCP models.

Planar electric double layer for a restricted primitive model electrolyte at low temperatures

Monte Carlo simulation and the modified Poisson-Boltzmann theory are used to investigate the planar electric double layer for a restricted primitive model electrolyte at low temperatures. Capacitance as a function of temperature at low surface charge is determined for 1:1, 2:2, 2:1, and 3:1 electrolytes. Negative adsorption can occur for 1:1 electrolytes at low surface charge with low electrolyte concentration. The 1:1 electrolyte diffuse layer potential as a function of surface charge displays a maximum at low densities. At high densities, the diffuse layer potential is negative with a negative slope. The Gouy-Chapman-Stern theory fails in this low-temperature regime, whereas the modified Poisson-Boltzmann theory is fairly successful in this regard.

Further Test of Third Order + Second-Order Perturbation DFT Approach: Hard Core Repulsive Yukawa Fluid Subjected to Diverse External Fields

Grand canonical Monte Carlo simulation is used to investigate density profiles of hard-core repulsive Yukawa (HCRY) model fluid under the influence of various external fields and radial distribution function (RDF) of the bulk HCRY system. The aim of these extensive simulations is to provide exact data for purely repulsive interaction potential against which the validity of a third order + second-order perturbation DFT approach can be tested. It is found that a semiempirical parametrized bridge function due to Malijevsky and Labik performs very well for the RDF of the bulk HCRY fluid. Incorporation of a bulk second-order direct correlation function (DCF) of the HCRY fluid based on the Malijevsky-Labik bridge function into the third order + second-order perturbation DFT approach yields the resulting theoretical predictions for the density profiles of inhomogeneous HCRY fluid that are in a very good agreement with the simulation data, an exception being somewhat larger deviations appearing for the structure of the fluid around the center of a hard spherical cavity. Both theory and simulation predict layering transition and gas-liquid coexistence phenomena occurring with the HCRY model fluid under confined conditions. For the case of an inverse sixth-power repulsive potential under the influence of a flat stationary wall defined by an inverse twelfth-power repulsive potential, the present third order + second-order perturbation DFT approach is found to be superior to several existing weighted density approximations (WDA) and partitioned WDA.

Electric Double Layers with Electrolyte Mixtures: Integral Equations Theories and Simulations

A study of a planar electric double layer (EDL) in the presence of mixtures of electrolyte is presented. In particular, results from the Hyper-Netted-Chain/Mean-Spherical-Approximation (HNC/MSA) theory are compared with Monte Carlo (MC) simulations. In this way, the charge inversion induced by mixtures of multivalent and monovalent counterions is probed. Since overcharging phenomena in nature emerge under such conditions, the role of ion-ion correlations in the EDL appears as a crucial point in this kind of study. Unlike previous related works, a realistic hydrated ion size is used in the HNC/MSA calculations and simulations. In this way, a qualitative agreement between the results obtained from the theory and MC simulations is found. However, some discrepancies arise when the charge inversion is expected to be more noticeable, namely at high surface charges and/or elevated concentrations of multivalent electrolytes. Such differences are explained in terms of an overestimation of the charge inversion by the integral equation (IE) formalism.

Simulation of electric double layers undergoing charge inversion: mixtures of mono- and multivalent ions

In this paper, the electric double layer (EDL) of a charged plane in the presence of mixtures of 1:1 and 3:1 electrolytes has been investigated through Monte Carlo (MC) simulations using a nonrestrictive primitive model of EDL. In particular, the charge inversion in colloids (attributable to an accumulation of counterions on the surface) can be better understood by means of the simulations performed here. Moreover, two mechanisms proposed for charge inversion are probed: The formation of a strongly correlated layer (SCL) of multivalent counterions and excluded volume effects (to which we will also refer as ion size correlations). Our results are in agreement with the behavior found experimentally for some model colloids with increasing the concentration of monovalent salt in the presence of trivalent ions, which clearly supports the relevance of ion size correlations. In contrast, certain disagreement with predictions of SCL theories is reported.

Effect of counterion species on colloidal crystal

The effect of counterion species on the colloidal crystal structure in a dispersion was carefully investigated as a function of the degree of neutralization (alpha) by the ultra-small-angle X-ray scattering technique. The nearest neighbor interparticle distance (2D(exp)) first increased with decreasing alpha, and then decreased after passing through the maximum. This behavior was confirmed for K(+), Li(+), Ca(2+), TMA(+) (tetramethylammonium) as a counterion, and Na(+) in our previous report (Harada, T.; Matsuoka, H.; Ikeda, T.; Yamaoka, H. Langmuir 2000, 16, 1612). However, the alpha value of the maximum position (alpha(max)) largely depended on the counterion species, and it was in the order K(+) < Na(+) < TMA(+) approximately Li(+). This behavior was well characterized by the specific features of each ion: the alpha(max) map could be well superimposed in the Stokes radius-crystal ion radius relationship of counterions. The alpha(max) dependence on Stokes radius was very similar to that of the B coefficient by Jones and Dole except in the case of Ca(2+). In principle, the smaller the value for B, the smaller alpha(max), indicating that a water structure breaker such as K(+) can more easily destroy the colloidal crystal structure. In other words, the effect of the counterion species on colloidal crystal stability follows the Hofmeister series. Including Ca(2+), the relationship was linear for the alpha(max) values plotted as a function of the limiting equivalent conductivity of small ions. A counterion with larger conductivity would be a stronger breaker for the colloidal crystal structure.