Direct determination of phase behavior of square-well fluids

We have combined Gibbs ensemble Monte Carlo simulations with the aggregation volume-biased method in conjunction with the Gibbs-Duhem method to provide the first direct estimates for the vapor-solid, vapor-liquid, and liquid-solid phase coexistences of square-well fluids with three different ranges of attraction. Our results are consistent with the previous simulations and verify the notion that the vapor-liquid coexistence behavior becomes metastable for cases where the attraction well becomes smaller than 1.25 times the particle diameter. In these cases no triple point is found.     

Interactions between spherical colloids mediated by a liquid crystal: a molecular simulation and mesoscale study

Monte Carlo simulations and dynamic field theory (DyFT) are used to study the interactions between dilute spherical particles, dispersed in nematic and isotropic phases of a liquid crystal. A recently developed simulation method (expanded ensemble density of states) was used to determine the potential of mean force (PMF) between the two spheres as a function of their separation and size. The PMF was also calculated by a dynamic field theory that describes the evolution of the local tensor order parameter. Both methods reveal an overall attraction between the colloids in the nematic phase; in the isotropic phase, the overall attraction between the colloids is much weaker, whereas the repulsion at short range is stronger. In addition, both methods predict a new topology of the disclination lines, which arises when the particles approach each other. The theory is found to describe the results of simulations remarkably well, down to length scales comparable to the size of the molecules. At separations corresponding to the width of individual molecular layers on the particles' surface, the two methods yield different defect structures. We attribute this difference to the neglect of density inhomogeneities in the DyFT. We also investigate the effects of the size of spherical colloids on their interactions.     

Depletion Interactions Produced by Nonadsorbing Charged and Uncharged Spheroids

The effect of macromolecule shape on the depletion attraction between two hard spherical particles in a solution with nonadsorbing hard spheroidal macromolecules of arbitrary size and aspect ratio was investigated using a modified form of the force-balance model of J. Y. Walz and A. Sharma (1994, J. Colloid Interface Sci. 168, 495). The macromolecules were represented as general spheroids, which could be either charged or uncharged. For the uncharged case, a set of analytical expressions describing the depletion attraction, valid for particles much larger than the characteristic macromolecule size, was developed. Comparisons with the case of spherical macromolecules were made under the condition of either constant macromolecule number density, rho(b), or constant volume fraction, phi. It was found that increasing the spheroidal macromolecule aspect ratio (major axis length/minor axis length) decreases the depletion attraction at constant rho(b), but increases the interaction at constant phi. In the latter case, the interaction produced by prolate macromolecules is greater than that produced by oblate macromolecules of equal axis lengths, while the opposite is true at constant rho(b). A simple scaling analysis is used to explain these trends. Surface charge is found to increase both the range and the magnitude of the depletion attraction; however, the general trends are the same as those found in the uncharged systems. Finally, the effect of the depletion attraction produced by spherical and spheroidal macromolecules on the stability of a dispersion of charged particles was examined. It was found that charged spheroids at concentrations of order 1% volume can produce secondary energy wells of sufficient magnitude to induce flocculation in a dispersion of charged spherical particles. Copyright 2000 Academic Press.     

Haloing, flocculation, and bridging in colloid-nanoparticle suspensions

Integral equation theory with a hybrid closure approximation is employed to study the equilibrium structure of highly size asymmetric mixtures of spherical colloids and nanoparticles. Nonequilibrium contact aggregation and bridging gel formation is also qualitatively discussed. The effect of size asymmetry, nanoparticle volume fraction and charge, and the spatial range, strength, and functional form of colloid-nanoparticle and colloid-colloid attractions in determining the potential-of-mean force (PMF) between the large spheres is systematically explored. For hard, neutral particles with weak colloid-nanoparticle attraction qualitatively distinct forms of the PMF are predicted: (i) a contact depletion attraction, (ii) a repulsive form associated with thermodynamically stable "nanoparticle haloing," and (iii) repulsive at contact but with a strong and tight bridging minimum. As the interfacial cohesion strengthens and becomes shorter range the PMF acquires a deep and tight bridging minimum. At sufficiently high nanoparticle volume fractions, a repulsive barrier then emerges which can provide kinetic stabilization. The charging of nanoparticles can greatly reduce the volume fractions where significant changes of the PMF occur. For direct and interfacial van der Waals attractions, the large qualitative consequences of changing the absolute magnitude of nanoparticle and colloid diameters at fixed size asymmetry ratio are also studied. The theoretical results are compared with recent experimental and simulation studies. Calculations of the real and Fourier space mixture structure at nonzero colloid volume fractions reveal complex spatial reorganization of the nanoparticles due to many body correlations.     

Depletion-induced phase separation in colloid-polymer mixtures

Phase separation can be induced in a colloidal dispersion by adding non-adsorbing polymers. Depletion of polymer around the colloidal particles induces an effective attraction, leading to demixing at sufficient polymer concentration. This communication reviews theoretical and experimental work carried out on the polymer-mediated attraction between spherical colloids and the resulting phase separation of the polymer-colloid mixture. Theoretical studies have mainly focused on the limits where polymers are small or large as compared to the colloidal size. Recently, however, theories are being developed that cover a wider colloid-polymer size ratio range. In practical systems, size polydispersity and polyelectrolytes (instead of neutral polymers) and/or charges on the colloidal surfaces play a role in polymer-colloid mixtures. The limited amount of theoretical work performed on this is also discussed. Finally, an overview is given on experimental investigations with respect to phase behavior and results obtained with techniques enabling measurement of the depletion-induced interaction potential, the structure factor, the depletion layer thickness and the interfacial tension between the demixed phases of a colloid-polymer mixture.     

Many body effects on the phase separation and structure of dense polymer-particle melts

Liquid state theory is employed to study phase transitions and structure of dense mixtures of hard nanoparticles and flexible chains (polymer nanocomposites). Calculations are performed for the first time over the entire compositional range from the polymer melt to the hard sphere fluid. The focus is on polymers that adsorb on nanoparticles. Many body correlation effects are fully accounted for in the determination of the spinodal phase separation instabilities. The nanoparticle volume fraction at demixing is determined as a function of interfacial cohesion strength (or inverse temperature) for several interaction ranges and nanoparticle sizes. Both upper and lower critical temperature demixing transitions are predicted, separated by a miscibility window. The phase diagrams are highly asymmetric, with the entropic depletion-like lower critical temperature occurring at a nanoparticle volume fraction of approximately 10%, and a bridging-induced upper critical temperature at approximately 95% filler loading. The phase boundaries are sensitive to both the spatial range of interfacial cohesion and nanoparticle size. Nonmonotonic variations of the bridging (polymer-particle complex formation) demixing boundary on attraction range are predicted. Moreover, phase separation due to many body bridging effects occurs for systems that are fully stable at a second order virial level. Real and Fourier space pair correlations are examined over the entire volume fraction regime with an emphasis on identifying strong correlation effects. Special attention is paid to the structure near phase separation and the minimum in the potential of mean force as the demixing boundaries are approached. The possibility that nonequilibrium kinetic gelation or nanoparticle cluster formation preempts equilibrium phase separation is discussed.     

Monte Carlo simulation of structure and nanoscale interactions in polymer nanocomposites

Off-lattice Monte Carlo simulations in the canonical ensemble are used to study polymer-particle interactions in nanocomposite materials. Specifically, nanoscale interactions between long polymer chains (N=550) and strongly adsorbing colloidal particles of comparable size to the polymer coils are quantified and their influence on nanocomposite structure and dynamics investigated. In this work, polymer-particle interactions are computed from the integrated force-distance curve on a pair of particles approaching each other in an isotropic polymer medium. Two distinct contributions to the polymer-particle interaction potential are identified: a damped oscillatory component that is due to chain density fluctuations and a steric repulsive component that arises from polymer confinement between the surfaces of approaching particles. Significantly, in systems where particles are in a dense polymer melt, the latter effect is found to be much stronger than the attractive polymer bridging effect. The polymer-particle interaction potential and the van der Waals potential between particles determine the equilibrium particle structure. Under thermodynamic equilibrium, particle aggregation is observed and there exists a fully developed polymer-particle network at a particle volume fraction of 11.3%. Near-surface polymer chain configurations deduced from our simulations are in good agreement with results from previous simulation studies.     

A Method for Directly Counting and Quantitatively Comparing Aggregated Structures during Cluster Formation

Assembling of a few particles into a cluster commonly occurs in many systems. However, it is still challenging to precisely control particle assembling, due to the various amorphous structures induced by thermal fluctuations during cluster formation. Although these structures may have very different degrees of aggregation, a quantitative method is lacking to describe them, and how these structures evolve remains unclear. Therefore a significant step towards precise control of particle self-assembly is to describe and analyze various aggregation structures during cluster formation quantitatively. In this work, we are motivated to propose a method to directly count and quantitatively compare different aggregated structures. We also present several case studies to evaluate how the aggregated structures during cluster formation are affected by external controlling factors, e.g., different interaction ranges, interaction strengths, or anisotropy of attraction.     

Monte Carlo computer simulations and electron microscopy of colloidal cluster formation via emulsion droplet evaporation

We consider a theoretical model for a binary mixture of colloidal particles and spherical emulsion droplets. The hard sphere colloids interact via additional short-ranged attraction and long-ranged repulsion. The droplet-colloid interaction is an attractive well at the droplet surface, which induces the Pickering effect. The droplet-droplet interaction is a hard-core interaction. The droplets shrink in time, which models the evaporation of the dispersed (oil) phase, and we use Monte Carlo simulations for the dynamics. In the experiments, polystyrene particles were assembled using toluene droplets as templates. The arrangement of the particles on the surface of the droplets was analyzed with cryogenic field emission scanning electron microscopy. Before evaporation of the oil, the particle distribution on the droplet surface was found to be disordered in experiments, and the simulations reproduce this effect. After complete evaporation, ordered colloidal clusters are formed that are stable against thermal fluctuations. Both in the simulations and with field emission scanning electron microscopy, we find stable packings that range from doublets, triplets, and tetrahedra to complex polyhedra of colloids. The simulated cluster structures and size distribution agree well with the experimental results. We also simulate hierarchical assembly in a mixture of tetrahedral clusters and droplets, and find supercluster structures with morphologies that are more complex than those of clusters of single particles.     

Colloidal systems with a short-range attraction and long-range repulsion: Phase diagrams,structures, and dynamics

Colloidal systems with both a short-range attraction and long-range repulsion (SALR) have rich phases compared with the traditional hard sphere systems or sticky hard sphere systems. The competition between the short-range attraction and long-range repulsion results in the frustrated phase separation, which leads to the formation of intermediate range order (IRO) structures and introduces new phases to both equilibrium and nonequilibrium phase diagrams, such as clustered fluid, cluster percolated fluid, Wigner glass, and cluster glass. One hallmark feature of many SALR systems is the appearance of the IRO peak in the interparticle structure factor, which is associated with different types of IRO structures. The relationship between the IRO peak and the clustered fluid state has been careful investigated. Not surprisingly, the morphology of clusters in solutions can be affected and controlled by the SALR potential. And the effect of the SALR potential on the dynamic properties is also reviewed here. Even though much progress has been made in understanding SALR systems, many future works are still needed to have quantitative comparisons between experiments and simulations/theories and understand the differences from different experimental systems. Owing to the large parameter space available for SALR systems, many exciting features of SALR systems are not fully explored yet. Because proteins in low-salinity solutions have SALR interactions, the understanding of SALR systems can greatly help understand protein behavior in concentrated solutions or crowded conditions.     

Phase diagram of mixtures of hard colloidal spheres and discs: a free-volume scaled-particle approach

Phase diagrams of mixtures of colloidal hard spheres with hard discs are calculated by means of the free-volume theory. The free-volume fraction available to the discs is determined from scaled-particle theory. The calculations show that depletion induced phase separation should occur at low disc concentrations in systems now experimentally available. The gas-liquid equilibrium of the spheres becomes stable at comparable size ratios as with bimodal mixtures of spheres or mixtures of rods and spheres. Introducing finite thickness of the platelets gives rise to a significant lowering of the fluid branch of the binodal.     

Influence of initial conditions on homogeneous nucleation kinetics in a closed system

The formation of nuclei of a new phase from the supersaturated mother phase in a closed system is studied. The depletion of the mother phase due to phase transition is taken into account. Basic kinetic equations describing such process are solved numerically to determine the number density of nuclei of newly forming phase and nucleation rate. It is shown that in contrary to the standard nucleation model, when the depletion of the mother phase is not taken into account, the initial size distribution of the clusters affects considerably the nucleation process at higher supersaturations. Our model starts with the equilibrium size distribution of clusters up to various cluster sizes in the undercritical region. At lower supersaturation the formation of nuclei is similar to the standard model because of the low depletion of the mother phase. At higher supersaturation, the depletion of the mother phase plays an important role and some extremal value appears at the size distribution of nuclei, which is not observed in the standard model. The extremum in the size distribution is not a consequence of the coalescence process itself, but it is caused rather by the depletion of the mother phase during the phase transformation.     

Effect of polymer size and chain length on depletion interactions between two colloids

A density functional theory based on the weighted density has been developed to investigate the depletion interactions between two colloids immersed in a bath of the binary polymer mixtures, where the colloids are modeled as hard spheres and the polymers as freely jointed tangent hard-sphere chain mixtures. The theoretical calculations for the depletion forces between two colloids induced by the polymer are in good agreement with the computer simulations. The effects of polymer packing fraction, degree of polymerization, polymer/polymer size ratio, colloid/polymer size ratio on the depletion interactions, and colloid-colloid second virial coefficient B2 due to polymer-mediated interactions have been studied. With increasing the polymer packing fraction, the depletion interaction becomes more long ranged and the attractive interaction near the colloid becomes deeper. The effect of degree polymerization shows that the long chain gives a more stable dispersion for colloids rather than the short chain. The strong effective colloid-colloid attraction appears for the large colloid/polymer and polymer/polymer size ratio. The location of maximum repulsion Rmax is found to appear Rmax approximately sigmac+Rg2 for the low polymer packing fraction and this is shifted to smaller separation Rmax approximately sigmac+sigmap2 with increasing the polymer packing fraction, where sigmap2 and Rg2 are the small-particle diameter and the radius of gyration of the polymer with the small-particle diameter, respectively.     

Polymer-induced depletion interaction between weakly attractive plates

Lattice Monte Carlo simulations have been employed to calculate depletion interaction of excluded volume chains in a weakly attractive slit, particularly in the region around the critical point of adsorption. The simulations were performed under full equilibrium conditions where a dilute solution in a slit was in contact with the reservoir. The free energy of confinement deltaA, the force f, and the relative pressurepI/pE on the slit walls were calculated as a function of slit width D and the attraction strength epsilon. The depletion region in the pressure profile pI/pE vs D is reduced by an increase in the attraction potential epsilon in a manner resembling the influence of polymer concentration. At the critical point of adsorption epsilonc the depletion interaction vanishes both in the pressure pI/pE and in the intraslit concentration profile phiI(x). The parameters used to assess the stability of colloidal dispersions such as the depletion potential W(D) (an integral of the net pressure deltap) reach a unique value at the critical condition. A monotonic repulsive profilepI vs D was found for chains trapped in the slit at restricted equilibrium. The mean dimensions (R2) of chains compressed in attractive slits feature a distinct minimum at intermediate slit widths.     

Self-diffusion of reversibly aggregating spheres

Reversible diffusion limited cluster aggregation of hard spheres with rigid bonds was simulated and the self-diffusion coefficient was determined for equilibrated systems. The effect of increasing attraction strength was determined for systems at different volume fractions and different interaction ranges. It was found that the slowing down of the diffusion coefficient due to crowding is decoupled from that due to cluster formation. The diffusion coefficient could be calculated from the cluster size distribution and became zero only at infinite attraction strength when permanent gels are formed. It is concluded that so-called attractive glasses are not formed at finite interaction strength.     

Salt-induced protein phase transitions in drying drops

Protein phase transitions in drying sessile drops of protein-salt-water colloidal systems were studied by means of optical and atom-force microscopy. The following sequence of events was observed during drop drying: attachment of a drop to a glass support; redistribution of colloidal phase due to hydrodynamic centrifugal stream; protein ring formation around the edge; formation of protein spatial structures inside a protein ring that pass into gel in the middle of the drop; salt crystallization in the shrinking gel. It was assumed that rapid drying of a protein ring over the circle of high colloidal volume fraction and low strength of interparticle attraction leads to formation of colloidal glass, whereas gel forms only in the middle of the drop at very low protein volume fraction and strong attraction between the particles. Before gelation, colloidal particles form fractal clusters. In dried drops of salt-free protein solutions, no visual protein structures were observed. Structural evolution of protein in sessile drying drops of protein-salt aqueous colloidal solutions is discussed on the basis of experimental data.     

Interfacial colloidal crystallization via tunable hydrogel depletants

We demonstrate an approach using temperature-dependent hydrogel depletants to thermoreversibly tune colloidal attraction and interfacial colloidal crystallization. Total internal reflection and video microscopy are used to measure temperature-dependent depletion potentials between approximately 2 microm silica colloids and surfaces as mediated by approximately 0.2 microm poly-N-isopropylacrylamide (PNIPAM) hydrogel particles. Measured depletion potentials are modeled using the Asakura-Oosawa theory while treating PNIPAM depletants as swellable hard spheres. Monte Carlo simulations using the measured potentials predict reversible, quasi-2D crystallization and melting at approximately 27 degrees C in quantitative agreement with video microscopy images of measured microstructures (i.e., radial distribution functions) over the temperature range of interest (20-29 degrees C). Additional measurements of short-time self-diffusivities display excellent agreement with predicted diffusivities by considering multibody hydrodynamic interactions and using a swellable hard sphere model for the PNIPAM solution viscosity. Our findings demonstrate the ability to quantitatively measure, model, and manipulate kT-scale depletion attraction and phase behavior as a means of formally engineering interfacial colloidal crystallization.     

Charged micelle depletion attraction and interfacial colloidal phase behavior

Ensemble total internal reflection microscopy (TIRM) is used to directly measure the evolution of colloid-surface depletion attraction with increasing sodium dodecyl sulfate (SDS) concentration near the critical micelle concentration (CMC). Measured potentials are well described by a modified Asakura-Oosawa (AO) depletion potential in addition to electrostatic and van der Waals contributions. The modified AO potential includes effects of electrostatic interactions between micelles and surfaces via effective depletant dimensions in an excluded volume term and partitioning in an osmotic pressure term. Directly measured colloid-surface depletion potentials are used in Monte Carlo (MC) simulations to capture video microscopy (VM) measurements of micelle-mediated quasi-two-dimensional phase behavior including fluid, crystal, and gel microstructures. Our findings provide information to develop more rigorous and analytically simple models of depletion attraction in charged micellar systems.     

Liquid-vapor phase diagram and cluster formation of two-dimensional ionic fluids

Direct molecular dynamics simulations on interfaces at constant temperature are performed to obtain the liquid-vapor phase diagram of the two-dimensional soft primitive model, an equimolar mixture of equal size spheres carrying opposite charges. Constant temperature and pressure simulations are also carried out to check consistency with interface simulations results. In addition, an analysis of the cluster formation of mixtures of particles with charge asymmetry in the range 1:1 to 1:36 at low and high densities is performed. The number of free ions, when plotted as a function of the positive ion charge, Z(+), has an oscillatory behavior and is independent of the density. The formation of aggregates is analyzed in terms of the attraction and repulsion between ions.