Multiple-image treatment of induced charges in Monte Carlo simulations of electrolytes near a spherical dielectric interface

Hydrodynamic Effects on the Collapse Kinetics of Flexible Polyelectrolytes

Understanding the collapse kinetics of polyelectrolytes (PEs) is crucial for comprehending various biological and industrial phenomena. Despite occurring in an aqueous environment, previous computational studies have overlooked the influence of hydrodynamic interactions (HIs) facilitated by fluid motion. Here, we directly compute the Navier-Stokes equation to investigate the collapse kinetics of a highly charged flexible PE. Our findings reveal that HI accelerates PE collapse induced by hydrophobicity and multivalent salt. In the case of hydrophobicity, HI induces long-range collective motion of monomers, accelerating the coarsening of local clusters through either Brownian-coagulation-like or evaporation-condensation-like processes, depending on the strength of hydrophobicity with respect to electrostatic interaction. Regarding multivalent salt, HI does not affect the condensation dynamics of multivalent ions but facilitates quicker movement of local dipolar clusters along the PE, thereby expediting the collapse process. These results provide valuable insights into the underlying mechanisms of HI in PE collapse kinetics.

From predictive modelling to machine learning and reverse engineering of colloidal self-assembly

An overwhelming diversity of colloidal building blocks with distinct sizes, materials and tunable interaction potentials are now available for colloidal self-assembly. The application space for materials composed of these building blocks is vast. To make progress in the rational design of new self-assembled materials, it is desirable to guide the experimental synthesis efforts by computational modelling. Here, we discuss computer simulation methods and strategies used for the design of soft materials created through bottom-up self-assembly of colloids and nanoparticles. We describe simulation techniques for investigating the self-assembly behaviour of colloidal suspensions, including crystal structure prediction methods, phase diagram calculations and enhanced sampling techniques, as well as their limitations. We also discuss the recent surge of interest in machine learning and reverse-engineering methods. Although their implementation in the colloidal realm is still in its infancy, we anticipate that these data-science tools offer new paradigms in understanding, predicting and (inverse) design of novel colloidal materials.

On the physics of both surface overcharging and charge reversal at heterophase interfaces

A series of Monte Carlo simulations are employed to reveal the physics of both surface overcharging and charge reversal at a negatively charged dielectric interface exposed to a bulk solution containing a +2:−1 electrolyte in the absence and presence of a monovalent salt.

Image-charge forces in thin interlayers due to surface charges in electrolyte

The surface forces arising in wetting films of nonpolar liquids or in thin air interlayers between an electrolyte and a nonpolar medium in the case of discrete charging of the dielectric-electrolyte interface are considered. The contributions of polarization effects to the distribution of the electrostatic potential in the three contacting media were calculated. Within the Debye-Hückel approximation, the analytical solutions were derived for the disjoining pressure in thin films, for the case of either dilute or relatively concentrated electrolyte solutions in the aforementioned systems. Analysis of the analytical and numerical results demonstrated that for dilute solutions the contribution of image forces to the disjoining pressure may significantly exceed the van der Waals forces for films from a few to tens of nanometers thick.

Self-energy-modified Poisson-Nernst-Planck equations: WKB approximation and finite-difference approaches

We propose a modified Poisson-Nernst-Planck (PNP) model to investigate charge transport in electrolytes of inhomogeneous dielectric environment. The model includes the ionic polarization due to the dielectric inhomogeneity and the ion-ion correlation. This is achieved by the self energy of test ions through solving a generalized Debye-Hückel (DH) equation. We develop numerical methods for the system composed of the PNP and DH equations. Particularly, toward the numerical challenge of solving the high-dimensional DH equation, we developed an analytical WKB approximation and a numerical approach based on the selective inversion of sparse matrices. The model and numerical methods are validated by simulating the charge diffusion in electrolytes between two electrodes, for which effects of dielectrics and correlation are investigated by comparing the results with the prediction by the classical PNP theory. We find that, at the length scale of the interface separation comparable to the Bjerrum length, the results of the modified equations are significantly different from the classical PNP predictions mostly due to the dielectric effect. It is also shown that when the ion self energy is in weak or mediate strength, the WKB approximation presents a high accuracy, compared to precise finite-difference results.

Electrostatic interaction in the presence of dielectric interfaces and polarization-induced like-charge attraction

Electrostatic polarization is important in many nano- and micro-scale physical systems such as colloidal suspensions, biopolymers, and nanomaterials assembly. The calculation of polarization potential requires an efficient algorithm for solving 3D Poisson's equation. We have developed a useful image charge method to rapid evaluation of the Green's function of the Poisson's equation in the presence of spherical dielectric discontinuities. This paper presents an extensive study of this method by giving a convergence analysis and developing a coarse-graining algorithm. The use of the coarse graining could reduce the number of image charges to around a dozen, by 1-2 orders of magnitude. We use the algorithm to investigate the interaction force between likely charged spheres in different dielectric environments. We find the size and charge asymmetry leads to an attraction between like charges, in agreement with existing results. Furthermore, we study three-body interactions and find, in the presence of an external interface, that the interaction force depends on the curvature of the interface and performs a nonmonotonic electrostatic force.

Simulation of charged systems in heterogeneous dielectric media via a true energy functional

For charged systems in heterogeneous dielectric media, a key obstacle for molecular dynamics (MD) simulations is the need to solve the Poisson equation in the media. This obstacle can be bypassed using MD methods that treat the local polarization charge density as a dynamic variable, but such approaches require access to a true free energy functional, one that evaluates to the equilibrium electrostatic energy at its minimum. In this Letter, we derive the needed functional. As an application, we develop a Car-Parrinello MD method for the simulation of free charges present near a spherical emulsion droplet separating two immiscible liquids with different dielectric constants. Our results show the presence of nonmonotonic ionic profiles in the dielectric with a lower dielectric constant.

Computational evidence of two driving mechanisms for overcharging in an electric double layer near the point of zero charge

We have adopted an ensemble Monte Carlo simulation method to systematically verify two physical driving mechanisms responsible for overcharging which refers to the adsorption of an effective charge onto a like-charged planar surface around the point of zero charge within the primitive model of mixed electrolytes with varying salt concentrations. One is electrostatic in character dominated by dielectric images and the other is purely entropic in origin by ionic size asymmetry effects, of which the former has never been reported both theoretically and experimentally and the latter could be interpreted satisfactorily in terms of available theoretical approaches. The electrostatically driven mechanism is found to critically depend on the ionic sizes while the entropically driven mechanism occurs with almost the same efficiency in a relative wide range of surface charge density. Depending on the delicate interplay between charge and steric correlations, the two distinct driving mechanisms may cooperatively give rise to a more pronounced overcharging process.

Competitive adsorption and ordered packing of counterions near highly charged surfaces: From mean-field theory to Monte Carlo simulations

Competitive adsorption of counterions of multiple species to charged surfaces is studied by a size-effect-included mean-field theory and Monte Carlo (MC) simulations. The mean-field electrostatic free-energy functional of ionic concentrations, constrained by Poisson's equation, is numerically minimized by an augmented Lagrangian multiplier method. Unrestricted primitive models and canonical ensemble MC simulations with the Metropolis criterion are used to predict the ionic distributions around a charged surface. It is found that, for a low surface charge density, the adsorption of ions with a higher valence is preferable, agreeing with existing studies. For a highly charged surface, both the mean-field theory and the MC simulations demonstrate that the counterions bind tightly around the charged surface, resulting in a stratification of counterions of different species. The competition between mixed entropy and electrostatic energetics leads to a compromise that the ionic species with a higher valence-to-volume ratio has a larger probability to form the first layer of stratification. In particular, the MC simulations confirm the crucial role of ionic valence-to-volume ratios in the competitive adsorption to charged surfaces that had been previously predicted by the mean-field theory. The charge inversion for ionic systems with salt is predicted by the MC simulations but not by the mean-field theory. This work provides a better understanding of competitive adsorption of counterions to charged surfaces and calls for further studies on the ionic size effect with application to large-scale biomolecular modeling.

A fast algorithm for treating dielectric discontinuities in charged spherical colloids

Electrostatic interactions between multiple colloids in ionic fluids are attracting much attention in studies of biological and soft matter systems. The evaluation of the polarization surface charges due to the spherical dielectric discontinuities poses a challenging problem to highly efficient computer simulations. In this paper, we propose a new method for fast calculating the electric field of spaced spheres using the multiple reflection expansion. The method uses a technique of recursive reflections among the spherical interfaces based on a formula of the multiple image representation, resulting in a simple, accurate and close-form expression of the surface polarization charges. Numerical calculations of the electric potential energies of charged spheres demonstrate the method is highly accurate with small number of reflections, and thus attractive for the use in practical simulations of related problems such as colloid suspension and macromolecular interactions.