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

Phase diagram of the hard-sphere potential model in three and four dimensions using a pseudo-hard-sphere potential

The hard-sphere potential has become a cornerstone in the study of both molecular and complex fluids. Despite its mathematical simplicity, its implementation in fixed time step molecular simulations remains a formidable challenge due to the discontinuity at contact. To avoid the issues associated with the ill-defined force at contact, a continuous potential has recently been proposed-here referred to as the pseudo-hard-sphere potential (pHS) [Báez et al., J. Chem, Phys. 149, 164907 (2018)]. This potential is constructed to match the second virial coefficient of the hard-sphere potential and is expected to mimic its thermodynamic properties. However, this hypothesis has only been partially validated within the fluid region of the phase diagram for hard-sphere dispersions in two and three dimensions. In this contribution, we examine the ability of the continuous pHS potential to reproduce the equation of state of a hard-sphere fluid, not only in the fluid phase but also across the fluid-solid coexistence region. Our focus is primarily on the phase diagram of hard-sphere systems in three and four dimensions; however, we also report on the feasibility of the pHS to reproduce the long time dynamics of a three-dimensional colloidal dispersion. We compare the thermodynamic properties obtained from Brownian dynamics simulations of the pHS potential with those derived from refined event-driven simulations of the corresponding hard-sphere potential. Furthermore, we provide a comparative analysis with theoretical equations of state based on both mean-field and integral equation approximations.

Coupled dynamics in binary mixtures of model colloidal Yukawa systems

The dynamical behavior of binary mixtures consisting of highly charged colloidal particles is studied by means of Brownian dynamics simulations. We investigate differently sized, but identically charged particles with nearly identical interactions between all species in highly dilute suspensions. Different short-time self-diffusion coefficients induce, mediated by electrostatic interactions, a coupling of both self and collective dynamics of differently sized particles: the long-time self-diffusion coefficients of a larger species are increased by the presence of a more mobile, smaller species and vice versa. Similar coupling effects are observed in collective dynamics where both time constant and functional form of intermediate scattering functions' initial decay are influenced by the presence of a differently sized species. We provide a systematic analysis of coupling effects in dependence on the ratio of sizes, number densities, and the strength of electrostatic interactions.

Rescaled mode-coupling scheme for the quantitative description of experimentally observed colloid dynamics

We describe experimentally observed collective dynamics in colloidal suspensions of model hard-sphere particles using a modified mode coupling theory (MCT). This rescaled MCT is capable of describing quantitatively the wave-vector and time-dependent diffusion in these systems. Intermediate scattering functions of liquidlike structured dispersions are determined by means of static and dynamic light-scattering experiments. The structure and short-time dynamics of the systems can be described quantitatively employing a multicomponent Percus-Yevick ansatz for the partial structure factors and an effective, one-component description of hydrodynamic interactions based on the semianalytical δγ expansion. Combined with a recently proposed empirical modification of MCT in which memory functions are calculated using effective structure factors at rescaled number densities, the scheme is able to model the collective dynamics over the entire accessible time and wave-vector range and predicts the volume-fraction-dependence of long-time self-diffusion coefficients and the zero-shear viscosity quantitatively. This highlights the potential of MCT as a practical tool for the quantitative analysis and prediction of experimental observations.

Quasicrystal of Binary Hard Spheres on a Plane Stabilized by Configurational Entropy

Because of their aperiodic nature, quasicrystals are one of the least understood phases in statistical physics. One significant complication they present in comparison to their periodic counterparts is the fact that any quasicrystal can be realized as an exponentially large number of different tilings, resulting in a significant contribution to the quasicrystal entropy. Here, we use free-energy calculations to demonstrate that it is this configurational entropy which stabilizes a dodecagonal quasicrystal in a binary mixture of hard spheres on a plane. Our calculations also allow us to quantitatively confirm that in this system all tiling realizations are essentially equally likely, with free-energy differences less than 0.0001k_{B}T per particle-an observation that could be related to the observation of only random tilings in soft-matter quasicrystals. Owing to the simplicity of the model and its available counterparts in colloidal experiments, we believe that this system is an excellent candidate to achieve the long-awaited quasicrystal self-assembly on the micron scale.

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.

Modeling the structure and thermodynamics of multicomponent and polydisperse hard-sphere dispersions with continuous potentials

A model system of identical particles interacting via a hard-sphere potential is essential in condensed matter physics; it helps to understand in and out of equilibrium phenomena in complex fluids, such as colloidal dispersions. Yet, most of the fixed time-step algorithms to study the transport properties of those systems have drawbacks due to the mathematical nature of the interparticle potential. Because of this, mapping a hard-sphere potential onto a soft potential has been recently proposed [Báez et al., J. Chem. Phys. 149, 164907 (2018)]. More specifically, using the second virial coefficient criterion, one can set a route to estimate the parameters of the soft potential that accurately reproduces the thermodynamic properties of a monocomponent hard-sphere system. However, real colloidal dispersions are multicomponent or polydisperse, making it important to find an efficient way to extend the potential model for dealing with such kind of many-body systems. In this paper, we report on the extension and applicability of the second virial coefficient criterion to build a description that correctly captures the phenomenology of both multicomponent and polydisperse hard-sphere dispersions. To assess the accuracy of the continuous potentials, we compare the structure of soft polydisperse systems with their hard-core counterpart. We also contrast the structural and thermodynamic properties of soft binary mixtures with those obtained through mean-field approximations and the Ornstein-Zernike equation for the two-component hard-sphere dispersion.

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.

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.

Structure and dynamics of a two-dimensional colloid of liquid droplets

Droplet arrays in thin, freely suspended liquid-crystalline smectic A films can form two-dimensional (2D) colloids. The droplets interact repulsively, arranging locally in a more or less hexagonal arrangement with only short-range spatial and orientational correlations and local lattice cell parameters that depend on droplet size. In contrast to quasi-2D colloids described earlier, there is no 3D bulk liquid subphase that affects the hydrodynamics. Although the films are surrounded by air, the droplet dynamics are genuinely 2D, the mobility of each droplet in its six-neighbor cage being determined by the ratio of cage and droplet sizes, rather than by the droplet size as in quasi-2D colloids. These experimental observations are described well by Saffman's model of a diffusing particle in a finite 2D membrane. The experiments were performed in microgravity, on the International Space Station.

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.

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.