Depletion force in a bidisperse granular layer

Cooling of a granular gas mixture in microgravity

Granular gases are fascinating non-equilibrium systems with interesting features such as spontaneous clustering and non-Gaussian velocity distributions. Mixtures of different components represent a much more natural composition than monodisperse ensembles but attracted comparably little attention so far. We present the observation and characterization of a mixture of rod-like particles with different sizes and masses in a drop tower experiment. Kinetic energy decay rates during granular cooling and collision rates were determined and Haff's law for homogeneous granular cooling was confirmed. Thereby, energy equipartition between the mixture components and between individual degrees of freedom is violated. Heavier particles keep a slightly higher average kinetic energy than lighter ones. Experimental results are supported by numerical simulations.

Bidisperse Nanospheres Jammed on a Liquid Surface

Jammed packings of bidisperse nanospheres were assembled on a nonvolatile liquid surface and visualized to the single-particle scale by using an in situ scanning electron microscopy method. The PEGylated silica nanospheres, mixed at different number fractions and size ratios, had large enough in-plane mobilities prior to jamming to form uniform monolayers reproducibly. From the collected nanometer-resolution images, local order and degree of mixing were assessed by standard metrics. For equimolar mixtures, a large-to-small size ratio of about 1.5 minimized correlated metrics for local orientational and positional order, as previously predicted in simulations of bidisperse disk jamming. Despite monolayer uniformity, structural and depletion interactions caused spheres of a similar size to cluster, a feature evident at size ratios above 2. Uniform nanoparticle monolayers of high packing disorder are sought in many liquid interface technologies, and these experiments outlined key design principles, buttressing extensive theory/simulation literature on the topic.

Demixing and tetratic ordering in some binary mixtures of hard superellipses

We examine the fluid phase behavior of binary mixtures of hard superellipses using the scaled particle theory. The superellipse is a general two-dimensional convex object that can be tuned between the elliptical and rectangular shapes continuously at a given aspect ratio. We find that the shape of the particle affects strongly the stability of isotropic, nematic, and tetratic phases in the mixture even if the side lengths of both species are fixed. While the isotropic-isotropic demixing transition can be ruled out using the scaled particle theory, the first order isotropic-nematic and the nematic-nematic demixing transition can be stabilized with strong fractionation between the components. It is observed that the demixing tendency is strongest in small rectangle-large ellipse mixtures. Interestingly, it is possible to stabilize the tetratic order at lower densities in the mixture of hard squares and rectangles where the long rectangles form a nematic phase, while the squares stay in the tetratic order.

Small-scale demixing in confluent biological tissues

Surface tension governed by differential adhesion can drive fluid particle mixtures to sort into separate regions, i.e., demix. Does the same phenomenon occur in confluent biological tissues? We begin to answer this question for epithelial monolayers with a combination of theory via a vertex model and experiments on keratinocyte monolayers. Vertex models are distinct from particle models in that the interactions between the cells are shape-based, as opposed to distance-dependent. We investigate whether a disparity in cell shape or size alone is sufficient to drive demixing in bidisperse vertex model fluid mixtures. Surprisingly, we observe that both types of bidisperse systems robustly mix on large lengthscales. On the other hand, shape disparity generates slight demixing over a few cell diameters, a phenomenon we term micro-demixing. This result can be understood by examining the differential energy barriers for neighbor exchanges (T1 transitions). Experiments with mixtures of wild-type and E-cadherin-deficient keratinocytes on a substrate are consistent with the predicted phenomenon of micro-demixing, which biology may exploit to create subtle patterning. The robustness of mixing at large scales, however, suggests that despite some differences in cell shape and size, progenitor cells can readily mix throughout a developing tissue until acquiring means of recognizing cells of different types.

Flow-assisted droplet assembly in a 3D microfluidic channel

Self-assembly of soft matter, such as droplets or colloids, has become a promising scheme to engineer novel materials, model living matter, and explore non-equilibrium statistical mechanics. In this article, we present detailed numerical simulations of few non-Brownian droplets in various flow conditions, specifically, focusing on their self-assembly within a short distance in a three-dimensional (3D) microfluidic channel, cf. [Shen et al., Adv. Sci., 2016, 3(6), 1600012]. Contrary to quasi two-dimensional (q2D) systems, where dipolar interaction is the key mechanism for droplet rearrangement, droplets in 3D confinement produce much less disturbance to the underlying flow, thus experiencing weaker dipolar interactions. Using confined simple shear and Poiseuille flows as reference flows, we show that the droplet dynamics is mostly affected by the shear-induced cross-stream migration, which favors chain structures if the droplets are under an attractive depletion force. For more compact clusters, such as three droplets in a triangular shape, our results suggest that an inhomogeneous cross-sectional inflow profile is further required. Overall, the accelerated self-assembly of a small-size droplet cluster results from the combined effects of strong depletion forces, confinement-mediated shear alignments, and fine-tuned inflow conditions. The deterministic nature of the flow-assisted self-assembly implies the possibility of large throughputs, though calibration of all different effects to directly produce large droplet crystals is generally difficult.

Clustering-Induced Attraction in Granular Mixtures of Rods and Spheres

Depletion-induced aggregation of rods enhanced by clustering is observed to produce a novel model of attractive pairs of rods separated by a line of spheres in a quasi-2D, vertically-shaken, granular gas of rods and spheres. We show that the stability of these peculiar granular aggregates increases as a function of shaking intensity. Velocity distributions of spheres inside and outside of a pair of rods trapping a line of spheres show a clear suppression of the momentum acquired by the trapped spheres. The condensed phase formed between the rods is caused by a clustering instability of the trapped spheres, enhanced by a vertical guidance produced by the confining rods. The liberated area corresponding to direct excluded-volume pairs and indirect depletion-aggregated pairs is measured as a function of time. The stability of rod pairs mediated by spheres reveals an attraction comparable in strength to the one purely induced by depletion forces.

Nucleation, aggregation, annealing, and disintegration of granular clusters

The processes of nucleation, aggregation, annealing, and disintegration of clusters of non-Brownian paramagnetic beads in a vibrofluidized system are experimentally investigated. The interaction among the beads is induced by a magnetic seed composed of two dipoles allocated outside the container cell. We observe a clearly differentiated nucleation stage, whose evolution (nucleation time versus acceleration strength) follows a power law. Thereafter, the beads aggregate forming 2D disordered clusters around the nucleus. Both processes (nucleation and aggregation) are determined by the competition between magnetic forces and the drag produced by a thermal bath created by glass particles. Once the agglomerates reach a final state (shape and length), they are annealed by increasing and decreasing the granular temperature. We found that the fractal dimension and the lacunarity index clearly describe the structural variations of the clusters. Our discussion on this phenomenon is addressed, making a rough analogy with the glass transition in a super-cooled liquid. Finally, we study the disintegration of the clusters as a function of time and the density of the surrounding gas. The question is not if, but how they disintegrate upon removing the external field; we find that the disintegration follows an exponential decay.

Horizontal segregation of mono-layer granules coordinated by vertical motion

We experimentally investigate the segregation of a binary mixture of spherical beads confined between two horizontal vertically vibrating plates. The two kinds of beads are of equal diameter and mass but have different restitution coefficients. Segregation occurs in particular ranges of vibration amplitude and frequency. We find that the collisions between beads at an angle to the horizontal plane induce an effective horizontal repulsive force. When one or both bead types bounce up and down in synchronization, the effective repulsive force between the two types of beads is likely to be larger than that found within a single bead type, resulting in the mixture segregating. Non-horizontal collisions also play a role in stabilizing the segregation state by transferring the horizontal kinetic energy back into vertical motion.

Depletion forces drive polymer-like self-assembly in vibrofluidized granular materials

Ranging from nano- to granular-scales, control of particle assembly can be achieved by limiting the available free space, for example by increasing the concentration of particles ("crowding") or through their restriction to 2D environments. It is unclear, however, if self-assembly principles governing thermally-equilibrated molecules can also apply to mechanically-excited macroscopic particles in non-equilibrium steady-state. Here we show that low densities of vibrofluidized steel rods, when crowded by high densities of spheres and confined to quasi-2D planes, can self-assemble into linear polymer-like structures. Our 2D Monte Carlo simulations show similar finite sized aggregates in thermally equilibrated binary mixtures. Using theory and simulations, we demonstrate how depletion interactions create oriented "binding" forces between rigid rods to form these "living polymers." Unlike rod-sphere mixtures in 3D that can demonstrate well-defined equilibrium phases, our mixtures confined to 2D lack these transitions because lower dimensionality favors the formation of linear aggregates, thus suppressing a true phase transition. The qualitative and quantitative agreement between equilibrium and granular patterning for these mixtures suggests that entropy maximization is the determining driving force for bundling. Furthermore, this study uncovers a previously unknown patterning behavior at both the granular and nanoscales, and may provide insights into the role of crowding at interfaces in molecular assembly.

Self-organized patterns in collections of chains

Results are presented of an experimental investigation into patterned segregation in thin layers of poppy seeds and short lengths of metal chains subjected to vibration. Critical phenomena are uncovered and both continuous and discontinuous transitions are observed. A phase diagram for the behavior is mapped out and a tricritical point that separates hysteretic from continuous segregation is identified. Remarkable similarities are found between the observed behavior in this driven granular system and phase separation phenomena in mixtures where the dynamics of the constituent components are markedly different.

Segregation induced by phase synchronization in a bidisperse granular layer

We propose an alternative segregation mechanism where the species-dependent interactions are dynamically induced by the phase synchronization of beads. Based on this scenario, we report an alternative segregation among beads of different restitution coefficients by molecular dynamics simulations. Since the beads are of equal size and mass, this is not related to the Brazilian-nut effect, nor can it be explained by the depletion force. Instead, this phenomenon derives from the phase synchronization, a concept which helps us determine the criteria for segregation and the phase boundaries that agree excellently with the simulation results.

Third-order thermodynamic perturbation theory for effective potentials that model complex fluids

We have performed Monte Carlo simulations to obtain the thermodynamic properties of fluids with two kinds of hard-core plus attractive-tail or oscillatory potentials. One of them is the square-well potential with small well width. The other is a model potential with oscillatory and decaying tail. Both model potentials are suitable for modeling the effective potential arising in complex fluids and fluid mixtures with extremely-large-size asymmetry, as is the case of the solvent-induced depletion interactions in colloidal dispersions. For the former potential, the compressibility factor, the excess energy, the constant-volume excess heat capacity, and the chemical potential have been obtained. For the second model potential only the first two of these quantities have been obtained. The simulations cover the whole density range for the fluid phase and several temperatures. These simulation data have been used to test the performance of a third-order thermodynamic perturbation theory (TPT) recently developed by one of us [S. Zhou, Phys. Rev. E 74, 031119 (2006)] as compared with the well-known second-order TPT based on the macroscopic compressibility approximation due to Barker and Henderson. It is found that the first of these theories provides much better accuracy than the second one for all thermodynamic properties analyzed for the two effective potential models.