Nanoparticle adsorption on a weak polyelectrolyte. Stiffness, pH, charge mobility, and ionic concentration effects investigated by Monte Carlo simulations

Monte Carlo simulations have been used to study two different models for a weak linear polyelectrolyte in the presence of nanoparticles: (i) a rodlike and (ii) a flexible polyelectrolytes. The use of simulated annealing has made it possible to simulate a polyelectrolyte chain in the presence of several nanoparticles by improving conformation sampling and avoiding multiple minima problems when dense conformations are produced. Nanoparticle distributions along the polymer backbone were analyzed versus the ionic concentration, polyelectrolyte stiffness, and nanoparticle surface charge. Titration curves were calculated and the influences of the ionic concentration, solution pH, and number of adsorbed nanoparticles on the acid/base polyelectrolyte properties have been systematically investigated. The subtle balance of attractive and repulsive interactions has been discussed, and some characteristic conformations are presented. The comparison of the two limit models provides a good representation of the stiffness influence on the complex formation. In some conditions, overcharging was obtained and presented with respect to both the polyelectrolyte and nanoparticle as the central element. Finally, the charge mobility influence along the polyelectrolyte backbone was investigated by considering annealed and quenched polyelectrolyte chains.     

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.     

Molecular dynamics simulations of dendritic polyelectrolytes with flexible spacers in salt free solution

We present the results of molecular dynamics simulations of dendritic polyelectrolytes in dilute salt-free solutions. The dendritic polyelectrolytes are modeled as an ensemble of regular-branched bead-spring chains of neutral and charged Lennard-Jones particles with explicit counterions. A wide range of molecular variables of the dendritic polyelectrolytes such as generation number, spacer length, and charge density were considered in the simulations. The effect of dendrimer size on relaxation time, the conformation of spacers, and the size dependence of the dendrimer on molecular variables are discussed and compared with a Flory type theory. The osmotic coefficients of the dilute dendritic polyelectrolyte solutions, as well as the profiles of monomers and counterions, are calculated directly from the simulations. Our simulation results show that the inner spacers of the dendrimers are extensively stretched, and the size dependence on the molecular weight deviates from the scaling prediction that assumes a Gaussian elasticity of the spacer.     

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.     

Electrostatic-interaction-induced phase separation in solutions of flexible-chain polyelectrolytes

A model of a polyelectrolyte solution has been formulated on the basis of the formalism of the thermodynamic perturbation theory. Macromolecules have been described in terms of the model of a flexible chain with an excluded volume and a variable electrical charge. During construction of the thermodynamic perturbation theory, a set of three independent subsystems—polyelectrolyte macromolecules placed in a structureless charge background of counterions, counterions placed in a structureless charge background of macromolecules, and Coulomb gas ions of a low-molecular-mass salt—has been taken as the reference system. In the framework of this model, liquid-liquid phase separation due to strong correlation-induced attraction has been predicted. The behavior of the degree of ionization over a wide monomer concentration range, including the region of phase separation either in a salt-free solution or in the presence of univalent ions of a low-molecular-mass salt in the solution, has been studied. It has been shown that macromolecules in the coexisting phases should have different degrees of ionization. The occurrence of phase separation under normal conditions in the case when dimethylformamide is taken as a solvent and the nonoccurrence of this phase separation in the case of aqueous solutions of flexible-chain polyelectrolytes are predicted.     

Influence of weak groups on polyelectrolyte mobilities

The ionization of dissociable groups in weak polyelectrolytes does not occur in a homogenous fashion. Monomer connectivity imposes constraints on the localization of the dissociated (charged) monomers that affect the local electric potential. As a result, the mean bare charge along a weak polyelectrolyte can vary depending on the proximity to topological features (e.g. presence of crosslinks or dangling ends). Using reaction‐ensemble Monte‐Carlo simulations we calculate the dissociation inhomogeneities for a few selected PE configurations, linear, rod‐like, flexible four‐arm star, and a star with stiff arms. An ensemble preaverage is used to obtain the annealed bare charge profile for these different polymer configurations. Using molecular dynamics simulations within a Lattice‐Boltzman fluid, we investigate how the electrophoretic mobility is affected by the bare charge inhomogeneities arising from the annealed weak polyelectrolytes. Surprisingly, the mobility obtained for the situations corresponding to the predicted charge profile for annealed weak polyelectrolytes are not significantly different than the mobility obtained when all the monomers have an identical charge (under the constraint that the total polyelectrolyte bare charge is the same). This is also true for the stiff rod‐like variants where conformational changes induced from the localization of the monomer charges are negligible. In salty solutions, we find that counterions are affected by the electric potential modulations induced by the topological features. Since the counterions crowd in regions where the electric potential caused by the dissociated monomers is highest, they wash‐out the bare charge inhomogeneities and contribute to a more uniform effective backbone charge.     

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.     

Electrophoresis of a spherical dispersion of polyelectrolytes in a salt-free solution

The electrophoretic behavior of a spherical dispersion of polyelectrolytes of arbitrary concentration is analyzed theoretically under a salt-free condition, that is, the liquid phase contains only counterions which come from the dissociation of the functional groups of polyelectrolytes. We show that, in general, the surface potential of a polyelectrolyte increases nonlinearly with its surface charge. A linear relation exists between them, however, when the latter is sufficiently small; and the more dilute the concentration of polyelectrolytes, the broader the range in which they are linearly correlated. If the amount of surface charge is sufficiently large, counterion condensation occurs, and the rate of increase of surface potential as the amount of surface charge increases declined. Also, it leads to an inverse in the perturbed potential near the surface of a polyelectrolyte, and its mobility decreases accordingly. For a fixed amount of surface charge, the lower the concentration of polyelectrolytes and/or the lower the valence of counterions, the higher the surface potential. The qualitative behavior of the mobility of a polyelectrolyte as the amount of its surface charge varies is similar to that of its surface charge.     

Theory of counter-ion condensation on flexible polyelectrolytes: adsorption mechanism

A new model is presented for counterion distribution around flexible polyelectrolytes by considering (i) free energy of the polyelectrolyte chain, (ii) translational entropy of adsorbed counterions, (iii) adsorption energy, (iv) translational entropy of unadsorbed counterions, (v) fluctuations of dissociated ions, and (vi) correlation among ion-pairs formed by adsorbed counterions on the polymer. The effective charge and size of the polymer are calculated self-consistently. The degree of ionization f of the polymer decreases continuously with 1/epsilonT (epsilon and T are the dielectric constant of the solvent and temperature, respectively), depending sensitively on local dielectric heterogeneity. Further, f decreases with an increase in salt concentration, monomer concentration, or chain flexibility. The polymer size, accompanying the changes in f, depends nonmonotonically on 1/epsilonT. The predictions of the model are consistent with all trends observed previously in simulations and are distinctly different from the Manning argument for rodlike chains.     

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.     

Supramolecular Assembly of Poly(propyleneimine) Dendrimers Driven By Simple Monovalent Counterions

The self‐assembly of semiglobular, positively charged poly(propyleneimine) (PPI) dendrimers with small monovalent counterions (e.g., Cl^?) in water/acetone mixtures was investigated. We showed that PPI dendrimers can assemble into hollow, spherical, single‐layered blackberry‐type structures mediated by the presence of monovalent counterions. The effects on the assembly of changing the solvent polarity and adjusting the pH were further investigated to confirm the presence of electrostatic interactions and hydrogen bonding as the driving forces. Results showed that PPI dendrimers form stable, hollow spheres in 5–20 % v/v acetone/water and that the size of the spheres decreases monotonically as the solvent polarity and/or the charge on the dendrimers (i.e., lower solution pH) increases. This is the first example to show that small monovalent counterions can trigger attraction among PPI dendrimers (or broadly defined polyelectrolytes) that is strong enough to bring them together to form large, stable supramolecular assemblies, which indicates that these organic macroions have similar solution behavior to more‐well‐defined inorganic molecular macroions.     

Conformational effect on small angle neutron scattering behavior of interacting polyelectrolyte solutions: a perspective of integral equation theory

We present small angle neutron scattering (SANS) measurements of deuterium oxide (D(2)O) solutions of linear and star sodium poly(styrene sulfonate) (NaPSS) as a function of polyelectrolyte concentration. Emphasis is on understanding the dependence of their SANS coherent scattering cross section I(Q) on the molecular architecture of single polyelectrolyte. The key finding is that for a given concentration, star polyelectrolytes exhibit more pronounced characteristic peaks in I(Q), and the position of the first peak occurs at a smaller Q compared to their linear counterparts. Based on a model of integral equation theory, we first compare the SANS experimental I(Q) of salt-free polyelectrolyte solutions with that predicted theoretically. Having seen their satisfactory qualitative agreement, the dependence of counterion association behavior on polyelectrolyte geometry and concentration is further explored. Our predictions reveal that the ionic environment of polyelectrolyte exhibits a strong dependence on polyelectrolyte geometry at lower polyelectrolyte concentration. However, when both linear and star polyelectrolytes exceed their overlap concentrations, the spatial distribution of counterion is found to be essentially insensitive to polyelectrolyte geometry due to the steric effect.     

Proteins and polyelectrolytes: A charged relationship

We review the interaction of charged polymeric systems with proteins. In solutions of low ionic strength there are many examples of proteins attracted to polyelectrolytes even if both systems carry the same overall charge. This attractive interaction is widespread, having been observed for single polyelectrolyte chains as well as for polyelectrolytes grafted to surfaces (polyelectrolyte brushes) and charged polymeric networks. In all cases, adding salt weakens the interaction considerably. We discuss the suggestion that the attractive force at low salinity originates from the asymmetry of interaction between charged polymer segments and charged patches on the surface of the protein globule. This can be explained if the attractive force is mainly due to a counterion release force, i.e., the polyelectrolyte chains become the multivalent counterions for the patches of opposite charge localized on the surface of the proteins. We review a selection of simple models that lead to semi-quantitative estimates of this force as the function of salt concentration.     

Linear polyelectrolytes in tetravalent salt solutions

The effect of adding tetravalent salt of different sizes to a solution of linear and flexible polyelectrolytes is investigated by molecular dynamics simulations. Upon the addition of salt, a chain reexpansion takes place, following a well-known collapsed conformation. The degrees of collapse and reexpansion increase with ion size. In the solution, tetravalent counterions replace monovalent ones and condense onto the chains. The condensation for small ions displays a profile different from that for large ones. In a high-salt region, ions can form layering orders around a polyelectrolyte and locally overcompensate the charge inside. Consequently, the integrated charge distribution reveals an oscillatory behavior away from a chain. By studying the radial distribution function between monomers on different polyelectrolytes, like-charge attraction between chains is demonstrated. This attraction is a prerequisite to chain aggregation or precipitation. The results show a strong dependence of salt concentration and ion size on the properties of polyelectrolyte solutions.     

Spectral collocation methods for polymer brushes

We provide an in-depth study of pseudo-spectral numerical methods associated with modeling the self-assembly of molten mixed polymer brushes in the framework of self-consistent field theory (SCFT). SCFT of molten polymer brushes has proved numerically challenging in the past because of sharp features that arise in the self-consistent pressure field at the grafting surface due to the chain end tethering constraint. We show that this pressure anomaly can be reduced by smearing the grafting points over a narrow zone normal to the surface in an incompressible model, and/or by switching to a compressible model for the molten brush. In both cases, we use results obtained from a source (delta function) distribution of grafting points as a reference. At the grafting surface, we consider both Neumann and Dirichlet conditions, where the latter is paired with a masking method to mimic a confining surface. When only the density profiles and relative free energies of two comparison phases are of interest, either source or smeared distributions of grafting points can be used, but a smeared distribution of grafting points exhibits faster convergence with respect to the number of chain contour steps. Absolute free energies converge only within the smeared model. In addition, when a sine basis is used with the masking method and a smeared distribution, fewer iterations are necessary to converge the SCFT fields for the compressible model. The numerical methods described here and investigated in one-dimension will provide an enabling platform for computationally more demanding three-dimensional SCFT studies of a broad range of mixed polymer brush systems.     

Diffusion of cationic polyelectrolytes into cellulosic fibers

The penetration of cationic polyelectrolytes into anionic cellulosic fibers was evaluated with fluorescent imaging techniques in order to clarify the mechanism and time scales for the diffusion process. The bulk charge of the cellulosic fibers indirectly creates a driving force for diffusion into the porous fiber wall, which is entropic in nature due to a release of counterions as the polyelectrolyte adsorbs. The individual bulk charges in the fiber cell wall also interact with the diffusing polyelectrolyte, such that the polyelectrolyte diffuses to the first available charge and consequently adsorbs and remains fixed. Thus, subsequent polyelectrolyte chains must first diffuse through the adsorbed polyelectrolyte layer before adsorbing to the next available bulk charges. This behavior differs from earlier suggested diffusion mechanisms, by which polyelectrolytes were assumed to first adsorb to the outermost surface and then reptate into the pore structure. The time scales for polyelectrolyte diffusion were highly dependent on the flexibility of the chain, which was estimated from calculations of the persistence length. The persistence length ultimately depended on the charge density and electrolyte concentration. The charge density of the polyelectrolyte had a greater influence on the time scales for diffusion. High charge density polyelectrolytes were observed to diffuse on a time scale of months, whereas the diffusion of low charge density polyelectrolytes was measured on the order of hours. An influence of the chain length, that is, steric interactions due the persistence length of the polyelectrolyte and to the tortuosity of the porous structure of the fiber wall, could only be noted for low charge density polyelectrolytes. Increasing the electrolyte concentration increased the chain flexibility by screening the electrostatic contribution to the persistence length, in turn inducing a faster diffusion process. However, a significant change in the diffusion behavior was observed at high electrolyte concentrations, at which the interaction between the polyelectrolyte charges and the fiber charges was almost completely screened.     

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.     

Density functional theory for planar electric double layers: closing the gap between simple and polyelectrolytes

We report a nonlocal density functional theory (NLDFT) for polyelectrolyte solutions within the primitive model; i.e., the solvent is represented by a continuous dielectric medium, and the small ions and polyions by single and tangentially connected charged hard spheres, respectively. The excess Helmholtz energy functional is derived from a modified fundamental measure theory for hard-sphere repulsion, an extended first-order thermodynamic perturbation theory for chain connectivity, and a quadratic functional Taylor expansion for electrostatic correlations. With the direct and cavity correlation functions of the corresponding monomeric systems as inputs, the NLDFT predicts the segment-level microscopic structures and adsorption isotherms of polyelectrolytes at oppositely charged surfaces in good agreement with molecular simulations. In particular, it faithfully reproduces the layering structures of polyions, charge inversion, and overcharging that cannot be captured by alternative methods including the polyelectrolyte Poisson-Boltzmann equation and an earlier version of DFT. The NLDFT has also been used to investigate the influences of the small ion valence, polyion chain length, and size disparity between polyion segments and counterions on the microscopic structure, mean electrostatic potential, and overcharging in planar electric double layers containing polyelectrolytes.     

Ionizable polyelectrolyte brushes: brush height and electrosteric interaction

Semi-analytical scaling theory is used to describe quenched and annealed (weakly charged, ionizable, charge-regulating) polyelectrolyte brushes in electrolyte solutions of arbitrary salt concentration. An Alexander-De Gennes box model with homogeneous distribution of polymer segments and the free ends located at the edge of the brush is assumed, as is local electroneutrality in the brush. For annealed polyelectrolyte and in the low-salt regime, the theory predicts that for sufficiently dense brushes, the salt concentration has a small influence on brush height, while the brush expands with increasing grafting density, in agreement with experiment. Expressions are presented for the interaction free energy of compressed ionizable and quenched polyelectrolyte brushes (proportional to the force between particles or curved surfaces). In all cases, the required prefactors are explicitly stated. The theory is compared directly with published experiments on the influence of salt concentration, pH, and grafting density on the thickness and interaction force of polystyrene sulfonate (quenched) and poly(meth)acrylic acid (annealed) brushes. In general, trends are well reproduced but significant deviations remain.     

Explicit ions condensation around strongly charged polyelectrolytes and spherical macroions: the influence of salt concentration and chain linear charge density. Monte Carlo simulations

The condensation of monovalent counterions and trivalent salt particles around strong rigid and flexible polyelectrolyte chains as well as spherical macroions is investigated by Monte Carlo simulations. The results are compared with the condensation theory proposed by Manning. Considering flexible polyelectrolyte chains, the presence of trivalent salt is found to play an important role by promoting chain collapse. The attraction of counterions and salt particles near the polyelectrolyte chains is found to be strongly dependent on the chain linear charge density with a more important condensation at high values. When trivalent salt is added in a solution containing monovalent salt, the trivalent cations progressively replace the monovalent counterions. Ion condensation around flexible chains is also found to be more efficient compared with rigid rods due to monomer rearrangement around counterions and salt cations. In the case of spherical macroions, it is found that a fraction of their bare charge is neutralized by counterions and salt cations. The decrease of the Debye length, and thus the increase of salt concentration, promotes the attraction of counterions and salt particles at the macroion surface. Excluded volume effects are also found to significantly influence the condensation process, which is found to be more important by decreasing the ion size.