Molecular dynamics studies of heterogeneous dynamics and dynamic crossover in supercooled atomic liquids

Entropic timescales of dynamic heterogeneity in supercooled liquid

Non-Gaussian displacement distributions are universal predictors of dynamic heterogeneity in slowly varying environments. Here, we explore heterogeneous dynamics in supercooled liquid using molecular dynamics simulations and show the efficiency of the relative-entropy based measure, negentropy, in quantifying dynamic heterogeneity over the widely used non-Gaussian parameter. Our analysis shows that the heterogeneity quantified by the negentropy is significantly different from the one obtained using the conventional moment-based definition that considers deviation from Gaussianity up to lower-order moments. We extract the timescales of dynamic heterogeneity using the two methods and show that the differential changes diverge as the system experiences strong intermittency near the glass transition. Further, we interpret the entropic timescales and discuss the general implications of our work.

Calculated solvent reorganization entropies, free energies, and fluctuations of the energy gaps for intramolecular electron transfer and excitation of the solvatochromic dye B30

Intramolecular electron transfer between two biphenyl groups linked by an androstane spacer and excitation of the pyridinium-N-phenolate betaine dye B30 to the first excited singlet state are studied by quantum/classical molecular-dynamics simulations at temperatures between 150 and 300 K in solvents with a range of polarities. Temperature dependences of the solvent reorganization energies, free energies, entropies, and the inhomogeneous broadening of B30's absorption band are examined. The variances of fluctuations of the energy gap between the reactant and product states do not have the direct proportionality to temperature that often is assumed to hold. An explanation for the observations is suggested.

Mesoscopic two-point collective dynamics of glass-forming liquids

The collective density-density and hydrostatic pressure-pressure correlations of glass-forming liquids are spatiotemporally mapped out using molecular dynamics simulations. It is shown that the sharp rise of structural relaxation time below the Arrhenius temperature coincides with the emergence of slow, nonhydrodynamic collective dynamics on mesoscopic scales. The observed long-range, nonhydrodynamic mode is independent of wave numbers and closely coupled to the local structural dynamics. Below the Arrhenius temperature, it dominates the slow collective dynamics on length scales immediately beyond the first structural peak in contrast to the well-known behavior at high temperatures. These results highlight a key connection between the qualitative change in mesoscopic two-point collective dynamics and the dynamic crossover phenomenon.

Automated characterization of spatial and dynamical heterogeneity in supercooled liquids via implementation of machine learning

A computational approach by an implementation of the principle component analysis (PCA) withK-means and Gaussian mixture (GM) clustering methods from machine learning algorithms to identify structural and dynamical heterogeneities of supercooled liquids is developed. In this method, a collection of the average weighted coordination numbers (WCNs‾) of particles calculated from particles' positions are used as an order parameter to build a low-dimensional representation of feature (structural) space forK-means clustering to sort the particles in the system into few meso-states using PCA. Nano-domains or aggregated clusters are also formed in configurational (real) space from a direct mapping using associated meso-states' particle identities with some misclassified interfacial particles. These classification uncertainties can be improved by a co-learning strategy which utilizes the probabilistic GM clustering and the information transfer between the structural space and configurational space iteratively until convergence. A final classification of meso-states in structural space and domains in configurational space are stable over long times and measured to have dynamical heterogeneities. Armed with such a classification protocol, various studies over the thermodynamic and dynamical properties of these domains indicate that the observed heterogeneity is the result of liquid-liquid phase separation after quenching to a supercooled state.

Orientational wetting and dynamical correlations toward glass transition on the surface of imidazolium-based ionic liquids

Surface induces many fascinating physical phenomena, such as dynamic acceleration, surface anchoring, and orientational wetting, and, thus, is of great interest to study. Here, we report classic molecular dynamics simulations on the free-standing surface of imidazolium-based ionic liquids (ILs) [C4mim][PF6] and [C10mim][PF6]. On [C10mim][PF6] surface, a significant orientational wetting is observed, with the wetting strength showing a diverging tendency. Depth of the wetting was captured from the density and orientational order profile by a static length, which remarkably increases below the temperature Tstat upon cooling down. The dynamical correlation length that measures the distance of surface-dynamics acceleration into the bulk was characterized via the spatial-dependent mobility. The translational correlation exhibits a similar drastic increment at Tstat, while the rotational correlation drastically increases at a lower temperature Trot. We connect these results to the dynamics in bulk liquids, by finding Tstat and Trot that correspond to the onset temperatures where the liquids become cooperative for translational and rotational relaxation, respectively. This signifies the importance of collective dynamics in the bulk on the orientational wetting and surface dynamics in the ILs.

Construction of a quantitative relation between structural relaxation and dynamic heterogeneity by vibrational dynamics in glass-forming liquids and polymers

The structural relaxation slows down drastically upon approaching the glass transition, accompanied by the significant growth of dynamic heterogeneity. The fundamental question of elusiveness and interest is whether there exists an underlying quantitative relationship between structural relaxation and dynamic heterogeneity. Here, we reveal that b̃ which is related to the reduced mean square displacements to overcome the energy barriers of activated jumps, instead of the kinetic fragility m, is the genuine key parameter connecting dynamic heterogeneity with structural relaxation for varying types of glass formers. Furthermore, based on the dependence of dynamic heterogeneity on the Debye-Waller factor we obtained a direct quantitative relation between dynamic heterogeneity and structural relaxation is built for different glass-forming liquids. More importantly, a scaling collapse of structural relaxation and dynamic heterogeneity is achieved by the important parameter b̃. These results are of fundamental and critical importance for developing a unified theory of glassy dynamics.

A Dynamically Correlated Network Model for the Collective Dynamics in Glass-Forming Molecular Liquids and Polymers

The non-Arrhenius behavior of segmental dynamics in glass-forming liquids is one of the most profound mysteries in soft matter physics. In this article, we propose a dynamically correlated network (DCN) model to understand the growing behavior of dynamically correlated regions during cooling, which leads to the viscous slowdown of supercooled liquids. The fundamental concept of the model is that the cooperative region of collective motions has a network structure that consists of string-like parts, and networks of various sizes interpenetrate each other. Each segment undergoes dynamical coupling with its neighboring segments via a finite binding energy. Monte Carlo simulations showed that the fractal dimension of the DCNs generated at different temperatures increased and their size distribution became broader with decreasing temperature. The segmental relaxation time was evaluated based on a power law with four different exponents for the activation energy of rearrangement with respect to the DCN size. The results of the present DCN model are consistent with the experimental results for various materials of molecular and polymeric liquids.

Unifying fluctuation-dissipation temperatures of slow-evolving nonequilibrium systems from the perspective of inherent structures

For nonequilibrium systems, how to define temperature is one of the key and difficult issues to solve. Although effective temperatures have been proposed and studied to this end, it still remains elusive what they actually are. Here, we focus on the fluctuation-dissipation temperatures and report that such effective temperatures of slow-evolving systems represent characteristic temperatures of their equilibrium counterparts. By calculating the fluctuation-dissipation relation of inherent structures, we obtain a temperature-like quantity T IS For monocomponent crystal-formers, T IS agrees well with the crystallization temperature T c, while it matches with the onset temperature T on for glass-formers. It also agrees with effective temperatures of typical nonequilibrium systems, such as aging glasses, quasi-static shear flows, and quasi-static self-propelled flows. From the unique perspective of inherent structures, our study reveals the nature of effective temperatures and the underlying connections between nonequilibrium and equilibrium systems and confirms the equivalence between T on and T c.

Jump-precursor state emerges below the crossover temperature in supercooled o-terphenyl

In a supercooled liquid, the crossover temperature T_{c} separates a high-temperature region of diffusive dynamics from a low-temperature region of activated dynamics. A molecular-dynamics simulation of all-atom, flexible o-terphenyl [Eastwood et al., J. Phys. Chem. B 117, 12898 (2013)10.1021/jp402102w] is analyzed with advanced statistical methods to reveal the molecular features associated with this crossover. The simulations extend to an α-relaxation time of 14 μs (272.5 K), two orders of magnitude slower than at T_{c} (290 K). At T_{c} and below, a distinct state emerges that immediately precedes an orientational jump. Compared to the initial, tightly caged state, this jump-precursor state has a looser cage, with solid-angular excursions of 0.054-0.0125 × 4π sr. At T_{c} (290 K), rate heterogeneity is already the dominant cause of stretched relaxation. Exchange within the distribution of rates is faster than α relaxation at T_{c}, but becomes equal to it at the lowest temperature simulated (272.5 K). The results trend toward a recent experimental observation near the glass transition (243 K) [Kaur et al., Phys. Rev. E 98, 040603(R) (2018)10.1103/PhysRevE.98.040603], which saw exchange substantially slower than α relaxation. Overall, the dynamic crossover comprises multiple phenomena: the development of heterogeneity, an increasing jump size, an emerging jump-precursor state, and a lengthening exchange time. The crossover is neither sharp, nor a simple superposition of the high- and low-temperature regimes; it is a broad region that contains unique and complex phenomena.

Connecting relaxation time to a dynamical length scale in athermal active glass formers

Supercooled liquids display dynamics that are inherently heterogeneous in space. This essentially means that at temperatures below the melting point, particle dynamics in certain regions of the liquid can be orders of magnitude faster than other regions. Often dubbed dynamical heterogeneity, this behavior has fascinated researchers involved in the study of glass transition for over two decades. A fundamentally important question in all glass transition studies is whether one can connect the growing relaxation time to a concomitantly growing length scale. In this paper, we go beyond the realm of ordinary glass forming liquids and study the origin of a growing dynamical length scale ξ in a self-propelled "active" glass former. This length scale, which is constructed using structural correlations, agrees well with the average size of the clusters of slow-moving particles that are formed as the liquid becomes spatially heterogeneous. We further report that the concomitantly growing α-relaxation time exhibits a simple scaling law, τ_{α}∼exp(μξ/T_{eff}), with μ as an effective chemical potential, T_{eff} as the effective temperature, and μξ as the growing free energy barrier for cluster rearrangements. The findings of our study are valid over four decades of persistence times, and hence they could be very useful in understanding the slow dynamics of a generic active liquid such as an active colloidal suspension, or a self-propelled granular medium.

Breakdown of the Stokes-Einstein Relation in Supercooled Water/Methanol Binary Mixtures: Explanation Using the Translational Jump-Diffusion Approach

A recent experiment has directly checked the validity of the Stokes-Einstein (SE) relation for pure water, pure methanol, and their binary mixtures of three different compositions at different temperatures. The effect of composition on the nature of breakdown of the SE relation is interesting. While in the majority of the systems, an increasing SE breakdown is observed with decreasing temperature, the breakdown is already significant at higher temperatures for the equimolar mixture. Violations of the SE relation in pure supercooled water at different temperatures and pressures have been previously explained using the translational jump-diffusion (TJD) approach, which provides a fundamental molecular basis, by directly connecting the SE breakdown with jump-diffusion of the molecules. We have used the same TJD approach for explaining the SE breakdown for the methanol/water binary mixtures of compositions studied in the experiment over a wide range of temperatures between 220 K and 300 K. We have understood that the jump-diffusion is the key responsible factor for the SE breakdown. The maximum jump-diffusion contribution gives rise to the early SE breakdown observed for the equimolar mixture observed in the experiment. This study, therefore, provides molecular insight into the SE breakdown for the supercooled water/methanol binary mixture, as found in the experiment.

Elucidation of the Nature of Structural Relaxation in Glassy d-Sorbitol

The nature and origin of the glass transition is one of the great unsolved problems of condensed-matter science. With the rapid increase of viscosity upon cooling the liquid near the glass-transition temperature, a range of dynamical motifs are observed, revealing the sheer complexity of interactions between the amorphous units. Yet, the causal link between those motifs and the solidification process remains unclear. Here, we apply a novel approach for exploring nontrivial interactions between structural units in d-sorbitol, a canonical example of a hydrogen-bonded organic glass, by introducing a dihedral-rearrangement-indicator analysis to shed light on relaxation processes and dynamical heterogeneity, which are known for their association with the stability of a glass. We find that both α- and β-relaxation processes are governed by cooperative and heterogeneous changes in hydrogen-bond dynamics that can be described by spatial and dihedral-angle-rearrangement indicators. The methodology and findings are of general applicability to other glass-forming systems.

Understanding the Origin of the Breakdown of the Stokes-Einstein Relation in Supercooled Water at Different Temperature-Pressure Conditions

A recent experiment has measured the viscosity of water down to approximately 244 K and up to 300 MPa. The correct viscosity and translational diffusivity data at various temperature-pressure (T-P) state points allowed for checking the validity of the Stokes-Einstein (SE) relation, which accounts for the coupling between translational self-diffusion and medium viscosity. The diffusion-viscosity decoupling increases with decreasing temperature, but the increasing pressure reduces the extent of the decoupling. Earlier simulation studies explained the breakdown of the SE relation in terms of the location of the Widom line, emanating from the liquid-liquid critical point (LLCP). Although these studies made a significant contribution to the current understanding of the above phenomena, a detailed molecular picture is still lacking. Recently, our group has explained the diffusion-viscosity decoupling from a jump-diffusion perspective. The jump-diffusion coefficient, emanating from the jump translation of water molecules, is calculated using a quantitative approach for different temperatures at ambient pressure. It has been observed that jump-diffusion is the key factor for diffusion-viscosity decoupling in supercooled water. The same method is adopted in the present work to estimate the jump-diffusion coefficient for different T-P state points and, thereby, explains the role of jump-diffusion for the different extents of the SE relation breakdown at different pressures. The residual diffusion coefficient, the other component of the total diffusion that originates from small step displacement and that is calculated by subtracting the jump-diffusion coefficient from the total diffusion, is seen to be fairly coupled to the viscosity at the entire range of temperature and pressure. Furthermore, we have calculated the average number of H-bonds per water molecule and the tetrahedral order for different T-P state points and investigated an approximate correlation between the average local structure and the contribution of the jump-diffusion to the total diffusion of water. This study, therefore, puts forward a new perspective for explaining the SE relation breakdown in supercooled water under different pressure conditions.

Molecular Dynamics Simulation of the Capillary Leveling of a Glass-Forming Liquid

Motivated by recent experimental studies probing (i) the existence of a mobile layer at the free surface of glasses and (ii) the capillary leveling of polymer nanofilms, we study the evolution of square-wave patterns at the free surface of a generic glass-forming binary Lennard-Jones mixture over a wide temperature range, by means of molecular dynamics simulations. The pattern's amplitude is monitored, and the associated decay rate is extracted. The evolution of the latter as a function of temperature exhibits a crossover between two distinct behaviors, over a temperature range typically bounded by the glass-transition temperature and the mode-coupling critical temperature. Layer-resolved analysis of the film particles' mean-squared displacements further shows that diffusion at the surface is considerably faster than in the bulk, below the glass-transition temperature. The diffusion coefficient of the surface particles is larger than its bulk counterpart by a factor that reaches 10⁵ at the lowest temperature studied. This factor decreases upon heating, in agreement with recent experimental studies.

A comparative study on the dynamic heterogeneity of supercooled polymers under nanoconfinement

Dynamic heterogeneity (DH) is a universal property of glass transition phenomena. In this work, we perform a comparative analysis of DH for pure polymer and polymer/nanoparticle composite systems in both film and bulk states via molecular dynamics simulations. We find that the dynamic gradient and the faster average dynamics due to the presence of a free surface are two leading factors, resulting from a nanoconfinement effect, which influence different parts of DH in a film system. The dynamic gradient results from differences in dynamics at different distances from the mobile surface, which induces a large deviation from the Gaussian distribution for the displacement distribution in the film. At the same time, the maximum string size which describes the region size for cooperative motion (dynamic correlation) can also be influenced by the dynamic gradient, although this influence is much weaker than that on the displacement distribution. On the other hand, reflecting temporal fluctuations of dynamics or temporal parts of DH, characteristic peak times of the non-Gaussian parameter and string size, and the ratio between persistent times and exchange times which describe the dynamic exchange properties, are mainly influenced by the faster dynamics on average. Our results demonstrate that measuring different properties (dynamic distribution, dynamic correlation or dynamic exchange) place an emphasis on distinct temporal and spatial parts of DH. It is necessary to use combinational measurements of these properties to give a complete picture of DH in nanoconfinement environments.

Evidence of a two-dimensional glass transition in graphene: Insights from molecular simulations

Liquids exhibit a sudden increase in viscosity when cooled fast enough, avoiding thermodynamically predicted route of crystallization. This phenomenon, known as glass transition, leads to the formation of non-periodic structures known as glasses. Extensive studies have been conducted on model materials to understand glass transition in two dimensions. However, despite the synthesis of disordered/amorphous single-atom thick structures of carbon, little attention has been given to glass transition in realistic two-dimensional materials such as graphene. Herein, using molecular dynamics simulation, we demonstrate the existence of glass transition in graphene leading to a realistic two-dimensional glassy structure, namely glassy graphene. We show that the resulting glassy structure exhibits excellent agreement with experimentally realized disordered graphene. Interestingly, this glassy graphene exhibits a wrinkled but stable structure, with reduced thermal vibration in comparison to its crystalline counterpart. We suggest that the topological disorder induced by glass transition governs the unique properties of this structure.

Structural Anomalies in Ionic Liquids near the Glass Transition Revealed by Pulse EPR

Unusual physical and chemical properties of ionic liquids (ILs) open up prospects for various applications. We report the first observation of density/rigidity heterogeneities in a series of ILs near the glass transition temperature ( T_g) by means of pulse electron paramagnetic resonance (EPR). Unprecedented suppression of molecular mobility is evidenced near the glass transition, which is assigned to unusual structural rearrangements of ILs on the nanometer scale. Indeed, pulse and continuous wave EPR clearly indicate the occurrence of heterogeneities near T_g, which exist in a rather broad temperature range of ∼50 K. The two types of local environments are evidenced, being drastically different by their stiffness. The more rigid one suppresses molecular mobility, whereas the softer one instead promotes diffusive molecular rotation. Such properties of ILs near T_g are of general importance; moreover, the observed density/rigidity heterogeneities controlled by temperature might be considered as a new type of tunable reaction nanoenvironment.

Elucidation of transport mechanism and enhanced alkali ion transference numbers in mixed alkali metal–organic ionic molten salts

This work reveals how structure facilitates diffusion of the Li/Na ion in ionic liquids with the high Li/Na concentration.

Influence of molecular-weight polydispersity on the glass transition of polymers

It is well known that the polymer glass transition temperature T_{g} is dependent on molecular weight, but the role of molecular-weight polydispersity on T_{g} is unclear. Using molecular-dynamics simulations, we clarify that for polymers with the same number-average molecular weight, the molecular-weight distribution profile (either in Schulz-Zimm form or in bimodal form) has very little influence on the glass transition temperature T_{g}, the average segment dynamics (monomer motion, bond orientation relaxation, and torsion transition), and the relaxation-time spectrum, which are related to the local nature of the glass transition. By analyzing monomer motions in different chains, we find that the motion distribution of monomers is altered by molecular-weight polydispersity. Molecular-weight polydispersity dramatically enhances the dynamic heterogeneity of monomer diffusive motions after breaking out of the "cage," but it has a weak influence on the dynamic heterogeneity of the short time scales and the transient spatial correlation between temporarily localized monomers. The stringlike cooperative motion is also not influenced by molecular-weight polydispersity, supporting the idea that stringlike collective motion is not strongly correlated with chain connectivity.

Variation of the dynamic susceptibility along an isochrone

Koperwas et al. showed in a recent paper [Phys. Rev. Lett. 111, 125701 (2013)] that the dynamic susceptibility χ4 as estimated by dielectric measurements for certain glass-forming liquids decreases substantially with increasing pressure along a curve of constant relaxation time. This observation is at odds with other measures of dynamics being invariant and seems to pose a problem for theories of glass formation. We show that this variation is in fact consistent with predictions for liquids with hidden scale invariance: Measures of dynamics at constant volume are invariant along isochrones, called isomorphs in such liquids, but contributions to fluctuations from long-wavelength fluctuations can vary. This is related to the known noninvariance of the isothermal bulk modulus. Considering the version of χ4 defined for the NVT ensemble, data from simulations of a binary Lennard-Jones liquid show in fact a slight increase with increasing density. This is a true departure from the formal invariance expected for this quantity.

Divergent dynamics and the Kauzmann temperature in glass forming systems

In the last decade the challenging analysis of previtreous behavior of relaxation time (τ(T)) in ultraviscous low molecular weight liquids led to the conceptual shift of the glass transition physics toward theories not predicting a “finite-temperature” divergence. This “breakthrough” experimental finding was strengthened by the discovery that “dynamic” (i.e. from τ(T) fitting) and “thermodynamic” estimations of the “ideal glass” (Kauzmann) temperature do not match, what in fact questioned its existence. In this report, due to the novel way of analysis based on the transformation of τ(T) experimental data to the activation energy temperature index form, the clear prevalence of the “finite-temperature” divergence is proved. The obtained “dynamic” singular temperatures clearly coincide with “thermodynamic” estimations of the Kauzmann temperature, thus solving also the second mystery. The comprehensive picture was obtained due to the analysis of 55 experimental data-sets, ranging from low molecular weight liquids and polymers to liquid crystal and plastic crystals.

Perspective: The glass transition

We provide here a brief perspective on the glass transition field. It is an assessment, written from the point of view of theory, of where the field is and where it seems to be heading. We first give an overview of the main phenomenological characteristics, or "stylised facts," of the glass transition problem, i.e., the central observations that a theory of the physics of glass formation should aim to explain in a unified manner. We describe recent developments, with a particular focus on real space properties, including dynamical heterogeneity and facilitation, the search for underlying spatial or structural correlations, and the relation between the thermal glass transition and athermal jamming. We then discuss briefly how competing theories of the glass transition have adapted and evolved to account for such real space issues. We consider in detail two conceptual and methodological approaches put forward recently, that aim to access the fundamental critical phenomenon underlying the glass transition, be it thermodynamic or dynamic in origin, by means of biasing of ensembles, of configurations in the thermodynamic case, or of trajectories in the dynamic case. We end with a short outlook.

A universal description of ultraslow glass dynamics

The dynamics of glass is of importance in materials science but its nature has not yet been fully understood. Here we report that a verification of the temperature dependencies of the primary relaxation time or viscosity in the ultraslowing/ultraviscous domain of glass-forming systems can be carried out via the analysis of the inverse of the Dyre-Olsen temperature index. The subsequent analysis of experimental data indicates the possibility of the self-consistent description of glass-forming low-molecular-weight liquids, polymers, liquid crystals, orientationally disordered crystals and Ising spin-glass-like systems, as well as the prevalence of equations associated with the 'finite temperature divergence'. All these lead to a new formula for the configurational entropy in glass-forming systems. Furthermore, a link to the dominated local symmetry for a given glass former is identified here. Results obtained show a new relationship between the glass transition and critical phenomena.

Bond orientational order in liquids: Towards a unified description of water-like anomalies, liquid-liquid transition, glass transition, and crystallization: Bond orientational order in liquids

There are at least three fundamental states of matter, depending upon temperature and pressure: gas, liquid, and solid (crystal). These states are separated by first-order phase transitions between them. In both gas and liquid phases a complete translational and rotational symmetry exist, whereas in a solid phase both symmetries are broken. In intermediate phases between liquid and solid, which include liquid crystal and plastic crystal phases, only one of the two symmetries is preserved. Among the fundamental states of matter, the liquid state is the most poorly understood. We argue that it is crucial for a better understanding of liquids to recognize that a liquid generally has the tendency to have a local structural order and its presence is intrinsic and universal to any liquid. Such structural ordering is a consequence of many-body correlations, more specifically, bond angle correlations, which we believe are crucial for the description of the liquid state. We show that this physical picture may naturally explain difficult unsolved problems associated with the liquid state, such as anomalies of water-type liquids (water, Si, Ge, ...), liquid-liquid transition, liquid-glass transition, crystallization and quasicrystal formation, in a unified manner. In other words, we need a new order parameter representing a low local free-energy configuration, which is a bond orientational order parameter in many cases, in addition to a density order parameter for the physical description of these phenomena. Here we review our two-order-parameter model of liquid and consider how transient local structural ordering is linked to all of the above-mentioned phenomena. The relationship between these phenomena is also discussed.

Driving kinetically constrained models into nonequilibrium steady states: Structural and slow transport properties

Complex fluids in shear flow and biased dynamics in crowded environments exhibit counterintuitive features which are difficult to address both at a theoretical level and by molecular dynamic simulations. To understand some of these features we study a schematic model of a highly viscous liquid, the two-dimensional Kob-Andersen kinetically constrained model, driven into nonequilibrium steady states by a uniform non-Hamiltonian force. We present a detailed numerical analysis of the microscopic behavior of the model, including transversal and longitudinal spatial correlations and dynamic heterogeneities. In particular, we show that at high particle density the transition from positive to negative resistance regimes in the current vs field relation can be explained via the emergence of nontrivial structures that intermittently trap the particles and slow down the dynamics. We relate such spatial structures to the current vs field relation in the different transport regimes.

Finite-size effects in the dynamics of glass-forming liquids

We present a comprehensive theoretical study of finite-size effects in the relaxation dynamics of glass-forming liquids. Our analysis is motivated by recent theoretical progress regarding the understanding of relevant correlation length scales in liquids approaching the glass transition. We obtain predictions both from general theoretical arguments and from a variety of specific perspectives: mode-coupling theory, kinetically constrained and defect models, and random first-order transition theory. In the last approach, we predict in particular a nonmonotonic evolution of finite-size effects across the mode-coupling crossover due to the competition between mode-coupling and activated relaxation. We study the role of competing relaxation mechanisms in giving rise to nonmonotonic finite-size effects by devising a kinetically constrained model where the proximity to the mode-coupling singularity can be continuously tuned by changing the lattice topology. We use our theoretical findings to interpret the results of extensive molecular dynamics studies of four model liquids with distinct structures and kinetic fragilities. While the less fragile model only displays modest finite-size effects, we find a more significant size dependence evolving with temperature for more fragile models, such as Lennard-Jones particles and soft spheres. Finally, for a binary mixture of harmonic spheres we observe the predicted nonmonotonic temperature evolution of finite-size effects near the fitted mode-coupling singularity, suggesting that the crossover from mode-coupling to activated dynamics is more pronounced for this model. Finally, we discuss the close connection between our results and the recent report of a nonmonotonic temperature evolution of a dynamic length scale near the mode-coupling crossover in harmonic spheres.

Facilitated spin models of dissipative quantum glasses

We introduce a class of dissipative quantum spin models with local interactions and without quenched disorder that show glassy behavior. These models are the quantum analogs of the classical facilitated spin models. Just like their classical counterparts, quantum facilitated models display complex glassy dynamics despite the fact that their stationary state is essentially trivial. In these systems, dynamical arrest is a consequence of kinetic constraints and not of static ordering. These models display a quantum version of dynamic heterogeneity: the dynamics toward relaxation is spatially correlated despite the absence of static correlations. Associated dynamical fluctuation phenomena such as decoupling of time scales is also observed. Moreover, we find that close to the classical limit, quantum fluctuations can enhance glassiness, as recently reported for quantum liquids.

The physics of the colloidal glass transition

As one increases the concentration of a colloidal suspension, the system exhibits a dramatic increase in viscosity. Beyond a certain concentration, the system is said to be a colloidal glass; structurally, the system resembles a liquid, yet motions within the suspension are slow enough that it can be considered essentially frozen. For several decades, colloids have served as a valuable model system for understanding the glass transition in molecular systems. The spatial and temporal scales involved allow these systems to be studied by a wide variety of experimental techniques. The focus of this review is the current state of understanding of the colloidal glass transition, with an emphasis on experimental observations. A brief introduction is given to important experimental techniques used to study the glass transition in colloids. We describe features of colloidal systems near and in glassy states, including increases in viscosity and relaxation times, dynamical heterogeneity and ageing, among others. We also compare and contrast the glass transition in colloids to that in molecular liquids. Other glassy systems are briefly discussed, as well as recently developed synthesis techniques that will keep these systems rich with interesting physics for years to come.

Time dependent stretching of aging dynamics in a generalized hydrodynamic model for supercooled liquids

The nonequilibrium dynamics and aging behavior of a supercooled liquid is investigated from an analysis of the correlation of density fluctuations at two different times. The dynamic correlation functions are computed by solving numerically the equations of nonlinear fluctuating hydrodynamics. The aging time dependence follows a modified stretched exponential form with a relaxation time which is dependent on the aging time. This is similar to the behavior seen in the aging data of dielectric response functions of a typical glass forming liquid.

Cooperative rearrangement regions and dynamical heterogeneities in colloidal glasses with attractive versus repulsive interactions

Water-lutidine mixtures permit the interparticle potentials of colloidal particles suspended therein to be tuned, in situ, from repulsive to attractive. We employ these systems to directly elucidate the effects of interparticle potential on glass dynamics in experimental samples composed of the same particles at the same packing fractions. Cooperative rearrangement regions (CRRs) and heterogeneous dynamics are observed in both types of glasses. Compared to repulsive glasses, the attractive glass dynamics are found to be heterogeneous over a wider range of time and length scales, and its CRRs involve more particles. Additionally, the CRRs are observed to be stringlike structures in repulsive glasses and compact structures in attractive glasses. Thus, the experiments demonstrate explicitly that glassy dynamics can depend on the sign of the interparticle interaction.