Evidence of growing spatial correlations at the glass transition from nonlinear response experiments

A "hump" in the third-order dielectric response of a highly polar liquid: Now you see it, now you don't

We have measured the linear and nonlinear dielectric responses of S-methoxypropylene carbonate, a highly polar glass-former, for which it has been reported that the "hump," which is typical of third harmonic susceptibilities, disappears across a 5 K temperature change. To understand this unusual feature, we have measured the responses to high amplitude ac and dc electric fields at the fundamental frequency. The static limits of these results are entered into a model aimed at reproducing nonlinear dielectric susceptibility spectra using the concept of a fictive electric field. This model reproduces the "hump" in the third-harmonic response and its seeming disappearance. It is revealed that the "hump" is predominantly the result of reduced time constants, a consequence of the energy the sample absorbs from the electric field. At elevated temperatures, the "hump" only appears to vanish because its reduced amplitude submerges below the extraordinarily high level of polarization saturation of this liquid.

Prigogine-Defay ratio of glassy freezing scales with liquid fragility

A detailed study of published experimental data for a variety of materials on the incremental variation of heat capacity, thermal expansion, and compressibility at glassy freezing reveals a striking dependence of the Prigogine-Defay ratio R on the fragility index m. At high m, R approaches values of ∼1, the Ehrenfest expectation for second-order continuous phase transitions, while R reaches values >20 for low fragilities. We explain this correlation by the degree of separation of the glassy freezing temperature from a hidden phase transition into an ideal low-temperature glass.

Glassy dynamics in a glass-forming liquid: A first-principles study of toluene

The microscopic understanding of the dramatic increase in viscosity of liquids when cooled toward the glass transition is a major unresolved issue in condensed matter physics. Here, we use machine learning methods to accelerate molecular dynamics simulations with first-principles accuracy for the glass-former toluene. We show that the increase in viscosity is intimately linked to the increasing number of dynamically correlated molecules N^{*}. While certain hallmark features of glassy dynamics, like physical aging, are linked to N^{*} as well, others, like relaxation stretching, are not.

Soft pinning: Experimental validation of static correlations in supercooled molecular glass-forming liquids

Enormous enhancement in the viscosity of a liquid near its glass transition is a hallmark of glass transition. Within a class of theoretical frameworks, it is connected to growing many-body static correlations near the transition, often called "amorphous ordering." At the same time, some theories do not invoke the existence of such a static length scale in the problem. Thus, proving the existence and possible estimation of the static length scales of amorphous order in different glass-forming liquids is very important to validate or falsify the predictions of these theories and unravel the true physics of glass formation. Experiments on molecular glass-forming liquids become pivotal in this scenario as the viscosity grows several folds (∼ 10 14 ), and simulations or colloidal glass experiments fail to access these required long-time scales. Here we design an experiment to extract the static length scales in molecular liquids using dilute amounts of another large molecule as a pinning site. Results from dielectric relaxation experiments on supercooled Glycerol with different pinning concentrations of Sorbitol and Glucose, as well as the simulations on a few model glass-forming liquids with pinning sites, indicate the versatility of the proposed method, opening possible new avenues to study the physics of glass transition in other molecular liquids.

Nonlinear response theory for Markov processes. IV. The asymmetric double-well potential model revisited

The dielectric response of noninteracting dipoles is discussed in the framework of the classical model of stochastic reorientations in an asymmetric double-well potential (ADWP). In the nonlinear regime, this model exhibits some pecularities in the static response. We find that the saturation behavior of the symmetric double-well potential model does not follow the Langevin function and only in the linear regime are the standard results recovered. If a finite asymmetry is assumed, then the nonlinear susceptibilities are found to change the sign at a number of characteristic temperatures that depend on the magnitude of the asymmetry, as has been observed earlier for the third-order and fifth-order responses. If the kinetics of the barrier crossing in the ADWP model is described as a two-state model, then we can give analytical expressions for the values of the characteristic temperatures. The results for the response obtained from a (numerical) solution of the Fokker-Planck equation for the Brownian motion in a model ADWP behaves very similarly to the two-state model for high barriers. For small barriers no clear-cut timescale separation between the barrier crossing process and the intrawell relaxation exists and the model exhibits a number of timescales. In this case, the frequency-dependent linear susceptibility at low temperatures is dominated by the fast intrawell transitions and at higher temperatures by the barrier crossing kinetics. We find that for nonlinear susceptibilities the latter process appears to be more important and the intrawell transitions play only a role at the lowest temperatures.

Dynamical susceptibilities near ideal glass transitions

Building on the recently derived inhomogeneous mode-coupling theory, we extend the generalized mode-coupling theory of supercooled liquids to inhomogeneous environments. This provides a first-principles-based, systematic, and rigorous way of deriving high-point dynamical susceptibilities from variations of the many-body dynamic structure factors with respect to their conjugate field. This framework allows for a fully microscopic possibility to probe for collective relaxation mechanisms in supercooled liquids near the mode-coupling glass transition. The behavior of these dynamical susceptibilities is then studied in the context of simplified self-consistent relaxation models.

Temperature Dependence of Structural Relaxation in Glass-Forming Liquids and Polymers

Understanding the microscopic mechanism of the transition of glass remains one of the most challenging topics in Condensed Matter Physics. What controls the sharp slowing down of molecular motion upon approaching the glass transition temperature Tg, whether there is an underlying thermodynamic transition at some finite temperature below Tg, what the role of cooperativity and heterogeneity are, and many other questions continue to be topics of active discussions. This review focuses on the mechanisms that control the steepness of the temperature dependence of structural relaxation (fragility) in glass-forming liquids. We present a brief overview of the basic theoretical models and their experimental tests, analyzing their predictions for fragility and emphasizing the successes and failures of the models. Special attention is focused on the connection of fast dynamics on picosecond time scales to the behavior of structural relaxation on much longer time scales. A separate section discusses the specific case of polymeric glass-forming liquids, which usually have extremely high fragility. We emphasize the apparent difference between the glass transitions in polymers and small molecules. We also discuss the possible role of quantum effects in the glass transition of light molecules and highlight the recent discovery of the unusually low fragility of water. At the end, we formulate the major challenges and questions remaining in this field.

Glassy dynamics in polyalcohols: intermolecular simplicity vs. intramolecular complexity

Using depolarized light scattering, we have recently shown that structural relaxation in a broad range of supercooled liquids follows, to good approximation, a generic line shape with high-frequency power law ω^-1/2. We now continue this study by investigating a systematic series of polyalcohols (PAs), frequently used as model-systems in glass-science, i.a., because the width of their respective dielectric loss spectra varies strongly along the series. Our results reveal that the microscopic origin of the observed relaxation behavior varies significantly between different PAs: while short-chained PAs like glycerol rotate as more or less rigid entities and their light scattering spectra follow the generic shape, long-chained PAs like sorbitol display pronounced intramolecular dynamic contributions on the time scale of structural relaxation, leading to systematic deviations from the generic shape. Based on these findings we discuss an important limitation for observing the generic shape in a supercooled liquid: the dynamics that is probed needs to reflect the intermolecular dynamic heterogeneity, and must not be superimposed by effects of intramolecular dynamic heterogeneity.

Prolonged Association between Water Molecules under Hydrophobic Nanoconfinement

We present an investigation of the dynamics of water confined among rigid carbon rods and between parallel graphene sheets with molecular dynamics simulations. Diffusion coefficients, activation energy of diffusion, and residence-time correlation functions as a function of confinement geometry reveal a retardation of water dynamics under hydrophobic confinement compared to bulk water. In fact, water under various confinements possesses longer associations with its neighbors and exhibits diffusion dynamics characteristic of a lower temperature. Analysis of the residence-time correlation functions reveals long and short residence times, which we relate to the diffusion coefficient and activation energy of diffusion, respectively. Additional investigations reveal how the level of confining surface hydrophobicity affects water dynamics, further broadening our understanding of water diffusion inside diverse media. Overall, this study sheds light on the physical origin of retarded water dynamics under hydrophobic confinement and the close relationship between residence times and diffusion behavior.

Amorphous Order and Nonlinear Susceptibilities in Glassy Materials

We review 15 years of theoretical and experimental work on the nonlinear response of glassy systems. We argue that an anomalous growth of the peak value of nonlinear susceptibilities is a signature of growing "amorphous order" in the system, with spin-glasses as a case in point. Experimental results on supercooled liquids are fully compatible with the random first-order transition (RFOT) prediction of compact "glassites" of increasing volume as temperature is decreased, or as the system ages. We clarify why such a behavior is hard to explain within purely kinetic theories of glass formation, despite recent claims to the contrary.

Quantifying dielectric permittivities in the nonlinear regime

In contrast to the static dielectric permittivity,ε, associated with linear response, its high-field counterpart,ε_E, is not a material specific quantity, but rather depends on the experimental method used to determine the nonlinear dielectric effect (NDE). Here, we defineε_Ein a manner consistent with how high field permittivities are typically derived from a capacitance measurement using high voltages. Based upon characterizing the materials nonlinear behavior via its third order susceptibility,χ₃, the relations between a givenχ₃and the observableε_Eis calculated for six different experimental or theoretical approaches to NDEs in the static limit. It is argued that the quantityχ₃is superior overε_Eor the Piekara factor, (ε_E-ε)/E², because it facilitates an unambiguous comparison among different experimental techniques and it provides a more robust connection between experiment and theory.

Modeling non-linear dielectric susceptibilities of supercooled molecular liquids

Advances in high-precision dielectric spectroscopy have enabled access to non-linear susceptibilities of polar molecular liquids. The observed non-monotonic behavior has been claimed to provide strong support for theories of dynamic arrest based on the thermodynamic amorphous order. Here, we approach this question from the perspective of dynamic facilitation, an alternative view focusing on emergent kinetic constraints underlying the dynamic arrest of a liquid approaching its glass transition. We derive explicit expressions for the frequency-dependent higher-order dielectric susceptibilities exhibiting a non-monotonic shape, the height of which increases as temperature is lowered. We demonstrate excellent agreement with the experimental data for glycerol, challenging the idea that non-linear response functions reveal correlated relaxation in supercooled liquids.

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.

On the fluctuation-dissipation relation in non-equilibrium and non-Hamiltonian systems

We review generalized fluctuation-dissipation relations, which are valid under general conditions even in "nonstandard systems," e.g., out of equilibrium and/or without a Hamiltonian structure. The response functions can be expressed in terms of suitable correlation functions computed in the unperturbed dynamics. In these relations, typically, one has nontrivial contributions due to the form of the stationary probability distribution; such terms take into account the interaction among the relevant degrees of freedom in the system. We illustrate the general formalism with some examples in nonstandard cases, including driven granular media, systems with a multiscale structure, active matter, and systems showing anomalous diffusion.

Emergence and percolation of rigid domains during the colloidal glass transition

Using video microscopy, we measure local spatial constraints in disordered binary colloidal samples, ranging from dilute fluids to jammed glasses, and probe their spatial and temporal correlations to local dynamics during the glass transition. We observe the emergence of significant correlations between constraints and local dynamics within the Lindemann criterion, which coincides with the onset of glassy dynamics in supercooled liquids. Rigid domains in fluids are identified based on local constraints and demonstrate a percolation transition near the glass transition, accompanied by the emergence of dynamical heterogeneities. Our results show that spatial constraint instead of the geometry of amorphous structures is the key that connects the complex spatial-temporal correlations in disordered materials.

Ionic conductivity of deep eutectic solvents: the role of orientational dynamics and glassy freezing

Dielectric spectroscopy reveals that the ionic conductivity of deep eutectic solvents is closely coupled to their reorientational dipolar relaxation dynamics.

Transient Nonlinear Response of Dynamically Decoupled Ionic Conductors

The present study demonstrates that large electric fields progressively enhance the conductivity of ionic systems up to timescales corresponding to those on which their structural rearrangements take place. Yet, in many ionic materials, some regarded as candidates for electrical energy storage applications, the structural relaxation process can be tremendously slower than (or highly decoupled from) the charge fluctuations. Consequently, nonlinear dielectric spectroscopy may be employed to access rheological information in dynamically decoupled ionic conductors, whereas the combination of large electric power density and good mechanical stability, both technologically highly desired, imposes specific experimental constraints to reliably determine the steady-state conductivity of such materials.

Cooperativity in plastic crystals

A statistical mechanical model previously adopted for the analysis of the α-relaxation in structural glass formers is rederived within a general theoretical framework originally developed for systems approaching the ideal glassy state. The interplay between nonexponentiality and cooperativity is reconsidered in the light of energy landscape concepts. The method is used to estimate the cooperativity in orientationally disordered crystals, either from the analysis of literature data on linear dielectric response or from the enthalpy relaxation function obtained by temperature-modulated calorimetry. Knowledge of the specific heat step due to the freezing of the configurational or conformational modes at the glass transition is needed in order to properly account for the extent to which the relaxing system deviates from equilibrium during the rearrangement processes. A number of plastic crystals have been analyzed, and relatively higher cooperativities are found in the presence of hydrogen bonding interaction.

Unifying different interpretations of the nonlinear response in glass-forming liquids

This work aims at reconsidering several interpretations coexisting in the recent literature concerning nonlinear susceptibilities in supercooled liquids. We present experimental results on glycerol and propylene carbonate, showing that the three independent cubic susceptibilities have very similar frequency and temperature dependences, for both their amplitudes and phases. This strongly suggests a unique physical mechanism responsible for the growth of these nonlinear susceptibilities. We show that the framework proposed by two of us [J.-P. Bouchaud and G. Biroli, Phys. Rev. B 72, 064204 (2005)PRBMDO1098-012110.1103/PhysRevB.72.064204], where the growth of nonlinear susceptibilities is intimately related to the growth of glassy domains, accounts for all the salient experimental features. We then review several complementary and/or alternative models and show that the notion of cooperatively rearranging glassy domains is a key (implicit or explicit) ingredient to all of them. This paves the way for future experiments, which should deepen our understanding of glasses.

Pinning susceptibility: a novel method to study growth of amorphous order in glass-forming liquids

The existence and growth of amorphous order in supercooled liquids approaching glass transition is a subject of intense research. Even after decades of work, there is still no clear consensus on the molecular mechanisms that lead to a rapid slowing down of liquid dynamics approaching this putative transition. The existence of a correlation length associated with amorphous order has recently been postulated and has also been estimated using multi-point correlation functions which cannot be calculated easily in experiments. Thus the study of growing amorphous order remains mostly restricted to systems like colloidal glasses and simulations of model glass-forming liquids. In this work, we propose an experimentally realizable yet simple susceptibility to study the growth of amorphous order. We then demonstrate the validity of this approach for a few well-studied model supercooled liquids and obtain results which are consistent with other conventional methods.

Nonlinear dielectric effects in liquids: a guided tour

Dielectric relaxation measurements probe how the polarization of a material responds to the application of an external electric field, providing information on structure and dynamics of the sample. In the limit of small fields and thus linear response, such experiments reveal the properties of the material in the same thermodynamic state it would have in the absence of the external field. At sufficiently high fields, reversible changes in enthalpy and entropy of the system occur even at constant temperature, and these will in turn alter the polarization responses. The resulting nonlinear dielectric effects feature field induced suppressions (saturation) and enhancements (chemical effect) of the amplitudes, as well as time constant shifts towards faster (energy absorption) and slower (entropy reduction) dynamics. This review focuses on the effects of high electric fields that are reversible and observed at constant temperature for single component glass-forming liquids. The experimental challenges involved in nonlinear dielectric experiments, the approaches to separating and identifying the different sources of nonlinear behavior, and the current understanding of how high electric fields affect dielectric materials will be discussed. Covering studies from Debye's initial approach to the present state-of-the-art, it will be emphasized what insight can be gained from the nonlinear responses that are not available from dielectric relaxation results obtained in the linear regime.

Study of dynamical heterogeneities in colloidal nanoclay suspensions approaching dynamical arrest

The dynamics of aqueous Laponite clay suspensions slow down with increasing sample waiting time (t w ). This behavior, and the material fragility that results, closely resemble the dynamical slowdown in fragile supercooled liquids with decreasing temperature, and are typically ascribed to the increasing sizes of distinct dynamical heterogeneities in the sample. In this article, we characterize the dynamical heterogeneities in Laponite suspensions by invoking the three-point dynamic susceptibility formalism. The average time-dependent two-point intensity autocorrelation and its sensitivity to t w are probed in dynamic light scattering experiments. Distributions of relaxation time scales, deduced from the Kohlrausch-Williams-Watts equation, are seen to widen with increasing t w . The calculated three-point dynamic susceptibility of Laponite suspensions exhibits a peak, with the peak height increasing with evolving t w at fixed volume fraction or with increasing volume fraction at fixed t w , thereby signifying the slowdown of the sample dynamics. The number of dynamically correlated particles, calculated from the peak-height, is seen to initially increase rapidly with increasing t w , before eventually slowing down close to the non-ergodic transition point. This observation is in agreement with published reports on supercooled liquids and hard sphere colloidal suspensions and offers a unique insight into the colloidal glass transition of Laponite suspensions.

Nonlinear response theory for Markov processes. II. Fifth-order response functions

The nonlinear response of stochastic models obeying a master equation is calculated up to fifth order in the external field, thus extending the third-order results obtained earlier [G. Diezemann, Phys. Rev. E 85, 051502 (2012)PLEEE81539-375510.1103/PhysRevE.85.051502]. For sinusoidal fields the 5ω component of the susceptibility is computed for the model of dipole reorientations in an asymmetric double well potential and for a trap model with a Gaussian density of states. For most realizations of the models a hump is found in the higher-order susceptibilities. In particular, for the asymmetric double well potential model there are two characteristic temperature regimes showing the occurrence of such a hump as compared to a single characteristic regime in the case of the third-order response. In the case of the trap model the results strongly depend on the variable coupled to the field. As for the third-order response, the low-frequency limit of the susceptibility plays a crucial role with respect to the occurrence of a hump. The findings are discussed in light of recent experimental results obtained for supercooled liquids. The differences found for the third-order and the fifth-order response indicate that nonlinear response functions might serve as a powerful tool to discriminate among the large number of existing models for glassy relaxation.

Divergence of the third harmonic stress response to oscillatory strain approaching the glass transition

The leading nonlinear stress response in a periodically strained concentrated colloidal dispersion is studied experimentally and by theory. A thermosensitive microgel dispersion serves as well-characterized glass-forming model, where the stress response at the first higher harmonic frequency (3ω for strain at frequency ω) is investigated in the limit of small amplitude. The intrinsic nonlinearity at the third harmonic exhibits a scaling behavior which has a maximum in an intermediate frequency window and diverges when approaching the glass transition. It captures the (in-) stability of the transient elastic structure. Elastic stresses in-phase with the third power of the strain dominate the scaling. Our results qualitatively differ from previously derived scaling behavior in dielectric spectroscopy of supercooled molecular liquids. This might indicate a dependence of the nonlinear response on the symmetry of the external driving under time reversal.

External dc bias-field effects in the nonlinear ac stationary response of dipolar particles in a mean-field potential

External dc bias-field effects on the nonlinear dielectric relaxation and dynamic Kerr effect of a system of permanent dipoles in a uniaxial mean-field potential are studied via the rotational Brownian motion model postulated in terms of the infinite hierarchy of differential-recurrence equations for the statistical moments f_{n}(t)=〈P_{n}〉(t) (the expectation value of the Legendre polynomials P_{n}). By solving these equations, the nonlinear dielectric and Kerr-effect ac stationary responses are evaluated for arbitrary dc field strength via perturbation theory in the ac field. Simple analytic equations based on the large separation of the time scales of the fast intrawell and slow overbarrier (interwell) relaxation processes are also derived.

Length scales in glass-forming liquids and related systems: a review

The central problem in the study of glass-forming liquids and other glassy systems is the understanding of the complex structural relaxation and rapid growth of relaxation times seen on approaching the glass transition. A central conceptual question is whether one can identify one or more growing length scale(s) associated with this behavior. Given the diversity of molecular glass-formers and a vast body of experimental, computational and theoretical work addressing glassy behavior, a number of ideas and observations pertaining to growing length scales have been presented over the past few decades, but there is as yet no consensus view on this question. In this review, we will summarize the salient results and the state of our understanding of length scales associated with dynamical slow down. After a review of slow dynamics and the glass transition, pertinent theories of the glass transition will be summarized and a survey of ideas relating to length scales in glassy systems will be presented. A number of studies have focused on the emergence of preferred packing arrangements and discussed their role in glassy dynamics. More recently, a central object of attention has been the study of spatially correlated, heterogeneous dynamics and the associated length scale, studied in computer simulations and theoretical analysis such as inhomogeneous mode coupling theory. A number of static length scales have been proposed and studied recently, such as the mosaic length scale discussed in the random first-order transition theory and the related point-to-set correlation length. We will discuss these, elaborating on key results, along with a critical appraisal of the state of the art. Finally we will discuss length scales in driven soft matter, granular fluids and amorphous solids, and give a brief description of length scales in aging systems. Possible relations of these length scales with those in glass-forming liquids will be discussed.

Questioning the relationship between the χ4 susceptibility and the dynamical correlation length in a glass former

Clusters of fast and slow correlated particles, identified as dynamical heterogeneities (DHs), constitute a central aspect of glassy dynamics. A key factor of the glass transition scenario is a significant increase of the cluster size ξ4 as the transition is approached. In need of easy-to-compute tools to measure ξ4, the dynamical susceptibility χ4 was introduced recently, and used in various experimental studies to probe DHs. Here, we investigate DHs in dense microgel suspensions using image correlation analysis, and compute both χ4 and the four-point correlation function G4. The spatial decrease of G4 provides a direct access to ξ4, which is found to grow significantly with increasing volume fraction. However, this increase is not captured by χ4. We show that the assumptions that validate the connection between χ4 and ξ4 are not fulfilled in our experiments.

Field induced anisotropic cooperativity in a magnetic colloidal glass

The translational dynamics of a repulsive colloidal glass-former is probed by time-resolved X-ray Photon Correlation Spectroscopy. In this dense dispersion of charge-stabilized and magnetic nanoparticles, the interaction potential can be tuned, from quasi-isotropic to anisotropic by applying an external magnetic field. This powerful control parameter finely tunes the anisotropy of the intricate energy landscape in the colloidal glass-former, which is seen here as a new tunable model-system to probe the dynamical heterogeneities at the approach of the glass transition. Both structural and dynamical anisotropies are reported on interparticle lengthscales associated with highly anisotropic cooperativity, almost two orders of magnitude larger in the field direction than in the perpendicular direction and in zero field.

Unified Criterion for Temperature-Induced and Strain-Driven Glass Transitions in Metallic Glass

In a model metallic glass, we study the relaxation dynamics in both the linear and the nonlinear response regimes by numerical simulations of dynamical mechanical spectroscopy and analyze the atomic displacement statistics. We find that the primary (α) relaxation always takes place when the most probable atomic displacement reaches a critical fraction (~20%) of the average interatomic distance, irrespective of whether the relaxation is induced by temperature (linear response) or by mechanical strain (nonlinear response). Such a unified scenario, analogous to the well-known Lindemann criterion for crystal melting, provides insight into the structural origin of the strain-induced glass-liquid transition.