Role of polydispersity in anomalous interactions in electrostatically levitated colloidal systems

Anisotropic colloidal interactions & assembly in AC electric fields

We match experimental and simulated configurations of anisotropic epoxy colloidal particles in high frequency AC electric fields by identifying analytical potentials for dipole-field and dipole-dipole interactions. We report an inverse Monte Carlo simulation algorithm to determine optimal fits of analytical potentials by matching simulated and experimental distribution functions for non-uniform liquid, liquid crystal, and crystal microstructures in varying amplitude electric fields. Two potentials that include accurate particle volume and dimensions along with a concentration dependent prefactor quantitatively capture experimental observations. At low concentrations, an effective ellipsoidal point dipole potential works well, whereas a novel stretched point dipole potential is found to be suitable at all concentrations, field amplitudes, and degrees of ordering. The simplicity, accuracy, and adjustability of the stretched point dipole potential suggest it can be applied to model field mediated microstructures and assembly of systematically varying anisotropic particle shapes.

Critical Casimir interactions of colloids in micellar critical solutions

We study the temperature-dependence of critical Casimir interactions in a critical micellar solution of the nonionic surfactant C12E5 dissolved in water. Experimentally, this is achieved with total internal reflection microscopy (TIRM) which measures the interaction between a single particle and a flat wall. For comparison, we also studied the pair interactions of a two dimensional layer of colloidal particles in the identical micellar system which yields good agreement with the TIRM results. Although, at the surfactant concentration considered here, the fluid forms a dynamical network of wormlike micelles whose structure is considerably more complex than that of simple critical molecular fluids, the temperature-dependence of the measured interactions is - for surface-to-surface distances above 160 nm - in excellent quantitative agreement with theory. Below 160 nm, deviations arise which we attribute to the adsorption of micelles to the interacting surfaces.

Effective structure of a system with continuous polydispersity

In a system of N particles, with continuous size polydispersity, there exists an N(N - 1) number of partial structure factors, making it analytically less tractable. A common practice is to treat the system as an effective one component system, which is known to exhibit an artificial softening of the structure. The aim of this study is to describe the system in terms of M pseudospecies such that we can avoid this artificial softening but, at the same time, have a value of M ≪ N. We use potential energy and pair excess entropy to estimate an optimum number of species, M₀. We then define the maximum width of polydispersity, Δσ₀, that can be treated as a monodisperse system. We show that M₀ depends on the degree and type of polydispersity and also on the nature of the interaction potential, whereas Δσ₀ weakly depends on the type of polydispersity but shows a stronger dependence on the type of interaction potential. Systems with a softer interaction potential have a higher tolerance with respect to polydispersity. Interestingly, M₀ is independent of system size, making this study more relevant for bigger systems. Our study reveals that even 1% polydispersity cannot be treated as an effective monodisperse system. Thus, while studying the role of polydispersity by using the structure of an effective one component system, care must be taken in decoupling the role of polydispersity from that of the artificial softening of the structure.

Domain size polydispersity effects on the structural and dynamical properties in lipid monolayers with phase coexistence

In lipid monolayers with phase coexistence, domains of the liquid-condensed phase always present size polydispersity. However, very few theoretical works consider size distribution effects on the monolayer properties. Because of the difference in surface densities, domains have excess dipolar density with respect to the surrounding liquid expanded phase, originating a dipolar inter-domain interaction. This interaction depends on the domain area, and hence the presence of a domain size distribution is associated with interaction polydispersity. Inter-domain interactions are fundamental to understanding the structure and dynamics of the monolayer. For this reason, it is expected that polydispersity significantly alters monolayer properties. By means of Brownian dynamics simulations, we study the radial distribution function (RDF), the average mean square displacement and the average time-dependent self-diffusion coefficient, D(t), of lipid monolayers with normally distributed size domains. For this purpose, we vary the relevant system parameters, polydispersity and interaction strength, within a range of experimental interest. We also analyze the consequences of using a monodisperse model to determine the interaction strength from an experimental RDF. We find that polydispersity strongly affects the value of the interaction strength, which is greatly underestimated if polydispersity is not considered. However, within a certain range of parameters, the RDF obtained from a polydisperse model can be well approximated by that of a monodisperse model, by suitably fitting the interaction strength, even for 40% polydispersities. For small interaction strengths or small polydispersities, the polydisperse systems obtained from fitting the experimental RDF have an average mean square displacement and D(t) in good agreement with that of the monodisperse system.

Reconfigurable multi-scale colloidal assembly on excluded volume patterns

The ability to create multi-scale, periodic colloidal assemblies with unique properties is important to emerging applications. Dynamically manipulating colloidal structures via tunable kT-scale attraction can provide the opportunity to create particle-based nano- and microstructured materials that are reconfigurable. Here, we report a novel tactic to obtain reconfigurable, multi-scale, periodic colloidal assemblies by combining thermoresponsive depletant particles and patterned topographical features that, together, reversibly mediate local kT-scale depletion interactions. This method is demonstrated in optical microscopy experiments to produce colloidal microstructures that reconfigure between well-defined ordered structures and disordered fluid states as a function of temperature and pattern feature depth. These results are well described by Monte Carlo simulations using theoretical depletion potentials that include patterned excluded volume. Ultimately, the approach reported here can be extended to control the size, shape, orientation, and microstructure of colloidal assemblies on multiple lengths scales and on arbitrary pre-defined pattern templates.

Colloidal crystal grain boundary formation and motion

The ability to assemble nano- and micro- sized colloidal components into highly ordered configurations is often cited as the basis for developing advanced materials. However, the dynamics of stochastic grain boundary formation and motion have not been quantified, which limits the ability to control and anneal polycrystallinity in colloidal based materials. Here we use optical microscopy, Brownian Dynamic simulations, and a new dynamic analysis to study grain boundary motion in quasi-2D colloidal bicrystals formed within inhomogeneous AC electric fields. We introduce "low-dimensional" models using reaction coordinates for condensation and global order that capture first passage times between critical configurations at each applied voltage. The resulting models reveal that equal sized domains at a maximum misorientation angle show relaxation dominated by friction limited grain boundary diffusion; and in contrast, asymmetrically sized domains with less misorientation display much faster grain boundary migration due to significant thermodynamic driving forces. By quantifying such dynamics vs. compression (voltage), kinetic bottlenecks associated with slow grain boundary relaxation are understood, which can be used to guide the temporal assembly of defect-free single domain colloidal crystals.

Depletion-mediated potentials and phase behavior for micelles, macromolecules, nanoparticles, and hydrogel particles

We report a simple depletion potential that captures measured potentials and phase behavior for micrometer-sized colloids in the presence of unadsorbing charged micelles, charged nanoparticles, nonionic macromolecules, and nonionic hydrogel particles. Total internal reflection microscopy (TIRM) is used to measure net potentials between colloids and surfaces, and video microscopy (VM) is used to measure quasi-2D phase behavior in the same material systems. A modified Asakura-Oosawa (AO) depletion potential is developed to accurately quantify particle-wall potentials and interfacial crystallization via particle-particle potentials in Monte Carlo (MC) simulations. The modified AO potential includes effective depletant sizes, accurate osmotic equations of state, and partition coefficients. Partition coefficients are used as the sole adjustable fitting parameter, although an approach to their theoretical prediction from depletant density profiles is also presented. Our results demonstrate a model that accurately captures depletion interactions and phase behavior in a variety of material systems.

kT-scale colloidal interactions in high-frequency inhomogeneous AC electric fields. II. Concentrated ensembles

We report nonintrusive optical microscopy measurements of ensembles of polystyrene colloids in inhomogeneous AC electric fields as a function of field frequency and particle size. By using an inverse Monte Carlo (MC) simulation analysis of time-averaged particle microstructures, we sensitively measure induced dipole-dipole interactions on the kT energy scale. Measurements are reported for frequencies when the particle polarizability is greater and less than the medium, as well as the crossover between these conditions when dipole-dipole interactions vanish. By using measured single dipole-field interactions and associated parameters from Part I as input in the inverse analysis, the dipole-dipole interactions in this work are accurately modeled with no adjustable parameters for conditions away from the crossover frequency (i.e., |f(CM)| > 0). As dipolar interactions vanish at the crossover, a single frequency-dependent parameter is introduced to account for microstructures that appear to result from weak AC electro-osmotic flow induced interactions. By connecting quantitative measures of equilibrium microstructures and kT-scale dipole-field and dipole-dipole interactions, our findings provide a basis to understand colloidal assembly in inhomogeneous AC electric fields.

Like-charge interactions between colloidal particles are asymmetric with respect to sign

Two-dimensional dispersions of colloidal particles with a range of surface chemistries and electrostatic potentials are characterized under a series of solution ionic strengths. A combination of optical imaging techniques are employed to monitor both the colloid structure and the electrostatic surface potential of individual particles in situ. We find that like-charge multiparticle interactions can be tuned from exclusively repulsive to long-range attractive by changing the particle surface composition. This behavior is strongly asymmetric with respect to the sign of the surface potential. Collective long-range attractive interactions are only observed among negatively charged particles.

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.

Imaging energy landscapes with concentrated diffusing colloidal probes

The ability to locally interrogate interactions between particles and energetically patterned surfaces provides essential information to design, control, and optimize template directed self-assembly processes. Although numerous techniques are capable of characterizing local physicochemical surface properties, no current method resolves interactions between colloids and patterned surfaces on the order of the thermal energy kT, which is the inherent energy scale of equilibrium self-assembly processes. Here, the authors describe video microscopy measurements and an inverse Monte Carlo analysis of diffusing colloidal probes as a means to image three dimensional free energy and potential energy landscapes due to physically patterned surfaces. In addition, they also develop a consistent analysis of self-diffusion in inhomogeneous fluids of concentrated diffusing probes on energy landscapes, which is important to the temporal imaging process and to self-assembly kinetics. Extension of the concepts developed in this work suggests a general strategy to image multidimensional and multiscale physical, chemical, and biological surfaces using a variety of diffusing probes (i.e., molecules, macromolecules, nanoparticles, and colloids).

Interactions from diffraction data: historical and comprehensive overview of simulation assisted methods

A large part of statistical mechanics is concerned with the determination of condensed matter structure on the basis of known microscopic interactions. An increasing emphasis has been put on the opposite situation in the last decades as well, where structural data, e.g. pair-distance statistics, are known from diffraction experiments, and one looks for the corresponding interaction functions. The solution of this inverse problem was searched for within the integral equation theories of condensed matter in the early investigations, but before long computer simulation assisted methods were suggested. The interest in this field showed an increasing trend after some attempts appeared in the late 1980s. Several methods were published in the 1990s, and one-two methods appear annually nowadays. In this paper a comprehensive and historical overview is given on the solution of the inverse problem with simulation assisted methods. Emphasis is put on the theoretical grounds of the methods, on the choice of possible input structural functions, on the numerically local or global schemes of the potential modifications, on some advantages and limits of the different methods and on the scientific impact of the methods.

Evanescent wave excited luminescence from levitated quantum dot modified colloids

Evanescent wave excited luminescence of quantum dot modified polystyrene (QDPS) colloids is investigated to measure potential energy profiles of QDPS colloids electrostatically levitated above a planar glass surface. Luminescence is characterized for three different-sized PS colloids modified with three different-sized QDs using confocal microscopy, emission spectra, flow cytometry, and temporal measurements of levitated and deposited colloids. Colloid-surface potential energy profiles constructed from scattering and luminescence intensity data display excellent agreement with each other, theoretical predictions, and independently measured parameters. QDPS luminescence intensity is indirectly confirmed to have an exponential dependence on height similar to conventional colloidal evanescent wave scattering. Our findings indicate that evanescent wave excited QDPS luminescence could enable total internal reflection microscopy measurements of index-matched hard spheres, multiple specific biomolecular interactions via spectral multiplexing, enhanced morphology-dependent resonance modes, and integrated evanescent wave-video-confocal microscopy experiments not possible with scattering.

Diffusing colloidal probes of protein and synthetic macromolecule interactions

A new approach is described for measuring kT and nanometer scale protein-protein and protein-synthetic macromolecule interactions. The utility of this method is demonstrated by measuring interactions of bovine serum albumin (BSA) and copolymers with exposed polyethyleneoxide (PEO) moieties adsorbed to hydrophobically modified colloids and surfaces. Total internal reflection and video microscopy are used to track three-dimensional trajectories of many single diffusing colloids that are analyzed to yield interaction potentials, mean-square displacements, and colloid-surface association lifetimes. A criterion is developed to identify colloids as being levitated, associated, or deposited based on energetic, spatial, statistical, and temporal information. Whereas levitation and deposition occur for strongly repulsive or attractive potentials, association is exponentially sensitive to weak interactions influenced by adsorbed layer architectures and surface heterogeneity. Systematic experiments reveal how BSA orientation and PEO molecular weight produce adsorbed layers that either conceal or expose substrate heterogeneities to generate a continuum of colloid-surface association lifetimes. These measurements provide simultaneous access to a broad range of information that consistently indicates purely repulsive BSA and PEO interactions and a role for surface heterogeneity in colloid-surface association. The demonstrated capability to measure nonspecific protein interactions provides a basis for future measurements of specific protein interactions.

Mapping patterned potential energy landscapes with diffusing colloidal probes

We report a new method for mapping patterned surfaces based on monitoring the interactions of freely diffusing colloidal probes with pattern features to generate measured potential energy landscapes. Evanescent wave scattering and video microscopy are used to track 3D center positions of nominal 2 microm silica colloids as they diffuse over 5-20-nm-thick patterned gold films. An analysis of ensemble-averaged particle height histograms on different pattern features using Boltzmann's equation produces local electrostatic and van der Waals potentials in good agreement with independent measurements and predictions. Absolute separation is obtained from theoretical fits to measured potential-energy profiles and direct measurement by depositing silica colloids onto gold surfaces via electrophoretic deposition. As colloidal probe and pattern feature dimensions become comparable, potential energy profiles suffer some distortion due to the increased probability of probes sampling pattern feature edges. An analysis of interfacial colloidal probe diffusion in conjunction with potential energy measurements demonstrates a consistent interpretation of dissipative and conservative forces in these measurements. Future extensions of this work should produce similar approaches for interrogating physical, chemical, and biomolecular heterogeneous/patterned surfaces and structures with diffusing colloidal probes.

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.

Like-charge attraction in confinement: myth or truth?

It is general wisdom that like-charged colloidal particles repel each other when suspended in liquids. This is in perfect agreement with mean field theories being developed more than 60 years ago. Accordingly, it was a big surprise when several groups independently reported long-ranged attractive components in the pair potential () of equally charged colloids. This so-called (LCA) was only observed in thin sample cells while the pair-interaction in unconfined suspensions has been experimentally confirmed to be entirely repulsive. Despite considerable experimental and theoretical efforts, LCA remains one of the most challenging mysteries in colloidal science. We experimentally reinvestigate the pair-potential () of charged colloidal particles with digital video microscopy and demonstrate that optical distortions in the particle's images lead to slightly erroneous particle positions. If not properly taken into account, this artefact pretends a minimum in () which was in the past misleadingly interpreted as LCA. After correcting optical distortions we obtain entirely repulsive pair interactions which show good agreement with linearized mean field theories.

Stacking fault structure in shear-induced colloidal crystallization

We report measurements of the spatial distribution of stacking faults in colloidal crystals formed by means of an oscillatory shear field at a particle volume fraction of 52% in a system where the pair potential interactions are mildly repulsive. Stacking faults are directly visualized via confocal laser scanning microscopy. Consistent with previous scattering studies, shear orders the initially amorphous colloids into close-packed planes parallel to the shearing surface. Upon increasing the strain amplitude, the close-packed direction of the (111) crystal plane shifts from an orientation parallel to the vorticity direction to parallel the flow direction. The quality of the layer ordering, as characterized by the mean stacking parameter, decreases with strain amplitude. In addition, we directly observe the three-dimensional structure of stacking faults in sheared crystals. We observe and quantify spatial heterogeneity in the stacking fault arrangement in both the flow-vorticity plane and the gradient direction, particularly at high strain amplitudes (gamma> or =3). At these conditions, layer ordering persists in the flow-vorticity plane only over scales of approximately 5-10 particle diameters. This heterogeneity is one component of the random layer ordering deduced from previous scattering studies. In addition, in the gradient direction, the stacking registry shows that crystals with intermediate global mean stacking probability are comprised of short sequences of face-centered cubic and hexagonal close-packed layers with a stacking that includes a component that is nonrandom and alternating in character.

Anomalous potentials from inverse analyses of interfacial polydisperse attractive colloidal fluids

This paper investigates effects of using monodisperse inverse analyses to extract particle-particle and particle-surface potentials from simulated interfacial colloidal fluids of polydisperse attractive particles. Effects of polydispersity are investigated as functions of particle concentration and attractive well depth and range for van der Waals and depletion potentials. Forward Monte Carlo simulations are used to generate particle distribution functions for polydisperse interfacial colloidal fluids from which inverted potentials are obtained using an inverse Ornstein-Zernike analysis and an inverse Monte Carlo simulation method. Attractive potentials are successfully recovered for monodisperse colloidal fluids, but polydispersity that is unaccounted for in inverse analyses produces (1) apparent softening of strong forces, (2) anomalous repulsive and attractive interactions, and (3) aphysical particle overlaps. This investigation provides insights into the role of polydispersity in altering the equilibrium structure and corresponding inverted potentials of attractive colloidal fluids near surfaces. These findings should assist the design and interpretation of optical microscopy experiments involving interfacial colloidal fluids similar to the simulated experiments reported here.