An exact method to obtain effective electrostatic interactions from computer simulations: The case of effective charge amplification

Effective interactions, structure, and pressure in charge-stabilized colloidal suspensions: Critical assessment of charge renormalization methods

Charge-stabilized colloidal suspensions display a rich variety of microstructural and thermodynamic properties, which are determined by electro-steric interactions between all ionic species. The large size asymmetry between molecular-scale microions and colloidal macroions allows the microion degrees of freedom to be integrated out, leading to an effective one-component model of microion-dressed colloidal quasi-particles. For highly charged colloids with strong macroion-microion correlations, nonlinear effects can be incorporated into effective interactions by means of charge renormalization methods. Here, we compare and partially extend several practical mean-field methods of calculating renormalized colloidal interaction parameters, including effective charges and screening constants, as functions of concentration and ionic strength. Within the one-component description, we compute structural and thermodynamic properties from the effective interactions and assess the accuracy of the different methods by comparing predictions with elaborate primitive-model simulations [P. Linse, J. Chem. Phys. 113, 4359 (2000)]. We also compare various prescriptions for the osmotic pressure of suspensions in Donnan equilibrium with a salt ion reservoir and analyze instances where the macroion effective charge becomes larger than the bare one. The methods assessed include single-center cell, jellium, and multi-center mean-field theories. The strengths and weaknesses of the various methods are critically assessed, with the aim of guiding optimal and accurate implementations.

Anomalous Underscreening in the Restricted Primitive Model

Underscreening is a collective term for charge correlations in electrolytes decaying slower than the Debye length. Anomalous underscreening refers to phenomenology that cannot be attributed alone to steric interactions. Experiments with concentrated electrolytes and ionic fluids report anomalous underscreening, which so far has not been observed in simulation. We present Molecular Dynamics simulation results exhibiting anomalous underscreening that can be connected to cluster formation. A theory that accounts for ion pairing confirms the trend. Our results challenge the classic understanding of dense electrolytes impacting the design of technologies for energy storage and conversion.

A self-consistent Ornstein-Zernike jellium for highly charged colloids (microgels) in suspensions with added salt

This work discusses a jellium scheme, built within the framework of the multicomponent Ornstein-Zernike (OZ) equation, which is capable of describing the collective structure of suspensions of highly charged colloids with added salt, even in the presence of finite-size multivalent microions. This approach uses a suitable approximation to decouple the microion-microion correlations from the macroion-microion profiles, which in combination with the methodology from the dressed ion theory (DIT) gives a full account of the electrostatic effective potential among the colloids. The main advantages of the present contribution reside in its ability to manage the short-range potentials and non-linear correlations among the microions, as well as its realistic characterization of the ionic clouds surrounding each macroion. The structure factors predicted by this jellium scheme are contrasted with previously reported experimental results for microgel suspensions with monovalent salts (2019Phys. Rev. E100032602), thus validating its high accuracy in these situations. The present theoretical analysis is then extended to microgel suspensions with multivalent salts, which reveals the prominent influence of the counterion valence on the makeup of the effective potentials. Although the induced differences may be difficult to identify through the mesoscopic structure, our results suggest that the microgel collapsing transition may be used to enhance these distinct effects, thus giving a feasible experimental probe for these phenomena.

Primitive model electrolytes in the near and far field: Decay lengths from DFT and simulations

Inspired by recent experimental observations of anomalously large decay lengths in concentrated electrolytes, we revisit the Restricted Primitive Model (RPM) for an aqueous electrolyte. We investigate the asymptotic decay lengths of the one-body ionic density profiles for the RPM in contact with a planar electrode using classical Density Functional Theory (DFT) and compare these with the decay lengths of the corresponding two-body correlation functions in bulk systems, obtained in previous Integral Equation Theory (IET) studies. Extensive Molecular Dynamics (MD) simulations are employed to complement the DFT and IET predictions. Our DFT calculations incorporate electrostatic interactions between the ions using three different (existing) approaches: one is based on the simplest mean-field treatment of Coulomb interactions (MFC), while the other two employ the Mean Spherical Approximation (MSA). The MSAc invokes only the MSA bulk direct correlation function, whereas the MSAu also incorporates the MSA bulk internal energy. Although MSAu yields profiles that are in excellent agreement with MD simulations in the near field, in the far field, we observe that the decay lengths are consistent between IET, MSAc, and MD simulations, whereas those from MFC and MSAu deviate significantly. Using DFT, we calculated the solvation force, which relates directly to surface force experiments. We find that its decay length is neither qualitatively nor quantitatively close to the large decay lengths measured in experiments and conclude that the latter cannot be accounted for by the primitive model. The anomalously large decay lengths found in surface force measurements require an explanation that lies beyond primitive models.

Volume transition effects on the correlations and effective interactions among highly charged microgels

Recent experimental studies have demonstrated the huge influence that the volume phase transition (VPT) has on the collective structure of highly charged thermo-responsive microgels in aqueous solution with low concentrations of added monovalent salt, thus opening a promising new route for controlling the overall properties of practical colloidal suspensions. We present here an analysis of this structure based on the effective electrostatic potential obtained with the exact methodology of the dressed ion theory (DIT). Starting with a description at the primitive model level, we determine the correlations among the components of our model system (macroions plus monovalent anions and cations) by utilizing the two-density integral equation theory, thus allowing us to consider realistic values for the microgel charges. The resulting microgel structure factors show a good agreement with the reported light scattering measurements, whereas the microscopic pair distributions reveal that in this regime the shrunken states promote an enhanced counterion absorption into the microgels. This packing of counterions inside the microgels induces strongly non-linear correlations among the microions, and in turn provokes a substantial weakening of the microgel-microgel correlations. The ensuing effective interactions are then obtained by contracting the description to the level in which only the macroions are present. We find not only that the magnitude and reach of the corresponding pair potentials are markedly inhibited in the shrunken states, but also that their general form diverges from the conventional screened Coulomb shape. This makes it necessary to rethink the concepts of effective charge and screening length.

Accounting for effective interactions among charged microgels

We introduce a theoretical approach to describe structural correlations among charged permeable spheres at finite particle concentrations. This theory explicitly accounts for correlations among microions and between microions and macroions and allows for the proposal of an effective interaction among macroions that successfully captures structural correlations observed in poly-N-isopropyl acrylamide microgel systems. In our description the bare charge is fixed and independent of the microgel size, the microgel concentration, and the ionic strength, which contrasts with results obtained using linear response approximations, where the bare charge needs to be adapted to properly account for microgel correlations obtained at different conditions.

Nanoparticles confined to a spherical surface in the presence of an external field: Interaction forces and induced microstructure

The structural response of a set of charged nanoparticles confined to move on a spherical surface under the influence of an external field is studied by Brownian Dynamics (BD) simulations and by an integral equation approach (IEA). Considering an identical nanoparticle as the source of the external field, we analyze the force exerted by the N confined particles on the external one, as well as the corresponding potential energy, focusing on their dependence on the distance of the external particle to the center of the sphere r₀. The connection of the force and potential to the equilibrium local distribution of the adsorbed particles, that is, the microstructure within the spherical monolayer induced by the external nanoparticle, which is also dependent on r₀, is elucidated by this analysis. It is found that the external particle needs to surmount a considerable potential barrier when moving toward the spherical surface, although much smaller than the one generated by a uniform surface distribution with an equivalent amount of charge. This is understood in terms of the correlation hole within the confined monolayer induced by the external particle. Another interesting conclusion is that the IEA provides an accurate, almost quantitative, description of the main features observed in the BD results, yet it is much less computationally demanding. The connection of these results with the overall chemical equilibrium of charged surfactant nanoparticles in the context of Pickering emulsions is also briefly discussed.

Evidence of electrostatic-enhanced depletion attraction in the structural properties and phase behavior of binary charged colloidal suspensions

In this paper we study the structure and phase behavior of binary mixtures of charged particles at low ionic strength. Due to the large size asymmetry between both species, light scattering measurements give us access only to the partial static structure factor that corresponds to the big particles. We observe that the addition of small charged colloids produces a decrease of the main peak of the measured static structure factor and a shift to larger scattering vector values. This finding is in agreement with theory based on integral equations with the Hypernetted-Chain Closure (HNC) relation. The effective interaction between two big particles due to the presence of small particles is obtained by a HNC inversion scheme and used in numerical simulations that adequately reproduce the experiments. We find that the presence of small particles induces an electrostatic depletion screening among the big colloids, creating around them an exclusion zone for the small charged colloids greater than that caused in the case of neutral small colloids, which in turn augments the depletion effect.

Colloidal Stability in Asymmetric Electrolytes: Modifications of the Schulze-Hardy Rule

The Schulze-Hardy rule suggests a strong dependence of the critical coagulation concentration (CCC) on the ionic valence. This rule is addressed theoretically and confronted with recent experimental results. The commonly presented derivation of this rule assumes symmetric electrolytes and highly charged particles. Both assumptions are incorrect. Symmetric electrolytes containing multivalent ions are hardly soluble, and experiments are normally carried out with the well-soluble salts of asymmetric electrolytes containing monovalent and multivalent ions. In this situation, however, the behavior is completely different whether the multivalent ions represent the counterions or co-ions. When these ions represent the counterions, meaning that the multivalent ions have the opposite sign than the charge of the particle, they adsorb strongly to the particles. Thereby, they progressively reduce the magnitude of the surface charge with increasing valence. In fact, this dependence of the charge density on the counterion valence is mainly responsible for the decrease of the CCC with the valence. In the co-ion case, where the multivalent ions have the same sign as the charge of the particle, the multivalent ions are repelled from the particles, and the surfaces remain highly charged. In this case, the inverse Schulze-Hardy rule normally applies, whereby the CCC varies inversely proportional to the co-ion valence.

Electrical double layer properties of spherical oxide nanoparticles

The accurate characterization of the electrical double layer properties of nanoparticles is of fundamental importance for optimizing their physicochemical properties for specific biotechnological and biomedical applications. In this article, we use classical solvation density functional theory and a surface complexation model to investigate the effects of the pH and the nanoparticle size on the structural and electrostatic properties of an electrolyte solution surrounding a spherical silica oxide nanoparticle. The formulation has been particularly useful for identifying dominant interactions governing the ionic driving force at a variety of pH levels and nanoparticle sizes. As a result of the energetic interplay displayed between electrostatic potential, ion-ion correlation and particle crowding effects on the nanoparticle surface titration, rich, non-trivial ion density profiles and mean electrostatic potential behavior have been found.

Effective charges and virial pressure of concentrated macroion solutions

The stability of colloidal suspensions is crucial in a wide variety of processes, including the fabrication of photonic materials and scaffolds for biological assemblies. The ionic strength of the electrolyte that suspends charged colloids is widely used to control the physical properties of colloidal suspensions. The extensively used two-body Derjaguin-Landau-Verwey-Overbeek (DLVO) approach allows for a quantitative analysis of the effective electrostatic forces between colloidal particles. DLVO relates the ionic double layers, which enclose the particles, to their effective electrostatic repulsion. Nevertheless, the double layer is distorted at high macroion volume fractions. Therefore, DLVO cannot describe the many-body effects that arise in concentrated suspensions. We show that this problem can be largely resolved by identifying effective point charges for the macroions using cell theory. This extrapolated point charge (EPC) method assigns effective point charges in a consistent way, taking into account the excluded volume of highly charged macroions at any concentration, and thereby naturally accounting for high volume fractions in both salt-free and added-salt conditions. We provide an analytical expression for the effective pair potential and validate the EPC method by comparing molecular dynamics simulations of macroions and monovalent microions that interact via Coulombic potentials to simulations of macroions interacting via the derived EPC effective potential. The simulations reproduce the macroion-macroion spatial correlation and the virial pressure obtained with the EPC model. Our findings provide a route to relate the physical properties such as pressure in systems of screened Coulomb particles to experimental measurements.

Excluded volume and ion-ion correlation effects on the ionic atmosphere around B-DNA: theory, simulations, and experiments

The ionic atmosphere around a nucleic acid regulates its stability in aqueous salt solutions. One major source of complexity in biological activities involving nucleic acids arises from the strong influence of the surrounding ions and water molecules on their structural and thermodynamic properties. Here, we implement a classical density functional theory for cylindrical polyelectrolytes embedded in aqueous electrolytes containing explicit (neutral hard sphere) water molecules at experimental solvent concentrations. Our approach allows us to include ion correlations as well as solvent and ion excluded volume effects for studying the structural and thermodynamic properties of highly charged cylindrical polyelectrolytes. Several models of size and charge asymmetric mixtures of aqueous electrolytes at physiological concentrations are studied. Our results are in good agreement with Monte Carlo simulations. Our numerical calculations display significant differences in the ion density profiles for the different aqueous electrolyte models studied. However, similar results regarding the excess number of ions adsorbed to the B-DNA molecule are predicted by our theoretical approach for different aqueous electrolyte models. These findings suggest that ion counting experimental data should not be used alone to validate the performance of aqueous DNA-electrolyte models.

Specific ion effects induced by mono-valent salts in like charged aggregates in water

While salt mediated association between similarly charged poly-electrolytes occurs in a broad range of biological and colloidal systems, the effects of mono-valent salts remains little known experimentally. In this communication we systematically study influences of assorted mono-valent salts on structures of and interactions in two dimensional ordered bundles of charged fibrils assembled in water using Small Angle X-ray Scattering (SAXS). By quantitatively analyzing the scattering peak features, we discern two competing effects with opposite influences due to partitioning of salts in the aqueous complex. While electrostatic effects from salts residing between the fibrils suppress attraction between fibrils and expand the bundles, it is compensated by external osmotic pressure from peripheral salts in the aqueous media. The balance between the two effects varies for different salts and gives rise to ion-specific equilibrium behavior as well as structure of ordered bundles in salty water. The specific ions effects in like charged aggregates can be attributed to preferential distribution of ions inside or outside the bundles, correlated to the ranking of ions in Hofmeister series for macromolecules. Unlike conventional studies on Hofmeister effects by thermodynamic measurements relying on modeling for data interpretation, our study is based directly on structural analysis and is model-insensitive.

Static and dynamic screening effects in the electrostatic self-assembly of nano-particles

3D structure of self-assembled and oppositely charged nanoparticles.

Large counterions boost the solubility and renormalized charge of suspended nanoparticles

Colloidal particles are ubiquitous in biology and in everyday products such as milk, cosmetics, lubricants, paints, or drugs. The stability and aggregation of colloidal suspensions are of paramount importance in nature and in diverse nanotechnological applications, including the fabrication of photonic materials and scaffolds for biological assemblies, gene therapy, diagnostics, targeted drug delivery, and molecular labeling. Electrolyte solutions have been extensively used to stabilize and direct the assembly of colloidal particles. In electrolytes, the effective electrostatic interactions among the suspended colloids can be changed over various length scales by tuning the ionic concentration. However, a major limitation is gelation or flocculation at high salt concentrations. This is explained by classical theories, which show that the electrostatic repulsion among charged colloids is significantly reduced at high electrolyte concentrations. As a result, these screened colloidal particles are expected to aggregate due to short-range attractive interactions or dispersion forces as the salt concentration increases. We discuss here a robust, tunable mechanism for colloidal stability by which large counterions prevent highly charged nanoparticles from aggregating in salt solutions with concentrations up to 1 M. Large counterions are shown to generate a thicker ionic cloud in the proximity of each charged colloid, which strengthens short-range repulsions among colloidal particles and also increases the corresponding renormalized colloidal charge perceived at larger separation distances. These effects thus provide a reliable stabilization mechanism in a broad range of biological and synthetic colloidal suspensions.