Strong and weak adsorptions of polyelectrolyte chains onto oppositely charged spheres

We investigate the complexation of long thin polyelectrolyte (PE) chains with oppositely charged spheres. In the limit of strong adsorption, when strongly charged PE chains adapt a definite wrapped conformation on the sphere surface, we analytically solve the linear Poisson-Boltzmann equation and calculate the electrostatic potential and the energy of the complex. We discuss some biological applications of the obtained results. For weak adsorption, when a flexible weakly charged PE chain is localized next to the sphere in solution, we solve the Edwards equation for PE conformations in the Hulthen potential, which is used as an approximation for the screened Debye-Huckel potential of the sphere. We predict the critical conditions for PE adsorption. We find that the critical sphere charge density exhibits a distinctively different dependence on the Debye screening length than for PE adsorption onto a flat surface. We compare our findings with experimental measurements on complexation of various PEs with oppositely charged colloidal particles. We also present some numerical results of the coupled Poisson-Boltzmann and self-consistent field equation for PE adsorption in an assembly of oppositely charged spheres.     

Electrostatic screening and energy barriers of ions in low-dielectric membranes

We present exact solutions of the linear Poisson-Boltzmann equation for several problems relevant for ion translocation across low-dielectric membranes. Our results are obtained for a finite Debye screening length, and they generalize the classical results for pure Coulombic electrostatics (Parsegian, A. Nature (London) 1969, 221, 844). We calculate the electrostatic self-energy of an ion in the middle of a low-dielectric slab, its energy inside a cylindrical high-dielectric pore, and its energy inside a high-dielectric spherical jacket. We consider also the influence of negative charges distributed on the walls of the cylindrical pore. We show that ion self-energy barriers are considerably reduced due to screening of electrolyte. We compare our results with some numerical results for screened electrostatics of ion channels and wide pores.     

Electrokinetic motion of a charged colloidal sphere in a spherical cavity with magnetic fields

The magnetohydrodynamic (MHD) effects on the translation and rotation of a charged colloidal sphere situated at the center of a spherical cavity filled with an arbitrary electrolyte solution when a constant magnetic field is imposed are analyzed at the quasisteady state. The electric double layers adjacent to the solid surfaces may have an arbitrary thickness relative to the particle and cavity radii. Through the use of a perturbation method to the leading order, the Stokes equations modified with the electric∕Lorentz force term are dealt by using a generalized reciprocal theorem. Using the equilibrium double-layer potential distribution in the fluid phase from solving the linearized Poisson-Boltzmann equation, we obtain explicit formulas for the translational and angular velocities of the colloidal sphere produced by the MHD effects valid for all values of the particle-to-cavity size ratio. For the limiting case of an infinitely large cavity with an uncharged wall, our result reduces to the relevant solution for an unbounded spherical particle available in the literature. The boundary effect on the MHD motion of the spherical particle is a qualitatively and quantitatively sensible function of the parameters a∕b and κa, where a and b are the radii of the particle and cavity, respectively, and κ is the reciprocal of the Debye screening length. In general, the proximity of the cavity wall reduces the MHD migration but intensifies the MHD rotation of the particle.     

A new and simpler approach for the solution of the electrostatic potential differential equation. Enhanced solution for planar, cylindrical and annular geometries

The motivation of the study performed in this project is focused on deriving a more effective, accurate, and mathematically friendly solution for the prediction of the electrostatic potential, commonly used on electrokinetic research and its related applications. In this contribution, based on the Debye-Hückel approximation, a new solution strategy for the differential equations of the electrostatic potential is proposed. In fact, a simple predictor-corrector calculation is developed to achieve more accurate predictions of electrostatic potential profiles. Furthermore, in this study the authors introduce the correction function f(AO) to the inverse Debye length, lambda. The f(AO) function improves the Debye-Hückel approximation and it is a recursive function of the electrical potential. Once the inverse Debye length, lambda, has been corrected by the f(AO) function and introduced in the simplified solution of the Poisson-Boltzmann equation (i.e., the linear approximation, due to Debye and Hückel), the electrostatic potential outcome little differs from the numerical solution of the complete (nonlinear) differential equation. This new approach embraces different geometries of interest, such as planar, cylindrical, and annular, with excellent results in all the cases and for a wide range of electrostatic potential values. This new predicting semi-analytical technique can be a useful tool on electrical field applications such as the separation of a mixture of macromolecules and the removal of contaminants in soil cleaning processes. Illustrative results are presented for the geometries identified above.     

Electrophoresis and electric conduction in a salt-free suspension of charged particles

The electrophoresis and electric conduction of a suspension of charged spherical particles in a salt-free solution are analyzed by using a unit cell model. The linearized Poisson-Boltzmann equation (valid for the cases of relatively low surface charge density or high volume fraction of the particles) and Laplace equation are solved for the equilibrium electric potential profile and its perturbation caused by the imposed electric field, respectively, in the fluid containing the counterions only around the particle, and the ionic continuity equation and modified Stokes equations are solved for the electrochemical potential energy and fluid flow fields, respectively. Explicit analytical formulas for the electrophoretic mobility of the particles and effective electric conductivity of the suspension are obtained, and the particle interaction effects on these transport properties are significant and interesting. The scaled zeta potential, electrophoretic mobility, and effective electric conductivity increase monotonically with an increase in the scaled surface charge density of the particles and in general decrease with an increase in the particle volume fraction, keeping each other parameter unchanged. Under the Debye-Hückel approximation, the dependence of the electrophoretic mobility normalized with the surface charge density on the ratio of the particle radius to the Debye screening length and particle volume fraction in a salt-free suspension is same as that in a salt-containing suspension, but the variation of the effective electric conductivity with the particle volume fraction in a salt-free suspension is found to be quite different from that in a suspension containing added electrolyte.     

Electric double layer effect on observable characteristics of the tunnel current through a bridged electrochemical contact

Scanning tunneling microscopy and electrical conductivity of redox molecules in conducting media (aqueous or other media) acquire increasing importance both as novel single-molecule science and with a view on molecular scale functional elements. Such configurations require full and independent electrochemical potential control of both electrodes involved. We provide here a general formalism for the electric current through a redox group in an electrochemical tunnel contact. The formalism applies broadly in the limits of both weak and strong coupling of the redox group with the enclosing metal electrodes. Simple approximate expressions better suited for experimental data analysis are also derived. Particular attention is given to the effects of the Debye screening of the electric potential in the narrow tunneling gap based on the limit of the linearized Poisson-Boltzmann equation. The current/overpotential relation shows a maximum at a position which depends on the ionic strength. It is shown, in particular, that the dependence of the maximum position on the bias voltage may be nonmonotonous. Approximate expressions for the limiting value of the slope of the current/overpotential dependence and the width of the maximum on the bias voltage are also given and found to depend strongly on both the Debye screening and the position of the redox group in the tunnel gap, with diagnostic value in experimental data analysis.     

Effective Debye length in closed nanoscopic systems: a competition between two length scales

The Poisson-Boltzmann equation (PBE) is widely employed in fields where the thermal motion of free ions is relevant, in particular in situations involving electrolytes in the vicinity of charged surfaces. The applications of this non-linear differential equation usually concern open systems (in osmotic equilibrium with an electrolyte reservoir, a semi-grand canonical ensemble), while solutions for closed systems (where the number of ions is fixed, a canonical ensemble) are either not appropriately distinguished from the former or are dismissed as a numerical calculation exercise. We consider herein the PBE for a confined, symmetric, univalent electrolyte and quantify how, in addition to the Debye length, its solution also depends on a second length scale, which embodies the contribution of ions by the surface (which may be significant in high surface-to-volume ratio micro- or nanofluidic capillaries). We thus establish that there are four distinct regimes for such systems, corresponding to the limits of the two parameters. We also show how the PBE in this case can be formulated in a familiar way by simply replacing the traditional Debye length by an effective Debye length, the value of which is obtained numerically from conservation conditions. But we also show that a simple expression for the value of the effective Debye length, obtained within a crude approximation, remains accurate even as the system size is reduced to nanoscopic dimensions, and well beyond the validity range typically associated with the solution of the PBE.     

Electrostatics of DNA complexes with cationic lipid membranes

We present the exact solutions of the linear Poisson-Boltzmann equation for several problems relevant to electrostatics of DNA complexes with cationic lipids. We calculate the electrostatic potential and electrostatic energy for lamellar and inverted hexagonal phases, concentrating on the effects of dielectric boundaries. We compare our results for the complex energy with the known results of numerical solution of the nonlinear Poisson-Boltzmann equation. Using the solution for the lamellar phase, we calculate the compressibility modulus and compare our findings with the experimental data available. Also, we treat charge-charge interactions across, along, and between two low-dielectric membranes. We obtain an estimate for the strength of electrostatic interactions of one-dimensional DNA smectic layers across the lipid membrane. We discuss in the end some aspects of two-dimensional DNA condensation and DNA-DNA attraction in the DNA-lipid lamellar phase in the presence of di- and trivalent cations. We analyze the equilibrium DNA-DNA separations in lamellar complexes using the recently developed theory of electrostatic interactions of DNA helical charge motifs.     

Transport properties of long straight nano-channels in electrolyte solutions: a systematic approach

The principle of local thermodynamic equilibrium is systematically employed for obtaining various transport properties of long straight nano-channels. The concept of virtual solution is used to describe situations of non-negligible overlap of diffuse parts of electric double layers (EDLs) in nano-channels. Generic expressions for a variety of transport properties of long straight nano-channels are obtained in terms of quasi-equilibrium distribution coefficients of ions and functionals of quasi-equilibrium distribution of electrostatic potential. Further, the Poisson-Boltzmann approach is used to specify these expressions for long straight slit-like nano-channels. In the approximation of non-overlapped diffuse parts of double electric layers in nano-channels, simple analytical expressions are obtained for the apparent electrophoretic mobilities of (trace) analytes of arbitrary charge as well as for the salt reflection coefficient (osmotic pressure), salt diffusion permeability and electro-viscosity (electrokinetic energy conversion). The approximate solutions are compared with the results of rigorous solution of non-linearized Poisson-Boltzmann equation, and the accuracy of approximation is shown to be typically excellent when the nano-channel half-height exceeds ca.3 Debye screening lengths. Due to non-negligible electrostatic adsorption of ions by nano-channels, the apparent electrophoretic mobilities of counter-ionic analytes in nano-channels are smaller than in micro-channels whereas those of co-ionic analytes are larger. This dependence on the charge is useful for the separation of analytes of close electrophoretic mobilities. The osmotic pressure is shown to be positive, negative or pass through maxima as a function of applied salt-concentration difference within a fairly narrow range of ratios of nano-channel height to the Debye screening length. The diffusion permeability of charged nano-channels to single salts is demonstrated (for the first time) to be typically larger than that of neutral nano-channels of the same dimensions due to electrical facilitation of salt diffusion.     

Three-dimensional discrete charge effects on the electrostatic interaction between two soft particles

An expression for the electrostatic interaction energy between two parallel plate-like soft particles (i.e., hard particles covered with an ion-penetrable surface layer of polyelectrolytes) in an electrolyte solution is derived by using the linearized Poisson-Boltzmann equation. This expression is based on a discrete charge model in which the surface layer consists of a cubic lattice of fixed point charges. We show that the deviation of the results of the discrete charge model from those of the conventional smeared charge model becomes significant as the ratio of the lattice spacing to the Debye length becomes large. As this ratio decreases, on the other hand, the discrete charge model approaches a smeared charge model, leading to the Donnan-potential regulated interaction model.     

Comparison of the Nernst-Planck model and the Poisson-Boltzmann model for electroosmotic flows in microchannels

In the analysis of electroosmotic flows, the internal electric potential is usually modeled by the Poisson-Boltzmann equation. The Poisson-Boltzmann equation is derived from the assumption of thermodynamic equilibrium where the ionic distributions are not affected by fluid flows. Although this is a reasonable assumption for steady electroosmotic flows through straight microchannels, there are some important cases where convective transport of ions has nontrivial effects. In these cases, it is necessary to adopt the Nernst-Planck equation instead of the Poisson-Boltzmann equation to model the internal electric field. In the present work, the predictions of the Nernst-Planck equation are compared with those of the Poisson-Boltzmann equation for electroosmotic flows in various microchannels where the convective transport of ions is not negligible.     

Magnetohydrodynamic effects on a charged colloidal sphere with arbitrary double-layer thickness

An analytical study is presented for the magnetohydrodynamic (MHD) effects on a translating and rotating colloidal sphere in an arbitrary electrolyte solution prescribed with a general flow field and a uniform magnetic field at a steady state. The electric double layer surrounding the charged particle may have an arbitrary thickness relative to the particle radius. Through the use of a simple perturbation method, the Stokes equations modified with an electric force term, including the Lorentz force contribution, are dealt by using a generalized reciprocal theorem. Using the equilibrium double-layer potential distribution from solving the linearized Poisson-Boltzmann equation, we obtain closed-form formulas for the translational and angular velocities of the spherical particle induced by the MHD effects to the leading order. It is found that the MHD effects on the particle movement associated with the translation and rotation of the particle and the ambient fluid are monotonically increasing functions of κa, where κ is the Debye screening parameter and a is the particle radius. Any pure rotational Stokes flow of the electrolyte solution in the presence of the magnetic field exerts no MHD effect on the particle directly in the case of a very thick double layer (κa→0). The MHD effect caused by the pure straining flow of the electrolyte solution can drive the particle to rotate, but it makes no contribution to the translation of the particle.     

Electrostatic stretching of a charged vesicle

We present a closed-form solution of electrostatic potential self-induced by a uniformly charged micro/nanovesicle and the corresponding elastic deformation of the vesicle membrane due to Maxwell stress. At equilibrium, the electrostatic force induced on both sides of the membrane is balanced by the elastic force of the stretched membrane. We develop differential and integral solutions of the coupled Poisson-Boltzmann system for a spherical vesicle and demonstrate that the integral solution is relatively flexible in formulating asymmetric configurations. Analytical results are formulated in terms of vesicle size, Debye length, and the surface charge density. The membrane stretching is characterized by the dimensionless group that defines the relative strength of the net electric force with respect to the membrane stiffness. We found that the self-induced electrostatic interaction will lead to a pre-stressed membrane although the small displacement is often negligible compared with the vesicle size. Quantitative analysis also reveals that the electric force can assist the vesicle in recovering its opening pore.     

Screening of charged spheroidal colloidal particles

We study the effective screened electrostatic potential created by a spheroidal colloidal particle immersed in an electrolyte, within the mean field approximation, using Poisson-Boltzmann equation in its linear and nonlinear forms, and also beyond the mean field by means of Monte Carlo computer simulation. The anisotropic shape of the particle has a strong effect on the screened potential, even at large distances (compared to the Debye length) from it. To quantify this anisotropy effect, we focus our study on the dependence of the potential on the position of the observation point with respect with the orientation of the spheroidal particle. For several different boundary conditions (constant potential, or constant surface charge) we find that, at large distance, the potential is higher in the direction of the large axis of the spheroidal particle.     

Line of charges in electrolyte solution near a half-space II. Electric field of a single charge

The electric potential of a single charge in electrolyte solution near a dielectric or a semiconducting half-space is determined in closed form when the electrostatics is described by the linear Debye-Hückel (D-H) equation. The electric potential strongly depends on the Debye length of the solution, the substrate-to-solution dielectric constant ratio, and the Debye length of the semiconductor. The technique of Hankel transforms is shown to be a useful tool in solving such axially symmetric boundary value problems for the D-H equation.     

A simple polarizable continuum solvation model for electrolyte solutions

We propose a Debye-Hu?ckel-like screening model (DESMO) that generalizes the familiar conductor-like screening model (COSMO) to solvents with non-zero ionic strength and furthermore provides a numerical generalization of the Debye-Hu?ckel model that is applicable to non-spherical solute cavities. The numerical implementation of DESMO is based upon the switching/Gaussian (SWIG) method for smooth cavity discretization, which we have recently introduced in the context of polarizable continuum models (PCMs). This approach guarantees that the potential energy is a smooth function of the solute geometry and analytic gradients for DESMO are reported here. The SWIG formalism also facilitates analytic implementation of two other PCMs that are based on a screened Coulomb potential: the "integral equation formalism" (IEF-PCM) and the "surface and simulation of volume polarization for electrostatics" [SS(V)PE] method. Fully analytic implementations of these screened PCMs are reported here for the first time. Numerical results, for model systems where an exact solution of the linearized Poisson-Boltzmann equation is available, demonstrate that these screened PCMs are highly accurate. In realistic test cases, they are as accurate as the best available three-dimensional finite-difference methods. In polar solvents, DESMO is nearly as accurate as more sophisticated screened PCMs, but is significantly simpler and more efficient.     

Electrical double layer interactions for spherical charge regulating colloidal particles

A new way of linearising the Poisson-Boltzmann equation is presented which is capable of producing accurate answers for moderately high surface charge and surface potentials. Unlike the Debye-Hückel approximation, the linearisation is based on the fact that electric potential variation, rather than the magnitude of the potential itself, remains small between two charged plates at close separation. The method can be applied quite generally to any type of surface. The cases of constant surface charge, the constant surfaceppotential and charge regulating systems are considered as examples. The results for constant charge and potentials as high as 5 μC/cm² and 200 mV are within 7% of the exact answers for plate separations of one Debye length or less. In contrast, the Debye-Hückel approximation produces answers that are at least an order of magnitude too large, for the same problems. The application of the method to the problem of calculating the electrostatic forces between charge regulated spherical particles is also addressed in some detail.     

Double layer interaction between two plates with polyelectrolyte brushes

When polyelectrolyte chains are grafted to colloidal particles, the electric field between particles is affected by the charges of the chains. In some previous theoretical attempts, the charge density of the polyelectrolyte chains per unit length was considered constant, and its effect was accounted for by introducing an additional constant charge density into the unidimensional Poisson-Boltzmann equation, which was evaluated assuming that it is uniformly distributed in the polyelectrolyte volume of the brush. In this paper, a more detailed model is employed for the calculation of the electrical potential between two plates on which polyelectrolyte brushes are present. In this model, the polyelectrolyte chain is viewed as a rigid cylinder, on the surface of which charges are generated through the dissociation of ionizable sites and adsorption of the cations of the electrolyte. To each of the chains an atmosphere is attached which for simplicity is assumed cylindrical. In the brush region, the electrical potential is described by a two-dimensional Poisson-Boltzmann equation, while in the region free of polyelectrolyte chains by a unidimensional Poisson-Boltzmann equation. Such a model is physically suitable when the charges of the chains are sufficiently large for the repulsion they generate to ensure that the chains are fully extended. Such cases are quite frequent, because relatively low charges lead to an almost complete extension of the chains. In this paper, both the plate surface and the surface of the cylinders are considered charged. The effects of electrolyte concentration, pH, brush thickness and chain coverage density on the repulsion between plates are examined.     

A new approach for efficient simulation of Coulomb interactions in ionic fluids

We propose a simplified version of local molecular field (LMF) theory to treat Coulomb interactions in simulations of ionic fluids. LMF theory relies on splitting the Coulomb potential into a short-ranged part that combines with other short-ranged core interactions and is simulated explicitly. The averaged effects of the remaining long-ranged part are taken into account through a self-consistently determined effective external field. The theory contains an adjustable length parameter sigma that specifies the cutoff distance for the short-ranged interaction. This can be chosen to minimize the errors resulting from the mean-field treatment of the complementary long-ranged part. Here we suggest that in many cases an accurate approximation to the effective field can be obtained directly from the equilibrium charge density given by the Debye theory of screening, thus eliminating the need for a self-consistent treatment. In the limit sigma-->0, this assumption reduces to the classical Debye approximation. We examine the numerical performance of this approximation for a simple model of a symmetric ionic mixture. Our results for thermodynamic and structural properties of uniform ionic mixtures agree well with similar results of Ewald simulations of the full ionic system. In addition, we have used the simplified theory in a grand-canonical simulation of a nonuniform ionic mixture where an ion has been fixed at the origin. Simulations using short-ranged truncations of the Coulomb interactions alone do not satisfy the exact condition of complete screening of the fixed ion, but this condition is recovered when the effective field is taken into account. We argue that this simplified approach can also be used in the simulations of more complex nonuniform systems.