A simple model for dynamic heterogeneity in glass-forming liquids

Rajib K. Pandit,Horacio E. Castillo
DOI: https://doi.org/10.1103/PhysRevLett.131.218202
2023-10-21
Abstract:Liquids near the glass transition exhibit dynamical heterogeneity, i.e. local relaxation rates fluctuate strongly over space and time. Here we introduce a simple continuum model that allows for quantitative predictions for the correlators describing these fluctuations. We find remarkable agreement of the model predictions for the dynamic susceptibility $\chi_4(t)$ with numerical results for a binary hard-sphere liquid and for a Kob-Andersen Lennard-Jones mixture. Under this model, the lifetime $\tau_{\rm ex}$ of the heterogeneities has little effect on the position $t = t_4 \sim \tau_{\alpha}$ of the peak of $\chi_4(t)$, but it controls the decay of $\chi_4(t)$ after the peak, and we show how to estimate it from this decay.
Soft Condensed Matter,Disordered Systems and Neural Networks
What problem does this paper attempt to address?
The problem that this paper attempts to solve is about the quantitative description and understanding of dynamic heterogeneity in liquids near the glass - transition. Specifically, the authors aim to explain and predict the correlation functions that describe the fluctuations of these dynamic heterogeneities by introducing a simple continuous model. The following is a detailed summary of the key problems that the paper attempts to solve: 1. **Quantification of Dynamic Heterogeneity**: - Liquids exhibit strong dynamic heterogeneity when approaching the glass - transition, that is, the local relaxation rate fluctuates violently in space and time. The author hopes that by introducing a simple model, it is possible to quantitatively predict the correlation functions that describe these fluctuations. 2. **Model Validation**: - The authors use this model to predict the dynamic susceptibility \(\chi_4(t)\) and compare it with the numerical simulation results of binary hard - sphere liquids and Kob - Andersen Lennard - Jones mixtures, showing the consistency between the model prediction and the numerical results. 3. **Influence of Heterogeneity Lifetime**: - The model shows that the lifetime of heterogeneity \(\tau_{\text{ex}}\) has a small influence on the peak position \(t = t_4\sim\tau_\alpha\) of \(\chi_4(t)\), but controls the decay of \(\chi_4(t)\) after the peak. The author shows how to estimate \(\tau_{\text{ex}}\) from this decay. 4. **Intuitive Description of Dynamic Heterogeneity**: - The model is directly based on the intuitive description of dynamic heterogeneity, using the local relaxation rate \(\gamma(\vec{r},t)\equiv1/\tau_{\vec{r}}(t)\) as the basic variable, thus expressing the multi - point correlation function as a two - point - and - a - half correlation function \(s(\vec{q},t,t')\), which helps to better understand the size and lifetime of heterogeneous regions. 5. **Connection between Numerical Simulation and Experimental Results**: - The model provides a more direct way to connect the numerical simulation results with the experimental results, especially in the study of dynamic heterogeneity. 6. **New Simplified Method**: - The model provides a new and simpler method to extract the estimated value of heterogeneity lifetime from the numerical simulation data. ### Summary of Mathematical Formulas - The dynamic susceptibility \(\chi_4(t)\) is defined as: \[ \chi_4(t)=\lim_{q\rightarrow0}S_4(\vec{q},t) \] - The local relaxation function \(C_{\vec{r}}(t)\) is defined as: \[ C_{\vec{r}}(t)\equiv\rho^{-1}\sum_{n = 1}^N w_n(t)\delta(\vec{r}_n(0)-\vec{r}) \] - The four - point dynamic structure factor \(S_4(\vec{q},t)\) is defined as: \[ S_4(\vec{q},t)\equiv\rho^2N^{-1}\int d\vec{r}e^{-i\vec{q}\cdot\vec{r}}\langle[C_{\vec{r}}(t)-C(t)][C_{\vec{0}}(t)-C(t)]\rangle \] - The collective contribution \(C_{\text{coll}}^{\vec{r}}(t)\) is defined as: \[ C_{\text{coll}}^{\vec{r}}(t)=\rho^{-1}\rho(\vec{r},0)C[\phi_{\vec{r}}(t)-\phi_{\vec{r}}(0)] \] - The exchange time \(\tau_{\text{ex}}\) is defined as: \[ \tau_{\text{ex}}\equiv\frac{\int_0^\infty t\chi_{\phi}^2(t)dt}{\int_0^\infty}