1
|
Moch K, Gainaru C, Böhmer R. Nonlinear susceptibilities and higher-order responses related to physical aging: Wiener-Volterra approach and extended Tool-Narayanaswamy-Moynihan models. J Chem Phys 2024; 161:014502. [PMID: 38949281 DOI: 10.1063/5.0207122] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/06/2024] [Accepted: 06/11/2024] [Indexed: 07/02/2024] Open
Abstract
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.
Collapse
Affiliation(s)
- Kevin Moch
- Fakultät Physik, Technische Universität Dortmund, D-44221 Dortmund, Germany
| | - Catalin Gainaru
- Chemical Sciences Division, Oak Ridge National Laboratory, Oak Ridge, Tennessee 37831, USA
| | - Roland Böhmer
- Fakultät Physik, Technische Universität Dortmund, D-44221 Dortmund, Germany
| |
Collapse
|
2
|
Bruns M, Varnik F. Rejuvenation in Deep Thermal Cycling of a Generic Model Glass: A Study of Per-Particle Energy Distribution. MATERIALS 2022; 15:ma15030829. [PMID: 35160779 PMCID: PMC8836976 DOI: 10.3390/ma15030829] [Citation(s) in RCA: 4] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Subscribe] [Scholar Register] [Received: 12/23/2021] [Revised: 01/20/2022] [Accepted: 01/20/2022] [Indexed: 02/04/2023]
Abstract
We investigate the effect of low temperature (cryogenic) thermal cycling on a generic model glass and observe signature of rejuvenation in terms of per-particle potential energy distributions. Most importantly, these distributions become broader and its average values successively increase when applying consecutive thermal cycles. We show that linear dimension plays a key role for these effects to become visible, since we do only observe a weak effect for a cubic system of roughly one hundred particle diameter but observe strong changes for a rule-type geometry with the longest length being two thousand particle diameters. A consistent interpretation of this new finding is provided in terms of a competition between relaxation processes, which are inherent to glassy systems, and excitation due to thermal treatment. In line with our previous report (Bruns et al., PRR 3, 013234 (2021)), it is shown that, depending on the parameters of thermal cycling, rejuvenation can be either too weak to be detected or strong enough for a clear observation.
Collapse
|
3
|
Liu C, Fan Y. Emergent Fractal Energy Landscape as the Origin of Stress-Accelerated Dynamics in Amorphous Solids. PHYSICAL REVIEW LETTERS 2021; 127:215502. [PMID: 34860096 DOI: 10.1103/physrevlett.127.215502] [Citation(s) in RCA: 4] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Received: 04/09/2021] [Revised: 08/23/2021] [Accepted: 10/05/2021] [Indexed: 06/13/2023]
Abstract
The ageing dynamics in a multiplicity of metastable glasses are investigated at various thermomechanical conditions. By using data analytics to deconvolute the integral effects of environmental factors (e.g., energy level, temperature, stress), and by directly scrutinizing the minimum energy pathways for local excitations, we demonstrate external shear would make the system's energy landscape surprisingly fractal and create an emergent low-barrier mode with highly tortuous pathways, leading to an accelerated relaxation. This finding marks a departure from the classic picture of shear-induced simple bias of energy landscape. The insights and implications of this study are also discussed.
Collapse
Affiliation(s)
- Chaoyi Liu
- Department of Mechanical Engineering, University of Michigan, Ann Arbor, Michigan 48109, USA
| | - Yue Fan
- Department of Mechanical Engineering, University of Michigan, Ann Arbor, Michigan 48109, USA
| |
Collapse
|
4
|
Structural relaxation and crystallization in supercooled water from 170 to 260 K. Proc Natl Acad Sci U S A 2021; 118:2022884118. [PMID: 33790015 DOI: 10.1073/pnas.2022884118] [Citation(s) in RCA: 9] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/18/2022] Open
Abstract
The origin of water's anomalous properties has been debated for decades. Resolution of the problem is hindered by a lack of experimental data in a crucial region of temperatures, T, and pressures where supercooled water rapidly crystallizes-a region often referred to as "no man's land." A recently developed technique where water is heated and cooled at rates greater than 109 K/s now enables experiments in this region. Here, it is used to investigate the structural relaxation and crystallization of deeply supercooled water for 170 K < T < 260 K. Water's relaxation toward a new equilibrium structure depends on its initial structure with hyperquenched glassy water (HQW) typically relaxing more quickly than low-density amorphous solid water (LDA). For HQW and T > 230 K, simple exponential relaxation kinetics is observed. For HQW at lower temperatures, increasingly nonexponential relaxation is observed, which is consistent with the dynamics expected on a rough potential energy landscape. For LDA, approximately exponential relaxation is observed for T > 230 K and T < 200 K, with nonexponential relaxation only at intermediate temperatures. At all temperatures, water's structure can be reproduced by a linear combination of two, local structural motifs, and we show that a simple model accounts for the complex kinetics within this context. The relaxation time, τ rel , is always shorter than the crystallization time, τ xtal For HQW, the ratio, τ xtal /τ rel , goes through a minimum at ∼198 K where the ratio is about 60.
Collapse
|
5
|
Mehri S, Ingebrigtsen TS, Dyre JC. Single-parameter aging in a binary Lennard-Jones system. J Chem Phys 2021; 154:094504. [PMID: 33685153 DOI: 10.1063/5.0039250] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
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.
Collapse
Affiliation(s)
- Saeed Mehri
- Glass and Time, IMFUFA, Department of Science and Environment, Roskilde University, P.O. Box 260, DK-4000 Roskilde, Denmark
| | - Trond S Ingebrigtsen
- Glass and Time, IMFUFA, Department of Science and Environment, Roskilde University, P.O. Box 260, DK-4000 Roskilde, Denmark
| | - Jeppe C Dyre
- Glass and Time, IMFUFA, Department of Science and Environment, Roskilde University, P.O. Box 260, DK-4000 Roskilde, Denmark
| |
Collapse
|
6
|
Janssen LMC. Active glasses. JOURNAL OF PHYSICS. CONDENSED MATTER : AN INSTITUTE OF PHYSICS JOURNAL 2019; 31:503002. [PMID: 31469099 DOI: 10.1088/1361-648x/ab3e90] [Citation(s) in RCA: 15] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 06/10/2023]
Abstract
Active glassy matter has recently emerged as a novel class of non-equilibrium soft matter, combining energy-driven, active particle movement with dense and disordered glass-like behavior. Here we review the state-of-the-art in this field from an experimental, numerical, and theoretical perspective. We consider both non-living and living active glassy systems, and discuss how several hallmarks of glassy dynamics (dynamical slowdown, fragility, dynamical heterogeneity, violation of the Stokes-Einstein relation, and aging) are manifested in such materials. We start by reviewing the recent experimental evidence in this area of research, followed by an overview of the main numerical simulation studies and physical theories of active glassy matter. We conclude by outlining several open questions and possible directions for future work.
Collapse
Affiliation(s)
- Liesbeth M C Janssen
- Theory of Polymers and Soft Matter, Department of Applied Physics, Eindhoven University of Technology, PO Box 513, 5600MB Eindhoven, The Netherlands
| |
Collapse
|
7
|
Diezemann G. Nonlinear response theory for Markov processes. II. Fifth-order response functions. Phys Rev E 2017; 96:022150. [PMID: 28950644 DOI: 10.1103/physreve.96.022150] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/03/2017] [Indexed: 06/07/2023]
Abstract
The nonlinear response of stochastic models obeying a master equation is calculated up to fifth order in the external field, thus extending the third-order results obtained earlier [G. Diezemann, Phys. Rev. E 85, 051502 (2012)PLEEE81539-375510.1103/PhysRevE.85.051502]. For sinusoidal fields the 5ω component of the susceptibility is computed for the model of dipole reorientations in an asymmetric double well potential and for a trap model with a Gaussian density of states. For most realizations of the models a hump is found in the higher-order susceptibilities. In particular, for the asymmetric double well potential model there are two characteristic temperature regimes showing the occurrence of such a hump as compared to a single characteristic regime in the case of the third-order response. In the case of the trap model the results strongly depend on the variable coupled to the field. As for the third-order response, the low-frequency limit of the susceptibility plays a crucial role with respect to the occurrence of a hump. The findings are discussed in light of recent experimental results obtained for supercooled liquids. The differences found for the third-order and the fifth-order response indicate that nonlinear response functions might serve as a powerful tool to discriminate among the large number of existing models for glassy relaxation.
Collapse
Affiliation(s)
- Gregor Diezemann
- Institut für Physikalische Chemie, Universität Mainz, Duesbergweg 10-14, 55128 Mainz, Germany
| |
Collapse
|
8
|
Energy landscape-driven non-equilibrium evolution of inherent structure in disordered material. Nat Commun 2017; 8:15417. [PMID: 28524879 PMCID: PMC5454540 DOI: 10.1038/ncomms15417] [Citation(s) in RCA: 22] [Impact Index Per Article: 3.1] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/07/2016] [Accepted: 03/28/2017] [Indexed: 11/08/2022] Open
Abstract
Complex states in glasses can be neatly expressed by the potential energy landscape (PEL). However, because PEL is highly multi-dimensional it is difficult to describe how the system moves around in PEL. Here we demonstrate that it is possible to predict the evolution of macroscopic state in a metallic glass, such as ageing and rejuvenation, through a set of simple equations describing excitations in the PEL. The key to this simplification is the realization that the step of activation from the initial state to the saddle point in PEL and the following step of relaxation to the final state are essentially decoupled. The model shows that the interplay between activation and relaxation in PEL is the key driving force that simultaneously explains both the equilibrium of supercooled liquid and the thermal hysteresis observed in experiments. It further predicts anomalous peaks in truncated thermal scanning, validated by independent molecular dynamics simulation.
Collapse
|
9
|
Dasriaux M, Castagnet S, Thilly L, Chocinski-Arnault L, Boyer SAE. Evolution of the amorphous fraction of PEEK during annealing at atmospheric and high pressure above the glass transition temperature. J Appl Polym Sci 2013. [DOI: 10.1002/app.39297] [Citation(s) in RCA: 11] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/11/2022]
Affiliation(s)
- Marion Dasriaux
- Institut P PRIME (UPR 3346 CNRS - ISAE ENSMA - Université de Poitiers), Chasseneuil-Futuroscope; France
| | - Sylvie Castagnet
- Institut P PRIME (UPR 3346 CNRS - ISAE ENSMA - Université de Poitiers), Chasseneuil-Futuroscope; France
| | - Ludovic Thilly
- Institut P PRIME (UPR 3346 CNRS - ISAE ENSMA - Université de Poitiers), Chasseneuil-Futuroscope; France
| | | | - Séverine A. E. Boyer
- Institut P PRIME (UPR 3346 CNRS - ISAE ENSMA - Université de Poitiers), Chasseneuil-Futuroscope; France
| |
Collapse
|
10
|
Diezemann G. Higher-order correlation functions and nonlinear response functions in a Gaussian trap model. J Chem Phys 2013; 138:12A505. [DOI: 10.1063/1.4769254] [Citation(s) in RCA: 18] [Impact Index Per Article: 1.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
|
11
|
Gnan N, Maggi C, Parisi G, Sciortino F. Generalized fluctuation-dissipation relation and effective temperature upon heating a deeply supercooled liquid. PHYSICAL REVIEW LETTERS 2013; 110:035701. [PMID: 23373934 DOI: 10.1103/physrevlett.110.035701] [Citation(s) in RCA: 5] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Received: 10/12/2012] [Indexed: 06/01/2023]
Abstract
We show that a generalized fluctuation-dissipation relation applies upon instantaneously increasing the temperature of a deeply supercooled liquid. This has the same two-step shape of the relation found upon cooling the liquid, but with opposite violation, indicating an effective temperature that is lower than bath temperature. We show that the effective temperature exhibits some sensible time dependence and that it retains its connection with the partitioned phase space visited in aging. We underline the potential relevance of our numerical results for experimental studies of the fluctuation-dissipation relation in glassy systems.
Collapse
Affiliation(s)
- Nicoletta Gnan
- Dipartimento di Fisica, Università di Roma Sapienza, Roma, Italy.
| | | | | | | |
Collapse
|
12
|
Vollmayr-Lee K, Bjorkquist R, Chambers LM. Microscopic picture of aging in SiO2. PHYSICAL REVIEW LETTERS 2013; 110:017801. [PMID: 23383839 DOI: 10.1103/physrevlett.110.017801] [Citation(s) in RCA: 6] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 09/18/2012] [Indexed: 06/01/2023]
Abstract
We investigate the aging dynamics of amorphous SiO(2) via molecular dynamics simulations of a quench from a high temperature T(i) to a lower temperature T(f). We obtain a microscopic picture of aging dynamics by analyzing single particle trajectories, identifying jump events when a particle escapes the cage formed by its neighbors, and determining how these jumps depend on the waiting time t(w), the time elapsed since the temperature quench to T(f). We find that the only t(w)-dependent microscopic quantity is the number of jumping particles per unit time, which decreases with age. Similar to previous studies for fragile glass formers, we show here for the strong glass former SiO(2) that neither the distribution of jump lengths nor the distribution of times spent in the cage are t(w) dependent. We conclude that the microscopic aging dynamics is surprisingly similar for fragile and strong glass formers.
Collapse
Affiliation(s)
- Katharina Vollmayr-Lee
- Department of Physics and Astronomy, Bucknell University, Lewisburg, Pennsylvania 17837, USA.
| | | | | |
Collapse
|
13
|
Rehwald C, Heuer A. How coupled elementary units determine the dynamics of macroscopic glass-forming systems. PHYSICAL REVIEW. E, STATISTICAL, NONLINEAR, AND SOFT MATTER PHYSICS 2012; 86:051504. [PMID: 23214786 DOI: 10.1103/physreve.86.051504] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 09/27/2012] [Indexed: 06/01/2023]
Abstract
We investigate the dynamics of a binary mixture Lennard-Jones system of different system sizes with respect to the importance of the properties of the underlying potential energy landscape (PEL). We show that the dynamics of small systems can be very well described within the continuous time random walk formalism, which is determined solely by PEL parameters. Finite size analysis shows that the diffusivity of large and small systems are very similar. This suggests that the PEL parameters of the small system also determine the local dynamics in large systems. The structural relaxation time, however, displays significant finite size effects. Furthermore, using a nonequilibrium configuration of a large system, we find that causal connections exist between nearby regions of the system. These findings can be described by the coupled landscape model for which a macroscopic system is described by a superposition of elementary systems, each described by its PEL. A minimum coupling is introduced which accounts for the finite size behavior. The coupling strength, as the single adjustable parameter, becomes smaller closer to the glass transition.
Collapse
Affiliation(s)
- Christian Rehwald
- Institut für Physikalische Chemie, Westfälische Wilhelms-Universität Münster, Corrensstrasse 28/30, 48149 Münster, Germany
| | | |
Collapse
|
14
|
Diezemann G. Nonlinear response theory for Markov processes: simple models for glassy relaxation. PHYSICAL REVIEW. E, STATISTICAL, NONLINEAR, AND SOFT MATTER PHYSICS 2012; 85:051502. [PMID: 23004762 DOI: 10.1103/physreve.85.051502] [Citation(s) in RCA: 16] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Received: 03/08/2012] [Indexed: 06/01/2023]
Abstract
The theory of nonlinear response for Markov processes obeying a master equation is formulated in terms of time-dependent perturbation theory for the Green's functions and general expressions for the response functions up to third order in the external field are given. The nonlinear response is calculated for a model of dipole reorientations in an asymmetric double well potential, a standard model in the field of dielectric spectroscopy. The static nonlinear response is finite with the exception of a certain temperature T_{0} determined by the value of the asymmetry. In a narrow temperature range around T_{0}, the modulus of the frequency-dependent cubic response shows a peak at a frequency on the order of the relaxation rate and it vanishes for both low frequencies and high frequencies. At temperatures at which the static response is finite (lower and higher than T_{0}), the modulus is found to decay monotonously from the static limit to zero at high frequencies. In addition, results of calculations for a trap model with a Gaussian density of states are presented. In this case, the cubic response depends on the specific dynamical variable considered and also on the way the external field is coupled to the kinetics of the model. In particular, a set of different dynamical variables that gives rise to identical shapes of the linear susceptibility and only to different temperature dependencies of the relaxation times is considered. It is found that the frequency dependence of the nonlinear response functions, however, strongly depends on the particular choice of the variables. The results are discussed in the context of recent theoretical and experimental findings regarding the nonlinear response of supercooled liquids and glasses.
Collapse
Affiliation(s)
- Gregor Diezemann
- Institut für Physikalische Chemie, Universität Mainz, Jakob-Welder-Weg 11, 55128 Mainz, Germany
| |
Collapse
|
15
|
Diezemann G, Heuer A. Memory effects in the relaxation of the Gaussian trap model. PHYSICAL REVIEW. E, STATISTICAL, NONLINEAR, AND SOFT MATTER PHYSICS 2011; 83:031505. [PMID: 21517505 DOI: 10.1103/physreve.83.031505] [Citation(s) in RCA: 11] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Received: 11/03/2010] [Indexed: 05/30/2023]
Abstract
We investigate the memory effect in a simple model for glassy relaxation, a trap model with a Gaussian density of states. In this model, thermal equilibrium is reached at all finite temperatures and we therefore can consider jumps from low to high temperatures in addition to the quenches usually considered in aging studies. We show that the evolution of the energy following the Kovacs protocol can approximately be expressed as a difference of two monotonously decaying functions and thus show the existence of a so-called Kovacs hump whenever these functions are not single exponentials. It is well established that the Kovacs effect also occurs in the linear response regime, and we show that most of the gross features do not change dramatically when large temperature jumps are considered. However, there is one distinguishing feature that only exists beyond the linear regime, which we discuss in detail. For the memory experiment with inverted temperatures, i.e., jumping up and then down again, we find a very similar behavior apart from an opposite sign of the hump.
Collapse
Affiliation(s)
- Gregor Diezemann
- Institut für Physikalische Chemie, Universität Mainz, Welderweg 11, DE-55099 Mainz, Federal Republic of Germany
| | | |
Collapse
|