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.

Calculating particle pair potentials from fluid-state pair correlations: Iterative ornstein-zernike inversion.

An iterative Monte Carlo inversion method for the calculation of particle pair potentials from given particle pair correlations is proposed in this article. The new method, which is best referred to as Iterative Ornstein-Zernike Inversion, represents a generalization and an improvement of the established Iterative Boltzmann Inversion technique (Reith, Pütz and Müller-Plathe, J. Comput. Chem. 2003, 24, 1624). Our modification of Iterative Boltzmann Inversion consists of replacing the potential of mean force as an approximant for the pair potential with another, generally more accurate approximant that is based on a trial bridge function in the Ornstein-Zernike integral equation formalism. As an input, the new method requires the particle pair correlations both in real space and in the Fourier conjugate wavenumber space. An accelerated iteration method is included in the discussion, by which the required number of iterations can be greatly reduced below that of the simple Picard iteration that underlies most common implementations of Iterative Boltzmann Inversion. Comprehensive tests with various pair potentials show that the new method generally surpasses the Iterative Boltzmann Inversion method in terms of reliability of the numerical solution for the particle pair potential. © 2018 Wiley Periodicals, Inc.

Short- and long-time diffusion and dynamic scaling in suspensions of charged colloidal particles.

We report on a comprehensive theory-simulation-experimental study of collective and self-diffusion in concentrated suspensions of charge-stabilized colloidal spheres. In theory and simulation, the spheres are assumed to interact directly by a hard-core plus screened Coulomb effective pair potential. The intermediate scattering function, fc(q, t), is calculated by elaborate accelerated Stokesian dynamics (ASD) simulations for Brownian systems where many-particle hydrodynamic interactions (HIs) are fully accounted for, using a novel extrapolation scheme to a macroscopically large system size valid for all correlation times. The study spans the correlation time range from the colloidal short-time to the long-time regime. Additionally, Brownian Dynamics (BD) simulation and mode-coupling theory (MCT) results of fc(q, t) are generated where HIs are neglected. Using these results, the influence of HIs on collective and self-diffusion and the accuracy of the MCT method are quantified. It is shown that HIs enhance collective and self-diffusion at intermediate and long times. At short times self-diffusion, and for wavenumbers outside the structure factor peak region also collective diffusion, are slowed down by HIs. MCT significantly overestimates the slowing influence of dynamic particle caging. The dynamic scattering functions obtained in the ASD simulations are in overall good agreement with our dynamic light scattering (DLS) results for a concentration series of charged silica spheres in an organic solvent mixture, in the experimental time window and wavenumber range. From the simulation data for the time derivative of the width function associated with fc(q, t), there is indication of long-time exponential decay of fc(q, t), for wavenumbers around the location of the static structure factor principal peak. The experimental scattering functions in the probed time range are consistent with a time-wavenumber factorization scaling behavior of fc(q, t) that was first reported by Segrè and Pusey [Phys. Rev. Lett. 77, 771 (1996)] for suspensions of hard spheres. Our BD simulation and MCT results predict a significant violation of exact factorization scaling which, however, is approximately restored according to the ASD results when HIs are accounted for, consistent with the experimental findings for fc(q, t). Our study of collective diffusion is amended by simulation and theoretical results for the self-intermediate scattering function, fs(q, t), and its non-Gaussian parameter α2(t) and for the particle mean squared displacement W(t) and its time derivative. Since self-diffusion properties are not assessed in standard DLS measurements, a method to deduce W(t) approximately from fc(q, t) is theoretically validated.

Structural correlations in diffusiophoretic colloidal mixtures with nonreciprocal interactions.

Nonreciprocal effective interaction forces can occur between mesoscopic particles in colloidal suspensions that are driven out of equilibrium. These forces violate Newton's third law actio = reactio on coarse-grained length and time scales. Here we explore the statistical mechanics of Brownian particles with nonreciprocal effective interactions. Our model system is a binary fluid mixture of spherically symmetric, diffusiophoretic mesoscopic particles, and we focus on the time-averaged particle pair- and triplet-correlation functions. Based on the many-body Smoluchowski equation we develop a microscopic statistical theory for the particle correlations and test it by computer simulations. For model systems in two and three spatial dimensions, we show that nonreciprocity induces distinct nonequilibrium pair correlations. Our predictions can be tested in experiments with chemotactic colloidal suspensions.

Classical Liquids in Fractal Dimension.

We introduce fractal liquids by generalizing classical liquids of integer dimensions d=1,2,3 to a noninteger dimension dl. The particles composing the liquid are fractal objects and their configuration space is also fractal, with the same dimension. Realizations of our generic model system include microphase separated binary liquids in porous media, and highly branched liquid droplets confined to a fractal polymer backbone in a gel. Here, we study the thermodynamics and pair correlations of fractal liquids by computer simulation and semianalytical statistical mechanics. Our results are based on a model where fractal hard spheres move on a near-critical percolating lattice cluster. The predictions of the fractal Percus-Yevick liquid integral equation compare well with our simulation results.

Glass transition of charged particles in two-dimensional confinement.

The glass transition of mesoscopic charged particles in two-dimensional confinement is studied by mode-coupling theory. We consider two types of effective interactions between the particles, corresponding to two different models for the distribution of surrounding ions that are integrated out in coarse-grained descriptions. In the first model, a planar monolayer of charged particles is immersed in an unbounded isotropic bath of ions, giving rise to an isotropically screened Debye-Hückel (Yukawa)-type effective interaction. The second, experimentally more relevant system is a monolayer of negatively charged particles that levitate atop a flat horizontal electrode, as frequently encountered in laboratory experiments with complex (dusty) plasmas. A steady plasma current toward the electrode gives rise to an anisotropic effective interaction potential between the particles, with an algebraically long-ranged in-plane decay. In a comprehensive parameter scan that covers the typical range of experimentally accessible plasma conditions, we calculate and compare the mode-coupling predictions for the glass transition in both kinds of systems.

Rotational self-diffusion in suspensions of charged particles: simulations and revised Beenakker-Mazur and pairwise additivity methods.

We present a comprehensive joint theory-simulation study of rotational self-diffusion in suspensions of charged particles whose interactions are modeled by the generic hard-sphere plus repulsive Yukawa (HSY) pair potential. Elaborate, high-precision simulation results for the short-time rotational self-diffusion coefficient, D(r), are discussed covering a broad range of fluid-phase state points in the HSY model phase diagram. The salient trends in the behavior of D(r) as a function of reduced potential strength and range, and particle concentration, are systematically explored and physically explained. The simulation results are further used to assess the performance of two semi-analytic theoretical methods for calculating D(r). The first theoretical method is a revised version of the classical Beenakker-Mazur method (BM) adapted to rotational diffusion which includes a highly improved treatment of the salient many-particle hydrodynamic interactions. The second method is an easy-to-implement pairwise additivity (PA) method in which the hydrodynamic interactions are treated on a full two-body level with lubrication corrections included. The static pair correlation functions required as the only input to both theoretical methods are calculated using the accurate Rogers-Young integral equation scheme. While the revised BM method reproduces the general trends of the simulation results, it significantly underestimates D(r). In contrast, the PA method agrees well with the simulation results for D(r) even for intermediately concentrated systems. A simple improvement of the PA method is presented which is applicable for large concentrations.

Short-time diffusion in concentrated bidisperse hard-sphere suspensions.

Diffusion in bidisperse Brownian hard-sphere suspensions is studied by Stokesian Dynamics (SD) computer simulations and a semi-analytical theoretical scheme for colloidal short-time dynamics, based on Beenakker and Mazur's method [Physica A 120, 388-410 (1983); 126, 349-370 (1984)]. Two species of hard spheres are suspended in an overdamped viscous solvent that mediates the salient hydrodynamic interactions among all particles. In a comprehensive parameter scan that covers various packing fractions and suspension compositions, we employ numerically accurate SD simulations to compute the initial diffusive relaxation of density modulations at the Brownian time scale, quantified by the partial hydrodynamic functions. A revised version of Beenakker and Mazur's δγ-scheme for monodisperse suspensions is found to exhibit surprisingly good accuracy, when simple rescaling laws are invoked in its application to mixtures. The so-modified δγ scheme predicts hydrodynamic functions in very good agreement with our SD simulation results, for all densities from the very dilute limit up to packing fractions as high as 40%.

Highly asymmetric electrolytes in the primitive model: hypernetted chain solution in arbitrary spatial dimensions.

The pair-correlation functions for fluid ionic mixtures in arbitrary spatial dimensions are computed in hypernetted chain (HNC) approximation. In the primitive model (PM), all ions are approximated as nonoverlapping hyperspheres with Coulomb interactions. Our spectral HNC solver is based on a Fourier-Bessel transform introduced by Talman (J. Comput. Phys. 1978, 29, 35), with logarithmically spaced computational grids. Numeric efficiency for arbitrary spatial dimensions is a commonly exploited virtue of this transform method. Here, we highlight another advantage of logarithmic grids, consisting in efficient sampling of pair-correlation functions for highly asymmetric ionic mixtures. For three-dimensional fluids, ion size and charge-ratios larger than 1000 can be treated, corresponding to hitherto computationally not accessed micrometer-sized colloidal spheres in 1-1 electrolyte. Effective colloidal charge numbers are extracted from our PM results. For moderately large ion size and charge-asymmetries, we present molecular dynamics simulation results that agree well with the approximate HNC pair correlations.

Sterically stabilized colloids with tunable repulsions.

When studying tunable electrostatic repulsions in aqueous suspensions of charged colloids, irreversible colloid aggregation or gelation may occur at high salt concentrations. For many commonly used synthetic colloids, such as polystyrene and silica particles, the reason for coagulation is the presence of unbalanced, strongly attractive, and short-ranged van der Waals (VDW) forces. Here, we present an aqueous polystyrene model colloid that is sterically stabilized against VDW attractions. We show that the synthesis procedure, based on a neutral initiator couple and a nonionic surfactant, introduces surface charges that can be further increased by the addition of charged comonomer methacrylic acid. Thus, the interactions between the polystyrene spheres can be conveniently tuned from hard-sphere-like to charge-stabilized with long-ranged electrostatic repulsions described by a Yukawa-type pair potential. The particle size, grafting density, core-shell structure, and surface charge are characterized by light and neutron scattering. Using X-ray and neutron scattering in combination with an accurate analytic integral equation scheme for the colloidal static structure factor, we deduce effective particle charges for colloid volume fractions ≥0.1 and salt concentrations in the range of 1.5 to 50 mM.

Structure and short-time dynamics in concentrated suspensions of charged colloids.

We report a comprehensive joint experimental-theoretical study of the equilibrium pair-structure and short-time diffusion in aqueous suspensions of highly charged poly-acrylate (PA) spheres in the colloidal fluid phase. Low-polydispersity PA sphere systems with two different hard-core radii, R(0) = 542 and 1117 Å, are explored over a wide range of concentrations and salinities using static and dynamic light scattering (DLS), small angle x-ray scattering, and x-ray photon correlation spectroscopy (XPCS). The measured static and dynamic scattering functions are analyzed using state-of-the-art theoretical methods. For all samples, the measured static structure factor, S(Q), is in good agreement with results by an analytical integral equation method for particles interacting by a repulsive screened Coulomb plus hard-core pair potential. In our DLS and XPCS measurements, we have determined the short-time diffusion function D(Q) = D(0) H(Q)∕S(Q), comprising the free diffusion coefficient D(0) and the hydrodynamic function H(Q). The latter is calculated analytically using a self-part corrected version of the δγ-scheme by Beenakker and Mazur which accounts approximately for many-body hydrodynamic interactions (HIs). Except for low-salinity systems at the highest investigated volume fraction φ ≈ 0.32, the theoretical predictions for H(Q) are in excellent agreement with the experimental data. In particular, the increase in the collective diffusion coefficient D(c) = D(Q → 0), and the decrease of the self-diffusion coefficient, D(s) = D(Q → ∞), with increasing φ is well described. In accord with the theoretical prediction, the peak value, H(Q(m)), of H(Q) relates to the nearest neighbor cage size ∼2π∕Q(m), for which concentration scaling relations are discussed. The peak values H(Q(m)) are globally bound from below by the corresponding neutral hard-spheres peak values, and from above by the limiting peak values for low-salinity charge-stabilized systems. HIs usually slow short-time diffusion on colloidal length scales, except for the cage diffusion coefficient, D(cge) = D(Q(m)), in dilute low-salinity systems where a speed up of the system dynamics and corresponding peak values of H(Q(m)) > 1 are observed experimentally and theoretically.

Short-time rheology and diffusion in suspensions of Yukawa-type colloidal particles.

A comprehensive study is presented on the short-time dynamics in suspensions of charged colloidal spheres. The explored parameter space covers the major part of the fluid-state regime, with colloid concentrations extending up to the freezing transition. The particles are assumed to interact directly by a hard-core plus screened Coulomb potential, and indirectly by solvent-mediated hydrodynamic interactions. By comparison with accurate accelerated Stokesian Dynamics (ASD) simulations of the hydrodynamic function H(q), and the high-frequency viscosity η(∞), we investigate the accuracy of two fast and easy-to-implement analytical schemes. The first scheme, referred to as the pairwise additive (PA) scheme, uses exact two-body hydrodynamic mobility tensors. It is in good agreement with the ASD simulations of H(q) and η(∞), for smaller volume fractions up to about 10% and 20%, respectively. The second scheme is a hybrid method combining the virtues of the δγ scheme by Beenakker and Mazur with those of the PA scheme. It leads to predictions in good agreement with the simulation data, for all considered concentrations, combining thus precision with computational efficiency. The hybrid method is used to test the accuracy of a generalized Stokes-Einstein (GSE) relation proposed by Kholodenko and Douglas, showing its severe violation in low salinity systems. For hard spheres, however, this GSE relation applies decently well.

Pair structure of the hard-sphere Yukawa fluid: an improved analytic method versus simulations, Rogers-Young scheme, and experiment.

We present a comprehensive study of the equilibrium pair structure in fluids of nonoverlapping spheres interacting by a repulsive Yukawa-like pair potential, with special focus on suspensions of charged colloidal particles. The accuracy of several integral equation schemes for the static structure factor, S(q), and radial distribution function, g(r), is investigated in comparison to computer simulation results and static light scattering data on charge-stabilized silica spheres. In particular, we show that an improved version of the so-called penetrating-background corrected rescaled mean spherical approximation (PB-RMSA) by Snook and Hayter [Langmuir 8, 2880 (1992)], referred to as the modified PB-RMSA (MPB-RMSA), gives pair structure functions which are in general in very good agreement with Monte Carlo simulations and results from the accurate but nonanalytical and therefore computationally more expensive Rogers-Young integral equation scheme. The MPB-RMSA preserves the analytic simplicity of the standard rescaled mean spherical (RMSA) solution. The combination of high accuracy and fast evaluation makes the MPB-RMSA ideally suited for extensive parameter scans and experimental data evaluation, and for providing the static input to dynamic theories. We discuss the results of extensive parameter scans probing the concentration scaling of the pair structure of strongly correlated Yukawa particles, and we determine the liquid-solid coexistence line using the Hansen-Verlet freezing rule.