Effect of Hydrodynamic Interactions on Self-Diffusion of Quasi-Two-Dimensional Colloidal Hard Spheres

Rotational diffusion of colloidal microspheres near flat walls

Recently, colloids with an off-center fluorescent core and homogeneous composition have been developed to measure the rotational diffusivity of microparticles using 3D confocal microscopy in refractive index-matched suspensions. Here, we show that the same particles may be imaged using a standard fluorescence microscope to yield their rotational diffusion coefficients. Trajectories of the off-center core may be combined with known expressions for the correlation decay of particle orientations to determine an effective rotational diffusivity. For sedimented particles, we also find the rotational diffusivity about axes perpendicular and parallel to the interface by adding some bright field illumination and simultaneously tracking both the core and the particle. Trajectories for particles of different sizes yield excellent agreement with hydrodynamic models of rotational diffusion near flat walls, taking the sedimentation-diffusion equilibrium into account. Finally, we explore the rotational diffusivity of particles in crowded two-dimensional monolayers, finding a different reduction of the rotational motion about the two axes depending on the colloidal microstructure.

Evidence for Enhanced Tracer Diffusion in Densely Packed Interfacial Assemblies of Hairy Nanoparticles

Nearly monodisperse nanoparticle (NP) spheres attached to a nonvolatile ionic liquid surface were tracked by in situ scanning electron microscopy to obtain the tracer diffusion coefficient Dtr as a function of the areal fraction ϕ. The in situ technique resolved both tracer (gold) and background (silica) particles for ∼1-2 min, highlighting their mechanisms of diffusion, which were strongly dependent on ϕ. Structure and dynamics at low and moderate ϕ paralleled those reported for larger colloidal spheres, showing an increase in order and a decrease in Dtr by over 4 orders of magnitude. However, ligand interactions were more important near jamming, leading to different caging and jamming dynamics for smaller NPs. The normalized Dtr at ultrahigh ϕ depended on particle diameter and ligand molecular weight. Increasing the PEG molecular weight by a factor of 4 increased Dtr by 2 orders of magnitude at ultrahigh ϕ, indicating stronger ligand lubrication for smaller particles.

Universal Relation between Entropy and Kinetics

Relating thermodynamic and kinetic properties is a conceptual challenge with many practical benefits. Here, based on first principles, we derive a rigorous inequality relating the entropy and the dynamic propagator of particle configurations. It is universal and applicable to steady states arbitrarily far from thermodynamic equilibrium. Applying the general relation to diffusive dynamics yields a relation between the entropy and the (normal or anomalous) diffusion coefficient. The relation can be used to obtain useful bounds for the late-time diffusion coefficient from the calculated steady-state entropy or, conversely, to estimate the entropy based on measured diffusion coefficients. We demonstrate the validity and usefulness of the relation through several examples and discuss its broad range of applications, in particular, for systems far from equilibrium.

Collisions Enhance Self-Diffusion in Odd-Diffusive Systems

It is generally believed that collisions of particles reduce the self-diffusion coefficient. Here we show that in odd-diffusive systems, which are characterized by diffusion tensors with antisymmetric elements, collisions surprisingly can enhance the self-diffusion. In these systems, due to an inherent curving effect, the motion of particles is facilitated, instead of hindered by collisions leading to a mutual rolling effect. Using a geometric model, we analytically predict the enhancement of the self-diffusion coefficient with increasing density. This counterintuitive behavior is demonstrated in the archetypal odd-diffusive system of Brownian particles under Lorentz force. We validate our findings by many-body Brownian dynamics simulations in dilute systems.

Observing capture with a colloidal model membrane channel

We use video microscopy to study the full capture process for colloidal particles transported through microfluidic channels by a pressure-driven flow. In particular, we obtain trajectories for particles as they move from the bulk into confinement, using these to map in detail the spatial velocity and concentration fields for a range of different flow velocities. Importantly, by changing the height profiles of our microfluidic devices, we consider systems for which flow profiles in the channel are the same, but flow fields in the reservoir differ with respect to the quasi-2D monolayer of particles. We find that velocity fields and profiles show qualitative agreement with numerical computations of pressure-driven fluid flow through the systems in the absence of particles, implying that in the regimes studied here particle-particle interactions do not strongly perturb the flow. Analysis of the particle flux through the channel indicates that changing the reservoir geometry leads to a change between long-range attraction of the particles to the pore and diffusion-to-capture-like behaviour, with concentration fields that show qualitative changes based on device geometry. Our results not only provide insight into design considerations for microfluidic devices, but also a foundation for experimental elucidation of the concept of a capture radius. This long standing problem plays a key role in transport models for biological channels and nanopore sensors.

Are strongly confined colloids good models for two dimensional liquids?

Quasi-two-dimensional (quasi-2D) colloidal hard-sphere suspensions confined in a slit geometry are widely used as two-dimensional (2D) model systems in experiments that probe the glassy relaxation dynamics of 2D systems. However, the question to what extent these quasi-2D systems indeed represent 2D systems is rarely brought up. Here, we use computer simulations that take into account hydrodynamic interactions to show that dense quasi-2D colloidal bi-disperse hard-sphere suspensions exhibit much more rapid diffusion and relaxation than their 2D counterparts at the same area fraction. This difference is induced by the additional vertical space in the quasi-2D samples in which the small colloids can move out of the 2D plane, therefore allowing overlap between particles in the projected trajectories. Surprisingly, this difference in the dynamics can be accounted for if, instead of using the surface density, one characterizes the systems by means of a suitable structural quantity related to the radial distribution function. This implies that in the two geometries, the relevant physics for glass formation is essentially identical. Our results provide not only practical implications on 2D colloidal experiments but also interesting insights into the 3D-to-2D crossover in glass-forming systems.

Hydrodynamics induce superdiffusive jumps of passive tracers along critical paths of random networks and colloidal gels

We present a numerical study on the effect of hydrodynamic interactions (HI) on the diffusion of inert point tracer particles in several fixed random structures. As expected, the diffusion is hampered by the extra hydrodynamic friction introduced by the obstacle network. However, a non-trivial effect due to HI appears in the analysis of the van-Hove displacement probability close to the percolation threshold, where tracers diffuse through critical fractal paths. We show that the tracer dynamics can be split up into short and long jumps, the latter being ruled by either exponential or Gaussian van Hove distribution tails. While at short time HI slow down the tracer diffusion, at long times, hydrodynamic interactions with the obstacles increase the probability of longer jumps, which circumvent the traps of the labyrinth more easily. Notably, the relation between the anomalous diffusion exponent and the fractal dimension of the critical (intricate) paths is greater than one, which implies that the long-time (long-jump) diffusion is mildly superdiffuse. A possible reason for such a hastening of the diffusion along the network corridors is the hydrodynamically induced mobility anisotropy, which favours displacements parallel to the walls, an effect which has already been experimentally observed in collagen gels.

Spatial Crossover Between Far-From-Equilibrium and Near-Equilibrium Dynamics in Locally Driven Suspensions

We examine the response of a quasi-two-dimensional colloidal suspension to a localized circular driving induced by optical tweezers. This approach allows us to resolve over 3 orders of magnitude in the Péclet number (Pe) and provide a direct observation of a sharp spatial crossover from far- to near-thermal-equilibrium regions of the suspension. In particular, particles migrate from high to low Pe regions and form strongly inhomogeneous steady-state density profiles with an emerging length scale that does not depend on the particle density and is set by Pe≈1. We show that the phenomenological two phase fluid constitutive model is in line with our results.

Accelerated annealing of colloidal crystal monolayers by means of cyclically applied electric fields

External fields are commonly applied to accelerate colloidal crystallization; however, accelerated self-assembly kinetics can negatively impact the quality of crystal structures. We show that cyclically applied electric fields can produce high quality colloidal crystals by annealing local disorder. We find that the optimal off-duration for maximum annealing is approximately one-half of the characteristic melting half lifetime of the crystalline phase. Local six-fold bond orientational order grows more rapidly than global scattering peaks, indicating that local restructuring leads global annealing. Molecular dynamics simulations of cyclically activated systems show that the ratio of optimal off-duration for maximum annealing and crystal melting time is insensitive to particle interaction details. This research provides a quantitative relationship describing how the cyclic application of fields produces high quality colloidal crystals by cycling at the fundamental time scale for local defect rearrangements; such understanding of dynamics and kinetics can be applied for reconfigurable colloidal assembly.

Diffusion of colloidal particles in model porous media

Using video microscopy and simulations, we study the long-time diffusion of colloidal tracers in a wide range of model porous media composed of frozen colloidal matrices with different structures. We found that the diffusion coefficient of a tracer can be quantitatively determined by the structures of porous media. In particular, a universal scaling relation exists between the dimensionless diffusion coefficient of the tracer and the structural entropy of the system. This universal scaling relation is an extension of the scaling law previously discovered for the diffusion of colloidal particles in fluctuating media.

Pattern detection in colloidal assembly: A mosaic of analysis techniques

Characterization of the morphology, identification of patterns and quantification of order encountered in colloidal assemblies is essential for several reasons. First of all, it is useful to compare different self-assembly methods and assess the influence of different process parameters on the final colloidal pattern. In addition, casting light on the structures formed by colloidal particles can help to get better insight into colloidal interactions and understand phase transitions. Finally, the growing interest in colloidal assemblies in materials science for practical applications going from optoelectronics to biosensing imposes a thorough characterization of the morphology of colloidal assemblies because of the intimate relationship between morphology and physical properties (e.g. optical and mechanical) of a material. Several image analysis techniques developed to investigate images (acquired via scanning electron microscopy, digital video microscopy and other imaging methods) provide variegated and complementary information on the colloidal structures under scrutiny. However, understanding how to use such image analysis tools to get information on the characteristics of the colloidal assemblies may represent a non-trivial task, because it requires the combination of approaches drawn from diverse disciplines such as image processing, computational geometry and computational topology and their application to a primarily physico-chemical process. Moreover, the lack of a systematic description of such analysis tools makes it difficult to select the ones more suitable for the features of the colloidal assembly under examination. In this review we provide a methodical and extensive description of real-space image analysis tools by explaining their principles and their application to the investigation of two-dimensional colloidal assemblies with different morphological characteristics.

Long-time self-diffusion in quasi-two-dimensional colloidal fluids of paramagnetic particles

The effect of hydrodynamic interactions (HI) on the long-time self-diffusion in quasi-two-dimensional fluids of paramagnetic colloidal particles is investigated using a combination of experiments and Brownian dynamics (BD) simulations. In the BD simulations, the direct interactions (DI) between the particles consist of a short-ranged repulsive part and a long-ranged part that is proportional to 1/r^{3}, with r the interparticle distance. By studying the equation of state, the simulations allow for the identification of the regime where the properties of the fluid are fully controlled by the long-ranged interactions, and the thermodynamic state solely depends on the dimensionless interaction strength Γ. In this regime, the radial distribution functions from the simulations are in quantitative agreement with those from the experiments for different fluid area fractions. This agreement confirms that the DI in the experiments and simulations are identical, which thus allows us to isolate the role of HI, as these are not taken into account in the BD simulations. Experiment and simulation fall onto a master curve with respect to the Γ dependence of D_{L}^{★}=D_{L}/(D_{0}Γ^{1/2}), with D_{0} the self-diffusion coefficient at infinite dilution and D_{L} the long-time self-diffusion coefficient. Our results thus show that, although HI affect the short-time self-diffusion, for a quasi-two-dimensional system with 1/r^{3} long-ranged DI, the reduced quantity D_{L}^{★} is effectively not affected by HI. Interestingly, this is in agreement with prior work on quasi-two-dimensional colloidal hard spheres.

Universal Scaling Law for Colloidal Diffusion in Complex Media

Using video microscopy and simulations, we study the diffusion of probe particles in a wide range of complex backgrounds, both crystalline and disordered, in quasi-2D colloidal systems. The dimensionless diffusion coefficients D^{*} from different systems collapse to a single master curve when plotted as a function of the structural entropy of the backgrounds, confirming the universal relation between diffusion dynamics and the structure of the medium. A new scaling equation is proposed with consideration for the viscous friction from the solvent, which is absent in previous theoretical models. This new universal law quantitatively predicts the diffusion coefficients from different systems over several orders of magnitude of D^{*}, with a single common fitting parameter.

Excess entropy and long-time diffusion in colloidal fluids with short-range interparticle attraction

Liquid structure and dynamics are experimentally investigated in colloidal suspensions with short-range depletion attraction. The colloidal fluid samples consist of hard-sphere colloidal particles suspended along with rodlike depletants based on surfactant micelles. The spheres have a range of surface chemistries, diameters, and packing fractions, and the rodlike micelle length depends on the temperature. Thus, the combination of hard-spheres and depletants generates a sample wherein short-range interparticle attraction can be temperature-tuned in situ. Video optical microscopy and particle tracking techniques are employed to measure particle trajectories from which structural and dynamical quantities are derived, including the particle pair correlation function [g(r)], mean square displacement, long-time diffusion coefficient, and the sample two-body excess entropy (S₂). The samples with stronger short-range attractions exhibit more order, as characterized by g(r) and S₂. The stronger short-range attractions are also observed to lead to slower long-time diffusion and more heterogeneous dynamics at intermediate time scales. Finally, the excess entropy scaling law prediction, i.e., the exponential relationship between two-body excess entropy and long-time diffusivity, is observed across the full range of samples.

Bond-orientational order and Frank's constant in two-dimensional colloidal hard spheres

Recently, the full phase behaviour of 2D colloidal hard spheres was experimentally established, and found to involve a first order liquid to hexatic transition and a continuous hexatic to crystal transition (Thorneywork et al 2017 Phys. Rev. Lett. 118 158001). Here, we expand upon this work by considering the behaviour of the bond-orientational correlation time and Frank's constant in the region of these phase transitions. We also consider the excess entropy, as calculated from the radial distribution functions, for a wide range of area fractions covering the liquid, hexatic and crystal phases. In all cases, the behaviour of these quantities further corroborates the previously reported melting scenario of 2D colloidal hard spheres.

2D Stokesian Approach to Modeling Flow Induced Deformation of Particle Laden Interfaces

A Stokesian dynamics simulation of the effect of surface Couette flow on the microstructure of particles irreversibly adsorbed to an interface is presented. Rather than modeling both bulk phases, the interface, and particles in a full 3D simulation, known interfacial interactions between adsorbed particles are used to create a 2D model from a top down perspective. This novel methodology is easy to implement and computationally inexpensive, which makes it favorable to simulate behavior of particles under applied flow at fluid-fluid interfaces. The methodology is used to examine microstructure deformation of monodisperse, rigid spherical colloids with repulsive interactions when a surface Couette flow is imposed. Simulation results compare favorably to experimental results taken from literature, showing that interparticle forces must be 1 order of magnitude greater than viscous drag for microstructure to transition from aligned particle strings to rotation of local hexagonal domains. Additionally, it is demonstrated that hydrodynamic interactions between particles play a significant role in the magnitude of these microstructure deformations.

Experimental Evidence for a Structural-Dynamical Transition in Trajectory Space

Among the key insights into the glass transition has been the identification of a nonequilibrium phase transition in trajectory space which reveals phase coexistence between the normal supercooled liquid (active phase) and a glassy state (inactive phase). Here, we present evidence that such a transition occurs in experiments. In colloidal hard spheres, we find a non-Gaussian distribution of trajectories leaning towards those rich in locally favored structures (LFSs), associated with the emergence of slow dynamics. This we interpret as evidence for a nonequilibrium transition to an inactive LFS-rich phase. Reweighting trajectories reveals a first-order phase transition in trajectory space between a normal liquid and a LFS-rich phase. We also find evidence for a purely dynamical transition in trajectory space.

Dynamic heterogeneities and non-Gaussian behavior in two-dimensional randomly confined colloidal fluids

A binary mixture of superparamagnetic colloidal particles is confined between glass plates such that the large particles become fixed and provide a two-dimensional disordered matrix for the still mobile small particles, which form a fluid. By varying fluid and matrix area fractions and tuning the interactions between the superparamagnetic particles via an external magnetic field, different regions of the state diagram are explored. The mobile particles exhibit delocalized dynamics at small matrix area fractions and localized motion at high matrix area fractions, and the localization transition is rounded by the soft interactions [T. O. E. Skinner et al., Phys. Rev. Lett. 111, 128301 (2013)PRLTAO0031-900710.1103/PhysRevLett.111.128301]. Expanding on previous work, we find the dynamics of the tracers to be strongly heterogeneous and show that molecular dynamics simulations of an ideal gas confined in a fixed matrix exhibit similar behavior. The simulations show how these soft interactions make the dynamics more heterogeneous compared to the disordered Lorentz gas and lead to strong non-Gaussian fluctuations.

Two-Dimensional Melting of Colloidal Hard Spheres

We study the melting of quasi-two-dimensional colloidal hard spheres by considering a tilted monolayer of particles in sedimentation-diffusion equilibrium. In particular, we measure the equation of state from the density profiles and use time-dependent and height-resolved correlation functions to identify the liquid, hexatic, and crystal phases. We find that the liquid-hexatic transition is first order and that the hexatic-crystal transition is continuous. Furthermore, we directly measure the width of the liquid-hexatic coexistence gap from the fluctuations of the corresponding interface, and thereby experimentally establish the full phase behavior of hard disks.

Combining Diffusion NMR and Small-Angle Neutron Scattering Enables Precise Measurements of Polymer Chain Compression in a Crowded Environment

The effect of particles on the behavior of polymers in solution is important in a number of important phenomena such as the effect of "crowding" proteins in cells, colloid-polymer mixtures, and nanoparticle "fillers" in polymer solutions and melts. In this Letter, we study the effect of spherical inert nanoparticles (which we refer to as "crowders") on the diffusion coefficient and radius of gyration of polymers in solution using pulsed-field-gradient NMR and small-angle neutron scattering (SANS), respectively. The diffusion coefficients exhibit a plateau below a characteristic polymer concentration, which we identify as the overlap threshold concentration c^{⋆}. Above c^{⋆}, in a crossover region between the dilute and semidilute regimes, the (long-time) self-diffusion coefficients are found, universally, to decrease exponentially with polymer concentration at all crowder packing fractions, consistent with a structural basis for the long-time dynamics. The radius of gyration obtained from SANS in the crossover regime changes linearly with an increase in polymer concentration, and must be extrapolated to c^{⋆} in order to obtain the radius of gyration of an individual polymer chain. When the polymer radius of gyration and crowder size are comparable, the polymer size is very weakly affected by the presence of crowders, consistent with recent computer simulations. There is significant chain compression, however, when the crowder size is much smaller than the polymer radius gyration.

Self-diffusion in two-dimensional binary colloidal hard-sphere fluids

We present a systematic experimental study of the dynamic behavior of monodisperse and bidisperse two-dimensional colloidal hard-sphere fluids. We consider the diffusive behavior of the two types of particles for systems with a variety of compositions and total area fractions. In particular, we measure the short- and long-time diffusion coefficients for both species independently. We find that the short-time self-diffusion coefficients show an approximately linear dependence on the area fraction and that the long-time self-diffusion coefficients are well described by an expression dependent upon only the area fraction and contact value of the radial distribution function. Finally, we consider the effect of composition change and find some variation in the long-time self-diffusion coefficients, which we ascribe to the complex packing effects exhibited by binary systems.

Colloidal diffusion in confined geometries

Colloidal diffusion in confined geometries is analysed at the level of anisotropic pair densities.

Short-time dynamics of monomers and dimers in quasi-two-dimensional colloidal mixtures

We report on the short-time dynamics in colloidal mixtures made up of monomers and dimers highly confined between two glass plates. At low concentrations, the experimental measurements of colloidal motion agree well with the solution of the Navier-Stokes equation at low Reynolds numbers; the latter takes into account the increase in the drag force on a colloidal particle due to wall-particle hydrodynamic forces. More importantly, we find that the ratio of the short-time diffusion coefficient of the monomer and that of the center of mass of the dimmer is almost independent of both the dimer molar fraction, x_{d}, and the total packing fraction, ϕ, up to ϕ≈0.5. At higher concentrations, this ratio displays a small but systematic increase. A similar physical scenario is observed for the ratio between the parallel and the perpendicular components of the short-time diffusion coefficients of the dimer. This dynamical behavior is corroborated by means of molecular dynamics computer simulations that include explicitly the particle-particle hydrodynamic forces induced by the solvent. Our results suggest that the effects of colloid-colloid hydrodynamic interactions on the short-time diffusion coefficients are almost identical and factorable in both species.

On the Gaussian approximation in colloidal hard sphere fluids

We study the behaviour of the self-intermediate scattering function and self-van Hove correlation function for quasi-two-dimensional colloidal hard sphere fluids at a range of area fractions. We compute these functions first directly from the particle coordinates and secondly from the mean squared displacement via the Gaussian approximation. This allows us to test the validity of this approximation over a range of length and time scales, where we find that the Gaussian approximation holds if the hydrodynamic limits are appropriately probed. Surprisingly, only small deviations from Gaussian behaviour are seen at intermediate times, even for dense fluids. We next consider these deviations from Gaussian behaviour firstly via the non-Gaussian parameter and secondly by considering the relaxation times of the intermediate scattering function. From these measurements we develop a scaling relation in order to directly determine the combinations of wavevectors and times at which the non-Gaussian behavior is seen. This allows for the clear identification of the hydrodynamic regimes and thus provides new insight into the crossover between long- and short-time self-diffusion.