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     

Virial expansion for charged colloids and electrolytes

Using a field-theoretic approach, we derive the first few coefficients of the exact low-density (“virial”) expansion of a binary mixture of positively and negatively charged hard spheres (two-component hard-core plasma, TCPHC). Our calculations are nonperturbative with respect to the diameters d₊ and d_- and charge valences q₊ and q_- of positive and negative ions. Consequently, our closed-form expressions for the coefficients of the free energy and activity can be used to treat dilute salt solutions, where typically d ₊∼d _- and q ₊∼q _-, as well as colloidal suspensions, where the difference in size and valence between macroions and counterions can be very large. We show how to map the TCPHC on a one-component hard-core plasma (OCPHC) in the colloidal limit of large size and valence ratio, in which case the counterions effectively form a neutralizing background. A sizable discrepancy with the standard OCPHC with uniform, rigid background is detected, which can be traced back to the fact that the counterions cannot penetrate the colloids. For the case of electrolyte solutions, we show how to obtain the cationic and anionic radii as independent parameters from experimental data for the activity coefficient. Received 6 September 2001 / Received in final form 20 May 2002 Published online 24 September 2002     

Effect of colloidal charge discretization in the primitive model

The effect of fixed discrete colloidal charges in the primitive model is investigated for spherical macroions. Instead of considering a central bare charge, as it is traditionally done, we distribute discrete charges randomly on the sphere. We use molecular dynamics simulations to study this effect on various properties such as overcharging, counterion distribution and diffusion. In the vicinity of the colloid surface the electrostatic potential may considerably differ from the one obtained with a central charge. In the strong Coulomb coupling, we showed that the colloidal charge discretization qualitatively influences the counterion distribution and leads to a strong colloidal charge-counterion pair association. However, we found that charge inversion still persists even if strong pair association is observed. Received 30 June 2000 and Received in final form 28 November 2000     

Electrophoresis of a charge-inverted macroion complex: Molecular-dynamics study

We have performed molecular-dynamics simulations to study the effect of an external electric field on a macroion in the solution of multivalent Z : 1 salt. To obtain plausible hydrodynamics of the medium, we explicitly make the simulation of many neutral particles along with ions. In a weak electric field, the macroion drifts together with the strongly adsorbed multivalent counterions along the electric field, in the direction proving inversion of the charge sign. The reversed mobility of the macroion is insensitive to the external field, and increases with salt ionic strength. The reversed mobility takes a maximal value at intermediate counterion valence. The motion of the macroion complex does not induce any flow of the neutral solvent away from the macroion, which reveals screening of hydrodynamic interactions at short distances in electrolyte solutions. A very large electric field, comparable to the macroion unscreened field, disrupts charge inversion by stripping the adsorbed counterions off the macroion. Received 5 December 2001 and Received in final form 10 April 2002     

Simple model for attraction between like-charged polyions

We present a simple model for the possible mechanism of appearance of attraction between like charged polyions inside a polyelectrolyte solution. The attraction is found to be short ranged, and exists only in the presence of multivalent counterions. It is produced by the correlations in layers of condensed counterions surrounding each polyion and is only weakly temperature dependent. We find the attraction to be maximum at zero temperature and dimish as the temperature is raised. The attraction is only possible if the number of condensed counterions exceeds the threshold, , where is the valence of counterions and Z is the polyion charge. Received 10 March 1999 and Received in final form 20 April 1999     

Monte Carlo simulations of star-branched polyelectrolyte micelles

The concentration profiles of monomers and counterions in star-branched polyelectrolyte micelles are calculated through Monte Carlo simulations, using the freely jointed chain model. We have investigated the onset of different regimes corresponding to the spherical and Manning condensation of counterions as a function of the strength of the Coulomb coupling. The Monte Carlo results are in fair agreement with the predictions of Self-Consistent-Field analytical models. We have simulated a real system of diblock copolymer micelles of (sodium-polystyrene-sulfonate)(NaPSS)-(polyethylene-propylene)(PEP) with f = 54 hydrophilic branches of N = 251 monomers at room temperature in salt-free solution. The calculated form factor compares nicely with our neutron scattering data. Received 18 July 2002 and Received in final form 11 October 2002 RID="a" ID="a"e-mail: roger@drecam.saclay.cea.fr     

Adsorption of hydrophobic polyelectrolytes onto oppositely charged surfaces

We present a scaling theory for the adsorption of a weakly charged polyelectrolyte chain in a poor solvent onto an oppositely charged surface. Depending on the fraction of charged monomers and on the solvent quality for uncharged monomers, the globule in the bulk of the solution has either a spherical conformation or a necklace structure. At sufficiently high surface charge density, a chain in the globular conformation adsorbs in a flat pancake conformation due to the Coulombic attraction to the oppositely charged surface. Different adsorption regimes are predicted depending on two screening lengths (the Debye screening length monitored by the salt concentration and the Gouy-Chapman length monitored by the surface charge density), on the degree of ionization of the polymer and on the solvent strength. At low bulk ionic strength, an increase in the surface charge density may induce a transition from an adsorbed necklace structure to a uniform pancake due to the enhanced screening of the intra-chain Coulombic repulsion by the counterions localized near the surface. Received 12 April 2001     

Correlation effects,image charge effects and finite size in the macro-ion-electrolyte system: A field-theoretic approach

We consider a model of a macro-ion surrounded by small ions of an electrolyte solution. The finite size of ionic charge distributions of ions, and image charge effects are considered. From such a model it is possible to construct a statistical field theory with a single fluctuating field and derive physical interpretations for both the mean field and two-point correlation function. For point-like charges, at the level of a Gaussian (or saddle point) approximation, we recover the standard Poisson-Boltzmann equation. However, to include ionic correlation effects, as well as image charge effects of individual ions, we must go beyond this. From the field theory considered, it is possible to construct self-consistent approximations. We consider the simplest of these, namely the Hartree approximation. The Hartree equations take the form of two coupled equations. One is a modified Poisson-Boltzmann equation; the other describes both image charge effects on the individual ions, as well as correlations. Such equations are difficult to solve numerically, so we develop an (a WKB-like) approximation for obtaining approximate solutions. This, we apply to a uniformly charged rod in univalent electrolyte solution, for point like ions, as well as for extended spherically symmetric distributions of ionic charge on electrolyte ions. The solutions show how correlation effects and image charge effects modify the Poisson-Boltzmann result. Finite-size charge distributions of the ions reduce both the effects of correlations and image charge effects. For point charges, we test the WKB approximation by calculating a leading-order correction from the exact Hartree result, showing that the WKB-like approximation works reasonably well in describing the full solution to the Hartree equations. From these solutions, we also calculate an effective charge compensation parameter in an analytical formula for the interaction of two charged cylinders. Electronic supplementary material  Supplementary material in the form of a doc file available from the journal web page at and are accessible for authorised users.     

Surfaces with quenched and annealed disordered charge distributions

We consider surfaces with disordered charge distribution. The disorder can be caused by mobile charges, as for example in mixed lipid bilayers, or by weakly charged surfaces where charge regulation takes place (e.g. carboxyl groups). Using Monte-Carlo simulation methods we find for quenched as well as annealed disordered charge distributions counterion densities close to the surface that are significantly larger than for ordered regularly spaced surface ions. Our field-theoretic results agree well with results obtained from Monte-Carlo simulations of the system. Furthermore, we obtain expressions for the effective interaction between charged colloids and charged rods close to a charged surface and discuss the effect of the surface-ion mobility and polarization charges on the interaction. In general, polarization effects as well as surface-ion mobility lead to a weakening of the effective interaction between charged objects.     

Static and dynamic properties of supercooled thin polymer films

The dynamic and static properties of a supercooled (non-entangled) polymer melt are investigated via molecular-dynamics (MD) simulations. The system is confined between two completely smooth and purely repulsive walls. The wall-to-wall separation (film thickness), D, is varied from about 3 to about 14 times the bulk radius of gyration. Despite the geometric confinement, the supercooled films exhibit many qualitative features which were also observed in the bulk and could be analyzed in terms of mode-coupling theory (MCT). Examples are the two-step relaxation of the incoherent intermediate scattering function, the time-temperature superposition property of the late time α-process and the space-time factorization of the scattering function on the intermediate time scale of the MCT β-process. An analysis of the temperature dependence of the α-relaxation time suggests that the critical temperature, T _c, of MCT decreases with D. If the confinement is not too strong ( D≥10monomer diameter), the static structure factor of the film coincides with that of the bulk when compared for the same distance, T - T _c(D), to the critical temperature. This suggests that T - T _c(D) is an important temperature scale of our model both in the bulk and in the films. Received 12 September 2001     

The distribution of counterions around synthetic rod-like polyelectrolytes in solution

An investigation of the radial distribution of the counterions of a synthetic rodlike polyelectrolyte in aqueous solution is presented. The cationic polyelectrolyte used here has a poly(p-phenylene) backbone. For typical molecular weights the macroion comprises approximately one persistence length (ca. 20 nm) and effects of finite stiffness may be disregarded. Each repeating unit bears four charges which leads to a charge parameter of ξ = 6.65. The distribution of the iodide counterions around this highly charged macroion is studied by small-angle X-ray scattering (SAXS) in dilute aqueous solution. These investigations are supplemented by measurements using anomalous small-angle X-ray scattering (ASAXS) that furnishes additional information about the contrast of the macroion. Data taken at high scattering angles give indication for contributions caused by the longitudinal fluctuations of the counterions. After correction for this effect the experimental results are compared to intensities calculated by use of the Poisson-Boltzmann (PB)-cell model. It is found that the PB-cell model describes the corrected data at intermediate and high scattering angles. Deviations at low scattering angle are attributed to the mutual interaction of the rod-like polyelectrolyte that can be described in terms of an effective structure factor. Data taken at lowest scattering angles point to a weak attraction between the rod-like macroions. Received: 27 July 2001 and Received in final form 27 March 2002     

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.     

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.     

Enthalpy recovery of a glass-forming liquid constrained in a nanoporous matrix: Negative pressure effects

The T _g of organic liquids confined to nanoporous matrices and that of thin polymer films can decrease dramatically from the bulk value. One possible explanation for this phenomenon is the development of hydrostatic tension during vitrification under confinement that results in a concomitant increase in the free volume. Here we present experimental evidence and modeling results for ortho-terphenyl (o-TP) confined in pores as small as 11.6 nm that indicate that, although there is an important hydrostatic tension in the liquid in the pores, it does not develop until near the reduced T _g of the constrained material --well below the bulk T _g. Enthalpy recovery for the o-TP in the nanopores exhibits accelerated physical aging relative to the bulk, as well as a leveling off of the fictive temperature at equilibrium values greater than the aging temperature. An adaptation of the structural recovery model that incorporates vitrification under isochoric conditions is able to provide a quantitative explanation for the apparently anomalous aging observed in nanopore confined liquids and in thin polymeric films. The results strongly support the existence of an intrinsic size effect as the cause of the reduced T _g. Received 3 September 2001     

Flexible linear polyelectrolytes in multivalent salt solutions: Solubility conditions

Single- and double-stranded DNA and many biological and synthetic polyelectrolytes undergo two structural transitions upon increasing the concentration of multivalent salt or molecules. First, the expanded-stretched chains in low monovalent salt solutions collapse into nearly neutral compact structures when the density of multivalent salt approaches that of the monomers. With further addition of multivalent salt the chains redissolve acquiring expanded-coiled conformations. We study the redissolution transition using a two-state model (F.J. Solis, M. Olvera de la Cruz, J. Chem. Phys. 112, 2030 (2000)). The redissolution occurs when there is a high degree of screening of the electrostatic interactions between monomers, thus reducing the energy of the expanded state. The transition is determined by the chemical potential of the multivalent ions in the solution, μ and the inverse screening length, κ. The transition point also depends on the charge distribution along the chain but is nearly independent of the molecular weight and degree of flexibility of the polyelectrolytes. We generate a diagram of μversusκ² where we find two regions of expanded conformations, one with charged chains and the other with overcharged (inverted charge) chains, separated by a collapsed nearly neutral conformation region. The collapse and redissolution transitions occur when the trajectory of the properties of the salt crosses the boundaries between these regions. We find that in most cases the redissolution occurs within the same expanded branch from which the chain precipitates. Received 15 May 2000 and Received in final form 28 June 2000     

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

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

Direct observation of counterion organization in F-actin polyelectrolyte bundles

Attractions between like-charged polyelectrolytes have been observed in a variety of systems (W.M. Gelbart, R.F. Bruinsma, P.A. Pincus, V.A. Parsegian, Phys. Today 53, September issue, 38 (2000)). Recent biological examples include DNA, filamentous viruses, and F-actin. Theoretical investigations on idealized systems indicate that counterion correlations play a central role, but no experiments that specifically probe such correlations have been performed. Using synchrotron X-ray diffraction, we have directly observed the organization of multivalent ions on cytoskeletal filamentous actin (a well-defined biological polyelectrolyte) and found an unanticipated symmetry-breaking collective counterion mechanism for generating attractions. Surprisingly, the counterions do not form a lattice that simply follows actins helical symmetry; rather, the counterions organize into frozen ripples parallel to the actin filaments and form structures reminiscent of charge density waves. Moreover, these 1D counterion charge density waves form a coupled mode with twist deformations of the oppositely charged actin filaments. This counterion organization is not sensitive to thermal fluctuations in temperature range accessible to protein-based polyelectrolyte systems. Moreover, the counterion density waves are pinned to the spatial periodicity of charges on the actin filament even if the global filament charge density is varied, indicating the importance of charge periodicity on the polyelectrolyte substrate.     

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.     

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

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

High frequency sound and the boson peak in amorphous silica

We present the results of extensive molecular dynamics computer simulations in which the high frequency dynamics of silica, i.e. for frequencies ν > 0.5 THz, is investigated in the viscous liquid state as well as in the glass state. We characterize the properties of high frequency sound modes by analyzing J _l(q,ν) and J _t(q,ν), the longitudinal and transverse current correlation function, respectively. For wave-vectors q > 0.4 ?^-1 the spectra are sitting on top of a flat background. The dynamic structure factor S(q,ν) exhibits for q > 0.23 ?^-1 a boson peak which is located nearly independent of q around 1.7 THz and for which the intensity scales approximately linearly with temperature. We show that the low frequency part of the boson peak is mainly due to the elastic scattering of transverse acoustic modes with frequencies around 1 THz. The strength of this scattering depends on q and is largest around q = 1.7 ?^-1, the location of the first sharp diffraction peak in the static structure factor. By studying S(q,ν) for different system sizes we show that strong finite size effects are present in the low frequency part of the boson peak in that for small systems part of its intensity is missing. We discuss the consequences of these finite size effects for the structural relaxation. Received 27 June 2000 and Received in final form 9 January 2001