Local micromechanics in a mean-field model of glasses reveal key properties of its non-equilibrium RSB phase

技术摘要：平均场玻璃模型中的局部微观力学

问题与动机
玻璃是一种以无序和缺乏长程有序为特征的非平衡态固体。虽然连续介质尺度的性质（如弹性）已得到充分理解，但微观起源的非平衡现象——例如非德拜（non-Debye）振动谱、屈服转变以及机械记忆——仍然具有挑战性。最近在有限维计算机玻璃研究中的进展识别出了“准局域化”低频模式以及振动态密度（VDoS）中普遍存在的 $\sim\omega^4$ 尾部，$D_G(\omega)$。

最近提出的一种无序哈密顿平均场模型 [47, 48] 展示了统一这些玻璃特性的潜力。该模型在特定的参数范围内具有平衡态、零温复制对称破缺（RSB）转变。在接近 RSB 转变的复制对称（RS）相中，其解可以用一个有效、自洽的随机势能来描述，该势能作用于不耦合的自由度。然而，该模型在进入非平衡 RSB 相（紧随瞬时淬火之后）后的行为，及其与有效势能描述之间的联系，仍需进一步研究。具体而言，需要为一个缺乏空间结构的完全连接模型定义“微观力学”响应函数，并建立局部可观测物理量与全局易受性之间的精确关系。

方法论
作者分析了一个由以下哈密顿量定义的平均场玻璃模型：
$$H = \frac{1}{2}\sum_{i=1}^N k_i x_i^2 + \frac{1}{4!}\sum_{i=1}^N x_i^4 + \sum_{i<j} J_{ij}x_i x_j - f(t)\sum_{i=1}^N x_i$$
其中 $k_i$ 是从某种分布中抽取的初始刚度（在所研究的案例中，$p(k) \propto k^{-1}$ 在 $[0,1]$ 上），而 $J_{ij}$ 是高斯随机耦合。系统在非平衡 RSB 相中通过最小化 $H$（瞬时淬火）进行研究，而不进行热化。

其核心方法论创新是在这种平均场语境下定义了力单极子响应（force monopole response）。类比于有限维玻璃中的力偶极子（它们屏蔽了长波声子），作者施加一个局部力单极子 $f_i$（作用在第 $i$ 个振子上的单位力），并测量线性响应 $z_i = M^{-1} \cdot f_i$，其中 $M$ 是 Hessian 矩阵。他们定义了局部单极子刚度 $\kappa_i$：
$$\kappa_i = \frac{f_i \cdot M^{-1} \cdot f_i}{f_i \cdot M^{-2} \cdot f_i} = \frac{(M^{-1})_{ii}}{(M^{-2})_{ii}}$$
利用随机矩阵理论（RMT）技术，作者推导出了一个将这些局部刚度与全局线性（$\chi$）和自旋玻璃（$\chi_{SG}$）易受性联系起来的精确平均场恒等式。随后，他们利用数值模拟验证了这些关系，重建了完整的 VDoS，并基于边际稳定性判据识别出了一个“稳定扇区（stabilized sector）”的振子。

主要贡献与结果

1. 精确平均场恒等式： 作者推导出了一个将局部单极子刚度 $\kappa_i$ 与全局易受性联系起来的严格关系：
   $$\kappa_i = \frac{(k_i - J^2\chi) + \frac{1}{2}x_i^2}{1 + J^2\chi_{SG}}$$
   该恒等式表明，完全局部的性质是由全局量决定的，这是模型完全连接特性的体现。至关重要的是，为了使 $\kappa_i$ 有限，自旋玻璃易受性 $\chi_{SG} = \int \omega^{-4} D_G(\omega) d\omega$ 必须是有限的。这要求 VDoS 具有伪能隙（pseudogap），其尾部指数大于 3，这与模型观察到的 $\sim\omega^4$ 尾部一致。

2. 微观力学局域化与伪能隙： 数值分析显示，单极子响应 $z_i$ 是强局域化的（参与数 $\langle N_e \rangle \approx 4$，与系统尺寸 $N$ 无关），这镜像了有限维玻璃中的准局域化模式。此外，单极子刚度分布 $p(\kappa)$ 具有伪能隙（当 $\kappa \to 0$ 时消失），表明了边际稳定性。

3. 振动谱重建： 作者展示了仅通过单极子刚度 $\{\kappa_i\}$ 的统计特性即可重建完整的非声子 VDoS $D_G(\omega)$。通过利用推导出的刚度统计特性求解解析函数的自洽方程，他们恢复出的 $D_G(\omega)$ 与直接进行 Hessian 矩阵对角化得到的结果具有极佳的定量一致性。

4. 特征尺度与玻色峰： 研究识别出了一个控制非声子 VDoS 的单一特征频率尺度。

   • 作者定义了一个“稳定扇区”的振子，其中初始刚度 $k_i < J^2\chi$。这些振子依赖于非谐位移项（$\frac{1}{2}x_i^2$）来维持正刚度。
   • 该稳定扇区的平均刚度 $\langle \kappa \rangle_{st}$ 随内在玻色峰频率的平方缩放：$\langle \kappa \rangle_{st} \sim \omega_p^2$。
   • 此外，由力涨落 $\Delta f$ 推导出的有效随机势能的特征频率尺度满足 $\omega_p \sim (\Delta f)^{1/3}$。

5. 介观与全局弹性： 介观刚度（由 $\langle \kappa \rangle_{st}$ 表征）与全局剪切模量 $G$ 的比值随无序度 $J$ 的变化显著。具体而言，$\langle \kappa \rangle_{st}$ 的变化范围比 $G$ 更宽，这一行为与在有限维计算机玻璃中的观察结果一致。

意义与主张
本文声称，该平均场模型的统计微观力学（特别是力单极子响应）为平衡态 RS 相中由有效随机势能描述的集体自由度提供了具体的实现。研究结果表明：

• 有效势能（针对平衡态 RS 相推导）的数学结构在深入非平衡 RSB 相时仍然有效。
• “稳定扇区”的振子（即 $k_i < J^2\chi$ 的振子）构成了负责 $\omega^4$ 尾部和玻色峰的软玻璃模式的物理体现。
• 该模型提供了一个统一的框架，连接了平均场理论与有限维玻璃物理，其中局部微观力学可观测量（单极子刚度）直接编码了全局谱性质和易受性极限。

作者总结道，这些结果强化了如下理论建议：该模型的能量景观本质上是玻璃态的，且有效势能图像为理解其非平衡性质提供了一种统一的描述。他们指出，未来的工作应探索稀疏连接性、非谐特性以及强驱动动力学。