Dynamic phase coexistence in glass-forming liquids

Comparing Microscopic and Macroscopic Dynamics in a Paradigmatic Model of Glass-Forming Molecular Liquid

Glass transition is a most intriguing and long-standing open issue in the field of molecular liquids. From a macroscopic perspective, glass-forming systems display a dramatic slowing-down of the dynamics, with the inverse diffusion coefficient and the structural relaxation times increasing by orders of magnitude upon even modest supercooling. At the microscopic level, single-molecule motion becomes strongly intermittent, and can be conveniently described in terms of “cage-jump” events. In this work, we investigate a paradigmatic glass-forming liquid, the Kob–Andersen Lennard–Jones model, by means of Molecular Dynamics simulations, and compare the macroscopic and microscopic descriptions of its dynamics on approaching the glass-transition. We find that clear changes in the relations between macroscopic timescales and cage-jump quantities occur at the crossover temperature where Mode Coupling-like description starts failing. In fact, Continuous Time Random Walk and lattice model predictions based on cage-jump statistics are also violated below the crossover temperature, suggesting the onset of a qualitative change in cage-jump motion. Interestingly, we show that a fully microscopic relation linking cage-jump time- and length-scales instead holds throughout the investigated temperature range.

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.

Biological physics by high-speed atomic force microscopy

While many fields have contributed to biological physics, nanotechnology offers a new scale of observation. High-speed atomic force microscopy (HS-AFM) provides nanometre structural information and dynamics with subsecond resolution of biological systems. Moreover, HS-AFM allows us to measure piconewton forces within microseconds giving access to unexplored, fast biophysical processes. Thus, HS-AFM provides a tool to nourish biological physics through the observation of emergent physical phenomena in biological systems. In this review, we present an overview of the contribution of HS-AFM, both in imaging and force spectroscopy modes, to the field of biological physics. We focus on examples in which HS-AFM observations on membrane remodelling, molecular motors or the unfolding of proteins have stimulated the development of novel theories or the emergence of new concepts. We finally provide expected applications and developments of HS-AFM that we believe will continue contributing to our understanding of nature, by serving to the dialogue between biology and physics. This article is part of a discussion meeting issue 'Dynamic in situ microscopy relating structure and function'.

Breakdown of the Stokes-Einstein Equation for Solutions of Water in Oil Reverse Micelles

An experimental study is presented for the reverse micellar system of 15% by mass polydisperse hexaethylene glycol monodecylether (C₁₀E₆) in cyclohexane with varying amounts of added water up to 4% by mass. Measurements of viscosity and self-diffusion coefficients were taken as a function of temperature between 10 and 45 °C at varying sample water loads but fixed C₁₀E₆/cyclohexane composition. The results were used to inspect the validity of the Stokes-Einstein equation for this system. Unreasonably small reverse average micelle radii and aggregation numbers were obtained with the Stokes-Einstein equation, but reasonable values for these quantities were obtained using the ratio of surfactant-to-cyclohexane self-diffusion coefficients. While bulk viscosity increased with increasing water load, a concurrent expected decrease of self-diffusion coefficient was only observed for the surfactant and water but not for cyclohexane, which showed independence of water load. Moreover, a spread of self-diffusion coefficients was observed for the protons associated with the ethylene oxide repeat unit in samples with polydisperse C₁₀E₆ but not in a sample with monodisperse C₁₀E₆. These findings were interpreted by the presence of reverse micelle to reverse micelle hopping motions that with higher water load become increasingly selective toward C₁₀E₆ molecules with short ethylene oxide repeat units, while those with long ethylene oxide repeat units remain trapped within the reverse micelle because of the increased hydrogen bonding interactions with the water inside the growing core of the reverse micelle. Despite the observed breakdown of the Stokes-Einstein equation, the temperature dependence of the viscosities and self-diffusion coefficients was found to follow Arrhenius behavior over the investigated range of temperatures.

Connection of large amplitude angular jump motions with temporal heterogeneity in aqueous solutions

It has now been established that large angular jumps do take place when a rotating water molecule exchanges its hydrogen bond (H-bond) identity. This motion differs from the small angular diffusional steps occurring within short time intervals which define the 'Debye diffusion model' of water dynamics. We intend to investigate whether these two processes do eventually complement each other. In this present investigation the orientational dynamics of water in its mixture with a small hydrophobic molecule 1,2-dimethoxy ethane (DME) is studied microscopically using the all-atom classical molecular dynamics (MD) simulation technique. We found that the reorientational motions of water molecules are governed by continuous making and breaking of intermolecular H-bonds with their partners. We characterise these H-bond reorientation motions with the description of the "large amplitude angular jump model" and explore the coupling between the rotational and translational motions. By following the trajectories of each molecule in the solutions we describe the orientational dynamics of liquid water with a 'continuous time random walk' (CTRW) approach. Finally, we explore the diffusivity distribution through the jump properties of the water molecules, which successfully leads to the inherent transient heterogeneity of the solutions. We observe that the heterogeneity increases with increasing DME content in the mixtures. Our study correlates the coupling between rotational and translational motions of water molecules in the mixtures.

A proposal for the structure of high- and low-density fluctuations in liquid water

Based on recent experimental data that can be interpreted as indicating the presence of specific structures in liquid water, we build and optimize two structural models which we compare with the available experimental data. To represent the proposed high-density liquid structures, we use a model consisting of chains of water molecules, and for low-density liquid, we investigate fused dodecahedra as templates for tetrahedral fluctuations. The computed infrared spectra of the models are in very good agreement with the extracted experimental spectra for the two components, while the extracted structures from molecular dynamics (MD) simulations give spectra that are intermediate between the experimentally derived spectra. Computed x-ray absorption and emission spectra as well as the O-O radial distribution functions of the proposed structures are not contradicted by experiment. The stability of the proposed dodecahedral template structures is investigated in MD simulations by seeding the starting structure, and remnants found to persist on an ∼30 ps time scale. We discuss the possible significance of such seeds in simulations and whether they can be viable candidates as templates for structural fluctuations below the compressibility minimum of liquid water.

Glassy dynamics predicted by mutual role of free and activation volumes

Broadband Dielectric Spectroscopy (BDS) at elevated pressures and Positron Annihilation Lifetime Spectroscopy (PALS) are employed to elucidate the importance of the ratio of activation and free volumes during vitrification. We show that this ratio has a linear correlation with the structural relaxation of glass forming liquids in a wide temperature range hence engendering it as a vital input in the description of the dynamic glass transition.

Statistics of small length scale density fluctuations in supercooled viscous liquids

Many successful theories of liquids near the melting temperature assume that small length scale density fluctuations follow Gaussian statistics. This paper presents a numerical investigation of density fluctuations in the supercooled viscous regime using an enhanced sampling method. Five model systems are investigated: the single component Lennard-Jones liquid, the Kob-Andersen binary mixture, the Wahnström binary mixture, the Lewis-Wahnström model of ortho-terphenyl, and the TIP4P/Ice model of water. The results show that the Gaussian approximation persists to a good degree into the supercooled viscous regime; however, it is less accurate at low temperatures. The analysis suggests that non-Gaussian fluctuations are related to crystalline configurations. Implications for theories of the glass transition are discussed.

Heterogeneous Activation, Local Structure, and Softness in Supercooled Colloidal Liquids

We experimentally characterize heterogeneous nonexponential relaxation in bidisperse supercooled colloidal liquids utilizing a recent concept called "softness" [Phys. Rev. Lett. 114, 108001 (2015)PRLTAO0031-900710.1103/PhysRevLett.114.108001]. Particle trajectory and structure data enable classification of particles into subgroups with different local environments and propensities to hop. We determine residence times t_{R} between particle hops and show that t_{R} derived from particles in the same softness subgroup are exponentially distributed. Using the mean residence time t[over ¯]_{R} for each softness subgroup, and a Kramers' reaction rate model, we estimate the activation energy barriers E_{b} for particle hops, and show that both t[over ¯]_{R} and E_{b} are monotonic functions of softness. Finally, we derive information about the combinations of large and small particle neighbors that determine particle softness, and we explicitly show that multiple exponential relaxation channels in the supercooled liquid give rise to its nonexponential behavior.

Non-Gaussianity of the van Hove function and dynamic-heterogeneity length scale

Non-Gaussian nature of the probability distribution of particles' displacements in the supercooled temperature regime in glass-forming liquids are believed to be one of the major hallmarks of glass transition. It has already been established that this probability distribution, which is also known as the van Hove function, shows universal exponential tail. The origin of such an exponential tail in the distribution function is attributed to the hopping motion of particles observed in the supercooled regime. The non-Gaussian nature can also be explained if one assumes that the system has heterogeneous dynamics in space and time. Thus exponential tail is the manifestation of dynamic heterogeneity. In this work we directly show that non-Gaussianity of the distribution of particles' displacements occur over the dynamic heterogeneity length scale and the dynamical behavior course grained over this length scale becomes homogeneous. We study the non-Gaussianity of the van Hove function by systematically coarse graining at different length scales and extract the length scale of dynamic heterogeneity at which the shape of the van Hove function crosses over from non-Gaussian to Gaussian. The obtained dynamic heterogeneity scale is found to be in very good agreement with the scale obtained from other conventional methods.

Single molecule translocation in smectics illustrates the challenge for time-mapping in simulations on multiple scales

Understanding the connections between the characteristic dynamical time scales associated with a coarse-grained (CG) and a detailed representation is central to the applicability of the coarse-graining methods to understand molecular processes. The process of coarse graining leads to an accelerated dynamics, owing to the smoothening of the underlying free-energy landscapes. Often a single time-mapping factor is used to relate the time scales associated with the two representations. We critically examine this idea using a model system ideally suited for this purpose. Single molecular transport properties are studied via molecular dynamics simulations of the CG and atomistic representations of a liquid crystalline, azobenzene containing mesogen, simulated in the smectic and the isotropic phases. The out-of-plane dynamics in the smectic phase occurs via molecular hops from one smectic layer to the next. Hopping can occur via two mechanisms, with and without significant reorientation. The out-of-plane transport can be understood as a superposition of two (one associated with each mode of transport) independent continuous time random walks for which a single time-mapping factor would be rather inadequate. A comparison of the free-energy surfaces, relevant to the out-of-plane transport, qualitatively supports the above observations. Thus, this work underlines the need for building CG models that exhibit both structural and dynamical consistency to the underlying atomistic model.

Structural-dynamical transition in the Wahnström mixture

In trajectory space, dynamical heterogeneities in glass-forming liquids correspond to the emergence of a dynamical phase transition between an active phase poor in local structure and an inactive phase which is rich in local structure. We support this scenario with the study of a model additive mixture of Lennard-Jones particles, quantifying how the choice of the relevant structural and dynamical observable affects the transition in trajectory space. We find that the low mobility, structure-rich phase is dominated by icosahedral order. Applying a non-equilibrium rheological protocol, we connect local order to the emergence of mechanical rigidity.

From Nanoscale to Microscale: Crossover in the Diffusion Dynamics within Two Pyrrolidinium-Based Ionic Liquids

Knowledge of the ion motion in room temperature ionic liquids (RTILs) is critical for their applications in a number of fields, from lithium batteries to dye-sensitized solar cells. Experiments on a limited number of RTILs have shown that on macroscopic time scales the ions typically undergo conventional, Gaussian diffusion. On shorter time scales, however, non-Gaussian behavior has been observed, similar to supercooled fluids, concentrated colloidal suspensions, and more complex systems. Here we characterize the diffusive motion of ionic liquids based on the N-butyl-N-methylpyrrolidinium (PYR₁₄) cation and bis(trifluoro methanesulfonyl)imide (TFSI) or bis(fluorosulfonyl)imide (FSI) anions. A combination of pulsed gradient spin-echo (PGSE) NMR experiments and molecular dynamics (MD) simulations demonstrates a crossover from subdiffusive behavior to conventional Gaussian diffusion at ∼10 ns. The deconvolution of molecular displacements into a continuous spectrum of diffusivities shows that the short-time behavior is related to the effects of molecular caging. For PYR₁₄FSI, we identify the change of short-range ion-counterion associations as one possible mechanism triggering long-range displacements.

Length-Scale Dependence of the Stokes-Einstein and Adam-Gibbs Relations in Model Glass Formers

The Adam-Gibbs (AG) relation connects the dynamics of a glass-forming liquid to its thermodynamics via the configurational entropy and is of fundamental importance in descriptions of glassy behavior. The breakdown of the Stokes-Einstein relation (SEB) between the diffusion coefficient and the viscosity (or structural relaxation times) in glass formers raises the question as to which dynamical quantity the AG relation describes. By performing molecular dynamics simulations, we show that the AG relation is valid over the widest temperature range for the diffusion coefficient and not for the viscosity or relaxation times. Studying relaxation times defined at a given wavelength, we find that SEB and the deviation from the AG relation occur below a temperature at which the correlation length of dynamical heterogeneity equals the wavelength probed.

Accelerated crystallization of colloidal glass by mechanical oscillation

Crystallization of a hard-sphere colloidal glass by mechanical oscillation is investigated, and accelerated crystallization is found at a specific frequency. The crystallization frequency increases as attractive force between particles increases, indicating that interparticle interaction affects the crystallization frequency. Time scale of the mechanical oscillation is different from that of the slow relaxation, and notable relationship with the low-frequency mode is not observed. The experimental results are not explained by the previously proposed model for crystallization by oscillatory shear. Conversely, we speculate that activations of the fast relaxation and particle motion in crystal-like clusters are possible causes of the observations.

Cage Size and Jump Precursors in Glass-Forming Liquids: Experiment and Simulations

Glassy dynamics is intermittent, as particles suddenly jump out of the cage formed by their neighbors, and heterogeneous, as these jumps are not uniformly distributed across the system. Relating these features of the dynamics to the diverse local environments explored by the particles is essential to rationalize the relaxation process. Here we investigate this issue characterizing the local environment of a particle with the amplitude of its short time vibrational motion, as determined by segmenting in cages and jumps the particle trajectories. Both simulations of supercooled liquids and experiments on colloidal suspensions show that particles in large cages are likely to jump after a small time-lag, and that, on average, the cage enlarges shortly before the particle jumps. At large time-lags, the cage has essentially a constant size, which is smaller for longer-lasting cages. Finally, we clarify how this coupling between cage size and duration controls the average behavior and opens the way to a better understanding of the relaxation process in glass-forming liquids.

Differential Variance Analysis: a direct method to quantify and visualize dynamic heterogeneities

Many amorphous materials show spatially heterogenous dynamics, as different regions of the same system relax at different rates. Such a signature, known as Dynamic Heterogeneity, has been crucial to understand the nature of the jamming transition in simple model systems and is currently considered very promising to characterize more complex fluids of industrial and biological relevance. Unfortunately, measurements of dynamic heterogeneities typically require sophisticated experimental set-ups and are performed by few specialized groups. It is now possible to quantitatively characterize the relaxation process and the emergence of dynamic heterogeneities using a straightforward method, here validated on video microscopy data of hard-sphere colloidal glasses. We call this method Differential Variance Analysis (DVA), since it focuses on the variance of the differential frames, obtained subtracting images at different time-lags. Moreover, direct visualization of dynamic heterogeneities naturally appears in the differential frames, when the time-lag is set to the one corresponding to the maximum dynamic susceptibility. This approach opens the way to effectively characterize and tailor a wide variety of soft materials, from complex formulated products to biological tissues.

Hydrodynamics in kinetically constrained lattice-gas models

Kinetically constrained models are lattice-gas models that are used for describing glassy systems. By construction, their equilibrium state is trivial and there are no equal-time correlations between the occupancy of different sites. We drive such models out of equilibrium by connecting them to two reservoirs of different densities, and we measure the response of the system to this perturbation. We find that under the proper coarse-graining, the behavior of these models may be expressed by a nonlinear diffusion equation, with a model- and density-dependent diffusion coefficient. We find a simple approximation for the diffusion coefficient, and we show that the relatively mild discrepancy between the approximation and our numerical results arises due to non-negligible correlations that appear as the system is driven out of equilibrium, even when the density gradient is infinitesimally small. Similar correlations appear when such kinetically constrained models are driven out of equilibrium by applying a uniform external force. We suggest that these correlations are the reason for the same discrepancy between the approximate diffusion coefficient and the numerical results for a broader group of models-nongradient lattice-gas models-for which kinetically constrained models are arguably the simplest example thereof.

Escape rate and diffusion of a Stochastically Driven particle

The dynamical properties of a tracer repeatedly colliding with heat bath particles can be described within a Langevin framework provided that the tracer is more massive than the bath particles, and that the collisions are frequent. Here we consider the escape of a particle from a potential well, and the diffusion coefficient in a periodic potential, without making these assumptions. We have thus investigated the dynamical properties of a Stochastically Driven particle that moves under the influence of the confining potential in between successive collisions with the heat bath. In the overdamped limit, both the escape rate and the diffusion coefficient coincide with those of a Langevin particle. Conversely, in the underdamped limit the two dynamics have a different temperature dependence. In particular, at low temperature the Stochastically Driven particle has a smaller escape rate, but a larger diffusion coefficient.

Single-particle dynamics near the glass transition of a metallic glass

The single-particle dynamics of the glass-forming Cu_{50}Zr_{50} alloy, from the supercooled liquid well above the glass-transition temperature, T_{g} to the glassy state, is studied by using the molecular dynamics simulations. When the liquid is cooled below 1.2T_{g}, the dynamics heterogeneity characterized by the cage-jump motion becomes increasingly pronounced. The analyses based on the continuous time random walk method indicate that the liquid falls out of equilibrium in the present simulation time scale when it is cooled into the regime below 1.02T_{g}. However, we find that the jump length and the jump rate do not display the non-equilibrium behaviors even in the glassy state below T_{g}, which allows us to study the intrinsic dynamic characteristics through T_{g}. The mean waiting time between two successive jumps has a rapid growth following the Vogel-Fulcher-Tammann law as the non-equilibrium regime is approached, in analogy with the temperature behaviors of transport properties for fragile supercooled liquids. In contrast, the jump rate maintains the Arrhenius decay and the jump length has even a weaker temperature dependence when the liquid is cooled into glassy state. We find that a pronounced enhancement of the spatial correlation of jumps occurs accompanied by the glass transition: the string-like cooperative jumps dominate the fast motion instead of the uncorrelated and individual jumps. Our work offers an insight into the equilibrium effect of the single-particle dynamics in glass transition.

Decoupling of viscosity and relaxation processes in supercooled water: a molecular dynamics study with the TIP4P/2005f model

We show that the TIP4P/2005f water model describes accurately the experimental viscosity and self-diffusion over a large temperature range. We then show the decoupling of viscosity and structural relaxation time in supercooled water.

Active cage model of glassy dynamics

We build up a phenomenological picture in terms of the effective dynamics of a tracer confined in a cage experiencing random hops to capture some characteristics of glassy systems. This minimal description exhibits scale invariance properties for the small-displacement distribution that echo experimental observations. We predict the existence of exponential tails as a crossover between two Gaussian regimes. Moreover, we demonstrate that the onset of glassy behavior is controlled only by two dimensionless numbers: the number of hops occurring during the relaxation of the particle within a local cage and the ratio of the hopping length to the cage size.

Field-induced superdiffusion and dynamical heterogeneity

By analyzing two kinetically constrained models of supercooled liquids we show that the anomalous transport of a driven tracer observed in supercooled liquids is another facet of the phenomenon of dynamical heterogeneity. We focus on the Fredrickson-Andersen and the Bertin-Bouchaud-Lequeux models. By numerical simulations and analytical arguments we demonstrate that the violation of the Stokes-Einstein relation and the field-induced superdiffusion observed during a long preasymptotic regime have the same physical origin: while a fraction of probes do not move, others jump repeatedly because they are close to local mobile regions. The anomalous fluctuations observed out of equilibrium in the presence of a pulling force ε,σ_{x}^{2}(t)=〈x_{ε}^{2}(t)〉-〈x_{ε}(t)〉^{2}∼t^{3/2}, which are accompanied by the asymptotic decay α_{ε}(t)∼t^{-1/2} of the non-Gaussian parameter from nontrivial values to zero, are due to the splitting of the probes population in the two (mobile and immobile) groups and to dynamical correlations, a mechanism expected to happen generically in supercooled liquids.

The heterogeneous energy landscape expression of KWW relaxation

Here we show a heterogeneous energy landscape approach to describing the Kohlrausch-Williams-Watts (KWW) relaxation function. For a homogeneous dynamic process, the distribution of free energy landscape is first proposed, revealing the significance of rugged fluctuations. In view of the heterogeneous relaxation given in two dynamic phases and the transmission coefficient in a rate process, we obtain a general characteristic relaxation time distribution equation for the KWW function in a closed, analytic form. Analyses of numerical computation show excellent accuracy, both in time and frequency domains, in the convergent performance of the heterogeneous energy landscape expression and shunning the catastrophic truncations reported in the previous work. The stretched exponential β, closely associated to temperature and apparent correlation with one dynamic phase, reveals a threshold value of 1/2 defining different behavior of the probability density functions. Our work may contribute, for example, to in-depth comprehension of the dynamic mechanism of glass transition, which cannot be provided by existing approaches.

Particle jumps in structural glasses

Particles in structural glasses rattle around temporary equilibrium positions, that seldom change through a process which is much faster than the relaxation time, known as particle jump. Since the relaxation of the system is due to the accumulation of many such jumps, it could be possible to connect the single particle short time motion to the macroscopic relaxation by understanding the features of the jump dynamics. Here we review recent results in this research direction, clarifying the features of particle jumps that have been understood and those that are still under investigation, and examining the role of particle jumps in different theories of the glass transition.

Glassy dynamics of a polymer monolayer on a heterogeneous disordered substrate

We present molecular dynamics simulations of a polymer monolayer on randomly functionalized surfaces that are characterized by different fractions of weakly and strongly attractive sites. We show that the dynamics slow-down upon cooling resembles that of a strong glass-forming liquid. Indeed, the mean-square displacements show an increasingly lasting subdiffusive behaviour before the diffusive regime, with signs of Fickian yet not Gaussian diffusion, and the dynamic correlation functions exhibit a stretched exponential decay. The glassy dynamics of this relatively dilute system is dominated by the interaction of the polymer with the substrate and becomes more marked when the substrate composition is heterogeneous. Accordingly, the estimated glass transition temperature shows a non-monotic dependence on surface composition, in agreement with previous results for the activation energy and with an analysis of the potential energy landscape experienced by the polymer beads. Our findings are relevant to the description of polymer-surface adhesion and friction and the development of polymer nanocomposites with tailored structural and mechanical properties.

Spatial correlations of elementary relaxation events in glass-forming liquids

The dynamical facilitation scenario, by which localized relaxation events promote nearby relaxation events in an avalanche process, has been suggested as the key mechanism connecting the microscopic and the macroscopic dynamics of structural glasses. Here we investigate the statistical features of this process via numerical simulations of a model structural glass. First we show that the relaxation dynamics of the system occurs through particle jumps that are irreversible, and that cannot be decomposed in smaller irreversible events. Then we show that each jump does actually trigger an avalanche. The characteristics of this avalanche change upon cooling, suggesting that the relaxation dynamics crossovers from a noise dominated regime, where jumps do not trigger other relaxation events, to a regime dominated by the facilitation process, where a jump triggers more relaxation events.