Perspective: Coulomb fluids—Weak coupling, strong coupling, in between and beyond

Charge regulation in ionic solutions: thermal fluctuations and Kirkwood-Schumaker interactions

We study the behavior of two macroions with dissociable charge groups, regulated by local variables such as pH and electrostatic potential, immersed in a monovalent salt solution, considering cases where the net charge can either change sign or remain of the same sign depending on these local parameters. The charge regulation in both cases is described by the proper free-energy function for each of the macroions, while the coupling between the charges is evaluated on the approximate Debye-Hückel level. The charge correlation functions and the ensuing charge fluctuation forces are calculated analytically and numerically. Strong attraction between like-charged macroions is found close to the point of zero charge, specifically due to asymmetric, anticorrelated charge fluctuations of the macroion charges. The general theory is then implemented for a system of two proteinlike macroions, generalizing the form and magnitude of the Kirkwood-Schumaker interaction.

Continuity of states between the cholesteric → line hexatic transition and the condensation transition in DNA solutions

A new method of finely temperature-tuning osmotic pressure allows one to identify the cholesteric → line hexatic transition of oriented or unoriented long-fragment DNA bundles in monovalent salt solutions as first order, with a small but finite volume discontinuity. This transition is similar to the osmotic pressure-induced expanded → condensed DNA transition in polyvalent salt solutions at small enough polyvalent salt concentrations. Therefore there exists a continuity of states between the two. This finding, together with the corresponding empirical equation of state, effectively relates the phase diagram of DNA solutions for monovalent salts to that for polyvalent salts and sheds some light on the complicated interactions between DNA molecules at high densities.

Experimental observation of the ion-ion correlation effects on charge inversion and strong adhesion between mica surfaces in aqueous electrolyte solutions

Direct force measurements between two mica surfaces in aqueous electrolyte solutions over broad ranges of LaCl3 concentrations and pH values were carried out with a surface forces apparatus. Charge inversion on mica surfaces is detected once the LaCl3 concentration reaches a critical value. With the continual increase of LaCl3 concentrations, the mica surface will be overscreened by the counterions. It is demonstrated that the two mica surfaces may experience the jump-in contact even at high LaCl3 concentrations, which is seldom seen in monovalent salt solutions. The strong adhesion cannot be attributed to the van der Waals force alone, but should include the ion-ion correlation forces. Through adjusting the pH values in LaCl3 solutions, the ion-ion correlation force can be evaluated quantitatively. These results provide important insight into the fundamental understanding in the role of ion-ion correlations in ion screening mechanism and interactions between charged objects.

Theory of non-equilibrium critical phenomena in three-dimensional condensed systems of charged mobile nanoparticles

A study of 3d electrostatic self-assembly (SA) in systems of charged nanoparticles (NPs) is one of the most difficult theoretical problems. In particular, the limiting case of negligible or very low polar media (e.g. salt) concentration, where the long-range NP interactions cannot be reduced to commonly used effective short-range (Yukawa) potentials, remains unstudied. Moreover, the present study has demonstrated that unlike the Debye-Hückel theory, a complete screening of the charges in SA kinetics (dynamic SA) is not always possible. Generally speaking, one has to take into account implicitly how each NP interacts with all other NPs (the true long-range interactions). Traditional theoretical methods allow us to monitor such electrostatic 3d system kinetics only for very short times, which is far from sufficient for understanding the dynamic SA. In this paper, combining an integrated analytical approach (the non-linear integro-differential kinetic equation for correlation functions) and reverse Monte Carlo in the 3d case, we have obtained a self-consistent solution of this challenging problem. We demonstrate, in particular, the existence of critical points and critical phenomena in the non-equilibrium kinetics in a 3d system of oppositely charged mobile NPs.

Nonlocal and nonlinear electrostatics of a dipolar Coulomb fluid

We study a model Coulomb fluid consisting of dipolar solvent molecules of finite extent which generalizes the point-like dipolar Poisson-Boltzmann model (DPB) previously introduced by Coalson and Duncan (1996 J. Phys. Chem. 100 2612) and Abrashkin et al (2007 Phys. Rev. Lett. 99 077801). We formulate a nonlocal Poisson-Boltzmann equation (NLPB) and study both linear and nonlinear dielectric response in this model for the case of a single plane geometry. Our results shed light on the relevance of nonlocal versus nonlinear effects in continuum models of material electrostatics.

Electronic structure, dielectric response, and surface charge distribution of RGD (1FUV) peptide

Long and short range molecular interactions govern molecular recognition and self-assembly of biological macromolecules. Microscopic parameters in the theories of these molecular interactions are either phenomenological or need to be calculated within a microscopic theory. We report a unified methodology for the ab initio quantum mechanical (QM) calculation that yields all the microscopic parameters, namely the partial charges as well as the frequency-dependent dielectric response function, that can then be taken as input for macroscopic theories of electrostatic, polar, and van der Waals-London dispersion intermolecular forces. We apply this methodology to obtain the electronic structure of the cyclic tripeptide RGD-4C (1FUV). This ab initio unified methodology yields the relevant parameters entering the long range interactions of biological macromolecules, providing accurate data for the partial charge distribution and the frequency-dependent dielectric response function of this peptide. These microscopic parameters determine the range and strength of the intricate intermolecular interactions between potential docking sites of the RGD-4C ligand and its integrin receptor.

Field-theoretic description of charge regulation interaction

In order to find the exact form of the electrostatic interaction between two proteins with dissociable charge groups in aqueous solution, we have studied a model system composed of two macroscopic surfaces with charge dissociation sites immersed in a counterion-only ionic solution. Field-theoretic representation of the grand canonical partition function is derived and evaluated within the mean-field approximation, giving the Poisson-Boltzmann theory with the Ninham-Parsegian boundary condition. Gaussian fluctuations around the mean field are then analyzed in the lowest-order correction that we calculate analytically and exactly, using the path integral representation for the partition function of a harmonic oscillator with time-dependent frequency. The first-order (one loop) free-energy correction gives the interaction free energy that reduces to the zero-frequency van der Waals form in the appropriate limit but in general gives rise to a monopolar fluctuation term due to charge fluctuation at the dissociation sites. Our formulation opens up the possibility to investigate the Kirkwood-Shumaker interaction in more general contexts where their original derivation fails.

Global analysis of the ground-state wrapping conformation of a charged polymer on an oppositely charged nano-sphere

We investigate the wrapping conformations of a single, strongly adsorbed polymer chain on an oppositely charged nano-sphere by employing a reduced (dimensionless) representation of a primitive chain-sphere model. This enables us to determine the global behavior of the chain conformation in a wide range of values for the system parameters including the chain contour length, its linear charge density and persistence length as well as the nano-sphere charge and radius, and also the salt concentration in the bathing solution. The structural behavior of a charged chain-sphere complex can be described in terms of a few distinct conformational symmetry classes separated by continuous or discontinuous transition lines which are determined by means of appropriately defined (order) parameters. Our results can be applied to a wide class of strongly coupled polymer-sphere complexes including, for instance, complexes that comprise a mechanically flexible or semiflexible polymer chain or an extremely short or long chain and, as a special case, include the biologically relevant example of DNA-histone complexes.

Effect of magnesium ions on dielectric relaxation in semidilute DNA aqueous solutions

The effect of magnesium ion Mg(2+) on the dielectric relaxation of semidilute DNA aqueous solutions has been studied by means of dielectric spectroscopy in the 100 Hz-100 MHz frequency range. de Gennes-Pfeuty-Dobrynin semidilute solution correlation length is the pertinent fundamental length scale for sufficiently low concentration of added salt, describing the collective properties of Mg-DNA solutions. No relaxation fingerprint of the DNA denaturation bubbles, leading to exposed hydrophobic core scaling, was detected at low DNA concentrations, thus indicating an increased stability of the double-stranded conformation in Mg-DNA solutions as compared to the case of Na-DNA solutions. Some changes are detected in the behavior of the fundamental length scale pertaining to the single molecule DNA properties, reflecting modified electrostatic screening effects of the Odijk-Skolnick-Fixman type. All results consistently demonstrate that Mg(2+) ions interact with DNA in a similar way as Na(1+) ions do, their effect being mostly describable through an enhanced screening.