Ion condensation in salt‐free dilute polyelectrolyte solutions

Inter-DNA attraction mediated by divalent counterions

Can nonspecifically bound divalent counterions induce attraction between DNA strands? Here, we present experimental evidence demonstrating attraction between short DNA strands mediated by Mg2+ ions. Solution small angle x-ray scattering data collected as a function of DNA concentration enable model independent extraction of the second virial coefficient. As the [Mg2+] increases, this coefficient turns from positive to negative reflecting the transition from repulsive to attractive inter-DNA interaction. This surprising observation is corroborated by independent light scattering experiments. The dependence of the observed attraction on experimental parameters including DNA length provides valuable clues to its origin.

Direct observation of cations and polynucleotides explains polyion adsorption to like-charged surfaces

We show an experimental approach for directly observing the condensation of polynucleotides and their electrolyte counterions at a liquid/solid interface. X-ray standing waves (XSW) generated by Bragg diffraction from a d = 20 nm Si/Mo multilayer substrate are used to measure the distinct distribution profiles of the polyanions and simple cations along the surface normal direction with subnanometer resolution. The 1D spatial sensitivity of this approach is enhanced by observing the XSW induced fluorescence modulations over multiple orders of Bragg peaks. We study the interesting divalent cation driven adsorption of anionic polynucleotides to anionic surfaces by exposing a hydroxyl-terminated silica surface to an aqueous solution with ZnCl2 and mercurated poly-uridylic acid (a synthetic RNA molecule). The in situ long-period XSW measurements are used to follow the evolution of both the Zn and Hg distribution profiles during the adsorption process. The conditions and physical mechanisms that govern the observed divalent cation adsorption and subsequent polynucleotide adsorption to an anionic surface are explained by a thermodynamic model that incorporates nonlinear electrostatic effects.

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

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

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.

Structure and dynamics of water surrounding the poly(methacrylic acid): A molecular dynamics study

All-atom molecular dynamics simulations are used to study a single chain of poly(methacrylic acid) in aqueous solutions at various degrees of charge density. Through a combination of analysis on the radial distribution functions of water and snapshots of the equilibrated structure, we observe that local arrangements of water molecules, surrounding the functional groups of COO- and COOH in the chain, behave differently and correlated well to the resulting chain conformation behavior. In general, due to strong attractive interactions between water and charged COO- via the formation of hydrogen bonds, water molecules tend to form shell-like layers around the COO- groups. Furthermore, water molecules often act as a bridging agent between two neighboring COO- groups. These bridged water molecules are observed to stabilize the rodlike chain conformation that the highly charged chain reveals, as they significantly limit torsional and bending degrees of the backbone monomers. In addition, they display different dynamic properties from the bulk water. Both the resulting oxygen and hydrogen spectra are greatly shifted due to the presence of strong H-bonded interactions.

Scaling and universality in the counterion-condensation transition at charged cylinders

Counterions at charged rodlike polymers exhibit a condensation transition at a critical temperature (or, equivalently, at a critical linear charge density for polymers), which dramatically influences various static and dynamic properties of charged polymer solutions. We address the critical and universal aspects of this transition for counterions at a single charged cylinder in two and three spatial dimensions using numerical and analytical methods. By introducing a Monte Carlo sampling method in logarithmic radial scale, we are able to numerically simulate the critical limit of infinite system size (corresponding to the infinite-dilution limit) within tractable equilibration times. The critical exponents are determined for the inverse moments of the counterionic density profile (which play the role of the order parameters and represent the mean inverse localization length of counterions) both within mean-field theory and within Monte Carlo simulations. In three dimensions (3D), we demonstrate that correlation effects (neglected within mean-field theory) lead to an excessive accumulation of counterions near the charged cylinder below the critical temperature (i.e., in the condensation phase), while surprisingly, the critical region exhibits universal critical exponents in accordance with mean-field theory. Also in contrast with the typical trend in bulk critical phenomena, where fluctuations become more enhanced in lower dimensions, we demonstrate, using both numerical and analytical approaches, that mean-field theory becomes exact for the two-dimensional (2D) counterion-cylinder system at all temperatures (Manning parameters), when the number of counterions tends to infinity. For a finite number of particles, however, the 2D problem displays a series of peculiar singular points (with diverging heat capacity), which reflect successive delocalization events of individual counterions from the central cylinder. In both 2D and 3D, the heat capacity shows a universal jump at the critical point and the internal energy develops a pronounced peak. The asymptotic behavior of the energy peak location is used to determine the critical temperature, which is also found to be in agreement with the mean-field prediction.

Polynucleotide adsorption to negatively charged surfaces in divalent salt solutions

Polynucleotide adsorption to negatively charged surfaces via divalent ions is extensively used in the study of biological systems. We analyze here the adsorption mechanism via a self-consistent mean-field model that includes the pH effect on the surface-charge density and the interactions between divalent ions and surface groups. The adsorption is driven by the cooperative effect of divalent metal ion condensation along polynucleotides and their reaction with the surface groups. Although the apparent reaction constants are enhanced by the presence of polynucleotides, the difference between reaction constants of different divalent ions at the ideal condition explains why not all divalent cations mediate DNA adsorption onto anionic surfaces. Calculated divalent salt concentration and pH value variations on polynucleotide adsorption are consistent with atomic force microscope results. Here we use long-period x-ray standing waves to study the adsorption of mercurated-polyuridylic acid in a ZnCl2 aqueous solution onto a negatively charged hydroxyl-terminated silica surface. These in situ x-ray measurements, which simultaneously reveal the Hg and Zn distribution profiles along the surface normal direction, are in good agreement with our model. The model also provides the effects of polyelectrolyte line-charge density and monovalent salt on adsorption.

Optimal cell approach to osmotic properties of finite stiff-chain polyelectrolytes

We propose a self-consistent geometry optimized cell model approach to study osmotic properties of stiff-chain polyelectrolyte solutions. In contrast with the usual monotonic Poisson-Boltzmann prediction, the cell model predicts the correct nonmonotonic dependence of the osmotic coefficient on concentration. A lower degree of polymerization is found to reduce significantly the counterion condensation in a typical dilute strong polyelectrolyte. The results agree quantitatively with simulations of a corresponding many-body bulk system up to a dense semidilute regime.

Effective Charges of Polyelectrolytes in a Salt-Free Solution Based on Counterion Chemical Potential

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

Lattice Monte Carlo simulations of three-dimensional charged polymer chains

The configurational properties of strongly charged polyelectrolytes accompanied by neutralizing counterions in dilute solutions are simulated using the cooperative motion algorithm on the face-centered-cubic lattice. The full Coulomb potential and the excluded volume condition between different ions/beads are taken into account and the reduced temperature T* is considered the main, variable parameter. The calculations that have been carried out for solutions of both single and several chains indicate a few regions of their behavior: (1) for T*--> infinity, it corresponds to that of neutral, self-avoiding polymers under good solvent conditions; (2) for T* approximately 1, due to the electrostatic interactions being effectively stronger, the chains are more outstretched compared to their size at other temperatures; (3) for T* well below one, the counterion condensation becomes more and more dominant, which gradually leads to strongly collapsed chains; and (4) at the lowest temperatures the chains and counterions assume low-energy configurations in the form of neutral, compact aggregates.

Attraction of like-charged macroions in the strong-coupling limit

Like-charged macroions attract each other as a result of strong electrostatic correlations in the presence of multivalent counterions or at low temperatures. We investigate the effective electrostatic interaction between i) two like-charged rods and ii) two like-charged spheres using the recently introduced strong-coupling theory, which becomes asymptotically exact in the limit of large coupling parameter (i.e. for large counterion valency, low temperature, or high surface charge density on macroions). In contrast to previous applications of the strong-coupling theory, we deal with curved surfaces and an additional parameter, referred to as Manning parameter, is introduced, which measures the ratio between the radius of curvature of macroions to the Gouy-Chapman length. This parameter, together with the size of the confining box enclosing the two macroions and their neutralizing counterions, controls the counterion-condensation process that directly affects the effective interactions. For sufficiently large Manning parameters (weakly-curved surfaces), we find a strong long-ranged attraction between two macroions that form a closely-packed bound state with small surface-to-surface separation of the order of the counterion diameter in agreement with recent simulations results. For small Manning parameters (highly-curved surfaces), on the other hand, the equilibrium separation increases and the macroions unbind from each other as the confinement volume increases to infinity. This occurs via a continuous universal unbinding transition for two charged rods at a threshold Manning parameter of Epsilon c = 2/3, while the transition is strongly discontinuous for spheres because of a pronounced potential barrier at intermediate distances. Unlike the cylindrical case, the attractive forces between spheres disappear slowly for increasing confinement volume due to the complete de-condensation of counterions. Scaling arguments suggest that for moderate values of coupling parameter, strong-coupling predictions remain valid for sufficiently small surface-to-surface separations.

Electrical conductivity of polyelectrolyte solutions in the presence of added salt: The role of the solvent quality factor in light of a scaling approach

The effects of added salt on the electrical conductivity behavior of a polyelectrolyte solution are described in light of the scaling approach recently proposed by Dobrynin and Rubinstein [Macromolecules 28, 1859 (1995); 32, 915 (1999)], taking into account the influence of the solvent quality factor. The coupling between the conformation of the chain and the local charge distribution, giving rise to different conductometric behaviors, has been investigated under different conditions, in a wide concentration range of added salt. The polyion equivalent conductances lambda(p) have been evaluated in different concentration regimes for a hydrophilic polyion in good solvent condition and compared with the experimental values obtained from electrical conductivity measurements. The agreement is rather good in the wide range of concentration of the added salt investigated. In the case of poor solvent conditions, we find the appropriate expressions for the electrical conductivity when the polyion chain consists into collapsed beads alternating with stretched segments in the framework of the necklace globule model.

Monte carlo study of coulombic criticality in polyelectrolytes

The role of charges in determining the water solubility of polyelectrolytes, a question of considerable relevance to biology, is currently unresolved. We use computer simulations to study the purely Coulombic phase separation of flexible polyelectrolytes with monovalent counterions in an athermal solvent. In agreement with recent theories we find that the critical temperature for this transition increases with chain length, but that the critical density remains unchanged. We therefore stress that the phase behavior of polyelectrolytes is qualitatively different from uncharged polymers, where the critical density decreases towards zero for long chains.

Conformational instability of rodlike polyelectrolytes due to counterion fluctuations

The effective elasticity of highly charged stiff polyelectrolytes is studied in the presence of counterions, with and without added salt. The rigid polymer conformations may become unstable due to an effective attraction induced by counterion density fluctuations. Instabilities at the longest, or intermediate length scales, may signal collapse to globule, or necklace states, respectively. In the presence of added salt, a generalized electrostatic persistence length is obtained, which has a nontrivial dependence on the Debye screening length. It is also found that the onset of conformational instability is a reentrant phenomenon as a function of polyelectrolyte length for the unscreened case, and the Debye length or salt concentration for the screened case. This may be relevant in understanding the experimentally observed reentrant condensation of DNA.

Attractive interactions between rodlike polyelectrolytes: polarization, crystallization, and packing

We study the attractive interactions between rodlike charged polymers in solution that appear in the presence of multivalence counterions. The counterions condensed to the rods exhibit both a strong transversal polarization and a longitudinal crystalline arrangement. At short distances between the rods, the fraction of condensed counterions increases, and the majority of these occupy the region between the rods, where they minimize their repulsive interactions by arranging themselves into packing structures. The attractive interaction is strongest for multivalent counterions. Our model takes into account the hard-core volume of the condensed counterions, and their angular distribution around the rods. The hard-core constraint strongly suppresses longitudinal charge fluctuations.

Counterion condensation revisited

We review some of the characteristic properties of the structure of polyelectrolyte solutions: the condensed layer of counterions that forms abruptly at a critical threshold charge density on the polymer chain; the more diffuse Debye-Hückel cloud, which is spatially distinct from the condensed layer; and the entropic release of counterions from the condensed layer as a driving force for the binding of oppositely charged ligands. We present a reminder of the basis of our current understanding in a variety of experiments, simulations, and theories; and we attempt as well to clarify some misunderstandings. We present a new analysis of a lattice model that suggests why the limiting laws for polyelectrolyte thermodynamics have proved to be accurate despite the neglect of polymer-polymer interactions in their original derivation. We sketch recent progress in constructing a potential between counterion and polyion. A counterion located in the interface between condensed layer and Debye cloud is repelled from the polyion, creating a sharp boundary between the two counterion populations.