Anatomy of cage formation in a two-dimensional glass-forming liquid

The cooperative migration dynamics of particles correlates to the nature of hexatic-isotropic phase transition in 2D systems of corner-rounded hexagons

In the two-dimensional (2D) melting transition of colloidal systems, the hexatic-isotropic (H-I) transition can be either first-order or continuous. However, how particle dynamics differs at the single-particle level during these two different melting transitions remains to be disclosed. In this work, by Brownian dynamics (BD) simulations, we have systematically studied the dynamic behavior of corner-rounded hexagons during the H-I transition, for a range of corner-roundness ζ = 0.40 to 0.99 that covers the crossover from the continuous to first-order nature of H-I transition. The results show that hexagons with ζ ≤ 0.5 display a continuous H-I transition, whereas those with ζ ≥ 0.6 demonstrate a first-order H-I transition. Dynamic analysis shows different evolution pathways of the dominant cluster formed by migrating particles, which results in a droplet-like cluster structure for ζ = 0.40 hexagons and a tree-like cluster structure for ζ = 0.99 hexagons. Further investigations on the hopping activities of particles suggest a cooperative origin of migrating clusters. Our work provides a new aspect to understand the dependence of the nature of H-I transition on the roundness of hexagons through particle dynamic behavior.

Anisotropic-Isotropic Transition of Cages at the Glass Transition

Characterizing the local structural evolution is an essential step in understanding the nature of glass transition. In this work, we probe the evolution of Voronoi cell geometry in simple glass models by simulations and colloid experiments, and find that the individual particle cages deform anisotropically in supercooled liquid and isotropically in glass. We introduce an anisotropy parameter k for each Voronoi cell, whose mean value exhibits a sharp change at the mode-coupling glass transition ϕ_{c}. Moreover, a power law of packing fraction ϕ∝q_{1}^{d} is discovered in the supercooled liquid regime with d>D, in contrast to d=D in the glass regime, where q_{1} is the first peak position of structure factor, and D is the space dimension. This power law is qualitatively explained by the change of k. The active motions in supercooled liquid are spatially correlated with long axes rather than short axes of Voronoi cells. In addition, the dynamic slowing down approaching the glass transition can be well characterized through a modified free-volume model based on k. These findings reveal that the structural parameter k is effective in identifying the structure-dynamics correlations and the glass transition in these systems.

Different glassy characteristics are related to either caging or dynamical heterogeneity

Despite the enormous theoretical and application interests, a fundamental understanding of the glassy dynamics remains elusive. The static properties of glassy and ordinary liquids are similar, but their dynamics are dramatically different. What leads to this difference is the central puzzle of the field. Even the primary defining glassy characteristics, their implications, and if they are related to a single mechanism remain unclear. This lack of clarity is a severe hindrance to theoretical progress. Here, we combine analytical arguments and simulations of various systems in different dimensions and address these questions. Our results suggest that the myriad of glassy features are manifestations of two distinct mechanisms. Particle caging controls the mean, and coexisting slow- and fast-moving regions govern the distribution of particle displacements. All the other glassy characteristics are manifestations of these two mechanisms; thus, the Fickian yet non-Gaussian nature of glassy liquids is not surprising. We discover a crossover, from stretched exponential to a power law, in the behavior of the overlap function. This crossover is prominent in simulation data and forms the basis of our analyses. Our results have crucial implications on how the glassy dynamics data are analyzed, challenge some recent suggestions on the mechanisms governing glassy dynamics, and impose strict constraints that a correct theory of glasses must have.

Learning molecular dynamics: predicting the dynamics of glasses by a machine learning simulator

Many-body dynamics of atoms such as glass dynamics is generally governed by complex (and sometimes unknown) physics laws. This challenges the construction of atom dynamics simulations that both (i) capture the physics laws and (ii) run with little computation cost. Here, based on graph neural network (GNN), we introduce an observation-based graph network (OGN) framework to "bypass all physics laws" to simulate complex glass dynamics solely from their static structure. By taking the example of molecular dynamics (MD) simulations, we successfully apply the OGN to predict atom trajectories evolving up to a few hundred timesteps and ranging over different families of complex atomistic systems, which implies that the atom dynamics is largely encoded in their static structure in disordered phases and, furthermore, allows us to explore the capacity of OGN simulations that is potentially generic to many-body dynamics. Importantly, unlike traditional numerical simulations, the OGN simulations bypass the numerical constraint of small integration timestep by a multiplier of ≥5 to conserve energy and momentum until hundreds of timesteps, thus leapfrogging the execution speed of MD simulations for a modest timescale.

Transition between Different Diffusion Modes of Individual Lipids during the Membrane-Specific Action of As-CATH4 Peptides

The cell membrane permeabilization ability of immune defense antimicrobial peptides (AMPs) is widely applied in biomedicine. Although the mechanisms of peptide-membrane interactions have been widely investigated, analyses at the molecular level are still lacking. Herein, the membrane-specific action of a native AMP, As-CATH4, is investigated using a single-lipid tracking method in combination with live cell and model membrane assays conducted at different scales. The peptide-membrane interaction process is characterized by analyzing single-lipid diffusion behaviors. As-CATH4 exhibits potent antimicrobial activity through bacterial membrane permeabilization, with moderate cytotoxicity against mammalian cells. In-plane diffusion analyses of individual lipids show that the lipid molecules exhibit non-Gaussian and heterogeneous diffusion behaviors in both pristine and peptide-treated membranes, which can be decomposed into two Gaussian subgroups corresponding to normal- and slow-diffusive lipids. Assessment of the temporal evolution of these two diffusion modes of lipids reveal that the peptide action states of As-CATH4 include surface binding, transmembrane defect formation, and dynamic equilibrium. The action mechanisms of As-CATH4 at varying concentrations and against different membranes are distinguished. This work resolves the simultaneous mixed diffusion mechanisms of single lipids in biomimetic cell membranes, especially during dynamic membrane permeabilization by AMPs.

Comparative study of the dynamics of colloidal glass and gel

We investigate and compare the difference in the dynamics of two arrested states: colloidal glass and colloidal gel. Real-space experiments reveal two distinct nonergodicity origins for their slow dynamics, namely, cage effects for the glass and attractive bondings for the gel. Such distinct origins lead to a faster decay of the correlation function and a smaller nonergodicity parameter of the glass than those of the gel. We also find that the gel exhibits stronger dynamical heterogeneity compared with the glass due to the greater correlated motions in the gel. Moreover, a logarithmic decay in the correlation function is observed as the two nonergodicity origins merge, consistent with the mode coupling theory.

Non-Gaussian Diffusion of Individual Lipids Unveils the Unique Peptide-Membrane Interaction Dynamics

The dynamics of protein (or peptide)-membrane interactions plays a central role in cellular functions; however, the underlying mechanisms remain unclear. Herein, through analyzing the diffusion of individual lipids in a bilayer membrane during the membrane actions of typical peptides (e.g., pore-forming peptide melittin and cell-penetrating peptide TAT) at varying concentrations, the spatial heterogeneity as well as temporal dynamics of lipid motions were investigated which showed close correlation with the peptide action mechanism. Specifically, the spatial heterogeneity of lipid diffusion was characterized by the non-Gaussianity of lipid trajectories, which was further decomposed into two basic diffusion modes; moreover, the temporal evolution of the Gaussian fitting parameters provided quantitative information on the varying metastable interaction states between peptides and the membrane (e.g., peptide landing, membrane insertion, and equilibrium). Generally, this work gives an insight into the correlation between single-lipid diffusion and function-realization of membrane-active peptides.

A fractal structural feature related to dynamic crossover in metallic glass-forming liquids

The dynamic crossover in supercooled liquids initially predicted by model coupling theory has been widely accepted, but its underlying structural origin is still an open issue for glass-forming liquids. By molecular dynamics simulations of binary CuZr liquids, the present work verifies that high pressure could enhance this crossover, facilitating the studies on the structural features at the crossover temperature T_c. We discover that the topological connectivity of icosahedral clusters is responsible for this dynamic crossover, rather than all clusters. T_c is the temperature at which the connectivity degree between these clusters reaches a maximum and the dynamic heterogeneity begins to keep stable. Below T_c, the fractal topological structures appear in the medium-range order scale. The icosahedral clusters with a certain connectivity pattern can be regarded as a fractal structural unit. By employing the established fractal analysis method, the fractal dimension D of the icosahedral network is calculated. Our results indicate that the D value increases monotonically with increasing pressure and the fractal behavior of the icosahedral network is an inherent feature of metallic glasses. We also find similar fractal behavior in clusters with high local five-fold symmetry. Our findings shed light on the origin of a dynamic crossover in the deep supercooled region of metallic glasses and also demonstrate the important role of icosahedral clusters in uncovering the fractal behavior of metallic glass.

Nonmonotonic Dynamical Correlations beneath the Surface of Glass-Forming Liquids

Collective motion over increasing length scales is a signature of the vitrification process of liquids. We demonstrate how distinct static and dynamic length scales govern the dynamics of vitrifying films. In contrast to a monotonically growing static correlation length, the dynamical correlation length that measures the extent of surface-dynamics acceleration into the bulk displays a striking nonmonotonic temperature evolution that is robust also against changes in detailed interatomic interaction. This nonmonotonic change defines a crossover temperature T_{*} that is distinct from the critical temperature T_{c} of mode-coupling theory. We connect this nonmonotonic change to a morphological change of cooperative rearrangement regions of fast particles, and to the point where the decoupling of fast-particle motion from the bulk relaxation is most sensitive to fluctuations. We propose a rigorous definition of this new crossover temperature T_{*} within a recent extension of mode-coupling theory, the stochastic β-relaxation theory.

Surface melting of a colloidal glass

Despite their technological relevance, a full microscopic understanding of glasses is still lacking. This applies even more to their surfaces whose properties largely differ from that of the bulk material. Here, we experimentally investigate the surface of a two-dimensional glass as a function of the effective temperature. To yield a free surface, we use an attractive colloidal suspension of micron-sized particles interacting via tunable critical Casimir forces. Similar to crystals, we observe surface melting of the glass, i.e., the formation of a liquid film at the surface well below the glass temperature. Underneath, however, we find an unexpected region with bulk density but much faster particle dynamics. It results from connected clusters of highly mobile particles which are formed near the surface and deeply percolate into the underlying material. Because its thickness can reach several tens of particle diameters, this layer may elucidate the poorly understood properties of thin glassy films which find use in many technical applications.

Understanding enhanced rotational dynamics of active probes in rod suspensions

Active Brownian particles (APs) have recently been shown to exhibit enhanced rotational diffusion (ERD) in complex fluids. Here, we experimentally observe ERD and numerically corroborate its microscopic origin for a quasi-two-dimensional suspension of colloidal rods. At high density, the rods form small rafts, wherein they perform small-amplitude, high-frequency longitudinal displacements. Activity couples AP-rod contacts to reorientation, with the variance therein leading to ERD. This is captured by a local, rather than a global relaxation time, as used in previous phenomenological modeling. Our result should prove relevant to the microrheological characterization of complex fluids and furthering our understanding of the dynamics of microorganisms in such media.

Local Plastic Response and Slow Heterogeneous Dynamics of Supercooled Liquids

We demonstrate, via numerical simulations, that the relaxation dynamics of supercooled liquids correlates well with a plastic length scale measuring a particle's response to impulsive localized perturbations and weakly to measures of local elasticity. We find that the particle averaged plastic length scale vanishes linearly in temperature and controls the super-Arrhenius temperature dependence of the relaxation time. Furthermore, we show that the plastic length scale of individual particles correlates with their typical displacement at the relaxation time. In contrast, the local elastic response only correlates with the dynamics on the vibrational timescale.

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.

Controlling the volume fraction of glass-forming colloidal suspensions using thermosensitive host "mesogels"

The key parameter controlling the glass transition of colloidal suspensions is φ, the fraction of the sample volume occupied by the particles. Unfortunately, changing φ by varying an external parameter, e.g., temperature T as in molecular glass formers, is not possible, unless one uses thermosensitive colloidal particles, such as the popular poly(N-isopropylacrylamide) (PNiPAM) microgels. These, however, have several drawbacks, including high deformability, osmotic deswelling, and interpenetration, which complicate their use as a model system to study the colloidal glass transition. Here, we propose a new system consisting of a colloidal suspension of non-deformable spherical silica nanoparticles, in which PNiPAM hydrogel spheres of ∼100-200μm size are suspended. These non-colloidal "mesogels" allow for controlling the sample volume effectively available to the silica nanoparticles and hence their φ, thanks to the T-induced change in mesogels' volume. Using optical microscopy, we first show that the mesogels retain their ability to change size with T when suspended in Ludox suspensions, similarly as in water. We then show that their size is independent of the sample thermal history such that a well-defined, reversible relationship between T and φ may be established. Finally, we use space-resolved dynamic light scattering to demonstrate that, upon varying T, our system exhibits a broad range of dynamical behaviors across the glass transition and beyond, comparable with those exhibited by a series of distinct silica nanoparticle suspensions of various φ.

Local symmetry predictors of mechanical stability in glasses

The mechanical properties of crystals are controlled by the translational symmetry of their structures. But for glasses with a disordered structure, the link between the symmetry of local particle arrangements and stability is not well established. In this contribution, we provide experimental verification that the centrosymmetry of nearest-neighbor polyhedra in a glass strongly correlates with the local mechanical stability. We examine the distribution of local stability and local centrosymmetry in a glass during aging and deformation using microbeam x-ray scattering. These measurements reveal the underlying relationship between particle-level structure and larger-scale behavior and demonstrate that spatially connected, coordinated local transformations to lower symmetry structures are fundamental to these phenomena. While glassy structures lack obvious global symmetry breaking, local structural symmetry is a critical factor in predicting stability.

Translation-rotation decoupling of tracers reflects medium-range crystalline order in two-dimensional colloid glasses

The dynamic heterogeneity and the translation-rotation decoupling are the dynamic signatures of glasses and supercooled liquids. Whether and how the dynamic heterogeneity would relate to the local structure of glasses has been a puzzle for decades. In this work we perform molecular dynamics simulations for tracers in both two-dimensional polydisperse colloids (2DPC) and two-dimensional binary colloids (2DBC). In 2DPC glasses, hexatic local structures develop at low enough temperatures and grow quickly along with the dynamic correlation length of the 2DPC, which is well known as the medium-range crystalline order (MRCO). In 2DBC glasses, on the other hand, any explicit local structure has not been reported to grow significantly with the dynamic correlation length at low temperatures. We introduce two different types of tracers into colloidal systems: A diamond tracer that resembles the MRCO of 2DPC glasses and a square tracer that is dissimilar to any local structure of glasses. The translation-rotation decoupling of the diamond tracer in 2DPC glasses is much more significant than that of the square tracer in the same 2DPC glasses. On the other hand, such a tracer shape-dependence of the decoupling is not observed in 2DBC glasses where the local hexatic structure does not develop significantly. We introduce a shape-dependency parameter of the decoupling and find that the shape-dependency parameter grows along with the dynamic correlation length in 2DPC glasses but not in 2DBC glasses. This illustrates that the dynamic heterogeneity and the translation-rotation decoupling of tracers could reveal the local structure that develops in glasses.

Generalized mode-coupling theory for mixtures of Brownian particles

Generalized mode-coupling theory (GMCT) has recently emerged as a promising first-principles theory to study the poorly understood dynamics of glass-forming materials. Formulated as a hierarchical extension of standard mode-coupling theory (MCT), it is able to systematically improve its predictions by including the exact dynamics of higher-order correlation functions into its hierarchy. However, in contrast to Newtonian dynamics, a fully generalized version of the theory based on Brownian dynamics is still lacking. To close this gap, we provide a detailed derivation of GMCT for colloidal mixtures obeying a many-body Smoluchowski equation. We demonstrate that a hierarchy of coupled equations can again be established and show that these, consistent with standard MCT, are identical to the ones obtained from Newtonian GMCT when taking the overdamped limit. Consequently, the nontrivial similarity between Brownian and Newtonian MCT is maintained for our multicomponent GMCT. As a proof of principle, we also solve the generalized mode-coupling equations for the binary Kob-Andersen Lennard-Jones mixture undergoing Brownian dynamics and confirm the improved predictive power of the theory upon using more levels of the GMCT hierarchy of equations.

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.

Enhanced Heterogeneous Diffusion of Nanoparticles in Semiflexible Networks

The transport of nanoparticles in semiflexible networks, which form diverse principal structural components throughout living systems, is important in biology and biomedical applications. By combining large-scale molecular simulations as well as theoretical analysis, we demonstrate here that nanoparticles in polymer networks with semiflexible strands possess enhanced heterogeneous diffusion characterized by more evident hopping dynamics. Particularly, the hopping energy barrier approximates to linear dependence on confinement parameters in the regime of moderate rigidity, in contrast to the quadratic dependence of both its soft and hard counterparts. This nonmonotonic feature can be attributed to the competition between the conformation entropy and the bending energy regulated by the chain rigidity, captured by developing an analytical model of a hopping energy barrier. Moreover, these theoretical results agree reasonably well with previous experiments. The findings bear significance in unraveling the fundamental physics of substance transport confined in network-topological environments and would provide an explanation for the dynamics diversity of nanoparticles within various networks, biological or synthetic.

Autonomously Probing Viscoelasticity in Disordered Suspensions

Recent experiments show a strong rotational diffusion enhancement for self-propelled microrheological probes in colloidal glasses. Here, we provide microscopic understanding using simulations with a frictional probe-medium coupling that converts active translation into rotation. Diffusive enhancement emerges from the medium's disordered structure and peaks at a second-order transition in the number of contacts. Our results reproduce the salient features of the colloidal glass experiment and support an effective description that is applicable to a broader class of viscoelastic suspensions.