Chandra N. Patra*ab
aTheoretical Chemistry Section, Chemistry Group, Bhabha Atomic Research Centre, Mumbai 400 085, India. E-mail: chandra@barc.gov.in
bHomi Bhabha National Institute, Mumbai, India
First published on 18th February 2015
A systematic study on the structure of mixed electrolytes with arbitrary size and charge asymmetry around a charged nanoparticle is carried out using density functional theory and Monte Carlo simulation. A primitive model representation, with the spherical macroion surrounded by small ions in a continuum dielectric as the solvent, is used. A weighted density approximation is used to evaluate the hard-sphere correlation, whereas the ionic part is calculated using perturbation expansion around the bulk density. A canonical ensemble Monte Carlo simulation on the same system is also performed for comparison. Parametric variations on component ratios of the electrolyte, ionic concentrations, surface charge densities and macroion as well as small ion sizes show interesting phenomena of overcharging and charge reversal. The theoretical predictions are found to be in good agreement with the simulation results concerning the density as well as mean electrostatic potential profiles. The present study shows distinctive evidence of size and charge correlations around the interface, in a fully asymmetric situation.
Macromolecular solutions denote another class of substance where interfacial interactions play a vital role, and the quality and quantity of the same have been used to predict a number of interfacial phenomena encompassing electrochemistry, biology, materials science, etc.7 Interfacial interactions also provide a detailed understanding for the development of properties of technologically important polyelectrolytes used in agricultural engineering and food science, paints, varnishes, etc.8 This particular subject has found increasing attention over the last decade due to the feasibility of measurements of subtle length and time scales.9 Thus, optical tweezer measurements10 in colloidal solutions provide direct evidence of attraction in like-charged systems. Similarly, confocal microscopy11 in bilayers and biological systems gives direct imaging of the topology of layers and their interactions.
A colloidal solution, along with supporting electrolytes, leads to the formation of electric double layers (EDLs),12,13 where the first layer (the so-called Helmholtz layer) remains within the size of the counterion and the diffuse layer extends beyond that, depending on the electrostatic interactions between the component ions and also with the interface. The geometry of the system merely determines the geometry of the respective double layer, i.e., a spherical double layer (SDL) for the case of a colloidal solution.14 However, when the colloidal macroion becomes sufficiently large, the symmetry of the system tends to a planar one, resulting in planar double layers.15 Similarly, biomacromolecular solutions, as well as DNA, represent cylindrical double layers.16 EDLs have been subjected to rigorous theoretical studies, mainly because of a wide variety of parametric variations leading to different types of double layers. Most of these studies are restricted to model systems as the applicability of the theoretical calculations can be compared to that of computer simulations. These studies also advance proper understanding of the suitability or refinement of different models used in double layer theories.
Theoretical study on EDLs started with the point-ion Poisson–Boltzmann (PB) description,17 however, the introduction of correlations has emerged only in the last decades, thereby, the modified Poisson–Boltzmann theory,18 integral equation theory (IET),19 density functional theory (DFT)20 and their different variants have taken their course in the last few years to study EDL in detail. All these studies not only improved fundamental understanding but also proposed suitable applications of the same in even more complex systems. Planar21–26 and cylindrical27–30 EDLs have been studied extensively because of their applicability to electrochemical cells, synthetic polyelectrolytes, and DNA, however, SDLs have attracted a great deal of attention with the advent of mechanisms of fabrication of nanomaterials through suitable modifications and the assembling of spherical nanoparticles. In most of these studies, the restricted primitive model (RPM)12 is used, where the electrode is a structureless large hard sphere with uniform surface charge density and small ions are charged hard spheres with equal diameter embedded in a solvent of dielectric continuum. This model is extended to include solvent as a component as the solvent primitive model (SPM).31 Including different sized small ions beyond the RPM has been worked out in great detail recently,32–36 although the PB solution for the system has been known over the years.37 Thus, IET involving the hypernetted chain/mean spherical approximation (HNC/MSA) and Monte Carlo (MC) simulations has been used quite extensively by Guerrero-García et al.34,36,38–40 to study the effect of size and charge asymmetry on SDLs. The other contributions that require a mention here are reference hypernetted chain approximations,32 reference fluid density functional approaches,33 and the exclusion volume corrected PB theory.35 The important inferences that counterions do not always determine the properties of the double layers and the presence of potential of zero charge in size-asymmetric SDLs as compared to its presence in certain cases in size-symmetric ones have been revealed in all these works.
Due to its ensuing success in handling classical correlations, DFT has been applied to size-asymmetric SDLs in recent times. Thus, Kim et al.41 applied perturbative DFT, with the hard-sphere component of the excess free energy from the Percus–Yevick equation and the ionic part through a perturbation expansion around a uniform fluid. Another perturbative DFT approach involving approximating the hard-sphere part through fundamental measure theory has also been applied42 to the same system within the RPM, SPM and solvent RPM. In both cases, the theoretical results seem to reproduce the MC data quite well. In a recent work,43 we carried out a comprehensive study using DFT and MC simulations of a fully asymmetric SDL, i.e., with both charge- and size-asymmetry of its components. The theoretical predictions compare quite quantitatively with that of the simulation results and provide interesting observations of overcharging (OC)44 and charge reversal (CR)28 phenomena.
Mixed salts containing ionic components having different valance and sizes are quite important45 in many macromolecular and biological phenomena,46 and hence have been the subject of recent research.42 The structure of a SDL containing mixed electrolytes within the RPM has been studied in detail47 and important observations of CR with the gradual addition of multivalent ions have been reported.48 Planar49,50 as well as cylindrical EDL51 systems containing mixed electrolytes within the RPM have also been studied. Since it is now quite well known that the presence of multivalent ions has a major influence on size-asymmetric SDLs, it would be worthwhile to work on such basic electrolyte systems containing arbitrary charge and size asymmetry. The theoretical formalism and MC simulation method remain the same as reported in the binary electrolyte case.43 Thus, the partially perturbative DFT used here considers a weighted density approach (WDA) for the hard-sphere part and the perturbation around the uniform fluid for the residual ionic component. The presence of multivalent ions in the charge and size asymmetric electrolyte is studied in terms of density as well as the mean electrostatic potential profiles of the SDL under different parametric conditions, viz. the concentration of the multivalent ions, bulk concentration of electrolyte, surface charge density, size of the macroion, and size of the small ions. In Section II, we present the model and formalisms, the results and discussion part is presented in Section III and the concluding remarks in Section IV.
(1) |
(2) |
Therefore, the macroion–ion interaction potential is written as
(3) |
Although the above equations are written in CGS form, SI units are actually used for different parameters as mentioned throughout the text. In the present work, the diameter ratio for the ions in the mixed electrolyte is always taken as 1:2:3, σ1 = σ, σ2 = 2σ, σ3 = 3σ, with the smallest ionic diameter (σ) as 0.2125 nm. The multivalent ion is always taken as the smallest one, which follows from the analysis of Rakitin and Pack52 about the use of bare ionic diameters without a hydration shell. It is worth noting that the diameter ratio chosen here is only a representative one and the theoretical formalism is equally applicable for any arbitrary diameter ratio. Also because of the presence of three different sized small ions in the PM electrolyte, there appears three Helmholtz planes, viz. the inner Helmholtz plane (IHP) and two outer Helmholtz planes (OHP) at distances R + σα/2, α = 1, 2, 3, respectively, as compared to only a single Helmholtz plane at R + σ/2 for size-symmetric (σ1 = σ2 = σ3 = σ) SDLs.
(4) |
ρα(r) = ρ0αexp{−β0zαψ(r) + c(1)hsα(r;[{ρα}]) − c(1)hsα([{ρ0α}]) + c(1)elα(r;[{ρα}]) − c(1)elα([{ρ0α}])}, | (5) |
(6) |
Although eqn (5) for the density profile is formally exact, it still requires approximations for c(1)hsα and c(1)elα for the nonuniform SDL system. In a recent work, we adopted a partially perturbative approach to calculate these two quantities for a SDL with arbitrary size and charge asymmetry. In this approach, the quantity c(1)hsα(x; [{ρα}]) is calculated using the weighted density approximation of Denton and Ashcroft (DA)54 as
c(1)hsα(r;[{ρα}]) = (1)hsα((α)(r)), | (7) |
(8) |
The direct correlation functions (2)hsαβ and (2)elαβ for the PM electrolyte are taken from the analytical expressions within the mean spherical approximation as given by Blum55 and Hiroike56 for arbitrary ionic size and charge asymmetry. Once the approximate first order correlations functions are obtained through eqn (7) and (8), the density and the mean electrostatic potential profiles are calculated from eqn (5) and (6) using a Picard iteration. The convergence criterion is chosen in such a way that after nth iteration, the norm defined as attains a small value (<10−6), where N is the total number of mesh points.
The increased charge inversion is also reflected in the mean electrostatic potential (MEP) profiles for systems of varying [Mg2+]:[Na+] concentrations. As depicted in Fig. 2, the contact value of the MEP at the first OHP continuously decreases as the concentration of divalent ions increases. The faster decay of the MEP near the interface is accompanied by an increase in the same a little away from the interface, indicating a change in sign. There is an increase in the depth of the inversion, although the width of the inversion layer remains more or less constant. The density and the MEP profiles predicted by the present DFT are in quantitative agreement with the MC data over the entire range of concentrations as seen by the variations in [Mg2+]:[Na+].
An interesting situation emerges as the charge on the macroion turns negative. Fig. 3 depicts the density profiles of counterions (Na+ and Mg2+) and coions (Cl−) as the ratio of [Mg2+]:[Na+] increases around a macroion of R = 1.5 nm with Q = −0.102 C m−2 (ZM = −18). The substantial increase in coion concentration in the second OHP leads to considerable OC phenomena as more and more divalent counterion is added. With a small amount of divalent salt (MgCl2) (Fig. 3(a)), the contact density of Mg2+ at the macroion surface is increased substantially accompanied by the substantial decrease of Na+. This is quite obvious because of the higher valency of Mg2+ as compared to Na+. However, when more Mg2+ is added, the contact density starts decreasing because of screening of the surface charge on the macroion. The is also reflected in the increase in CR effect with Mg2+, which shifted to a shorter distance. This is quite evident from the MEP profiles (Fig. 4), where the contact density at the IHP continuously increases with an increase in [Mg2+]:[Na+] ratio. In effect, size correlations due to volume exclusion have also become equally important as charge correlations in this system. Once again the theoretical prediction and MC data go hand in hand.
The effect of bulk ion–ion correlations of small ions is studied by increasing the bulk concentration of NaCl from 0.01 M to 2 M, while maintaining [Mg2+]:[Na+] = 1:2. As can be seen in Fig. 5, the volume exclusion becomes extremely important, leading to an increase in oscillations in all the ionic density profiles. The double layer becomes more compact and the diffuse layer is narrowed down due to an increase in ion–ion correlations. This leads to an effective screening of macroion charge resulting in a CR effect as is depicted in Fig. 6. The CR is observed twice because of synergy in size and charge correlation. The “periodicity” of oscillations in the density profiles relates to the diameter of the individual ionic components. The layering and the charge inversion becomes quite distinct, although these originate from the same physical phenomena. There is a substantial increase in the depth of the inversion layer, however its width decreases due to stronger electrostatic interactions. Quite surprisingly, the DFT predictions show quantitative matching with MC simulation. This is in sharp contrast to our earlier work on size-asymmetric electrolytes.43
Fig. 5 Ionic density profiles around a spherical macroion of R = 1.5 nm and Q = 0.102 C m−2 for NaCl/MgCl2 salt with [Mg2+]:[Na+] = 1:2, at various concentrations of: (a) 0.01 M, (b) 0.1 M, (c) 1 M, and (d) 2 M. The key is the same as in Fig. 1. |
To study the implications of variation in the surface charge density on the SDL formed from 1 M NaCl/MgCl2 with [Mg2+]:[Na+] = 1:2 around a macroion of R = 1.5 nm, Q is continuously varied from 0.102 to 0.408 C m−2 (corresponding variation in total charge on the nanoparticle from ZM = 18 to ZM = 72). This is depicted in Fig. 7, where there is considerable accumulation of counterions and depletion of coions from the interface as Q increases. The stronger electrostatic attraction between the counterion and the macroion also leads to effective screening of the macroion charge resulting in charge inversion. Once the inversion sets in, the width as well as the depth of the double layer remains constant. The same conclusions can be drawn from the MEP profiles as shown in Fig. 8.
Fig. 7 Ionic density profiles for 1 M NaCl/MgCl2 salt with [Mg2+]:[Na+] = 1:2, around a spherical macroion of R = 1.5 nm at varying surface charge densities: (a) Q = 0.102 C m−2, (b) Q = 0.204 C m−2, (c) Q = 0.306 C m−2, and (d) Q = 0.408 C m−2. The key is the same as in Fig. 1. |
The charge correlations in SDL are better understood in terms of the size of the macroion, along with its surface charge density. This is shown in Fig. 9 for a 1 M NaCl/MgCl2 mixed electrolyte with [Mg2+]:[Na+] = 1:2 at Q = 0.102 C m−2 for four different macroion radii, viz. R = 0.5 nm, 1 nm, 1.5 nm, and 6 nm. The total charge on the macroion for these four sizes corresponds to ZM = 2, 8, 18, and 288 respectively. As we move on to the larger macroion, there is an increase in absolute charge (ZM) on the macroion resulting in an enhancement of counterions and a consequent depletion of coions at the macroion surface. The diffuse layer widens with the faster decay of counterions, accompanied by a sharp growth of the coion density profiles at a distance away from the surface. This also results in an increase in the CR effect. As evident from the MEP profile in Fig. 10, the depth of the CR layer increases, whereas the width remains more or less constant. Once again, the agreement between the theory and simulation results is excellent.
Fig. 9 Ionic density profiles for 1 M NaCl/MgCl2 salt with [Mg2+]:[Na+] = 1:2, around a spherical macroion of Q = 0.102 C m−2 at different macroion radii: (a) R = 0.5 nm, (b) R = 1 nm, (c) R = 1.5 nm, and (d) R = 6 nm. The key is the same as in Fig. 1. |
In order to study the effect of valence on the behavior of the diffuse layer, the added electrolyte is changed to AlCl3 instead of MgCl2. Thus, Fig. 11 depicts the ionic density and the MEP profiles of a NaCl/AlCl3 mixed electrolyte with [Al3+]:[Na+] = 1:2 with 1 M NaCl. The coion and counterion densities remain almost the same, although there is a slight increase in Cl− and Al3+ at the interface (cf. Fig. 1(c)). This also results in a stronger charge inversion for the higher valence coion (Al3+) because of the excess of counterions effectively screening the surface charge. This is also clear from the MEP profile, where the inversion width essentially gets reduced (cf. Fig. 2).
The importance of size correlations is also seen with a change in ionic sizes (σ) of the small ions. Fig. 12 shows the ionic density profiles for 1 M NaCl with added MgCl2 for [Mg2+]:[Na+] = 1:2 at four different small ion diameters, viz. σ = 0.1 nm, 0.15 nm, 0.2 nm, and 0.25 nm. For the smallest σ = 0.1 nm, the counterion (Cl−) density profile decays monotonically at the expense of the growth of the coion (Na+ and Mg2+) density profiles. However, the diffuse layer gets compacted for σ = 0.15 nm because of stronger size correlations. As the ionic size is increased to σ = 0.2 nm, the size correlation converts to charge correlation leading to charge inversion at a distance away from the interface. At the largest ionic sizes studied here for σ = 0.25 nm, the double layer become more structured with the presence of another layer for counterions (Cl−) and monovalent coions (Na+). This is due to the volume exclusion resulting from the larger ionic sizes. In fact, this system, along with the concentration variation (cf. Fig. 5), shows clear evidence of a distinctive contribution from size and charge correlations. This step by step effect of correlations can also be visualized in the MEP profiles shown in Fig. 13, which also indicates double charge inversions for the case of the largest ionic size (σ = 0.25 nm). The increase of the amount of coions at all the Helmholtz planes is an indication of the increase of the OC effect with an increase in small ion sizes. The double layer also becomes thinner and deeper as is evident from both the density as well as the MEP profiles. The DFT predictions are in quantitative agreement with the MC profiles for all the ionic diameters studied in the present work.
Fig. 12 Ionic density profiles for 1 M NaCl/MgCl2 salt with [Mg2+]:[Na+] = 1:2, around a spherical macroion of R = 1.5 nm and Q = 0.102 C m−2 with small ion diameters of (a) σ = 0.1 nm, (b) σ = 0.15 nm, (c) σ = 0.2 nm, and (d) σ = 0.25 nm. The key is the same as in Fig. 1. |
The importance of OC and CR effects are clearly reflected with the addition of multivalent ions in PM electrolytes. The present work provides direct support of the important fact that coions become extremely important in a SDL formed from a PM electrolyte, unlike RPM, where counterions determine the very nature of the SDL formed. The interplay between charge and size correlations has important implications on the diffuse layer characteristics in terms of compactness as well as charge inversions. Thus, ionic correlations become important with an increase in [Mg2+]:[Na+] ratio, bulk concentration, the size of macroion, and the surface charge density on the macroion. Similarly, size correlations take precedence with an increase in the size of the small ions of the mixed electrolytes, as revealed in charge inversions and zeta potentials. Thus, the addition of divalent coions (Mg2+) to the 1:1 (NaCl) electrolyte system causes a decrease of the layering of counterions (Cl−). The formation of multiple layers at large bulk concentration and in the case of the largest small ion is indicative of the distinctive contributions of the charge and size correlations in a SDL formed from a PM electrolyte solution. The effective screening of macroion surface charge also leads to important conclusions regarding depth and width of inversion layers.
The present primitive model used to study the mixed electrolyte system is a completely idealized situation, since at these sizes, non-electrostatic effects such as ionic polarizability will have sufficient contributions.59 The interaction between hydration shells around individual nanoparticles as well as that of other nanoparticles is also important. The actual atomistic structure of the nanoparticle itself will be another factor determining the SDL formed. A detailed insight into such situations may emerge from the use of more civilized models60 as well as considering explicit interactions from the solvent and ions.61
The present work can be equally important in cylindrical geometry, as that will provide important information concerning DNA conformation in different supporting electrolyte solutions,62 which gives a direct mechanism to control gene expression. The addition of solvent and its associated effects are another important area that require immediate attention.63 Other work should be to predict the electrokinetic charge of the colloidal macroion in electrophoresis since this is related to the MEP at the slipping plane.64 The present theory can be directly applicable to calculate the force between the large colloidal particles immersed in a mixed electrolyte solution and provide important implications about the nature of forces involved.65 At present, work along these directions is in progress and will be reported in the future.
This journal is © The Royal Society of Chemistry 2015 |