Interactions between charged surfaces immersed in a polyelectrolyte solution

With grand canonical simulations invoking a configurationally weighted scheme, we have calculated interactions between charged surfaces immersed in a polyelectrolyte solution. In contrast to previous simulations of such systems, we have imposed full equilibrium conditions (i.e., we have included diffusive equilibrium with a bulk solution). This has a profound impact on the resulting interactions: even at modest surface charge densities, oppositely charged chains will, at sufficiently large separations, adsorb strongly enough to overcompensate for the nominal surface charge. This phenomenon, known as charge inversion, generates a double-layer repulsion and a free-energy barrier. Simpler canonical approaches, where the chains are assumed to neutralize the surfaces perfectly, will not capture this stabilizing barrier. The barrier height increases with the length of the polyions. Interestingly enough, the separation at which the repulsion becomes attractive is independent of chain length. The short chains here are unable to reach across from one surface to the other. We therefore conclude that the transition to an attractive regime is not provided by the formation of such "intersurface" bridges. With long chains and at large separations, charge inversion displays decaying oscillatory behavior (i.e., the apparent surface charge switches sign once again). This is due to polyion packing effects. We have also investigated responses to salt addition and changes in polyelectrolyte concentration. Our results are in qualitative and semiquantitative agreement with experimental findings, although it should be noted that our chains are comparatively short, and the experimental surface charge density is poorly established.     

Molecular dynamics simulations of polyelectrolyte adsorption

We have performed molecular dynamics simulations of polyelectrolyte adsorption at oppositely charged surfaces from dilute polyelectrolyte solutions. In our simulations, polyelectrolytes were modeled by chains of charged Lennard-Jones particles with explicit counterions. We have studied the effects of the surface charge density, surface charge distribution, solvent quality for the polymer backbone, strength of the short-range interactions between polymers and substrates on the polymer surface coverage, and the thickness of the adsorbed layer. The polymer surface coverage monotonically increases with increasing surface charge density for almost all studied systems except for the system of hydrophilic polyelectrolytes adsorbing at hydrophilic surfaces. In this case the polymer surface coverage saturates at high surface charge densities. This is due to additional monomer-monomer repulsion between adsorbed polymer chains, which becomes important in dense polymeric layers. These interactions also preclude surface overcharging by hydrophilic polyelectrolytes at high surface charge densities. The thickness of the adsorbed layer shows monotonic dependence on the surface charge density for the systems of hydrophobic polyelectrolytes for both hydrophobic and hydrophilic surfaces. Thickness is a decreasing function of the surface charge density in the case of hydrophilic surfaces while it increases with the surface charge density for hydrophobic substrates. Qualitatively different behavior is observed for the thickness of the adsorbed layer of hydrophilic polyelectrolytes at hydrophilic surfaces. In this case, thickness first decreases with increasing surface charge density, then it begins to increase.     

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)].     

Forces between interfaces in concentrated nanoparticle suspensions and polyelectrolyte solutions

This article provides an overview of interactions between charged interfaces across concentrated suspensions of charged nanoparticles or solutions of polyelectrolytes. These systems bear many similarities. We distinguish the like-charged and oppositely charged situations. In the like-charged situation, a layered structure in the proximity of the interface is formed. This structure induces a strongly repulsive energy profile at shorter distances, which originates from a gap that is free of nanoparticles or polyelectrolytes. At larger distances, the profile becomes oscillatory. This energy profile can be quantified with a simple model, which distinguishes the near-field region and the far-field region. The parameters entering the model show characteristic scaling relations. In the oppositely charged situation, a saturated, tightly bound layer at the interface forms. This layer leads to a charge reversal of the interface and induces a similar layered structure as in the like-charged case.     

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

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

The interplay of diffusional and electrophoretic transport mechanisms of charged solutes in the liquid film surrounding charged nonporous adsorbent particles employed in finite bath adsorption systems

A model that describes the diffusive and electrophoretic mass transport of the cation and anion species of a buffer electrolyte and of a charged adsorbate in the liquid film surrounding nonporous adsorbent particles in a finite bath adsorption system, in which adsorption of the charged adsorbate onto the charged surface of the nonporous particles occurs, is constructed and solved. The dynamic behavior of the mechanisms of this model explicitly demonstrates (a) the interplay between the diffusive and electrophoretic molar fluxes of the charged adsorbate and of the species of the buffer electrolyte in the liquid film surrounding the nonporous adsorbent particles, (b) the significant effect that the functioning of the electrical double layer has on the transport of the charged species and on the adsorption of the charged adsorbate, and (c) the substantial effect that the dynamic behavior of the surface charge density has on the functioning of the electrical double layer. It is found that at equilibrium, the value of the concentration of the charged adsorbate in the fluid layer adjacent to the surface of the adsorbent particles is significantly greater than the value of the concentration of the adsorbate in the finite bath, while, of course, the net molar flux of the charged adsorbate in the liquid film is equal to zero at equilibrium. This result is very different than that obtained from the conventional model that is currently used to describe the transport of a charged adsorbate in the liquid film for systems involving the adsorption of a charged adsorbate onto the charged surface of nonporous adsorbent particles; the conventional model (i) does not consider the existence of an electrical double layer, (ii) assumes that the transport of the charged adsorbate occurs only by diffusion in the liquid film, and (iii) causes at equilibrium the value of the charged adsorbate in the liquid layer adjacent to the surface of the particles to become equal to the value of the concentration of the charged adsorbate in the liquid of the finite bath. Furthermore, it was found that a maximum can occur in the dynamic behavior of the concentration of the adsorbate in the adsorbed phase when the value of the free molecular diffusion coefficient of the adsorbate is relatively large, because the increased magnitude of the synergistic interplay between the diffusive and electrophoretic molar fluxes of the adsorbate in the liquid film allows the adsorbate to accumulate (to be entrapped) in the liquid layer adjacent to the surface of the adsorbent particles faster than the concentrations of the electrolyte species, whose net molar fluxes are significantly hindered due to their opposing diffusive and electrophoretic molar fluxes, can adjust to account for the change in the surface charge density of the particles that arises from the adsorption of the charged adsorbate. The results presented in this work also have significant implications in finite bath adsorption systems involving the adsorption of a charged adsorbate onto the surface of the pores of charged porous adsorbent particles, because the diffusion and the electrophoretic migration of the charged solutes (cations, anions, and charged adsorbate) in the pores of the adsorbent particles will depend on the dynamic concentration profiles of the charged solutes in the liquid film surrounding the charged porous adsorbent particles. The results of the present work are also used to illustrate how the functioning of the electrical double layer could contribute to the development of inner radial humps (concentration rings) in the concentration of the adsorbate in the adsorbed phase of charged porous adsorbent particles.     

Diffusioosmosis and electroosmosis in a capillary slit with surface charge layers

This study analytically examines the steady diffusioosmotic and electroosmotic flows of an electrolyte solution in a fine capillary slit with each of its inside walls covered by a layer of adsorbed polyelectrolytes. In this solvent-permeable and ion-penetrable surface charge layer, idealized polyelectrolyte segments are assumed to distribute at a uniform density. The electric double layer and the surface charge layer may have arbitrary thicknesses relative to the gap width between the slit walls. The electrostatic potential distribution on a cross section of the slit is obtained by solving the linearized Poisson–Boltzmann equation, which applies to the case of low potentials or low fixed-charge densities. Explicit formulas for the fluid velocity profile due to the imposed electrolyte concentration gradient or electric field through the slit are derived as the solution of a modified Navier–Stokes/Brinkman equation. The results demonstrate that the structure of the surface charge layer can lead to an augmented or a diminished electrokinetic flow (even a reversal in direction of the flow) relative to that in a capillary with bare walls, depending on the characteristics of the capillary, of the surface charge layer, and of the electrolyte solution. For the diffusioosmotic flow with an induced electric field, competition between electroosmosis and chemiosmosis can result in more than one reversal in direction of the flow over a range of the Donnan potential of the adsorbed polyelectrolyte in the capillary.     

Cylindrical cell model for the electrostatic free energy of polyelectrolyte complexes

Associative phase separation (complex coacervation) in a mixture of oppositely charged polyelectrolytes can lead to different types of (inter-)polyelectrolyte complexes (soluble micelles, macroscopic precipitation). In a previous report [Langmuir 2004, 20, 2785-2791], we presented a model for the electrostatic free energy change when (weakly charged) polyelectrolyte forms a homogeneous complex phase. The influence of ionization of the polymer on the electrostatic free energy of the complex was incorporated but the influence of complex density neglected. In the present effort, cylindrical cells are assumed around each polyelectrolyte chain in the complex, and on the basis of the Poisson-Boltzmann equation, the electrostatic free energy is calculated as a function of the complex density. After combination with Flory-Huggins mixing free energy terms and minimization of the total free energy, the equilibrium complex density is obtained, for a given ratio of polycations to polyanions in the complex. The analysis is used in an example calculation ofpolyelectrolyte film formation by alternatingly applying a polycation and a polyanion solution. The calculation suggests that the often observed exponential growth of a polyelectrolyte film when the polymer is weakly charged has a thermodynamic origin: the polyelectrolyte complex shifts repeatedly between two equilibrium states of different densities and compositions. However, when the polyelectrolytes are strongly charged the difference in the compositions between the two equilibrium states is very small, and exponential growth by an absorption mechanism is no longer possible.     

Ionic exclusion phase transition in neutral and weakly charged cylindrical nanopores

A field theoretic variational approach is introduced to study ion penetration into water-filled cylindrical nanopores in equilibrium with a bulk reservoir [S. Buyukdagli, M. Manghi, and J. Palmeri, Phys. Rev. Lett. 105, 158103 (2010)]. It is shown that an ion located in a neutral pore undergoes two opposing mechanisms: (i) a deformation of its surrounding ionic cloud of opposite charge, with respect to the reservoir, which increases the surface tension and tends to exclude ions from the pore, and (ii) an attractive contribution to the ion self-energy due to the increased screening with ion penetration of the repulsive image forces associated with the dielectric jump between the solvent and the pore wall. For pore radii around 1 nm and bulk concentrations lower than 0.2 mol/l, this mechanism leads to a first-order phase transition, similar to capillary "evaporation," from an ionic-penetration state to an ionic-exclusion state. The discontinuous phase transition exists within the biological concentration range (～0.15 mol/l) for small enough membrane dielectric constants (ε(m) < 5). In the case of a weakly charged pore, counterion penetration exhibits a nonmonotonic behavior and is characterized by two regimes: at low reservoir concentrations or small pore radii, coions are excluded and counterions enter the pore to enforce electroneutrality; dielectric repulsion (image forces) remain strong and the counterion partition coefficient decreases with increasing reservoir concentration up to a characteristic value. For larger reservoir concentrations, image forces are screened and the partition coefficient of counterions increases with the reservoir concentration, as in the neutral pore case. Large surface charge densities (>2 × 10(-3) e/nm(2)) suppress the discontinuous transition by reducing the energy barrier for ion penetration and shifting the critical point toward very small pore sizes and reservoir concentrations. Our variational method is also compared to a previous self-consistent approach and yields important quantitative corrections. The role of the curvature of dielectric interfaces is highlighted by comparing ionic penetration into slit and cylindrical pores. Finally, a charge regulation model is introduced in order to explain the key effect of pH on ionic exclusion and explain the origin of observed time-dependent nanopore electric conductivity fluctuations and their correlation with those of the pore surface charge.     

Constant-NTμ simulations: Free energy difference method for excess adsorption

We describe a new method to simulate confined fluids in equilibrium with a bulk. The equilibrium is first located at a low density, where conventional methods relying upon real or virtual particle exchanges are reliable. Thereafter, the chemical potential of the fluid is increased by the same amount in both systems using a variation of a recently developed free energy difference method in an isobaric (isotension) ensemble. The method, illustrated here for a simple fluid in a planar slit, is shown to be reliable up to high densities. As the method does not rely upon particle exchanges it is ideally suited to the simulation of equilibria in complex fluids, e.g., molecular liquids and polymers. © 1996 by John Wiley & Sons, Inc.     

Molecular structure of single DNA complexes with positively charged dendronized polymers

Positively charged dendronized polymers with protonated amine groups at the periphery and different dendron generations are cylindrically shaped nanoobjects whose radii and linear charge densities can be varied systematically. These polyelectrolytes have been complexed with DNA and subsequently adsorbed on precoated mica substrates. The analysis of scanning force microscopy data indicates that DNA wraps around the dendronized polymers. The calculated pitch is 2.30 +/- 0.27 and 2.16 +/- 0.27 nm for DNA wrapped around dendronized polymers of generation two and four, respectively. The complex with the second generation has been shown to be negatively charged, which is consistent with the theory of spontaneous overcharging of macro-ion complexes, when the electrostatic contribution to the free energy dominates over the elastic energy. The complexes may be of interest for the development of nonviral gene delivery systems.     

Adsorption of heterogeneously charged nanoparticles on a variably charged surface by the extended surface complexation approach: charge regulation, chemical heterogeneity, and surface complexation

Adsorption of randomly branched polyelectrolytes, "hairy" particles and internally structured macromolecules, collectively denoted as heterogeneously charged nanoparticles, on charged surfaces is important in many technological and natural processes. In this paper, we will focus on (1) the charge regulation of both the nanoparticle and the surface and (2) the surface complexation between the particle functional groups and the surface sites and will theoretically study the adsorption using the extended surface complexation approach. The model explicitly considers the electrochemical potential of a nanoparticle with an average (smeared-out) structure and charge both in bulk solution and on the surface to obtain the equilibrium adsorption. The chemical heterogeneity of the particle is described by a distribution of the protonation constant. Detailed analysis of the chemical potential of the adsorbed nanoparticle reveals that the pH and salt dependence of the adsorption can be largely explained by the balance between an energy gain resulting from the particle and surface charge regulation and the surface complexation and an energy loss from the unfavorable interparticle electrostatic repulsion close to the surface. This conclusion is also supported by the strong impacts that the chemical heterogeneity of the particle functional groups, the magnitude of the surface complexation, the number of the functional groups, and the size of the particle have on the adsorption.     

Monte Carlo simulation and molecular theory of tethered polyelectrolytes

We investigate the structure of end-tethered polyelectrolytes using Monte Carlo simulations and molecular theory. In the Monte Carlo calculations we explicitly take into account counterions and polymer configurations and calculate electrostatic interaction using Ewald summation. Rosenbluth biasing, distance biasing, and the use of a lattice are all used to speed up Monte Carlo calculation, enabling the efficient simulation of the polyelectrolyte layer. The molecular theory explicitly incorporates the chain conformations and the possibility of counterion condensation. Using both Monte Carlo simulation and theory, we examine the effect of grafting density, surface charge density, charge strength, and polymer chain length on the distribution of the polyelectrolyte monomers and counterions. For all grafting densities examined, a sharp decrease in brush height is observed in the strongly charged regime using both Monte Carlo simulation and theory. The decrease in layer thickness is due to counterion condensation within the layer. The height of the polymer layer increases slightly upon charging the grafting surface. The molecular theory describes the structure of the polyelectrolyte layer well in all the different regimes that we have studied.     

Theory of repulsive charged colloids in slit-pores

Using classical density functional theory (DFT) we analyze the structure of the density profiles and solvation pressures of negatively charged colloids confined in slit pores. The considered model, which was already successfully employed to study a real colloidal (silica) suspension [S. H. L. Klapp et al., Phys. Rev. Lett. 100, 118303 (2008)], involves only the macroions which interact via the effective Derjaguin-Landau-Verwey-Overbeek (DLVO) potential supplemented by a hard core interaction. The solvent enters implicitly via the screening length of the DLVO interaction. The free energy functional describing the colloidal suspension consists of a hard sphere contribution obtained from fundamental measure theory and a long range contribution which is treated using two types of approximations. One of them is the mean field approximation (MFA) and the remaining is based on Rosenfeld's perturbative method for constructing the Helmholtz energy functional. These theoretical calculations are carried out at different bulk densities and wall separations to compare finally to grand canonical Monte Carlo simulations. We also consider the impact of charged walls. Our results show that the perturbative DFT method yields generally qualitatively consistent and, for some systems, also quantitatively reliable results. In MFA, on the other hand, the neglect of charge-induced correlations leads to a breakdown of this approach in a broad range of densities.     

Charged Semi-Permeable Shell with Encapsulated Polyions: Concentration Profile,Surface Potential,and Electrostatic Pressure

Summary: We study theoretically the electrostatic equilibrium for a charged shell filled with a suspension of polyions (e.g., colloids, polyelectrolytes, etc.) and immersed in an infinite salt-free reservoir. The shell is impermeable for polyions, but allows free diffusion of counterions. From the solution of the linearized Poisson-Boltzmann equation we obtain the distribution of the potential and concentration profiles for polyions. We then derive explicit formulas for the electrostatic pressure exerted by the shell. If the overall charge of the filled shell has the same sign as the surface alone the pressure on the shell increases with increase of the surface charge density. Otherwise the surface charge density suppresses the electro-osmotic pressure due to the electrostatic attraction between the oppositely charged polyions and shell.     

Electrokinetic flow in a capillary with a charge-regulating surface polymer layer

An analytical study of the steady electrokinetic flow in a long uniform capillary tube or slit is presented. The inside wall of the capillary is covered by a layer of adsorbed or covalently bound charge-regulating polymer in equilibrium with the ambient electrolyte solution. In this solvent-permeable and ion-penetrable surface polyelectrolyte layer, ionogenic functional groups and frictional segments are assumed to distribute at uniform densities. The electrical potential and space charge density distributions in the cross section of the capillary are obtained by solving the linearized Poisson-Boltzmann equation. The fluid velocity profile due to the application of an electric field and a pressure gradient through the capillary is obtained from the analytical solution of a modified Navier-Stokes/Brinkman equation. Explicit formulas for the electroosmotic velocity, the average fluid velocity and electric current density on the cross section, and the streaming potential in the capillary are also derived. The results demonstrate that the direction of the electroosmotic flow and the magnitudes of the fluid velocity and electric current density are dominated by the fixed charge density inside the surface polymer layer, which is determined by the regulation characteristics such as the dissociation equilibrium constants of the ionogenic functional groups in the surface layer and the concentration of the potential-determining ions in the bulk solution.     

Mixed systems of hydrophobically modified polyelectrolytes: controlling rheology by charge and hydrophobe stoichiometry and interaction strength

Rheology and phase separation were investigated for aqueous mixtures of two oppositely charged hydrophobically modified polyelectrolytes. The typical phase separation, normally seen for oppositely charged polymer mixtures, is dramatically reduced by the presence of hydrophobic modification, and phase separation is only detected close to the point of charge neutralization. While the two polyelectrolytes separately can give high viscosities and a gel-like behavior, a pronounced maximum in viscosity and storage modulus with the mixing ratio of the polyelectrolytes is observed; the maximum is located between the points of charge and hydrophobe stoichiometry and reflects a combination of hydrophobic and electrostatic association. Lowering the charge density of the anionic polymer leads to a strengthened association at first, but at lower charge densities there is a weakened association due to the onset of phase separation. The strength of the electrostatic interaction was modified by adding salt. Increased ionic strength can lead to phase separation and to increased or decreased viscosity depending on the polyelectrolyte mixing ratio.     

Efficiently accounting for ion correlations in electrokinetic nanofluidic devices using density functional theory

The electrokinetic behavior of nanofluidic devices is dominated by the electrical double layers at the device walls. Therefore, accurate, predictive models of double layers are essential for device design and optimization. In this paper, we demonstrate that density functional theory (DFT) of electrolytes is an accurate and computationally efficient method for computing finite ion size effects and the resulting ion-ion correlations that are neglected in classical double layer theories such as Poisson-Boltzmann. Because DFT is derived from liquid-theory thermodynamic principles, it is ideal for nanofluidic systems with small spatial dimensions, high surface charge densities, high ion concentrations, and/or large ions. Ion-ion correlations are expected to be important in these regimes, leading to nonlinear phenomena such as charge inversion, wherein more counterions adsorb at the wall than is necessary to neutralize its surface charge, leading to a second layer of co-ions. We show that DFT, unlike other theories that do not include ion-ion correlations, can predict charge inversion and other nonlinear phenomena that lead to qualitatively different current densities and ion velocities for both pressure-driven and electro-osmotic flows. We therefore propose that DFT can be a valuable modeling and design tool for nanofluidic devices as they become smaller and more highly charged.     

Equilibrium domains on heterogeneously charged surfaces

Recent experiments have shown that salt solutions containing surfaces with two oppositely charged species show stable, possibly equilibrium, structures with finite domain sizes. The short-range interactions between the two species would normally result in phase separation that is driven by the line tension with macroscopically large domains of each species. In this paper, we show that, when at least one of the charged species is mobile, finite domains can occur in equilibrium. The domain size is determined by a competition of the electrostatic free energy that promotes charge mixing and small domains, with the line tension that promotes macroscopic phase separation. We calculate the equilibrium patch size as a function of the surface charge and the concentration of dissolved monovalent salts in the bulk phase. An important finding is the prediction of a first-order transition from finite patches to macroscopic phase separation of the two charge species as the salt concentration is increased.     

Dilational surface visco-elasticity of polyelectrolyte/surfactant solutions: formation of heterogeneous adsorption layers

Recent application of the methods of surface dilational rheology to solutions of the complexes between synthetic polyelectrolytes and oppositely charged surfactants (PSC) gave a possibility to determine some steps of the adsorption layer formation and to discover an abrupt transition connected with the formation of microaggregates at the liquid surface. The kinetic dependencies of the dynamic surface elasticity are always monotonous at low surfactant concentrations but can have one or two local maxima in the range beyond the critical aggregation concentration. The first maximum is accompanied by the generation of higher harmonics of induced surface tension oscillations and caused by heterogeneities in the adsorption layer. The formation of a multilayered structure at the surface for some systems leads to the second maximum in the dynamic surface elasticity. The hydrophobicity and charge density of a polymer chain influence strongly the surface structure, resulting in a variety of dynamic surface properties of PSC solutions. Optical methods and atomic force microscopy give additional information for the systems under consideration. Experimental results and existing theoretical frameworks are reviewed with emphasis on the general features of all studied PSC systems.