Flexible confinement leads to multiple relaxation regimes in glassy colloidal liquids

Computer simulations and mode-coupling theory of glass-forming confined hard-sphere fluids

We present mode-coupling theory (MCT) results for densely packed hard-sphere fluids confined between two parallel walls and compare them quantitatively to computer simulations. The numerical solution of MCT is calculated using the full system of matrix-valued integro-differential equations. We investigate several dynamical properties of supercooled liquids including scattering functions, frequency-dependent susceptibilities, and mean-square displacements. Close to the glass transition, we find quantitative agreement between the coherent scattering function predicted from theory and that evaluated from simulations, which enables us to make quantitative statements on caging and relaxation dynamics of the confined hard-sphere fluid.

The rheology of confined colloidal hard disks

Colloids may be treated as “big atoms” so that they are good models for atomic and molecular systems. Colloidal hard disks are, therefore, good models for 2d materials, and although their phase behavior is well characterized, rheology has received relatively little attention. Here, we exploit a novel, particle-resolved, experimental setup and complementary computer simulations to measure the shear rheology of quasi-hard-disk colloids in extreme confinement. In particular, we confine quasi-2d hard disks in a circular “corral” comprised of 27 particles held in optical traps. Confinement and shear suppress hexagonal ordering that would occur in the bulk and create a layered fluid. We measure the rheology of our system by balancing drag and driving forces on each layer. Given the extreme confinement, it is remarkable that our system exhibits rheological behavior very similar to unconfined 2d and 3d hard particle systems, characterized by a dynamic yield stress and shear-thinning of comparable magnitude. By quantifying particle motion perpendicular to shear, we show that particles become more tightly confined to their layers with no concomitant increase in density upon increasing the shear rate. Shear thinning is, therefore, a consequence of a reduction in dissipation due to weakening in interactions between layers as the shear rate increases. We reproduce our experiments with Brownian dynamics simulations with Hydrodynamic Interactions (HI) included at the level of the Rotne–Prager tensor. That the inclusion of HI is necessary to reproduce our experiments is evidence of their importance in transmission of momentum through the system.

Isomorphs in nanoconfined liquids

We study in this paper the possible existence of Roskilde-simple liquids and their isomorphs in a rough-wall nanoconfinement. Isomorphs are curves in the thermodynamic phase diagram along which structure and dynamics are invariant in suitable nondimensionalized units. Two model liquids using molecular dynamics computer simulations are considered: the single-component Lennard-Jones (LJ) liquid and the Kob-Andersen binary LJ mixture, both of which in the bulk phases are known to have good isomorphs. Nanoconfinement is implemented by adopting a slit-pore geometry with fcc crystalline walls; this implies inhomogenous density profiles both parallel and perpendicular to the confining walls. Despite this fact and consistent with an earlier study [Ingebrigtsen et al., Phys. Rev. Lett., 2013, 111, 235901] we find that these two nanoconfined liquids have isomorphs to a good approximation. More specifically, we show good invariance along the isomorphs of inhomogenous density profiles, mean-square displacements, and higher-order structures probed using the topological cluster classification algorithm. Our study thus provides an alternative framework for understanding nanoconfined liquids.

Ordering and Dynamics of Interacting Colloidal Particles under Soft Confinement

Confinement can induce substantial changes in the physical properties of macromolecules in suspension. Soft confinement is a particular class of restriction where the boundaries that constraint the particles in a region of the space are not well-defined. This scenario leads to a broader structural and dynamical behavior than observed in systems enclosed between rigid walls. In this contribution, we study the ordering and diffusive properties of a two-dimensional colloidal model system subjected to a one-dimensional parabolic trap. Increasing the trap strength makes it possible to go through weak to strong confinement, allowing a dimensional transition from two- to one-dimension. The non-monotonic response of the static and dynamical properties to the gradual dimensionality change affects the system phase behavior. We find that the particle dynamics are connected to the structural transitions induced by the parabolic trap. In particular, at low and intermediate confinement regimes, complex structural and dynamical scenarios arise, where the softness of the external potential induces melting and freezing, resulting in faster and slower particle diffusion, respectively. Besides, at strong confinements, colloids move basically along one direction, and the whole system behaves structurally and dynamically similar to a one-dimensional colloidal system.

Tagged-particle dynamics in confined colloidal liquids

We present numerical results for the tagged-particle dynamics by solving the mode-coupling theory in confined geometry for colloidal liquids (cMCT). We show that neither the microscopic dynamics nor the type of intermediate scattering function qualitatively changes the asymptotic dynamics in vicinity of the glass transition. In particular, we find similar characteristics of confinement in the low-frequency susceptibility spectrum which we interpret as footprints of parallel relaxation. We derive predictions for the localization length and the scaling of the diffusion coefficient in the supercooled regime and discover a pronounced nonmonotonic dependence on the confinement length. For dilute liquids in the hydrodynamic limit we calculate an analytical expression for the intermediate scattering functions, which is in perfect agreement with event-driven Brownian dynamics simulations. From this, we derive an expression for persistent anticorrelations in the velocity autocorrelation function (VACF) for confined motion. Using numerical results of the cMCT equations for the VACF we also identify a crossover between different scalings corresponding to a transition from unconfined to confined behavior.

Dynamical properties of densely packed confined hard-sphere fluids

Numerical solutions of the mode-coupling theory (MCT) equations for a hard-sphere fluid confined between two parallel hard walls are elaborated. The governing equations feature multiple parallel relaxation channels which significantly complicate their numerical integration. We investigate the intermediate scattering functions and the susceptibility spectra close to structural arrest and compare to an asymptotic analysis of the MCT equations. We corroborate that the data converge in the β-scaling regime to two asymptotic power laws, viz. the critical decay and the von Schweidler law. The numerical results reveal a nonmonotonic dependence of the power-law exponents on the slab width and a nontrivial kink in the low-frequency susceptibility spectra. We also find qualitative agreement of these theoretical results to event-driven molecular dynamics simulations of polydisperse hard-sphere systems. In particular, the nontrivial dependence of the dynamical properties on the slab width is well reproduced.

Understanding, predicting, and tuning the fragility of vitrimeric polymers

Fragility is an empirical property that describes how abruptly a glass-forming material solidifies upon supercooling. The degree of fragility carries important implications for the functionality and processability of a material, as well as for our fundamental understanding of the glass transition. However, the microstructural properties underlying fragility still remain poorly understood. Here, we explain the microstructure-fragility link in vitrimeric networks, a novel type of high-performance polymers with unique bond-swapping functionality and unusual glass-forming behavior. Our results are gained from coarse-grained computer simulations and first-principles mode-coupling theory (MCT) of star-polymer vitrimers. We first demonstrate that the vitrimer fragility can be tuned over an unprecedentedly broad range, from fragile to strong and even superstrong behavior, by decreasing the bulk density. Remarkably, this entire phenomenology can be reproduced by microscopic MCT, thus challenging the conventional belief that existing first-principles theories cannot account for nonfragile behaviors. Our MCT analysis allows us to rationally identify the microstructural origin of the fragile-to-superstrong crossover, which is rooted in the sensitivity of the static structure factor to temperature variations. On the molecular scale, this behavior stems from a change in dominant length scales, switching from repulsive excluded-volume interactions to intrachain attractions as the vitrimer density decreases. Finally, we develop a simplified schematic MCT model which corroborates our microscopically founded conclusions and which unites our findings with earlier MCT studies. Our work sheds additional light on the elusive structure-fragility link in glass-forming matter and provides a first-principles-based platform for designing amorphous materials with an on-demand dynamic response.

A Journey Through the Landscapes of Small Particles in Binary Colloidal Assemblies: Unveiling Structural Transitions from Isolated Particles to Clusters upon Variation in Composition

Two-dimensional (2D) amorphous binary colloidal assemblies composed of particles of two different sizes are characterized by the loss of hexagonal close-packing for larger particles, occurring when the size ratio between small (S) and large (L) particles dSdL exceeds a certain threshold value. For moderately low particle number ratios NSNL large particles still retain a denser arrangement with transitions from hexagonal symmetry to the coexistence of different types of symmetries as NSNL progressively departs from 0 to higher values. On the other hand, small particles reveal sparser arrangements: shape identification and quantification of structural transitions in small particle arrangements appear particularly challenging. In this article, we investigate their shapes and transitions for amorphous binary colloidal particles assembled at the air/water interface. For the quantitative characterization of the evolution in particle arrangements for NSNL variable between 0.5 and 2, we develop an innovative procedure for morphological analysis, combining Minkowski functionals, Voronoi diagrams and ad hoc techniques to recognize and classify specific features. Such a powerful approach has revealed a wide variety of landscapes featuring isolated particles, dimers, chains, small clusters evolving with the colloidal suspension composition. Our method can be applied to the analysis of spatial configurations of sparse colloidal patterns obtained in different conditions.

Magic number colloidal clusters as minimum free energy structures

Clusters in systems as diverse as metal atoms, virus proteins, noble gases, and nucleons have properties that depend sensitively on the number of constituent particles. Certain numbers are termed ‘magic’ because they grant the system with closed shells and exceptional stability. To this point, magic number clusters have been exclusively found with attractive interactions as present between atoms. Here we show that magic number clusters exist in a confined soft matter system with negligible interactions. Colloidal particles in an emulsion droplet spontaneously organize into a series of clusters with precisely defined shell structures. Crucially, free energy calculations demonstrate that colloidal clusters with magic numbers possess higher thermodynamic stability than those off magic numbers. A complex kinetic pathway is responsible for the efficiency of this system in finding its minimum free energy configuration. Targeting similar magic number states is a strategy towards unique configurations in finite self-organizing systems across the scales.

Magic number cluster with closed shells and increased stability often result from potential energy minimization between attractive atoms or particles. Here, Wang et al. show that such magic number clusters can also result from entropy maximization in colloidal systems with negligible interactions.

The race to the bottom: approaching the ideal glass?

Key to resolving the scientific challenge of the glass transition is to understand the origin of the massive increase in viscosity of liquids cooled below their melting temperature (avoiding crystallisation). A number of competing and often mutually exclusive theoretical approaches have been advanced to describe this phenomenon. Some posit a bona fide thermodynamic phase to an 'ideal glass', an amorphous state with exceptionally low entropy. Other approaches are built around the concept of the glass transition as a primarily dynamic phenomenon. These fundamentally different interpretations give equally good descriptions of the data available, so it is hard to determine which-if any-is correct. Recently however this situation has begun to change. A consensus has emerged that one powerful means to resolve this longstanding question is to approach the putative thermodynamic transition sufficiently closely, and a number of techniques have emerged to meet this challenge. Here we review the results of some of these new techniques and discuss the implications for the existence-or otherwise-of the thermodynamic transition to an ideal glass.

Nonergodicity parameters of confined hard-sphere glasses

Within a recently developed mode-coupling theory for fluids confined to a slit we elaborate numerical results for the long-time limits of suitably generalized intermediate scattering functions. The theory requires as input the density profile perpendicular to the plates, which we obtain from density functional theory within the fundamental-measure framework, as well as symmetry-adapted static structure factors, which can be calculated relying on the inhomogeneous Percus-Yevick closure. Our calculations for the nonergodicity parameters for both the collective as well as for the self motion are in qualitative agreement with our extensive event-driven molecular dynamics simulations for the intermediate scattering functions for slightly polydisperse hard-sphere systems at high packing fraction. We show that the variation of the nonergodicity parameters as a function of the wavenumber correlates with the in-plane static structure factors, while subtle effects become apparent in the structure factors and relaxation times of higher mode indices. A criterion to predict the multiple reentrant from the variation of the in-plane static structure is presented.