Glassy dynamics in asymmetric binary mixtures of hard spheres

Anisotropic remixing of a phase separated binary colloidal system with particles of different sizes in an external modulation

We explore the phase behavior of a binary colloidal system under external spatially periodic modulation. We perform Monte Carlo simulations on a binary mixture of big and small repulsive Lennard-Jones particles with a diameter ratio of 2:1. We characterize structure by isotropic and anisotropic pair correlation functions, cluster size distribution, bond angle distribution, order parameter, and specific heat. We observe the demixing of the species in the absence of external modulation. However, the mixing of the species gets enhanced with increasing potential strength along with the alignment of the particles transverse to the modulation. The de-mixing order parameter shows discontinuity with increasing modulation strength, characterizing a first order phase transition. The peak in specific heat increases linearly with the size of the system. We also look into the dynamical behavior of the system via computing Mean Square Displacement (MSD) along both parallel and perpendicular directions to the modulation. We observe a decrease in the diffusion coefficient for both types of particles as we increase the strength of the modulation.

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.

Cooperative activated hopping dynamics in binary glass-forming liquids: effects of the size ratio, composition, and interparticle interactions

Slow dynamics in supercooled and glassy liquids is an important research topic in soft matter physics. Compared to the traditionally focused one-component systems, glassy dynamics in mixture systems adds in a rich set of new complexities, which are fundamentally interesting and also relevant for many technological applications. In this paper, we apply the recently developed self-consistent cooperative hopping theory (SCCHT) to systematically investigate the effects of the size ratio, composition and interparticle interactions on the cooperative activated hopping dynamics of matrix (in larger size) and penetrant (in smaller size) particles in varied binary sphere mixture model systems, with a specific focus on ultrahigh mixture packing fractions that mimic the deeply supercooled glass transition conditions for molecular/polymeric mixture materials. Analysis shows that in these high activation barrier cases, the long-range elastic distortion associated with a matrix particle hopping over its cage confinement always generates an elastic barrier of a nonnegligible magnitude, although the ratio between the elastic barrier and local barrier contribution is sensitively dependent on all three mixture-specific system factors considered in this work. SCCHT predicts two general scenarios of penetrant-matrix cooperative activated hopping dynamics: matrix/penetrant co-hopping (regime 1) or the penetrant mean barrier hopping time shorter than that of the matrix (regime 2). Increasing the penetrant-to-matrix size ratio or the penetrant-matrix cross-attraction strength is found to universally enlarge the composition window of regime 1. Diverse dynamical properties characterising different aspects of the cooperative activated hopping process, including the penetrant and matrix transient localization lengths, penetrant and matrix hopping jump distances, different types of local and elastic activated barriers, and matrix long-time diffusivity, relaxation time and dynamic fragility are quantitatively studied against a wide range of variations over the three system factors. Of particular interest is the universal "anti-plasticization" phenomenon achievable for sufficiently strong cross-attractive interactions. The prospects this work opens for the exploration of a wide variety of polymer-based mixture materials are briefly discussed at the end.

Non-equilibrium view of the amorphous solidification of liquids with competing interactions

The interplay between short-range attractions and long-range repulsions (SALR) characterizes the so-called liquids with competing interactions, which are known to exhibit a variety of equilibrium and non-equilibrium phases. The theoretical description of the phenomenology associated with glassy or gel states in these systems has to take into account both the presence of thermodynamic instabilities (such as those defining the spinodal line and the so called λ line) and the limited capability to describe genuine non-equilibrium processes from first principles. Here, we report the first application of the non-equilibrium self-consistent generalized Langevin equation theory to the description of the dynamical arrest processes that occur in SALR systems after being instantaneously quenched into a state point in the regions of thermodynamic instability. The physical scenario predicted by this theory reveals an amazing interplay between the thermodynamically driven instabilities, favoring equilibrium macro- and micro-phase separation, and the kinetic arrest mechanisms, favoring non-equilibrium amorphous solidification of the liquid into an unexpected variety of glass and gel states.

Phase separation and dynamical arrest of protein solutions dominated by short-range attractions

The interplay of liquid-liquid phase separation (LLPS) and dynamical arrest can lead to the formation of gels and glasses, which is relevant for such diverse fields as condensed matter physics, materials science, food engineering, and the pharmaceutical industry. In this context, protein solutions exhibit remarkable equilibrium and non-equilibrium behaviors. In the regime where attractive and repulsive forces compete, it has been demonstrated, for example, that the location of the dynamical arrest line seems to be independent of ionic strength, so that the arrest lines at different ionic screening lengths overlap, in contrast to the LLPS coexistence curves, which strongly depend on the salt concentration. In this work, we show that the same phenomenology can also be observed when the electrostatic repulsions are largely screened, and the range and strength of the attractions are varied. In particular, using lysozyme in brine as a model system, the metastable gas-liquid binodal and the dynamical arrest line as well as the second virial coefficient have been determined for various solution conditions by cloud-point measurements, optical microscopy, centrifugation experiments, and light scattering. With the aim of understanding this new experimental phenomenology, we apply the non-equilibrium self-consistent generalized Langevin equation theory to a simple model system with only excluded volume plus short-range attractions, to study the dependence of the predicted arrest lines on the range of the attractive interaction. The theoretical predictions find a good qualitative agreement with experiments when the range of the attraction is not too small compared with the size of the protein.

Molecular dynamics simulations of binary sphere mixtures

Explicit simulations of fluid mixtures of highly size-dispersed particles are constrained by numerical challenges associated with identifying pair-interaction neighbors. Recent algorithmic developments have ameliorated these difficulties to an extent, permitting more efficient simulations of systems with many large and small particles of disperse sizes. We leverage these capabilities to perform molecular dynamics simulations of binary sphere mixtures with elastically stiff particles approaching the hard sphere limit and particle size ratios of up to 50, approaching the colloidal limit. The systems considered consist of 500 large particles and up to nearly 3.6×10^{6} small particles with total particle volume fractions up to 0.51. Our simulations confirm qualitative predictions for correlations between large particles previously obtained analytically and for simulations employing effective depletion interactions, but also reveal additional insights into the near-contact structure that result from the explicit treatment of the small particle solvent. No spontaneous crystal nucleation was observed during the simulations, suggesting that nucleation rates in the fluid-solid coexistence region are too small to observe crystal nucleation for feasible simulation system sizes and timescales.

Finite-size effects in the diffusion dynamics of a glass-forming binary mixture with large size ratio

Extensive molecular dynamics computer simulations of an equimolar, glass-forming AB mixture with a large size ratio are presented. While the large A particles show a glass transition around the critical density of mode-coupling theory ρ c, the small B particles remain mobile with a relatively weak decrease in their self-diffusion coefficient DB with increasing density. Surprisingly, around ρ c, the self-diffusion coefficient of species A, DA, also starts to show a rather weak dependence on density. We show that this is due to finite-size effects that can be understood from the analysis of the collective interdiffusion dynamics.

Rheological response of a glass-forming liquid having large bidispersity

Using extensive numerical simulations, we investigate the flow behaviour of a model glass-forming binary mixture whose constituent particles have a large size ratio. The rheological response to applied shear is studied in the regime where the larger species are spatially predominant. We demonstrate that the macroscopic rigidity that emerges with increasing density occurs in the regime where the larger species undergo a glass transition while the smaller species continue to be highly diffusive. We analyse the interplay between the timescale imposed by the shear and the quiescent relaxation dynamics of the two species to provide a microscopic insight into the observed rheological response. Finally, by tuning the composition of the mixture, we illustrate that the systematic insertion of the smaller particles affects the rheology by lowering of viscosity of the system.

"Inner clocks" of glass-forming liquids

Providing a physically sound explanation of aging phenomena in non-equilibrium amorphous materialsis a challenging problem in modern statistical thermodynamics. The slow evolution of physical propertiesafter quenches of control parameters is empirically well interpreted via the concept of material time (orinternal clock), based on the Tool-Narayanaswamy-Moynihan (TNM) model. Yet, the fundamental reasonsof its striking success remain unclear. We propose a microscopic rationale behind the material time onthe basis of the linear laws of irreversible thermodynamics and its extension that treats the correspondingkinetic coefficients as state functions of a slowly evolving material state. Our interpretation is based onthe recognition that the same mathematical structure governs both the Tool model and the recently devel-oped non-equilibrium extension of the self-consistent generalized Langevin equation theory (NE-SCGLE),guided by the universal principles of Onsager's theory of irreversible processes. This identification opensthe way for a generalization of the material-time concept to aging systems where several relaxation modeswith very different equilibration processes must be considered, and partially frozen glasses manifest theappearance of partial ergodicity breaking, and hence materials with multiple very distinct inner clocks.

Efficient event-driven simulations of hard spheres

 Hard spheres are arguably one of the most fundamental model systems in soft matter physics, and hence a common topic of simulation studies. Event-driven simulation methods provide an efficient method for studying the phase behavior and dynamics of hard spheres under a wide range of different conditions. Here, we examine the impact of several optimization strategies for speeding up event-driven molecular dynamics of hard spheres and present a light-weight simulation code that outperforms existing simulation codes over a large range of system sizes and packing fractions. The presented differences in simulation speed, typically a factor of five to ten, save significantly on both CPU time and energy consumption and may be a crucial factor for studying slow processes such as crystal nucleation and glassy dynamics.

Connecting Packing Efficiency of Binary Hard Sphere Systems to Their Intermediate Range Structure

Using computed x-ray tomography we determine the three dimensional (3D) structure of binary hard sphere mixtures as a function of composition and size ratio of the particles q. Using a recently introduced four-point correlation function we reveal that this 3D structure has on intermediate and large length scales a surprisingly regular order, the symmetry of which depends on q. The related structural correlation length has a minimum at the composition at which the packing fraction is highest. At this composition also the number of different local particle arrangements has a maximum, indicating that efficient packing of particles is associated with a structure that is locally maximally disordered.

Non-equilibrium relaxation and aging in the dynamics of a dipolar fluid quenched towards the glass transition

The recently developed non-equilibrium self-consistent generalized Langevin equation theory of the dynamics of liquids of non-spherically interacting particles [2016J. Phys. Chem. B1207975] is applied to the description of the irreversible relaxation of a thermally and mechanically quenched dipolar fluid. Specifically, we consider a dipolar hard-sphere liquid quenched (att_w= 0) from full equilibrium conditions towards different ergodic-non-ergodic transitions. Qualitatively different scenarios are predicted by the theory for the time evolution of the system after the quench (t_w> 0), that depend on both the kind of transition approached and the specific features of the protocol of preparation. Each of these scenarios is characterized by the kinetics displayed by a set of structural correlations, and also by the development of two characteristic times describing the relaxation of the translational and rotational dynamics, allowing us to highlight the crossover from equilibration to aging in the system and leading to the prediction of different underlying mechanisms and relaxation laws for the dynamics at each of the glass transitions explored.

Theory of the effect of external stress on the activated dynamics and transport of dilute penetrants in supercooled liquids and glasses

We generalize the self-consistent cooperative hopping theory for a dilute spherical penetrant or tracer activated dynamics in dense metastable hard sphere fluids and glasses to address the effect of external stress, the consequences of which are systematically established as a function of matrix packing fraction and penetrant-to-matrix size ratio. All relaxation processes speed up under stress, but the difference between the penetrant and matrix hopping (alpha relaxation) times decreases significantly with stress corresponding to less time scale decoupling. A dynamic crossover occurs at a critical "slaving onset" stress beyond which the matrix activated hopping relaxation time controls the penetrant hopping time. This characteristic stress increases (decreases) exponentially with packing fraction (size ratio) and can be well below the absolute yield stress of the matrix. Below the slaving onset, the penetrant hopping time is predicted to vary exponentially with stress, differing from the power law dependence of the pure matrix alpha time due to system-specificity of the stress-induced changes in the penetrant local cage and elastic barriers. An exponential growth of the penetrant alpha relaxation time with size ratio under stress is predicted, and at a fixed matrix packing fraction, the exponential relation between penetrant hopping time and stress for different size ratios can be collapsed onto a master curve. Direct connections between the short- and long-time activated penetrant dynamics and between the penetrant (or matrix) alpha relaxation time and matrix thermodynamic dimensionless compressibility are also predicted. The presented results should be testable in future experiments and simulations.

Effect of softness on glass melting and re-entrant solidification in mixtures of soft and hard colloids

We present a systematic investigation of the structure and dynamic properties of model soft-hard colloidal mixtures. Results of a coarse-grained theoretical model are contrasted with rheological data, where the soft and hard colloids are mimicked by large star polymers with high functionality as the soft component and smaller stars with ultrahigh functionality as the hard one. Previous work by us revealed the recovery of the ergodicity of glassy soft star solutions and subsequent arrested phase separation and re-entrant solid transition upon progressive addition of small hard depletants. Here, we use different components to show that a small variation in softness has a significant impact on the state diagram of such mixtures. In particular, we establish that rendering the soft component more penetrable and modifying the size ratio bring about a remarkable shift in both the phase separation region and the glass-melting line so that the region of restored ergodicity can be notably enhanced and extended to much higher star polymer concentrations than for pure systems. We further rationalize our findings by analyzing the features of the depletion interaction induced by the smaller component that result from the interplay between the size ratio and the softness of the large component. These results demonstrate the great sensitivity of the phase behavior of entropic mixtures to small changes in the molecular architecture of the soft stars and point to the importance of accounting for details of the internal microstructure of soft colloidal particles for tailoring the flow properties of soft composites.

Reorientational dynamics in highly asymmetric binary low-molecular mixtures-A quantitative comparison of dielectric and NMR spectroscopy results

Previously, we scrutinized the dielectric spectra of a binary glass former made by a low-molecular high-T_g component 2-(m-tertbutylphenyl)-2'-tertbutyl-9,9'-spirobi[9H]fluorene (m-TPTS; T_g = 350 K) and low-T_g tripropyl phosphate (TPP; T_g = 134 K) [Körber et al., Phys. Chem. Chem. Phys. 23, 7200 (2021)]. Here, we analyze nuclear magnetic resonance (NMR) spectra and stimulated echo decays of deuterated m-TPTS-d₄ (²H) and TPP (³¹P) and attempt to understand the dielectric spectra in terms of component specific dynamics. The high-T_g component (α₁) shows relaxation similar to that of neat systems, yet with some broadening upon mixing. This correlates with high-frequency broadening of the dielectric spectra. The low-T_g component (α₂) exhibits highly stretched relaxations and strong dynamic heterogeneities indicated by "two-phase" spectra, reflecting varying fractions of fast and slow liquid-like reorienting molecules. Missing for the high-T_g component, such two-phase spectra are identified down to w_TPP = 0.04, indicating that isotropic reorientation prevails in the rigid high-T_g matrix stretching from close to T_g ^TPP to T_g1 w_TPP. This correlates with low-frequency broadening of the dielectric spectra. Two T_g values are defined: T_g1 (w_TPP) displays a plasticizer effect, whereas T_g2 (w_TPP) passes through a maximum, signaling extreme separation of the component dynamics at low w_TPP. We suggest understanding the latter counter-intuitive feature by referring to a crossover from "single glass" to "double glass" scenario revealed by recent MD simulations. Analyses reveal that a second population of TPP molecules exists, which is associated with the dynamics of the high-T_g component. However, the fractions are lower than suggested by the dielectric spectra. We discuss this discrepancy considering the role of collective dynamics probed by dielectric but not by NMR spectroscopy.

Waiting-time dependent non-equilibrium phase diagram of simple glass- and gel-forming liquids

Under numerous circumstances, many soft and hard materials are present in a puzzling wealth of non-equilibrium amorphous states, whose properties are not stationary and depend on preparation. They are often summarized in unconventional "phase diagrams" that exhibit new "phases" and/or "transitions" in which time, however, is an essential variable. This work proposes a solution to the problem of theoretically defining and predicting these non-equilibrium phases and their time-evolving phase diagrams, given the underlying molecular interactions. We demonstrate that these non-equilibrium phases and the corresponding non-stationary (i.e., aging) phase diagrams can indeed be defined and predicted using the kinetic perspective of a novel non-equilibrium statistical mechanical theory of irreversible processes. This is illustrated with the theoretical description of the transient process of dynamic arrest into non-equilibrium amorphous solid phases of an instantaneously quenched simple model fluid involving repulsive hard-sphere plus attractive square well pair interactions.

Particle sizing for flowing colloidal suspensions using flow-differential dynamic microscopy

Particle size is a key variable in understanding the behaviour of the particulate products that underpin much of our modern lives. Typically obtained from suspensions at rest, measuring the particle size under flowing conditions would enable advances for in-line testing during manufacture and high-throughput testing during development. However, samples are often turbid, multiply scattering light and preventing the direct use of common sizing techniques. Differential dynamic microscopy (DDM) is a powerful technique for analysing video microscopy of such samples, measuring diffusion and hence particle size without the need to resolve individual particles while free of substantial user input. However, when applying DDM to a flowing sample, diffusive dynamics are rapidly dominated by flow effects, preventing particle sizing. Here, we develop "flow-DDM", a novel analysis scheme that combines optimised imaging conditions, a drift-velocity correction and modelling of the impact of flow. Flow-DDM allows a decoupling of flow from diffusive motion that facilitates successful particle size measurements at flow speeds an order of magnitude higher than for DDM. We demonstrate the generality of the technique by applying flow-DDM to two separate microscopy methods and flow geometries.

Multiple glass transitions and higher-order replica symmetry breaking of binary mixtures

We extend the replica liquid theory in order to describe the multiple glass transitions of binary mixtures with large size disparities, by taking into account the two-step replica symmetry breaking (2RSB). We determine the glass phase diagram of the mixture of large and small particles in the large-dimension limit where the mean-field theory becomes exact. When the size ratio of particles is beyond a critical value, the theory predicts three distinct glass phases; (i) the one-step replica symmetery breaking (1RSB) double glass where both components vitrify simultaneously, (ii) the 1RSB single glass where only large particles are frozen while small particles remain mobile, and (iii) a glass phase called the 2RSB double glass where both components vitrify simultaneously but with an energy landscape topography distinct from the 1RSB double glass.

Activated penetrant dynamics in glass forming liquids: size effects, decoupling, slaving, collective elasticity and correlation with matrix compressibility

We employ the microscopic self-consistent cooperative hopping theory of penetrant activated dynamics in glass forming viscous liquids and colloidal suspensions to address new questions over a wide range of high matrix packing fractions and penetrant-to-matrix particle size ratios. The focus is on the mean activated relaxation time of smaller tracers in a hard sphere fluid of larger particle matrices. This quantity also determines the penetrant diffusion constant and connects directly with the structural relaxation time probed in an incoherent dynamic structure factor measurement. The timescale of the non-activated fast dissipative process is also studied and is predicted to follow power laws with the contact value of the penetrant-matrix pair correlation function and the penetrant-matrix size ratio. For long time penetrant relaxation, in the relatively lower packing fraction metastable regime the local cage barriers are dominant and matrix collective elasticity effects unimportant. As packing fraction and/or penetrant size grows, much higher barriers emerge and the collective elasticity associated with the correlated matrix dynamic displacement that facilitates penetrant hopping becomes important. This results in a non-monotonic variation with packing fraction of the degree of decoupling between the matrix and penetrant alpha relaxation times. The conditions required for penetrant hopping to become slaved to the matrix alpha process are determined, which depend mainly on the penetrant to matrix particle size ratio. By analyzing the absolute and relative importance of the cage and elastic barriers we establish a mechanistic understanding of the origin of the predicted exponential growth of the penetrant hopping time with size ratio predicted at very high packing fractions. A dynamics-thermodynamics power law connection between the penetrant activation barrier and the matrix dimensionless compressibility is established as a prediction of theory, with different scaling exponents depending on whether matrix collective elasticity effects are important. Quantitative comparisons with simulations of the penetrant relaxation time, diffusion constant, and transient localization length of tracers in dense colloidal suspensions and cold viscous liquids reveal good agreements. Multiple new predictions are made that are testable via future experiments and simulations. Extension of the theoretical approach to more complex systems of high experimental interest (nonspherical molecules, semiflexible polymers, crosslinked networks) interacting via variable hard or soft repulsions and/or short range attractions is possible, including under external deformation.

On the Size Effect of Additives in Amorphous Shape Memory Polymers

Small additive molecules often enhance structural relaxation in polymers. We explore this effect in a thermoplastic shape memory polymer via molecular dynamics simulations. The additive-to-monomer size ratio is shown to play a key role here. While the effect of additive-concentration on the rate of shape recovery is found to be monotonic in the investigated range, a non-monotonic dependence on the size-ratio emerges at temperatures close to the glass transition. This work thus identifies the additives' size to be a qualitatively novel parameter for controlling the recovery process in polymer-based shape memory materials.

Observation of soft glassy behavior in a magnetic colloid exposed to an external magnetic field

We provide the first experimental evidence for soft glassy behavior in a sterically stabilized magnetic colloid (ferrofluid) of relatively low volume fraction (φ = 0.037) when a uniform magnetic field is applied at a sufficiently high rate (fast quench). Fast magnetic-field quenches favor structural arrest of field-induced aggregates, owing to insufficient time to settle into lower energy states, thereby pushing the system to a frustrated metastable configuration like a repulsive glass. Brownian dynamics simulations are used to show that the polydisperse ferrofluid (as in experiments) forms thick ropes aligned along the field direction, while a monodisperse ferrofluid does not. The simulations show that there is practically no ordering of the thin, monodisperse chains, while the thick, polydisperse ropes show positional ordering with a typical center-center separation between the particles in different ropes of about 0.39 μm. As a consequence of structural arrest, the ferrofluid exhibits aging with broken time-translational invariance, a hallmark of glassy dynamics. The superposition of strain and creep compliance curves obtained from rheological measurements at different waiting times in the effective time domain corroborates the soft glassy behavior when exposed to a magnetic field applied at a fast ramp rate.

Rheology of colloidal and metallic glass formers

Colloidal hard-sphere suspensions are convenient experimental models to understand soft matter, and also by analogy the structural-relaxation behavior of atomic or small-molecular fluids. We discuss this analogy for the flow and deformation behavior close to the glass transition. Based on a mapping of temperature to effective hard-sphere packing, the stress–strain curves of typical bulk metallic glass formers can be quantitatively compared with those of hard-sphere suspensions. Experiments on colloids give access to the microscopic structure under deformation on a single-particle level, providing insight into the yielding mechanisms that are likely also relevant for metallic glasses. We discuss the influence of higher-order angular signals in connection with non-affine particle rearrangements close to yielding. The results are qualitatively explained on the basis of the mode-coupling theory. We further illustrate the analogy of pre-strain dependence of the linear-elastic moduli using data on PS-PNiPAM suspensions.

Multiple character of non-monotonic size-dependence for relaxation dynamics in polymer-particle and binary mixtures

Adding plasticizers is a well-known procedure to reduce the glass transition temperature in polymers. It has been recently shown that this effect shows a non-monotonic dependence on the size of additive molecules (2019 J. Chem. Phys. 150 024903). In this work, we demonstrate that, as the size of the additive molecules is changed at fixed concentration, multiple extrema emerge in the dependence of the system's relaxation time on the size ratio. The effect occurs on all relevant length scales including single monomer dynamics, decay of Rouse modes and relaxation of the chain's end-to-end vector. A qualitatively similar trend is found within mode-coupling theoretical results for a binary hard-sphere mixture. An interpretation of the effect in terms of local packing efficiency and coupling between the dynamics of minority and majority species is provided.

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.

Glass-transition asymptotics in two theories of glassy dynamics: Self-consistent generalized Langevin equation and mode-coupling theory

We contrast the generic features of structural relaxation close to the idealized glass transition that are predicted by the self-consistent generalized Langevin equation theory (SCGLE) against those that are predicted by the mode-coupling theory of the glass transition (MCT). We present an asymptotic solution close to conditions of kinetic arrest that is valid for both theories, despite the different starting points that are adopted in deriving them. This in particular provides the same level of understanding of the asymptotic dynamics in the SCGLE as was previously done only for MCT. We discuss similarities and different predictions of the two theories for kinetic arrest in standard glass-forming models, as exemplified through the hard-sphere system. Qualitative differences are found for models where a decoupling of relaxation modes is predicted, such as the generalized Gaussian core model, or binary hard-sphere mixtures of particles with very disparate sizes. These differences, which arise in the distinct treatment of the memory kernels associated to self- and collective motion of particles, lead to distinct scenarios that are predicted by each theory for partially arrested states and in the vicinity of higher-order glass-transition singularities.

Effect of size disparity on the structure and dynamics of the small component in concentrated binary colloidal mixtures

We determined, using confocal microscopy, the structure and dynamics of the small component in concentrated binary colloidal mixtures with moderate and large size ratios and different compositions of Polymethyl methacrylate particles. We show that when increasing the content of small spheres at fixed total volume fraction, a transition in the local environment of the small particles is observed, from a mixed environment of other small and large particles to a local environment of only small particles. The transition is rather abrupt for moderate size ratios, while it becomes particularly broad for large size ratios. This can be associated with the improved ability of the small particles to pack in between the large particles for larger size ratios. The dynamics reflect the transition with an increase of the mobility observed at intermediate mixing. This increase becomes particularly pronounced for large size ratios, leading to diffusive dynamics of the small particles, in agreement with predictions of theories of the glass transition in binary hard-sphere mixtures. The composition at which the fastest dynamics are observed is apparently independent of the size ratio.