Molecular dynamics simulations of end‐grafted polymers with charged side chains

We report molecular dynamics simulations on bottle‐brush polyelectrolytes end‐grafted to a planar surface. For each bottle‐brush polyelectrolyte, flexible charged side chains are anchored to one neutral main chain. The effects of the counterion valence and the grafting density on the density profiles and the structural characteristics of the brush were studied in this work. It is found that the electrostatic repulsion between charged monomers in the side chains leads an extended conformation of the brush in a solution containing monovalent counterions, while strong electrostatic binding of multivalent counterions to the side chains has a significant contribution to the collapse of the brush. For the trivalent case, the distribution of end monomers in the main chains becomes broader upon decreasing the grafting density, as compared with the monovalent case. However, the position of the distribution for the monovalent case is relatively insensitive to the change of the grafting density. Additionally, with increased counterion valence, enhanced electrostatic correlation between counterions and charged side chains also weakens the diffusive ability of counterions. © 2011 Wiley Periodicals, Inc. J Polym Sci Part B: Polym Phys, 2011     

Conformational transitions of flexible hydrophobic polyelectrolytes in solutions of monovalent and multivalent salts and their mixtures

Conformations of cationic polyelectrolytes (PEs), a weak poly(2-vinylpyridine) (P2VP) and a strong poly(N-methyl-2-vinylpyridinium iodide) (qP2VP), adsorbed on mica from saline solutions in the presence of counterions of different valences are studied using in situ atomic force microscopy (AFM). Quantitative characteristics of chain conformations are analyzed using AFM images of the adsorbed molecules. The results of the statistical analysis of the chain contour reveal collapse of the PE coils when ionic strength is in a range from tens to hundreds of millimoles per kilogram and re-expansion of the coils with a further increase of ionic strength up to a region of the saturated saline solutions. The competition between monovalent and multivalent counterions simultaneously present in solutions strongly affects conformations of PE chains even at a very small fraction of multivalent counterions. Shrinkage of PE coils is steeper for multivalent counterions than for monovalent counterions. However, the re-expansion is only incremental in the presence of multivalent counterions. Extended adsorbed coils at low salt concentrations and at very high concentrations of monovalent salt exhibit conformation corresponding to a 2D coil with 0.95 fraction of bound segments (segments in "trains") in the regime of diluted surface concentration of the PE. Shrunken coils in the intermediate range of ionic strength resemble 3D-globules with 0.8 fraction of trains. The incrementally re-expanded PE coils at a high ionic strength remain unchanged at higher multivalent salt concentrations up to the solubility limit of the salt. The formation of a strong PE complex with multivalent counterions at high ionic strength is not well understood yet. A speculative explanation of the observed experimental result is based on possible stabilization of the complex due to hydrophobic interactions of the backbone.     

Anionic polyelectrolyte adsorption on mica mediated by multivalent cations: a solution to DNA imaging by atomic force microscopy under high ionic strengths

Adsorption of DNA molecules on mica, a highly negatively charged surface, mediated by divalent or trivalent cations is considered. By analyzing atomic force microscope (AFM) images of DNA molecules adsorbed on mica, phase diagrams of DNA molecules interacting with a mica surface are established in terms of concentrations of monovalent salt (NaCl) and divalent (MgCl2) or multivalent (spermidine, cobalt hexamine) salts. These diagrams show two transitions between nonadsorption and adsorption. The first one arises when the concentration of multivalent counterions is larger than a limit value, which is not sensitive to the monovalent salt concentration. The second transition is due to the binding competition between monovalent and multivalent counterions. In addition, we develop a model of polyelectrolyte adsorption on like-charged surfaces with multivalent counterions. This model shows that the correlations of the multivalent counterions at the interface between DNA and mica play a critical role. Furthermore, it appears that DNA adsorption takes place when the energy gain in counterion correlations overcomes an energy barrier. This barrier is induced by the entropy loss in confining DNA in a thin adsorbed layer, the entropy loss in the interpenetration of the clouds of mica and DNA counterions, and the electrostatic repulsion between DNA and mica. The analysis of the experimental results provides an estimation of this energy barrier. We then discuss some important issues, including DNA adsorption under physiological conditions.     

Chain stiffness, salt valency, and concentration influences on titration curves of polyelectrolytes: Monte Carlo simulations

Monte Carlo simulations have been used to study two different models of a weak linear polyelectrolyte surrounded by explicit counterions and salt particles: (i) a rigid rod and (ii) a flexible chain. We focused on the influence of the pH, chain stiffness, salt concentration, and valency on the polyelectrolyte titration process and conformational properties. It is shown that chain acid-base properties and conformational properties are strongly modified when multivalent salt concentration variation ranges below the charge equivalence. Increasing chain stiffness allows to minimize intramolecular electrostatic monomer interactions hence improving the deprotonation process. The presence of di and trivalent salt cations clearly promotes the chain degree of ionization but has only a limited effect at very low salt concentration ranges. Moreover, folded structures of fully charged chains are only observed when multivalent salt at a concentration equal or above charge equivalence is considered. Long-range electrostatic potential is found to influence the distribution of charges along and around the polyelectrolyte backbones hence resulting in a higher degree of ionization and a lower attraction of counterions and salt particles at the chain extremities.     

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.     

Surface Patterns of a Tetrahedral Polyelectrolyte Brush Induced by Grafting Density and Charge Fraction

A tetrahedral polyelectrolyte brush in the presence of trivalent counterions is researched under the condition of good solution by mea ns of molecular dynamics simulati ons.Grafting density and charge fraction are varied to gen erate a series of surface patter ns.Lateral microphase separation happens and various interesting pinned patches appear at appropriate charge fraction and grafting density.Through a careful analysis on the brush thickness,the pair correlation functions,the distributions of net charge,and the four states of trivalent counterions in the brush,we find that the ordered surface patterns and special properties are induced by the pure electrostatic correlation effect of trivalent ions even in the good solvent.Furthermore,the dependences of electrostatic correlation on the charge fraction of tethered chains are evaluated for fixed grafting den sity.Also,our results can serve as a guide for precise control over the stimuli-responsive materials rational and self-assembly of nanoparticles.     

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.     

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.     

Conformation of poly(styrene sulfonate) layers physisorbed from high salt solution studied by force measurements on two different length scales

The conformation of poly(styrene sulfonate) (PSS) layers physisorbed from 1 M NaCl is determined by force measurements and imaging on two length scales. With colloidal probe technique steric forces as predicted for neutral grafted brushes are observed. On decrease and increase of the NaCl concentration, the grafting density remains constant, yet the brush thickness swells and shrinks reversibly with the salt concentration with an exponent of -0.3. At low salt conditions, the brush length amounts to 30% of the contour length, a behavior known for polyelectrolyte brushes and attributed to the entropy of the counterions trapped in the brush. Between a PSS layer and a pure colloidal silica sphere, the same steric forces are observed, and additionally at large separations (beyond the range of the steric repulsion) an electrostatic force is found. A negatively charged AFM tip penetrates the brush--a repulsive electrostatic force between the tip and surface is found, and single chains can be imaged. Thus, with the nanometer-sized AFM tip, the flatly adsorbed fraction of the PSS chains is seen, whereas the micrometer-sized colloidal probe interacts with the fraction of the chains penetrating into solution.     

Overcharging of nanoparticles in electrolyte solutions

Monte Carlo simulations are performed to investigate the effects of salt concentration, valence and size of small ions, surface charge density, and Bjerrum length on the overcharging of isolated spherical nanoparticles within the framework of a primitive model. It is found that charge inversion is most probable in solutions containing multivalent counterions at high salt concentrations. The maximum strength of overcharging occurs near the nanoparticle surface where counterions and coions have identical local concentrations. The simulation results also suggest that both counterion size and electrostatic correlations play major roles for the occurrence of overcharging.     

Synthesis and swelling behavior of pH-responsive polybase brushes

We synthesize polybase brushes and investigate their swelling behavior. Poly(2-(dimethylamino)ethyl methacrylate)) (PDMAEMA) brushes are prepared by the "grafting from" method using surface-initiated Atom Transfer Radical Polymerization to obtain dense brushes with relatively monodisperse chains (PDI = 1.35). In situ quaternization reaction can be performed to obtain poly(2-(trimethylamino)ethyl methacrylate)) (PTMAEMA) brushes. We determine the swollen thickness of the brushes using ellipsometry and neutron reflectivity techniques. Brushes are submitted to different solvent conditions to be investigated as neutral brushes and weak and strong polyelectrolyte brushes. The swelling of the brushes is systematically compared to scaling models. It should be pointed out that the scaling analysis of different types of brushes (neutral polymer and weak and strong polyelectrolyte brushes) is performed with identical samples. The scaling behavior of the PDMAEMA brush in methanol and the PTMAEMA brush in water is in good agreement with the predicted scaling laws for a neutral polymer brush in a good solvent and a polyelectrolyte brush in the osmotic regime. The salt-induced contraction of the quaternized brush is observed for high salt concentration, in agreement with the predicted transition between the regimes of the osmotic brush and the salted brush. From the crossover concentration, we calculate the effective charge ratio of the brush following the Manning counterion condensation. We also use PDMAEMA brushes as pH-responsive polybase brushes. The swelling behavior of the polybase brush is intermediate with respect to the behavior of the neutral polymer brush in a good solvent and the behavior of the quenched polyelectrolyte brush, as expected. The effective charge ratio of the PDMAEMA brush is determined as a function of pH using the scaling law of the polyelectrolyte brush in the osmotic regime.     

Effects of counter ions of clay platelets on the swelling behavior of nanocomposite gels

The effects of replacing the native Na(+) counter ions associated with the clay platelets by various other cations on the swelling behavior of nanocomposite (NC) gels consisting of an organic (polymer)/inorganic (clay) network were investigated. The negative surface charge of the clay platelet conferred an ionic nature on the NC gels making them a type of polyelectrolyte gel; consequently, the swelling behavior of the NC gels was strongly influenced by the valence of the co-existing counter ions. NC gels containing monovalent cations such as Na(+), K(+) and Li(+) exhibited large swellings and subsequent deswelling in water after attaining maximum degrees of swelling. In contrast, introduction of multivalent cations such as Ca(2+), Mg(2+), and Al(3+) into NC gels depressed markedly both the swelling and subsequent deswelling. The decreased swelling and suppressed deswelling with multivalent ions were strongly influenced by the initial gel state and result from the formation of additional cross-links through ionic interactions between the clay platelets and the multivalent cations. Also, the similar swelling behaviors were observed for all NC gels with different clay concentration. Further, reversible absorption/desorption and selective absorption of multivalent cations were observed for the NC gels examined.     

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.     

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.     

High elongation of polyelectrolyte chains in the osmotic limit of spherical polyelectrolyte brushes: a study by cryogenic transmission electron microscopy

We investigate the conformation of long polyelectrolyte chains attached to colloidal latex particles by cryogenic transmission electron microscopy (cryo-TEM). The dense grafting of the polyelectrolyte chains ("polyelectrolyte brush") leads to a confinement of the counterions and a concomitantly high osmotic pressure within the polyelectrolyte layer attached to the core particles. Cryo-TEM has provided first model-independent direct proof for the strong stretching of the polyelectrolyte chains by direct visualization. If salt is added, cryo-TEM clearly shows how chains collapse because of the strong screening of the electrostatic interaction. Moreover, the analysis of interacting particles by cryo-TEM shows that the polyelectrolyte chains retract at close contact. Hence, we demonstrate how cryo-TEM can be used to analyze directly the spatial structure of polyelectrolyte brushes in situ.     

On the monomer density of grafted polyelectrolyte brushes and their interactions

Most of the modern theories of grafted polyelectrolyte brushes are valid only for moderate stretching of the polyelectrolyte. However, particularly at low ionic strength and high grafting densities, even a moderate charge of the polyelectrolyte can generate a strong stretching. A simple mean field model for strongly stretched grafted polyelectrolyte brushes is suggested, based on an approximate calculation of the partition function of a polyelectrolyte chain. It is shown that the average Boltzmann factor of a possible chain configuration can be approximated by the Boltzmann factor of a configuration with a constant monomer distribution, for which the free energy can be readily obtained. The monomer density in the brush and the interaction between two surfaces with grafted polyelectrolyte brushes could be calculated as a statistical average over all possible configurations. Some simple analytical results are derived, and their accuracy is examined. The dependence of the brush thickness on the electrolyte concentration is investigated, and it is shown that the trapping of a fraction of counterions in the brush influences strongly the thickness of the brush. When two surfaces with grafted polyelectrolyte brushes approach each other more rapidly than the ion diffusion parallel to the surface, the trapping of the counterions between the brushes can affect the interactions by orders of magnitude.     

Salt effect on the nanostructure of strong polyelectrolyte brushes in amphiphilic diblock copolymer monolayers on the water surface

The nanostructure of a spread monolayer of diblock copolymers of poly(hydrogenated isoprene) and poly(styrenesulfonate) at the air/water interface were studied by in situ X-ray reflectivity as a function of the brush density and salt concentration. When the monolayer was compressed beyond the "critical brush density", its nanostructure changed from a flat, adsorbed "carpet" layer to a "carpet + brush" structure. The critical brush density was found to be about 0.12 nm-2, independent the brush length and salt concentration under a low-salt condition. The brush formation behavior was considered to be controlled by an electrostatic interaction between polyelectrolyte chains rather than by a steric hindrance. This might be because the distance between the chains at the critical point is rather long and also because of the effect of the salt on the critical brush density. The critical brush density increased at higher added salt concentration beyond 1 M. As a result, we found a new structure transition behavior of the polymer brushes between carpet-only and carpet + brush structures, which was induced by salt addition. Finally, we succeeded in the controlled release of salt ions from the salted brush layer by changing the brush density by compression of the monolayer.     

Lattice Monte Carlo simulations of a charged polymer chain: effect of valence and concentration of the added salt

The configurational properties of a single polyelectrolyte chain accompanied by counterions and added salt are simulated using the cooperative motion algorithm on the face-centered cubic lattice. In particular, a greater emphasis is put on the effect of valence z(s) and concentration of the added positive (negative) salt ions n(s) on the polymer behavior. This is achieved by inspecting two families of systems with widely varying numbers n(s) of monovalent (z(s)=1) or multivalent (z(s)=4) salt ions at two fixed reduced temperatures T*=0.5, 1. The calculations indicate that especially at the lower temperature the addition of some amount of multivalent salt has a tremendous impact on chain conformations compared to the situation with monovalent salt. Even for relatively low concentrations of the former, the mean radius of gyration (1/2) and the mean end-to-end distance (1/2) decrease sharply, i.e., the polymer exists in strongly collapsed forms. This reduction of polymer size is also accompanied by a drop in the system inner energy e* and the effective mean charge per monomer q*. The analysis of various pair-correlation functions g(ab)(r) indicates that the latter effect-caused by condensation of ions onto the chain-is dominated by the multivalent ones. Furthermore, it is found that for z(s)=4, the uncondensed salt ions tend to group themselves into small clusters.     

Packaging of a flexible polyelectrolyte inside a viral capsid: effect of salt concentration and salt valence

The effect of salt on the location and structure of a flexible polyelectrolyte confined inside a viral capsid and the Donnan equilibrium of the salt across the capsid have been examined using a coarse-grained model solved by Monte Carlo simulations. The polyelectrolyte was represented by a linear jointed chain of charged beads, and the capsid was represented by a spherical shell with embedded charges. At low salt concentration, the polyelectrolyte was strongly adsorbed onto the inner capsid surface, whereas at high salt concentration it was located preferentially in the central part of the capsid. Under the condition of equal Debye screening length, the electrostatic screening increased as the valence of the polyelectrolyte counterion was increased. The distribution of the small cations and anions was unequal across the capsid. An excess of polyelectrolyte counterions occurred inside the capsid, and the excess increased with the salt concentration. A simplified representation of the small ions through the use of the screened Coulomb potential provided only a qualitatively correct picture; the electrostatic screening originating from the small ions was exaggerated.     

Poisson-Boltzmann theory of pH-sensitive (annealing) polyelectrolyte brush

We present a self-consistent field analytical theory of a polymer brush formed by weakly charged pH-sensitive (annealing) polyelectrolytes tethered to a solid-liquid interface and immersed in buffer solution of low molecular weight salt. We use the Poisson-Boltzmann framework, applied by us previously to polyelectrolyte (PE) brushes with quenched charge (Zhulina, E. B.; Borisov, O. V. J. Chem. Phys. 1997, 107, 5952). This approach allows for detailed analysis of the internal structure of annealing PE brush in terms of polymer density distribution, profiles of electrostatic potential and of local degree of chain ionization as a function of buffer ionic strength and pH without any assumptions on mobile ion distribution imposed in earlier scaling-type models. The presented analytical theory recovers all major asymptotic dependences for average brush properties predicted earlier. In particular, a nonmonotonic dependence of brush thickness on ionic strength and grafting density is confirmed and specified with accuracy of numerical coefficients including crossover regions. Moreover, the theory predicts qualitatively new effects, such as, e.g., disproportionation of tethered polyions into weakly charged concentrated proximal and strongly charged sparse distal brush domains at low salt and moderate grating densities. The presented results allow us to quantify responsive features of annealing PE brushes whose large-scale and local conformational properties can be manipulated by external stimuli.