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

Charge Regulation Triggers Condensation of Short Oligopeptides to Polyelectrolytes

Electrostatic interactions between charged macromolecules are ubiquitous in biological systems, and they are important also in materials design. Attraction between oppositely charged molecules is often interpreted as if the molecules had a fixed charge, which is not affected by their interaction. Less commonly, charge regulation is invoked to interpret such interactions, i.e., a change of the charge state in response to a change of the local environment. Although some theoretical and simulation studies suggest that charge regulation plays an important role in intermolecular interactions, experimental evidence supporting such a view is very scarce. In the current study, we used a model system, composed of a long polyanion interacting with cationic oligolysines, containing up to 8 lysine residues. We showed using both simulations and experiments that while these lysines are only weakly charged in the absence of the polyanion, they charge up and condense on the polycations if the pH is close to the pKa of the lysine side chains. We show that the lysines coexist in two distinct populations within the same solution: (1) practically nonionized and free in solution; (2) highly ionized and condensed on the polyanion. Using this model system, we demonstrate under what conditions charge regulation plays a significant role in the interactions of oppositely charged macromolecules and generalize our findings beyond the specific system used here.

Cluster Formation in Solutions of Polyelectrolyte Rings

We use molecular dynamics simulations to explore concentrated solutions of semiflexible polyelectrolyte ring polymers, akin to the DNA mini-circles, with counterions of different valences. We find that the assembly of rings into nanoscopic cylindrical stacks is a generic feature of the systems, but the morphology and dynamics of such a cluster can be steered by the counterion conditions. In general, a small addition of trivalent ions can stabilize the emergence of clusters due to the counterion condensation, which mitigates the repulsion between the like-charged rings. Stoichiometric addition of trivalent ions can even lead to phase separation of the polyelectrolyte ring phase due to the ion-bridging effects promoting otherwise entropically driven clustering. On the other hand, monovalent counterions cause the formation of stacks to be re-entrant with density. The clusters are stable within a certain window of concentration, while above the window the polyelectrolytes undergo an osmotic collapse, disfavoring ordering. The cluster phase exhibits characteristic cluster glass dynamics with arrest of collective degrees of freedom but not the self-ones. On the other hand, the collapsed phase shows arrest on both the collective and single level, suggesting an incipient glass-to-glass transition, from a cluster glass of ring clusters to a simple glass of rings.

Interaction between two polyelectrolytes in monovalent aqueous salt solutions

We use the recently developed soft-potential-enhanced Poisson-Boltzmann (SPB) theory to study the interaction between two parallel polyelectrolytes (PEs) in monovalent ionic solutions in the weak-coupling regime. The SPB theory is fitted to ion distributions from coarse-grained molecular dynamics (MD) simulations and benchmarked against all-atom MD modelling for poly(diallyldimethylammonium) (PDADMA). We show that the SPB theory is able to accurately capture the interactions between two PEs at distances beyond the PE radius. For PDADMA positional correlations between the charged groups lead to locally asymmetric PE charge and ion distributions. This gives rise to small deviations from the SPB prediction that appear as short-range oscillations in the potential of mean force. Our results suggest that the SPB theory can be an efficient way to model interactions in chemically specific complex PE systems.

Counterion condensation theory for finite polyelectrolyte and salt concentrations

In the present work we analyze the physical fundamentals of Manning's counterion condensation using his charged line model in a simple salt solution. We extend the theory for the cases of finite saline concentration and polymeric concentration tending to zero and the case of both finite concentrations. To find the equilibrium between the phases of free and condensed counterions, besides minimizing the free energy, we deduce an auxiliary equation to determine the two characteristic parameters of the theory, the fraction of condensed counterions and the volume of condensation. We compare the obtained results in the present work for only one infinite charged line with the ones of counterion condensation theory by Schurr and Fujimoto. We find that the linear density of critical charge depends on the concentration of added salt and takes values higher than one, instead of the unitary value predicted by Manning. We obtain the equations by the activity and osmotic coefficients in function of the critical charge density. We compare them with the corresponding equations by Manning for these parameters. We extend the counterion condensation theory to solutions of linear polyelectrolytes for finite saline and polymeric concentrations using a cell model. We modify the electrostatic contribution to the Gibbs energy adding, to the traditional one calculated by Manning, the energy excess due to the macroion present in a cylindrical cell. We apply the theory to obtain the osmotic coefficient and we compare our results with experimental data of DNA osmotic coefficient and with theoretical adjustment using the Poisson-Boltzmann equation.

Modified Poisson–Boltzmann theory for polyelectrolytes in monovalent salt solutions with finite-size ions

We present a soft-potential-enhanced Poisson–Boltzmann (SPB) theory to efficiently capture ion distributions and electrostatic potential around rodlike charged macromolecules. The SPB model is calibrated with a coarse-grained particle-based model for polyelectrolytes (PEs) in monovalent salt solutions as well as compared to a full atomistic molecular dynamics simulation with the explicit solvent. We demonstrate that our modification enables the SPB theory to accurately predict monovalent ion distributions around a rodlike PE in a wide range of ion and charge distribution conditions in the weak-coupling regime. These include excess salt concentrations up to 1M and ion sizes ranging from small ions, such as Na+ or Cl−, to softer and larger ions with a size comparable to the PE diameter. The work provides a simple way to implement an enhancement that effectively captures the influence of ion size and species into the PB theory in the context of PEs in aqueous salt solutions.

Explicit characterization of counterion dynamics around a flexible polyelectrolyte

The article presents a comprehensive study of counterion dynamics around a generic linear polyelectrolyte chain with the help of coarse-grained computer simulations. The ion-chain coupling is discussed in the form of binding time, mean-square displacement (MSD) relative to the chain, local ion transport coefficient, and spatiotemporal correlations in the effective charge. We have shown that a counterion exhibits subdiffusive behavior 〈δR^{2}〉∼t^{δ}, δ≈0.9 w.r.t. chain's center of mass. The MSD of ions perpendicularly outward from the chain segment exhibits a smaller subdiffusive exponent compared to the one along the chain backbone. The effective diffusivity of ion is the lowest in chain's close proximity, extending up to the length-scale of radius of gyration R_{g}. Beyond R_{g} at larger distances, they attain diffusivity of free ion with a smooth cross-over from the adsorbed regime to the free ion regime. We have shown that the effective diffusivity drastically decreases for the multivalent ions, while the crossover length scale remains the same. Conversely, with increasing salt concentration the coupling-length scale reduces, while the diffusivity remains unaltered. The effective diffusivity of adsorbed-ion reveals an exponential reduction with electrostatic interaction strength. We further corroborate this from the binding time of ions on the chain, which also grows exponentially with the coupling strength of the ion-polymer duo. Moreover, the binding time of ions exhibits a weak dependence with salt concentration for the monovalent salt, while for multivalent salts the binding time decreases dramatically with concentration. Our work also elucidates fluctuations in the effective charge per site, where it exhibits strong negative correlations at short length-scales.

Active particles with polar alignment in ring-shaped confinement

We study steady-state properties of active, nonchiral and chiral Brownian particles with polar alignment and steric interactions confined within a ring-shaped confinement (annulus) in two dimensions. Exploring possible interplays between polar interparticle alignment, geometric confinement and the surface curvature, being incorporated here on minimal levels, we report a surface-population reversal effect, whereby active particles migrate from the outer concave boundary of the annulus to accumulate on its inner convex boundary. This contrasts the conventional picture, implying stronger accumulation of active particles on concave boundaries relative to the convex ones. The population reversal is caused by both particle alignment and surface curvature, disappearing when either of these factors is absent. We explore the ensuing consequences for the chirality-induced current and swim pressure of active particles and analyze possible roles of system parameters, such as the mean number density of particles and particle self-propulsion, chirality, and alignment strengths.

Poisson-Boltzmann equation with a random field for charged fluids

The classical Poisson-Boltzmann equation (CPBE), which is a mean field theory by averaging the ion fluctuation, has been widely used to study ion distributions in charged fluids. In this study, we derive a modified Poisson-Boltzmann equation with a random field from the field theory and recover the ion fluctuation through a multiplicative noise added in the CPBE. The Poisson-Boltzmann equation with a random field (RFPBE) captures the effect of the ion fluctuation and gives different ion distributions in the charged fluids compared to the CPBE. To solve the RFPBE, we propose a Monte Carlo method based on the path integral representation. Numerical results show that the effect of the ion fluctuation strengthens the ion diffusion into the domain and intends to distribute the ions in the fluid uniformly. The final ion distribution in the fluid is determined by the competition between the ion fluctuation and the electrostatic forces exerted by the boundaries. The RFPBE is general and feasible for high dimensional systems by taking the advantage of the Monte Carlo method. We use the RFPBE to study a two dimensional system as an example, in which the effect of ion fluctuation is clearly captured.

Charged nanorods at heterogeneously charged surfaces

We study the spatial and orientational distribution of charged nanorods (rodlike counterions) as well as the effective interaction mediated by them between two plane-parallel surfaces that carry fixed (quenched) heterogeneous charge distributions. The nanorods are assumed to have an internal charge distribution, specified by a multivalent monopolar moment and a finite quadrupolar moment, and the quenched surface charge is assumed to be randomly distributed with equal mean and variance on the two surfaces. While equally charged surfaces are known to repel within the traditional mean-field theories, the presence of multivalent counterions has been shown to cause attractive interactions between uniformly charged surfaces due to the prevalence of strong electrostatic couplings that grow rapidly with the counterion valency. We show that the combined effects due to electrostatic correlations (caused by the coupling between the mean surface field and the multivalent, monopolar, charge valency of counterions) as well as the disorder-induced interactions (caused by the coupling between the surface disorder field and the quadrupolar moment of counterions) lead to much stronger attractive interactions between two randomly charged surfaces. The interaction profile turns out to be a nonmonotonic function of the intersurface separation, displaying an attractive minimum at relatively small separations, where the ensuing attraction can exceed the maximum strong-coupling attraction (produced by multivalent monopolar counterions between uniformly charged surfaces) by more than an order of magnitude.

Rodlike Counterions near Charged Cylinders: Counterion Condensation and Intercylinder Interaction

We study a system composed of like-charged cylinders and dumbbell-like counterions, with the focus laid on the role of the internal structure of counterions, using Monte Carlo simulations. The dumbbell ions are found to exhibit novel counterion condensation behavior governed by their length. Effective electrostatic interactions mediated between charged parallel cylinders also turn out significantly different from the case of pointlike ions, as a result of the complex interplay between the spatially separated charge distribution in the dumbbell counterions, their orientation, and the curvature of the charged cylinder. We show that at a weak-to-moderate electrostatic coupling strength, where effective like-charge interactions are usually found to be repulsive, the intercylinder interaction can become attractive and display a distinctive sensitivity to the cylinder curvature and dumbbell size, proving the significant effect of ion structure.

Specific ion effects for polyelectrolytes in aqueous and non-aqueous media: the importance of the ion solvation behavior

We present the results of atomistic molecular dynamics simulations with regard to specific ion effects in water, methanol and N,N-dimethylacetamide (DMAc). As a reference system, we introduce rigid and rod-like models of polyanions and polycations in combination with alkali metal cations and halide anions as counterions. Pronounced specific ion effects can be observed in terms of the individual anion and cation condensation behavior. The outcomes of our simulations thus reveal significant deviations from standard electrostatic mean-field theories. A detailed investigation of the individual energy contributions shows that ion-dipole interactions play a pivotal role in rationalizing the findings. The corresponding deviations in terms of the cation and anion distribution can be brought into agreement with the donor and acceptor numbers of the solvents, which thus highlights the importance of solvent-ion interactions in addition to electrostatic attraction.

Ability of the Poisson-Boltzmann equation to capture molecular dynamics predicted ion distribution around polyelectrolytes

Here, we examine polyelectrolyte (PE) and ion chemistry specificity in ion condensation via all-atom molecular dynamics (MD) simulations and assess the ability of the Poisson-Boltzmann (PB) equation to describe the ion distribution predicted by the MD simulations. The PB model enables the extraction of parameters characterizing ion condensation. We find that the modified PB equation which contains the effective PE radius and the energy of the ion-specific interaction as empirical fitting parameters describes ion distribution accurately at large distances but close to the PE, especially when strongly localized charge or specific ion binding sites are present, the mean field description of PB fails. However, the PB model captures the MD predicted ion condensation in terms of the Manning radius and fraction of condensed counterions for all the examined PEs and ion species. We show that the condensed ion layer thickness in our MD simulations collapses on a single master curve for all the examined simple, monovalent ions (Na⁺, Br⁺, Cs⁺, Cl^-, and Br^-) and PEs when plotted against the Manning parameter (and consequently the PE line charge density). The significance of this finding is that, contrary to the Manning radius extracted from the mean field PB model, the condensed layer thickness in the all atom detail MD modelling does not depend on the PE chemistry or counterion type. Furthermore, the fraction of condensed counterions in the MD simulations exceeds the PB theory prediction. The findings contribute toward understanding and modelling ion distribution around PEs and other charged macromolecules in aqueous solutions, such as DNA, functionalized nanotubes, and viruses.

Local pH and effective pKA of weak polyelectrolytes - insights from computer simulations

In this work we study the titration behavior of weak polyelectrolytes by computer simulations. We analyze the local pH near the chains at various conditions and provide molecular-level insight which complements the recent experimental determination of this quantity. Next, we analyze the non-ideal titration behaviour of weak polyelectrolytes in solution, calculate the effective ionization constant and compare the simulation results with theoretical predictions. In contrast with the universal behaviour with respect to chain length, we find non-universality and deviations from theory with respect to polymer concentration and permittivity of the solvent. The latter we explain in terms of counterion condensation and ion correlation effects, which lead to reversal of the non-ideal titration behaviour at very low permittivities. We discuss the impact of these findings on the interpretation of experimental results.

50th Anniversary Perspective: A Perspective on Polyelectrolyte Solutions

From the beginning of life with the information-containing polymers until the present era of a plethora of water-based materials in health care industry and biotechnology, polyelectrolytes are ubiquitous with a broad range of structural and functional properties. The main attribute of polyelectrolyte solutions is that all molecules are strongly correlated both topologically and electrostatically in their neutralizing background of charged ions in highly polarizable solvent. These strong correlations and the necessary use of numerous variables in experiments on polyelectrolytes have presented immense challenges toward fundamental understanding of the various behaviors of charged polymeric systems. This Perspective presents the author's subjective summary of several conceptual advances and the remaining persistent challenges in the contexts of charge and size of polymers, structures in homogeneous solutions, thermodynamic instability and phase transitions, structural evolution with oppositely charged polymers, dynamics in polyelectrolyte solutions, kinetics of phase separation, mobility of charged macromolecules between compartments, and implications to biological systems.

Hidden Criticality of Counterion Condensation Near a Charged Cylinder

Counterion condensation onto a charged cylinder, known as the Manning transition, has received a great deal of attention since it is essential to understand the properties of polyelectrolytes in ionic solutions. However, the current understanding is still far from complete and poses a puzzling question: While the strong-coupling theory valid at large ionic correlations suggests a discontinuous nature of the counterion condensation, the mean-field theory always predicts a continuous transition at the same critical point. This naturally leads to a question how one can reconcile the mean-field theory with the strong-coupling prediction. Here, we study the counterion condensation transition on a charged cylinder via Monte Carlo simulations. Varying the cylinder radius systematically in relation to the system size, we find that in addition to the Manning transition, there exists a novel transition where all counterions are bound to the cylinder and the heat capacity shows a drop at a finite Manning parameter. A finite-size scaling analysis is carried out to confirm the criticality of the complete condensation transition, yielding the same critical exponents with the Manning transition. We show that the existence of the complete condensation is essential to explain how the condensation nature alters from continuous to discontinuous transition.

Scaling and criticality of the Manning transition

We consider a system consisting of a charged cylinder in the presence of neutralizing counterions. This system is well known to exhibit the Manning transition of counterion condensation onto the charged cylinder. We study the criticality and the scaling properties of the Manning transition, analyzing involved thermodynamic quantities such as condensed fraction, its fluctuation, and heat capacity. Through the Monte Carlo simulations and finite-size scaling analysis, we find that near the transition point the examined quantities exhibit scale-invariant behaviors with specific exponents, which provides an evidence that the Manning transition is a critical phenomenon having a scale-invariant property, analogous to bulk phase transitions. Furthermore, we numerically confirm that such scaling properties are not affected by the coupling strength.

Potential of mean force and transient states in polyelectrolyte pair complexation

The pair association between two polyelectrolytes (PEs) of the same size but opposite charge is systematically studied in terms of the potential of mean force (PMF) along their center-of-mass reaction coordinate via coarse-grained, implicit-solvent, explicit-salt computer simulations. The focus is set on the onset and the intermediate transient stages of complexation. At conditions above the counterion-condensation threshold, the PE association process exhibits a distinct sliding-rod-like behavior where the polymer chains approach each other by first stretching out at a critical distance close to their contour length, then "shaking hand" and sliding along each other in a parallel fashion, before eventually folding into a neutral complex. The essential part of the PMF for highly charged PEs can be very well described by a simple theory based on sliding charged "Debye-Hückel" rods with renormalized charges in addition to an explicit entropy contribution owing to the release of condensed counterions. Interestingly, at the onset of complex formation, the mean force between the PE chains is found to be discontinuous, reflecting a bimodal structural behavior that arises from the coexistence of interconnected-rod and isolated-coil states. These two microstates of the PE complex are balanced by subtle counterion release effects and separated by a free-energy barrier due to unfavorable stretching entropy.

Wang-Landau Reaction Ensemble Method: Simulation of Weak Polyelectrolytes and General Acid-Base Reactions

We present a novel method for the study of weak polyelectrolytes and general acid-base reactions in molecular dynamics and Monte Carlo simulations. The approach combines the advantages of the reaction ensemble and the Wang-Landau sampling method. Deprotonation and protonation reactions are simulated explicitly with the help of the reaction ensemble method, while the accurate sampling of the corresponding phase space is achieved by the Wang-Landau approach. The combination of both techniques provides a sufficient statistical accuracy such that meaningful estimates for the density of states and the partition sum can be obtained. With regard to these estimates, several thermodynamic observables like the heat capacity or reaction free energies can be calculated. We demonstrate that the computation times for the calculation of titration curves with a high statistical accuracy can be significantly decreased when compared to the original reaction ensemble method. The applicability of our approach is validated by the study of weak polyelectrolytes and their thermodynamic properties.

Poisson-Boltzmann thermodynamics of counterions confined by curved hard walls

We consider a set of identical mobile pointlike charges (counterions) confined to a domain with curved hard walls carrying a uniform fixed surface charge density, the system as a whole being electroneutral. Three domain geometries are considered: a pair of parallel plates, the cylinder, and the sphere. The particle system in thermal equilibrium is assumed to be described by the nonlinear Poisson-Boltzmann theory. While the effectively one-dimensional plates and the two-dimensional cylinder have already been solved, the three-dimensional sphere problem is not integrable. It is shown that the contact density of particles at the charged surface is determined by a first-order Abel differential equation of the second kind which is a counterpart of Enig's equation in the critical theory of gravitation and combustion or explosion. This equation enables us to construct the exact series solutions of the contact density in the regions of small and large surface charge densities. The formalism provides, within the mean-field Poisson-Boltzmann framework, the complete thermodynamics of counterions inside a charged sphere (salt-free system).

Counterion density profile around a charged disk: From the weak to the strong association regime

We present a comprehensive study of the two-dimensional one-component plasma in the cell model with charged boundaries. Starting from weak couplings through a convenient approximation of the interacting potential we were able to obtain an analytic formulation to the problem deriving the partition function, density profile, contact densities, and integrated profiles that compared well with the numerical data from Monte Carlo simulations. Additionally, we derived the exact solution for the special cases of Ξ=1,2,3,⋯, finding a correspondence between those from weak couplings and the latter. Furthermore, we investigated the strong-coupling regime taking into consideration the Wigner formulation. Elaborating on this, we obtained the profile to leading order, computed the contact density values as compared to those derived in an earlier work on the contact theorem. We formulated adequately the strong-coupling regime for this system that differed from previous formulations. Ultimately, we calculated the first-order corrections and compared them against numerical results from our simulations with very good agreement; these results compared equally well in the planar limit, whose results are well known.

Electrostatic correlations near charged planar surfaces

Electrostatic correlation effects near charged planar surfaces immersed in a symmetric electrolytes solution are systematically studied by numerically solving the nonlinear six-dimensional electrostatic self-consistent equations. We compare our numerical results with widely accepted mean-field (MF) theory results, and find that the MF theory remains quantitatively accurate only in weakly charged regimes, whereas in strongly charged regimes, the MF predictions deviate drastically due to the electrostatic correlation effects. We also observe a first-order like phase-transition corresponding to the counterion condensation phenomenon in strongly charged regimes, and compute the phase diagram numerically within a wide parameter range. Finally, we investigate the interactions between two likely-charged planar surfaces, which repulse each other as MF theory predicts in weakly charged regimes. However, our results show that they attract each other above a certain distance in strongly charged regimes due to significant electrostatic correlations.

Low-frequency collective exchange mode in the dielectric spectrum of salt-free dilute polyelectrolyte solutions

We study the frequency-dependent dielectric response of flexible polyelectrolyte chains in aqueous salt-free dilute solution as a function of polymer concentration and length by means of Brownian dynamics simulations including hydrodynamic interactions. We show that condensed and uncondensed counter ions are characterized by different dielectric response spectra and thus confirm long-held scaling assumptions used for the interpretation of experimental data. In addition to the well-known relaxation modes of condensed and uncondensed counter ions we observe for a single polyelectrolyte chain, we find in many-chain simulations a novel spectral feature at low frequencies that essentially corresponds to an exchange of counter ions between neighboring polyelectrolytes. Our results suggest that the experimental low-frequency dielectric mode might not be due to condensed counter ions, as is commonly assumed, but rather due to a collective many-chain process. This could resolve a long-standing puzzle in the comparison of experimental results and scaling predictions.

Ensemble equivalence for counterion condensation on a two-dimensional charged disk

We study the counterion condensation on a two-dimensional charged disk in the limit of infinite dilution, and compare the energy-temperature relation obtained from the canonical free energy and microcanonical entropy. The microcanonical entropy is piecewise linear in energy, and is shown to be concave for all energies. As a result, even though the interactions are long-ranged, the energy-temperature relation and hence the counterion condensation transition points are identical in both the ensembles.

Wigner-crystal formulation of strong-coupling theory for counterions near planar charged interfaces

We present a new analytical approach to the strong electrostatic coupling regime (SC) that can be achieved equivalently at low temperatures, high charges, low dielectric permittivity, etc. Two geometries are analyzed in detail: one charged wall first, and then two parallel walls at small distances that can be likely or oppositely charged. In all cases, only one type of mobile counterions is present, and ensures electroneutrality (salt-free case). The method is based on a systematic expansion around the ground state formed by the two-dimensional Wigner crystal(s) of counterions at the plate(s). The leading SC order stems from a single-particle theory, and coincides with the virial SC approach that has been much studied in the last 10 years. The first correction has the functional form of the virial SC prediction, but the prefactor is different. The present theory is free of divergences and the obtained results, both for symmetrically and asymmetrically charged plates, are in excellent agreement with available data of Monte Carlo simulations under strong and intermediate Coulombic couplings. All results obtained represent relevant improvements over the virial SC estimates. The present SC theory starting from the Wigner crystal and therefore coined Wigner SC, sheds light on anomalous phenomena like the counterion mediated like-charge attraction, and the opposite-charge repulsion.

A submicron device to rectify a square-wave angular velocity

We study a system composed of two thick dielectric disks separated by a thin layer of an electrolyte solution. Initially both plates have the same surface charge distribution. The surface charge distribution has no rotational symmetry. We show that the top plate experiences a torque [Formula: see text]([Formula: see text]) if it rotates about its axis by an angle [Formula: see text] . The torque can be controlled by varying the electrolyte concentration, the separation and the surface charge density of the plates. For a specific example of charged rods attached to the plates, we find [Formula: see text]([Formula: see text]) [Formula: see text] sin(4[Formula: see text]) . We also study the dynamics of the system. We consider the case where the angular velocity of the bottom disk is a square-wave signal. We find that the average angular velocity of the top disk is not zero.

Counter-ions at charged walls: two-dimensional systems

We study equilibrium statistical mechanics of classical point counter-ions, formulated on 2D Euclidean space with logarithmic Coulomb interactions (infinite number of particles) or on the cylinder surface (finite particle numbers), in the vicinity of a single uniformly charged line (one single double layer), or between two such lines (interacting double layers). The weak-coupling Poisson-Boltzmann theory, which applies when the coupling constant [Formula: see text] is small, is briefly recapitulated (the coupling constant is defined as [Formula: see text] [Formula: see text] [Formula: see text] e (2) , where [Formula: see text] is the inverse temperature, and e the counter-ion charge). The opposite limit ( [Formula: see text] [Formula: see text] ∞ is treated by using a recent method based on an exact expansion around the ground-state Wigner crystal of counter-ions. These two limiting results are compared at intermediary values of the coupling constant [Formula: see text] = 2[Formula: see text] ([Formula: see text] = 1, 2, 3) , to exact results derived within a 1D lattice representation of 2D Coulomb systems in terms of anti-commuting field variables. The models (density profile, pressure) are solved exactly for any particles numbers N at [Formula: see text] = 2 and up to relatively large finite N at [Formula: see text] = 4 and 6. For the one-line geometry, the decay of the density profile at asymptotic distance from the line undergoes a fundamental change with respect to the mean-field behavior at [Formula: see text] = 6 . The like-charge attraction regime, possible for large [Formula: see text] but precluded at mean-field level, survives for [Formula: see text] = 4 and 6, but disappears at [Formula: see text] = 2 .

Dressed counterions: strong electrostatic coupling in the presence of salt

We reformulate the theory of strong electrostatic coupling in order to describe an asymmetric electrolyte solution of monovalent salt ions and polyvalent counterions using field-theoretical techniques and Monte Carlo simulations. The theory is based on an asymmetric treatment of the different components of the electrolyte solution. The weak coupling Debye-Hückel approach is used in order to describe the monovalent salt ions while a strong coupling approach is used to tackle the polyvalent counterions. This combined weak-strong coupling approach effectively leads to dressed interactions between polyvalent counterions and thus directly affects the correlation attraction mediated by polyvalent counterions between like-charged objects. The general theory is specifically applied to a system composed of two uniformly charged plane-parallel surfaces in the presence of salt and polyvalent counterions. In the strong coupling limit for polyvalent counterions, the comparison with Monte Carlo simulations shows good agreement for large enough values of the electrostatic coupling parameter. We delineate two limiting laws that in fact encompass all the Monte Carlo data.

Salt-induced counterion-mobility anomaly in polyelectrolyte electrophoresis

We study the electrokinetics of a single polyelectrolyte chain in salt solution using hydrodynamic simulations. The salt-dependent chain mobility compares well with experimental DNA data. The mobility of condensed counterions exhibits a salt-dependent change of sign, an anomaly that is also reflected in the counterion excess conductivity. Using Green's function techniques this anomaly is explained by electrostatic screening of the hydrodynamic interactions between the chain and counterions.

Charge density coordination and dynamics in a rodlike polyelectrolyte

We study static and dynamic correlation functions of segmental charge density fluctuation on a rodlike polymer in ionic solutions. The polymer is described by an effective Hamiltonian which incorporates the Coulomb interaction between fluctuating charges screened by ionic fluid environment. We analytically calculate the correlation functions and discuss how charge distribution and dynamics are affected by counterion valency and concentration in electrolyte. We find that the charge correlation exhibits an underdamped oscillation with a wavelength comparable to the counterion size and with the amplitude increasing with the counterion valency. The dynamic correlation of charge density is shown to decay with a characteristic time varying with the counterion valency and concentration. The multivalency gives rise to faster decay of the correlation than that given by one dimensional diffusion.

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.

Brownian dynamics of double-stranded DNA in periodic systems with discrete salt

Numerical models of mesoscale DNA dynamics relevant to in vivo scenarios require methods that incorporate important features of the intracellular environment, while maintaining computational tractability. Because the explicit inclusion of ions leads to electrostatic calculations that scale as the square of the number of charged particles, such models typically handle these calculations using low-potential, mean-field approaches, rather than by considering the discrete interactions of ions. This allows approximation of the long-range, screened self-repulsion of DNA, but is unable to capture detailed electrostatic phenomena, such as short-range attractions mediated by ion-ion correlations. Here, we develop a dynamical model of explicitly double-stranded, sequence-specific DNA in a bulk environment consisting of other polyions and explicitly represented counterions and coions. DNA is represented as two interwound chains of charged Stokes spheres, and ions as free, monovalently charged Stokes spheres. Brownian dynamics simulations performed at salt concentrations of 0.1, 1, 10, and 100 mM demonstrate this model captures anticipated behaviors of the system, including increasing compaction of the polyion by the ionic atmosphere with increasing ionic strength. The decay of the distance dependence of the ion concentrations as one moves away from the polyion approaches their equilibrium values in quantitative agreement with predictions of Poisson-Boltzmann theory. The simulation results also demonstrate quantitative agreement with experimental measurements of the persistence length of B-DNA, which increases significantly at low ionic strengths. The model also captures behaviors intimating the importance of explicitly representing ionic and polyionic structure. These include penetration of the polyion interior by both coions and counterions, and counterion-mediated accumulation of coions near the surface of the polyion. Such phenomena are likely to play an important role in the formation of alternative DNA secondary structures, suggesting the present methods will prove valuable to dynamic models of superhelical stress-induced DNA structural transitions.

Counterion Condensation on Charged Spheres, Cylinders, and Planes

We use the framework of counterion condensation theory, in which deviations from linear electrostatics are ascribed to charge renormalization caused by collapse of counterions from the ion atmosphere, to explore the possibility of condensation on charged spheres, cylinders, and planes immersed in dilute solutions of simple salt. In the limit of zero concentration of salt, we obtain Zimm-Le Bret behavior: a sphere condenses none of its counterions regardless of surface charge density, a cylinder with charge density above a threshold value condenses a fraction of its counterions, and a plane of any charge density condenses all of its counterions. The response in dilute but nonzero salt concentrations is different. Spheres, cylinders, and planes all exhibit critical surface charge densities separating a regime of counterion condensation from states with no condensed counterions. The critical charge densities depend on salt concentration, except for the case of a thin cylinder, which exhibits the invariant criticality familiar from polyelectrolyte theory.