Role of polydispersity in anomalous interactions in electrostatically levitated colloidal systems

In this paper, we investigate the effects of using inverse analyses developed for monodisperse particles to extract particle-particle and particle-surface potentials from simulated interfacial colloidal configurations having finite-size polydispersity. Forward Monte Carlo simulations are used to generate three-dimensional equilibrium configurations of log normal-distributed polydisperse particles confined by gravity near an underlying surface. Particles remain levitated above the substrate and stabilized against aggregation by repulsive electrostatic Derjaguin-Landau-Verwey-Overbeek pair potentials. An inverse Ornstein-Zernike analysis and an inverse Monte Carlo simulation method are used to obtain interactions from simulated distribution functions as a function of polydispersity (sigma), relative range of repulsion (kappa a), and projected interfacial concentration (rho). Both inverse analyses successfully recover input potentials for all monodisperse cases, but fail for polydispersities often encountered in experiments. For different conditions (sigma, kappa a, and rho), our results indicate softened short-range repulsion, anomalous long-range attraction, and apparent particle overlaps, which are similar to commonly reported observations in optical microscopy measurements of quasi-two-dimensional interfacial colloidal ensembles. By demonstrating signatures of, and limitations due to, polydispersity when extracting pair potentials from measured distribution functions, our specific goal is to provide a basis to objectively interpret and resolve the effects of polydispersity in optical microscopy experiments.     

Self-diffusion in submonolayer colloidal fluids near a wall

Theoretical expressions are developed to describe self-diffusion in submonolayer colloidal fluids that require only equilibrium structural information as input. Submonolayer colloidal fluids are defined for the purpose of this work to occur when gravity confines colloids near a planar wall surface so that they behave thermodynamically as two dimensional fluids. Expressions for self-diffusion are generalized to consider different colloid and surface interaction potentials and interfacial concentrations from infinite dilution to near fluid-solid coexistence. The accuracy of these expressions is demonstrated by comparing self-diffusion coefficients predicted from Monte Carlo simulated equilibrium particle configurations with standard measures of self-diffusion evaluated from Stokesian Dynamics simulated particle trajectories. It is shown that diffusivities predicted for simulated equilibrium fluid structures via multibody hydrodynamic resistance tensors and particle distribution functions display excellent agreement with values computed from mean squared displacements and autocorrelation functions of simulated tracer particles. Results are obtained for short and long time self-diffusion both parallel and normal to underlying planar wall surfaces in fluids composed of particles having either repulsive electrostatic or attractive van der Waals interactions. The demonstrated accuracy of these expressions for self-diffusion should allow their direct application to experiments involving submonolayer colloidal fluids having a range of interaction potentials and interfacial concentrations.     

Measurement and interpretation of particle-particle and particle-wall interactions in levitated colloidal ensembles

This paper reports measurements of particle-wall and particle-particle interactions in levitated colloidal ensembles using integrated total internal reflection microscopy (TIRM) and video microscopy (VM) techniques. In levitated colloidal ensembles with area fractions of phiA = 0.03-0.25, ensemble TIRM measured height distribution functions are used to interpret particle-wall interactions, and VM measured pair distribution functions are used to interpret particle-particle interactions using inverse Ornstein-Zernike (OZ) and three-dimensional inverse Monte Carlo (MC) analyses. An inconsistent finding is the observation of an anomalous long-range particle-particle attraction and recovery of the expected Derjaguin-Landau-Verwey-Overbeek (DLVO) particle-wall interactions for all concentrations examined. Because particle-wall and particle-particle potentials are expected to be consistent in several respects, the analytical and experimental methods employed in this investigation are examined for possible sources of error. Comparison of inverse OZ and three-dimensional inverse MC analyses are used to address uncertainties related to dimensionality, effects of particle concentration, and assumptions of the OZ theory and closure relations. The possible influence of charge heterogeneity and particle size polydispersity on measured distribution functions is discussed with regard to inconsistent particle-wall and particle-particle potentials. Ultimately, achieving a consistent understanding of particle-wall and particle-particle interactions in interfacial and confined colloidal systems is essential to numerous complex fluid and advanced material technologies.     

Interfacial colloidal sedimentation equilibrium. I. Intensity based confocal microscopy

This paper reports confocal microscopy measurements of inhomogeneous colloidal sedimentation equilibrium profiles near planar wall surfaces for conditions when colloid dimensions are comparable to the characteristic gravitational length scale. The intensity based confocal method developed in this work enables real-space measurements of one-dimensional density profiles of Brownian colloids without identifying many single colloid centers in large imaging volumes. Measured sedimentation equilibrium profiles for single-phase interfacial fluids and for coexisting inhomogeneous fluid and solid phases are in agreement with a perturbation theory and Monte Carlo simulations within the local density approximation. Monte Carlo simulated colloid scale density profiles display some minor differences with confocal images in terms of microstructural transitions involving the onset of interfacial crystallization and the precise elevation of the fluid-solid interface. These discrepancies are attributed to polydispersity unaccounted for in the analyses, sensitivity of the perturbation theory to the effective hard sphere size, and the influence of ensemble, system size, and box shape in Monte Carlo simulations involving anisotropic/inhomogeneous solids. Successful demonstration of intensity based confocal microscopy provides a basis for future measurements of three-dimensional colloidal interactions, dynamics, and structure near surfaces.     

Interfacial colloidal sedimentation equilibrium. II. Closure-based density functional theory

In Part I [R. E. Beckham and M. A. Bevan, J. Chem. Phys. 127, 164708 (2007)], results were presented for the sedimentation equilibrium of concentrated colloidal dispersions using confocal scanning laser microscopy experiments, Monte Carlo (MC) simulations, and a local density approximation perturbation theory. In this paper, we extended the modeling effort on those systems to include nonlocal density functional theory (DFT), which is capable of predicting the microstructure of the sediment at length scales comparable to the colloidal particle dimension. Specifically, we use a closure-based DFT formulation to predict interfacial colloidal sedimentation equilibrium density profiles. The colloid-colloid and colloid-surface interactions were modeled with DLVO screened electrostatic potentials using parameters taken directly from the experimental work. The DFT profiles were compared to the experimental and MC results from Part I. Good agreement was found for relatively dilute interfacial colloidal fluids, but agreement was less satisfactory as interfacial layering became more pronounced for conditions approaching the onset of interfacial crystallization. We also applied DFT in an inverse sense using the measured colloid density profile to extract the underlying colloid-surface potential; this can be thought of as a microscopic analog to the well-known procedure of using the macroscopic (coarse-grained) density profile to extract the osmotic equation of state. For the dilute interfacial fluid, the inverse DFT calculations reproduced the true colloid-surface potential to within 0.5kT at all elevations.     

Attraction-driven disorder in a hard-core colloidal monolayer

Monte Carlo simulation techniques were employed to explore the effect of short-range attraction on the orientational ordering in a two-dimensional assembly of monodisperse spherical particles. We find that if the range of square-well attraction is approximately 15% of the particle diameter, the dense attractive fluid shows the same ordering behavior as the same density fluid composed of purely repulsive hard spheres. Fluids with an attraction range larger than 15% show an enhanced tendency to crystallization, while disorder occurs for fluids with an attractive range shorter than 15% of the particle diameter. A possible link with the existence of "repulsive" and "attractive" states in dense colloidal systems is discussed.     

Effect of quenched size polydispersity on the fluid-solid transition in charged colloidal suspensions

We study the effect of quenched size polydispersity on the phase behavior of charged colloidal suspensions using free-energy calculations in Monte Carlo simulations. The colloids are assumed to interact with a hard-core repulsive Yukawa (screened-Coulomb) interaction with constant surface potential, so that the particles are polydisperse both in size and charge. In addition, we take the size distribution to be fixed in both the fluid and crystal phase (no size fractionation is allowed). We study the fluid-solid transition for various screening lengths and surface potentials, finding that upon increasing the size polydispersity the freezing transition shifts toward higher packing fractions and the density discontinuity between the two coexisting phases diminishes. Our results provide support for a terminal polydispersity above which the freezing transition disappears.     

Structure of polydisperse inverse ferrofluids: theory and computer simulation

By using theoretical analysis and molecular dynamics simulations, we investigate the structure of colloidal crystals formed by nonmagnetic microparticles (or magnetic holes) suspended in ferrofluids (called inverse ferrofluids), by taking into account the effect of polydispersity in size of the nonmagnetic microparticles. Such polydispersity often exists in real situations. We obtain an analytical expression for the interaction energy of monodisperse, bidisperse, and polydisperse inverse ferrofluids. Body-centered tetragonal (bct) lattices are shown to possess the lowest energy when compared with other sorts of lattices and thus serve as the ground state of the systems. Also, the effect of microparticle size distributions (namely, polydispersity in size) plays an important role in the formation of various kinds of structural configurations. Thus, it seems possible to fabricate colloidal crystals by choosing appropriate polydispersity in size.     

How do interactions control droplet size during nanoprecipitation?

Nanoprecipitation provides colloidal dispersions through successive recombination events between nanometric objects. In the present article, we explain why the nanoprecipitation pathways induced through solvent-shifting – the Ouzo effect –, are fascinating study-cases. Indeed, they allow to address the question of how the interactions between the colloidal particles control the dynamics of the process, thus the particle size distribution. Experimental monitoring of the precipitation dynamics demonstrates that the colloidal dispersion polydispersity decreases over time as the droplets coalesce. Monte Carlo simulations within the Smoluchowski framework agree quantitatively with these observations, and show how the interactions between the particles naturally force the system to become nearly monodisperse. The mechanistic understanding gained from the solvent-shifting experiments is also relevant to other nanoprecipitation processes.     

Investigation of particle interactions in concentrated colloidal suspensions using frequency domain photon migration: monodisperse systems

Frequency domain photon migration (FDPM) measurements were conducted to assess particle interactions of concentrated, monodisperse, polystyrene samples obtained directly from industry by using multiple scattering light. The angle-integrated static structure factor, S(q), and static structure factor at small wave vector q, S(0), were obtained from FDPM measurements at high volume fractions ranging from 0.05 to 0.3, and were compared with those obtained from the monodisperse hard sphere Percus-Yevick (HSPY) model. Effects of different colloid sizes on structure factor evaluated at two different wavelengths were also investigated. Results show that the monodisperse HSPY model is suitable for accounting for particle interactions and local microstructures in these colloidal suspensions. Upon using the HSPY model, particle sizes of polystyrene suspensions were recovered from FDPM measurements at high volume fractions (up to 0.3), which agree well with the DLS measurement of diluted sample ( approximately 0.001). The study of polydispersity effect shows that the FDPM method can be successfully used for recovering the mean particle size of polydisperse colloidal suspension with low polydispersity (<16%) under the assumption of monodisperse hard sphere systems.     

Liquid-gas coexistence and critical point shifts in size-disperse fluids

Specialized Monte Carlo simulations and the moment free energy (MFE) method are employed to study liquid-gas phase equilibria in size-disperse fluids. The investigation is made subject to the constraint of fixed polydispersity, i.e., the form of the "parent" density distribution rho(0)(sigma) of the particle diameters sigma, is prescribed. This is the experimentally realistic scenario for, e.g., colloidal dispersions. The simulations are used to obtain the cloud and shadow curve properties of a Lennard-Jones fluid having diameters distributed according to a Schulz form with a large (delta approximately 40%) degree of polydispersity. Good qualitative accord is found with the results from a MFE method study of a corresponding van der Waals model that incorporates size dispersity both in the hard core reference and the attractive parts of the free energy. The results show that polydispersity engenders considerable broadening of the coexistence region between the cloud curves. The principal effect of fractionation in this region is a common overall scaling of the particle sizes and typical interparticle distances, and we discuss why this effect is rather specific to systems with Schulz diameter distributions. Next, by studying a family of such systems with distributions of various widths, we estimate the dependence of the critical point parameters on delta. In contrast to a previous theoretical prediction, size dispersity is found to raise the critical temperature above its monodisperse value. Unusually for a polydisperse system, the critical point is found to lie at or very close to the extremum of the coexistence region in all cases. We outline an argument showing that such behavior will occur whenever polydispersity affects only the range, rather than the strength of the interparticle interactions.     

Two-dimensional Monte Carlo simulations of a polydisperse colloidal dispersion composed of ferromagnetic particles for the case of no external magnetic field

We have investigated the aggregation phenomena in a polydisperse colloidal dispersion composed of ferromagnetic particles by means of the cluster-moving Monte Carlo method. The results have been compared with those for a monodisperse system. The internal structures of aggregates have been analyzed in terms of the radial distribution function in order to clarify the quantitative differences in the internal structures of clusters. In addition, the cluster size distribution and angular distribution function have been investigated. The results obtained in the present study are summarized as follows. In a monodisperse system, open necklacelike clusters are formed and they extend with increasing strength of the magnetic particle-particle interaction. In a polydisperse system with a small standard deviation in the particle size distribution, sigma=0.2, larger necklacelike clusters are formed and some looplike clusters can also be observed. In a polydisperse system with a larger standard deviation, sigma=0.35, clumplike clusters are formed for a weak magnetic particle-particle interaction. For a stronger magnetic interaction, larger size clusters that exhibit a complicated network structure are formed. These complicated cluster formations found in a polydisperse system are mainly due to the effect of the presence of larger particles.     

Particle Penetration into Polydisperse Polymer Brushes: A Theoretical Analysis

Based on theoretical analysis, the effect of polydispersity on particle penetration into polydisperse polymer brushes is investigated. Three different polydispersities representing sharp, moderate, and extremely wide chain length distributions are chosen, since the corresponding explicit expressions of brush density at these polydispersities are available. To simplify the discussion, this study is restricted to spherical particles of small size which ensure that the particle insertion only causes local conformational perturbations. By analyzing the particle distribution, it is found that polydispersity always facilitates particle penetration. This prediction is confirmed by analyzing the surface fluctuations of the brushes. Interestingly, uniform scaling relations are observed for particles penetrating into monodisperse and moderately polydisperse brushes. The uniformity predicted by monodisperse and moderately polydisperse brushes originates from the same asymptotic behavior of their densities approaching the brush edge. This indicates that polydispersity brings significant influence only at high polydispersities.     

Structure formation and fractionation in systems of polydisperse attractive rods

Fractionation effects and the formation of structured domains are investigated in polydisperse systems of attractive spherocylinders with the help of Monte Carlo simulations. For sufficiently high attractive interaction strength and pressure, the large rods in the system aggregate and form a highly ordered hexatic monolayer that coexists with an isotropic fluid of smaller rods. Fractionation diminishes with decreasing interaction strength but is still observed for hard rod systems, in which the large rods form a nematic droplet rather than a monolayer. Results for polydisperse systems are accompanied by phase diagrams for monodisperse systems of attractive spherocylinders. Here, the phase behavior is shown as a function of rod length and pressure.     

Size fractionation in a phase-separated colloidal fluid

Phase separation of a polydisperse colloidal dispersion implies size fractionation. An application of this effect is given by size-selective purification procedures associated with the colloidal synthesis of so-called monodisperse nanoparticles. We used electron microscopy to determine detailed particle size distributions of coexisting colloidal fluid phases containing highly polydisperse iron oxide nanoparticles with a log-normal distribution (sigma = 0.54 for the total system). Analysis of N approximately 10000 particles per phase yields the first five statistical moments of the distributions. Within experimental error, the interdependence of the statistical moments is in quantitative agreement with the "universal law of fractionation" proposed by Evans, Fairhurst, and Poon [Phys. Rev. Lett. 1998, 81, 1326], even though the theory was derived in the limit of slight polydispersity.     

Polymer Polydispersity Effect on Depletion Interaction between Colloidal Particles

The effect of polymer polydispersity on the polymer‐induced interaction between colloidal particles due to non‐adsorbing ideal chains is investigated. An analytical theory is developed for the polymer‐segment density between two plates and in the space surrounding two spheres by extending a recently proposed superposition approximation to include polymer polydispersity. Monte Carlo computer simulations were made to test the validity of the analytical theory. The polymer densities predicted by the superposition approximation are in reasonable agreement with simulation results for the polydisperse case. The simulations show that depletion leads to a size fractionation of the polymers. It is shown that size polydispersity has a small effect on the interaction between two parallel plates but a more significant effect on the interaction between two spheres. The range of the potential increases and the contact potential drops with increasing polydispersity.

Polymer‐segment density as a function of y for three values of x, as indicated, in the space surrounding two colloidal spheres with radius R = R_g0 and h = 0.48R_g0. Symbols are the MC results: polydisperse polymer (○; z = 1) and monodisperse polymer (•) samples. Curves are the predictions of the product‐function approximation for monodisperse polymer (solid lines) and polydisperse polymer (z = 1, dashed lines).     

Structural and dynamical analysis of monodisperse and polydisperse colloidal systems

We present a semigrand ensemble Monte Carlo and Brownian dynamics simulation study of structural and dynamical properties of polydisperse soft spheres interacting via purely repulsive power-law potentials with a varying degree of "softness." Comparisons focus on crystal and amorphous phases at their coexistence points. It is shown through detailed structural analysis that as potential interactions soften, the "quality of crystallinity" of both monodisperse and polydisperse systems deteriorates. In general, polydisperse crystalline phases are characterized by a more ordered structure than the corresponding monodisperse ones (i.e., for the same potential softness). This counter-intuitive feature originates partly from the fact that particles of different sizes may be accommodated more flexibly in a crystal structure and from the reality that coexistence (osmotic) pressure is substantially higher for polydisperse systems. These trends diminish for softer potentials. Potential softness eventually produces substitutionally disordered crystals. However, substitutional order is apparent for the hard-spherelike interactions. Diffusionwise, crystals appear quite robust with a slight difference in the vibrational amplitudes of small and large particles. This difference, again, diminishes with potential softness. Overcrowding in amorphous polydisperse suspensions causes "delayed" diffusion at intermediate times.     

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.     

Investigation of structure factors in dense colloidal suspensions using frequency domain photon migration: polydisperse systems

Frequency domain photon migration (FDPM) technique was employed to investigate the structure factors of dense, polydisperse colloidal suspensions. The angle-integrated structure factors, [S(q)], extracted from FDPM measurements of scattering properties at volume fractions ranging from 0.05 to 0.4, were compared with the values predicted from the polydisperse hard sphere Percus-Yevick (HSPY) model, as well as decoupling approximation (DA) and local monodisperse approximation (LMA) models that incorporated independently measured particle size information. Results show that the polydisperse HSPY model is the most suitable for accounting for particle interactions which predominantly arise from volume exclusion effects. Furthermore, the influence of size polydispersity upon [S(q)] is most significant at high volume fractions. The static structure factors at small wave vector q, S(0), were also assessed from dual wavelength FDPM measurements by using the small wave number approximation as well as the local monodisperse approximation. The measured S(0) agrees well with the values predicted by the polydisperse HSPY model.