Temperature Accelerated Dynamics in Glass-Forming Materials

Structural Transitions in Glassy Atactic Polystyrene Using Transition-State Theory

Transition pathways on the energy landscape of atactic polystyrene (aPS) glassy specimens are probed below its glass-transition temperature. Each of these transitions is considered an elementary structural relaxation event, whose corresponding rate constant is calculated by applying multidimensional transition-state theory. Initially, a wide spectrum of first-order saddle points surrounding local minima on the energy landscape is discovered by a stabilized hybrid eigenmode-following method. Then, (minimal-energy) “reaction paths” to the adjacent minima are constructed by a quadratic descent method. The heights of the free energy, the potential energy, and the entropy barriers are estimated for every connected triplet of transition state and minima. The resulting distribution of free energy barriers is asymmetric and extremely broad, extending to very high barrier heights (over 50 k_BT); the corresponding distribution of rate constants extends over 30 orders of magnitude, with well-defined peaks at the time scales corresponding to the subglass relaxations of polystyrene. Analysis of the curvature along the reaction paths reveals a multitude of different rearrangement mechanisms; some of them bearing multiple distinct phases. Finally, connections to theoretical models of the glass phenomenology allows for the prediction, based on first-principles, of the “ideal” glass-transition temperature entering the Vogel–Fulcher–Tammann (VFT) equation describing the super-Arrhenius temperature dependence of glassy dynamics. Our predictions of the time scales of the subglass relaxations and the VFT temperature are in favorable agreement with available experimental literature data for systems of similar molecular weight under the same conditions.

Computer simulations of glasses: the potential energy landscape

We review the current state of research on glasses, discussing the theoretical background and computational models employed to describe them. This article focuses on the use of the potential energy landscape (PEL) paradigm to account for the phenomenology of glassy systems, and the way in which it can be applied in simulations and the interpretation of their results. This article provides a broad overview of the rich phenomenology of glasses, followed by a summary of the theoretical frameworks developed to describe this phenomonology. We discuss the background of the PEL in detail, the onerous task of how to generate computer models of glasses, various methods of analysing numerical simulations, and the literature on the most commonly used model systems. Finally, we tackle the problem of how to distinguish a good glass former from a good crystal former from an analysis of the PEL. In summarising the state of the potential energy landscape picture, we develop the foundations for new theoretical methods that allow the ab initio prediction of the glass-forming ability of new materials by analysis of the PEL.

Temporal disconnectivity of the energy landscape in glassy systems

An alternative graphical representation of the potential energy landscape (PEL) has been developed and applied to a binary Lennard-Jones glassy system, providing insight into the unique topology of the system's potential energy hypersurface. With the help of this representation one is able to monitor the different explored basins of the PEL, as well as how--and mainly when--subsets of basins communicate with each other via transitions in such a way that details of the prior temporal history have been erased, i.e., local equilibration between the basins in each subset has been achieved. In this way, apart from detailed information about the structure of the PEL, the system's temporal evolution on the PEL is described. In order to gather all necessary information about the identities of two or more basins that are connected with each other, we consider two different approaches. The first one is based on consideration of the time needed for two basins to mutually equilibrate their populations according to the transition rate between them, in the absence of any effect induced by the rest of the landscape. The second approach is based on an analytical solution of the master equation that explicitly takes into account the entire explored landscape. It is shown that both approaches lead to the same result concerning the topology of the PEL and dynamical evolution on it. Moreover, a "temporal disconnectivity graph" is introduced to represent a lumped system stemming from the initial one. The lumped system is obtained via a specially designed algorithm [N. Lempesis, D. G. Tsalikis, G. C. Boulougouris, and D. N. Theodorou, J. Chem. Phys. 135, 204507 (2011)]. The temporal disconnectivity graph provides useful information about both the lumped and the initial systems, including the definition of "metabasins" as collections of basins that communicate with each other via transitions that are fast relative to the observation time. Finally, the two examined approaches are compared to an "on the fly" molecular dynamics-based algorithm [D. G. Tsalikis, N. Lempesis, G. C. Boulougouris, and D. N. Theodorou, J. Chem. Theory Comput. 6, 1307 (2010)].