Search

The ionic composition and pair correlations in fluid phases of realistically salt-free charged colloidal sphere suspensions are calculated in the primitive model. We obtain the number densities of all ionic species in suspension, including low-molecular weight microions, and colloidal macroions with acidic surface groups, from a self-consistent solution of a coupled physicochemical set of nonlinear algebraic equations and non-mean-field liquid integral equations. Here, we study suspensions of colloidal spheres with sulfonate or silanol surface groups, suspended in demineralized water that is saturated with carbon dioxide under standard atmosphere. The only input required for our theoretical scheme are the acidic dissociation constants pKa, and effective sphere diameters of all involved ions. Our method allows for an ab initio calculation of colloidal bare and effective charges, at high numerical efficiency. [ABSTRACT FROM AUTHOR]

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.

A nonequilibrium method is developed to sample evaporation of a liquid across a planar interface in a stationary scenario by molecular dynamics. The method does not rely on particle insertions which are challenging when they are used to maintain mass conservation. Its algorithm has a low complexity and is well suited for massively parallel simulations that may yield results with an excellent statistical accuracy. Spatially resolved classical profiles, e.g., for temperature, density, and force, are sampled with a high resolution for a varying hydrodynamic velocity of the evaporation flow. Relatively large systems are simulated, allowing for a detailed study of velocity distribution functions. Varying the hydrodynamic velocity from zero to the speed of sound, it is found that the evaporation flux increases asymptotically, reaching about 90% of its maximum value when the hydrodynamic velocity is about half of its maximum value. A deviation from the Maxwell distribution is identified for the transversal particle velocity near the interface which selectively hinders the migration of individual particles from liquid to vapor with its potential well, allowing only the faster ones to escape. The vapor region in the vicinity of the interface exhibits a spread between the transversal and longitudinal temperature, but equipartition is reattained through particle interactions such that Maxwell distributions are found at a certain distance from the interface. A detailed discussion of the atomistic mechanisms during evaporation is provided, facilitating understanding of this ubiquitous process.

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, f c (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 f c (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 f c (q, t), there is indication of long-time exponential decay of f c (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 f c (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 f c (q, t). Our study of collective diffusion is amended by simulation and theoretical results for the self-intermediate scattering function, f s (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 f c (q, t) is theoretically validated.

Molecular dynamics simulations are reported for the evaporation of a liquid into vacuum, where a Lennard-Jones type fluid with truncated and shifted potential at 2.5σ is considered. Vacuum is enforced locally by particle deletion and the liquid is thermostated in its bulk so that heat flows to the planar interface driving stationary evaporation. The length of the non-thermostated transition region between the bulk liquid and the interface Ln is under study. First, it is found for the reduced bulk liquid temperature Tl/Tc = 0.74 (Tc is the critical temperature) that by increasing Ln from 5.2σ to 208σ the interface temperature Ti drops by 17% and the evaporation flux decreases by a factor of 4.4. From a series of simulations for increasing values of Ln, an asymptotic value Ti (∞) of the interface temperature for Ln → ∞ can be estimated which is 21% lower than the bulk liquid temperature Tl. Second, it is found that the evaporation flux is solely determined by the interface temperature Ti, independent on Tl or Ln. Combining these two findings, the evaporation coefficient α of a liquid thermostated on a macroscopic scale is estimated to be α ≈ 0.14 for Tl/Tc = 0.74.

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%.

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.

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., (© 2011 American Institute of Physics)

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.