Light-scattering spectroscopy of the liquid-glass transition in CaKNO3 and in the molecular glass Salol: Extended-mode-coupling-theory analysis

Glassy dynamics of a binary Voronoi fluid: a mode-coupling analysis

The binary Voronoi mixture is a fluid model whose interactions are derived from the Voronoi-Laguerre tessellation of the configurations of the system. The resulting interactions are local and many-body. Here we perform molecular-dynamics (MD) simulations of an equimolar mixture that is weakly polydisperse and additive. For the first time we study the structural relaxation of this mixture in the supercooled-liquid regime. From the simulations we determine the time- and temperature-dependent coherent and incoherent scattering functions for a large range of wave vectors, as well as the mean-square displacements of both particle species. We perform a detailed analysis of the dynamics by comparing the MD results with the first-principles-based idealized mode-coupling theory (MCT). To this end, we employ two approaches: fits to the asymptotic predictions of the theory, and fit-parameter-free binary MCT calculations based on static-structure-factor input from the simulations. We find that many-body interactions of the Voronoi mixture do not lead to strong qualitative differences relative to similar analyses carried out for simple liquids with pair-wise interactions. For instance, the fits give an exponent parameter λ ≈ 0.746 comparable to typical values found for simple liquids, the wavevector dependence of the Kohlrausch relaxation time is in good qualitative agreement with literature results for polydisperse hard spheres, and the MCT calculations based on static input overestimate the critical temperature, albeit only by a factor of about 1.2. This overestimation appears to be weak relative to other well-studied supercooled-liquid models such as the binary Kob-Andersen Lennard-Jones mixture. Overall, the agreement between MCT and simulation suggests that it is possible to predict several microscopic dynamic properties with qualitative, and in some cases near-quantitative, accuracy based solely on static two-point structural correlations, even though the system itself is inherently governed by many-body interactions.

Generalized mode-coupling theory of the glass transition. I. Numerical results for Percus-Yevick hard spheres

Mode-coupling theory (MCT) constitutes one of the few first-principles-based approaches to describe the physics of the glass transition, but the theory's inherent approximations compromise its accuracy in the activated glassy regime. Here, we show that microscopic generalized mode-coupling theory (GMCT), a recently proposed hierarchical framework to systematically improve upon MCT, provides a promising pathway toward a more accurate first-principles description of glassy dynamics. We present a comprehensive numerical analysis for Percus-Yevick hard spheres by performing explicitly wavenumber- and time-dependent GMCT calculations up to sixth order. Specifically, we calculate the location of the critical point, the associated non-ergodicity parameters, and the time-dependent dynamics of the density correlators at both absolute and reduced packing fractions, and we test several universal scaling relations in the α- and β-relaxation regimes. It is found that higher-order GMCT can successfully remedy some of MCT's pathologies, including an underestimation of the critical glass transition density and an overestimation of the hard-sphere fragility. Furthermore, we numerically demonstrate that the celebrated scaling laws of MCT are preserved in GMCT and that the predicted critical exponents manifestly improve as more levels are incorporated in the GMCT hierarchy. Although formally the GMCT equations should be solved up to infinite order to reach full convergence, our finite-order GMCT calculations unambiguously reveal a uniform convergence pattern for the dynamics. We thus argue that GMCT can provide a feasible and controlled means to bypass MCT's main uncontrolled approximation, offering hope for the future development of a quantitative first-principles theory of the glass transition.

Characteristic length scales of the secondary relaxations in glass-forming glycerol

We investigate the secondary relaxations and their link to the main structural relaxation in glass-forming liquids using glycerol as a model system. We analyze the incoherent neutron scattering signal dependence on the scattering momentum transfer, Q , in order to obtain the characteristic length scale for different secondary relaxations. Such a capability of neutron scattering makes it somewhat unique and highly complementary to the traditional techniques of glass physics, such as light scattering and broadband dielectric spectroscopy, which provide information on the time scale, but not the length scales, of relaxation processes. The choice of suitable neutron scattering techniques depends on the time scale of the relaxation of interest. We use neutron backscattering to identify the characteristic length scale of 0.7 Å for the faster secondary relaxation described in the framework of the mode-coupling theory (MCT). Neutron spin-echo is employed to probe the slower secondary relaxation of the excess wing type at a low temperature ( ∼ 1.13T g . The characteristic length scale for this excess wing dynamics is approximately 4.7 Å. Besides the Q -dependence, the direct coupling of neutron scattering signal to density fluctuation makes this technique indispensable for measuring the length scale of the microscopic relaxation dynamics.

Excess wing in glass-forming glycerol and LiCl-glycerol mixtures detected by neutron scattering

The relaxational dynamics in glass-forming glycerol and glycerol mixed with LiCl is investigated using different neutron scattering techniques. The performed neutron spin echo experiments, which extend up to relatively long relaxation time scales of the order of 10 ns, should allow for the detection of contributions from the so-called excess wing. This phenomenon, whose microscopic origin is controversially discussed, arises in a variety of glass formers and, until now, was almost exclusively investigated by dielectric spectroscopy and light scattering. Here we show that the relaxational process causing the excess wing can also be detected by neutron scattering, which directly couples to density fluctuations.

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.

Relaxational Dynamics and Strength in Supercooled Liquids from Impulsive Stimulated Thermal Scattering

Impulsive stimulated thermal scattering (ISTS), a time-domain light scattering technique, provides a more than 6-decade time range from sub-ns to many ms. It permits characterization of the structural relaxation dynamics and determination of the relaxation strength or Debye-Waller factor in supercooled liquids, and thus allows testing of the mode coupling theory of the liquidglass transition. ISTS experiments were performed on glass formers salol, butylbenzene, and the molten salt [Ca(N03)]0.4[KNO3]0.6. The relaxational dynamics and the Debye-Waller factorfq=0 were obtained. A square-root anomaly was observed in fq=0 (T) at a crossover temperature Tc for all three materials, consistent with the prediction of mode coupling theory.

The role of nucleation in vitrification of supercooled liquids

Low-frequency Raman and differential scanning calorimetry (DSC) investigations were carried out during a structural transformation of supercooled liquid salol (phenyl salicylate) in a wide temperature range. DSC experiments indicate that in the supercooled liquid salol at temperature ~40 K above the glass transition temperature metastable nuclei start to form. During subsequent cooling the nuclei become an important element of the glass structure, and thereby are considered as a measure of the intermediate range order in this glass. It was shown that the crystalline structure of the metastable nuclei differ from that of the stable nuclei. Low-frequency Raman spectra of the glassy salol show a broad band in the spectral range from 14.5 to 17.2 cm(-1); the so called 'Boson peak', which can be interpreted in terms of its relationship to the formation of structured clusters, with typical sizes in the nanometer range (critical radii).

Glassy dynamics and domains: explicit results for the East model

A general matrix-based scheme for analyzing the long-time dynamics in kinetically constrained models such as the East model is presented. The treatment developed here is motivated by the expectation that slowly relaxing spin domains of arbitrary size govern the highly cooperative events that lead to spin relaxation at long times. To account for the role of large spin domains in the dynamics, a complete basis expressed in terms of domains of all sizes is introduced. It is first demonstrated that accounting for single domains of all possible sizes leads to a simple analytical result for the two-time single-spin correlation function in the East model that is in excellent quantitative agreement with simulation data for equilibrium spin-up density values c > or = 0.6. It is then shown that including also two neighboring domains leads to a closed expression that describes the slow relaxation of the system down to c approximately 0.3. Ingredients of generalizing the method to lower values of c are also provided, as well as to other models. The main advantage of this approach is that it gives explicit analytical results and that it requires neither an arbitrary closure for the memory kernel nor the construction of an irreducible memory kernel. It also allows one to calculate quantities that measure heterogeneity in the same framework, as is illustrated on the neighbor-pair correlation function, the average relaxation time, and the width of the distribution of relaxation times.

Slow dynamics of salol: a pressure- and temperature-dependent light scattering study

We study the slow dynamics of salol by varying both temperature and pressure using photon correlation spectroscopy and pressure-volume-temperature measurements, and compare the behavior of the structural relaxation time with equations derived within the Adam-Gibbs entropy theory and the Cohen-Grest free volume theory. We find that pressure-dependent data are crucial to assess the validity of these model equations. Our analysis supports the entropy-based equation, and estimates the configurational entropy of salol at ambient pressure approximately 70% of the excess entropy. Finally, we investigate the evolution of the shape of the structural relaxation process, and find that a time-temperature-pressure superposition principle holds over the range investigated.

Brillouin scattering study of salol: exploring the effects of rotation-translation coupling

Brillouin scattering in liquids composed of optically and mechanically anisotropic molecules is affected by coupling between rotational and translational dynamics. While this effect has been extensively studied in depolarized (VH) scattering where it produces the "Rytov dip," recent theoretical analyses by Pick, Franosch show that it should also produce observable effects in polarized (VV) scattering [Eur. Phys. J. B 31, 217 (2003)]; 31, 229 (2003)]]. To test this theory, we carried out Brillouin scattering studies of the molecular glassformer salol in the temperature range 210-380 K, including VH-backscattering, VH-90 degrees, and VV-90 degrees spectra. The data were analyzed consistently to determine the effects of rotation-translation coupling on both the polarized and depolarized spectra. A previously unanticipated feature predicted by these authors was observed: a narrow negative region in the q -dependent part of the 90 degrees VV spectra, which we designate as the "VV dip." It is an analog of the Rytov dip observed at high temperatures in the 90 degrees VH spectra, which is also accurately described by this theory. Analysis of the 90 degrees VV spectra was carried out both with and without inclusion of translation-rotation coupling in order to determine quantitatively the role this coupling plays.

Universality of the dynamic crossover in glass-forming liquids: a "magic" relaxation time

Analysis of experimental data on the structural relaxation time tau(alpha) in various glass formers revealed its universality at the critical temperature T(c) of the mode-coupling theory. In most glass formers studied ln tau(alpha)(T(c))=-(6.5-7.5). Possible reasons for such a universality are discussed.

Anomaly of the nonergodicity parameter and crossover to white noise in the fast relaxation spectrum of a simple glass former

We present quasielastic light scattering and dielectric spectra of the glass former alpha-picoline. At high temperatures the evolution of the susceptibility minimum is well described by the mode coupling theory (MCT). Below the critical temperature T(c) the simple scaling laws of MCT fail due to the appearance of the excess wing of the alpha process, which shows a universal evolution as a function of log(10)tau(alpha). Taking this into account, however, we observe the predicted cusplike anomaly of the nonergodicity parameter as well as a crossover to "white noise."

Microscopic dynamics of molecular liquids and glasses: role of orientations and translation-rotation coupling

We investigate the dynamics of a fluid of dipolar hard spheres in its liquid and glassy phases, with emphasis on the microscopic time or frequency regime. This system shows rather different glass transition scenarios related to its rich equilibrium behavior, which ranges from a simple hard sphere fluid to long range ferroelectric orientational order. In the liquid phase close to the ideal glass transition line and in the glassy regime a medium range orientational order occurs leading to a softening of an orientational mode. To investigate the role of this mode we use the molecular mode-coupling equations to calculate the spectra straight phi"lm(q,omega) and chi"lm(q,omega). In the center of mass spectra straight phi"00(q,omega) and chi"00(q,omega) we found, besides a high frequency peak at omega(hf), a peak at omega(op), about one decade below omega(hf) x omega(op) has almost no q dependence and exhibits an "isotope" effect omega(op) proportional to I(-1/2), with I the moment of inertia. We give evidence that the existence of this peak is related to the occurrence of medium range orientational order. It is shown that some of these features also exist for schematic mode coupling models.

Calcium rubidium nitrate: mode-coupling beta scaling without factorization

The fast dynamics of viscous calcium rubidium nitrate is investigated by depolarized light scattering, neutron scattering, and dielectric loss. Fast beta relaxation evolves as in calcium potassium nitrate. The dynamic susceptibilities can be described by the asymptotic scaling law of mode-coupling theory with a shape parameter lambda=0.79; the temperature dependence of the amplitudes extrapolates to T(c) approximately equal 378 K. However, the frequencies of the minima of the three different spectroscopies never coincide, in conflict with the factorization prediction, indicating that the true asymptotic regime is unreachable.

Molecular mode-coupling theory applied to a liquid of diatomic molecules

We study the molecular mode-coupling theory for a liquid of diatomic molecules. The equations for the critical tensorial nonergodicity parameters F(m)(ll('))(q) and the critical amplitudes of the beta relaxation H(m)(ll('))(q) are solved up to a cutoff l(co)=2 without any further approximations. Here l,m are indices of spherical harmonics. Contrary to previous studies, where additional approximations were applied, we find in agreement with simulations that all molecular degrees of freedom vitrify at a single temperature T(c). The theoretical results for the nonergodicity parameters and the critical amplitudes are compared with those from simulations. The qualitative agreement is good for all molecular degrees of freedom. To study the influence of the cutoff on the nonergodicity parameter, we also calculate the nonergodicity parameters for an upper cutoff l(co)=4. In addition, we also propose a method for the calculation of the critical nonergodicity parameter from the liquid side of transition.

Molecular correlations in a supercooled liquid

We present static and dynamic properties of molecular correlation functions S(lmn,l(')m(')n('))(q-->,t) in a simulated supercooled liquid of water molecules, as a preliminary effort in the direction of solving the molecular mode-coupling theory (MMCT) equations for supercooled molecular liquids. The temperature and time dependence of various molecular correlation functions, calculated from 250 ns long molecular dynamics simulations, show the characteristic patterns predicted by MMCT and shed light on the driving mechanism responsible for the slowing down of the molecular dynamics. We also discuss the symmetry properties of the molecular correlation functions that can be predicted on the basis of the C(2v) symmetry of the molecule. Analysis of the molecular dynamics results for the static correlators S(lmn,l(')m(')n('))(q-->) reveals that additional relationships between correlators with different signs of n and n(') exist. We prove that for molecules with C(rv) symmetry this unexpected result becomes exact at least for high temperatures.

Orientational dynamics in supercooled liquids near T(c) and comparison with ideal mode-coupling theory

Orientational dynamics in supercooled salol and ortho-terphenyl were measured near their critical temperatures, T(c), with optical Kerr effect experiments spanning a very broad range of times. Above T(c), the decays are shown to be in excellent agreement with the master curve predicted by ideal mode-coupling theory when higher order terms are included. Between the critical decay and the von Schweidler power laws, the intermediate time range of the data can be modeled by a power law. This intermediate power law, located at 2<t<10 ps to 500 ps (depending on temperature), shows a significant temperature dependence with a power law exponent of approximately -1 below T(c).

Depolarized light scattering spectroscopy of Ca0.4K0.6(NO3)(1.4): a reexamination of the "knee"

The "knee" found in the depolarized light scattering spectra of Ca0.4K0.6(NO3)(1.4) at low temperatures by G. Li, W.M. Du, X.K. Chen, H.Z. Cummins, and N.J. Tao [Phys. Rev. A 45, 3867 (1992)] appears to have been an experimental artifact. The origin of this feature is analyzed, and its implications for the mode coupling theory of the liquid-glass transition are considered.