Suspensions of adhesive colloidal particles in sedimentation equilibrium in a planar pore

Interacting hard-sphere fluids in an external field

We present a method for studying equilibrium properties of interacting fluids in an arbitrary external field. The fluid is composed of monodisperse spherical particles with hard-core repulsion and additional interactions of arbitrary shape and limited range. Our method of analysis is exact in one dimension and provides demonstrably good approximations in higher dimensions. It can cope with homogeneous and inhomogeneous environments. We derive an equation for the pair distribution function. The solution, to be evaluated numerically, in general, or analytically for special cases, enters expressions for the entropy and free energy functionals. For some one-dimensional systems, our approach yields analytic solutions, reproducing available exact results from different approaches.

Interacting hard rods on a lattice: distribution of microstates and density functionals

We derive exact density functionals for systems of hard rods with first-neighbor interactions of arbitrary shape but limited range on a one-dimensional lattice. The size of all rods is the same integer unit of the lattice constant. The derivation, constructed from conditional probabilities in a Markov chain approach, yields the exact joint probability distribution for the positions of the rods as a functional of their density profile. For contact interaction ("sticky core model") between rods, we give a lattice fundamental measure form of the density functional and present explicit results for contact correlators, entropy, free energy, and chemical potential. Our treatment includes inhomogeneous couplings and external potentials.

Closure-based density functional theory applied to interfacial colloidal fluids

The behavior of dense colloidal fluids near surfaces can now be probed in great detail with experimental techniques like confocal microscopy. In fact, we are approaching a point where quantitative comparisons of experiment with particle-level theory, such as classical density functional theory (DFT), are appropriate. In a forward sense, we may use a known surface potential to predict a particle density distribution function from DFT; in an inverse sense, we may use an experimentally measured particle density distribution function to predict the underlying surface potential from DFT. In this paper, we tested the ability of the closure-based DFT of Zhou and Ruckenstein (J. Chem. Phys. 2000, 112, 8079-8082) to perform forward and inverse calculations on potential models commonly employed for colloidal particles and surfaces. To reduce sources of uncertainty in this initial study, Monte Carlo simulation results played the role of experimental data. The combination of Rogers-Young and modified-Verlet closures consistently performed well across the different potential models. For a reasonable range of choices of the density, temperature, and potential parameters, the inversion procedure yielded particle-surface potentials to an accuracy on the order of 0.1kT.

Nonequilibrium sedimentation of colloids on the particle scale

We investigate sedimentation of model hard-sphere-like colloidal dispersions confined in horizontal capillaries using laser scanning confocal microscopy, dynamical density functional theory, and Brownian dynamics computer simulations. For homogenized initial states we obtain quantitative agreement of the results from the respective approaches for the time evolution of the one-body density distribution and the osmotic pressure on the walls. We demonstrate that single-particle information can be obtained experimentally in systems that were initialized further out of equilibrium such that complex lateral patterns form.

Adsorption of a Binary Mixture of Adhesive Fluids in Planar Pores: A Monte Carlo Study

Grand canonical Monte Carlo simulation is used to study the adsorption of a binary sticky hard-sphere fluid mixture in planar pores. The wall-component 1 and wall-component 2 contact densities are determined to calculate the pressure as a function of the composition of the mixture and the separation between the walls. From these data dependence of the solvation force between the plates on pore width is estimated. The simulation results are compared with the predictions of the Percus-Yevick approximation for planar pores. The density profiles of particular components show interesting shapes stemming from the interplay between the steric effects and the competitive adhesion among all possible species pairs. It is shown that narrowing of the pore causes selective partitioning of individual components of the mixture between the bulk phase and the interior of the pore. The agreement between the two methods is better at wider pores and for the component comprised of weakly adhesive particles.

The density profile of hard sphere liquid system under gravity

The density profile of hard sphere liquid under gravity is calculated by using density functional theory and Monte Carlo simulation method. The two methods give consistent results for a wide range of parameters. Meanwhile, the validity range of the density functional theory is also established. The results are quite different from the barometric height distribution rho(z)=rho(0) exp(-zL(G)) in almost all cases studied, which indicates that the interaction between particles plays an important role in the density distribution under external fields. Moreover, the crystallizing phenomenon is also predicted at the bottom part of the system under strong gravitation.

Global and critical test of the perturbation density-functional theory based on extensive simulation of Lennard-Jones fluid near an interface and in confined systems

The structure of a Lennard-Jones (LJ) fluid subjected to diverse external fields maintaining the equilibrium with the bulk LJ fluid is studied on the basis of the third-order+second-order perturbation density-functional approximation (DFA). The chosen density and potential parameters for the bulk fluid correspond to the conditions situated at "dangerous" regions of the phase diagram, i.e., near the critical temperature or close to the gas-liquid coexistence curve. The accuracy of DFA predictions is tested against the results of a grand canonical ensemble Monte Carlo simulation. It is found that the DFA theory presented in this work performs successfully for the nonuniform LJ fluid only on the condition of high accuracy of the required bulk second-order direct correlation function. The present report further indicates that the proposed perturbation DFA is efficient and suitable for both supercritical and subcritical temperatures.

Sedimentation Equilibrium of Colloidal Suspensions in a Planar Pore Based on Density Functional Theory and the Hard-Core Attractive Yukawa Model

The sedimentation equilibrium of colloidal suspensions modeled by hard-core attractive Yukawa (HCAY) fluids in a planar pore is studied. The density profile of the HCAY fluid in a gravitational field and its distribution between the pore and uniform phases are investigated by a density functional theory (DFT) approach, which results from employing a recently proposed parameter-free version of the Lagrangian theorem-based density functional approximation (Zhou, S. Phys. Lett. A 2003, 319, 279) for hard-sphere fluids to the hard-core part of the HCAY fluid, and the second-order functional perturbation expansion approximation to the tail part as was done in a recent partitioned density functional approximation (Zhou, S. Phys. Rev. E: Stat. Phys., Plasmas, Fluids, Relat. Interdiscip. Top. 2003, 68, 061201). The resultant DFT approach is, thus, the first adjustable parameter-free DFT for HCAY fluids. The validity of the present DFT for HCAY fluids of reduced range parameter z(red) = 1.8 under various external potentials is established in the first of the papers cited previously. The present DFT for HCAY fluids can predict the radial distribution function for the bulk HCAY fluid accurately in the colloidal limit (large value of z(red)), and in the hard-sphere limit, its prediction for the density profile of the hard-sphere fluid in a gravitational field is in very good agreement with the existing simulation data. The dependence of the density profile and distribution coefficient on the magnitude of the interparticle attraction, gravitational field, and degree of confinement is investigated in detail by the present DFT approach. Intuitive and qualitative analyses are also compared with the quantitative DFT calculational results.