通过能量最小化探究热软球玻璃的临界慢化

使用混合分子动力学/SWAP 蒙特卡洛（MD/SMC）模拟，我们表明，随着样品变得越来越充分平衡，软球玻璃的 FIRE 能量最小化的终端弛豫时间 $\tau$ 表现出热起始现象。尽管随着平衡的进行以及阻塞密度 $\phi_{\rm J}$ 通过热起始而增加，$\tau(\phi)$ 可以减小几个数量级，但它始终标度为 $\tau(\phi) \sim (\phi_{\rm J} - \phi)^{-2} \sim [Z_{\rm iso} - Z_{\rm ms}(\tau)]^{-2}$，其中 $\phi_{\rm J}$ 是阻塞密度，$Z_{\rm ms}(\tau)$ 是满足最小局部力学稳定性判据（$Z \geq d+1$）的粒子的平均配位数，位于系统所遇最终势能景观（PEL）子盆地的顶部。这种标度允许我们将作为 $\phi$ 函数绘制时看起来非常不同的 $\tau$ 数据集坍缩，并解决一个密切相关的问题：在恒定 $\phi$ 和 $T$ 下，随着玻璃老化，致密热玻璃最通常占据的 PEL 盆地的特征如何演变？