Growing length and time scales in glass-forming liquids

The glass transition, whereby liquids transform into amorphous solids at low temperatures, is a subject of intense research despite decades of investigation. Explaining the enormous increase in relaxation times of a liquid upon supercooling is essential for understanding the glass transition. Although many theories, such as the Adam-Gibbs theory, have sought to relate growing relaxation times to length scales associated with spatial correlations in liquid structure or motion of molecules, the role of length scales in glassy dynamics is not well established. Recent studies of spatially correlated rearrangements of molecules leading to structural relaxation, termed ``spatially heterogeneous dynamics,'' provide fresh impetus in this direction. A powerful approach to extract length scales in critical phenomena is finite-size scaling, wherein a system is studied for sizes traversing the length scales of interest. We perform finite-size scaling for a realistic glass-former, using computer simulations, to evaluate the length scale associated with spatially heterogeneous dynamics, which grows as temperature decreases. However, relaxation times that also grow with decreasing temperature do not exhibit standard finite-size scaling with this length. We show that relaxation times are instead determined, for all studied system sizes and temperatures, by configurational entropy, in accordance with the Adam-Gibbs relation, but in disagreement with theoretical expectations based on spin-glass models that configurational entropy is not relevant at temperatures substantially above the critical temperature of mode-coupling theory. Our results provide new insights into the dynamics of glass-forming liquids and pose serious challenges to existing theoretical descriptions.

A chemical approach to understand fragilities of glass-forming liquids

Fragility is viewed as a measure of the loss of rigidity of a glass structure above its glass transition temperature. It is attributed to the weakness of directional bonding and to the presence of a high density of low-energy configurational states. An a priori fragility function of electronegativities and bond distances is proposed which quite remarkably reproduces the entire range of reported fragilities and demonstrates that the fragility of a melt is indeed encrypted in the chemistry of the parent material. It has also been shown that the use of fragility-modified activation barriers in the Arrhenius function account for the whole gamut of viscosity behavior of liquids. It is shown that fragility can be a universal scaling parameter to collapse all viscosity curves on to a master plot.

Characterization of the dynamics of glass-forming liquids from the properties of the potential energy landscape

We develop a framework for understanding the difference between strong and fragile behavior in the dynamics of glass-forming liquids from the properties of the potential energy landscape. Our approach is based on a master equation description of the activated jump dynamics among the local minima of the potential energy (the so-called inherent structures) that characterize the potential energy landscape of the system. We study the dynamics of a small atomic cluster using this description as well as molecular dynamics simulations and demonstrate the usefulness of our approach for this system. Many of the remarkable features of the complex dynamics of glassy systems emerge from the activated dynamics in the potential energy landscape of the atomic cluster. The dynamics of the system exhibits typical characteristics of a strong supercooled liquid when the system is allowed to explore the full configuration space. This behavior arises because the dynamics is dominated by a few lowest-lying minima of the potential energy and the potential energy barriers between these minima. When the system is constrained to explore only a limited region of the potential energy landscape that excludes the basins of attraction of a few lowest-lying minima, the dynamics is found to exhibit the characteristics of a fragile liquid.

Adam-Gibbs Relation for Glass-Forming Liquids in Two, Three, and Four Dimensions

The Adam-Gibbs relation between relaxation times and the configurational entropy has been tested extensively for glass formers using experimental data and computer simulation results. Although the form of the relation contains no dependence on the spatial dimensionality in the original formulation, subsequent derivations of the Adam-Gibbs relation allow for such a possibility. We test the Adam-Gibbs relation in two, three, and four spatial dimensions using computer simulations of model glass formers. We find that the relation is valid in three and four dimensions. But in two dimensions, the relation does not hold, and interestingly, no single alternate relation describes the results for the different model systems we study.

Growing Length Scales and Their Relation to Timescales in Glass-Forming Liquids

The question of whether the dramatic slowing down of the dynamics of glass-forming liquids near the structural glass transition is caused by the growth of one or more correlation lengths has received much attention in recent years. Several proposals have been made for both static and dynamic length scales that may be responsible for the growth of timescales as the glass transition is approached. These proposals are critically examined with emphasis on the dynamic length scale associated with spatial heterogeneity of local dynamics and the static point-to-set or mosaic length scale of the random first order transition theory of equilibrium glass transition. Available results for these length scales, obtained mostly from simulations, are summarized, and the relation of the growth of timescales near the glass transition with the growth of these length scales is examined. Some of the outstanding questions about length scales in glass-forming liquids are discussed, and studies in which these questions may be addressed are suggested.

Fragility correlates thermodynamic and kinetic properties of glass forming liquids

In our earlier communication we proposed a simple fragility determining function, (NBO]/(VmTg)-T-3), which we have now used to analyze several glass systems using available thermal data. A comparison with similar fragility determining function, Delta C-p/C-p(1), introduced by Chryssikos et al. in their investigation of lithium borate glasses has also been performed and found to be more convenient quantity for discussing fragilities. We now propose a new function which uses both Delta C-p and Delta T and which gives a numerical fragility parameter, F whose value lies between 0 and 1 for glass forming liquids. F can be calculated through the use of measured thermal parameters Delta C-p, C-p(1), T-g and T-m. Use of the new fragility values in reduced viscosity equation reproduces the whole range of viscosity curves of the Angell plot. The reduced viscosity equation can be directly compared with the Adam-Gibbs viscosity equation and a heat capacity function can be formulated which reproduces satisfactorily the Delta C-p versus In(T-r) curves and hence the configurational entropy. (C) 2014 Elsevier Ltd. All rights reserved.

Fragility correlates thermodynamic and kinetic properties of glass forming liquids

In our earlier communication we proposed a simple fragility determining function, (NBO]/(VmTg)-T-3), which we have now used to analyze several glass systems using available thermal data. A comparison with similar fragility determining function, Delta C-p/C-p(1), introduced by Chryssikos et al. in their investigation of lithium borate glasses has also been performed and found to be more convenient quantity for discussing fragilities. We now propose a new function which uses both Delta C-p and Delta T and which gives a numerical fragility parameter, F whose value lies between 0 and 1 for glass forming liquids. F can be calculated through the use of measured thermal parameters Delta C-p, C-p(1), T-g and T-m. Use of the new fragility values in reduced viscosity equation reproduces the whole range of viscosity curves of the Angell plot. The reduced viscosity equation can be directly compared with the Adam-Gibbs viscosity equation and a heat capacity function can be formulated which reproduces satisfactorily the Delta C-p versus In(T-r) curves and hence the configurational entropy. (C) 2014 Elsevier Ltd. All rights reserved.

Dynamics of Glass Forming Liquids with Randomly Pinned Particles

It is frequently assumed that in the limit of vanishing cooling rate, the glass transition phenomenon becomes a thermodynamic transition at a temperature T-K. However, with any finite cooling rate, the system falls out of equilibrium at temperatures near T-g(> T-K), implying that the very existence of the putative thermodynamic phase transition at T-K can be questioned. Recent studies of systems with randomly pinned particles have hinted that the thermodynamic glass transition may be observed for liquids with randomly pinned particles. This expectation is based on the results of approximate calculations that suggest that the thermodynamic glass transition temperature increases with increasing concentration of pinned particles and it may be possible to equilibrate the system at temperatures near the increased transition temperature. We test the validity of this prediction through extensive molecular dynamics simulations of two model glass-forming liquids in the presence of random pinning. We find that extrapolated thermodynamic transition temperature T-K does not show any sign of increasing with increasing pinning concentration. The main effect of pinning is found to be a rapid decrease in the kinetic fragility of the system with increasing pin concentration. Implications of these observations for current theories of the glass transition are discussed.

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.

Short-Time Beta Relaxation in Glass-Forming Liquids Is Cooperative in Nature

Temporal relaxation of density fluctuations in supercooled liquids near the glass transition occurs in multiple steps. Using molecular dynamics simulations for three model glass-forming liquids, we show that the short-time beta relaxation is cooperative in nature. Using finite-size scaling analysis, we extract a growing length scale associated with beta relaxation from the observed dependence of the beta relaxation time on the system size. We find, in qualitative agreement with the prediction of the inhomogeneous mode coupling theory, that the temperature dependence of this length scale is the same as that of the length scale that describes the spatial heterogeneity of local dynamics in the long-time alpha-relaxation regime.

Inorganic glasses, glass-forming liquids and amorphizing solids

Greaves, George; Sen, S., (2007) 'Inorganic glasses, glass-forming liquids and amorphizing solids', Advances in Physics 56(1) pp.1-166 RAE2008

The new insight into dynamic crossover in glass forming liquids from the apparent enthalpy analysis

One of the most intriguing phenomena in glass forming systems is the dynamic crossover (T(B)), occurring well above the glass temperature (T(g)). So far, it was estimated mainly from the linearized derivative analysis of the primary relaxation time τ(T) or viscosity η(T) experimental data, originally proposed by Stickel et al. [J. Chem. Phys. 104, 2043 (1996); J. Chem. Phys. 107, 1086 (1997)]. However, this formal procedure is based on the general validity of the Vogel-Fulcher-Tammann equation, which has been strongly questioned recently [T. Hecksher et al. Nature Phys. 4, 737 (2008); P. Lunkenheimer et al. Phys. Rev. E 81, 051504 (2010); J. C. Martinez-Garcia et al. J. Chem. Phys. 134, 024512 (2011)]. We present a qualitatively new way to identify the dynamic crossover based on the apparent enthalpy space (H(a)(') = dlnτ/d(1/T)) analysis via a new plot lnH(a)(') vs. 1∕T supported by the Savitzky-Golay filtering procedure for getting an insight into the noise-distorted high order derivatives. It is shown that depending on the ratio between the "virtual" fragility in the high temperature dynamic domain (m(high)) and the "real" fragility at T(g) (the low temperature dynamic domain, m = m(low)) glass formers can be splitted into two groups related to f < 1 and f > 1, (f = m(high)∕m(low)). The link of this phenomenon to the ratio between the apparent enthalpy and activation energy as well as the behavior of the configurational entropy is indicated.

Localized Excitations and the Morphology of Cooperatively Rearranging Regions in a Colloidal Glass-Forming Liquid

We develop a scheme based on a real space microscopic analysis of particle dynamics to ascertain the relevance of dynamical facilitation as a mechanism of structural relaxation in glass-forming liquids. By analyzing the spatial organization of localized excitations within clusters of mobile particles in a colloidal glass former and examining their partitioning into shell-like and corelike regions, we establish the existence of a crossover from a facilitation-dominated regime at low area fractions to a collective activated hopping-dominated one close to the glass transition. This crossover occurs in the vicinity of the area fraction at which the peak of the mobility transfer function exhibits a maximum and the morphology of cooperatively rearranging regions changes from stringlike to a compact form. Collectively, our findings suggest that dynamical facilitation is dominated by collective hopping close to the glass transition, thereby constituting a crucial step towards identifying the correct theoretical scenario for glass formation.

Glass forming ability and stability: Ternary Cu bearing Ti, Zr, Hf alloys

Four Cu bearing alloys of nominal composition Zr25Ti25Cu50, Zr34Ti16Cu50, Zr25Hf25Cu50 and Ti25Hf25Cu50 have been rapidly solidified in order to produce ribbons. All the alloys become amorphous after meltspinning. In the Zr34Ti16Cu50 alloy localized precipitation of cF24 Cu5Zr phase can be observed in the amorphous matrix. The alloys show a tendency of phase separation at the initial stages of crystallization. The difference in crystallization behavior of these alloys with Ni bearing ternary alloys can be explained by atomic size, binary heat of mixing and Mendeleev number. It has been observed that both Laves and Anti-Laves phase forming compositions are suitable for glass formation. The structures of the phases, precipitated during rapid solidification and crystallization can be viewed in terms of Bernal deltahedra and Frank-Kasper polyhedra.

Analysis of Dynamic Heterogeneity in a Glass Former from the Spatial Correlations of Mobility

The growth of characteristic length scales associated with dynamic heterogeneity in glass-forming liquids is investigated in an extensive computational study of a four-point, time-dependent structure factor defined from spatial correlations of mobility, for a model liquid for system sizes extending up to 351 232 particles, in constant-energy and constant-temperature ensembles. Our estimates for dynamic correlation lengths and susceptibilities are consistent with previous results from finite size scaling. We find scaling exponents that are inconsistent with predictions from inhomogeneous mode coupling theory and a recent simulation confirmation of these predictions.