Nonlinear phonon dispersion in disordered solids and non-Debye vibrational spectra

本研究证明，无序固体中的非线性声子色散源于无序诱导的介观长度尺度，并通过分析与模拟揭示，这种非线性软化与非声子振动均对玻色峰等非德拜反常现象有显著贡献，其相对重要性取决于材料的无序强度。

要点

想象一个固体物体，比如一块玻璃或一块金属，将其视为一个巨大而无形的管弦乐队。当你敲击它时，它不会静止不动，而是会振动。这些振动以波的形式在材料中传播，就像声波在空气中传播一样。在物理学中，我们将这些振动称为“声子”。

一个多世纪以来，科学家们一直使用一本名为德拜模型的经典规则书来预测这些振动的行为。该规则书指出：“如果你观察低频（缓慢）振动，它们会沿完美的直线传播，且振动的数量会以可预测的平滑方式增加。”

然而，当科学家观察无序固体（如窗玻璃，其没有晶体结构，与钻石不同）时，音乐变得混乱。振动不再遵循直线；它们发生弯曲，且振动的数量远远超过旧规则书的预测。这种额外的“噪声”在物理学中创造了一个著名的谜团，称为玻色峰。

长期以来，科学家们争论是什么导致了这种混乱。是因为波本身发生了弯曲（非线性色散）？还是因为无序产生了全新的、奇怪的振动类型，这些类型在完美晶体中并不存在（非声子模式）？

本文就像一部侦探故事，通过直接测量波并将两个嫌疑人区分开来，最终解决了这个案件。

侦探工作：一种新的聆听方式

最大的问题在于，在混乱的无序固体中，很难确切知道特定音高的波传播得有多快。这就像试图在拥挤嘈杂的房间里听清一把小提琴的声音。

作者开发了一种巧妙的新技术，称为**“强制波法”**。

• 类比：想象你在一个拥挤的房间里。与其等待有人开始唱歌，不如你轻轻推动房间里每个人，使其按照特定的波模式（就像体育场里的“人浪”）移动。然后你测量房间的反应。
• 通过在计算机上进行这种数学操作，他们能够迫使材料以特定模式振动，并精确测量随着音高升高，该波的速度如何变化。这使得他们能够以高精度描绘出波的“弯曲”路径。

发现：“隐藏的标尺”

他们发现，在无序固体中，波的弯曲并非仅仅因为碰到了原子的大小（就像在完美晶体中那样）。相反，它们开始弯曲是因为一个介观尺度。

• 类比：将完美晶体想象成排列完美的瓷砖网格。如果你走过它，当你踩到瓷砖边缘时就会绊倒。
• 在无序固体（如玻璃）中，没有瓷砖。然而，作者发现了一个“隐藏的标尺”（我们称之为 $\xi$），它比单个原子大得多。这个标尺代表了材料刚度开始随机波动的尺度。
• 发现：波平滑传播，直到它们变得足够大，能够“感知”到这个隐藏的标尺。一旦达到这个尺寸，它们就开始减速并弯曲（软化）。这个隐藏的标尺还控制了波被散射和损耗（衰减）的程度。玻璃越无序，这个隐藏的标尺就越大。

解开玻色峰之谜

一旦他们确切知道了波如何弯曲，他们就能仅根据这种弯曲计算出应该存在多少振动。然后，他们将这一计算结果与实际观测到的振动总数进行比较。

他们发现，“玻色峰”（额外的振动）实际上是两个不同来源的二重奏：

1. 弯曲的波（声子型）：波本身由于隐藏的标尺而发生弯曲，从而产生额外的振动。
2. 奇怪的局部抖动（非声子型）：由于材料混乱，其某些部分陷入了一种根本不传播的“抖动”。这些是局部的、被捕获的振动。

判决：

• 在高度无序的玻璃中（例如通过极快冷却制成的玻璃），“奇怪的局部抖动”是造成额外振动的主要元凶。
• 在稳定、现实实验室中的玻璃（我们实际使用的类型）中，“弯曲的波”和“奇怪的抖动”贡献几乎相等。

为何这很重要

该论文得出结论，长期以来，科学家们可能一直归咎于错误的因素。有些人认为额外的振动仅来自奇怪的局部抖动。另一些人则认为它们仅来自弯曲的波。

这项研究表明，两者都是关键角色。各自贡献的多少取决于玻璃的“混乱”程度。在我们实际实验室制造的玻璃中，你不能忽视其中任何一个；它们都显著地塑造了振动谱。

简而言之：这篇论文不仅找到了玻色峰，还绘制了一张新的地形图，表明玻璃的“混乱”创造了一个隐藏尺度，使波发生弯曲，而这种弯曲与被捕获的振动携手合作，共同创造了无序固体的独特声音。

技术摘要

技术摘要：无序固体中的非线性声子色散与非德拜振动谱

问题陈述
虽然晶体固体可由德拜经典模型很好地描述——其中低频声学声子遵循线性色散关系 $\omega(k) \propto k$，且振动态密度（VDoS）按 $D(\omega) \propto \omega^2$ 标度——但无序固体（玻璃）表现出显著的偏差，即非德拜异常。其中最突出的是“玻色峰”（BP），即振动模式超出德拜预测的过剩。这些异常的起源仍是争论的焦点。提出了两种主要机制：(1) 结构无序导致的声子色散关系 $\omega(k)$ 的软化非线性；(2) 无序结构固有的非声子型、准局域化振动模式的存在。当前理解的一个关键缺口在于无法定量分离并确定这两种机制对非德拜行为及玻色峰出现的相对贡献。此外，控制非晶体固体中声子色散和波衰减的无序诱导长度尺度的性质尚未被完全理解。

方法论
为解决这些挑战，作者对三类不同的三维无序固体进行了大规模计算机模拟：

1. 无序弹性网络：由松弛的胡克弹簧组成，具有不同的平均配位数（$z$），范围从临界麦克斯韦阈值（$z_c=6$）到远高于该阈值。
2. 多分散逆幂律（PIPL）玻璃：排斥系统，经历不同的热历史，由系统脱离平衡时的母体温度 $T_p$ 量化。
3. 二元玻璃：包括排斥性二元混合物以及具有可调吸引力（由参数 $r_c$ 控制）的“粘性球”玻璃。

研究采用了一种名为“强制波法”的新型计算技术来计算无序系统中的声子色散关系 $\omega(k)$。与晶体不同（晶体中声子是海森矩阵的精确本征函数），无序性使得无法为给定的波矢 $k$ 唯一地识别 $\omega$。强制波法通过对粒子施加类波力场并计算线性响应位移，为每个允许的波矢分配一个唯一的频率。这使得能够提取 $\omega(k)$ 及其相关的群速度。

作者还采用了核多项式方法（KPM）来计算超大系统尺寸（$N \sim 10^6$）下的 VDoS $D(\omega)$，从而能够解析低频尾部并分离声子型和非声子型贡献。为统一不同模型中无序的描述，计算了每个系统的无量纲无序量化器 $\chi$，定义为介观尺度上剪切模量涨落与其均值的比值。

主要贡献与结果

1. 发现介观无序诱导长度尺度（$\xi$）：
   作者证明，无序固体中声子色散关系的非线性软化并非由微观原子间距离 $a_0$（如晶体中那样）控制，而是由大得多的、涌现的介观长度尺度 $\xi \gg a_0$ 控制。色散关系的形式特征为 $\omega(k) = f(k\xi)ck$，其中一阶非线性项按 $\omega(k) \approx ck - c\Upsilon \xi^2 k^3$ 标度（其中 $\Upsilon = 1/96$）。

   • 标度关系：发现 $\xi$ 随无序程度标度变化：对于接近去阻塞转变的弹性网络，$\xi \sim 1/\sqrt{z-z_c}$；对于玻璃，$\xi$ 随 $T_p$ 和 $r_c$ 的增加而增加。
   • 统一无序度量：关键在于，作者表明 $\xi/a_0$ 与无量纲无序量化器 $\chi$ 呈线性标度关系（$\xi/a_0 \sim \chi$），从而统一了不同无序控制参数的影响。

2. 与波衰减的联系：
   相同的长度尺度 $\xi$ 被证明控制着波衰减率 $\Gamma(k)$ 的瑞利散射。作者推导出了小 $k$ 下的关系式 $\Gamma(k) \sim \xi^2 k^4$，这与非均匀弹性理论（HET）和瑞利 - 克莱门斯关系一致。这建立了非线性色散与波衰减之间的直接联系，两者均由同一介观无序诱导长度尺度控制。

3. 非德拜谱的定量分解：
   该论文定量分离了非线性声子色散（声子型）和非声子型振动对 VDoS 的贡献。

   • 声子型贡献：利用测得的 $\omega(k)$，作者计算了声子型 VDoS $D_{ph}(\omega)$，其偏离德拜定律的形式为 $D_{ph}(\omega) \approx A_D \omega^2 + A_w \omega^4 + \dots$。系数 $A_w$ 取决于 $\xi$。
   • 非声子型贡献：对于玻璃，存在无间隙的非声子型模式分布，具有通用尾部 $D_G(\omega) \approx A_g \omega^4$。
   • 相对重要性：比率 $A_g/A_w$ 决定了哪种机制主导非德拜行为的出现。结果表明，该比率对无序程度高度敏感。在高度无序（不稳定）的玻璃中，非声子型模式占主导地位（$A_g/A_w > 1$）。然而，在稳定、现实的实验室玻璃（低无序）中，来自非线性色散的声子型贡献往往占主导地位（$A_g/A_w < 1$）。这挑战了此前将 $\omega^4$ 尾部完全归因于非声子型模式的解释。

4. 玻色峰的组成：
   通过分析玻色峰频率 $\omega_{BP}$ 处的累积 VDoS，作者定义了一个比率 $R(\omega_{BP})$，代表总过剩模式中非声子型振动的分数。

   • 结果：发现 $R(\omega_{BP})$ 是接近 1 的显著分数（范围从 $\sim 0.2$ 到 $\sim 0.5$），跨越广泛的玻璃无序范围。
   • 结论：玻色峰并非仅由软化声子填充（如某些近期理论所建议），也非仅由非声子型模式填充。相反，它源于非线性声子软化与非声子型准局域化激发之间的显著混合。

意义
该论文声称通过建立一个全面的定量框架，将波色散、波衰减和振动谱通过单一介观长度尺度 $\xi$ 联系起来，从而在理解无序固体方面取得了“基础性进展”。通过解析声子型和非声子型机制的相对贡献，该工作阐明了玻色峰和非德拜异常的物理起源。研究结果表明，玻色峰的组成并非普适的，而是强烈依赖于玻璃的具体无序强度和热历史。所开发的分析框架和强制波法为实验振动谱提供了可检验的预测，有望解决关于非晶材料中低频激发性质的长期争议。