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.

Elliptical tracers in two-dimensional, homogeneous, isotropic fluid turbulence: The statistics of alignment, rotation, and nematic order

We study the statistical properties of orientation and rotation dynamics of elliptical tracer particles in two-dimensional, homogeneous, and isotropic turbulence by direct numerical simulations. We consider both the cases in which the turbulent flow is generated by forcing at large and intermediate length scales. We show that the two cases are qualitatively different. For large-scale forcing, the spatial distribution of particle orientations forms large-scale structures, which are absent for intermediate-scale forcing. The alignment with the local directions of the flow is much weaker in the latter case than in the former. For intermediate-scale forcing, the statistics of rotation rates depends weakly on the Reynolds number and on the aspect ratio of particles. In contrast with what is observed in three-dimensional turbulence, in two dimensions the mean-square rotation rate increases as the aspect ratio increases.

Relation between the Isentropic Index and the Gruneisen Parameter for Saturated Liquids

We investigate the isentropic index along the saturated vapor line as a correlating parameter with quantities both in the saturated liquid phase and the saturated vapor phase. The relation is established via temperatures such as T-hgmax and T* where the saturated vapor enthalpy and the product of saturation temperature and saturated liquid density attain a maximum, respectively. We obtain that the saturated vapor isentropic index is correlated with these temperatures but also with the saturated liquid Gruneisen parameters at T-hgmax. and T*.

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.

Organic alloys of room temperature liquids thiophenol and selenophenol

The first examples of organic alloys of two room temperature liquids, obtained and characterized via in situ cryo-crystallography, are presented. Thiophenol and selenophenol, which exhibit isostructurality and similar modes of S center dot center dot center dot S and Se center dot center dot center dot Se homo-chalcogen interactions along with weak and rare S-H center dot center dot center dot S and Se-H center dot center dot center dot Se hydrogen bonds, are shown to form solid solutions exhibiting Veggard's law-like trends.

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.