Counterion release from adsorbed highly charged polyelectrolyte: an electrooptical study

The adsorption of sodium poly(4-styrene sulfonate) on oppositely charged beta-FeOOH particles is studied by electrooptics. The focus of this paper is on the release of condensed counterions from adsorbed polyelectrolyte upon surface charge overcompensation. The fraction of condensed Na+ counterions on the adsorbed polyion surface is estimated according to the theory of Sens and Joanny and it is compared with the fraction of condensed counterions on nonadsorbed polyelectrolyte. The relaxation frequency of the electrooptical effect from the polymer-coated particle is found to depend on the polyelectrolyte molecular weight. This is attributed to polarization of the layer from condensed counterions on the polyion surface, being responsible for creation of the effect from particles covered with highly charged polyelectrolyte. The number of the adsorbed chains is calculated also assuming counterion condensation on the adsorbed polyelectrolyte and semiquantative agreement is found with the result obtained from the condensed counterion polarizability of the polymer-coated particle. Our findings are in line with theoretical predictions that the fraction of condensed counterions remains unchanged due to the adsorption of highly charged polyelectrolyte onto weakly charged substrate.     

Entropy and enthalpy of polyelectrolyte complexation: Langevin dynamics simulations

We report a systematic study by Langevin dynamics simulation on the energetics of complexation between two oppositely charged polyelectrolytes of same charge density in dilute solutions of a good solvent with counterions and salt ions explicitly included. The enthalpy of polyelectrolyte complexation is quantified by comparisons of the Coulomb energy before and after complexation. The entropy of polyelectrolyte complexation is determined directly from simulations and compared with that from a mean-field lattice model explicitly accounting for counterion adsorption. At weak Coulomb interaction strengths, e.g., in solvents of high dielectric constant or with weakly charged polyelectrolytes, complexation is driven by a negative enthalpy due to electrostatic attraction between two oppositely charged chains, with counterion release entropy playing only a subsidiary role. In the strong interaction regime, complexation is driven by a large counterion release entropy and opposed by a positive enthalpy change. The addition of salt reduces the enthalpy of polyelectrolyte complexation by screening electrostatic interaction at all Coulomb interaction strengths. The counterion release entropy also decreases in the presence of salt, but the reduction only becomes significant at higher Coulomb interaction strengths. More significantly, in the range of Coulomb interaction strengths appropriate for highly charged polymers in aqueous solutions, complexation enthalpy depends weakly on salt concentration and counterion release entropy exhibits a large variation as a function of salt concentration. Our study quantitatively establishes that polyelectrolyte complexation in highly charged Coulomb systems is of entropic origin.     

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

Effective charges of polyelectrolytes in a salt-free solution based on counterion chemical potential

The phenomenon of counterion condensation around a flexible polyelectrolyte chain with N monomers is investigated by Monte Carlo simulations in terms of the degree of ionization alpha, which is proportional to the effective charge. It is operationally defined as the ratio of observed to intrinsic counterion concentration, alpha = co/ci. The observed counterion concentration in the dilute polyelectrolyte solution is equivalent to an electrolyte solution of concentration co with the same counterion chemical potential. It can be determined directly by thermodynamic experiments such as ion-selective electrode. With the polyelectrolyte fixed at the center of the spherical Wigner-Seitz cell, the polymer conformation, counterion distribution, and chemical potential can be obtained. Our simulation shows that the degree of ionization rises as the polymer concentration decreases. This behavior is opposite to that calculated from the infinitely long charged rod model, which is often used to study counterion condensation. Moreover, we find that, for a specified line charge density, alpha decreases with an increment in chain length and chain flexibility. In fact, the degree of ionization is found to decline with increasing polymer fractal dimension, which can be tuned by varying bending modulus and solvent quality. Those results can be qualitatively explained by a simple model of two-phase approximation.     

Potential distribution around a polyelectrolyte-coated spherical particle in a salt-free medium

Simple analytic approximate expressions for the solution to the Poisson-Boltzmann equation around a spherical particle coated with an ion-penetrable polyelectrolyte layer in a salt-free medium containing counterions only are derived. The results of the calculation of the potential distribution using the approximate solution are found to be in good agreement with exact numerical results. It is shown that as in the case of a charged rigid particle, there is a certain critical value of the particle charge, separating two cases, that is, the low-particle-charge case and the high-particle-charge case. In the low-charge case the potential is essentially the same as if counterions were absent and thus the potential is proportional to the particle charge. In the high-charge case counterion condensation occurs in the polyelectrolyte layer region, so that the dependence of the potential on the particle charge is considerably suppressed.     

Molecular dynamics and neutron scattering study of the dependence of polyelectrolyte dendrimer conformation on counterion behavior

Atomistic molecular dynamics (MD) simulations and contrast variation small angle neutron scattering (SANS) have been combined to investigate the Generation-5 polyelectrolyte polyamidoamine starburst dendrimer. This work reveals the dendrimer conformational dependence on counterion association at different levels of molecular charge. The accuracy of the simulations is verified through satisfactory comparison between modeled results, such as excess intra-dendrimer scattering length density distribution and hydration level, and their experimental counterparts. While the counterion distributions are not directly measureable with SANS, the spatial distribution of the counterions and their dendrimer association are extracted from the validated MD equilibrium trajectories. It is found that the conformation of the charged dendrimer is strongly dependent on the counterion association. Sensitivity of the distribution of counterions around charged amines to the counterion valency is qualitatively explained by adopting Langmuir adsorption theory. Moreover, via extending the concept of electrical double layer for compact charged colloids, we define an effective radius of a charged dendrimer including the spatial distribution of counterions in its vicinity. Within the same framework, the correlation between the strength of intra-dendrimer electrostatic repulsion and the counterion valency and dynamics is also addressed.     

Electrophoretic evidence for a new type of counterion condensation

A systematic capillary electrophoresis study uncovered how polyelectrolyte effective charge density varies with backbone charge spacing and solvent dielectric constant. The study primarily focused on aliphatic ionenes, a special class of polyelectrolytes, which possess regularly spaced quaternary ammonium groups in the backbone. Complete ionization of functional units and good solvency in water or mixtures of water with lower dielectric constant solvents (methanol, acetonitrile) enabled continuous measurements of ionene effective charge density through the onset of counterion condensation. Ionenes with both uniform and alternating charge spacing were examined. As expected, effective charge density rose linearly with fixed charge density to a critical value, above which effective charge density remained constant. Deviating from expectation, the onset of condensation did not occur at a critical fixed charge density. Instead, condensation initiated at a constant critical Bjerrum length. The same onset condition was found for quaternized poly(vinyl pyridine)s. These experimental results suggest a new form of condensation, one driven by ion-pairing of polyelectrolyte with counterions. In support of this hypothesis, the onset of condensation appeared to correlate with counterion size. © 2004 Wiley Periodicals, Inc. J Polym Sci Part B: Polym Phys 42: 3616–3627, 2004     

Counterion condensation on charged spheres, cylinders, and planes

We use the framework of counterion condensation theory, in which deviations from linear electrostatics are ascribed to charge renormalization caused by collapse of counterions from the ion atmosphere, to explore the possibility of condensation on charged spheres, cylinders, and planes immersed in dilute solutions of simple salt. In the limit of zero concentration of salt, we obtain Zimm-Le Bret behavior: a sphere condenses none of its counterions regardless of surface charge density, a cylinder with charge density above a threshold value condenses a fraction of its counterions, and a plane of any charge density condenses all of its counterions. The response in dilute but nonzero salt concentrations is different. Spheres, cylinders, and planes all exhibit critical surface charge densities separating a regime of counterion condensation from states with no condensed counterions. The critical charge densities depend on salt concentration, except for the case of a thin cylinder, which exhibits the invariant criticality familiar from polyelectrolyte theory.     

Theory of polyelectrolyte adsorption onto surfaces patterned with charge and topography

Mean-field theory is used to derive criteria for the adsorption of a weakly charged polyelectrolyte molecule from salt solution onto surfaces patterned with charge and topography. For flat surfaces patterned with periodic arrays of charged patches, the adsorbed layer thickness predicted using mean-field theory and that found by Brownian dynamics simulations are in quantitative agreement in the strong-adsorption regime, which corresponds to sufficiently small kappa or sufficiently large |sigma(eff)q|, where kappa is the inverse Debye screening length, sigma(eff) is an effective surface charge density, and q is the charge on each segment of the polyelectrolyte. Qualitative agreement is obtained in the weak-adsorption regime, and for the case where surfaces are patterned with both charge and topography. For uniformly charged, sinusoidally corrugated surfaces, the theory predicts that the critical temperature required for adsorption can be greater than or less than the corresponding value for a flat surface depending on the relative values of kappa and the corrugation wave number. If the surface charge is also allowed to vary sinusoidally, then adsorption is predicted to occur only when the topography crests have a surface charge opposite to that of the polyelectrolyte. Surfaces patterned with rectangular indentations having charged bottoms which are separated by flat charged plateaus are investigated as well. Adsorption is predicted to occur even when the net surface charge is zero, provided that the plateaus have a charge opposite to that of the polyelectrolyte. If the charge on the plateaus and polyelectrolyte is the same, adsorption may still occur if electrostatic attraction from the indentation bottoms is sufficiently strong.     

The pH inside a swollen polyelectrolyte gel: poly(N-vinylimidazole)

The pH inside a dissolved polyelectrolyte coil or a swollen ionic polymer network is not accessible to direct measurement. It is here calculated through a simple model, based on Donnan equilibrium, counterion condensation (for charge density exceeding the critical value), and balance of mobile ions, without any assumption on the pKa of the ionizable groups. The data needed for the calculation with this model are polymer concentration, pH value in the initial solution, and pH value in the bath at equilibrium. All three can be determined experimentally by a batch method where the polymer is immersed in a different pot for each starting pH. The model is applied to a sample system, namely, chemically cross-linked poly(N-vinylimidazole) immersed in acidic baths of different pH values. The imidazole units are basic and become protonated by the acid, thus changing the pH of the initial bath. The model shows how the pH developed inside the swollen gel is several units higher than the pH of the bath at equilibrium, both with or without the correction for counterion condensation. Consequently, when the pKa of the polyelectrolyte is determined in the usual way (with the pH measured in the external bath), it gives an apparent value that is several units below the pKa determined from the actual pH inside the swollen gel at equilibrium. The inclusion of the counterion condensation decreases very slightly the polymer basicity. Surface effects and intramolecular association between protonated and unprotonated imidazole rings are discussed to explain the pKa behavior in the limit of low degree of ionization.     

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.     

Influence of the counterion size on swelling and collapse of polyelectrolyte gel

A theory that predicts the effect of the counterion size on the swelling and collapse of a weakly charged polyelectrolyte gel was developed. In addition to excluded-volume interactions between monomer units of the gel, the theory involves the counterion-monomer unit and counterion-counterion interactions in terms of the virial approximation. The character of interactions between different units in the system varies from repulsion to attraction depending on the type of solvent, counterion, and dielectric permittivity of the solvent. For solvents with a low permittivity, the effect of condensation of counterions resulting in the formation of ion pairs is taken into account.     

Monte Carlo simulations of charged dendrimer-linear polyelectrolyte complexes and explicit counterions

We study complexes composed of one dendrimer of generation G = 4 (G4 dendrimer) with N(t) = 32 charged terminal groups and an oppositely charged linear polyelectrolyte accompanied by neutralizing counterions in an athermal solvent using Monte Carlo simulations based on the bond fluctuation model. In our study both the full Coulomb potential and the excluded volume interactions are taken into account explicitly with the reduced temperature τ and the chain length N(ch) as the main simulation parameters. Our calculations indicate that there exist three temperature ranges that determine the behavior of such complexes. At τ(complex) stable charged dendrimer-linear polyelectrolyte complexes are formed first, which are subsequently accompanied by selective counterion localization within the complex interior at τ(loc) ≤ τ(complex), and counterion condensation as temperature is further decreased below τ(cond) < τ(loc). In particular, we observe that condensation takes place exclusively on the excess charges in the complex and thus no condensation is observed at the compensation point (N(ch) = N(t)), irrespective of τ. For N(ch) ≠ N(t) the complex is overally charged. Furthermore, we discuss the size and structure of the dendrimer and the linear polyelectrolyte within the complex, as well as spatial distributions of monomers and counterions. Conformations of the chain in the bound state are analysed in terms of loops, trains, and tails.     

New experiments for the quantification of counterion condensation

The condensation of counterions is an important aspect of charged macromolecules. Therefore an experimental characterization of the condensation of counterions is desirable. In this contribution two experimental techniques for the characterization of counterion condensation are introduced and compared: Anomalous Small Angle X-Ray Scattering (ASAXS) is able to probe the spatial distribution of counterions and electrophoresis nuclear magnetic resonance (NMR) measures counterion condensation via the effective charge obtained from the dynamic behaviour of molecules and complexes in an electric field.     

Langevin dynamics of semiflexible polyelectrolytes: rod-toroid-globule-coil structures and counterion distribution

We have investigated the nature of counterion condensation on uniformly charged semiflexible polyelectrolyte chains and the concomitant configurations by monitoring the role of chain stiffness, chain length, counterion valency, and the strength of electrostatic interaction. The counterion condensation is seen to follow the adsorption process and the effective polymer charge increases with chain stiffness. Size and shape, as calculated through the radius of gyration, effective persistence length, and hydrodynamic radius, are studied. Stable coil-like, globular, folded-chain, toroidal, and rodlike configurations are possible at suitable combinations of values of chain stiffness, chain length, electrostatic interaction strength, and the valency of counterion. For high strengths of electrostatic interactions, sufficiently stiff polyelectrolytes form toroids in the presence of multivalent counterions, whereas flexible polyelectrolytes form disordered globules. The kinetic features of the nucleation and growth of toroids are monitored. Several metastable structures are found to frustrate the formation of toroids. The generic pathway involves the nucleation of one primary loop somewhere along the chain contour, followed by a growth process where the rest of the chain is folded continuously on top of the primary loop. The dependence of the average radii of toroids on the chain length is found to be roughly linear, in disagreement with existing scaling arguments.     

In-situ investigation of the adsorption of globular model proteins on stimuli-responsive binary polyelectrolyte brushes

Binary brushes constituted from two incompatible polymers can be used in the form of ultrathin polymeric layers as a versatile tool for surface engineering to tune physicochemical surface characteristics such as wettability, surface charge, chemical composition, and morphology and furthermore to create responsive surface properties. Mixed brushes of oppositely charged weak polyelectrolytes represent a special case of responding surfaces that are sensitive to changes in the pH value of the aqueous environment and therefore represent interesting tools for biosurface engineering. The polyelectrolyte brushes used for this study were composed of two oppositely charged polyelelctrolytes poly(2-vinylpyridine) (P2VP) and poly(acrylic acid) (PAA). The in-situ properties and surface characteristics such as as surface charge, surface tension, and extent of swelling of these brush layers are functions of the pH value of the surrounding aqueous solution. To test the behavior of the mixed polylelctrolyte brushes in contact with biosystems, protein adsorption experiments with globular model proteins were performed at different pH values and salt concentrations (confinement of counterions) of the buffer solutions. The influence of the pH value, buffer salt concentration, and isoelectric points (IEP) of the brush and protein on the adsorbed amount and the interfacial tension during protein adsorption as well as the protein adsorption mechanism postulated in reference to recently developed theories of protein adsorption on polyelectrolyte brushes is discussed. In the salted regime, protein adsorption was found to be similar to the often-described adsorption at hydrophobic surfaces. However, in the osmotic regime the balance of electrostatic repulsion and a strong entropic driving force, "counterion release", was found to be the main influence on protein adsorption.     

Theory of competitive counterion adsorption on flexible polyelectrolytes: divalent salts

The counterion distribution around an isolated flexible polyelectrolyte in the presence of a divalent salt is evaluated using the adsorption model [M. Muthukumar, J. Chem. Phys. 120, 9343 (2004)] that considers the Bjerrum length, salt concentration, and local dielectric heterogeneity as physical variables in the system. Self-consistent calculations of effective charge and size of the polymer show that divalent counterions replace condensed monovalent counterions in competitive adsorption. The theory further predicts that at modest physical conditions for a flexible polyelectrolytes such as sodium polystyrene sulfonate in aqueous solutions polymer charge is compensated and reversed with increasing divalent salt. Consequently, the polyelectrolyte shrinks and reswells. Lower temperatures and higher degrees of dielectric heterogeneity between chain backbone and solvent enhance condensation of all species of ions. Complete diagrams of states for the effective charge calculated as functions of the Coulomb strength and salt concentration suggest that (a) overcharging requires a minimum Coulomb strength and (b) progressively higher presence of salt recharges the polymer due to either electrostatic screening (for low Coulomb strengths) or coion condensation (for high Coulomb strengths). Consideration of ion-bridging by divalent counterions leads to a first-order collapse of polyelectrolytes in modest presence of divalent salts and at higher Coulomb strengths. The authors' theoretical predictions are in agreement with the generic results from experiments and simulations.     

Evaluation of the adsorption trends of a low molecular-weight polyelectrolyte with a site-binding model

The site-binding model is very useful for describing the adsorption of ions and small ionized molecules. It has been slightly modified to include multi-site adsorption of larger molecules such as oligomers and low molecular weight polyelectrolytes. We describe alterations of the classical model and the results of calculations for adsorption of polyacrylic acid onto titanium dioxide as an example. The triple layer model is used to relate charge densities to interfacial potential profiles. Comparison between adsorption trends and the surface layer composition as a function of pH and ionic strength demonstrates the prominent influence of ions binding in the adsorption process. The site-binding model makes it easy to simulate the ions displacement associated with polyelectrolyte adsorption. Strongly bound electrolyte anions prevent polyacrylic acid from adsorbing, and, in contrast, electrostatic screening due to cation condensation makes it easier. Calculations of the pH change in the solution, due to adsorption, are also made by comparing ionization ratios of both the surface and polymer units in the adsorbed layer and before adsorption. Trends in electrokinetic potentials as a function of the solution's parameters are evaluated assuming the identity of the shearing surface and the inner boundary of the diffuse layer. All data compare well with experimental values. The very good agreement betwen the experiment and model calculations supports the fact that (small) polyelectrolyte molecules adsorb essentially flat on the surface.     

Simulation Study of the Conformational Properties of Diblock Polyelectrolytes in Salt Solutions

Coarse-grained molecular dynamics simulations are performed to understand the behavior of diblock polyelectrolytes in solutions of divalent salt by studying the conformations of chains over a wide range of salt concentrations. The polymer molecules are modeled as bead spring chains with different charged fractions and the counterions and salt ions are incorporated explicitly. Upon addition of a divalent salt, the salt cations replace the monovalent counterions, and the condensation of divalent salt cations onto the polyelectrolyte increases, and the chains favor to collapse. The condensation of ions changes with the salt concentration and depends on the charged fraction. Also, the degree of collapse at a given salt concentration changes with the increasing valency of the counterion due to the bridging effect. As a quantitative measure of the distribution of counterions around the polyelectrolyte chain, we study the radial distribution function between monomers on different polyelectrolytes and the counterions inside the counterion worm surrounding a polymer chain at different concentrations of the divalent salt. Our simulation results show a strong dependence of salt concentration on the conformational properties of diblock copolymers and indicate that it can tune the self-assembly behaviors of such charged polyelectrolyte block copolymers.