What do we learn from the shape of the dynamical susceptibility of glass-formers?

TL;DR

This paper analyzes the four-point correlation function in glass-formers, revealing universal features and how different relaxation mechanisms influence the dynamical susceptibility's shape and exponents.

Contribution

It provides an analytical and numerical study of the four-point susceptibility in various glass models, identifying universal behaviors and linking exponents to relaxation mechanisms.

Findings

The four-point susceptibility chi_4(t) exhibits a maximum at t=t^* related to relaxation time.

Exponents mu and lambda distinguish different relaxation mechanisms.

Numerical simulations agree with Mode-Coupling Theory predictions.

Abstract

We compute analytically and numerically the four-point correlation function that characterizes non-trivial cooperative dynamics in glassy systems within several models of glasses: elasto-plastic deformations, mode-coupling theory (MCT), collectively rearranging regions (CRR), diffusing defects and kinetically constrained models (KCM). Some features of the four-point susceptibility chi_4(t) are expected to be universal. at short times we expect an elastic regime characterized by a t or sqrt{t} growth. We find both in the beta, and the early alpha regime that chi_4 sim t^mu, where mu is directly related to the mechanism responsible for relaxation. This regime ends when a maximum of chi_4 is reached at a time t=t^* of the order of the relaxation time of the system. This maximum is followed by a fast decay to zero at large times. The height of the maximum also follows a power-law, chi_4(t^*) sim t^{*lambda}. The value of the exponents mu and lambda allows one to distinguish between different mechanisms. For example, freely diffusing defects in d=3 lead to mu=2 and lambda=1, whereas the CRR scenario rather predicts either mu=1 or a logarithmic behaviour depending on the nature of the nucleation events, and a logarithmic behaviour of chi_4(t^*). MCT leads to mu=b and lambda =1/gamma, where b and gamma are the standard MCT exponents. We compare our theoretical results with numerical simulations on a Lennard-Jones and a soft-sphere system. Within the limited time-scales accessible to numerical simulations, we find that the exponent mu is rather small, mu < 1, with a value in reasonable agreement with the MCT predictions.