Effect of Discrete Macroion Charge Distributions on Electric Double Layer of a Spherical Macroion

The effect of replacing the conventional uniform macroion surface charge density with discrete macroion charge distributions on the structure of electric double layer (EDL) of a spherical macroion has been investigated by Monte Carlo (MC) simulations. Two discrete models have been investigated in addition to the central macroion charge: point charges localized on the macroion surface and finite-sized charges protruding into the solution. Both models have been studied with fixed and mobile macroion charges. The radial functions of local densities and electrostatic potential in EDL, are calculated and compared to the results obtained for the central macroion charge distribution. It is concluded that the model of charge distribution significantly affects the EDL structure close to the macroion, while the effect is much weaker at larger distances. With point charges localized on the macroion surface, counterions become stronger accumulated to the macroion, as a result the absolute values of surface potential ?₀ and zeta ξ potential are decreased. With protruding charges, the excluded volume effect dominates over the increased correlation ability; hence the counterions are less accumulated near the macroions and the absolute values of ?₀ and ξ potentials are increased.     

Effect of discrete macroion charge distributions in solutions of like-charged macroions

The effect of replacing the conventional uniform macroion surface charge density with discrete macroion charge distributions on structural properties of aqueous solutions of like-charged macroions has been investigated by Monte Carlo simulations. Two discrete charge distributions have been considered: point charges localized on the macroion surface and finite-sized charges protruding into the solution. Both discrete charge distributions have been examined with fixed and mobile macroion charges. Different boundary conditions have been applied to examine various properties. With point charges localized on the macroion surface, counterions become stronger accumulated to the macroion and the effect increases with counterion valence. As a consequence, with mono- and divalent counterions the potential of mean force between two macroions becomes less repulsive and with trivalent counterions more attractive. With protruding charges, the excluded volume effect dominates over the increased correlation ability; hence the counterions are less accumulated near the macroions and the potential of mean force between two macroions becomes more repulsive/less attractive.     

Surface properties of fluids of charged platelike colloids

Surface properties of mixtures of charged platelike colloids and salt in contact with a charged planar wall are studied within density functional theory. The particles are modeled by hard cuboids with their edges constrained to be parallel to the Cartesian axes corresponding to the Zwanzig model [J. Chem. Phys. 39, 1714 (1963)] and the charges of the particles are concentrated at their centers. The density functional applied is an extension of a recently introduced functional for charged platelike colloids. It provides a qualitative approach because it does not determine the relation between the actual and the effective charges entering into the model. Technically motivated approximations, such as using the Zwanzig model, are expected not to influence the results qualitatively. Analytically and numerically calculated bulk and surface phase diagrams exhibit first-order wetting for sufficiently small macroion charges and isotropic bulk order as well as first-order drying for sufficiently large macroion charges and nematic bulk order. The asymptotic wetting and drying behaviors are investigated by means of effective interface potentials which turn out to be asymptotically the same as for a suitable neutral system governed by isotropic nonretarded dispersion forces. Wetting and drying points as well as predrying lines and the corresponding critical points have been located numerically. A crossover from monotonic to nonmonotonic electrostatic potential profiles upon varying the surface charge density has been observed. Nonmonotonic electrostatic potential profiles are equivalent to the occurrence of charge inversion. Due to the presence of both the Coulomb interactions and the hard-core repulsions, the surface potential and the surface charge do not vanish simultaneously, i.e., the point of zero charge and the isoelectric point of the surface do not coincide.     

Effects of image charges, interfacial charge discreteness, and surface roughness on the zeta potential of spherical electric double layers

We investigate the effects of image charges, interfacial charge discreteness, and surface roughness on spherical electric double layer structures in electrolyte solutions with divalent counterions in the setting of the primitive model. By using Monte Carlo simulations and the image charge method, the zeta potential profile and the integrated charge distribution function are computed for varying surface charge strengths and salt concentrations. Systematic comparisons were carried out between three distinct models for interfacial charges: (1) SURF1 with uniform surface charges, (2) SURF2 with discrete point charges on the interface, and (3) SURF3 with discrete interfacial charges and finite excluded volume. By comparing the integrated charge distribution function and the zeta potential profile, we argue that the potential at the distance of one ion diameter from the macroion surface is a suitable location to define the zeta potential. In SURF2 model, we find that image charge effects strongly enhance charge inversion for monovalent interfacial charges, and strongly suppress charge inversion for multivalent interfacial charges. For SURF3, the image charge effect becomes much smaller. Finally, with image charges in action, we find that excluded volumes (in SURF3) suppress charge inversion for monovalent interfacial charges and enhance charge inversion for multivalent interfacial charges. Overall, our results demonstrate that all these aspects, i.e., image charges, interfacial charge discreteness, their excluding volumes, have significant impacts on zeta potentials of electric double layers.     

Simulation of the electrical double layer of a macroion with different counterion charges

An electrical double layer of a spherical macroion with single-, double-, and triple-charged counterions in aqueous solution of 1: 1 background electrolyte at different concentrations are studied by the molecular dynamics method for models with discrete and continuous surface charge distribution. Radial profiles of ion partial densities and the electric potential distribution in the double layer are calculated. The degree of counterion binding with a macroion is determined. The effect of water permittivity on the structure of electrical double layer is studied.     

Discrete charge effects on the Donnan potential and surface potential of a soft particle

The Donnan potential and surface potential of soft particles (i.e., polyelectrolyte-coated hard particles) in an electrolyte solution play an essential role in their electric behaviors. These potentials are usually derived via a continuum model in which fixed charges inside the surface layer are distributed with a continuous charge density. In this paper, for a plate-like soft particle consisting of a cubic lattice of fixed point charges, on the basis of the linearized Poisson–Boltzmann equation, we derive expressions for the electric potential distribution in the regions inside and outside the surface layer. This expression is given in terms of a sum of the screened Coulomb potentials produced by the point charges within the surface layer. We show that the deviation of the results of the discrete charge model from those of the continuous charge model becomes significant as the ratio of the lattice spacing to the Debye length becomes large.     

A macroion electrokinetics algorithm

A numerical algorithm is presented for the standard model of macroion electrokinetics and certain generalizations of it. The macroion consists of a cylindrical section with identical, hemispheroidal endcaps, each piece having arbitrary length. The system of one macroion and adjoining salt solution is subjected to an arbitrary sequence of pulsed electrical fields and pulsed translational and rotational velocities. Numerical solutions are obtained for the time dependent electrostatic and mobile ion concentration fields and the solvent velocity. From these fields the dielectric response, force, and torque are calculated. Generalizations of the standard model include the diffusive motion of macroion surface charges, partial slip of solvent motion at the macroion surface, and a simple model for the reactive exchange of surface charge with solution ions. The primary illustrative application is to recent measurements of electric birefringence versus applied field frequency for poly-(tetrafluorothylene) colloidal particles, but a few results are presented for the dielectric response of DNA fragments and of spherical colloidal particles. The source code and additional details are provided as supplementary documentation.     

Investigation of the effect of finite-sized ions on the nanowire field-effect transistor in electrolyte concentration using a modified Poisson–Boltzmann model

One-dimensional (1D) nanowire field-effect transistors (FETs) have recently played a major role in sensing applications. Due to charging of the surface functional chemical groups with protonation and deprotonation, the transport properties of these nanowire transistors affect the aqueous environment, altering the electrical double layer (EDL) potential drops and charge distributions in the electrolyte concentration. In this work, we have implemented the simple modified Poisson–Boltzmann (MPB) theory in a 1D silicon nanowire FET, and the effect of the various finite sizes of ions in z:z symmetric electrolyte concentration was investigated. For a given ionic concentration and surface charge, the EDL potential drop, accumulation of charges and the charge distributions of NaCl and CsCl ions were studied. From the MPB model results with the nanowire FET, it was observed that the potential drop of the EDL depends on the size of the ions in the electrolyte. The study of various electrostatic investigations of finite-sized ions was successfully done by implementing the MPB model on a silicon nanowire FET. It can be used in both chemical and biological sensors.     

Free isotropic-nematic interfaces in fluids of charged platelike colloids

Bulk properties and free interfaces of mixtures of charged platelike colloids and salt are studied within the density-functional theory. The particles are modeled by hard cuboids with their edges constrained to be parallel to the Cartesian axes corresponding to the Zwanzig model. The charges of the particles are concentrated in their center. The density functional is derived by functional integration of an extension of the Debye-Hückel pair distribution function with respect to the interaction potential. For sufficiently small macroion charges, the bulk phase diagrams exhibit one isotropic and one nematic phase separated by a first-order phase transition. With increasing platelet charge, the isotropic and nematic binodals are shifted to higher densities. The Donnan potential between the coexisting isotropic and nematic phases is inferred from bulk structure calculations. Nonmonotonic density and nematic order parameter profiles are found at a free interface interpolating between the coexisting isotropic and nematic bulk phases. Moreover, electrically charged layers form at the free interface leading to monotonically varying electrostatic potential profiles. Both the widths of the free interfaces and the bulk correlation lengths are approximately given by the Debye length. For fixed salt density, the interfacial tension decreases upon increasing the macroion charge.     

Influence of Fixed Charge Distribution on Ion Mobility: Small Deviations of Equilibrium Potential from Donnan Value

The influence of the electrostatic interaction on the diffusion coefficients of ions is estimated in the approximation of the slightly varying equilibrium potential. Various assumptions have been considered for the fixed charge distribution, namely, the discrete charge distributed in the bulk, the surface charge, and the linear charge distributed around a macroion. Both the random and regular distributions of the fixed charge groups are being separately analyzed as applied to each assumption. Obtained dependencies show in some measure various types of the fixed charge heterogeneity contribute to the effect of the diffusion coefficient decrease.     

Effects of structural properties of the Stern layer on the electrophoretic migration of a highly charged spherical macroion

The electrophoretic migration of a highly charged spherical macroion suspended in an aqueous solution of NaCl is studied using the molecular dynamic method. The objective is to examine the effects of the colloidal surface charge density on the electrophoretic mobility (μ) of the spherical macroion. The bare charge and the size of the macroion are varied separately to induce changes in the colloidal surface charge density. Our results indicate that μ depends on colloidal surface charge density in a nonmonotonic manner, but that this relationship is independent of the way the surface charge density is varied. It is found that an increase in colloidal surface charge density may lead to the formation of new sublayers in the Stern layer. The μ profile is also found to have a local maximum for a bare charge at which a new sublayer is formed in the Stern layer, and a local minimum for a bare charge at which the outer sublayer becomes relatively dense. Finally, the electrophoretic flow caused by the migration of the spherical macroion is studied to find that one decisive factor causing the electrophoretic flow is the ability of the macroion to carry anions in the electrolyte solution.     

The Electrical Double Layer at the Bi–Ga Liquid Alloy in the Acetonitrile Electrolyte Solutions

The electrical double layer (EDL) characteristics of a Bi (0.25 at. %)–Ga liquid electrode in acetonitrile are studied. In acetonitrile, as in water, the electrode models Bi electrodes in their electrochemical properties. In contrast to aqueous solutions, in acetonitrile solutions it is possible to study the EDL structure not only at negative charges but also near the zero charge and at small positive charges. In acetonitrile, the electrode's potential of zero charge not distorted by specific ion adsorption and the corrected electrochemical work function are determined. The electrode is as lyophilic with respect to acetonitrile as Hg. Thus, the orientation of acetonitrile dipoles on Hg and Bi is identical at negative charges. The obtained data and data for a partially fused polycrystalline electrode in aqueous solutions are used to determine contributions of semimetal properties of a Bi electrode to the capacitance of the inner part of EDL and to a potential drop. The charge dependences of these contributions are found.     

Simulation of the electrical double layer of a spherical micelle of an anionic substance with regard to the solvent structure

Computer simulation methods are employed to consider the structure of the electrical double layer of a spherical micelle in aqueous surfactant solutions with allowance for the contribution of the solvent. Three micelle models were used in the calculations, namely, a macroion with discretely distributed charges and a continual solvent, a spherical model micelle with a coarse-grained representation of the solvent, and a spherical model micelle in an aqueous phase with an explicit account for water molecules. Based on these three models, the radial profiles of the local densities and electric potentials in the electrical double layer, as well as the degrees of binding single-, double-, and triple-charged counterions by the macroion in aqueous surfactant solutions, are calculated with regard to the Lennard-Jones and electrostatic interactions. The allowance for the molecular structure of the solvent leads to qualitatively different local dependences of the electric potential as compared to both the continual and coarse-grained representation of the solvent.     

Effect of Anion-Induced Cation Adsorption on the Electrical Double Layer Structure

Effect of supporting cations on the EDL structure is analyzed in the region of specific adsorption of anions. The analysis assumes that cations are pulled into the EDL dense part due to their interaction with specifically-adsorbed anions. The model describing this phenomenon is based on the introduction of a second, cationic IHP situated between the plane where centers of charges of specifically-adsorbed anions are localized and OHP. The adsorption of ions, electrode charge, PZC, and EDL capacitance are calculated using two adsorption isotherms for anions and cations.     

Electrophoretic properties of highly charged colloids: a hybrid molecular dynamics/lattice Boltzmann simulation study

Using computer simulations, the electrophoretic motion of a positively charged colloid (macroion) in an electrolyte solution is studied in the framework of the primitive model. In this model, the electrolyte is considered as a system of negatively and positively charged microions (counterions and coions, respectively) that are immersed into a structureless medium. Hydrodynamic interactions are fully taken into account by applying a hybrid simulation scheme, where the charged ions (i.e., macroion and electrolyte), propagated via molecular dynamics, are coupled to a lattice Boltzmann (LB) fluid. In a recent electrophoretic experiment by Martin-Molina et al. [J. Phys. Chem. B 106, 6881 (2002)], it was shown that, for multivalent salt ions, the mobility mu initially increases with charge density sigma, reaches a maximum, and then decreases with further increase of sigma. The aim of the present work is to elucidate the behavior of mu at high values of sigma. Even for the case of monovalent microions, a decrease of mu with sigma is found. A dynamic Stern layer is defined that includes all the counterions that move with the macroion while subjected to an external electrical field. The number of counterions in the Stern layer, q(0), is a crucial parameter for the behavior of mu at high values of sigma. In this case, the mobility mu depends primarily on the ratio q(0)/Q (with Q the valency of the macroion). The previous contention that the increase in the distortion of the electric double layer (EDL) with increasing sigma leads to the lowering of mu does not hold for high sigma. In fact, it is shown that the deformation of the EDL decreases with the increase of sigma. The role of hydrodynamic interactions is inferred from direct comparisons to Langevin simulations where the coupling to the LB fluid is switched off. Moreover, systems with divalent counterions are considered. In this case, at high values of sigma the phenomenon of charge inversion is found.     

Structure of spherical electric double layers: a density functional approach

A density functional theory is presented for the structure of spherical electric double layers within the restricted primitive model, where the macroion is considered as a hard sphere having uniform surface charge density, the small ions as charged hard spheres, and the solvent is taken as a dielectric continuum. The theory is partially perturbative as the hard-sphere contribution to the one-particle correlation function is evaluated using suitably averaged weighted density and the ionic part is obtained through a second-order functional Taylor expansion around the uniform fluid. The theory is in quantitative agreement with Monte Carlo simulation for the density profiles and the zeta potentials over a wide range of macroion sizes and electrolyte concentrations. The theory is able to provide interesting insights about the layering and the charge inversion phenomena occurring at the interface.     

A Monte Carlo study of spherical electrical double layer of macroions-polyelectrolytes systems in salt free solutions

A canonical Monte Carlo simulation is performed to investigate the microstructure and the electrical double layer (EDL) of polyelectrolytes around macroions in the bulk systems based on the primitive model. We explore the influences of particles size, chain length, and charge density of polyelectrolytes on the microscopic behavior of the macroions-polyelectrolytes systems. The simulation results show that the surface charge density and the chain length of the polyelectrolytes are two key factors that affect the microstructure of polyelectrolytes around the macroions and potential of mean force between the macroions as well as the zeta potential of the spherical EDL constructed by polyelectrolytes. The high surface charge density of a polyelectrolyte leads to the polyelectrolyte acting as a bridge for the aggregation of macroions, causing the presence of the attraction between macroions. The polyelectrolytes with a long chain length present a cooperativity effect for the adsorption of the polyelectrolytes on the surface of the macroions. Furthermore, the two key factors both induce the overcharge of the macroions. The longer the chain length and the higher surface charge density of the polyelectrolytes, the stronger is the overcharge.     

The influence of discrete surface charges on the force between charged surfaces

The force between two parallel charged flat surfaces, with discrete surface charges, has been calculated with Monte Carlo simulations for different values of the electrostatic coupling. For low electrostatic coupling (small counterion valence, small surface charge, high dielectric constant, and high temperature) the total force is dominated by the entropic contribution and can be described by mean field theory, independent of the character of the surface charges. For moderate electrostatic coupling, counterion correlation effects lead to a smaller repulsion than predicted by mean field theory. This correlation effect is strengthened by discrete surface charges and the repulsive force is further reduced. For large electrostatic coupling the total force for smeared out surface charges is known to be attractive due to counterion correlations. If discrete surface charges are considered the attractive force is weakened and can even be turned into a repulsive force. This is due to the counterions being strongly correlated to the discrete surface charges forming effective, oppositely directed, dipoles on the two walls.