Anisotropic colloidal particles in critical fluids

We consider anisotropic colloidal particles with dumbbell or lens shapes that are immersed in a critical binary fluid mixture. The orientation-dependent long-ranged universal interactions mediated by the critical solvent between a particle and a wall or between two particles are investigated for mesoscopic particle sizes small compared to the correlation length and interparticle distances. Exact results are obtained using a "small particle operator expansion." The amplitudes of the isotropic and anisotropic operators in the expansion depend on the size and aspect ratio of the dumbbell or lens and are determined by density profiles in the Ising model at the critical point in a wedge geometry with symmetry-breaking fixed-spin boundary conditions. Dumbbells and ellipsoids with a symmetry preserving surface are also considered.     

Implementation of Lees-Edwards periodic boundary conditions for direct numerical simulations of particle dispersions under shear flow

A general methodology is presented to perform direct numerical simulations of particle dispersions in a shear flow with Lees-Edwards periodic boundary conditions. The Navier-Stokes equation is solved in oblique coordinates to resolve the incompatibility of the fluid motions with the sheared geometry, and the force coupling between colloidal particles and the host fluid is imposed by using a smoothed profile method. The validity of the method is carefully examined by comparing the present numerical results with experimental viscosity data for particle dispersions in a wide range of volume fractions and shear rates including nonlinear shear-thinning regimes.     

The sticking probability of colloidal particles in polymer-induced flocculation

The interaction between colloidal particles in a solution containing polymer molecules is examined theoretically. In particular, the sticking probability of colloidal particles, which plays a significant role in the determination of the collision efficiency of an unstable dispersed system, is analyzed. We found that the sticking probability is extremely sensitive to the variation in the Hamaker constant of a colloidal particle. In the conventional analysis of the collision efficiency in polymer-induced flocculation, it is assumed that the adsorption of polymer is relatively fast than the collision of colloidal particles. Since the adsorption of polymer will significantly affect the properties of the surface of a colloidal particle, this assumption may lead to an appreciable deviation in the prediction of collision efficiency.     

Combining molecular dynamics with Lattice Boltzmann: a hybrid method for the simulation of (charged) colloidal systems

We present a hybrid method for the simulation of colloidal systems that combines molecular dynamics (MD) with the Lattice Boltzmann (LB) scheme. The LB method is used as a model for the solvent in order to take into account the hydrodynamic mass and momentum transport through the solvent. The colloidal particles are propagated via MD and they are coupled to the LB fluid by viscous forces. With respect to the LB fluid, the colloids are represented by uniformly distributed points on a sphere. Each such point [with a velocity V(r) at any off-lattice position r] is interacting with the neighboring eight LB nodes by a frictional force F = xi0(V(r)-u(r)), with xi0 being a friction coefficient and u(r) being the velocity of the fluid at the position r. Thermal fluctuations are introduced in the framework of fluctuating hydrodynamics. This coupling scheme has been proposed recently for polymer systems by Ahlrichs and Dunweg [J. Chem. Phys. 111, 8225 (1999)]. We investigate several properties of a single colloidal particle in a LB fluid, namely, the effective Stokes friction and long-time tails in the autocorrelation functions for the translational and rotational velocity. Moreover, a charged colloidal system is considered consisting of a macroion, counterions, and coions that are coupled to a LB fluid. We study the behavior of the ions in a constant electric field. In particular, an estimate of the effective charge of the macroion is yielded from the number of counterions that move with the macroion in the direction of the electric field.     

Particle Interactions in Diffusiophoresis and Electrophoresis of Colloidal Spheres with Thin but Polarized Double Layers

The diffusiophoretic and electrophoretic motions of two colloidal spheres in the solution of a symmetrically charged electrolyte are analyzed using a method of reflections. The particles are oriented arbitrarily with respect to the electrolyte gradient or the electric field, and they are allowed to differ in radius and in zeta potential. The thickness of the electric double layers surrounding the particles is assumed to be small relative to the radius of each particle and to the gap width between the particles, but the effect of polarization of the mobile ions in the diffuse layer is taken into account. A slip velocity of fluid and normal fluxes of solute ions at the outer edge of the thin double layer are used as the boundary conditions for the fluid phase outside the double layers. The method of reflections is based on an analysis of the electrochemical potential and fluid velocity disturbances produced by a single dielectric sphere placed in an arbitrarily varying electrolyte gradient or electric field. The solution for two-sphere interactions is obtained in expansion form correct to O(r(12)(-7)), where r(12) is the distance between the particle centers. Our analytical results are found to be in good agreement with the available numerical solutions obtained using a boundary collocation method. On the basis of a model of statistical mechanics, the results of two-sphere interactions are used to analytically determine the first-order effect of the volume fraction of particles of each type on the mean diffusiophoretic and eletrophoretic velocities in a bounded suspension. For a suspension of identical spheres, the mean diffusiophoretic velocity can be decreased or increased as the volume fraction of the particles is increased, while the mean electrophoretic velocity is reduced with the increase in the particle concentration. Generally speaking, the particle interaction effects can be quite significant in typical situations. Copyright 2000 Academic Press.     

Size-sorted anionic iron oxide nanomagnets as colloidal mediators for magnetic hyperthermia

Iron oxide colloidal nanomagnets generate heat when subjected to an alternating magnetic field. Their heating power, governed by the mechanisms of magnetic energy dissipation for single-domain particles (Brown and Néel relaxations), is highly sensitive to the crystal size, the material, and the solvent properties. This study was designed to distinguish between the contributions of Néel and Brownian mechanisms to heat generation. Anionic nanocrystals of maghemite and cobalt ferrite, differing by their magnetic anisotropy, were chemically synthesized and dispersed in an aqueous suspension by electrostatic stabilization. The particles were size-sorted by successive electrostatic phase separation steps. Parameters governing the efficiency of nanomagnets as heat mediators were varied independently; these comprised the particle size (from 5 to 16.5 nm), the solvent viscosity, magnetic anisotropy, and the magnetic field frequency and amplitude. The measured specific loss powers (SLPs) were in quantitative agreement with the results of a predictive model taking into account both Néel and Brown loss processes and the whole particle size distribution. By varying the carrier fluid viscosity, we found that Brownian friction within the carrier fluid was the main contributor to the heating power of cobalt ferrite particles. In contrast, Néel internal rotation of the magnetic moment accounted for most of the loss power of maghemite particles. Specific loss powers were varied by 3 orders of magnitude with increasing maghemite crystal size (from 4 to 1650 W/g at 700 kHz and 24.8 kA/m). This comprehensive parametric study provides the groundwork for the use of anionic colloidal nanocrystals to generate magnetically induced hyperthermia in various media, including complex systems and biological materials.     

Theoretical analysis of factors affecting the formation and stability of multilayered colloidal dispersions

A mathematical analysis of the major factors influencing the formation and stability of colloidal dispersions containing spherical particles surrounded by multilayered polymeric interfacial membranes formed by the layer-by-layer electrostatic deposition technique is carried out. The mathematical model assumes that (i) the colloidal dispersion initially consists of a mixture of electrically charged monodisperse spherical particles and oppositely charged polymer molecules, (ii) the adsorption of polymer molecules to the particle surfaces is diffusion-limited, and (iii) the dominant particle-particle collision mechanism is Brownian motion. This approach was used to produce stability maps that highlight conditions under which bridging flocculation, multilayer formation, or depletion flocculation occurs. The stability maps are derived from calculations of the critical polymer concentrations required to (i) saturate the particle surfaces (C(Sat)), (ii) ensure that polymer adsorption is faster than particle collisions (C(Ads)), and (iii) promote depletion flocculation (C(Dep)). In addition, the influence of interfacial properties on the stability of multilayer colloidal dispersions was assessed by calculating the colloidal interactions between the coated particles (i.e., van der Waals, electrostatic, steric, and depletion). These calculations indicated that the major factors are the interfacial charge and composition rather than the interfacial thickness. This article provides useful insights into the factors affecting the formation of stable multilayer colloidal dispersions.     

Effects of flow on measurements of interactions in colloidal suspensions

A hydrodynamic mechanism of interactions of colloidal particles is considered. The mechanism is based on the assumption of tiny background flows in the experimental cells during measurements by Grier et al. Both trivial (shear flow) and nontrivial (force propagation through viscous fluid) effects are taken into account for two colloidal particles near a wall bounding the solvent. Expressions for the radial (attractive or repulsive) forces and the polar torques are obtained. Quantitative estimates of the flow needed to produce the observed strength of attractive force are given; other necessary conditions are also considered. The following conclusion is made: the mechanism suggested most likely is not responsible for the attractive interactions observed in the experiments of Grier et al.; however, it may be applicable in other experimental realizations and should be kept in mind while conducting colloidal measurements of high sensitivity. Several distinctive features of the interactions due to this mechanism are identified.     

Mesoscopic dispersion of colloidal agglomerate in a complex fluid modelled by a hybrid fluid-particle model

The dispersion of the agglomerating fluid process involving colloids has been investigated at the mesoscale level by a discrete particle approach--the hybrid fluid-particle model (FPM). Dynamical processes occurring in the granulation of colloidal agglomerate in solvents are severely influenced by coupling between the dispersed microstructures and the global flow. On the mesoscale this coupling is further exacerbated by thermal fluctuations, particle-particle interactions between colloidal beds, and hydrodynamic interactions between colloidal beds and the solvent. Using the method of FPM, we have tackled the problem of dispersion of a colloidal slab being accelerated in a long box filled with a fluid. Our results show that the average size of the agglomerated fragments decreases with increasing shearing rate gamma, according to the power law A x gamma(k), where k is around 2. For larger values of gamma, the mean size of the agglomerate S(avg) increases slowly with gamma from the collisions between the aggregates and the longitudinal stretching induced by the flow. The proportionality constant A increases exponentially with the scaling factor of the attractive forces acting between the colloidal particles. The value of A shows a rather weak dependence on the solvent viscosity. But A increases proportionally with the scaling factor of the colloid-solvent dissipative interactions. Similar type of dependence can be found for the mixing induced by Rayleigh-Taylor instabilities involving the colloidal agglomerate and the solvent. Three types of fragmentation structures can be identified, which are called rupture, erosion, and shatter. They generate very complex structures with multiresolution character. The aggregation of colloidal beds is formed by the collisions between aggregates, which are influenced by the flow or by the cohesive forces for small dispersion energies. These results may be applied to enhance our understanding concerning the nonlinear complex interaction occurring in mesoscopic flows such as blood flow in small vessels.     

Thermodiffusion of interacting colloids. I. A statistical thermodynamics approach

The thermal diffusion coefficient of colloids consists of two additive contributions, one related to specific interactions between the surfaces of colloidal particles with solvent molecules, and a contribution due to interactions between the colloidal particles. In the present paper, the effect of intercolloidal particle interactions on their thermodiffusive behavior is discussed within a statistical thermodynamics framework. Transport coefficients are expressed in terms of the interaction potential between the colloidal spheres. A special feature of macromolecular systems is that this interaction potential is a potential of mean force, which is temperature dependent. It is shown that under certain conditions this implicit temperature dependence gives rise to negative Soret coefficients, that is, to diffusion of macromolecules to hot regions.     

Small-angle neutron scattering on a core-shell colloidal system: a contrast-variation study

Small-angle neutron scattering (SANS) measurements are reported on a sterically stabilized, core-shell colloidal system using contrast variation. Aqueous dispersions of polystyrene particles bearing grafted poly(ethylene glycol) (PEG) have been studied over a large range of particle concentrations and two different solvent conditions for the PEG polymer. SANS data are analyzed quantitatively by modeling the particles as core-shell colloids. In a good solvent and under particle contrast conditions, an effective hard-sphere interaction captures excluded-volume interactions up to high concentrations. Contrast variation, through isotopic substitution of both the core and solvent, expedite a detailed study of the PEG layer, both in the dilute limit and as a function of the particle concentration. Upon diminishing the solvent quality, subtle changes in the PEG layer translate into attractions among particles of moderate magnitude.     

Direct measurement of the binding force between microfabricated particles and a planar surface in aqueous solution by force-sensing piezoresistive cantilevers

We propose a force measurement method for evaluating the binding force between microscale flat surfaces in an aqueous solution. Using force-sensing piezoresistive cantilevers with sub-nanonewton force resolution, we have directly measured binding forces between SiO2-SiO2 microcontacts, which were created by gravity-driven random collision between microfabricated SiO2 cylindrical particles and a planar SiO2 substrate in a HCl solution. First, to examine our method we measured the pH dependence of the binding force. The binding forces were 12 and 5.8 nN at pH 1.0 and 2.0, respectively. As the pH increased, the binding force decreased and became zero at pH greater than 3.0. We confirmed that the bindings were based on the van der Waals' (VDW) force at pH 2.0 or less whereas a repulsive double-layer force acted between the surfaces at pH 3.0 or more. Second, the binding forces were categorized into a friction force or an adhesion force between the particles and the substrate. In the measurement, the friction force between the particle and the substrate was measured in the case when the particle slid on the substrate. On the contrary, the adhesion force was measured when the particle came off the substrate. Whether the particle slid or came off depended on the aspect ratio of the particle. We fabricated cylindrical particles with an aspect ratio of 0.03-2.0 and distinguished the friction force from the adhesion force by changing the aspect ratio of the particles. As a result, the friction force per unit contact area between SiO2-SiO2 flat surfaces was found to be 330 pN/microm2 +/- 20% when we used particles with a low aspect ratio (<0.1), and the adhesion force per unit contact area was 90 pN/microm2 +/- 20% for particles with a high aspect ratio (>0.4). For fluidic self-assembly that utilizes microscale surface contact in a liquid, our measurement method is an effective tool for studying and developing systems.     

Calculation of effective diffusion coefficient in a colloidal crystal by the finite-element method

The work is devoted to the calculation of effective diffusion coefficient of ions from the bulk solution to the electrode through a mask and the calculation of the distribution of the limiting current density over the electrode surface. A colloidal crystal, which is formed by orderly arranged monodispersed spherical particles, serves as a mask. It is shown that the diffusion of electroactive ions in the pores between spherical particles can be simulated by unit cells with rhombic, rectangular, or triangular cross-section. In the latter case, the cell side surface has no periodical boundaries. This simplifies significantly the numerical solution of the Laplace??s equation by the finite-element method. The effective diffusion coefficient in the bulk colloidal crystal is calculated at various values of its porosity. The calculated results agree well with the literature data. It is found that, for close-packed spherical particles, the relative effective diffusion coefficient in the bulk colloidal crystal is 0.16. The thicknesses of transient zones adjacent to the electrode surface and outer boundary of colloidal crystal and the effective diffusion coefficients for these zones are determined. The dependence of effective diffusion coefficient on the number of spherical particle layers in the colloidal crystal is obtained. The distribution of the limiting current density over the electrode surface is analyzed at various numbers of particle layers.     

Pair interaction potentials of colloids by extrapolation of confocal microscopy measurements of collective suspension structure

A method for measuring the pair interaction potential between colloidal particles by extrapolation measurement of collective structure to infinite dilution is presented and explored using simulation and experiment. The method is particularly well suited to systems in which the colloid is fluorescent and refractive index matched with the solvent. The method involves characterizing the potential of mean force between colloidal particles in suspension by measurement of the radial distribution function using 3D direct visualization. The potentials of mean force are extrapolated to infinite dilution to yield an estimate of the pair interaction potential, U(r). We use Monte Carlo simulation to test and establish our methodology as well as to explore the effects of polydispersity on the accuracy. We use poly-12-hydroxystearic acid-stabilized poly(methyl methacrylate) particles dispersed in the solvent dioctyl phthalate to test the method and assess its accuracy for three different repulsive systems for which the range has been manipulated by addition of electrolyte.     

The Correlation Functions of a Suspension of Large Particles in Amorphous Polybutadiene: A Molecular Dynamics Simulation Study

The correlation functions of a suspension of Lennard-Jones large particles in a model 1,4-cis-polybutadiene solvent have been investigated by molecular dynamics simulations. We present the effects of temperature, the degree of polymerization, and the solvent/large particle density on the effective interactions between the large particles and the solvent. It is found that as the temperature increases, the structure between solvent-solvent, solvent-large particle, and large-large particle decreases. Additionally, as the bulk density or the chain length is increased, the attractive part of the large-large interaction becomes weaker and small. We believe that in part, this is due to the polymer having collapsed onto itself and entangling the large particles and lessening their interaction until they are actually in contact. Increasing the length of the polymer also entangles the large particles. However, we believe that this may be a general feature that is a characteristic of a polymer solvent containing macroscopically large colloidal particles even though entanglement should be of less significance. Copyright 2000 Academic Press.     

Collisions, caging, thermodynamics, and jamming in the barrier hopping theory of glassy hard sphere fluids

An ultralocal limit of the microscopic single particle barrier hopping theory of glassy dynamics is proposed which allows explicit analytic expressions for the characteristic length scales, energy scales, and nonequilibrium free energy to be derived. All properties are shown to be controlled by a single coupling constant determined by the fluid density and contact value of the radial distribution function. This parameter quantifies an effective mean square force exerted on a tagged particle due to collisions with its surroundings. The analysis suggests a conceptual basis for previous surprising findings of multiple inter-relationships between characteristics of the transient localized state, the early stages of cage escape, non-Gaussian or dynamic heterogeneity effects, and the barrier hopping process that defines the alpha relaxation event. The underlying physical picture is also relevant to fluids of nonspherical molecules and sticky colloidal suspensions. The possibility of a unified view of liquid dynamics is suggested spanning the range from dense gases to the zero mobility jammed state.     

Measuring colloidal interactions with confocal microscopy

We use confocal laser scanning microscopy to measure interactions in colloidal suspensions. By inverting the radial distribution function, determined by tracking the particle coordinates, we obtain the effective interaction between the colloidal particles. Although this method can be applied to arbitrary colloidal interactions, here we demonstrate its efficacy with two well-known systems for which accurate theories are available: a colloid-polymer mixture and binary hard spheres. The high sensitivity of this method allows for the precise determination of complex interactions, as exemplified, for example, by the accurate resolution of the oscillatory effective potential of the binary hard sphere system. We argue that the method is particularly well suited for the determination of attractive forces.     

Nonequilibrium inertial dynamics of colloidal systems

We consider the properties of a one-dimensional fluid of Brownian inertial hard-core particles, whose microscopic dynamics is partially damped by a heat bath. Direct interactions among the particles are represented as binary, instantaneous elastic collisions. Collisions with the heat bath are accounted for by a Fokker-Planck collision operator, whereas direct collisions among the particles are treated by a well known method of kinetic theory, the revised Enskog theory. By means of a time multiple time-scale method we derive the evolution equation for the average density. Remarkably, for large values of the friction parameter and/or of the mass of the particles we obtain the same equation as the one derived within the dynamic density functional theory (DDF). In addition, at moderate values of the friction constant, the present method allows to study the inertial effects not accounted for by DDF method. Finally, a numerical test of these corrections is provided.     

Two-dimensional Monte Carlo simulations of a colloidal dispersion composed of rod-like ferromagnetic particles in the absence of an applied magnetic field

Influences of the magnetic interaction between particles and the aspect ratio of particles on aggregate structures in a colloidal dispersion composed of rod-like ferromagnetic particles were investigated by means of the cluster-moving Monte Carlo method. The internal structures of the aggregates obtained in simulations were analyzed in terms of the number density distribution of the clusters and radial distribution functions. The results show that as the magnetic interaction between particles increases, many small clusters such as anti-parallel particle pairs, raft-like clusters, and small loop-like clusters are formed; these gather to form larger aggregates. In the case of a relatively strong magnetic interaction between particles, solid-like rectangular clusters are formed when the aspect ratio is approximately 2, since the suitable distance between magnetic charges enables particles to form a fundamental structure of two normal anti-parallel particle pairs. As the aspect ratio increases beyond 2, many more stable raft-like clusters are formed, since the increase in distance between magnetic charges makes the two normal anti-parallel particle pair structures unstable. For a significantly larger aspect ratio, large network microstructures are produced by the formation of many chain-like and loop-like structures.     

熵驱大分子胶体粒子的柱状受限结晶

熵是物理化学的基本状态参量,在统计力学和热力学中处于核心位置.按照玻尔兹曼的微观解释,熵可以由孤立系统微观状态的数目(W)给出,即S=kBlnW,这里kB为玻尔兹曼常数[1,2].根据此公式,微观状态数越多,系统越混乱,熵越大,所以熵常被视作体系无序程度的度量.但熵增仅对应体系微观状态数的增加,与可观测的结构有序程度无关[3~5].在一些典型的软物质体系中,结构越有序熵反而越大,如胶体硬球在随机密堆积点的有序结晶[6]及描述各向异性棒状分子从各向同性相到向列相转变的Onsager原理[7].