Molecular dynamics simulations of end-grafted centipede-like polymers with stiff charged side chains

We use molecular dynamics simulations to investigate centipede-like polymers with stiff charged side chains, end-grafted to a planar wall. The effect of the grafting density and the Bjerrum length on the conformational behaviour of the brush is examined in detail. In addition, we make a comparison of centipede-like polyelectrolyte (CPE) brushes with neutral centipede-like polymer (NCP) and linear polyelectrolyte (LPE) brushes. At weak electrostatic interaction, the main chains of the CPE chains adopt a strongly stretched conformation, and the monomer density profiles of side chains exhibit a clear oscillatory behaviour. With increasing Bjerrum length, the CPE brush undergoes a collapse transition. Compared to the CPE brushes, the counterion condensation effect is stronger for the LPE brushes, regardless of whether the electrostatic interaction is weak or strong and of whether the grafting density is low or high. Additionally, it is shown that the architecture of the grafted chains makes a weak contribution to the counterion condensation at strong electrostatic interaction. We also find that the electrostatic repulsion between charged side chains can enhance the stiffness of the main chains and thus limit the range of movement of the free-end monomers.     

Oppositely charged polyelectrolytes grafted onto planar surface: Mean-field lattice theory

Equilibrium structures of planar polyelectrolyte brushes formed by grafted chains carrying charges of opposite sign are examined by employing mean-field lattice theory. Two brushes of different architecture are considered: one formed by grafted diblock copolymers with oppositely charged blocks and the other being a mixed brush composed of oppositely charged homopolymers. The systems display nontrivial intrinsically inhomogeneous brush structures originating from the chain connectivity and the electrostatic interaction among the segments. In addition, a coexistence of stretched and coiled chains inside the brush is observed. The influence of the charges of the blocks, the relative length of the oppositely charged blocks, and the ionic strength of the solution on the brush inhomogeneity and structural differences between the two types of brushes are discussed. Received 14 March 2001 and Received in final form 18 June 2001     

Effects of counterion fluctuations in a polyelectrolyte brush

We investigate the effect of counterion fluctuations in a single polyelectrolyte brush in the absence of added salt by systematically expanding the counterion free energy about Poisson-Boltzmann mean-field theory. We find that for strongly charged brushes, there is a collapse regime in which the brush height decreases with increasing charge on the polyelectrolyte chains. The transition to this collapsed regime is similar to the liquid-gas transition, which has a first-order line terminating at a critical point. We find that, for monovalent counterions, the transition is discontinuous in theta solvent, while for multivalent counterions, the transition is generally continuous. For collapsed brushes, the brush height is not independent of grafting density as it is for osmotic brushes, but scales linear with it.Received: 26 November 2003, Published online: 11 May 2004PACS: 61.25.Hq Macromolecular and polymer solutions; polymer melts; swelling - 61.20.Qg Structure of associated liquids: electrolytes, molten salts, etc.     

水合作用与pH调控聚电解质刷的构象转变

我们基于Flory-Huggins理论,建立理论模型研究水合作用与pH调控聚电解质刷的构象转变.理论模型考虑聚电解质链与水分子间的作用(聚电解质链的水合作用)、体系中的静电作用.研究发现,随着水合作用的改变,聚电解质刷出现由溶胀到塌缩的构象转变.由此表明了水合作用可在很大程度调节聚电解质刷的相变.通过分析pH的调控效应我们还发现,在碱性环境中(pH=8),聚电解质链单体的解离度增大,静电排斥会使得聚电解质刷溶胀.由此表明,聚电解质刷内水合作用与静电效应的耦合,将会共同决定聚电解质刷的构象转变特性.理论结果深刻揭示了水合作用的改变,会使得聚电解质刷体系发生相变,pH可在很大程度上改变其相变特性.     

Counterion condensation theory of attraction between like charges in the absence of multivalent counterions

There is abundant experimental evidence suggesting the existence of attractive interactions among identically charged polyelectrolytes in ordinary salt solutions. The presence of multivalent counterions is not required. We review the relevant literature in detail and conclude that it merits more attention than it has received. We discuss also some recent observations of a low ionic strength attraction of negatively charged DNA to the region of a negatively charged glass nanoslit where the floor of the nanoslit meets the walls, again in the absence of multivalent ions. On the theoretical side, it has become clear that purely electrostatic interactions require the presence of multivalent counterions if they are to generate like-charge attraction. Any theory of like-charge attraction in the absence of multivalent counterions must therefore contain a non-electrostatic component. We point out that counterion condensation theory, which has predicted like-charge polyelectrolyte attraction in an intermediate range of distances in ordinary 1:1 salt conditions, contains both electrostatic and non-electrostatic elements. The non-electrostatic component of the theory is the modeling constraint that the counterions fall into two explicit populations, condensed and uncondensed. As reviewed in the paper, this physically motivated constraint is supported by strong experimental evidence. We proceed to offer an explanation of the nanoslit observations by showing in an idealized model that the line of intersection of two intersecting planes is a virtual polyelectrolyte. Since we have previously developed a counterion condensation theory of attraction of two like-charged polyelectrolytes, our suggestion is that the DNA is attracted to the virtual polyelectrolytes that may be located in the nanoslit where floor meets walls. We present the detailed calculations needed to document this suggestion: an extension of previous theory to the case of polyelectrolytes with like but not identical charges; the demonstration of counterion condensation on a plane with bare charge density greater than an explicitly exhibited critical value; a calculation of the free energy of the plane; a calculation of the interaction of a line charge polyelectrolyte with a like-charged plane; and the detailed demonstration that the line of intersection of two planes is a virtual polyelectrolyte.     

中性/聚电解质高分子混合刷对蛋白质的吸附/解吸附

本文应用分子理论,研究中性(A)/聚电解质(B)高分子混合刷对蛋白质的吸附/解吸附特性.理论模型考虑蛋白质与中性高分子A的排斥、以及与聚电解质高分子B的静电吸引.研究发现,在pH=4～6、中性高分子A处于弱水合状态时,混合刷中A高分子链塌缩,B聚电解质链溶胀.由于蛋白质和B聚电解质链间的静电吸引,导致高分子混合刷对蛋白质的吸附.当A高分子水合性增强时,A高分子链溶胀,B聚电解质链塌缩.由于蛋白质与A高分子链间的排斥作用增强,与B聚电解质链间的静电吸引减弱,混合刷对蛋白质解吸附.     

Polyelectrolytes in the presence of multivalent ions: gelation versus segregation

We analyze solutions of strongly charged chains bridged by linkers such as multivalent ions. The gelation induced by the strong short range electrostatic attractions is dramatically suppressed by the long range electrostatic correlations due to the charge along the non-cross-linked monomers and ions. A modified Debye-Hückel approach of cross-linked clusters of charged chains is used to determine the mean field gelation transition self-consistently. Highly dilute polyelectrolyte solutions tend to segregate macroscopically. Semidilute solutions can form gels if the Bjerrum length l(B) and the distance between neighboring charged monomers along the chain b are both greater than the ion size a.     

带电聚合电介质表面吸附理论研究

应用聚合电介质吸附的定标理论,根据介质和表面电介质常数的比率,考虑多化合价吸附电介质之间强相关性作用,我们提出一种表面排斥电荷的近似定标理论方法,根据这种方法把电介质表面吸附层的相图分为本质上不同的两大类。从相图可知：当表面电荷密度低（或体带相反电荷离子密度高）,这时表面和体带相反电荷离子密度几乎相同;一旦表面电荷密度足够高,就使带相反电荷的离子在表面上浓缩。据此,可确定在这个区域内,低化合价聚合电介质形成一个相关的多链状态,当化合价足够高时,由于近邻链之间的更强排斥增强,使状态转变成单链。     

带电聚合电介质表面吸附理论研究

摘要：应用聚合电介质吸附的定标理论,根据介质和表面电介质常数的比率,考虑多化合价吸附电介质之间强相关性作用,我们提出一种表面排斥电荷的近似定标理论方法,根据这种方法把电介质表面吸附层的相图分为本质上不同的两大类。从相图可知：当表面电荷密度低（或体带相反电荷离子密度高）,这时表面和体带相反电荷离子密度几乎相同;一旦表面电荷密度足够高,就使带相反电荷的离子在表面上浓缩。据此,可确定在这个区域内,低化合价聚合电介质形成一个相关的多链状态,当化合价足够高时,由于近邻链之间的更强排斥增强,使状态转变成单链。     

Non-linear osmotic brush regime: Simulations and mean-field theory

We investigate polyelectrolyte brushes in the osmotic regime using both theoretical analysis and molecular dynamics simulation techniques. In the simulations at moderate Bjerrum length, we observe that the brush height varies weakly with grafting density, in contrast to the accepted scaling law, which predicts a brush thickness independent of the grafting density. We show that such behavior can be explained by considering lateral electrostatic effects (within the non-linear Poisson-Boltzmann theory) combined with the coupling between lateral and longitudinal degrees of freedom due to the conserved polymer volume (which are neglected in scaling arguments). We also take the non-linear elasticity of polyelectrolyte chains into consideration, which makes significant effects as chains are almost fully stretched in the osmotic regime. It is shown that all these factors lead to a non-monotonic behavior for the brush height as a function of the grafting density. At large grafting densities, the brush height increases with increasing the grafting density due to the volume constraint. At small grafting densities, we obtain a re-stretching of the chains for decreasing grafting density, which is caused by lateral electrostatic contributions and is controlled by the counterion-condensation process around polyelectrolyte chains. These results are obtained assuming all counterions to be trapped within the brush, which is valid for sufficiently long chains of large charge fraction.Received: 14 May 2003, Published online: 11 November 2003PACS: 61.25.Hq Macromolecular and polymer solutions; polymer melts; swelling - 36.20.-r Macromolecules and polymer molecules - 61.20.Qg Structure of associated liquids: electrolytes, molten salts, etc.     

Collapse of spherical polyelectrolyte brushes in the presence of multivalent counterions

We consider the interaction of multivalent counterions with spherical polyelectrolyte brushes (SPB). The SPB result if linear polyelectrolyte (PE) chains (contour length: 60 nm) are densely grafted to colloidal spheres of 116 nm in diameter. Dispersed in water, the surface layer consisting of chains of the strong PE poly(styrene sulfonic acid) (PSS) will swell. We demonstrate that successive addition of trivalent ions (La3+) leads to a collapse in which the surface layer is shrinking drastically. All findings are discussed on the base of a theoretical mean-field approach using the Donnan equilibrium. The ion exchange and a strong binding of trivalent ions by PE chains is followed up by a drop in the osmotic pressure inside the brush. This reduction is the driving force for the collapse. The strong ion-chain correlation is discussed with results obtained from molecular dynamics simulations.     

Ordering of urchin-like charged copolymer micelles: Electrostatic, packing and polyelectrolyte correlations

We report the results of small-angle neutron scattering (SANS) studies on aqueous solutions of spherical polyelectrolyte micelles formed by association of charged-neutral diblock copolymers. The neutral moieties are found to self-assemble into small dense spheres (cores of the micelles) whose sizes are independent of the polymer concentration c. In the dilute regime, c<c ^*, where c^* is the overlap concentration of the micelles, the conformation of the charged groups, which form the corona of the micelles, is found to be extended. A liquid-like order is observed over a wide concentration range spanning from the dilute regime to the concentrated regime. For c>c ^*, polyelectrolyte correlations appear at smaller spatial scales and coexist with the liquid-like order. These results suggest that for dense brushes, above c^*, the rod-like statistics of the charged chains begin to disappear due to contraction of corona arms or by interpenetration of coronae. For less dense brushes, the charged chains are found to be extended up to concentrations far above c^*, before the progressive development of polyelectrolyte correlations. Received 8 October 1999     

Variational charge renormalization in charged systems

We apply general variational techniques to the problem of the counterion distribution around highly charged objects where strong condensation of counterions takes place. Within a field-theoretic formulation using a fluctuating electrostatic potential, the concept of surface-charge renormalization is recovered within a simple one-parameter variational procedure. As a test, we reproduce the Poisson-Boltzmann surface potential for a single-charged planar surface both in the weak-charge and strong-charge regime. We then apply our techniques to non-planar geometries where closed-form solutions of the non-linear Poisson-Boltzmann equation are not available. In the cylindrical case, the Manning charge renormalization result is obtained in the limit of vanishing salt concentration. However, for intermediate salt concentrations a slow crossover to the non-charged-renormalized regime (at high salt) is found with a quasi-power-law behavior which helps to understand conflicting experimental and theoretical results for the electrostatic persistence length of polyelectrolytes. In the spherical geometry charge renormalization is only found at intermediate salt concentrations, in agreement with previous numerical results.Received: 16 December 2002, Published online: 22 July 2003PACS: 82.70.-y Disperse systems; complex fluids - 61.20.Ja Computer simulation of liquid structure - 61.20.Qg Structure of associated liquids: electrolytes, molten salts, etc.     

Counterion distribution in a spherical charged sparse brush

The counterion distribution within a spherical polyelectrolyte sparse brush was measured by small-angle X-ray scattering using contrast variation with different counterions by means of ion dialysis. The brush was made by self-association of charged diblock copolymers. Thanks to the contrast variation method, we were able to separate the signal due to the monomers and the signal due to the counterions. At a small length scale, it is demonstrated that the system behaves as independent charged rods whose counterion distribution follows the Poisson-Boltzmann model. Received 14 February 2001 and Received in final form 2 May 2001     

Effect of counterions on the swelling of spherical polyelectrolyte brushes

We investigate the swelling of colloidal spherical polyelectrolyte brushes in the presence of different counterions. The colloidal particles consist of a solid poly(styrene) core of ca. 100 nm diameter onto which linear polyelectrolyte chains are chemically grafted. Two types of polyelectrolyte chains have been used here: The cationic polyelectrolyte poly(2-(acryloyl)ethyltrimethylammonium chloride)) (PATAC) and the anionic poly(styrenesulfonate) (PSS). Both systems are dispersed in water and the degree of swelling of the surface layer is studied by dynamic light scattering. Adding more and more salt leads to a strong shrinking of the surface layer as expected for polyelectrolyte brushes. It is shown that data obtained at low ionic strength can be collapsed on suitable master curves for monovalent and divalent counterions, respectively. For some ions, however, high salt concentrations may lead to a re-swelling of the brush layer in case of the cationic systems. This points to specific interactions of the counterions with the PATAC chains. This strong specific interaction between the counterions and the attached polyelectrolyte may even lead to flocculation of the particles at intermediate salt concentration. Surprisingly, for iodide and magnesium counterions the solubility increases again if the salt concentration is raised to 1 mol/l. Hence, specific interaction leads to salting-out effects as well as to salting-in effects for these colloidal particles. All specific effects seen at high concentrations of added salt can be explained by the increase of the reduced excluded-volume parameter which is due to the adsorption of salt ions.     

Dynamical clustering of counterions on flexible polyelectrolytes

Molecular dynamics simulations are used to study the spatiotemporal dynamics of charge fluctuations around a polyelectrolyte molecule at charge densities above and below the classic counterion condensation threshold. Surprisingly, the counterions form weakly interacting clusters which exhibit slowly decaying short range orientational order. Local charge fluctuations create energy fluctuations at the order of k_(B)T that is sufficient to affect the polyelectrolyte interaction with an approaching ligand molecule. The predictions of the classical theory appear to be appropriate only over much longer time scales.     

Reentrant behavior of divalent-counterion-mediated DNA-DNA electrostatic interaction

The problem of DNA-DNA interaction mediated by divalent counterions is studied using computer simulation. Although divalent counterions cannot condense free DNA molecules in solution, we show that if DNA configurational entropy is restricted, divalent counterions can cause DNA reentrant condensation similar to that caused by tri- or tetravalent counterions. DNA-DNA interaction is strongly repulsive at small or large counterion concentration and is negligible or slightly attractive for a concentration in between. Implications of our results to experiments of DNA ejection from bacteriophages are discussed. The quantitative result serves to understand electrostatic effects in other experiments involving DNA and divalent counterions.     

The role of excluded volume and electrostatics from coarse-grain modeling of the interaction of gemini surfactants with like-charged membranes

Cationic systems composed of lipids and/or surfactants are of paramount importance in a variety of applications. Within these, gemini have attracted particular attention, mainly due to their improved aggregation properties and to the possibility of tuning offered by the presence of a spacer. In this work, a Monte Carlo simulation study with a coarse-grained model was employed to assess the interaction of cationic gemini surfactants with a like-charged model membrane. Separating the contribution of the excluded volume and that of the electrostatic effects in the organization of gemini–lipid membranes was the first goal of this work and the role of these factors was assessed varying the concentration, the spacer length and the headgroup charge of gemini surfactants. The results provide a new insight on the organization of lipid headgroups in the vicinity of gemini surfactants. It was found that the surfactant–lipid interaction is strongly affected by the surfactant spacer length, being controlled by an overall balance between excluded volume and surfactant–lipid and surfactant–surfactant electrostatic effects. It is also seen that the out-of-plane motion of the spacer has a significant effect upon membrane organization and counterion condensation. Good agreement was found with results previously obtained from atomistic simulation.     

多价抗衡离子诱导聚甲基丙烯酸N,N-二甲基氨基乙酯(PDMAEMA)刷构象转变的UCST行为

本文应用分子场理论,研究pH、[Fe(CN)_6]~(3-)诱导聚甲基丙烯酸N,N-二甲基氨基乙酯(PDMAEMA)刷的上临界溶解温度(UCST)构象转变与结构特性.理论模型考虑p H和[Fe(CN)_6]~(3-)对PDMAEMA刷体系的静电调控作用.研究发现,在不同[Fe(CN)_6]~(3-)浓度、不同p H条件下,PDMAEMA刷呈现了UCST构象转变行为.由于p H调节PDMAEMA单体质子化,[Fe(CN)_6]~(3-)通过与PDMAEMA带正电荷的单体结合,形成了在PDMAEMA链内以[Fe(CN)_6]~(3-)为中介的带电单体间的静电吸引结合.随着温度升高,[Fe(CN)_6]~(3-)与PDMAEMA带正电荷的单体结合被破坏,[Fe(CN)_6]~(3-)在链内凝聚导致的静电屏蔽效应减弱,PDMAEMA链内带电单体间的静电排斥增强,PDMAEMA刷的构象呈现了从塌缩到溶胀的UCST转变行为,并且在较高[Fe(CN)_6]~(3-)浓度条件下,PDMAEMA刷构象转变的UCST增高.在较低p H值条件下,较多的PDMAEMA单体被质子化,[Fe(CN)_6]~(3-)与PDMAEMA带正电单体的结合增强,PDMAEMA刷构象转变的UCST增大.基于pH和[Fe(CN)_6]~(3-)对PDMAEMA刷体系中的静电调控效应,可以预言,在较小p H和较大[Fe(CN)_6]~(3-)浓度条件下,PDMAEMA链在垂直培基表面沿着链方向形成结节状结构.这是由于以[Fe(CN)_6]~(3-)为中介的链内带电单体间的静电吸引作用增强,导致临近单体间汇聚结节.我们的理论结果符合实验观测,由此表明,pH调节PDMAEMA单体的带电状态,以及[Fe(CN)_6]~(3-)在PDMAEMA链内凝聚导致的静电屏蔽效应,决定着PDMAEMA刷的UCST构象转变和结构特性.     

Adsorption of polyelectrolytes on silica and gold surfaces

The results of a study that helps understand the mechanisms of adsorption of polyelectrolytes on particles, using numerical simulation methods, specifically the one known as dissipative particle dynamics are reported here. The adsorption of cationic polyelectrolytes of two different polymerisation degrees interacting with two types of surfaces, one made of gold and the other of silica, is predicted and compared. We find that a more negatively charged wall does not necessarily adsorb more cationic polyelectrolytes because the electrostatic repulsion between the wall and the polyelectrolytes is stronger. Additionally, intra-chain repulsion plays an important role, because the largest polyelectrolyte chains have larger excluded volume than the shorter ones. In regard to the adsorption dependence on the polyelectrolyte polymerisation degree, we find that the excluded volume drives the adsorption throughout the intra-chain electrostatic repulsion, because the SiO₂ surface is strongly negative. These results are expected to be useful for several nanotechnological applications of current interest, such as in gene therapy and in the improvement of drug delivering mechanisms.