🔬 材料科学

Where heat flows and where to change it: gauge-free sensitivity maps of thermal conductivity in disordered solids

本文开发并比较了无序固体热导率的两种局部描述方法，揭示了尽管由于局部软化效应，依赖于规范的原子投影与无规范的质量/刚度敏感度图在初始阶段呈负相关，但后者为优化材料设计以及理解非晶态和副晶态硅中的热输运限制提供了一种稳健且精确的指导。

原作者： D. A. Drabold, C. Ugwumadu

发布于 2026-11-06✓ Author reviewed ⓘ

📖 1 分钟阅读☕ 轻松阅读

CC BY 4.0

原作者： D. A. Drabold, C. Ugwumadu

原始论文采用 CC BY 4.0 许可（http://creativecommons.org/licenses/by/4.0/）。 ✨ 这是对下方论文的AI生成解释。它不是由作者撰写的。如需技术准确性，请参阅原始论文。阅读完整免责声明

热量是固体材料内部不断的旅行者，从热点向冷点移动；但在非晶态固体（如玻璃或用于计算机芯片的硅）这种混乱的世界里，这场旅程是混沌的。与原子位于完美重复晶格中、允许热量沿直线疾驰的晶体不同，非晶态固体是由没有长程有序性的原子构成的纠缠乱团。几十年来，科学家们已经知道如何计算一种材料所能携带的热量总量，这个数值被称为热导率。然而，对于试图设计更好电子设备或更高效能量转换器的工程师来说，仅知道总量是不够的。他们需要确切地知道热量在材料内部的具体流向，更重要的是，知道对材料结构的微小改变会在何处产生最大的影响。如果你想阻止热量从设备中逃逸，或者帮助它更快地逃逸，你必须知道哪些特定的原子在做“重活”，而哪些仅仅是旁观者。

这个问题之所以难以回答，是因为在如何定义单个原子在群体中的能量这一问题上，存在一个根本性的歧义。在计算热流时，研究人员必须决定如何将系统的总能量分配给各个原子。这种选择是任意的，就像在不清楚每个人具体吃了多少钱的情况下，决定如何分摊账单一样。分配账单的不同方式会导致不同的热流分布图，尽管总热量保持不变。这让科学家们感到困惑：他们绘制的热流图描述的是材料真实的物理特性，还是仅仅是他们用来进行计算的数学记账方式？

为了消除这种困惑，俄亥俄大学和洛斯阿拉莫斯国家实验室的研究人员将注意力转向了非晶硅——这种材料结构简单，但其复杂性足以代表广泛的无序固体。他们开发了两种截然不同的观察问题的方法。第一种方法是直接绘制热流图，描绘出材料中哪些区域携带了最多的能量。他们测试了当改变分配能量的任意方式时，这张地图会发生多大的变化。他们发现，虽然在单原子层面上地图看起来略有不同，但如果从大约几个原子宽度的微小区域来看，它会变得非常稳定。这证实了热流图是材料的一个真实物理特征，而非仅仅是数学上的人工产物。

第二种，也更令人惊讶的方法提出了一个不同的问题：如果你改变单个原子的质量，或者改变它与其相邻原子之间的键合刚度，总热流会如何变化？这个问题没有歧义，因为它不需要拆分能量；它只是测量系统对局部微调的敏感度。研究人员在他们的模型中为数千个原子计算了这种敏感度图。他们发现了一个惊人的响应符号反转现象：在热量流动最剧烈的区域，增加原子的质量往往会增加总热流；相反，在热量流动极少的区域，增加原子的质量实际上会降低总热流。

这种看似违反直觉的结果是由材料局部的“柔软度”驱动的。在非晶硅中，一些原子被松散地固定在原地，以较大的振幅振动，而另一些原子则受到严格约束，振动幅度很小。研究人员发现，热量流动最少的区域正是那些原子振动剧烈、松散且固定不牢的区域。在这些柔软的部位，由于原子本身运动幅度过大，增加质量会破坏允许能量传递的振动混合机制，从而有效地阻碍了流动。在热量自由流动的坚硬区域，原子被紧紧束缚，增加质量往往会提高导热性。然而，研究人员指出，虽然响应的符号与热流方向可靠地相反，但敏感度的量级并不总是与热流完全匹配；两者仅在平均意义上一致。

这项研究表明，这种柔软度与热流之间的关系是理解该材料的关键。研究人员展示了，如果考虑到一个原子如何振动，这种奇怪的敏感度符号反转现象就会基本消失。那些柔软、大振幅的区域充当了瓶颈，并不是因为它们直接阻挡了热量，而是因为它们通过混合振动模式的方式，阻止了能量的高效移动。当研究人员在更大的模型和不同的温度下进行测试时，他们发现这种效应在室温下最为显著，并且随着材料尺寸的增加而变得更加明显，这使得长波长振动能够参与其中。

研究结果对于设计材料具有直接的实际价值。敏感度图可以作为工程师的指南。例如，如果你想降低热流（以保持热电装置的高效），你不应该只是随机添加重原子。相反，你应该将它们放置在图中显示出最大敏感度的柔软、高振幅区域。研究人员通过模拟用较重的锗原子替换硅原子测试了这一点。他们发现，将这些重原子放置在特定预测的柔软点位，比随机放置或放置在坚硬的热流区域更能显著降低热流。

研究还观察了无序非晶区域与有序晶体之间的边界。他们发现，跨越该边界的热流几乎完全受无序侧控制。晶体侧过于坚硬且均匀，改变其原子对总热传递的影响微乎其微。瓶颈在于非晶区域，特别是进入无序材料仅几个原子层深处的那些柔软点位。这意味着，要改善或降低此类界面的热性能，必须专注于无序的一侧，而不是晶体一侧。

最终，这项工作提供了一种清晰、明确的方法来观察无序固体内部的情况。它区分了“热量在哪里流动”与“哪里发生改变会产生影响”这两个问题。这两个答案往往是相反的：承载热量最多的地方，是你能做出最小改变的地方；而发生微小改变能产生最大影响的地方，恰恰是热量流动最少的地方。通过关注原子局部的柔软度和振动幅度，科学家现在可以准确预测在哪里寻找控制材料热特性的方法。

技术摘要：无规范敏感度图谱在无序固体热导率中的应用

问题陈述
设计热管理材料不仅需要了解总热导率 ( $\kappa$ )，还需要识别热量在无序网络中何处流动，以及局部的改性会在何处改变这种传输。虽然 Allen–Feldman (AF) 谐波理论通过将热载体视为弥散子 (diffusons)、传播子 (propagons) 和定位子 (locons)，成功预测了玻璃的总 $\kappa$ ，但它仅提供一个全局数值。以往试图实现该物理量局域化的尝试（如位点投影热导率，SPTC, $\zeta(\ell)$ ）存在根本性的“规范 (gauge)”歧义问题：分配给单个原子的值取决于势能如何在原子间进行划分。这种任意性引发了一个疑问：局域图谱描述的是物理传输过程，还是仅仅是记账惯例？此外，在不重新计算整个系统的情况下，如何预测局部结构变化（例如质量取代或键强硬化）对总导热性的影响，目前仍不明确。

方法论
作者利用谐波 AF 理论，开发并测试了两种互补的非晶硅 (a-Si) 热传输局域描述方法。

位点投影热导率 ( $\zeta$ )： 作者分析了标准 $\kappa$ 在原子上的投影的鲁棒性。他们量化了 $\zeta$ 对局部能量选择（即“规范”）的依赖性，对比了标准的 AF 规范与基于机器学习势能（ACE, GAP-18）和 Keating 力场得到的位点能量规范。
无规范敏感度图谱 ( $\kappa_m$ 和 $s_b$ )： 为了消除规范歧义，作者推导了 $\kappa$ $κ$ 对局部参数（每个原子的质量 $m_\ell$ $m_{ℓ}$ 以及每条键的刚度）的精确梯度。他们定义了质量敏感度 $\kappa_m(\ell) = 2m_\ell \partial \kappa / \partial m_\ell$ $κ_{m} (ℓ) = 2 m_{ℓ} \partial κ / \partial m_{ℓ}$ 。这种方法无需进行能量划分。
- 精确分解： 质量敏感度图谱被精确分解为三个项：电流项（与 $\zeta$ 成正比）、特征值偏移项（经敏感度加权的局部态密度）以及旋转项（通过模式混合重新分布敏感度）。
- 局域软度代理： 作者利用每个原子的谐波振动振幅 ( $\delta R$ )（由动力学矩阵的伪逆导出）作为局域软度的结构代理。
计算模型： 研究采用了连续随机网络模型（包含 4096、512 和 32,768 个原子的 a-Si）以及准晶体模型（嵌入基质中的晶粒）。计算使用了 ACE、GAP-18 和 Keating 势能。对于较小系统，作者采用密集对角化法；对于较大系统（高达 10,000 个原子），则采用切比雪夫多项式展开法。

核心贡献

解决规范歧义： 本文证明，虽然 $\zeta(\ell)$ 在单个位点上具有中度的规范依赖性（不确定性约为空间跨度的 23%），但在 4 Å 半径范围内，该值变得稳健（不确定性约为 9.5%）。在超过几个埃的尺度上，位点的排序几乎是不变的。
推导精确敏感度： 作者首次提供了 AF 导热系数对局部质量和键的精确梯度，并建立了一个必须满足的欧拉求和规则 (Euler sum rule)。
发现反相关性： 一个核心发现是，热流图 ( $\zeta$ ) 与质量敏感度图 ( $\kappa_m$ ) 之间存在强烈的反相关关系。热量流向极少的区域（软、大振幅区域）表现出较大的正 $\kappa_m$ （增加质量会降低 $\kappa$ ），而热量流动最剧烈的区域则表现出负值或较小的 $\kappa_m$ （增加质量会提高 $\kappa$ ）。
机制识别： 这种反相关性是由局域软度 ( $\delta R$ ) 介导的。敏感度分解中的旋转项涉及相距 1–2 meV 的振动模式混合，这是该效应的主要驱动力。
应用于界面与设计： 这些图谱被应用于非晶/晶体界面及准晶硅，研究表明，无序基质而非晶粒控制着横向传输响应。

结果

$\zeta$ 的鲁棒性： 在不同势能（ACE, GAP-18, Keating）和网络尺寸下，热流图 $\zeta$ 在 4 Å 尺度上是一个稳健的结构描述符。不同规范或势能之间的 Spearman 相关系数很高（0.86–0.98）。
反相关机制： 在 a-Si 中，热流最高的 12% 体积与质量敏感度最高的 12% 体积之间的重叠部分仅为 0.27%。
- 由软度介导： 在室温 (300 K) 下，以局域振动振幅 $\delta R$ 为条件进行处理后，反相关性基本消失。这证实了“软”区域（大 $\delta R$ ）携带的热量很少，且通过混合附近的振动模式，改变那里的质量会显著改变 $\kappa$ 。
- 温度与尺寸依赖性： 反相关性随系统尺寸增大而增强（出现更多低频模式），并且在低温下最为显著。在 ~100 K 以下，即使在以 $\delta R$ 为条件后，仍存在残余的反相关性，这表明长波传播子 (propagons) 贡献了一个无法仅通过局域软度捕捉的部分。
准晶硅： 在含有纳米晶粒的模型中，晶粒单位原子的导热能力高于非晶基质。然而，敏感度图显示，非晶基质控制着对质量变化的响应。 $\kappa_m$ 的符号从正值（基质）变为负值（晶粒内部），但响应的量级受无序部分的控制。
设计启示：
- 取代物放置： 对于分散的重元素取代（如 Si 中的 Ge），敏感度图是一个有效的指南。将取代物放置在具有高正 $\kappa_m$ 的区域（软区域）比随机放置或放置在硬区域能更显著地降低 $\kappa$ 。
- 紧凑团簇： 对于紧凑的取代团簇，该图谱的预测能力会下降；其效果主要取决于取代原子的总数，表现为一个单一的重缺陷。
- 界面： 在 a-Si/c-Si 界面处，敏感度图识别出非晶侧（距离界面 4–6 Å 处）是局部扰动最能影响横向导热性的区域，而位点投影 $\zeta$ 在整个系列电阻层中保持均匀。

意义与主张
本文通过将关注点从“热量如何分配”（具有规范依赖性）转向“热传输如何响应变化”（无规范依赖性），解决了长期以来关于局域热传输描述的歧义。作者认为，敏感度图为设计问题提供了一个建设性的答案：哪里发生变化才会有意义？

这项工作确立了在无序固体中，热量携带最少的区域恰恰是热导率对局部扰动最敏感的区域。这一看似反直觉的结果可以通过软区域中振动模式的混合来解释。作者将敏感度图定位为一种实用的材料设计工具，特别是在优化热电材料（低 $\kappa$ ）或纳米级器件的热管理方面，并指出局域振动振幅 ( $\delta R$ ) 可以作为完整敏感度图的廉价且高效的计算代理。本研究受限于谐波网络和 $\Gamma$ 点，作者指出，物理寿命（非谐性）可能需要更大的胞单元来解析敏感度图而不产生人工展宽，但其在谐波框架内的定性反相关性和设计原则被认为是稳健的发现。

技术摘要：无规范敏感度图谱在无序固体热导率中的应用

类似论文