Role of anisotropy in electrodynamically induced colloidal aggregates

We investigate the assembly of spherical and anisotropic colloidal particles with the shape of peanuts when subjected to an external alternating electric field. By varying the strength and frequency of the applied field, we observe that both types of particles form clusters at low frequencies due to attractive electrohydrodynamic interactions or disperse into a liquidlike phase at high frequencies due to repulsive dipolar interactions. We characterize the observed structures via pair correlation functions and radius of gyration, and observe a clear difference in the ordering process between the isotropic and anisotropic colloids. Further on, we interpret the cluster formation kinetics in terms of dynamic scaling theory, and observe a faster aggregation of the anisotropic colloids with respect to the isotropic ones.     

计算平均力势的理论方法

A universal theoretical framework is proposed for calculating potential of mean force (PMF) between two solute particles immersed in a solvent bath, the present method overcomes all of drawbacks of previous methods. The only input required to implement the recipe is solvent density distribution profile around a single solute particle. The universal framework is applied to calculate the PMF between two large spherical particles immersed in small hard sphere solvent bath. Comparison between the present predictions and existing simulation data shows reliability of the present recipe. Effects of solvent-solute interaction detail, solvent bulk density, and solute size on the excess PMF are investigated. The resultant conclusion is that depletion of solvent component by the solute particle induces attractive excess PMF, while gathering of solvent component by the solute particle induces repulsive excess PMF, high solvent bulk density and large solute size can strengthen the tendency of attraction or repulsion. Relevance of transition from depletion attraction to gathering repulsion with the biomolecular interaction, i.e. hydrophobic attraction and hydration repulsion, is discussed.     

Density functional theory of freezing of a system of conjugated oligomers parameterised via Gay–Berne potential

Density functional theory (DFT) of freezing is used to study the isotropic–nematic, isotropic–smectic A and nematic–smectic A phase transitions in a system of large, semi-flexible conjugated oligomers parameterised within Gay–Berne (GB) potential. The pair correlation functions of the isotropic fluid, used as structural inputs in the DFT, are calculated by solving the Percus–Yevick integral equation theory. Large number of spherical harmonic coefficients of each orientation-dependent functions has been considered to ensure the numerical accuracy at different densities and temperatures for the system of these model GB ellipsoids having large aspect ratio (length-to-breadth ratio). We found that the system of GB ellipsoids parameterised for conjugated oligomers shows stable isotropic, nematic and smectic A phases. At low temperatures, on increasing the density, isotropic fluid makes a direct transition to smectic A phase. Nematic phase get stabilised in between the isotropic and smectic A phases on increasing the temperature. Using the transition parameter obtained through the DFT, we have plotted the temperature–density and pressure–temperature phase diagrams which are found to be qualitatively similar to the one obtained in simulations for the systems with low aspect ratio GB particles.     

Free isotropic-nematic interfaces in fluids of charged platelike colloids

Bulk properties and free interfaces of mixtures of charged platelike colloids and salt are studied within the density-functional theory. The particles are modeled by hard cuboids with their edges constrained to be parallel to the Cartesian axes corresponding to the Zwanzig model. The charges of the particles are concentrated in their center. The density functional is derived by functional integration of an extension of the Debye-Hückel pair distribution function with respect to the interaction potential. For sufficiently small macroion charges, the bulk phase diagrams exhibit one isotropic and one nematic phase separated by a first-order phase transition. With increasing platelet charge, the isotropic and nematic binodals are shifted to higher densities. The Donnan potential between the coexisting isotropic and nematic phases is inferred from bulk structure calculations. Nonmonotonic density and nematic order parameter profiles are found at a free interface interpolating between the coexisting isotropic and nematic bulk phases. Moreover, electrically charged layers form at the free interface leading to monotonically varying electrostatic potential profiles. Both the widths of the free interfaces and the bulk correlation lengths are approximately given by the Debye length. For fixed salt density, the interfacial tension decreases upon increasing the macroion charge.     

Nematic-fluid structure in wall-field geometry. II. The direct correlation function

An explicit expression for the wall-nematic direct correlation function (DCF) is obtained for any orientation of the wall with respect to an external orienting field. It is found that inside the surface of the wall, the DCF rapidly tends to a function of the nematogen orientation and depends only on parameters of the bulk fluid. We suggest that the wall-nematic DCF can be used as an ansatz for the colloid-nematic DCF in dilute nematic colloids. The reliability of this ansatz is investigated at different field strengths in both isotropic and nematic regions. Our calculations for spherical colloidal particles show that this approximation is valid for colloidal particles that are large, but well within the physically realistic size range. The ansatz could also be applied to nonspherical colloidal particles.     

Influence of the anisometry of magnetic particles on the isotropic–nematic phase transition

The influence of the shape anisotropy of magnetic particles on the isotropic–nematic phase transition was studied in ferronematics based on the nematic liquid crystal (LC) 4-(trans-4-n-hexylcyclohexyl)-isothiocyanato-benzene (6CHBT). The LC was doped with spherical or rod-like magnetic particles of different size and volume concentrations. The phase transition from isotropic to nematic phase was observed by polarising microscope as well as by capacitance measurements. The influence of the concentration and the shape anisotropy of the magnetic particles on the isotropic–nematic phase transition in LC are demonstrated here. The results are in a good agreement with recent theoretical predictions.     

Nanospheres in a nematic liquid crystal solvent: the influence of particle size

We use Monte Carlo simulations to investigate a simple lattice model for nematic liquid crystals containing nanospheres. The influence of particle size on the phase behaviour is studied using two different sized particles. The phase diagram is found to be topologically equivalent for both particle sizes, with a large biphasic region corresponding to coexistence between a rod-rich nematic and a rod-poor isotropic phase. For small spheres, the rod-rich nematic phase is stable for relatively large volume fractions of spheres (up to a maximum of about 16%). In contrast, the nematic phase for the system with larger spheres is constrained to a much narrower region of the phase diagram.     

Mesophase formation in a system of top-shaped hard molecules: density functional theory and Monte Carlo simulation

We present the phase diagram of a system of mesogenic top-shaped molecules based on the Parsons-Lee density functional theory and Monte Carlo simulation. The molecules are modeled as a hard spherocylinder with a hard sphere embedded in its center. The stability of five different phases is studied, namely, isotropic, nematic, smectic A, smectic C, and columnar phases. The positionally ordered phases are investigated only for the case of parallel alignment. It is found that the central spherical unit destabilizes the nematic with respect to the isotropic phase, while increasing the length of the cylinder has the opposite effect. Also, the central hard sphere has a strong destabilizing effect on the smectic A phase, due the inefficient packing of the molecules into layers. For large hard sphere units the smectic A phase is completely replaced by a smectic C structure. The columnar phase is first stabilized with increasing diameter of the central unit, but for very large hard sphere units it becomes less stable again. The density functional results are in good agreement with the simulations.     

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.     

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.     

Steric stabilization of spherical colloidal particles: Implicit and explicit solvent

We present the results of Monte Carlo simulations and density functional theory treatment of interactions between spherical colloidal brushes both in implicit (good) solvent and in an explicit polymeric solution. Overall, theory is seen to be in good agreement with simulations. We find that interactions between hard-sphere particles grafted with hard-sphere chains are always repulsive in implicit solvent. The range and steepness of the repulsive interaction is sensitive to the grafting density and the length of the grafted chains. When the brushes are immersed in an explicit solvent of hard-sphere chains, a weak mid-range attraction arises, provided the length of the free chains exceeds that of the grafted chains.     

A density-functional theory study of the confined soft ellipsoid fluid

A system of soft ellipsoid molecules confined between two planar walls is studied using classical density-functional theory. Both the isotropic and nematic phases are considered. The excess free energy is evaluated using two different Ans?tze and the intermolecular interaction is incorporated using two different direct correlation functions (DCF's). The first is a numerical DCF obtained from simulations of bulk soft ellipsoid fluids and the second is taken from the Parsons-Lee theory. In both the isotropic and nematic phases the numerical DCF gives density and order parameter profiles in reasonable agreement with simulation. The Parsons-Lee DCF also gives reasonable agreement in the isotropic phase but poor agreement in the nematic phase.     

Interpretation of conservative forces from Stokesian dynamic simulations of interfacial and confined colloids

This paper presents Stokesian dynamics simulations of experiments involving one or two charged colloids near either a single charged wall or confined between parallel charged walls. Equilibrium particle-particle and particle-wall interactions are interpreted from dynamic particle trajectories in simulations involving (1) a single particle levitated above a wall, (2) two particles below a wall, and (3) two particles confined between two parallel walls. By specifying only repulsive electrostatic Derjaguin-Landau-Verwey-Overbeek (DLVO) potentials and including multibody hydrodynamics, we successfully recover expected potentials in some cases, while anomalous attraction is observed in other cases. Attraction inferred in the latter simulations displays quantitative agreement with literature measurements when particle dynamics are interpreted using reported analyses. Because anomalous attraction is reproduced in simulations using only electrostatic repulsive DLVO potentials, our results reveal the one-dimensional analyses to be invalid for configurations that are inherently multidimensional via multibody hydrodynamics. Parameters related to experimental sampling of particle dynamics are also found to be critical for obtaining accurate potentials. We explain the anomalous attraction in each experiment using effective potentials, which can be employed in an a priori fashion to assist the confident design of future experiments involving interfacial and confined colloids. Ultimately, our findings reveal the importance of dimensionality and multibody hydrodynamics for understanding nonequilibrium dynamics of colloids near surfaces.     

Surface properties of fluids of charged platelike colloids

Surface properties of mixtures of charged platelike colloids and salt in contact with a charged planar wall are studied within density functional theory. The particles are modeled by hard cuboids with their edges constrained to be parallel to the Cartesian axes corresponding to the Zwanzig model [J. Chem. Phys. 39, 1714 (1963)] and the charges of the particles are concentrated at their centers. The density functional applied is an extension of a recently introduced functional for charged platelike colloids. It provides a qualitative approach because it does not determine the relation between the actual and the effective charges entering into the model. Technically motivated approximations, such as using the Zwanzig model, are expected not to influence the results qualitatively. Analytically and numerically calculated bulk and surface phase diagrams exhibit first-order wetting for sufficiently small macroion charges and isotropic bulk order as well as first-order drying for sufficiently large macroion charges and nematic bulk order. The asymptotic wetting and drying behaviors are investigated by means of effective interface potentials which turn out to be asymptotically the same as for a suitable neutral system governed by isotropic nonretarded dispersion forces. Wetting and drying points as well as predrying lines and the corresponding critical points have been located numerically. A crossover from monotonic to nonmonotonic electrostatic potential profiles upon varying the surface charge density has been observed. Nonmonotonic electrostatic potential profiles are equivalent to the occurrence of charge inversion. Due to the presence of both the Coulomb interactions and the hard-core repulsions, the surface potential and the surface charge do not vanish simultaneously, i.e., the point of zero charge and the isoelectric point of the surface do not coincide.     

Effect of charge inhomogeneity and mobility on colloid aggregation

The aggregation of inhomogeneously charged colloids with the same average charge is analyzed using Monte Carlo simulations. We find aggregation of colloids for sizes in the range 10-200 nm, which is similar to the range in which aggregation is observed in several experiments. The attraction arises from the strongly correlated electrostatic interactions associated with the increase in the counterion density in the region between the particles; this effect is enhanced by the discreteness and mobility of the surface charges. Larger colloids attract more strongly when their surface charges are discrete. We study the aggregation as functions of the surface charge density, counterion valence, and volume fraction.     

Percolation, phase separation, and gelation in fluids and mixtures of spheres and rods

The relationship between kinetic arrest, connectivity percolation, structure and phase separation in protein, nanoparticle, and colloidal suspensions is a rich and complex problem. Using a combination of integral equation theory, connectivity percolation methods, nai?ve mode coupling theory, and the activated dynamics nonlinear Langevin equation approach, we study this problem for isotropic one-component fluids of spheres and variable aspect ratio rigid rods, and also percolation in rod-sphere mixtures. The key control parameters are interparticle attraction strength and its (short) spatial range, total packing fraction, and mixture composition. For spherical particles, formation of a homogeneous one-phase kinetically stable and percolated physical gel is predicted to be possible, but depends on non-universal factors. On the other hand, the dynamic crossover to activated dynamics and physical bond formation, which signals discrete cluster formation below the percolation threshold, almost always occurs in the one phase region. Rods more easily gel in the homogeneous isotropic regime, but whether a percolation or kinetic arrest boundary is reached first upon increasing interparticle attraction depends sensitively on packing fraction, rod aspect ratio and attraction range. Overall, the connectivity percolation threshold is much more sensitive to attraction range than either the kinetic arrest or phase separation boundaries. Our results appear to be qualitatively consistent with recent experiments on polymer-colloid depletion systems and brush mediated attractive nanoparticle suspensions.     

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.     

Nanoparticles in nematic liquid crystals: interactions with nanochannels

A mesoscale theory for the tensor order parameter Q is used to investigate the structures that arise when spherical nanoparticles are suspended in confined nematic liquid crystals (NLCs). The NLC is "sandwiched" between a wall and a small channel. The potential of mean force is determined between particles and the bottom of the channels or between several particles. Our results suggest that strong NLC-mediated interactions between the particles and the sidewalls of the channels, on the order of hundreds of k(B)T, arise when the colloids are inside the channels. The magnitude of the channel-particle interactions is dictated by a combination of two factors, namely, the type of defect structures that develop when a nanoparticle is inside a channel, and the degree of ordering of the nematic in the region between the colloid and the nanochannel. The channel-particle interactions become stronger as the nanoparticle diameter becomes commensurate with the nanochannel width. Nanochannel geometry also affects the channel-particle interactions. Among the different geometries considered, a cylindrical channel seems to provide the strongest interactions. Our calculations suggest that small variations in geometry, such as removing the sharp edges of the channels, can lead to important reductions in channel-particle interactions. Our calculations for systems of several nanoparticles indicate that linear arrays of colloids with Saturn ring defects, which for some physical conditions are not stable in a bulk system, can be stabilized inside the nanochannels. These results suggest that nanochannels and NLCs could be used to direct the assembly of nanoparticles into ordered arrays with unusual morphologies.     

Dynamics of self-assembled chaining in magnetorheological fluids

The aggregation dynamics of paramagnetic spherical particles embedded in a viscous fluid is investigated via numerical simulations using a fully coupled three-dimensional model. Particles experience simultaneously Brownian motion, dipolar magnetic attraction, and multibody hydrodynamic interactions. When the dipole strength characterizing the ratio of magnetic attraction to random diffusion exceeds a critical value, particles join together forming supraparticle structures. As time evolves, particle/chain and chain/ chain interactions lead to a continuous increase of the cluster size. The mean length of particle chains has a power-law dependence with respect to time, as predicted by the theory of diffusion-limited aggregation. Both the exponent and the characteristic time scale agree very well with the experimental results of Promislow et al.     

Forces between cylindrical nanoparticles in a liquid crystal

Using classical density functional theory, the forces between two cylindrical nanoparticles in a liquid crystal solvent are calculated. Both the nematic and isotropic phases of the solvent are considered. In the nematic phase, the interaction is highly anisotropic. At short range, changes in the defect structure around the cylinders leads to a complex interaction between them. In the isotropic phase, an attractive interaction arises due to overlap between halos of ordered fluid adsorbed on the surfaces of the cylinders.