Ion-ion correlations in highly asymmetrical electrolytes

Integral equation theory for mixtures of spherical and patchy colloids. 2. Numerical results

Thermodynamic properties and structure of binary mixtures of patchy and spherical colloids are studied using a recently developed theory [Y. V. Kalyuzhnyi, et al., Soft Matter, 2020, 16, 3456]. The theory is based on a solution of the multidensity Ornstein-Zernike equation and provides completely analytical expressions for the structure factors of these systems and for all their major thermodynamical quantities. The considered mixtures are made up of particles of different size and with a different number of patches. A set of molecular simulation data has been generated to enable a systematic comparison and to access thus accuracy of the theoretical predictions. In general, the predictions of the theory appear to be in good agreement with computer simulation data. For the models with a lower number of patches (n_p = 1, 2) the theoretical results show very good accuracy. Less accurate are the predictions for the four-patch versions of the model. While theoretical results for the radial distribution functions are, generally, relatively accurate for all the models, results for thermodynamics deteriorate with increasing concentration of the spherical colloids. Possible ways to improve the theory are briefly outlined.

Assessment of the Wolf method using the Stillinger-Lovett sum rules: From strong electrolytes to weakly charged colloidal dispersions

The Ewald method has been the cornerstone in molecular simulations for modeling electrostatic interactions of charge-stabilized many-body systems. In the late 1990s, Wolf and collaborators developed an alternative route to describe the long-range nature of electrostatic interactions; from a computational perspective, this method provides a more efficient and straightforward way to implement long-range electrostatic interactions than the Ewald method. Despite these advantages, the validity of the Wolf potential to account for the electrostatic contribution in charged fluids remains controversial. To alleviate this situation, in this contribution, we implement the Wolf summation method to both electrolyte solutions and charged colloids with moderate size and charge asymmetries in order to assess the accuracy and validity of the method. To this end, we verify that the proper selection of parameters within the Wolf method leads to results that are in good agreement with those obtained through the standard Ewald method and the theory of integral equations of simple liquids within the so-called hypernetted chain approximation. Furthermore, we show that the results obtained with the original Wolf method do satisfy the moment conditions described by the Stillinger-Lovett sum rules, which are directly related to the local electroneutrality condition and the electrostatic screening in the Debye-Hückel regime. Hence, the fact that the solution provided by the Wolf method satisfies the first and second moments of Stillinger-Lovett proves, for the first time, the reliability of the method to correctly incorporate the electrostatic contribution in charge-stabilized fluids. This makes the Wolf method a powerful alternative compared to more demanding computational approaches.

Integral equation theory for a mixture of spherical and patchy colloids: analytical description

An analytic theory for the structure and thermodynamics of two-component mixtures of patchy and spherical colloids is developed. The theory is based on an analytical solution of the multidensity Ornstein-Zernike equation supplemented by the associative Percus-Yevick closure relations. We derive closed-form analytic expressions for the partial structure factors and thermodynamic properties using the energy route for the model with arbitrary number of patches and any hard-sphere size ratio of the particles. To assess the accuracy of the theoretical predictions we compare them against existing and newly generated set of computer simulation data. In our numerical calculations we consider the model with equal hard-sphere sizes and one patch. Very good agreement between results of the theory and simulation for the pair correlation functions, excess internal energy and pressure is observed for almost all values of the system density, temperature and composition studied. Only in the region of low concentrations of spherical colloids the theoretical results become less accurate.

Long-range forces and charge inversions in model charged colloidal dispersions at finite concentration

In charged colloidal dispersion systems the interest is in finding their stability conditions, phase transitions, and transport properties, either in bulk or confinement, among other physicochemical quantities, for which the knowledge of the dispersions' molecular structure and the associated macroion-macroion forces is crucial. To investigate these phenomena simple models have been proposed. Most of the theoretical and simulation studies on charged particles suspensions are at infinite dilution conditions. Hence, these studies have been focused on the electrolyte structure around one or two isolated central particle(s), where phenomena as charge reversal, charge inversion and surface charge amplification have been shown to be relevant. However, experimental studies at finite volume fraction exhibit interesting phenomenology which imply very long-range correlations. A simple, yet useful, model is the Colloidal Primitive Model, in which the colloidal dispersion is modeled as a mixture of size (and charge) asymmetrical hard spheres, at finite volume fraction. In this paper we review recent integral equations solutions for this model, where very long-range attractive-repulsive forces, as well as new long-range, giant charge inversions are reported. The calculated macroions radial distribution functions, charge distributions, and macroion-macroion forces are qualitatively consistent with existing experimental results, and Monte Carlo and molecular dynamics simulations.

Integral Equation Theory of Biomolecules and Electrolytes

The so-called three-dimensional version (3D-RISM) can be used to describe the interactions of solvent components (here we treat water and ions) with a chemical or biomolecular solute of arbitrary size and shape. Here we give an overview of the current status of such models, describing some aspects of “pure” electrolytes (water plus simple ions) and of ionophores, proteins and nucleic acids in the presence of water and salts. Here we focus primarily on interactions with water and dissolved salts; as a practical matter, the discussion is mostly limited to monovalent ions, since studies of divalent ions present many difficult problems that have not yet been addressed. This is not a comprehensive review, but covers a few recent examples that illustrate current issues.

Transport coefficients of aqueous dodecyltrimethylammonium bromide solutions: Comparison between experiments, analytical calculations and numerical simulations

We study dynamical properties of ionic species in aqueous solutions of dodecyltrimethylammonium bromide, for several concentrations below and above the critical micellar concentration (cmc). New experimental determinations of the electrical conductivity are given which are compared to results obtained from an analytical transport theory; transport coefficients of ions in these solutions above the cmc are also computed from Brownian dynamics simulations. Analytical calculations as well as the simulation treat the solution within the framework of the continuous solvent model. Above the cmc, three ionic species are considered: the monomer surfactant, the micelle and the counterion. The analytical transport theory describes the structural properties of the electrolyte solution within the mean spherical approximation and assumes that the dominant forces which determine the deviations of transport processes from the ideal behavior (i.e., without any interactions between ions) are hydrodynamic interactions and electrostatic relaxation forces. In the simulations, both direct interactions and hydrodynamic interactions between solutes are taken into account. The interaction potential is modeled by pairwise repulsive 1/r(12) interactions and Coulomb interactions. The input parameters of the simulation (radii and self-diffusion coefficients of ions at infinite dilution) are partially obtained from the analytical transport theory which fits the experimental determinations of the electrical conductivity. Both the electrical conductivity of the solution and the self-diffusion coefficients of each species computed from Brownian dynamics are compared to available experimental data. In every case, the influence of hydrodynamic interactions (HIs) on the transport coefficients is investigated. It is shown that HIs are crucial to obtain agreement with experiments. In particular, the self-diffusion coefficient of the micelle, which is the largest and most charged species in the present system, is enhanced when HIs are included whereas the diffusion coefficients of the monomer and the counterion are roughly not influenced by HIs.