Theory of the effect of small ions on the electrical admittance of macroion solutions

A general theory is given on the electrical admittance of a parallel plate capacitor filled with a conductive medium containing charged particles by considering both their electric atmospheres and the electric field due to charges induced on the electrodes. The theory is applied to a dilute macroion solution with added salt and it is found that the low frequency dielectric response of the solution reflects the motions of the macroion and bound small ions and the effect of motion of the ion atmosphere around the macroion is canceled out by the current due to the induced charges on the electrodes.     

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.     

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.     

Computer simulation of macroion layering in a wedge film

The layering of macroions confined to a wedge slit formed by two uncharged hard walls is studied using a canonical Monte Carlo method combined with a simulation cell that contains both wedge-shaped slit and bath regions. The macroion solution is modeled within a one-component fluid approach that in an effective way incorporates the double layer repulsion due to simple electrolyte ions as well as the discrete nature of an aqueous solvent. The layer formation under a wedge confinement is analyzed by carrying out separate simulation runs for a set of consecutive wedge segments designed to represent a single wedge slit. As the wedge thickness progressively increases, the sequence of regions along the wedge film with distinct features of macroion layering has been established. This sequence comprises (i) a wedge region of the thickness smaller than the macroion diameter that is free of macroions; (ii) a region with a one-dimensional macroion chain along the wedge corner at a wedge thickness of a one macroion diameter; (iii) a region comprising a low-ordered macroion monolayer that extends until the wedge thickness slightly above two macroion diameters; (iv) a region comprising a pair of well-defined two-dimensional configurations of macroions segregated on each of the wedge walls; and (v) a free-of-macroions wedge region between two surface monolayers that now originates from an electrostatic repulsion imposed by the surface macroions, which is followed by (vi) a well-defined macroion monolayer film between two surface monolayers, a less defined bilayer film, a three-layer film, and so on up to the bulk solution. Once formed, the macroion surface monolayers persist for all remaining wedge thicknesses up to the bulk, forming in such a way effective charged wedge boundaries. Such a formation of the macroion surface monolayers on the uncharged confining walls is related to the haloing mechanism for regulating the stability in colloidal suspensions [Tohver et al. Proc. Natl. Acad. Sci. U.S.A. 2001, 98, 8951] and is discussed as well. Finally, the estimated boundary of the free of macroion region between surface monolayers correlates well with the location of the boundary of the so-called "vacuum" phase that has been observed experimentally in an aqueous suspension of charged polystyrene spheres bounded by electrostatically repulsive glass walls [Pieransky et al. Phys. Rev. Lett. 1983, 50, 900].     

Combining molecular dynamics with Lattice Boltzmann: a hybrid method for the simulation of (charged) colloidal systems

We present a hybrid method for the simulation of colloidal systems that combines molecular dynamics (MD) with the Lattice Boltzmann (LB) scheme. The LB method is used as a model for the solvent in order to take into account the hydrodynamic mass and momentum transport through the solvent. The colloidal particles are propagated via MD and they are coupled to the LB fluid by viscous forces. With respect to the LB fluid, the colloids are represented by uniformly distributed points on a sphere. Each such point [with a velocity V(r) at any off-lattice position r] is interacting with the neighboring eight LB nodes by a frictional force F = xi0(V(r)-u(r)), with xi0 being a friction coefficient and u(r) being the velocity of the fluid at the position r. Thermal fluctuations are introduced in the framework of fluctuating hydrodynamics. This coupling scheme has been proposed recently for polymer systems by Ahlrichs and Dunweg [J. Chem. Phys. 111, 8225 (1999)]. We investigate several properties of a single colloidal particle in a LB fluid, namely, the effective Stokes friction and long-time tails in the autocorrelation functions for the translational and rotational velocity. Moreover, a charged colloidal system is considered consisting of a macroion, counterions, and coions that are coupled to a LB fluid. We study the behavior of the ions in a constant electric field. In particular, an estimate of the effective charge of the macroion is yielded from the number of counterions that move with the macroion in the direction of the electric field.     

Dielectric saturation of the ion hydration shell and interaction between two double helices of DNA in mono- and multivalent electrolyte solutions: foundations of the epsilon-modified Poisson-Boltzmann theory

Potentials of mean force between single Na+, Ca2+, and Mg2+ cations and a highly charged spherical macroion in SPC/E water have been determined using molecular dynamics simulations. Results are compared to the electrostatic energy calculations for the primitive polarization model (PPM) of hydrated cations describing the ion hydration shell as a dielectric sphere of low permittivity (Gavryushov, S.; Linse, P. J. Phys. Chem. B 2003, 107, 7135). Parameters of the ion dielectric sphere and radius of the macroion/water dielectric boundary were extracted by means of this comparison to approximate the short-range repulsion of ions near the interface. To explore the counterion distributions around a simplified model of DNA, the obtained PPM parameters for Na+ and Ca2+ have been substituted into the modified Poisson-Boltzmann (MPB) equations derived for the PPM and named the epsilon-MPB (epsilon-MPB) theory. epsilon-MPB results for DNA suggest that such polarization effects are important in the case of 2:1 electrolyte and highly charged macromolecules. The three-dimensional implementation of the epsilon-MPB theory was also applied to calculation of the energies of interaction between two parallel macromolecules of DNA in solutions of NaCl and CaCl2. Being compared to results of MPB calculations without the ion polarization effects, it suggests that the ion hydration shell polarization and inhomogeneous solvent permittivity might be essential factors in the experimentally known hydration forces acting between charged macromolecules and bilayers at separations of less than 20 A between their surfaces.     

Macroion solutions in the cell model studied by field theory and Monte Carlo simulations

Aqueous solutions of charged spherical macroions with variable dielectric permittivity and their associated counterions are examined within the cell model using a field theory and Monte Carlo simulations. The field theory is based on separation of fields into short- and long-wavelength terms, which are subjected to different statistical-mechanical treatments. The simulations were performed by using a new, accurate, and fast algorithm for numerical evaluation of the electrostatic polarization interaction. The field theory provides counterion distributions outside a macroion in good agreement with the simulation results over the full range from weak to strong electrostatic coupling. A low-dielectric macroion leads to a displacement of the counterions away from the macroion.     

Radiative and nonradiative decay rates of a molecule close to a metal particle of complex shape

We present a model to evaluate the radiative and nonradiative lifetimes of electronic excited states of a molecule close to a metal particle of complex shape and, possibly, in the presence of a solvent. The molecule is treated quantum mechanically at Hartree-Fock (HF) or density-functional theory (DFT) level. The metal/solvent is considered as a continuous body, characterized by its frequency dependent local dielectric constant. For simple metal shapes (planar infinite surface and spherical particle) a version of the polarizable continuum model based on the integral equation formalism has been used, while an alternative methodology has been implemented to treat metal particles of arbitrary shape. In both cases, equations have been numerically solved using a boundary element method. Excitation energies and nonradiative decay rates due to the energy transfer from the molecule to the metal are evaluated exploiting the linear response theory (TDHF or TDDFT where TD--time dependent). The radiative decay rate of the whole system (molecule + metal/solvent) is calculated, still using a continuum model, in terms of the response of the surrounding to the molecular transition. The model presented has been applied to the study of the radiative and nonradiative lifetimes of a lissamine molecule in solution (water) and close to gold spherical nanoparticles of different radius. In addition, the influence of the metal shape has been analyzed by performing calculations on a system composed by a coumarin-type molecule close to silver aggregates of complex shape.     

Ion size effects on the electrokinetics of salt-free concentrated suspensions in ac fields

We analyze the influence of finite ion size effects in the response of a salt-free concentrated suspension of spherical particles to an oscillating electric field. Salt-free suspensions are just composed of charged colloidal particles and the added counterions released by the particles to the solution that counterbalance their surface charge. In the frequency domain, we study the dynamic electrophoretic mobility of the particles and the dielectric response of the suspension. We find that the Maxwell-Wagner-O’Konski process associated with the counterions condensation layer is enhanced for moderate to high particle charges, yielding an increment of the mobility for such frequencies. We also find that the increment of the mobility grows with ion size and particle charge. All these facts show the importance of including ion size effects in any extension attempting to improve standard electrokinetic models.     

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.     

Molecular Dynamics Simulation of the Electrical Double Layer of Spherical Macroion

The structure of the electrical double layer (EDL) of a spherical macroion with a total charge of 60 elementary charges is studied by molecular dynamics methods. In calculations we used two models: continuous and discrete. In the continuous model, the total charge was concentrated in the center of the macroion; in the discrete model, elementary charges were randomly distributed over the surface of the macroion. The radial profiles of local densities and electric potential in EDL, as well as the degree of counterion binding by the macroion, are calculated with allowance for the Lennard-Jones and electrostatic interactions. It is established that the character of charge distribution significantly affects the EDL structure near the macroion, whereas its effect is much weaker at larger distances. The results obtained are compared with the experimental data on the surface potential and the diffuse part of EDL of sodium dodecyl sulfate micelles in aqueous solution, as well as on the micelle-bound charge. It is shown that even weak specific interaction between counterions and a macroion can substantially influence the structure of its EDL.     

Motion of a colloidal sphere with interfacial self-electrochemical reactions induced by a magnetic field

The motion of a spherical colloidal particle with spontaneous electrochemical reactions occurring on its surface in an ionic solution subjected to an applied magnetic field is analyzed for an arbitrary zeta potential distribution. The thickness of the electric double layer adjacent to the particle surface is assumed to be much less than the particle radius. The solutions of the Laplace equations governing the magnetic scalar potential and electric potential, respectively, lead to the magnetic flux and electric current density distributions in the particle and fluid phases of arbitrary magnetic permeabilities and electric conductivities. The Stokes equations modified with the Lorentz force contribution for the fluid motion are dealt by using a generalized reciprocal theorem, and closed-form formulas for the translational and angular velocities of the colloidal sphere induced by the magnetohydrodynamic effect are obtained. The dipole and quadrupole moments of the zeta potential distribution over the particle surface cause the particle translation and rotation, respectively. The induced velocities of the particle are unexpectedly significant, and their dependence on the characteristics of the particle-fluid system is physically different from that for electromagnetophoretic particles or phoretic swimmers.     

Nonlocal continuum electrostatic theory predicts surprisingly small energetic penalties for charge burial in proteins

We study the energetics of burying charges, ion pairs, and ionizable groups in a simple protein model using nonlocal continuum electrostatics. Our primary finding is that the nonlocal response leads to markedly reduced solvent screening, comparable to the use of application-specific protein dielectric constants. Employing the same parameters as used in other nonlocal studies, we find that for a sphere of radius 13.4 ? containing a single +1e charge, the nonlocal solvation free energy varies less than 18 kcal/mol as the charge moves from the surface to the center, whereas the difference in the local Poisson model is ～35 kcal/mol. Because an ion pair (salt bridge) generates a comparatively more rapidly varying Coulomb potential, energetics for salt bridges are even more significantly reduced in the nonlocal model. By varying the central parameter in nonlocal theory, which is an effective length scale associated with correlations between solvent molecules, nonlocal-model energetics can be varied from the standard local results to essentially zero; however, the existence of the reduction in charge-burial penalties is quite robust to variations in the protein dielectric constant and the correlation length. Finally, as a simple exploratory test of the implications of nonlocal response, we calculate glutamate pK(a) shifts and find that using standard protein parameters (ε(protein) = 2-4), nonlocal results match local-model predictions with much higher dielectric constants. Nonlocality may, therefore, be one factor in resolving discrepancies between measured protein dielectric constants and the model parameters often used to match titration experiments. Nonlocal models may hold significant promise to deepen our understanding of macromolecular electrostatics without substantially increasing computational complexity.     

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.     

Dielectric properties of sols of silver nanoparticles capped by alkyl carboxylate ligands

Sols of silver nanoparticles in toluene were studied by broadband dielectric spectroscopy (10⁻³–10⁵ Hz). The frequency dependences of the specific alternating current (ac) conductivity and the complex electric modulus were used to estimate the temperature/frequency intervals of long- and short-range charge transfer occurs, respectively. A considerable increase (by more than 30 °C) in the Vogel temperature T ₀ and the glass transition temperature T _g in sols compared with the pure solvent was found. It can be hypothesized that these cooperative effects reflect the initial stage of the superlattice formation. Although the dielectric characteristics of sols are generally controlled by the conductivity relaxation, the dielectric response was observed in the high-frequency range (1–10³ Hz) at low temperatures (from −50 to +10 °C). This response results from the presence of nanoparticles in solution. It is supposed that the relaxation is caused by the motion of ion impurities on the Ag nanoparticle surface within the carboxylate ligands shell. The dielectric properties of films strongly depend on both the characteristics of nanoparticles and the conditions of the film preparation. Like in sols, the direct current (dc) conductivity and the dielectric response of Ag nanoparticles in films are due to ion impurities.     

Transient phoretic migration of a permselective colloidal particle

We quantify the phoretic migration of a spherical cation-permselective colloidal particle immersed in a binary electrolyte under a time-dependent electric field. We invoke the thin-Debye-layer approximation, where the size of ionic Debye layer enveloping the particle is much smaller than the particle radius. The imposed electric field generates ion concentration gradients, or concentration polarization, in the bulk (electroneutral) electrolyte outside the Debye layer. The bulk ion concentration polarization--and consequently the particle's phoretic velocity--evolves on the time scale for ion diffusion around the particle, which can be on the order of milliseconds for typical colloidal dimensions. Notably, concentration polarization arises here solely due to the permselectivity of the particle; it does not require non-uniform ionic transport in the Debye layer (i.e., surface conduction). Thus, the phoretic transport of a permselective particle is significantly different to that of a inert, dielectric particle, since surface conduction is necessary to achieve bulk concentration polarization in the (more commonly studied) latter case. Calculations are presented for a permselective particle under oscillatory (ac) and suddenly applied electric fields. In the former case, the particle velocity possesses frequency-dependent components in phase and out of phase with the driving field; in the latter case, the particle approaches its terminal velocity with a long-time (algebraic) tail.     

Towards a modellisation of the solvation energy in multi-component solvents. The interesting case of a charged solute imbedded in a polymer-containing electrolyte solution

In this paper we propose a mean-field theory to calculate the solvation free energy of a charged solute imbedded in a complex multi-component solvent. We considered a solvent made up of a mixture of small (electrolyte solution) and large (polymer) components. The presence of macromolecules ensures reduced mixing entropy among the different solvent components, an effect due to polymer connectivity. The reduced entropy favours strong preferential distribution of a particular solvent even in the presence of weak preferential solute–solvent interactions. In addition, two energy terms must be considered: (a) the interaction between the solute electrostatic potential and the electrolyte solution and (b) the formation of a polymer–solute interface. Because of the different dielectric permittivity of the solvent components, the electrolyte and polymer distribution functions are strongly coupled: ions, indeed, are more solvated in regions of higher local dielectric permittivity arising from the inhomogeneous mixing of solvent and polymer. We combined together the different energy terms in the framework of the de Gennes free energy functional for polymer solutions along with a generalised Poisson–Boltzmann equation developed for inhomogeneous dielectric media. Moreover, the preferential electrolyte solvation in regions of greater polarity was considered by an extension of the Born equation. Setting the polymer dielectric permittivity smaller than the solvent one and making null the specific polymer–solute interactions, we calculated enhanced electrolyte concentration and reduced polymer concentration near the solute surface on raising the solute surface charge density. The theory shows also the breakdown of the widely used separation between electrostatic and surface tension-dependent contributions to solvation energy when non-ideal mixed solvents are considered. In fact, according to the model, the surface tension of such mixed solvents strongly depends on the solute surface charge density: at high potentials the interfacial tension may increase rather than decrease on raising the polymer volume fraction. The theoretical results have been compared with experimental data on polymer+electrolyte solution surface tension and with solubility data of colloidal particles. The comparison evidences the complex behaviour of multi-component solvents going well beyond the trivial weighted average of the dielectric permittivity and surface tension of the isolated chemical components. Deviations from the simple behaviour predicted by an average picture of multi-component solvents could be understood by developing more sophisticated, but still simple, approaches like that proposed in this paper.Contribution to the Jacopo Tomasi Honorary Issue. This paper is dedicated to Jacopo Tomasi. I learned much of the difficult art of transforming complex problems into simple models after reading his early works on solvation energy.