1
|
Shiraishi K, Berthier L. Characterizing the Slow Dynamics of the Swap Monte Carlo Algorithm. J Phys Chem B 2024; 128:12279-12291. [PMID: 39616495 DOI: 10.1021/acs.jpcb.4c06702] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/13/2024]
Abstract
The swap Monte Carlo algorithm introduces nonphysical dynamic rules to accelerate the exploration of the configuration space of supercooled liquids. Its success raises deep questions regarding the nature and physical origin of the slow dynamics of dense liquids and how it is affected by swap moves. We provide a detailed analysis of the slow dynamics generated by the swap Monte Carlo algorithm at very low temperatures in two glass-forming models. We find that the slowing down of the swap dynamics is qualitatively distinct from its local Monte Carlo counterpart, with considerably suppressed dynamic heterogeneity at both single-particle and collective levels. Our results suggest that local kinetic constraints are drastically reduced by swap moves, leading to nearly Gaussian and diffusive dynamics and weakly growing dynamic correlation length scales. The comparison between static and dynamic fluctuations shows that swap Monte Carlo is a nearly optimal local equilibrium algorithm, suggesting that further progress should necessarily involve collective or driven algorithms.
Collapse
Affiliation(s)
- Kumpei Shiraishi
- Laboratoire Charles Coulomb (L2C), Université de Montpellier, CNRS, 34095 Montpellier, France
| | - Ludovic Berthier
- Laboratoire Charles Coulomb (L2C), Université de Montpellier, CNRS, 34095 Montpellier, France
- Gulliver, UMR CNRS 7083, ESPCI Paris, PSL Research University, 75005 Paris, France
| |
Collapse
|
2
|
Sharma A, Liu C, Ozawa M. Selecting relevant structural features for glassy dynamics by information imbalance. J Chem Phys 2024; 161:184506. [PMID: 39530372 DOI: 10.1063/5.0235084] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/25/2024] [Accepted: 10/27/2024] [Indexed: 11/16/2024] Open
Abstract
We numerically investigate the identification of relevant structural features that contribute to the dynamical heterogeneity in a model glass-forming liquid. By employing the recently proposed information imbalance technique, we select these features from a range of physically motivated descriptors. This selection process is performed in a supervised manner (using both dynamical and structural data) and an unsupervised manner (using only structural data). We then apply the selected features to predict future dynamics using a machine learning technique. One of the advantages of the information imbalance technique is that it does not assume any model a priori, i.e., it is a non-parametric method. Finally, we discuss the potential applications of this approach in identifying the dominant mechanisms governing the glassy slow dynamics.
Collapse
Affiliation(s)
- Anand Sharma
- Indian Institute of Science Education and Research, Dr. Homi Bhabha Road, Pashan, Pune 411008, India
- CNRS, LIPhy, Univ. Grenoble Alpes, 38000 Grenoble, France
| | - Chen Liu
- Innovation and Research Division, Ge-Room, Inc., 93160 Noisy le Grand, France
| | - Misaki Ozawa
- CNRS, LIPhy, Univ. Grenoble Alpes, 38000 Grenoble, France
| |
Collapse
|
3
|
Parley JT, Sollich P. Ductile and brittle yielding of athermal amorphous solids: A mean-field paradigm beyond the random-field Ising model. Phys Rev E 2024; 110:045002. [PMID: 39562871 DOI: 10.1103/physreve.110.045002] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/01/2023] [Accepted: 08/30/2024] [Indexed: 11/21/2024]
Abstract
Amorphous solids can yield in either a ductile or brittle manner under strain: plastic deformation can set in gradually, or abruptly through a macroscopic stress drop. Developing a unified theory describing both ductile and brittle yielding constitutes a fundamental challenge of nonequilibrium statistical physics. Recently, it has been proposed that, in the absence of thermal effects, the nature of the yielding transition is controlled by physics akin to that of the quasistatically driven random field Ising model (RFIM), which has served as the paradigm for understanding the effect of quenched disorder in slowly driven systems with short-ranged interactions. However, this theoretical picture neglects both the dynamics of, and the elasticity-induced long-ranged interactions between, the mesoscopic material constituents. Here, we address these two aspects and provide a unified theory building on the Hébraud-Lequeux elastoplastic description. The first aspect is crucial to understanding the competition between the imposed deformation rate and the finite timescale of plastic rearrangements: We provide a dynamical description of the macroscopic stress drop, as well as predictions for the shifting of the brittle yield strain and the scaling of the peak susceptibility with inverse shear rate. The second is essential to capture properly the behavior in the limit of quasistatic driving, where avalanches of plasticity diverge with system size at any value of the strain. We fully characterise the avalanche behavior, which is radically different to that of the RFIM. In the quasistatic, infinite-size limit, we find that both models have mean-field Landau exponents, obscuring the effect of the interactions. We show, however, that the latter profoundly affect the behavior of finite systems approaching the spinodal-like brittle yield point, where we recover qualitatively the finite-size trends found in particle simulations. The interactions also modify the nature of the random critical point separating ductile and brittle yielding, where we predict critical behavior on top of the marginality present at any value of the strain. We finally discuss how all our predictions can be directly tested against particle simulations and eventually experiments, and make first steps in this direction.
Collapse
|
4
|
Semenov A, Baschnagel J. General Relations between Stress Fluctuations and Viscoelasticity in Amorphous Polymer and Glass-Forming Systems. Polymers (Basel) 2024; 16:2336. [PMID: 39204556 PMCID: PMC11359246 DOI: 10.3390/polym16162336] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/28/2024] [Revised: 08/10/2024] [Accepted: 08/14/2024] [Indexed: 09/04/2024] Open
Abstract
Mechanical stress governs the dynamics of viscoelastic polymer systems and supercooled glass-forming fluids. It was recently established that liquids with long terminal relaxation times are characterized by transiently frozen stress fields, which, moreover, exhibit long-range correlations contributing to the dynamically heterogeneous nature of such systems. Recent studies show that stress correlations and relaxation elastic moduli are intimately related in isotropic viscoelastic systems. However, the origin of these relations (involving spatially resolved material relaxation functions) is non-trivial: some relations are based on the fluctuation-dissipation theorem (FDT), while others involve approximations. Generalizing our recent results on 2D systems, we here rigorously derive three exact FDT relations (already established in our recent investigations and, partially, in classical studies) between spatio-temporal stress correlations and generalized relaxation moduli, and a couple of new exact relations. We also derive several new approximate relations valid in the hydrodynamic regime, taking into account the effects of thermal conductivity and composition fluctuations for arbitrary space dimension. One approximate relation was heuristically obtained in our previous studies and verified using our extended simulation data on two-dimensional (2D) glass-forming systems. As a result, we provide the means to obtain, in any spatial dimension, all stress-correlation functions in terms of relaxation moduli and vice versa. The new approximate relations are tested using simulation data on 2D systems of polydisperse Lennard-Jones particles.
Collapse
Affiliation(s)
- Alexander Semenov
- Institut Charles Sadron, CNRS–UPR 22, University of Strasbourg, 67034 Strasbourg, France
| | | |
Collapse
|
5
|
Gavazzoni C, Brito C, Wyart M. Testing Theories of the Glass Transition with the Same Liquid but Many Kinetic Rules. PHYSICAL REVIEW LETTERS 2024; 132:248201. [PMID: 38949336 DOI: 10.1103/physrevlett.132.248201] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Received: 08/02/2023] [Accepted: 04/23/2024] [Indexed: 07/02/2024]
Abstract
We study the glass transition by exploring a broad class of kinetic rules that can significantly modify the normal dynamics of supercooled liquids while maintaining thermal equilibrium. Beyond the usual dynamics of liquids, this class includes dynamics in which a fraction (1-f_{R}) of the particles can perform pairwise exchange or "swap moves," while a fraction f_{P} of the particles can move only along restricted directions. We find that (i) the location of the glass transition varies greatly but smoothly as f_{P} and f_{R} change and (ii) it is governed by a linear combination of f_{P} and f_{R}. (iii) Dynamical heterogeneities (DHs) are not governed by the static structure of the material; their magnitude correlates instead with the relaxation time. (iv) We show that a recent theory for temporal growth of DHs based on thermal avalanches holds quantitatively throughout the (f_{R},f_{P}) diagram. These observations are negative items for some existing theories of the glass transition, particularly those reliant on growing thermodynamic order or locally favored structure, and open new avenues to test other approaches, as we illustrate.
Collapse
Affiliation(s)
- Cristina Gavazzoni
- Instituto de Física, Universidade Federal do Rio Grande do Sul, Caixa Postal 15051, CEP 91501-970, Porto Alegre, Rio Grande do Sul, Brazil
| | - Carolina Brito
- Instituto de Física, Universidade Federal do Rio Grande do Sul, Caixa Postal 15051, CEP 91501-970, Porto Alegre, Rio Grande do Sul, Brazil
| | - Matthieu Wyart
- Institute of Physics, Ecole Polytechnique Federale de Lausanne, 729 BSP UNIL, 1015, Lausanne, Switzerland
| |
Collapse
|
6
|
Pica Ciamarra M, Ji W, Wyart M. Local vs. cooperative: Unraveling glass transition mechanisms with SEER. Proc Natl Acad Sci U S A 2024; 121:e2400611121. [PMID: 38787876 PMCID: PMC11145278 DOI: 10.1073/pnas.2400611121] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/11/2024] [Accepted: 04/17/2024] [Indexed: 05/26/2024] Open
Abstract
Which phenomenon slows down the dynamics in supercooled liquids and turns them into glasses is a long-standing question of condensed matter. Most popular theories posit that as the temperature decreases, many events must occur in a coordinated fashion on a growing length scale for relaxation to occur. Instead, other approaches consider that local barriers associated with the elementary rearrangement of a few particles or "excitations" govern the dynamics. To resolve this conundrum, our central result is to introduce an algorithm, Systematic Excitation ExtRaction, which can systematically extract hundreds of excitations and their energy from any given configuration. We also provide a measurement of the activation energy, characterizing the liquid dynamics, based on fast quenching and reheating. We use these two methods in a popular liquid model of polydisperse particles. Such polydisperse models are known to capture the hallmarks of the glass transition and can be equilibrated efficiently up to millisecond time scales. The analysis reveals that cooperative effects do not control the fragility of such liquids: the change of energy of local barriers determines the change of activation energy. More generally, these methods can now be used to measure the degree of cooperativity of any liquid model.
Collapse
Affiliation(s)
- Massimo Pica Ciamarra
- Division of Physics and Applied Physics, School of Physical and Mathematical Sciences, Nanyang Technological University, 637371, Singapore
- Consiglio Nazionale delle Ricerce, CNR-SPIN, NapoliI-80126, Italy
| | - Wencheng Ji
- Department of Physics of Complex Systems, Weizmann Institute of Science, Rehovot76100, Israel
| | - Matthieu Wyart
- Institute of Physics, École Polytechnique Fédérale de Lausanne, LausanneCH-1015, Switzerland
| |
Collapse
|
7
|
Wang H, Hu L, Xie W, Chang J, Zheng C, Li M, Wang Q, Liao H, Liu D, Wei B. Metastable Liquid Properties and Surface Flow Patterns of Ultrahigh Temperature Alloys Explored in Outer Space. Angew Chem Int Ed Engl 2024; 63:e202400312. [PMID: 38306324 DOI: 10.1002/anie.202400312] [Citation(s) in RCA: 1] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/05/2024] [Revised: 02/01/2024] [Accepted: 02/02/2024] [Indexed: 02/04/2024]
Abstract
The metastable liquid properties and chemical bonds beyond 2000 K remain a huge challenge for ground-based research on liquid materials chemistry. We show the strong undercooling capability, metastable liquid properties and surface wave patterns of refractory Nb-Si and Zr-V binary alloys explored in space environment. The floating droplet of Nb82.7Si17.3 eutectic alloy superheated up to 2338 K exhibited an extreme undercooling of 437 K, approaching the 0.2TE threshold for homogeneous nucleation of liquid-solid reaction. The microgravity state endowed alloy droplets with nearly perfect sphericity and thus ensured the high accuracy to determine metastable undercooled liquid properties. A special kind of swirling flow was induced for liquid alloy owing to Marangoni convection, which resulted in the spiral microstructures on Zr64V36 alloy surface during liquid-solid phase transition. The coupled impacts of surface nucleation and surface flow brought in a novel olivary shape for these binary alloys. Furthermore, the chemical bonds and atomic structures of high temperature liquids were revealed to understand the liquid properties in outer space circumstances.
Collapse
Affiliation(s)
- Haipeng Wang
- Shaanxi Key Laboratory of Space Materials Science and Technology, School of Physical Science and Technology, Northwestern Polytechnical University, Xi'an, 710072, China
| | - Liang Hu
- Shaanxi Key Laboratory of Space Materials Science and Technology, School of Physical Science and Technology, Northwestern Polytechnical University, Xi'an, 710072, China
| | - Wenjun Xie
- Shaanxi Key Laboratory of Space Materials Science and Technology, School of Physical Science and Technology, Northwestern Polytechnical University, Xi'an, 710072, China
| | - Jian Chang
- Shaanxi Key Laboratory of Space Materials Science and Technology, School of Physical Science and Technology, Northwestern Polytechnical University, Xi'an, 710072, China
| | - Chenhui Zheng
- Shaanxi Key Laboratory of Space Materials Science and Technology, School of Physical Science and Technology, Northwestern Polytechnical University, Xi'an, 710072, China
| | - Mingxing Li
- Shaanxi Key Laboratory of Space Materials Science and Technology, School of Physical Science and Technology, Northwestern Polytechnical University, Xi'an, 710072, China
| | - Qing Wang
- Shaanxi Key Laboratory of Space Materials Science and Technology, School of Physical Science and Technology, Northwestern Polytechnical University, Xi'an, 710072, China
| | - Hui Liao
- Shaanxi Key Laboratory of Space Materials Science and Technology, School of Physical Science and Technology, Northwestern Polytechnical University, Xi'an, 710072, China
| | - Dingnan Liu
- Shaanxi Key Laboratory of Space Materials Science and Technology, School of Physical Science and Technology, Northwestern Polytechnical University, Xi'an, 710072, China
| | - Bingbo Wei
- Shaanxi Key Laboratory of Space Materials Science and Technology, School of Physical Science and Technology, Northwestern Polytechnical University, Xi'an, 710072, China
| |
Collapse
|
8
|
Dyre JC. Solid-that-Flows Picture of Glass-Forming Liquids. J Phys Chem Lett 2024; 15:1603-1617. [PMID: 38306474 PMCID: PMC10875679 DOI: 10.1021/acs.jpclett.3c03308] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/24/2023] [Revised: 01/12/2024] [Accepted: 01/17/2024] [Indexed: 02/04/2024]
Abstract
This perspective article reviews arguments that glass-forming liquids are different from those of standard liquid-state theory, which typically have a viscosity in the mPa·s range and relaxation times on the order of picoseconds. These numbers grow dramatically and become 1012 - 1015 times larger for liquids cooled toward the glass transition. This translates into a qualitative difference, and below the "solidity length" which is roughly one micron at the glass transition, a glass-forming liquid behaves much like a solid. Recent numerical evidence for the solidity of ultraviscous liquids is reviewed, and experimental consequences are discussed in relation to dynamic heterogeneity, frequency-dependent linear-response functions, and the temperature dependence of the average relaxation time.
Collapse
Affiliation(s)
- Jeppe C Dyre
- "Glass and Time", IMFUFA, Dept. of Sciences, Roskilde University, P.O. Box 260, DK-4000 Roskilde, Denmark
| |
Collapse
|
9
|
Petersen CF, Harrowell P. Direct measurement of the structural change associated with amorphous solidification using static scattering of coherent radiation. J Chem Phys 2023; 159:244506. [PMID: 38156637 DOI: 10.1063/5.0177251] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/20/2023] [Accepted: 12/06/2023] [Indexed: 01/03/2024] Open
Abstract
In this paper, we demonstrate that the weak temperature dependence of the structure factor of supercooled liquids, a defining feature of the glass transition, is a consequence of the averaging of the scattering intensity due to angular averaging. We show that the speckle at individual wavevectors, calculated from a simulated glass former, exhibits a Debye-Waller factor with a sufficiently large temperature dependence to represent a structural order parameter capable of distinguishing liquid from glass. We also extract from the speckle intensities a quantity proportional to the variance of the local restraint, i.e., a direct experimental measure of the amplitude of structural heterogeneity.
Collapse
Affiliation(s)
- Charlotte F Petersen
- School of Chemistry, University of Sydney, Sydney, New South Wales 2006, Australia
- School of Chemistry, University of Melbourne, Melbourne, Victoria 3010, Australia
| | - Peter Harrowell
- School of Chemistry, University of Sydney, Sydney, New South Wales 2006, Australia
| |
Collapse
|
10
|
Richard D, Kapteijns G, Lerner E. Detecting low-energy quasilocalized excitations in computer glasses. Phys Rev E 2023; 108:044124. [PMID: 37978582 DOI: 10.1103/physreve.108.044124] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/03/2023] [Accepted: 09/15/2023] [Indexed: 11/19/2023]
Abstract
Soft, quasilocalized excitations (QLEs) are known to generically emerge in a broad class of disordered solids and to govern many facets of the physics of glasses, from wave attenuation to plastic instabilities. In view of this key role of QLEs, shedding light upon several open questions in glass physics depends on the availability of computational tools that allow one to study QLEs' statistical mechanics. The latter is a formidable task since harmonic analyses are typically contaminated by hybridizations of QLEs with phononic excitations at low frequencies, obscuring a clear picture of QLEs' abundance, typical frequencies, and other important micromechanical properties. Here we present an efficient algorithm to detect the field of quasilocalized excitations in structural computer glasses. The algorithm introduced takes a computer-glass sample as input and outputs a library of QLEs embedded in that sample. We demonstrate the power of the algorithm by reporting the spectrum of glassy excitations in two-dimensional computer glasses featuring a huge range of mechanical stability, which is inaccessible using conventional harmonic analyses due to phonon hybridizations. Future applications are discussed.
Collapse
Affiliation(s)
- David Richard
- Univ. Grenoble Alpes, CNRS, LIPhy, 38000 Grenoble, France
| | - Geert Kapteijns
- Institute for Theoretical Physics, University of Amsterdam, Science Park 904, 1098 XH Amsterdam, The Netherlands
| | - Edan Lerner
- Institute for Theoretical Physics, University of Amsterdam, Science Park 904, 1098 XH Amsterdam, The Netherlands
| |
Collapse
|
11
|
Jung G, Biroli G, Berthier L. Predicting Dynamic Heterogeneity in Glass-Forming Liquids by Physics-Inspired Machine Learning. PHYSICAL REVIEW LETTERS 2023; 130:238202. [PMID: 37354408 DOI: 10.1103/physrevlett.130.238202] [Citation(s) in RCA: 15] [Impact Index Per Article: 7.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Received: 11/03/2022] [Revised: 03/07/2023] [Accepted: 05/17/2023] [Indexed: 06/26/2023]
Abstract
We introduce GlassMLP, a machine learning framework using physics-inspired structural input to predict the long-time dynamics in deeply supercooled liquids. We apply this deep neural network to atomistic models in 2D and 3D. Its performance is better than the state of the art while being more parsimonious in terms of training data and fitting parameters. GlassMLP quantitatively predicts four-point dynamic correlations and the geometry of dynamic heterogeneity. Transferability across system sizes allows us to efficiently probe the temperature evolution of spatial dynamic correlations, revealing a profound change with temperature in the geometry of rearranging regions.
Collapse
Affiliation(s)
- Gerhard Jung
- Laboratoire Charles Coulomb (L2C), Université de Montpellier, CNRS, 34095 Montpellier, France
| | - Giulio Biroli
- Laboratoire de Physique de l'Ecole Normale Supérieure, ENS, Université PSL, CNRS, Sorbonne Université, Université de Paris, F-75005 Paris, France
| | - Ludovic Berthier
- Laboratoire Charles Coulomb (L2C), Université de Montpellier, CNRS, 34095 Montpellier, France
- Yusuf Hamied Department of Chemistry, University of Cambridge, Lensfield Road, Cambridge CB2 1EW, United Kingdom
| |
Collapse
|
12
|
Ozawa M, Biroli G. Elasticity, Facilitation, and Dynamic Heterogeneity in Glass-Forming Liquids. PHYSICAL REVIEW LETTERS 2023; 130:138201. [PMID: 37067329 DOI: 10.1103/physrevlett.130.138201] [Citation(s) in RCA: 6] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Received: 09/28/2022] [Accepted: 03/02/2023] [Indexed: 06/19/2023]
Abstract
We study the role of elasticity-induced facilitation on the dynamics of glass-forming liquids by a coarse-grained two-dimensional model in which local relaxation events, taking place by thermal activation, can trigger new relaxations by long-range elastically mediated interactions. By simulations and an analytical theory, we show that the model reproduces the main salient facts associated with dynamic heterogeneity and offers a mechanism to explain the emergence of dynamical correlations at the glass transition. We also discuss how it can be generalized and combined with current theories.
Collapse
Affiliation(s)
- Misaki Ozawa
- Laboratoire de Physique de l'Ecole Normale Supérieure, ENS, Université PSL, CNRS, Sorbonne Université, Université de Paris, F-75005 Paris, France
| | - Giulio Biroli
- Laboratoire de Physique de l'Ecole Normale Supérieure, ENS, Université PSL, CNRS, Sorbonne Université, Université de Paris, F-75005 Paris, France
| |
Collapse
|