Signatures of dynamical heterogeneity in the structure of glassy free-energy minima

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.

Non-hyperuniform metastable states around a disordered hyperuniform state of densely packed spheres: stochastic density functional theory at strong coupling

The disordered and hyperuniform structures of densely packed spheres near and at jamming are characterized by vanishing of long-wavelength density fluctuations, or equivalently by long-range power-law decay of the direct correlation function (DCF). We focus on previous simulation results that exhibit the degradation of hyperuniformity in jammed structures while maintaining the long-range nature of the DCF to a certain length scale. Here we demonstrate that the field-theoretic formulation of stochastic density functional theory is relevant to explore the degradation mechanism. The strong-coupling expansion method of stochastic density functional theory is developed to obtain the metastable chemical potential considering the intermittent fluctuations in dense packings. The metastable chemical potential yields the analytical form of the metastable DCF that has a short-range cutoff inside the sphere while retaining the long-range power-law behavior. It is confirmed that the metastable DCF provides the zero-wavevector limit of the structure factor in quantitative agreement with the previous simulation results of degraded hyperuniformity. We can also predict the emergence of soft modes localized at the particle scale by plugging this metastable DCF into the linearized Dean-Kawasaki equation, a stochastic density functional equation.

Glass forming liquids in a quenched random potential

The response of a model two-dimensional colloidal glass former to an externally imposed spatially random potential, which acts as a quenched disorder, is investigated using numerical simulations, motivated by recent experiments and also mean field predictions. The external potential induces the onset of the glassy dynamics at increasingly smaller field roughness, with increasing packing fraction of the particulate assembly, and the existence of aging processes within the glassy regime is also observed. Furthermore, along the axis of increasing field roughness, the dynamical slowdown is not correlated to the hexatic order within the supercooled regime.

Multi-scale dynamics at the glassy silica surface

Silica-based glass is a household name, providing insulation for windows to microelectronics. The debate over the types of motions thought to occur in or on SiO₂ glass well below the glass transition temperature continues. Here, we form glassy silica films by oxidizing the Si(100) surface (from 0.5 to 1.5 nm thick, to allow tunneling). We then employ scanning tunneling microscopy in situ to image and classify these motions at room temperature on a millisecond to hour time scale and 50-pm to 5-nm length scale. We observe two phenomena on different time scales. Within minutes, compact clusters with an average diameter of several SiO₂ glass-forming units (GFUs) hop between a few (mostly two) configurations, hop cooperatively (facilitation), and merge into larger clusters (aging) or split into smaller clusters (rejuvenation). Within seconds, Si-O-Si bridges connect two GFUs within a single cluster flip, providing a vibrational fine structure to the energy landscape. We assign the vibrational fine structure using electronic structure calculations. Calculations also show that our measured barrier height for whole cluster hopping at the glass surface (configurational dynamics) is consistent with the configurational entropy predicted by thermodynamic models of the glass transition and that the vibrational entropy for GFU flipping and configurational entropy for cluster hopping are comparable (on a per GFU basis).

Classical density functional theory, unconstrained crystallization, and polymorphic behavior

While in principle, classical density functional theory (cDFT) should be a powerful tool for the study of crystallization, in practice this has not so far been the case. Progress has been hampered by technical problems which have plagued the study of the crystalline systems using the most sophisticated fundamental measure theory models. In this paper, the reasons for the difficulties are examined and it is proposed that the tensor functionals currently favored are in fact numerically unstable. By reverting to an older, more heuristic model it is shown that all of the technical difficulties are eliminated. Application to a Lennard-Jones fluid results in a demonstration of power of cDFT to describe crystallization in a highly inhomogeneous system. First, we show that droplets attached to a slightly hydrophobic wall crystallize spontaneously upon being quenched. The resulting crystallites are clearly faceted structures and are predominantly HCP structures. In contrast, droplets in a fully periodic calculational cell remain stable to lower temperatures and eventually show the same spontaneous localization of the density into "atoms" but in an amorphous structure having many of the structural characteristics of a glass. A small change of the protocol leads, at the same temperature, to the formation of crystals, this time with the fcc structure typical of bulk Lennard-Jones solids. The fcc crystals have lower free energy than the amorphous structures which in turn are more stable than the liquid droplets. It is demonstrated that as the temperature is raised, the free energy differences between the structures decrease until the solid clusters become less stable than the liquid droplets and spontaneously melt. The presence of energy barriers separating the various structures is therefore clearly demonstrated.

Potential energy landscape and inherent dynamics of a hard-sphere fluid

Hard-sphere models exhibit many of the same kinds of supercooled-liquid behavior as more realistic models of liquids, but the highly nonanalytic character of their potentials makes it a challenge to think of that behavior in potential energy landscape terms. We show here that it is possible to calculate an important topological property of hard-sphere landscapes, the geodesic pathways through those landscapes, and to do so without artificially coarse-graining or softening the potential. We show, moreover, that the rapid growth of the lengths of those pathways with increasing packing fraction quantitatively predicts the precipitous decline in diffusion constants in a glass-forming hard-sphere mixture model. The geodesic paths themselves can be considered as defining the intrinsic dynamics of hard spheres, so it is also revealing to find that they (and therefore the features of the underlying potential energy landscape) correctly predict the occurrence of dynamic heterogeneity and nonzero values of the non-Gaussian parameter. The success of these landscape predictions for the dynamics of such a singular model emphasizes that there is more to potential energy landscapes than is revealed by looking at the minima and saddle points.