Experimental observation of structural crossover in binary mixtures of colloidal hard spheres

Nonadditive Interactions Unlock Small-Particle Mobility in Binary Colloidal Monolayers

We examine the organization and dynamics of binary colloidal monolayers composed of micron-scale silica particles interspersed with smaller-diameter silica particles that serve as minority component impurities. These binary monolayers are prepared at the surface of ionic liquid droplets over a range of size ratios (σ = 0.16-0.66) and are studied with low-dose minimally perturbative scanning electron microscopy (SEM). The high resolution of SEM imaging provides direct tracking of all particle coordinates over time, enabling a complete description of the microscopic state. In these bidisperse size mixtures, particle interactions are nonadditive because interfacial pinning to the droplet surface causes the equators of differently sized particles to lie in separate planes. By varying the size ratio, we control the extent of nonadditivity in order to achieve phase behavior inaccessible to additive 2D systems. Across the range of size ratios, we tune the system from a mobile small-particle phase (σ < 0.24) to an interstitial solid (0.24 < σ < 0.33) and furthermore to a disordered glass (σ > 0.33). These distinct phase regimes are classified through measurements of hexagonal ordering of the large-particle host lattice and the lattice's capacity for small-particle transport. Altogether, we explain these structural and dynamic trends by considering the combined influence of interparticle interactions and the colloidal packing geometry. Our measurements are reproduced in molecular dynamics simulations of 2D nonadditive disks, suggesting an efficient method for describing confined systems with reduced dimensionality representations.

Structural properties of additive binary hard-sphere mixtures. II. Asymptotic behavior and structural crossovers

The structural properties of additive binary hard-sphere mixtures are addressed as a follow-up of a previous paper [S. Pieprzyk et al., Phys. Rev. E 101, 012117 (2020)]2470-004510.1103/PhysRevE.101.012117. The so-called rational-function approximation method and an approach combining accurate molecular dynamics simulation data, the pole structure representation of the total correlation functions, and the Ornstein-Zernike equation are considered. The density, composition, and size-ratio dependencies of the leading poles of the Fourier transforms of the total correlation functions h_{ij}(r) of such mixtures are presented, those poles accounting for the asymptotic decay of h_{ij}(r) for large r. Structural crossovers, in which the asymptotic wavelength of the oscillations of the total correlation functions changes discontinuously, are investigated. The behavior of the structural crossover lines as the size ratio and densities of the two species are changed is also discussed.

Tetrahedrality Dictates Dynamics in Hard Sphere Mixtures

The link between local structure and dynamical slowdown in glassy fluids has been the focus of intense debate for the better part of a century. Nonetheless, a simple method to predict the dynamical behavior of a fluid purely from its local structural features is still missing. Here, we demonstrate that the diffusivity of perhaps the most fundamental family of glass formers-hard sphere mixtures-can be accurately predicted based on just the packing fraction and a simple order parameter measuring the tetrahedrality of the local structure. Essentially, we show that the number of tetrahedral clusters in a hard sphere mixture is directly linked to its global diffusivity. Moreover, the same order parameter is capable of locally pinpointing particles in the system with high and low mobility. We attribute the power of the local tetrahedrality for predicting local and global dynamics to the high stability of tetrahedral clusters, the most fundamental building and densest-packing building blocks for a disordered fluid.

Thermodynamic behavior of binary mixtures of hard spheres: Semianalytical solutions on a Husimi lattice built with cubes

We study binary mixtures of hard particles, which exclude up to their kth nearest neighbors (kNN) on the simple cubic lattice and have activities z_{k}. In the first model analyzed, point particles (0NN) are mixed with 1NN ones. The grand-canonical solution of this model on a Husimi lattice built with cubes unveils a phase diagram with a fluid and a solid phase separated by a continuous and a discontinuous transition line which meet at a tricritical point. A density anomaly, characterized by minima in isobaric curves of the total density of particles against z_{0} (or z_{1}), is also observed in this system. Overall, this scenario is identical to the one previously found for this model when defined on the square lattice. The second model investigated consists of the mixture of 1NN particles with 2NN ones. In this case, a very rich phase behavior is found in its Husimi lattice solution, with two solid phases-one associated with the ordering of 1NN particles (S1) and the other with the ordering of 2NN ones (S2)-beyond the fluid (F) phase. While the transitions between F-S2 and S1-S2 phases are always discontinuous, the F-S1 transition is continuous (discontinuous) for small (large) z_{2}. The critical and coexistence F-S1 lines meet at a tricritical point. Moreover, the coexistence F-S1,F-S2, and S1-S2 lines meet at a triple point. Density anomalies are absent in this case.

Active apolar doping determines routes to colloidal clusters and gels

Collections of polar active particles have been unable to form stable and long-living structures due to the presence of self-propulsion. We solve this timely issue by introducing the concept of “active doping” and show that a few light-activated apolar, i.e., non–self-propelling, units can be used to rapidly trigger the formation of solid clusters and gels composed of passive colloidal particles. Our active doping can be used to assemble disparate microscopic objects, including synthetic or biological ones, paving the way toward the extension of fundamental concepts of gel and glass formation to active out-of-equilibrium systems.

Collections of interacting active particles, self-propelling or not, have shown remarkable phenomena including the emergence of dynamic patterns across different length scales, from animal groups to vibrated grains, microtubules, bacteria, and chemical- or field-driven colloids. Burgeoning experimental and simulation activities are now exploring the possibility of realizing solid and stable structures from passive elements that are assembled by a few active dopants. Here we show that such an elusive task may be accomplished by using a small amount of apolar dopants, namely synthetic active but not self-propelling units. We use blue light to rapidly assemble 2D colloidal clusters and gels via nonequilibrium diffusiophoresis, where microscopic hematite dockers form long-living interstitial bonds that strongly glue passive silica microspheres. By varying the relative fraction of doping, we uncover a rich phase diagram including ordered and disordered clusters, space-filling gels, and bicontinuous structures formed by filamentary dockers percolating through a solid network of silica spheres. We characterize the slow relaxation and dynamic arrest of the different phases via correlation and scattering functions. Our findings provide a pathway toward the rapid engineering of mesoscopic gels and clusters via active colloidal doping.

Screening Lengths in Ionic Fluids

The decay of correlations in ionic fluids is a classical problem in soft matter physics that underpins applications ranging from controlling colloidal self-assembly to batteries and supercapacitors. The conventional wisdom, based on analyzing a solvent-free electrolyte model, suggests that all correlation functions between species decay with a common decay length in the asymptotic far field limit. Nonetheless, a solvent is present in many electrolyte systems. We show using an analytical theory and molecular dynamics simulations that multiple decay lengths can coexist in the asymptotic limit as well as at intermediate distances once a hard sphere solvent is considered. Our analysis provides an explanation for the recently observed discontinuous change in the structural force across a thin film of ionic liquid-solvent mixtures as the composition is varied, as well as reframes recent debates in the literature about the screening length in concentrated electrolytes.

Structural crossover in a model fluid exhibiting two length scales: Repercussions for quasicrystal formation

We investigate the liquid state structure of the two-dimensional model introduced by Barkan et al. [Phys. Rev. Lett. 113, 098304 (2014)10.1103/PhysRevLett.113.098304], which exhibits quasicrystalline and other unusual solid phases, focusing on the radial distribution function g(r) and its asymptotic decay r→∞. For this particular model system, we find that as the density is increased there is a structural crossover from damped oscillatory asymptotic decay with one wavelength to damped oscillatory asymptotic decay with another distinct wavelength. The ratio of these wavelengths is ≈1.932. Following the locus in the phase diagram of this structural crossover leads directly to the region where quasicrystals are found. We argue that identifying and following such a crossover line in the phase diagram towards higher densities where the solid phase(s) occur is a good strategy for finding quasicrystals in a wide variety of systems. We also show how the pole analysis of the asymptotic decay of equilibrium fluid correlations is intimately connected with the nonequilibrium growth or decay of small-amplitude density fluctuations in a bulk fluid.

Interfacial structure and structural forces in mixtures of ionic liquid with a polar solvent

Oscillatory and monotonic decay in mixtures of salt and solvent at interfaces with varying surface charge.

Switching the Structural Force in Ionic Liquid-Solvent Mixtures by Varying Composition

The structure and interactions in electrolytes at high concentration have implications from energy storage to biomolecular interactions. However, many experimental observations are yet to be explained in these mixtures, which are far beyond the regime of validity of mean-field models. Here, we study the structural forces in a mixture of ionic liquid and solvent that is miscible in all proportions at room temperature. Using the surface force balance to measure the force between macroscopic smooth surfaces across the liquid mixtures, we uncover an abrupt increase in the wavelength above a threshold ion concentration. Below the threshold concentration, the wavelength is determined by the size of the solvent molecule, whereas above the threshold, it is the diameter of a cation-anion pair that determines the wavelength.

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.

Equilibrium Grain Boundary Segregation and Clustering of Impurities in Colloidal Polycrystalline Monolayers

We investigate the segregation of impurities to grain boundaries in colloidal polycrystalline monolayers using video microscopy. A model colloidal alloy is prepared by embedding large spherical impurities in a polycrystalline monolayer of small host colloidal hard spheres, which stops grain growth at a finite grain size. The size ratio between the impurities and the host particles determines whether they behave as interstitial or substitutional impurities in the bulk crystal, akin to those in real alloys. We find that the partitioning of impurities between the grains and the grain boundaries is in very good agreement with the Langmuir-McLean adsorption model for equilibrium grain boundary segregation. This enables the direct measurement of the free energy of adsorption for the two types of impurities. Near saturation, we characterize the spatial distribution of the adsorbed impurities and find that it strongly depends on their interstitial or substitutional nature. This is because the relative importance of clustering and mixing due to nonadditivity is determined by geometrical constraints imposed by the crystalline host lattice.

Cavity averages for hard spheres in the presence of polydispersity and incomplete data

We develop a cavity-based method which allows to extract thermodynamic properties from position information in hard-sphere/disk systems. So far, there are available-volume and free-volume methods. We add a third one, which we call available volume after take-out, and which is shown to be mathematically equivalent to the others. In applications, where data sets are finite, all three methods show limitations, and they do this in different parameter ranges. We illustrate the principal equivalence and the limitations on data from molecular dynamics: In particular, we test robustness against missing data. We have in mind experimental limitations where there is a small polydispersity, say 4% in the particle radii, but individual radii cannot be determined. We observe that, depending on the used method, the errors in such a situation are easily 100% for the pressure and 10kT for the chemical potentials. Our work is meant as guideline to the experimentalists for choosing the right one of the three methods, in order to keep the outcome of experimental data analysis meaningful.

Understanding shape entropy through local dense packing

Many natural systems are structured by the ordering of repeated, distinct shapes. Understanding how this happens is difficult because shape affects structure in two ways. One is how the shape of a cell or nanoparticle, for example, affects its surface, chemical, or other intrinsic properties. The other is an emergent, entropic effect that arises from the geometry of the shape itself, which we term “shape entropy,” and is not well understood. In this paper, we determine how shape entropy affects structure. We quantify the mechanism and determine when shape entropy competes with intrinsic shape effects. Our results show that in a wide class of systems, shape affects bulk structure because crowded particles optimize their local packing.

Understanding shape entropy through local dense packing

Entropy drives the phase behavior of colloids ranging from dense suspensions of hard spheres or rods to dilute suspensions of hard spheres and depletants. Entropic ordering of anisotropic shapes into complex crystals, liquid crystals, and even quasicrystals was demonstrated recently in computer simulations and experiments. The ordering of shapes appears to arise from the emergence of directional entropic forces (DEFs) that align neighboring particles, but these forces have been neither rigorously defined nor quantified in generic systems. Here, we show quantitatively that shape drives the phase behavior of systems of anisotropic particles upon crowding through DEFs. We define DEFs in generic systems and compute them for several hard particle systems. We show they are on the order of a few times the thermal energy ([Formula: see text]) at the onset of ordering, placing DEFs on par with traditional depletion, van der Waals, and other intrinsic interactions. In experimental systems with these other interactions, we provide direct quantitative evidence that entropic effects of shape also contribute to self-assembly. We use DEFs to draw a distinction between self-assembly and packing behavior. We show that the mechanism that generates directional entropic forces is the maximization of entropy by optimizing local particle packing. We show that this mechanism occurs in a wide class of systems and we treat, in a unified way, the entropy-driven phase behavior of arbitrary shapes, incorporating the well-known works of Kirkwood, Onsager, and Asakura and Oosawa.

Magnetic assembly and annealing of colloidal lattices and superlattices

The ability to assemble mesoscopic colloidal lattices above a surface is important for fundamental studies related with nucleation and crystallization but also for a variety of technological applications in photonics and microengineering. Current techniques based on particle sedimentation above a lithographic template are limited by a slow deposition process and by the use of static templates, which make difficult to implement fast annealing procedures. Here it is demonstrated a method to realize and anneal a series of colloidal lattices displaying triangular, honeycomb, or kagome-like symmetry above a structure magnetic substrate. By using a binary mixture of particles, superlattices can be realized increasing further the variety and complexity of the colloidal patterns which can be produced.

Binary non-additive hard sphere mixtures: fluid demixing, asymptotic decay of correlations and free fluid interfaces

Using a fundamental measure density functional theory we investigate both bulk and inhomogeneous systems of the binary non-additive hard sphere model. For sufficiently large (positive) non-additivity the mixture phase separates into two fluid phases with different compositions. We calculate bulk fluid-fluid coexistence curves for a range of size ratios and non-additivity parameters and find that they compare well to simulation results from the literature. Using the Ornstein-Zernike equation, we investigate the asymptotic, [Formula: see text], decay of the partial pair correlation functions, g(ij)(r). At low densities a structural crossover occurs in the asymptotic decay between two different damped oscillatory modes with different wavelengths corresponding to the two intra-species hard-core diameters. On approaching the fluid-fluid critical point there is a Fisher-Widom crossover from exponentially damped oscillatory to monotonic asymptotic decay. Using the density functional we calculate the density profiles for the planar free fluid-fluid interface between coexisting fluid phases. We show that the type of asymptotic decay of g(ij)(r) not only determines the asymptotic decay of the interface profiles, but is also relevant for intermediate and even short-ranged behaviour. We also determine the surface tension of the free fluid interface, finding that it increases with non-additivity, and that on approaching the critical point mean-field scaling holds.

Lateral resolution enhancement in confocal microscopy by vectorial aperture engineering

This article reports the design and implementation of a lateral resolution-enhancement technique in confocal microscopy that can work, in principle, either in the reflection mode or in the fluorescence mode. Taking the difference between two images corresponding to two different vectorially (involving amplitude, phase, and polarization of light) engineered illumination pupils or apertures of a confocal microscope, high spatial frequency contents in the resultant image can be significantly enhanced. This can be realized by incorporating an extra vectorial beam-forming element into the illumination beam path of a conventional confocal microscope. The method of the proposed technique has been explained by giving it an analytical treatment supported by numerical simulation results. The technique has been implemented in a reflection mode confocal microscope and results obtained are presented.

Controlled structuring of binary hard-disk mixtures via a periodic, external potential

Ordering phenomena on surfaces or in monolayers can be successfully studied by model systems as binary hard-disk mixtures, the influence of a substrate being modeled by an external potential. For the field-free case the thermodynamic stability of space-filling lattice structures for binary hard-disk mixtures is studied by Monte Carlo computer simulations. As these structures prove to be thermodynamically stable only in high pressure environments, the phase behavior of an equimolar binary mixture with a diameter ratio of sigma_{B}/sigma_{A}=0.414 exposed to an external, one-dimensional, periodic potential is analyzed in detail. The underlying ordering mechanisms and the resulting order differ considerably, depending on which components of the mixture interact with the external potential. The simulations show that slight deviations in the concentration of large particles x_{A} or the diameter ratio sigma_{B}/sigma_{A} have no impact on the occurrence of the various field-induced phenomena as long as the mixture stays in the relevant regime of the packing fraction eta . Furthermore the importance of the commensurability of the external potential to the S1(AB) square lattice for the occurrence of the induced ordering is discussed.

Structure of binary colloidal systems confined in a quasi-one-dimensional channel

The structural properties of a binary colloidal quasi-one-dimensional system confined in a narrow channel are investigated through modified Monte Carlo simulations. Two species of particles with different magnetic moment interact through a repulsive dipole-dipole force are confined in a quasi-one-dimensional channel. The impact of three decisive parameters (the density of particles, the magnetic-moment ratio, and the fraction between the two species) on the transition from disordered phase to crystal-like phases and the transitions among the different mixed phases are summarized in a phase diagram.

Test particle limit for the pair structure of quenched-annealed fluid mixtures

A density functional route to the pair structure of quenched-annealed fluid mixtures is presented. The bulk two-body partial pair correlation functions of the mixture are identified with the one-body density distributions in an external potential that models a fixed test particle. Quenched-annealed (replica) density functional theory is used to calculate the inhomogeneous one-body density distributions. A closed theory is obtained by using an exact sum rule that equates two different expressions for the cross pair correlation function between unlike species. Results for binary quenched-annealed hard sphere mixtures are found to agree well with computer simulation data, improving over results from the replica Ornstein-Zernike equations using the direct correlation functions, obtained as second functional derivatives of the quenched-annealed excess free energy functional, as input. The proposed framework allows for the independent determination of the blocked part of the fluid-fluid partial pair correlation function.

Equilibrating nanoparticle monolayers using wetting films

Monolayers of bimodal gold nanoparticles on silicon are investigated by a combination of microscopy (dry monolayers) and x-ray diffraction (dry and wet monolayers). In the presence of an excess of small particles, the nanoscale packing structure closely resembles the small-particle-rich scenario of the structural crossover transition that has been predicted and also observed with micron-scale hard-sphere colloids. Structural morphology is monitored in situ during monolayer dissolution and reassembly within the thin liquid wetting film. This approach allows investigation of size and solvent effects on nanoparticles in quasi-two-dimensional confinement.

Disordering transitions and peak effect in polydisperse particle systems

We show numerically that in a binary system of Yukawa particles, a dispersity-driven disordering transition occurs. In the presence of quenched disorder this disordering transition coincides with a marked increase in the depinning threshold, known as a peak effect. We find that the addition of poorly pinned particles can increase the overall pinning in the sample by increasing the amount of topological disorder present. If the quenched disorder is strong enough to create a significant amount of topological disorder in the monodisperse system, addition of a poorly pinned species generates further disorder but does not produce a peak in the depinning force. Our results indicate that for binary mixtures, optimal pinning occurs for topological defect fraction densities from 0.2 to 0.25. For defect densities below this range, the system retains orientational order. We determine the effect of the pinning density, strength, and radius on the depinning peak and find that the peak effect is more pronounced in weakly pinning systems.

Entropy versus energy: the phase behavior of a hard-disk mixture in a periodic external potential

The phase behavior of a 50% binary hard-disk mixture with diameter ratio sigmaB/sigmaA=0.414 , which is exposed to a one-dimensional periodic potential, is examined via Monte Carlo simulations. We find an induced structural crossover in the modulated liquid. At higher densities, depending on the strength of the external potential, the system exhibits a tunable demixing transition, followed by fluid-solid coexistence of an equimolar mixture with the S1(AB) square lattice. We find a decoupled melting of the sublattices of the S1(AB) lattice. The melting of the small-component sublattice perpendicular to the external potential minima leads to fissuring in the large-component sublattice.

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.