Computing hydrodynamic interactions in confined doubly periodic geometries in linear time.

We develop a linearly scaling variant of the force coupling method [K. Yeo and M. R. Maxey, J. Fluid Mech. 649, 205-231 (2010)] for computing hydrodynamic interactions among particles confined to a doubly periodic geometry with either a single bottom wall or two walls (slit channel) in the aperiodic direction. Our spectrally accurate Stokes solver uses the fast Fourier transform in the periodic xy plane and Chebyshev polynomials in the aperiodic z direction normal to the wall(s). We decompose the problem into two problems. The first is a doubly periodic subproblem in the presence of particles (source terms) with free-space boundary conditions in the z direction, which we solve by borrowing ideas from a recent method for rapid evaluation of electrostatic interactions in doubly periodic geometries [Maxian et al., J. Chem. Phys. 154, 204107 (2021)]. The second is a correction subproblem to impose the boundary conditions on the wall(s). Instead of the traditional Gaussian kernel, we use the exponential of a semicircle kernel to model the source terms (body force) due to the presence of particles and provide optimum values for the kernel parameters that ensure a given hydrodynamic radius with at least two digits of accuracy and rotational and translational invariance. The computation time of our solver, which is implemented in graphical processing units, scales linearly with the number of particles, and allows computations with about a million particles in less than a second for a sedimented layer of colloidal microrollers. We find that in a slit channel, a driven dense suspension of microrollers maintains the same two-layer structure as above a single wall, but moves at a substantially lower collective speed due to increased confinement.

Efficient dynamic simulations of charged dielectric colloids through a novel hybrid method.

Modern particle-based simulations increasingly incorporate polarization charges arising from spatially nonuniform permittivity. For complex dielectric geometries, calculation of these induced many-body effects typically requires numerical solvers based upon boundary-element methods, which very significantly increase the required computational effort. For the special case of dielectric spheres, such as colloids or nanoparticles, we recently proposed a semianalytical spectrally accurate hybrid method that combines the method of moments, the image-charge method, and the fast multipole method. The hybrid method is efficient and accurate even when dielectric spheres are closely packed. Here, we extend the method to the evaluation of direct and induced electrostatic forces and demonstrate how this can be incorporated in molecular dynamics simulations. The choice of the relevant numerical parameters for molecular dynamics simulations is discussed in detail, as well as comparisons to the boundary-element method. As a concrete example, we examine the challenging case of binary crystal structures composed of close-packed dielectric colloidal spheres.

Ion Structure Near a Core-Shell Dielectric Nanoparticle.

A generalized image charge formulation is proposed for the Green's function of a core-shell dielectric nanoparticle for which theoretical and simulation investigations are rarely reported due to the difficulty of resolving the dielectric heterogeneity. Based on the formulation, an efficient and accurate algorithm is developed for calculating electrostatic polarization charges of mobile ions, allowing us to study related physical systems using the Monte Carlo algorithm. The computer simulations show that a fine-tuning of the shell thickness or the ion-interface correlation strength can greatly alter electric double-layer structures and capacitances, owing to the complicated interplay between dielectric boundary effects and ion-interface correlations.

Hybrid Monte Carlo and continuum modeling of electrolytes with concentration-induced dielectric variations.

The distribution of ions near a charged surface is an important quantity in many biological and material processes, and has been therefore investigated intensively. However, few theoretical and simulation approaches have included the influence of concentration-induced variations in the local dielectric permittivity of an underlying electrolyte solution. Such local variations have long been observed and known to affect the properties of ionic solution in the bulk and around the charged surface. We propose a hybrid computational model that combines Monte Carlo simulations with continuum electrostatic modeling to investigate such properties. A key component in our hybrid model is a semianalytical formula for the ion-ion interaction energy in a dielectrically inhomogeneous environment. This formula is obtained by solving for the Green's function Poisson's equation with ionic-concentration-dependent dielectric permittivity using a harmonic interpolation method and spherical harmonic series. We also construct a self-consistent continuum model of electrostatics to describe the effect of ionic-concentration-dependent dielectric permittivity and the resulting self-energy contribution. With extensive numerical simulations, we verify the convergence of our hybrid simulation scheme, show the qualitatively different structures of ionic distribution due to the concentration-induced dielectric variations, and compare our simulation results with the self-consistent continuum model. In particular, we study the differences between weakly and strongly charged surfaces and multivalencies of counterions. Our hybrid simulations conform particularly the depletion of ionic concentrations near a charged surface and also capture the charge inversion. We discuss several issues and possible further improvement of our approach for simulations of large charged systems.

Effects of image charges, interfacial charge discreteness, and surface roughness on the zeta potential of spherical electric double layers.

We investigate the effects of image charges, interfacial charge discreteness, and surface roughness on spherical electric double layer structures in electrolyte solutions with divalent counterions in the setting of the primitive model. By using Monte Carlo simulations and the image charge method, the zeta potential profile and the integrated charge distribution function are computed for varying surface charge strengths and salt concentrations. Systematic comparisons were carried out between three distinct models for interfacial charges: (1) SURF1 with uniform surface charges, (2) SURF2 with discrete point charges on the interface, and (3) SURF3 with discrete interfacial charges and finite excluded volume. By comparing the integrated charge distribution function and the zeta potential profile, we argue that the potential at the distance of one ion diameter from the macroion surface is a suitable location to define the zeta potential. In SURF2 model, we find that image charge effects strongly enhance charge inversion for monovalent interfacial charges, and strongly suppress charge inversion for multivalent interfacial charges. For SURF3, the image charge effect becomes much smaller. Finally, with image charges in action, we find that excluded volumes (in SURF3) suppress charge inversion for monovalent interfacial charges and enhance charge inversion for multivalent interfacial charges. Overall, our results demonstrate that all these aspects, i.e., image charges, interfacial charge discreteness, their excluding volumes, have significant impacts on zeta potentials of electric double layers.

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

The polarization-induced charges of a dielectric sphere are studied for charged colloidal systems in electrolyte solutions with a primitive model. The method of constructing multiple-image charges is used to approximate the polarization potential of a microion outside the sphere; it is based on a numerical discretization of the potential's analytical integral representation, and can systematically approximate the exact potential with desired accuracy by varyiation of the number of point images. Different aspects of the image effects are then investigated by Monte Carlo simulations for several colloidal systems, in both salt-free and salty environments. Furthermore, we studied the influence of discrete surface charges of different valences, and demonstrate that the polarization charges can significantly strengthen charge reversal for the colloid-microion complex, especially for multivalent interfacial ions.