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

Correlation between crystalline order and vitrification in colloidal monolayers

We investigate experimentally the relationship between local structure and dynamical arrest in a quasi-2d colloidal model system which approximates hard discs. We introduce polydispersity to the system to suppress crystallisation. Upon compression, the increase in structural relaxation time is accompanied by the emergence of local hexagonal symmetry. Examining the dynamical heterogeneity of the system, we identify three types of motion: 'zero-dimensional' corresponding to β-relaxation, 'one-dimensional' or stringlike motion and '2D' motion. The dynamic heterogeneity is correlated with the local order, that is to say locally hexagonal regions are more likely to be dynamically slow. However, we find that lengthscales corresponding to dynamic heterogeneity and local structure do not appear to scale together approaching the glass transition.

Simultaneous and accurate measurement of the dielectric constant at many frequencies spanning a wide range

We present an innovative technique which allows the simultaneous measurement of the dielectric constant of a material at many frequencies, spanning a four orders of magnitude range chosen between 10(-2) Hz and 10(4) Hz. The sensitivity and accuracy are comparable to those obtained using standard single frequency techniques. The technique is based on three new and simple features: (a) the precise real time correction of the amplification of a current amplifier, (b) the specific shape of the excitation signal and its frequency spectrum, and (c) the precise synchronization between the generation of the excitation signal and the acquisition of the dielectric response signal. This technique is useful in the case of relatively fast dynamical measurements when the knowledge of the time evolution of the dielectric constant is needed.

Third order dielectric time dependent susceptibilities and the "box model"

We present here the calculation of the "time dependent" moduli of the third order, dielectric, susceptibilities measured at the frequency ω of the applied electric field and at the frequency 3ω. The calculation is performed within the frame work of the "box model," both in the ideal case of a pure third order polarization and in the practical case studied by Samanta and Richert [J. Chem. Phys. 140, 247101 (2014)] where a first order polarization contribution is also detected. We show that, in the two cases, those two modules have a largely different dynamics and that the results can be easily compared with experiments thanks to the data gathered by Samanta and Richert [J. Chem. Phys. 140, 247101 (2014)]. This should provide a new test on the validity of the "box model."

Cooperativity and heterogeneity in plastic crystals studied by nonlinear dielectric spectroscopy

The glassy dynamics of plastic-crystalline cyclo-octanol and ortho-carborane, where only the molecular reorientational degrees of freedom freeze without long-range order, is investigated by nonlinear dielectric spectroscopy. Marked differences to canonical glass formers show up: While molecular cooperativity governs the glassy freezing, it leads to a much weaker slowing down of molecular dynamics than in supercooled liquids. Moreover, the observed nonlinear effects cannot be explained with the same heterogeneity scenario recently applied to canonical glass formers. This supports ideas that molecular relaxation in plastic crystals may be intrinsically nonexponential. Finally, no nonlinear effects were detected for the secondary processes in cyclo-octanol.

Probing cooperative liquid dynamics with the mean square displacement

Literature data for picosecond mean square displacements show that the anharmonicity explains only about half of the fragility (with different fractions for different glass formers). The other half must be ascribed to the Adam-Gibbs mechanism of a growing cooperatively rearranging region. One can measure both influences separately by a simultaneous measurement of liquid and crystal in the coexistence region.

Frequency dependence of dielectric saturation

Dielectric saturation originates from the upper bound to dipole orientation, reached when all dipoles are aligned "perfectly" with respect to the electric field. For Debye-type dynamics, it is well established that the saturation effect is diminished at high frequencies relative to its steady-state value. Here, it is argued that a similar frequency dependence of this nonlinear dielectric effect is expected also for dispersive dynamics, provided that the system is dynamically homogeneous. By contrast, more realistic relaxation time dispersions based upon heterogeneous dynamics display a strongly reduced frequency dependence of dielectric saturation. Calculations demonstrate this effect in terms of both the fundamental and the third harmonic frequency susceptibilities. The relations of signatures of nonlinearity in different Fourier components of the response are briefly discussed.

Cooperativity and the freezing of molecular motion at the glass transition

The slowing down of molecular dynamics when approaching the glass transition generally proceeds much stronger than expected for thermally activated motions. This strange phenomenon can be formally ascribed to a temperature-dependent activation energy E(T). In the present work, via measurements of the third-order nonlinear dielectric susceptibility, we deduce the increase of the number of correlated molecules N(corr) when approaching the glass transition and find a surprisingly simple correlation of E(T) and N(corr)(T). This provides strong evidence that the noncanonical temperature development of glassy dynamics is caused by a temperature-dependent energy barrier arising from the cooperative motion of ever larger numbers of molecules at low temperatures.

Dynamic length scales in glass-forming liquids: an inhomogeneous molecular dynamics simulation approach

In this work, we numerically investigate a new method for the characterization of growing length scales associated with spatially heterogeneous dynamics of glass-forming liquids. This approach, motivated by the formulation of the inhomogeneous mode-coupling theory (IMCT) [Biroli, G.; et al. Phys. Rev. Lett. 2006 97, 195701], utilizes inhomogeneous molecular dynamics simulations in which the system is perturbed by a spatially modulated external potential. We show that the response of the two-point correlation function to the external field allows one to probe dynamic correlations. We examine the critical properties shown by this function, in particular, the associated dynamic correlation length, that is found to be comparable to the one extracted from standardly employed four-point correlation functions. Our numerical results are in qualitative agreement with IMCT predictions but suggest that one has to take into account fluctuations not included in this mean-field approach to reach quantitative agreement. Advantages of our approach over the more conventional one based on four-point correlation functions are discussed.

Emergence of crystal-like atomic dynamics in glasses at the nanometer scale

The vibrational dynamics of a permanently densified silica glass is compared to the one of an α-quartz polycrystal, the silica polymorph of the same density and local structure. The combined use of inelastic x-ray scattering experiments and ab initio numerical calculations provides compelling evidence of a transition, in the glass, from the isotropic elastic response at long wavelengths to a microscopic regime as the wavelength decreases below a characteristic length ξ of a few nanometers, corresponding to about 20 interatomic distances. In the microscopic regime the glass vibrations closely resemble those of the polycrystal, with excitations related to the acoustic and optic modes of the crystal. A coherent description of the experimental results is obtained assuming that the elastic modulus of the glass presents spatial heterogeneities of an average size a ~ ξ/2 π.

Nonlinear dielectric response at the excess wing of glass-forming liquids

We present nonlinear dielectric measurements of glass-forming glycerol and propylene carbonate applying electrical fields up to 671 kV/cm. The measurements extend to sufficiently high frequencies to allow for the investigation of the nonlinear behavior in the regime of the so-far mysterious excess wing, showing up in the loss spectra of many glass formers as a second power law at high frequencies. Surprisingly, we find a complete lack of nonlinear behavior in the excess wing, in marked contrast to the α relaxation where, in agreement with previous reports, a strong increase of dielectric constant and loss is found.

Evidence of growing spatial correlations during the aging of glassy glycerol

We have measured, as a function of the age t(a), the aging of the nonlinear dielectric susceptibility χ(3) of glycerol below the glass transition. Whereas the linear susceptibility can be accurately accounted for in terms of an age dependent relaxation time τ(α)(t(a)), this scaling breaks down for χ(3), suggesting an increase of the amplitude of χ(3). This is a strong indication that the number N(corr) of molecules involved in relaxation events increases with t(a). For T=0.96×T(g), we find that N(corr) increases by ~10% when t(a) varies from 1 to 100 ks. This sheds new light on the relation between length scales and time scales in glasses.

Using mutual information to measure order in model glass formers

Whether or not there is growing static order accompanying the dynamical heterogeneity and increasing relaxation times seen in glassy systems is a matter of dispute. An obstacle to resolving this issue is that the order is expected to be amorphous and so not amenable to simple order parameters. We use mutual information to provide a general measurement of order that is sensitive to multiparticle correlations. We apply this to two glass-forming systems (two-dimensional binary mixtures of hard disks with different size ratios to give varying amounts of hexatic order) and show that there is little growth of amorphous order in the system without crystalline order. In both cases we measure the dynamical length with a four-point correlation function and find that it increases significantly faster than the static lengths in the system as density is increased. We further show that we can recover the known scaling of the dynamic correlation length in a kinetically constrained model, the two-vacancy-assisted-hopping triangular lattice gas.

Cooperative mercury motion in the ionic conductor Cu2HgI4

We present the observation of glasslike dynamic correlations of mobile mercury ions in the ionic conductor Cu2HgI4, detected in both NMR and nonlinear conductivity experiments. The results show that dynamic cooperativity appears in systems seemingly unrelated to glassy and soft arrested materials. A simple kinetic two-component model is proposed, which seems to provide a good description of the cooperative ionic dynamics.

Nonlinear response theory for Markov processes: simple models for glassy relaxation

The theory of nonlinear response for Markov processes obeying a master equation is formulated in terms of time-dependent perturbation theory for the Green's functions and general expressions for the response functions up to third order in the external field are given. The nonlinear response is calculated for a model of dipole reorientations in an asymmetric double well potential, a standard model in the field of dielectric spectroscopy. The static nonlinear response is finite with the exception of a certain temperature T_{0} determined by the value of the asymmetry. In a narrow temperature range around T_{0}, the modulus of the frequency-dependent cubic response shows a peak at a frequency on the order of the relaxation rate and it vanishes for both low frequencies and high frequencies. At temperatures at which the static response is finite (lower and higher than T_{0}), the modulus is found to decay monotonously from the static limit to zero at high frequencies. In addition, results of calculations for a trap model with a Gaussian density of states are presented. In this case, the cubic response depends on the specific dynamical variable considered and also on the way the external field is coupled to the kinetics of the model. In particular, a set of different dynamical variables that gives rise to identical shapes of the linear susceptibility and only to different temperature dependencies of the relaxation times is considered. It is found that the frequency dependence of the nonlinear response functions, however, strongly depends on the particular choice of the variables. The results are discussed in the context of recent theoretical and experimental findings regarding the nonlinear response of supercooled liquids and glasses.

A simple device for dielectric spectroscopy of polymers with temperature regulation close to 300 K based on a Peltier junction

We present a simple thermostat device for performing dielectric spectroscopy measurements on polymers close to their glass transition temperature. By using a vacuum chamber containing a Peltier junction with its regulator, we show that a very simple setup yields a temperature accuracy which is good enough for accurate studies of polymer dielectric properties. This technique is also more cost effective than standard setups using cryogenic fluids.

Third harmonics nonlinear susceptibility in supercooled liquids: a comparison to the box model

The box model, originally introduced to account for the nonresonant hole burning (NHB) dielectric experiments in supercooled liquids, is compared to the measurements of the third harmonics P(3) of the polarisation, reported recently in glycerol, close to the glass transition temperature T(g) [C. Crauste-Thibierge, C. Brun, F. Ladieu, D. L'Hôte, G. Biroli, and J.-P. Bouchaud, Phys. Rev. Lett. 104, 165703 (2010)]. In this model, each box is a distinct dynamical relaxing entity (hereafter called dynamical heterogeneity (DH)) which follows a Debye dynamics with its own relaxation time τ(dh). When it is submitted to a strong electric field, the model posits that a temperature increase δT(dh), depending on τ(dh), arises due to the dissipation of the electrical power. Each DH has thus its own temperature increase, on top of the temperature increase of the phonon bath δT(ph). Contrary to the "fast" hole burning experiments where δT(ph) is usually neglected, the P(3) measurements are, from a thermal point of view, fully in a stationary regime, which means that δT(ph) can no longer be neglected a priori. This is why the version of the box model that we study here takes δT(ph) into account, which implies that the δT(dh) of the DHs are all coupled together. The value of P(3), including both the "intrinsic" contribution of each DH as well as the "spurious" one coming from δT(ph), is computed within this box model and compared to the P(3) measurements for glycerol, in the same range of frequencies and temperatures T. Qualitatively, we find that this version of the box model shares with experiments some nontrivial features, e.g., the existence of a peak at finite frequency in the modulus of P(3) as well as its order of magnitude. Quantitatively, however, some experimental features are not accounted for by this model. We show that these differences between the model and the experiments do not come from δT(ph) but from the "intrinsic" contribution of the DHs. Finally, we show that the interferences between the 3ω response of the various DHs are the most important issue leading to the discrepancies between the box model prediction and the experiments. We argue that this could explain why the box model is quite successful to account for some kinds of nonlinear experiments (such as NHB) performed close to T(g), even if it does not completely account for all of them (such as the P(3) measurements). This conclusion is supported by an analytical argument which helps understanding how a "space-free" model as the box model is able to account for some of the experimental nonlinear features.

Dynamic heterogeneities, boson peak, and activation volume in glass-forming liquids

There are various arguments and models connecting the characteristic length associated with the boson peak vibrations ξ to the length scale of dynamical heterogeneity L(het). ξ is usually defined as the ratio of the transverse sound velocity to the boson peak frequency. Here we present pressure, temperature, and molecular weight dependencies of ξ, estimated using light scattering, in a few molecular and polymeric glass formers. These dependencies are compared with respective dependencies of the activation volume ΔV(#) in the same materials. Good agreement is found for the pressure and molecular weight dependencies of ξ and ΔV(#) measured at the glass transition temperature T(g). These results provide more evidence for a possible relationship between the sensitivity of structural relaxation to density (activation volume) and the heterogeneity volume. However, contrary to the expectations for L(het), ξ does not decrease with temperature above T(g) in most of the studied materials. The temperature dependence of ξ is compared to that of L(het) in glycerol and orthoterphenyl (OTP) estimated from literature data. The analysis shows a clear difference in the behavior of ξ(T) and ΔV(#)(T) at temperatures above T(g), although ΔV(#)(T)(1/3) and L(het)(T) have similar temperature dependence. Possible reasons for the observed difference are discussed.

Analysis of a growing dynamic length scale in a glass-forming binary hard-sphere mixture

We examine a length scale that characterizes the spatial extent of heterogeneous dynamics in a glass-forming binary hard-sphere mixture up to the mode-coupling volume fraction ϕ(c). First, we characterize the system's dynamics. Then, we utilize a method [Phys. Rev. Lett. 105, 217801 (2010)] to extract and analyze the ensemble-independent dynamic susceptibility χ(4)(t) and the dynamic correlation length ξ(t) for a range of times between the β and α relaxation times. We find that in this time range the dynamic correlation length follows a volume fraction-independent master curve ξ(t)~ln(t). For longer times, ξ(t) departs from this master curve and remains constant up to the largest time at which we can determine the length accurately. In addition to the previously established correlation τ(α)~exp[ξ(τ(α))] between the α relaxation time, τ(α), and the dynamic correlation length at this time, ξ(τ(α)), we also find a similar correlation for the diffusion coefficient D~exp[ξ(τ(α))(θ)] with θ≈0.6. We discuss the relevance of these findings for different theories of the glass transition.

Enthalpy space analysis of the evolution of the primary relaxation time in ultraslowing systems

For decades the Vogel-Fulcher-Tammann equation has dominated the description of dynamics of the non-Arrhenius behavior in glass forming systems. Recently, this dominance has been questioned. Hecksher et al. [Nat. Phys. 4, 737 (2008)], Elmatad et al. [J. Phys. Chem. B 113, 5563 (2009)], and Mauro et al. [Proc. Natl. Acad. Sci. U.S.A. 106, 19780 (2009)] indicated superiority of several equations showing no divergence at a finite (nonzero) temperature. This paper shows distortion-sensitive and derivative based empirical analysis of the validity of leading equations for portraying the previtreous evolution of primary relaxation time.

Glassy nature of water in an ultraconfining disordered material: the case of calcium-silicate-hydrate

We present the structural and dynamic nature of water ultraconfined in the quasi-two-dimensional nanopores of the highly disordered calcium-silicate-hydrate (C-S-H), the major binding phase in cement. Our approach is based on classical molecular simulations. We demonstrate that the C-S-H nanopore space is hydrophilic, particularly because of the nonbridging oxygen atoms on the disordered silicate chains which serve as hydrogen-bond acceptor sites, directionally orienting the hydrogen atoms of the interfacial water molecules toward the calcium-silicate layers. The water in this interlayer space adopts a unique multirange structure: a distorted tetrahedral coordination at short range up to 2.7 Å, a disordered structure similar to that of dense fluids and supercooled phases at intermediate range up to 4.2 Å, and persisting spatial correlations through dipole-dipole interactions up to 10 Å. A three-stage dynamics governs the mean square displacement (MSD) of water molecules, with a clear cage stage characteristic of the dynamics in supercooled liquids and glasses, consistent with its intermediate-range structure identified here. At the intermediate time scales corresponding to the β-relaxation of glassy materials, coincident with the cage stage in MSD, the non-Gaussian parameter indicates a significant heterogeneity in the translational dynamics. This dynamic heterogeneity is induced primarily because of the heterogeneity in the distribution of hydrogen bond strengths. The strongly attractive interactions of water molecules with the calcium silicate walls serve to constrain their motion. Our findings have important implications on describing the cohesion and mechanical behavior of cement from its setting to its aging.

Dynamic heterogeneity in a glass forming fluid: susceptibility, structure factor, and correlation length

We investigate the growth of dynamic heterogeneity in a glassy hard-sphere mixture for volume fractions up to and including the mode-coupling transition. We use an 80,000 particle system to test a new procedure to evaluate a dynamic correlation length ξ(t): we determine the ensemble independent dynamic susceptibility χ(4)(t) and use it to facilitate evaluation of ξ(t) from the small wave vector behavior of the four-point structure factor. We analyze relations between the α relaxation time τ(α), χ(4)(τ(α)), and ξ(τ(α)). We find that mode-coupling-like power laws provide a reasonable description of the data over a restricted range of volume fractions, but the power laws' exponents differ from those predicted by the inhomogeneous mode-coupling theory. We find ξ(τ(α))~ln(τ(α)) over the full range of volume fractions studied, which is consistent with an Adam-Gibbs-type relation.

Connection between slow and fast dynamics of molecular liquids around the glass transition

The mean-square displacement (MSD) was measured by neutron scattering at various temperatures and pressures for a number of molecular glass-forming liquids. The MSD is invariant along the glass-transition line at the pressure studied, thus establishing an "intrinsic" Lindemann criterion for any given liquid. A one-to-one connection between the MSD's temperature dependence and the liquid's fragility is found when the MSD is evaluated on a time scale of ∼4 ns , but does not hold when the MSD is evaluated at shorter times. The findings are discussed in terms of the elastic model and the role of relaxations, and the correlations between slow and fast dynamics are addressed.