Poisson-Boltzmann theory of the charge-induced adsorption of semi-flexible polyelectrolytes

A model is suggested for the structure of an adsorbed layer of a highly charged semi-flexible polyelectrolyte on a weakly charged surface of opposite charge sign. The adsorbed phase is thin, owing to the effective reversal of the charge sign of the surface upon adsorption, and ordered, owing to the high surface density of polyelectrolyte strands caused by the generally strong binding between polyelectrolyte and surface. The Poisson-Boltzmann equation for the electrostatic interaction between the array of adsorbed polyelectrolytes and the charged surface is solved for a cylindrical geometry, both numerically, using a finite element method, and analytically within the weak curvature limit under the assumption of excess monovalent salt. For small separations, repulsive surface polarization and counterion osmotic pressure effects dominate over the electrostatic attraction and the resulting electrostatic interaction curve shows a minimum at nonzero separations on the Angstrom scale. The equilibrium density of the adsorbed phase is obtained by minimizing the total free energy under the condition of equality of chemical potential and osmotic pressure of the polyelectrolyte in solution and in the adsorbed phase. For a wide range of ionic conditions and charge densities of the charged surface, the interstrand separation as predicted by the Poisson-Boltzmann model and the analytical theory closely agree. For low to moderate charge densities of the adsorbing surface, the interstrand spacing decreases as a function of the charge density of the charged surface. Above about 0.1 M excess monovalent salt, it is only weakly dependent on the ionic strength. At high charge densities of the adsorbing surface, the interstrand spacing increases with increasing ionic strength, in line with the experiments by Fang and Yang [J. Phys. Chem. B 101, 441 (1997)].     

Interaction between brush layers of bottle-brush polyelectrolytes: molecular dynamics simulations

Interactions between tethered layers composed of aggrecan (charged bottle-brush) macromolecules are responsible for the molecular origin of cartilage biomechanical behavior. To elucidate the role of the electrostatic forces in interaction between bottle-brush layers, we have performed molecular dynamics simulations of charged and neutral bottle-brush macromolecules tethered to substrates. In the case of charged bottle-brush layers, the disjoining pressure P between two brush layers in salt-free solutions increases with decreasing distance D between substrates as P ∝ D(-1.8). A stronger dependence of the disjoining pressure P on the surface separation D was observed for neutral bottle-brushes, P ∝ D(-4.6), in the same interval of disjoining pressures. These scaling laws for dependence of disjoining pressure P on distance D are due to bending energy of the bottle-brush macromolecules within compressed brush layers. The weaker distance dependence observed in polyelectrolyte bottle-brushes is due to interaction between counterion clouds surrounding the bottle-brush macromolecules preventing strong brush overlap.     

Molecular dynamics simulation of discontinuous volume phase transitions in highly-charged crosslinked polyelectrolyte networks with explicit counterions in good solvent

The volumetric properties of highly-charged defect-free polyelectrolyte networks with tetrafunctional crosslinks are studied through molecular dynamics simulations in the canonical ensemble. The network backbone monomers, which are monovalent, and the counterions, which are mono-, di-, or trivalent, are modeled explicitly in the simulations, but the solvent is treated implicitly as a dielectric medium of good solvation quality. The osmotic pressure of the network-solvent system is found to depend greatly on the strength of electrostatic interactions. Discontinuous volume phase transitions are observed when the electrostatic interactions are strong, and the onset of these transitions shifts to higher solvent dielectricity as the counterion valency increases. The roles of the various virial contributions to the osmotic pressure are examined. The network elasticity entropy is found to behave nearly classically. As the network contracts and collapses with increasing strength of electrostatic interactions, the loss of counterion entropy leads to increased counterion osmotic pressure contributions via two mechanisms. The reduction in available configurational space increases the counterion translational entropy contribution to the ideal part of the osmotic pressure, and the greater number of counterion-monomer contacts formed due to counterion condensation and confinement increases the counterion excluded-volume entropy contribution to the excess part of the osmotic pressure. These observations contrast the decrease in the single ideal-gas-like counterion translational entropy contribution to the osmotic pressure predicted by the counterion condensation-charge renormalization theory. An accompanying decrease in the total electrostatic energy balances the loss of counterion excluded-volume entropy as the polyelectrolyte networks collapse in low-dielectric solvents. This interplay between the electrostatic energy and the counterion excluded-volume entropy appears to be responsible for the discontinuous volume phase transitions that are observed in polyelectrolyte networks. The structure of the polyelectrolyte network is also found to be affine in the swollen state, with constituent chains nearly fully extended, and nonaffine in the collapsed state, with the chains adopting a Gaussian conformation.     

Ion induced lamellar-lamellar phase transition in charged surfactant systems

We propose a model for the liquid-liquid (L(alpha)-->L(alpha(') )) phase transition observed in osmotic pressure measurements of certain charged lamellae-forming amphiphiles. The model free energy combines mean-field electrostatic and phenomenological nonelectrostatic interactions, while the number of dissociated counterions is treated as a variable degree of freedom that is determined self-consistently. The model, therefore, joins two well-known theories: the Poisson-Boltzmann theory for ionic solutions between charged lamellae and the Langmuir-Frumkin-Davies adsorption isotherm modified to account for charged adsorbing species. Minimizing the appropriate free energy for each interlamellar spacing, we find the ionic density profiles and the resulting osmotic pressure. While in the simple Poisson-Boltzmann theory the osmotic pressure isotherms are always smooth, we observe a discontinuous liquid-liquid phase transition when the Poisson-Boltzmann theory is self-consistently augmented by the Langmuir-Frumkin-Davies adsorption. This phase transition depends on the area per amphiphilic head group, as well as on nonelectrostatic interactions of the counterions with the lamellae and interactions between counterion-bound and counterion-dissociated surfactants. Coupling the lateral phase transition in the bilayer plane with electrostatic interactions in the bulk, our results offer a qualitative explanation for the existence of the L(alpha)-->L(alpha(') ) phase transition of didodecyldimethylammonium bromide (DDABr), but the transition's apparent absence for the chloride and the iodide homologs. More quantitative comparisons with experiment require better understanding of the microscopic basis of the phenomenological model parameters.     

The bridging force between two plates by polyelectrolyte chains

The interaction between particles in a colloidal system can be significantly affected by their bridging by polyelectrolyte chains. In this paper, the bridging is investigated by using a self-consistent field approach which takes into account the van der Waals interactions between the segments of the polyelectrolyte molecules and the plates, as well as the electrostatic and volume exclusion interactions. A positive contribution to the force between two plates is generated by the van der Waals interactions between the segments and the plates. This positive (repulsive) contribution plays an important role in the force when the distances between the plates are small. With increasing van der Waals interaction strength between segments and plates, the force between the plates becomes more repulsive at small distances and more attractive at large distances. When the surfaces of the plates have a constant surface electrical potential and a charge sign opposite to that of the polyelectrolyte chains, the force between the two plates becomes less attractive as the bulk polyelectrolyte concentration increases. This behavior is due to a higher bulk counterion concentration dissociated from the polyelectrolyte molecules. At short distances, the force between plates is more repulsive for stiffer chains. A comparison between theoretical and experimental results regarding the contraction of the interlayer separation between the platelets of vermiculite clays against the concentration of poly(vinyl methyl ether) was made.     

Electrostatic Repulsion in Concentrated Disperse Systems

Electrostatic interactions are considered in the framework of the cell model to predict the osmotic pressure in concentrated disperse systems. A procedure was developed to represent the osmotic pressure as a function of two parameters, namely, the dispersed phase volume fraction and the electric potential attributed to the interface between the continuous and dispersed phases. The procedure is based on a general formula which was derived to express the electrostatic contribution to the osmotic pressure through the electric potential at the cell boundary. The potential of the cell boundary is predicted from the solution of the Poisson-Boltzmann problem which was specified for the cell model approach. The Poisson-Boltzmann problem is solved by a perturbation technique using a normalized interface potential as the perturbation parameter. Three leading terms were obtained in the expansion of the osmotic pressure in terms of the normalized interface potential. Two options for the formation of the interface electric potential are discussed in the analysis of the interface potential dependency on the volume fraction of the dispersed phase. The first one is associated with the difference between the individual ionic distribution coefficients characterizing the equilibrium ratio between the concentrations in the bulk of the constituent phases. The second one deals with preferential adsorption of the carriers having a given electric charge sign. The dependency of the osmotic pressure on the system parameters is discussed and interrelated with other relevant theories. Special discussion is presented concerning the theory's application for the study of hydrocarbon disperse systems, e.g., water-in-oil emulsions. Copyright 2001 Academic Press.     

Aqueous suspensions of natural swelling clay minerals. 1. Structure and electrostatic interactions

In this article, we present a general overview of the organization of colloidal charged clay particles in aqueous suspension by studying different natural samples with different structural charges and charge locations. Small-angle X-ray scattering experiments (SAXS) are first used to derive swelling laws that demonstrate the almost perfect exfoliation of clay sheets in suspension. Using a simple approach based on geometrical constraints, we show that these swelling laws can be fully modeled on the basis of morphological parameters only. The validity of this approach was further extended to other clay data from the literature, in particular, synthetic Laponite. For all of the investigated samples, experimental osmotic pressures can be properly described by a Poisson-Boltzmann approach for ionic strength up to 10(-3) M, which reveals that these systems are dominated by repulsive electrostatic interactions. However, a detailed analysis of the Poisson-Boltzmann treatment shows differences in the repulsive potential strength that are not directly linked to the structural charge of the minerals but rather to the charge location in the structure for tetrahedrally charged clays (beidellite and nontronites) undergoing stronger electrostatic repulsions than octahedrally charged samples (montmorillonites, laponite). Only minerals subjected to the strongest electrostatic repulsions present a true isotropic to nematic phase transition in their phase diagrams. The influence of ionic repulsions on the local order of clay platelets was then analyzed through a detailed investigation of the structure factors of the various clay samples. It appears that stronger electrostatic repulsions improve the liquidlike positional local order.     

Brownian dynamics simulations of polyelectrolyte adsorption onto topographically patterned surfaces

The effect of patterned surface topography on the adsorption of single polyelectrolyte molecules is explored using Brownian dynamics simulations. The polyelectrolyte is modeled as a free-draining, freely jointed bead-rod chain, and electrostatic interactions are incorporated using a screened Coulombic potential with excluded volume interactions accounted for by the repulsive part of a Lennard-Jones potential. Topography consisting of periodically spaced valleys of square cross section separated by flat hills is considered. Chain conformations are characterized for a wide range of valley widths, depths, and spacings as well as for several different types of surface charge distributions. Depending on the parameter values describing the topography, the chains are found to adopt conformations ranging from flat and extended to those associated with bridge-, brush-, or semi-bridge-like structures. The formation of these structures is rationalized on the basis of a free-energy model that takes into account the increase in free energy due to entropic confinement, excluded volume interactions, and chain stretching as well as the decrease in free energy due to bead-surface electrostatic attraction. The results of this work are expected to be useful in designing patterned surface topography to control the conformations of adsorbed polyelectrolyte molecules.     

Loose complexation of weakly charged microemulsion droplets and a polyelectrolyte

The effect of polyelectrolyte addition on the properties of an oil-in-water (O/W) microemulsion of weakly charged spherical micelles is studied. The 81 A radius O/W droplets in this system can be charged by the partial substitution of the nonionic surfactant by a cationic surfactant. The effect of the addition of poly(acrylic acid) (PAA), which is a charged pH-dependent polyelectrolyte, on the interactions between charged or noncharged droplets has been investigated using SANS. We have characterized the phase behavior of this pH-smart system as a function of the microemulsion and the polyelectrolyte concentration and the number of charges per droplet at three pH values: pH = 2, 4.5, and 12. In particular, an associative phase separation due to the bridging of the droplets by the neutral PAA chains through H-bonds is observed with extremely low PAA addition at low pH. At the opposite, an addition of PAA at pH = 4.5 generates a strong repulsive contribution between neutral droplets. Electrostatic bonds between charged droplets and PAA, controlled by the number of charges per droplet, are responsible for a pH drift and then for an associative phase separation similar to that observed at low pH. Finally, at high pH, the creation of electrostatic bonds between fully charged PAA and charged droplets liberates sufficiently counterions in solution at high droplet charge density to screen the electrostatic interactions and to allow an associative phase separation.     

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.     

Confinement free energy of flexible polyelectrolytes in spherical cavities

A weakly charged flexible polyelectrolyte chain in a neutral spherical cavity is analyzed by using self-consistent field theory within an explicit solvent model. Assuming the radial symmetry for the system, it is found that the confinement of the chain leads to creation of a charge density wave along with the development of a potential difference across the center of cavity and the surface. We show that the solvent entropy plays an important role in the free energy of the confined system. For a given radius of the spherical cavity and fixed charge density along the backbone of the chain, solvent and small ion entropies dominate over all other contributions when chain lengths are small. However, with the increase in chain length, chain conformational entropy and polymer-solvent interaction energy also become important. Our calculations reveal that energy due to electrostatic interactions plays a minor role in the free energy. Furthermore, we show that the total free energy under spherical confinement is not extensive in the number of monomers. Results for the osmotic pressure and mean activity coefficient for monovalent salt are presented. We demonstrate that fluctuations at one-loop level lower the free energy and corrections to the osmotic pressure and mean activity coefficient of the salt are discussed. Finite size corrections are shown to widen the range of validity of the fluctuation analysis.     

Correlation between conformation change of polyelectrolyte brushes and lubrication

We directly monitor the absolute separation profiles that function as film thickness between a single glass disk and the charged polyelectrolyte brushes decorated steel slider in water using a home-made slider-on-disk apparatus, which reflects the structural conformation variations and interactions of polymer brushes under externally applied pressure, in addition to probing the relative variation of friction forces under different applied loads and sliding velocities. We find that the polyelectrolyte brushes modified surfaces can sustain high pressure and have extremely low friction coefficients(around 0.006 at pressures of 0.13 MPa; 0.5-0.6 without brushes). The water-lubrication characteristics are correlated to the structural conformation changes of the polyelectrolyte brushes that are mainly governed by electrostatic interactions and the osmotic pressure of counterions inside the polymer chains, which can be used to support and distribute the normal pressure. The apparent thickness of the brush decreases with the increase of loading forces, an increase in the ionic strength causes the polymer chains collapse, and the friction forces increase. This fundamental research is of great importance to understand the mechanical and structural properties of polyelectrolyte brushes and their influences on the tribological behaviors, and helps to design friction/lubrication-controlled surface/interface by taking advantage of polyelectrolyte brushes.     

Density functional study on the structural and thermodynamic properties of aqueous DNA-electrolyte solution in the framework of cell model

A density functional theory (DFT) in the framework of cell model is proposed to calculate the structural and thermodynamic properties of aqueous DNA-electrolyte solution with finite DNA concentrations. The hard-sphere contribution to the excess Helmholtz energy functional is derived from the modified fundamental measure theory, and the electrostatic interaction is evaluated through a quadratic functional Taylor expansion around a uniform fluid. The electroneutrality in the cell leads to a variational equation with a constraint. Since the reference fluid is selected to be a bulk phase, the Lagrange multiplier proves to be the potential drop across the cell boundary (Donnan potential). The ion profiles and electrostatic potential profiles in the cell are calculated from the present DFT-cell model. Our DFT-cell model gives better prediction of ion profiles than the Poisson-Boltzmann (PB)- or modified PB-cell models when compared to the molecular simulation data. The effects of polyelectrolyte concentration, ion size, and added-salt concentration on the electrostatic potential difference between the DNA surface and the cell boundary are investigated. The expression of osmotic coefficient is derived from the general formula of grand potential. The osmotic coefficients predicted by the DFT are lower than the PB results and are closer to the simulation results and experimental data.     

Brownian dynamics simulations of polyelectrolyte adsorption in shear flow with hydrodynamic interaction

The adsorption of single polyelectrolyte molecules in shear flow is studied using Brownian dynamics simulations with hydrodynamic interaction (HI). Simulations are performed with bead-rod and bead-spring chains, and electrostatic interactions are incorporated through a screened Coulombic potential with excluded volume accounted for by the repulsive part of a Lennard-Jones potential. A correction to the Rotne-Prager-Yamakawa tensor is derived that accounts for the presence of a planar wall. The simulations show that migration away from an uncharged wall, which is due to bead-wall HI, is enhanced by increases in the strength of flow and intrachain electrostatic repulsion, consistent with kinetic theory predictions. When the wall and polyelectrolyte are oppositely charged, chain behavior depends on the strength of electrostatic screening. For strong screening, chains get depleted from a region close to the wall and the thickness of this depletion layer scales as N(1/3)Wi(2/3) at high Wi, where N is the chain length and Wi is the Weissenberg number. At intermediate screening, bead-wall electrostatic attraction competes with bead-wall HI, and it is found that there is a critical Weissenberg number for desorption which scales as N(-1/2)kappa(-3)(l(B)|sigmaq|)(3/2), where kappa is the inverse screening length, l(B) is the Bjerrum length, sigma is the surface charge density, and q is the bead charge. When the screening is weak, adsorbed chains are observed to align in the vorticity direction at low shear rates due to the effects of repulsive intramolecular interactions. At higher shear rates, the chains align in the flow direction. The simulation method and results of this work are expected to be useful for a number of applications in biophysics and materials science in which polyelectrolyte adsorption plays a key role.     

Maximum disjoining pressure in protein stabilized concentrated oil-in-water emulsions

A model for the prediction of the equilibrium profile of film thickness and continuous phase liquid holdup profile in a concentrated oil-in-water (O/W) emulsion is proposed. This model is employed to infer the maximum disjoining pressure in a concentrated corn oil-in-water emulsion stabilized by bovine serum albumin (BSA) from the experimental measurements of different proportions of oil, polyhedral O/W foam, and aqueous layers at different centrifugal accelerations. The inferred maximum disjoining pressures were found to be higher at higher concentrations of BSA, lower ionic strengths as well as at pH values farther away from pI. The predicted variations of disjoining pressure with film thickness for a concentrated O/W emulsion stabilized by BSA exhibited two maxima due to steric and electrostatic interactions, respectively. The experimental maximum disjoining pressures for toluene-in-water emulsion stabilized by BSA were found to be about two to three times the predicted maxima due to steric interactions but were two to three orders of magnitude higher than the maxima due to electrostatic interactions, thus indicating that steric interaction is the dominant stabilizing mechanism. The discrepancy between the experimental and predicted maximum disjoining pressures is believed to be mainly due to lack of information with regard to the thickness of the adsorbed protein layer at the oil—water interface.     

Modeling Disjoining Pressures in Submicrometer Liquid-Filled Cylindrical Geometries

This work develops models for calculating the disjoining pressures of a cylindrical fluid "plug", specifically in submicrometer cylindrical pores. This modeling produces closed-form, cylindrical-pore disjoining pressures for London/van der Waals and solute/pore-wall adsorption interactions, which are the slit-pore models with the characteristic pore size replaced by the radius and multiplied by 6, resulting in a 48-fold or more increase in magnitude. In addition, this work contains a numerical solution for electrostatic interactions. The result of the numerical solution was a 9-fold increase in the modeled disjoining pressure compared to that in the slit-pore model. The cylindrical models may apply to the chemical coating of the interior walls of cylindrical pores or to the thermodynamics within droplets after the breakup of a fluid coating a surface. However, the application used as the base case in this paper is the extension of transport and thermodynamic laws for porous media, previously developed with capillary pressure models, to fully saturated porous media with submicrometer-sized pores. As such, the models could apply to mass transport in ultrafiltration, nanofiltration, and reverse-osmosis membranes. Copyright 2001 Academic Press.     

Synthesis and properties in solution of rodlike polyelectrolytes

Efficient synthetic strategies are described for the preparation of rodlike polyelectrolytes based on the intrinsically rigid poly(p-phenylene). Uncharged precursors were first prepared via the Suzuki coupling and then characterized by different methods of polymer analysis. Finally, they were transformed into polyelectrolytes using macromolecular substitution reactions. Depending on the substitution pattern, the obtained polyelectrolytes are either soluble or insoluble in water. Using water-soluble derivatives, the Poisson-Boltzmann cell model was tested by osmotic measurements and small-angle X-ray scattering. It is shown that the cell model provides a good first approximation of the distribution of the counterions around the macroion but still underestimates their correlation. Moreover, the PPP polyelectrolytes show a very pronounced polyelectrolyte effect. Since the rodlike PPPs are very rigid in shape, this observation proves that the polyelectrolyte effect is caused by long-range intermolecular electrostatic repulsion of the dissolved macroions rather than due to conformational changes.     

Counterion localization in solutions of starlike polyelectrolytes and colloidal polyelectrolyte brushes: a self-consistent field theory

A quantitative analysis of the distribution of counterions in salt-free solutions of colloidal polyelectrolyte brushes and starlike polyelectrolytes is performed on the level of the Poisson-Boltzmann approximation. Exact numerical solutions are obtained for starlike polyelectrolyte molecules composed of f = 20, . . ., 50 arms with a fixed fractional charge alpha per segment by applying the self-consistent field method of Scheutjens and Fleer (SF-SCF). The Wigner-Seitz cell dimension defines the concentration of polyelectrolyte stars in the system. The numerical results are compared to predictions of an analytical mean field theory and related to experimental observations on the osmotic pressure in solutions of starlike polyelectrolytes and colloidal polyelectrolyte brushes.