New results from the contact theorem for the charge profile for symmetric electrolytes

In this paper the contact value of the charge profile at a charged interface is presented as the sum of the normal component of the Maxwell electrostatic tensor and a new electrostatic property defined by the integral from the product of the gradient of the electrical potential and the singlet distribution function of coions (ions with sign of the charge equal to that of the interface). On physical arguments, it is conjectured that this new property is a monotonic function of the electrical charge at the wall and is limited by the bulk electrolyte pressure for large electrical charges at the wall. Using the contact theorems for the density and the charge profiles, the exact expressions for the contact values of the profiles of coions and counterions are derived and some related general properties are discussed.     

Entropic effects in the electric double layer of model colloids with size-asymmetric monovalent ions

The structure of the electric double layer of charged nanoparticles and colloids in monovalent salts is crucial to determine their thermodynamics, solubility, and polyion adsorption. In this work, we explore the double layer structure and the possibility of charge reversal in relation to the size of both counterions and coions. We examine systems with various size-ratios between counterions and coions (ion size asymmetries) as well as different total ion volume fractions. Using Monte Carlo simulations and integral equations of a primitive-model electric double layer, we determine the highest charge neutralization and electrostatic screening near the electrified surface. Specifically, for two binary monovalent electrolytes with the same counterion properties but differing only in the coion's size surrounding a charged nanoparticle, the one with largest coion size is found to have the largest charge neutralization and screening. That is, in size-asymmetric double layers with a given counterion's size the excluded volume of the coions dictates the adsorption of the ionic charge close to the colloidal surface for monovalent salts. Furthermore, we demonstrate that charge reversal can occur at low surface charge densities, given a large enough total ion concentration, for systems of monovalent salts in a wide range of ion size asymmetries. In addition, we find a non-monotonic behavior for the corresponding maximum charge reversal, as a function of the colloidal bare charge. We also find that the reversal effect disappears for binary salts with large-size counterions and small-size coions at high surface charge densities. Lastly, we observe a good agreement between results from both Monte Carlo simulations and the integral equation theory across different colloidal charge densities and 1:1-electrolytes with different ion sizes.     

Dressed ion theory of size-asymmetric electrolytes: effective ionic charges and the decay length of screened Coulomb potential and pair correlations

The effects of ionic size asymmetry on long-range electrostatic interactions in electrolyte solutions are investigated within the primitive model. Using the formalism of dressed ion theory we analyze correlation functions from Monte Carlo simulations and the hypernetted chain approximation for size asymmetric 1:1 electrolytes. We obtain decay lengths of the screened Coulomb potential, effective charges of ions, and effective permittivity of the solution. It is found that the variation of these quantities with the degree of size asymmetry depends in a quite intricate manner on the interplay between the electrostatic coupling and excluded volume effects. In most cases the magnitude of the effective charge of the small ion species is larger than that of the large species; the difference increases with increasing size asymmetry. The effective charges of both species are larger (in absolute value) than the bare ionic charge, except for high asymmetry where the effective charge of the large ions can become smaller than the bare charge.     

Methodological problems in pressure profile calculations for lipid bilayers

From molecular dynamics simulations of a dipalmitoyl-phosphatidyl-choline (DPPC) lipid bilayer in the liquid crystalline phase, pressure profiles through the bilayer are calculated by different methods. These profiles allow us to address two central and unresolved problems in pressure profile calculations: The first problem is that the pressure profile is not uniquely defined since the expression for the local pressure involves an arbitrary choice of an integration contour. We have investigated two different choices leading to the Irving-Kirkwood (IK) and Harasima (H) expressions for the local pressure tensor. For these choices we find that the pressure profile is almost independent of the contour used, which indicates that the local pressure is well defined for a DPPC bilayer in the liquid crystalline phase. This may not be the case for other systems and we therefore suggest that both the IK and H profiles are calculated in order to test the uniqueness of the profile. The second problem is how to include electrostatic interactions in pressure profile calculations when the simulations are conducted without truncating the electrostatic potential, i.e., using the Ewald summation technique. Based on the H expression for the local pressure, we present a method for calculating the contribution to the lateral components of the local pressure tensor from electrostatic interactions evaluated by the Ewald summation technique. Pressure profiles calculated with an electrostatic potential truncation (cutoff) from simulations conducted with Ewald summation are shown to depend on the cutoff in a subtle manner which is attributed to the existence of long-ranged charge ordering in the system. However, the pressure profiles calculated with relatively long cutoffs are qualitatively similar to the Ewald profile for the DPPC bilayer studied here.     

Structure of adsorption layers of ionic surfactants at the solid/liquid interface

A large number of experimental results of different surfactant adsorption systems (mainly on the silicas) obtained from both equilibrium and kinetic studies under different conditions are interpreted by a model of small individual surface aggregates. The adsorption model is contrasted with the influences of various factors, including electrostatic interaction, hydrophobic interaction, concentrations, types of coions, types of counterions, surfactant structure, alkyl chain length, types of head groups, neutral electrolytes, pH, adsorbent structure, porosity, surface charge density, and surface polarity.Dedicated to Frau Professor Dr. Elsa Ullmann on the occasion of her 80th birthday     

Effective multipoles and Yukawa electrostatics in dressed molecule theory

In this paper we derive the multipolar expansion of the screened Coulomb potential in electrolyte solutions with molecular solvent. The solute and solvent molecules can have arbitrary sizes, shapes, and internal charge distributions. We use the exact statistical mechanical definition of renormalized charge distributions coming from "dressed molecule theory" to determine the effective multipoles of a molecule immersed in an electrolyte. The effects of many-body correlations are fully included in our formally exact theory. We restrict ourselves to sufficiently dilute solutions so the screened Coulomb potential decays for large distances like a Yukawa function, exp(-kappa r)/r, where r is the distance and 1/kappa is the decay length (it is normally different from the Debye length). The resulting "Yukawa electrostatics" differ in many respects from ordinary, unscreened electrostatics. The "Yukawa charge" of a molecule (the lowest order moment in the multipolar expansion) is in general not equal to its Coulombic charge and it is not the integral of the renormalized charge distribution of the molecule. Moreover, as shown in this paper, the multipolar expansion of the Yukawa potential does not correspond, contrary to the case of the Coulomb potential, to its asymptotic expansion for large r. As a consequence, the charge term in the multipolar expansion is not the leading term in the asymptotic expansion. Instead, for large r values, multipoles of all orders contribute to the leading asymptotic term. Thus, the electrostatic potential from, for example, an electroneutral solvent molecule in an electrolyte solution has generally the same range as that from an ion. The proper asymptotic expansion for electrostatic interactions in electrolytes is derived. It is briefly shown how the multipole expansion formalism can also be applied in the Poisson-Boltzmann approximation for primitive model electrolytes.     

Concentration polarization at Langmuir monolayer deposition: the role of indifferent electrolytes

Under dynamic conditions of the charged Langmuir monolayer deposition onto a substrate surface, ion concentration and electric potential profiles are induced in the subphase around the three-phase contact line. Such local changes in the subphase influence the deposition process, particularly the monolayer adhesion work and the maximum deposition rate. If indifferent electrolytes (not interacting chemically with interfacial groups) are present in the solutions, they can affect electric potential distributions and therefore the monolayer charge and the deposition process as a whole. With increasing deposition rate, the indifferent electrolyte counterions replace gradually the potential-determining counterions in a close vicinity to the contact line. This leads to increasing monolayer ionization and increasing electrostatic repulsion between the monolayer and substrate. When the deposition rate approaches the critical one, the charge of the monolayer increases dramatically and the stationary monolayer deposition becomes impossible. Such a significant increase of the monolayer charge is not observed in the absence of indifferent electrolytes.     

An exact formula for the contact value of the density profile of a system of charged hard spheres near a charged wall

An exact formula for the contact value of the density of a system of charged hard spheres near a charged hard wall is obtained by means of a general statistical mechanical argument. In addition, a formula for the contact value of the charge profile in the limit of large field is obtained. Comparison with the corresponding expressions in the Poisson-Boltzmann theory of Gouy and Chapman shows that these latter expressions become exact for large fields, independent of the density of the hard spheres.     

Contact conditions for the charge in the theory of the electrical double layer

In this paper, from the Born-Green-Yvon equations of the liquid-state theory, we derive a general expression for the charge-density contact value at charged interfaces. This relation is discussed, in particular, for symmetrical electrolytes. We emphasize an essential coupling between the electric properties and the density profile. Limiting behavior at small and large charges at the interface is discussed.     

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.     

Structural thermodynamics of lamellar cationic lipid-DNA complex: DNA compressibility modulus

We have studied theoretically the compressibility modulus B of DNA and complexation adsorption isotherms of DNA and lipids, as a function of DNA spacing d(DNA) and NaCl electrolyte concentration, respectively, in isoelectric states of lamellar DNA/cationic lipid (CL) self-assemblies. The electrostatic free energy derived from the Poisson-Boltzmann theory predicts partial agreement with measured B values for interhelical separations d(DNA)>33 A when use is made of a fit of hydration repulsion from bulk DNA hexagonal phases in solution. For lower interchain separations the prediction worsens due to the hydration interaction that overcomes the electrostatic contribution. An exact match of the system's counterion electrochemical potentials and the coions of salt in aqueous phase leads to the electrostatic part of the free energy that renders isotherms of d(DNA) versus ionic strength in qualitative consistency with general trends of available experimental data of CL-DNA complexes.     

Fast multipole method for three-dimensional systems with periodic boundary condition in two directions

We derived a new expression for the electrostatic interaction of three-dimensional charge-neutral systems with two-dimensional periodic boundary conditions (slab geometry) using a fast multipole method (FMM). Contributions from all the image cells are expressed as a sum of real and reciprocal space terms, and a self-interaction term. The reciprocal space contribution consists of two parts: zero and nonzero terms of the absolute value of the reciprocal lattice vector. To test the new expressions, electrostatic interactions were calculated for a randomly placed charge distribution in a cubic box and liquid water produced by molecular dynamics calculation. The accuracy could be controlled by the degree of expansion of the FMM. In the present expression, the computational complexity of the electrostatic interaction of N-particle systems is order N, which is superior to that of the conventional two-dimensional periodic Ewald method for a slab geometry and the particle mesh Ewald method with a large empty space at an interface of the unit cell. © 2020 Wiley Periodicals, Inc.     

Anisotropic solvent model of the lipid bilayer. 1. Parameterization of long-range electrostatics and first solvation shell effects

A new implicit solvation model was developed for calculating free energies of transfer of molecules from water to any solvent with defined bulk properties. The transfer energy was calculated as a sum of the first solvation shell energy and the long-range electrostatic contribution. The first term was proportional to solvent accessible surface area and solvation parameters (σ(i)) for different atom types. The electrostatic term was computed as a product of group dipole moments and dipolar solvation parameter (η) for neutral molecules or using a modified Born equation for ions. The regression coefficients in linear dependencies of solvation parameters σ(i) and η on dielectric constant, solvatochromic polarizability parameter π*, and hydrogen-bonding donor and acceptor capacities of solvents were optimized using 1269 experimental transfer energies from 19 organic solvents to water. The root-mean-square errors for neutral compounds and ions were 0.82 and 1.61 kcal/mol, respectively. Quantification of energy components demonstrates the dominant roles of hydrophobic effect for nonpolar atoms and of hydrogen-bonding for polar atoms. The estimated first solvation shell energy outweighs the long-range electrostatics for most compounds including ions. The simplicity and computational efficiency of the model allows its application for modeling of macromolecules in anisotropic environments, such as biological membranes.     

Electrolyte distribution around two like-charged rods: their effective attractive interaction and angular dependent charge reversal

A simple model for two like-charged parallel rods immersed in an electrolyte solution is considered. We derived the three point extension (TPE) of the hypernetted chain/mean spherical approximation (TPE-HNC/MSA) and Poisson-Boltzmann (TPE-PB) integral equations. We numerically solve these equations and compare them to our results of Monte Carlo (MC) simulations. The effective interaction force, F(T), the charge distribution profiles, rho(el)(x,y), and the angular dependent integrated charge function, P(theta), are calculated for this system. The analysis of F(T) is carried out in terms of the electrostatic and entropic (depletion) contributions, F(E) and F(C). We studied several cases of monovalent and divalent electrolytes, for which the ionic size and concentration are varied. We find good qualitative agreement between TPE-HNC/MSA and MC in all the cases studied. The rod-rod force is found to be attractive when immersed in large size, monovalent or divalent electrolytes. In general, the TPE-PB has poor agreement with the MC. For large monovalent and divalent electrolytes, we find angular dependent charge reversal charge inversion and polarizability. We discuss the intimate relationship between this angular dependent charge reversal and rod-rod attraction.     

Computation of the Madelung constant for hypercubic crystal structures in any dimension

A new method of computing the Madelung constants for hypercubic crystal structures in any dimension \(n\ge 2\) is given. It is shown for \(n\ge 3\) that the Madelung constant may be obtained in a simple, efficient and unambiguous way as the Hadamard finite part of the integral representation of the potential within the crystal which is divergent at any point charge location. Such a regularization method fails in the bidimensional case due to the logarithmic nature of singularities for the potential. In that case, a specific approach is proposed taking in account the scale invariance of the Poisson equation and the existence of a finite horizon for each point charge in the plane. Since a closed-form exact solution for the 2D electrostatic potential may be derived, one shows that the Madelung constant may be defined via an appropriate limit calculation as the mean value of potential energies of charges composing the unit cell.     

Equation of state of charged colloidal suspensions and its dependence on the thermodynamic route

The thermodynamic properties of highly charged colloidal suspensions in contact with a salt reservoir are investigated in the framework of the renormalized Jellium model (RJM). It is found that the equation of state is very sensitive to the particular thermodynamic route used to obtain it. Specifically, the osmotic pressure calculated within the RJM using the contact value theorem can be very different from the pressure calculated using the Kirkwood-Buff fluctuation relations. On the other hand, Monte Carlo simulations show that both the effective pair potentials and the correlation functions are accurately predicted by the RJM. It is suggested that the lack of self-consistency in the thermodynamics of the RJM is a result of neglected electrostatic correlations between the counterions and coions.     

Analytical solution of Poisson–Boltzmann equation for interacting plates of arbitrary potentials and same sign

Efficient calculation of electrostatic interactions in colloidal systems is becoming more important with the advent of such probing techniques as atomic force microscopy. Such practice requires solving the nonlinear Poisson–Boltzmann equation (PBE). Unfortunately, explicit analytical solutions are available only for the weakly charged surfaces. Analysis of arbitrarily charged surfaces is possible only through cumbersome numerical computations. A compact analytical solution of the one-dimensional PBE is presented for two plates interacting in symmetrical electrolytes. The plates can have arbitrary surface potentials at infinite separation as long they have the same sign. Such a condition covers a majority of the colloidal systems encountered. The solution leads to a simple relationship which permits determination of surface potentials, surface charge densities, and electrostatic pressures as a function of plate separation H for different charging scenarios. An analytical expression is also presented for the potential profile between the plates for a given separation. Comparison of these potential profiles with those obtained by numerical analysis shows the validity of the proposed solution.     

A finite expansion method for the calculation and interpretation of molecular electrostatic potentials

Because it is useful to have the molecular electrostatic potential as an element in a complex scheme to assess the toxicity of large molecules, efficient and reliable methods are needed for the calculation and characterization of these potentials. A multicenter multipole expansion of the molecular electron charge density calculated with a limited Gaussian basis set is shown here to have only a finite number of nonzero terms from which the molecular electrostatic potential can be calculated. The discrete contributions to the electrostatic potentials from the terms of this expansion provide a physically meaningful decomposition of the potential and a means for its characterization. With pyrrole as an example, the electrostatic potential calculated from this finite expansion of the electron density is compared to that obtained from exact calculations from the same wave function. Good agreement is obtained at distances greater than 1.5 A from any atom in the molecule. In contrast, rearrangement of the terms into an expansion corresponding only to Mulliken atomic charges and dipoles yields a decomposition that produces electrostatic potentials which agree less well with the exact potential. This discrepancy is attributable to the neglect of terms due to higher moments.     

Asymmetry of current–voltage characteristics: a bilayer model of a modified ion-exchange membrane

The asymmetry of the current–voltage characteristics of ion-exchange membranes is explained in terms of the model of a bilayer fine porous membrane with constant charge distributions over the thickness of layers. This model has previously been proposed for determining diffusion permeability of membranes. In the case of one uncharged (neutral) layer, a set of two implicit algebraic equations is derived for determining the total current–voltage characteristics (CVC) of a membrane. For the first time, implicit algebraic equations are obtained for calculating the limiting currents at different orientations of an anisotropic membrane in an electrodialysis cell and explicit expressions are derived for determining specific conductivity of the membrane from the slope of the ohmic region of a CVC under the approximation of “excluded coions.” The model may be successfully used for describing the CVCs of perfluorinated MF-4SC sulfonic cation-exchange membranes, the surface layers of which are modified with polyaniline or halloysite.     

A unified analytical treatment of the acid-dissociation equilibria of weakly acidic linear polyelectrolytes and the conjugate acids of weakly basic linear polyelectrolytes

 The influence of added sodium chloride concentration levels on the acid-dissociation equilibria of a weakly acidic linear polyelectrolyte and a conjugate acid of weakly basic linear polyelectrolyte has been investigated potentiometrically by use of polyacrylic acid (PAA) and poly(N-vinylimidazole) (PVIm) as examples of polyelectrolytes. Both equilibria are strongly influenced by the degree of dissociation of the polyacids as well as the concentration levels of sodium chloride due to an electrostatic effect originating from the negatively or positively charged polymer surfaces. These have been analyzed in a unified manner by taking accounts of two-phase properties of the charged linear polyions. Distribution of counterions and coions between a polyelectrolyte phase formed around the polymer skeleton and a bulk solution phase has been rationalized by a Donnan’s relation. Introduction of a volume term for the polyelectrolyte phase permits definition of averaged concentrations of mobile ions in the vicinity of the polyion molecules, which enables us to define hypothetical intrinsic acid-dissociation constants in the polyion domain. The intrinsic constants estimated by extrapolation of apparent acid-dissociation constants at zero-charge state are in good agreement with the acid-dissociation constants of the monomer analogs of the polymers, i.e., acetic acid for PAA and imidazole for PVIm, respectively. The difference between the apparent and intrinsic acid-dissociation constants for PVIm was much higher than that for PAA at defined degree of dissociation of the polyacids, even though the separations of the functionalities fixed on the linear polymers are approximately equal to each other. Received: 4 February 1997 Accepted: 26 May 1997