Spatially heterogeneous dynamics investigated via a time-dependent four-point density correlation function

• Published in: The Journal of chemical physics 2003-10, Vol.119 (14), p.7372-7387
• Main Authors: Lačević, N., Starr, F. W., Schrøder, T. B., Glotzer, S. C.
• Format: Article
• Language: English
• ISSN: 0021-9606; 1089-7690
• DOI: 10.1063/1.1605094

Relaxation in supercooled liquids above their glass transition and below the onset temperature of “slow” dynamics involves the correlated motion of neighboring particles. This correlated motion results in the appearance of spatially heterogeneous dynamics or “dynamical heterogeneity.” Traditional two-point time-dependent density correlation functions, while providing information about the transient “caging” of particles on cooling, are unable to provide sufficiently detailed information about correlated motion and dynamical heterogeneity. Here, we study a four-point, time-dependent density correlation function g4(r,t) and corresponding “structure factor” S4(q,t) which measure the spatial correlations between the local liquid density at two points in space, each at two different times, and so are sensitive to dynamical heterogeneity. We study g4(r,t) and S4(q,t) via molecular dynamics simulations of a binary Lennard-Jones mixture approaching the mode coupling temperature from above. We find that the correlations between particles measured by g4(r,t) and S4(q,t) become increasingly pronounced on cooling. The corresponding dynamical correlation length ξ4(t) extracted from the small-q behavior of S4(q,t) provides an estimate of the range of correlated particle motion. We find that ξ4(t) has a maximum as a function of time t, and that the value of the maximum of ξ4(t) increases steadily from less than one particle diameter to a value exceeding nine particle diameters in the temperature range approaching the mode coupling temperature from above. At the maximum, ξ4(t) and the α relaxation time τα are related by a power law. We also examine the individual contributions to g4(r,t), S4(q,t), and ξ4(t), as well as the corresponding order parameter Q(t) and generalized susceptibility χ4(t), arising from the self and distinct contributions to Q(t). These contributions elucidate key differences between domains of localized and delocalized particles.