Shape of dynamical heterogeneities and fractional Stokes-Einstein and Stokes-Einstein-Debye relations in quasi-two-dimensional suspensions of colloidal ellipsoids

Breakdown of the Stokes-Einstein relation in Stillinger-Weber silicon

We investigate the dynamical properties of liquid and supercooled liquid silicon, modeled using the Stillinger-Weber potential, to examine the validity of the Stokes-Einstein (SE) relation. Toward this end, we examine the relationship among various dynamical quantities, including (i) the macroscopic transport coefficients-self-diffusion coefficient D and viscosity η, (ii) relaxation time τα, and (iii) lengthscale dependent relaxation times τα(q) over a broad range of temperature T, pressure P, and density ρ covering both equilibrium and metastable liquid state points in the phase diagram. Our study shows a weak breakdown in the SE relation involving D and η, and the loci of the breakdown of the SE relation (SEB) are found in the high T liquid phase. The τα, when used as a proxy to η, shows a distinct breakdown in the SE relation, whose loci are found in the supercooled liquid phase. Interestingly, certain parts of the phase diagram show that the loci of onset of slow dynamics lie below the loci of SEB, suggesting a regime that exhibits Arrhenius but non-Fickian behavior. Computation of τα(q) enables us to extract the lengthscale associated with the Fickian to non-Fickian behavior using which we show that the breakdown of the SE relation occurs only below a specific lengthscale at a given temperature. Furthermore, we also compare the SEB loci with other features of the phase behavior, including the loci of compressibility maximum, density maximum, and diffusivity maximum.

Direction-dependent dynamics of colloidal particle pairs and the Stokes-Einstein relation in quasi-two-dimensional fluids

Hydrodynamic interactions are important for diverse fluids, especially those with low Reynolds number such as microbial and particle-laden suspensions, and proteins diffusing in membranes. Unfortunately, while far-field (asymptotic) hydrodynamic interactions are fully understood in two- and three-dimensions, near-field interactions are not, and thus our understanding of motions in dense fluid suspensions is still lacking. In this contribution, we experimentally explore the hydrodynamic correlations between particles in quasi-two-dimensional colloidal fluids in the near-field. Surprisingly, the measured displacement and relaxation of particle pairs in the body frame exhibit direction-dependent dynamics that can be connected quantitatively to the measured near-field hydrodynamic interactions. These findings, in turn, suggest a mechanism for how and when hydrodynamics can lead to a breakdown of the ubiquitous Stokes-Einstein relation (SER). We observe this breakdown, and we show that the direction-dependent breakdown of the SER is ameliorated along directions where hydrodynamic correlations are smallest. In total, the work uncovers significant ramifications of near-field hydrodynamics on transport and dynamic restructuring of fluids in two-dimensions.

Structure of saturated random-sequential-adsorption ellipse packings

Motivated by the recent observation of liquid glass in suspensions of ellipsoidal colloids, we examine the structure of (asymptotically) saturated RSA ellipse packings. We determine the packing fractions ϕ_{s}(α) to high precision, finding an empirical analytic formula that predicts ϕ_{s}(α) to within less than 0.1% for all α≤10. Then we explore how these packings' positional-orientational order varies with α. We find a transition from tip/side- to side/side-contact-dominated structure at α=α_{TS}≃2.4. At this aspect ratio, the peak value g_{max} of packings' positional-orientational pair correlation functions is minimal, and systems can be considered maximally locally disordered. For smaller (larger) α, g_{max} increases exponentially with deceasing (increasing) α. Local nematic order and structures comparable to the precursor domains observed in experiments gradually emerge as α increases beyond three. For α≳5, single-layer lamellae become more prominent and long-wavelength density fluctuations increase with α as packings gradually approach the rodlike limit.

Structural origin of excitations in a colloidal glass-former

Despite decades of intense research, whether the transformation of supercooled liquids into glass is a kinetic phenomenon or a thermodynamic phase transition remains unknown. Here, we analyzed optical microscopy experiments on 2D binary colloidal glass-forming liquids and investigated the structural links of a prominent kinetic theory of glass transition. We examined a possible structural origin for localized excitations, which are building blocks of the dynamical facilitation theory—a purely kinetic approach for the glass transition. To accomplish this, we utilize machine learning methods to identify a structural order parameter termed “softness” that has been found to be correlated with reorganization events in supercooled liquids. Both excitations and softness qualitatively capture the dynamical slowdown on approaching the glass transition and motivated us to explore spatial and temporal correlations between them. Our results show that excitations predominantly occur in regions with high softness and the appearance of these high softness regions precedes excitations, thus suggesting a causal connection between them. Thus, unifying dynamical and thermodynamical theories into a single structure-based framework may provide a route to understand the glass transition.

Are strongly confined colloids good models for two dimensional liquids?

Quasi-two-dimensional (quasi-2D) colloidal hard-sphere suspensions confined in a slit geometry are widely used as two-dimensional (2D) model systems in experiments that probe the glassy relaxation dynamics of 2D systems. However, the question to what extent these quasi-2D systems indeed represent 2D systems is rarely brought up. Here, we use computer simulations that take into account hydrodynamic interactions to show that dense quasi-2D colloidal bi-disperse hard-sphere suspensions exhibit much more rapid diffusion and relaxation than their 2D counterparts at the same area fraction. This difference is induced by the additional vertical space in the quasi-2D samples in which the small colloids can move out of the 2D plane, therefore allowing overlap between particles in the projected trajectories. Surprisingly, this difference in the dynamics can be accounted for if, instead of using the surface density, one characterizes the systems by means of a suitable structural quantity related to the radial distribution function. This implies that in the two geometries, the relevant physics for glass formation is essentially identical. Our results provide not only practical implications on 2D colloidal experiments but also interesting insights into the 3D-to-2D crossover in glass-forming systems.

Breakdown of the Stokes-Einstein relation in supercooled liquids: A cage-jump perspective

The breakdown of the Stokes-Einstein relation in supercooled liquids, which is the increase in the ratio τ_ατ_D between the two macroscopic times for structural relaxation and diffusion on decreasing the temperature, is commonly ascribed to dynamic heterogeneities, but a clear-cut microscopic interpretation is still lacking. Here, we tackle this issue exploiting the single-particle cage-jump framework to analyze molecular dynamics simulations of soft disk assemblies and supercooled water. We find that τ_ατ_D∝⟨t_p⟩⟨t_c⟩, where ⟨t_p⟩ and ⟨t_c⟩ are the cage-jump times characterizing slow and fast particles, respectively. We further clarify that this scaling does not arise from a simple term-by-term proportionality; rather, the relations τ_α∝⟨t_p⟩⟨Δr_J ²⟩ and τ_D∝⟨t_c⟩⟨Δr_J ²⟩ effectively connect the macroscopic and microscopic timescales, with the mean square jump length ⟨Δr_J ²⟩ shrinking on cooling. Our work provides a microscopic perspective on the Stokes-Einstein breakdown and generalizes previous results on lattice models to the case of more realistic glass-formers.

Emergence of cooperatively reorganizing cluster and super-Arrhenius dynamics of fragile supercooled liquids

In this paper, we develop a theory to calculate the structural relaxation time τ_{α} of fragile supercooled liquids. Using the information of the configurational entropy and structure, we calculate the number of dynamically free, metastable, and stable neighbors around a central particle. In supercooled liquids, the cooperatively reorganizing clusters (CRCs) in which the stable neighbors form "stable" nonchemical bonds with the central particle emerge. For an event of relaxation to take place, these bonds have to reorganize irreversibly; the energy involved in the processes is the effective activation energy of relaxation. The theory brings forth a temperature T_{a} and a temperature-dependent parameter ψ(T) which characterize slowing down of dynamics on cooling. It is shown that the value of ψ(T) is equal to 1 for T>T_{a}, indicating that the underlying microscopic mechanism of relaxation is dominated by the entropy-driven processes, while for T<T_{a}, ψ(T) decreases on cooling, indicating the emergence of the energy-driven processes. This crossover of ψ(T) from high to low temperatures explains the crossover seen in τ_{α}. The dynamics of systems that may have similar static structure but very different dynamics can be understood in terms of ψ(T). We present results for the Kob-Anderson model for three densities and show that the calculated values of τ_{α} are in excellent agreement with simulation values for all densities. We also show that when ψ(T), τ_{α}, and other quantities are plotted as a function of T/T_{a} (or T_{a}/T), the data collapse on master curves.

Violation of Stokes-Einstein and Stokes-Einstein-Debye relations in polymers at the gas-supercooled liquid coexistence

Molecular dynamics simulations are performed on a system of model linear polymers to look at the violations of Stokes-Einstein (SE) and Stokes-Einstein-Debye (SED) relations near the mode coupling theory transition temperatureT_cat three (one higher and two lower) densities. At low temperatures, both lower density systems show stable gas-supercooled-liquid coexistence whereas the higher density system is homogeneous. We show that monomer density relaxation exhibits SE violation for all three densities, whereas molecular density relaxation shows a weak violation of the SE relation nearT_cin both lower density systems. This study identifies disparity in monomer mobility and observation of jumplike motion in the typical monomer trajectories resulting in the SE violations. In addition to the SE violation, a weak SED violation is observed in the gas-supercooled-liquid coexisting domains of the lower densities. Both lower density systems also show a decoupling of translational and rotational dynamics in this polymer system.

Cooperatively rearranging regions change shape near the mode-coupling crossover for colloidal liquids on a sphere

The structure and dynamics of liquids on curved surfaces are often studied through the lens of frustration-based approaches to the glass transition. Competing glass transition theories, however, remain largely untested on such surfaces and moreover, studies hitherto have been entirely theoretical/numerical. Here we carry out single particle-resolved imaging of dynamics of bi-disperse colloidal liquids confined to the surface of a sphere. We find that mode-coupling theory well captures the slowing down of dynamics in the moderate to deeply supercooled regime. Strikingly, the morphology of cooperatively rearranging regions changed from string-like to compact near the mode-coupling crossover—a prediction unique to the random first-order theory of glasses. Further, we find that in the limit of strong curvature, Mermin–Wagner long-wavelength fluctuations are irrelevant and liquids on a sphere behave like three-dimensional liquids. A comparative evaluation of competing mechanisms is thus an essential step towards uncovering the true nature of the glass transition.

The static and dynamic behavior of condensed phases residing on curved surfaces can be fundamentally different from their counterparts in Euclidean space. Singh et al. test several competing glass theories on colloidal liquids confined to the surface of a sphere and show they behave like 3D bulk liquids.

Rotational and translational dynamics in dense fluids of patchy particles

We explore the effect of directionality on rotational and translational relaxation in glassy systems of patchy particles. Using molecular dynamics simulations, we analyze the impact of two distinct patch geometries, one that enhances the local icosahedral structure and the other one that does not strongly affect the local order. We find that in nearly all investigated cases, rotational relaxation takes place on a much faster time scale than translational relaxation. By comparing to a simplified dynamical Monte Carlo model, we illustrate that rotational diffusion can be qualitatively explained as purely local motion within a fixed environment, which is not coupled strongly to the cage-breaking dynamics required for translational relaxation. Nonetheless, icosahedral patch placement has a profound effect on the local structure of the system, resulting in a dramatic slowdown at low temperatures, which is strongest at an intermediate "optimal" patch size.

Long-wavelength fluctuations and anomalous dynamics in 2-dimensional liquids

In 2-dimensional systems at finite temperature, long-wavelength Mermin-Wagner fluctuations prevent the existence of translational long-range order. Their dynamical signature, which is the divergence of the vibrational amplitude with the system size, also affects disordered solids, and it washes out the transient solid-like response generally exhibited by liquids cooled below their melting temperatures. Through a combined numerical and experimental investigation, here we show that long-wavelength fluctuations are also relevant at high temperature, where the liquid dynamics do not reveal a transient solid-like response. In this regime, these fluctuations induce an unusual but ubiquitous decoupling between long-time diffusion coefficient D and structural relaxation time τ, where [Formula: see text], with [Formula: see text] Long-wavelength fluctuations have a negligible influence on the relaxation dynamics only at extremely high temperatures in molecular liquids or at extremely low densities in colloidal systems.

Correlations between short- and long-time relaxation in colloidal supercooled liquids and glasses

Spatiotemporal dynamics of short- and long-time structural relaxation are measured experimentally as a function of packing fraction, ϕ, in quasi-two-dimensional colloidal supercooled liquids and glasses. The relaxation times associated with long-time dynamic heterogeneity and short-time intracage motion are found to be strongly correlated and to grow by orders of magnitude with increasing ϕ toward dynamic arrest. We find that clusters of fast particles on the two timescales often overlap, and, interestingly, the distribution of minimum-spatial-separation between closest nonoverlapping clusters across the two timescales is revealed to be exponential with a decay length that increases with ϕ. In total, the experimental observations suggest short-time relaxation events are very often precursors to heterogeneous relaxation at longer timescales in glassy materials.

Dynamic Heterogeneities in Colloidal Supercooled Liquids: Experimental Tests of Inhomogeneous Mode Coupling Theory

The dynamics in supercooled liquids slow enormously upon approaching the glass transition, albeit without significant change of liquid structure. This empirical observation has stimulated development of many theoretical models which attempt to elucidate microscopic mechanisms in glasses and glass precursors. Here, quasi-two-dimensional colloidal supercooled liquids and glasses are employed to experimentally test predictions of widely used models: mode coupling theory (MCT) and its important extension, inhomogeneous MCT (IMCT). We measure two-point dynamic correlation functions in the glass forming liquids to determine structural relaxation times, τ_α, and mode coupling exponents, a, b, and γ; these parameters are then used to extract the mode coupling dynamic crossover packing area-fraction, ϕ _c. This information, along with our measurements of supercooled liquid spatiotemporal dynamics, permits characterization of dynamic heterogeneities in the samples and facilitates direct experimental tests of the scaling predictions of IMCT. The time scales at which dynamic heterogeneities are largest, and their spatial sizes, exhibit power law growth on approaching ϕ _c. Within experimental error, the exponents of the measured power laws are close to the predictions of IMCT.

Translational and Rotational Dynamical Heterogeneities in Granular Systems

We use x-ray tomography to investigate the translational and rotational dynamical heterogeneities of a three dimensional hard ellipsoid granular packing driven by oscillatory shear. We find that particles which translate quickly form clusters with a size distribution given by a power law with an exponent that is independent of the strain amplitude. Identical behavior is found for particles that are translating slowly, rotating quickly, or rotating slowly. The geometrical properties of these four different types of clusters are the same as those of random clusters. Different cluster types are considerably correlated or anticorrelated, indicating a significant coupling between translational and rotational degrees of freedom. Surprisingly, these clusters are formed already at time scales that are much shorter than the α-relaxation time, in stark contrast to the behavior found in glass-forming systems.

Measurements of growing surface tension of amorphous-amorphous interfaces on approaching the colloidal glass transition

There is mounting evidence indicating that relaxation dynamics in liquids approaching their glass transition not only become increasingly cooperative, but the relaxing regions also become more compact in shape. Of the many theories of the glass transition, only the random first-order theory—a thermodynamic framework—anticipates the surface tension of relaxing regions to play a role in deciding both their size and morphology. However, owing to the amorphous nature of the relaxing regions, even the identification of their interfaces has not been possible in experiments hitherto. Here, we devise a method to directly quantify the dynamics of amorphous–amorphous interfaces in bulk supercooled colloidal liquids. Our procedure also helped unveil a non-monotonic evolution in dynamical correlations with supercooling in bulk liquids. We measure the surface tension of the interfaces and show that it increases rapidly across the mode-coupling area fraction. Our experiments support a thermodynamic origin of the glass transition.

The existence of interfaces, separating distinct relaxing regions, has been predicted in glass  theory, but a direct proof remains challenging due to the amorphous nature of glasses. Ganapathi et al. identify and measure the surface tension of these interfaces in bulk supercooled colloidal liquids.

Formation of a transient amorphous solid in low density aqueous charged sphere suspensions

Colloidal glasses formed from hard spheres, nearly hard spheres, ellipsoids and platelets or their attractive variants, have been studied in great detail. Complementing and constraining theoretical approaches and simulations, the many different types of model systems have significantly advanced our understanding of the glass transition in general. Despite their early prediction, however, no experimental charged sphere glasses have been found at low density, where the competing process of crystallization prevails. We here report the formation of a transient amorphous solid formed from charged polymer spheres suspended in thoroughly deionized water at volume fractions of 0.0002-0.01. From optical experiments, we observe the presence of short-range order and an enhanced shear rigidity as compared to the stable polycrystalline solid of body centred cubic structure. On a density dependent time scale of hours to days, the amorphous solid transforms into this stable structure. We further present preliminary dynamic light scattering data showing the evolution of a second slow relaxation process possibly pointing to a dynamic heterogeneity known from other colloidal glasses and gels. We compare our findings to the predicted phase behaviour of charged sphere suspensions and discuss possible mechanisms for the formation of this peculiar type of colloidal glass.

Additivity, density fluctuations, and nonequilibrium thermodynamics for active Brownian particles

Using an additivity property, we study particle-number fluctuations in a system of interacting self-propelled particles, called active Brownian particles (ABPs), which consists of repulsive disks with random self-propulsion velocities. From a fluctuation-response relation, a direct consequence of additivity, we formulate a thermodynamic theory which captures the previously observed features of nonequilibrium phase transition in the ABPs from a homogeneous fluid phase to an inhomogeneous phase of coexisting gas and liquid. We substantiate the predictions of additivity by analytically calculating the subsystem particle-number distributions in the homogeneous fluid phase away from criticality where analytically obtained distributions are compatible with simulations in the ABPs.

A molecular dynamics simulations study on the relations between dynamical heterogeneity, structural relaxation, and self-diffusion in viscous liquids

We perform molecular dynamics simulations for viscous liquids to study the relations between dynamical heterogeneity, structural (α) relaxation, and self-diffusion. For atomistic models of supercooled water, polymer melts, and an ionic liquid, we characterize the space-time characteristics of dynamical heterogeneity by the degree of deviations from Gaussian displacement statistics (α2), the size of clusters comprising highly mobile particles (S(w)), and the length of strings consisting of cooperatively moving particles (L(w)). Comparison of our findings with previous simulation results for a large variety of viscous liquids, ranging from monoatomic liquids to silica melt, reveals a nearly universal decoupling between the time scales of maximum non-Gaussian parameter (τ(α2)) and the time constant of the α relaxation (τ(α)) upon cooling, explicitly, τ(α2) ∝τ(α)(3/4). Such uniform relation was not observed between the peak times of S(w) or L(w) and τ(α). On the other hand, the temperature-dependent time scale of maximum string length (τ(L)) follows the inverse of the self-diffusion coefficient (D) for various systems at sufficiently low temperatures, i.e., τ(L) ∝ D(-1). These observations are discussed in view of a breakdown of the Stokes-Einstein relation for the studied systems. It is found that the degree of deviation from this relation is correlated with the stretching of the α relaxation.

Tracer Shape and Local Media Structure Determine the Trend of Translation-Rotation Decoupling in Two-Dimensional Colloids

The translational diffusion of tracers in glass-forming materials often violates the Stokes-Einstein relation while their rotation follows the Debye-Stokes-Einstein relation faithfully, thus decoupling translational and rotational diffusion. In this Letter, we show by performing molecular dynamics simulations for two-dimensional (2D) colloids that the tracer shape and the local media structure are critical such that rotational diffusion is either suppressed or enhanced depending on the tracer shape. For square tracers dissimilar in structure to the local media structure of 2D colloids, the translation-rotation decoupling occurs and the rotational diffusion is enhanced relative to the translation. For sufficiently large diamond tracers similar in structure to the local media structure, tracers undergo rotational hopping motions and their rotation is suppressed relative to the translation. For distorted-diamond tracers, the decoupling is marginal. Translational diffusion does not change significantly with the tracer shape and obeys the Stokes-Einstein relation.