Narayanaswamy’s 1971 aging theory and material time

How to Distinguish Nonexponentiality and Nonlinearity in Isothermal Structural Relaxation of Glass-Forming Materials

The nonexponentiality and nonlinearity are two essential features of the structural relaxation in any glass-forming material, which seem to be inextricably bound together by the material time. It is shown that the temperature down-jump and up-jump experiments of the same magnitude ΔT = T0 - T to the same temperature T provide a clue for their separation. The isothermal structural relaxation can be quantified using the stabilization period on the logarithmic time scale log(tm/t0). It is described as the sum of the nonexponentiality term 1.181/ß and the nonlinearity term (σ/2.303)ΔT for the temperature down-jump, and as their difference for the temperature up-jump. The material parameter σ = -(∂lnτ/∂Tf)i quantifies variation of the relaxation time with structural changes at the inflection point of the relaxation curve and is formulated for the most widely used phenomenological models. The asymmetry of approach to equilibrium after the temperature down-jump and up-jump was first described by Kovacs in 1963. A detailed analysis of this asymmetry is provided, and a simple method for the estimation of the parameters characterizing the nonexponentiality (ß) and nonlinearity (σ) is proposed. The applicability of this method is tested using previously reported isothermal experimental data as well as calculated data for aging of polymers and other glass-forming materials. This concept illuminates differences in structural relaxation kinetics in a simple and consistent way that can be useful in the design of novel materials and the evaluation of their physical aging treatment.

Nonlinear susceptibilities and higher-order responses related to physical aging: Wiener-Volterra approach and extended Tool-Narayanaswamy-Moynihan models

Large-amplitude thermal excursions imposed on deeply supercooled liquids modulate the nonlinear time evolution of their structural rearrangements. The consequent aftereffects are treated within a Wiener-Volterra expansion in laboratory time that allows one to calculate the associated physical-aging and thermal response functions. These responses and the corresponding higher-harmonic susceptibilities are illustrated using calculations based on the Tool-Narayanaswamy-Moynihan (TNM) model. The conversion from laboratory to material time is thoroughly discussed. Similarities and differences to field-induced higher-harmonic susceptibilities are illustrated using Lissajous and Cole-Cole plots and discussed in terms of aging nonlinearity parameters. For the Lissajous plots, banana-type shapes emerge, while the Cole-Cole plots display cardioidic and other visually appealing patterns. For application beyond the regime in which conventional single-parameter aging concepts work, the Wiener-Volterra material-time-series is introduced as the central tool. Calculations and analyses within this general framework in conjunction with suitable choices of higher-order memory kernels and employing correspondingly extended TNM models yield at least qualitative agreement with recent large-perturbation physical aging experiments. Implications for differential scanning calorimetry and related methods are discussed. The introduced concepts and analyses provide a solid foundation for a generalized description of nonlinear thermal out-of-equilibrium dynamics of glass forming materials, differing from the nonlinear responses known from rheology and dielectric spectroscopy.

Length Scale of Molecular Motions Governing Glass Equilibration in Hyperquenched and Slow-Cooled Polystyrene

Polymer glasses attain thermodynamic equilibrium owing to structural relaxation at various length scales. Herein, calorimetry experiments were conducted to trace the macroscopic relaxation of slow-cooled (SC) and hyperquenched (HQ) polystyrene (PS) glasses and based on detailed comparisons with molecular dynamics probed by dye reorientation, we discussed the possible molecular process governing the equilibration of PS glasses near the glass transition temperatures (Tg). Both SC and HQ glasses equilibrate owing to the cooperative segment motion above a characteristic temperature (Tc) slightly lower than the Tg. In contrast, below the Tc, the localized backbone motion with an apparent activation energy of 290 ± 20 kJ/mol, involving approximately six repeating units, assists equilibrium recovery of PS glasses on the experimentally accessible time scales. The results possibly indicate the presence of an alternative mechanism other than the α-cooperative process controlling physical aging of materials in their deep glassy states.

Distance-as-time in physical aging

Although it has been known for half a century that the physical aging of glasses in experiments is described well by a linear thermal-history convolution integral over the so-called material time, the microscopic definition and interpretation of the material time remains a mystery. We propose that the material-time increase over a given time interval reflects the distance traveled by the system's particles. Different possible distance measures are discussed, starting from the standard mean-square displacement and its inherent-state version that excludes the vibrational contribution. The viewpoint adopted, which is inspired by and closely related to pioneering works of Cugliandolo and Kurchan from the 1990s, implies a "geometric reversibility" and a "unique-triangle property" characterizing the system's path in configuration space during aging. Both of these properties are inherited from equilibrium, and they are here confirmed by computer simulations of an aging binary Lennard-Jones system. Our simulations moreover show that the slow particles control the material time. This motivates a "dynamic-rigidity-percolation" picture of physical aging. The numerical data show that the material time is dominated by the slowest particles' inherent mean-square displacement, which is conveniently quantified by the inherent harmonic mean-square displacement. This distance measure collapses data for potential-energy aging well in the sense that the normalized relaxation functions following different temperature jumps are almost the same function of the material time. Finally, the standard Tool-Narayanaswamy linear material-time convolution-integral description of physical aging is derived from the assumption that when time is replaced by distance in the above sense, an aging system is described by the same expression as that of linear-response theory.

"Inner clocks" of glass-forming liquids

Providing a physically sound explanation of aging phenomena in non-equilibrium amorphous materialsis a challenging problem in modern statistical thermodynamics. The slow evolution of physical propertiesafter quenches of control parameters is empirically well interpreted via the concept of material time (orinternal clock), based on the Tool-Narayanaswamy-Moynihan (TNM) model. Yet, the fundamental reasonsof its striking success remain unclear. We propose a microscopic rationale behind the material time onthe basis of the linear laws of irreversible thermodynamics and its extension that treats the correspondingkinetic coefficients as state functions of a slowly evolving material state. Our interpretation is based onthe recognition that the same mathematical structure governs both the Tool model and the recently devel-oped non-equilibrium extension of the self-consistent generalized Langevin equation theory (NE-SCGLE),guided by the universal principles of Onsager's theory of irreversible processes. This identification opensthe way for a generalization of the material-time concept to aging systems where several relaxation modeswith very different equilibration processes must be considered, and partially frozen glasses manifest theappearance of partial ergodicity breaking, and hence materials with multiple very distinct inner clocks.

Temperature dependence of aging dynamics in highly non-equilibrium model polymer glasses

A central feature of the non-equilibrium glassy “state” is its tendency to age toward equilibrium, obeying signatures identified by Kovacs over 50 years ago. The origin of these signatures, their fate far from equilibrium and at high temperatures, and the underlying nature of the glassy “state” far from equilibrium remain unsettled. Here, we simulate physical aging of polymeric glasses, driven much farther from equilibrium and at much higher temperatures than possible in experimental melt-quenched glasses. While these glasses exhibit Kovacs’ signatures of glassy aging at sufficiently low temperatures, these signatures disappear above the onset TA of non-Arrhenius equilibrium dynamics, suggesting that TA demarcates an upper bound to genuinely glassy states. Aging times in glasses after temperature up-jumps are found to obey an Arrhenius law interpolating between equilibrium dynamics at TA and at the start of the temperature up-jump, providing a zero-parameter rule predicting their aging behavior and identifying another unrecognized centrality of TA to aging behavior. This differs qualitatively from behavior of our glasses produced by temperature down-jumps, which exhibit a fractional power law decoupling relation with equilibrium dynamics. While the Tool–Narayanaswamy–Moynihan model can predict the qualitative single-temperature behavior of these systems, we find that it fails to predict the disappearance of Kovacs signatures above TA and the temperature dependence of aging after large temperature up-jumps. These findings highlight a need for new theoretical insights into the aging behavior of glasses at ultra-high fictive temperatures and far from equilibrium.

Predicting nonlinear physical aging of glasses from equilibrium relaxation via the material time

The noncrystalline glassy state of matter plays a role in virtually all fields of materials science and offers complementary properties to those of the crystalline counterpart. The caveat of the glassy state is that it is out of equilibrium and therefore exhibits physical aging, i.e., material properties change over time. For half a century, the physical aging of glasses has been known to be described well by the material-time concept, although the existence of a material time has never been directly validated. We do this here by successfully predicting the aging of the molecular glass 4-vinyl-1,3-dioxolan-2-one from its linear relaxation behavior. This establishes the defining property of the material time. Via the fluctuation-dissipation theorem, our results imply that physical aging can be predicted from thermal-equilibrium fluctuation data, which is confirmed by computer simulations of a binary liquid mixture.

Single-parameter aging in a binary Lennard-Jones system

This paper studies physical aging by computer simulations of a 2:1 Kob-Andersen binary Lennard-Jones mixture, a system that is less prone to crystallization than the standard 4:1 composition. Starting from thermal-equilibrium states, the time evolution of the following four quantities is monitored by following up and down jumps in temperature: potential energy, virial, average squared force, and the Laplacian of the potential energy. Despite the fact that significantly larger temperature jumps are studied here than in typical similar experiments, to a good approximation, all four quantities conform to the single-parameter-aging scenario derived and validated for small jumps in experiments [T. Hecksher, N. B. Olsen, and J. C. Dyre, J. Chem. Phys. 142, 241103 (2015)]. As a further confirmation of single-parameter aging with a common material time for the four different quantities monitored, their relaxing parts are found to be almost identical for all temperature jumps.

Fast contribution to the activation energy of a glass-forming liquid

This paper presents physical-aging data for the silicone oil tetramethyl-tetraphenyl trisiloxane. The density and the high-frequency plateau shear modulus [Formula: see text] were monitored following temperature jumps starting from fully equilibrated conditions. Both quantities exhibit a fast change immediately after a temperature jump. Adopting the material-time formalism of Narayanaswamy, we determine from the dielectric loss at 0.178 Hz the time evolution of the aging-rate activation energy. The relative magnitude of the fast change of the activation energy differs from that of the density, but is identical to that of [Formula: see text] In fact, the activation energy is proportional to [Formula: see text] throughout the aging process, with minor deviations at the shortest times. This shows that for the silicone oil in question the dynamics are determined by [Formula: see text] in-as well as out of-equilibrium.

The very long-term physical aging of glassy polymers

The thermodynamic state of polymer glasses aged over 30 years reveals the existence of a metastable state with partial equilibrium recovery.

"Swarm relaxation": Equilibrating a large ensemble of computer simulations⋆

It is common practice in molecular dynamics and Monte Carlo computer simulations to run multiple, separately-initialized simulations in order to improve the sampling of independent microstates. Here we examine the utility of an extreme case of this strategy, in which we run a large ensemble of M independent simulations (a "swarm"), each of which is relaxed to equilibrium. We show that if M is of order [Formula: see text], we can monitor the swarm's relaxation to equilibrium, and confirm its attainment, within [Formula: see text], where [Formula: see text] is the equilibrium relaxation time. As soon as a swarm of this size attains equilibrium, the ensemble of M final microstates from each run is sufficient for the evaluation of most equilibrium properties without further sampling. This approach dramatically reduces the wall-clock time required, compared to a single long simulation, by a factor of several hundred, at the cost of an increase in the total computational effort by a small factor. It is also well suited to modern computing systems having thousands of processors, and is a viable strategy for simulation studies that need to produce high-precision results in a minimum of wall-clock time. We present results obtained by applying this approach to several test cases.

Mapping Isobaric Aging onto the Equilibrium Phase Diagram

The linear volume relaxation and the nonlinear volume aging of a glass-forming liquid are measured, directly compared, and used to extract the out-of-equilibrium relaxation time. This opens a window to investigate how the relaxation time depends on temperature, structure, and volume in parts of phase space that are not accessed by the equilibrium liquid. It is found that the temperature dependence of relaxation time is non-Arrhenius even in the isostructural case-challenging the Adam-Gibbs entropy model. Based on the presented data and the idea that aging happens through quasiequilibrium states, we suggest a mapping of the out-of-equilibrium states during isobaric aging to the equilibrium phase diagram. This mapping implies the existence of isostructural lines in the equilibrium phase diagram. The relaxation time is found to depend on the bath temperature, density, and a just single structural parameter, referred to as an effective temperature.

Model for the alpha and beta shear-mechanical properties of supercooled liquids and its comparison to squalane data

This paper presents data for supercooled squalane's frequency-dependent shear modulus covering frequencies from 10 mHz to 30 kHz and temperatures from 168 K to 190 K; measurements are also reported for the glass phase down to 146 K. The data reveal a strong mechanical beta process. A model is proposed for the shear response of the metastable equilibrium liquid phase of supercooled liquids. The model is an electrical equivalent-circuit characterized by additivity of the dynamic shear compliances of the alpha and beta processes. The nontrivial parts of the alpha and beta processes are each represented by a "Cole-Cole retardation element" defined as a series connection of a capacitor and a constant-phase element, resulting in the Cole-Cole compliance function well-known from dielectrics. The model, which assumes that the high-frequency decay of the alpha shear compliance loss varies with the angular frequency as ω^-1/2, has seven parameters. Assuming time-temperature superposition for the alpha and beta processes separately, the number of parameters varying with temperature is reduced to four. The model provides a better fit to the data than an equally parametrized Havriliak-Negami type model. From the temperature dependence of the best-fit model parameters, the following conclusions are drawn: (1) the alpha relaxation time conforms to the shoving model; (2) the beta relaxation loss-peak frequency is almost temperature independent; (3) the alpha compliance magnitude, which in the model equals the inverse of the instantaneous shear modulus, is only weakly temperature dependent; (4) the beta compliance magnitude decreases by a factor of three upon cooling in the temperature range studied. The final part of the paper briefly presents measurements of the dynamic adiabatic bulk modulus covering frequencies from 10 mHz to 10 kHz in the temperature range from 172 K to 200 K. The data are qualitatively similar to the shear modulus data by having a significant beta process. A single-order-parameter framework is suggested to rationalize these similarities.

Universal scaling in the aging of the strong glass former SiO2

We show that the aging dynamics of a strong glass former displays a strikingly simple scaling behavior, connecting the average dynamics with its fluctuations, namely, the dynamical heterogeneities. We perform molecular dynamics simulations of SiO2 with van Beest-Kramer-van Santen interactions, quenching the system from high to low temperature, and study the evolution of the system as a function of the waiting time tw measured from the instant of the quench. We find that both the aging behavior of the dynamic susceptibility χ4 and the aging behavior of the probability distribution P(fs,r) of the local incoherent intermediate scattering function fs,r can be described by simple scaling forms in terms of the global incoherent intermediate scattering function C. The scaling forms are the same that have been found to describe the aging of several fragile glass formers and that, in the case of P(fs,r), have been also predicted theoretically. A thorough study of the length scales involved highlights the importance of intermediate length scales. We also analyze directly the scaling dependence on particle type and on wavevector q and find that both the average and the fluctuations of the slow aging dynamics are controlled by a unique aging clock, which is not only independent of the wavevector q, but is also the same for O and Si atoms.