Numerical test of the Percus–Yevick approximation for continuum media of adhesive sphere model at percolation threshold

The generalized continuous multiple step (GCMS) potential: model systems and benchmarks

The generalized continuous multiple step (GCMS) potential is presented in this work. Its flexible form allows forrepulsiveand/orattractivecontributions to be encoded through adjustable energy and length scales. The GCMS interaction provides a continuous representation of square-well, square-shoulder potentials and their variants for implementation in computer simulations. A continuous and differentiable energy representation is required to derive forces in conventional simulation algorithms. Molecular dynamics simulations are of particular interest when considering the dynamic properties of a system. The GCMS potential can mimic other interactions with a judicious choice of parameters due to the versatile sigmoid form. In this study, our benchmarks for the GCMS representation include triangular, Yukawa, Franzese, and Lennard-Jones potentials. Comparisons made with published data on volumetric phase diagrams, liquid structure, and diffusivity from model systems are in excellent agreement.

Density functional theory for Baxter's sticky hard spheres in confinement

It has recently been shown that a free energy for Baxter's sticky hard-sphere fluid is uniquely defined within the framework of fundamental measure theory (FMT) for the inhomogeneous hard-sphere fluid, provided that it obeys scaled-particle theory and the Percus-Yevick (PY) result for the direct correlation function [H. Hansen-Goos and J. S. Wettlaufer, J. Chem. Phys. 134, 014506 (2011)]. Here, combining FMT weighted densities with a new vectorial weighted density, we regularize the divergences of the associated strongly confined limit. The free energy that emerges is exact in the zero-dimensional limit, leaves the underlying equation of state unaffected, and yields a direct correlation function distinct from that of PY. Comparison with simulation data for both the bulk pair-correlation function and the density profiles in confinement shows that the new theory is significantly more accurate than the PY-based results.

Percolation transition in fluids: scaling behavior of the spanning probability functions

We examine different spanning probability functions (wrapping and crossing) near the percolation threshold of a supercritical square-well fluid and determine the threshold values of these probabilities, which may be universal for all fluids. It is shown that for a continuous system, over a wide range of system size, the wrapping probabilities can be described by universal scaling functions, whereas the crossing probabilities do not show such universal behavior over the same range of system size. The obtained universal functions for the wrapping probabilities can be used for an estimation of the percolation threshold in fluids in general. The results for the crossing probabilities allow us then to characterize large clusters in real fluids.

Phase separation and percolation of reversibly aggregating spheres with a square-well attraction potential

Reversible aggregation of spheres is simulated using a novel method in which clusters of bound spheres diffuse collectively with a diffusion coefficient proportional to their radius. It is shown that the equilibrium state is the same as with other simulation techniques, but with the present method more realistic kinetics are obtained. The behavior as a function of volume fraction and interaction strength was tested for two different attraction ranges. The binodal and the percolation threshold were determined. The cluster structure and size distribution close to the percolation threshold were found to be consistent with the percolation model. Close to the binodal phase separation occurred through the growth of spherical dense domains, while for deep quenches a system spanning network is formed that coarsens with a rate that decreases with increasing attraction. We found no indication for arrest of the coarsening.

Dynamic signature for the equilibrium percolation threshold of attractive colloidal fluids

Short-time self-diffusivities, D(S)S, are computed for attractive colloidal fluids with van der Waals potentials to identify their percolation threshold and a previously unexplained dynamic transition. Our results show a discontinuous change in the slope of D(S)S vs temperature as the percolation threshold is crossed. Because D(S)S depends only on multibody hydrodynamic interactions, the percolation threshold of attractive colloidal fluids is shown to correspond to a transition in a dynamic property consistent with the linear viscoelastic measurements of Woutersen et al., J. Chem. Phys. 101, 542 (1994).

Specific ion-dependent attraction and phase behavior of polymer-coated colloids

This paper presents results demonstrating the role of temperature and specific ions in mediating attraction between polymer-coated colloids and determining their equilibrium phase behavior. In particular, theoretical predictions of continuum van der Waals attraction between poly(ethylene oxide)-poly(propylene oxide)-poly(ethylene oxide) (PEO-PPO-PEO)-coated polystyrene colloids are used to explain measured temperature and specific ion-dependent fluid-gel transitions in dispersions of these particles. Building on previous studies of PEO-PPO-PEO-coated polystyrene colloids dispersed in aqueous NaCl media, this work reports rheologically measured fluid-gel transitions as a function of temperature and NaCl/MgSO4 composition. Adhesive-sphere predictions of percolation thresholds are fit to measured fluid-gel data by allowing the adsorbed copolymer layer thickness as a single adjustable parameter. This allows the attraction between the PEO-PPO-PEO layers to be interpreted as a function of temperature and NaCl/MgSO4 composition. Quantitative predictions of a polymeric van der Waals attraction associated with the layer collapse in diminishing solvent conditions provides a simple mechanism for explaining the measured fluid-gel data as a dynamic percolation transition. Ultimately, this work identifies the importance of continuum polymeric van der Waals attraction for explaining specific ion-dependent phenomena.

Phase diagram of the adhesive hard sphere fluid

The phase behavior of the Baxter adhesive hard sphere fluid has been determined using specialized Monte Carlo simulations. We give a detailed account of the techniques used and present data for the fluid-fluid coexistence curve as well as parametrized fits for the supercritical equation of state and the percolation threshold. These properties are compared with the existing results of Percus-Yevick theory for this system.

Competition of percolation and phase separation in a fluid of adhesive hard spheres

Using a combination of Monte Carlo techniques, we locate the liquid-vapor critical point of adhesive hard spheres. We find that the critical point lies deep inside the gel region of the phase diagram. The (reduced) critical temperature and density are tau(c)=0.1133+/-0.0005 and rho(c)=0.508+/-0.01. We compare these results with the available theoretical predictions. Using a finite-size scaling analysis, we verify that the critical behavior of the adhesive hard sphere model is consistent with that of the 3D Ising universality class, the default for systems with short-range attractive forces.