🔬 physics

Chemical Disorder Suppresses Thermal Transport through Local Stress Redistribution in SrTiO3-Based High-Entropy Perovskites

通过将 MACE 教师模型蒸馏为用于大规模分子动力学的效率型 NEP，本研究揭示了局部应力无序（而非简单的质量或尺寸波动）是通过损失声子谱相干性来抑制高熵 SrTiO₃ 钙钛矿热传导的主要机制。

原作者： Yongheng Li, Maomao Liu, Hongwei Du, Baole Wei, Bin Wei, Bonan Zhu, Ziheng Lu, Yuanhua Lin

发布于 2026-09-07

📖 1 分钟阅读☕ 轻松阅读

原作者： Yongheng Li, Maomao Liu, Hongwei Du, Baole Wei, Bin Wei, Bonan Zhu, Ziheng Lu, Yuanhua Lin

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

热量是一种在材料中移动的能量形式，通常会将多余的热量从一处带到另一处。在热电领域，科学家们试图在阻止热量流动的同时，让电能自由流动。这是一个微妙的平衡过程。如果一种材料的导热性太好，它就无法用于将废热转化为有用的电能。几十年来，研究人员一直将一种被称为钙钛矿的晶体结构家族，特别是被称为钛酸锶的一种化合物，视为一个极具前景的起点。挑战在于如何在不破坏电流路径的情况下阻断热量的路径。一种强大的策略涉及“高熵”工程，这是一种通过在同一晶格位置混合多种不同类型的原子来进行的方法。其核心思想是，这种化学上的混乱应该能够散射携带热量的振动（即声子），并使其减速。然而，尽管科学家们知道这种混合确实有效，但他们并不清楚原子究竟是如何相互作用从而导致这种减速的，也不清楚哪种特定类型的原子紊乱起到了最大的作用。

一组研究人员现在利用一种新的计算机建模方法，揭开了这一谜团。他们专注于钛酸锶，这种材料中的原子排列在一种刚性的、重复的框架中。在这种结构中，一个中心的钛原子被六个氧原子包围，形成一个小八面体。这些八面体通过共用顶点来构建整个晶体。研究人员用包括铅、钙、钡和镧在内的其他元素混合物替换了外部的锶原子，创造了一个高度无序的化学环境。为了观察这个混乱混合物内部发生的情况，他们不能依赖标准的计算机模拟，因为对于如此复杂的系统，标准模拟往往过于缓慢或不够准确。相反，他们开发了一套复杂的流程。他们首先基于量子物理计算的数据，训练了一个高精度但运行缓慢的人工智能模型。接着，他们利用这个缓慢但精确的模型来教第二个更快速的模型如何预测原子的行为。这使得他们能够进行涉及数千个原子、跨越长时间周期的庞大模拟，捕捉决定热量传递的细微且转瞬即逝的运动。

结果揭示了不同原子如何影响热流的一个令人惊讶的层级关系。研究人员发现，仅仅向混合物中添加更多类型的原子并不总是会导致热导率的持续下降。虽然最初加入无序度会导致热传输急剧下降，但对于已经处于混乱状态的混合物再添加组分，其收益却递减。热导率在初期迅速下降，随后趋于平缓，并保持在热传导性能最差的理论极限值之上。这表明，仅仅存在许多不同的原子并不是问题的全部。随后，团队调查了为什么某些特定的原子在阻挡热量方面比其他原子更有效。他们发现，原子的尺寸比其重量更为重要。当他们用具有不同离子尺寸的钙取代锶时，热流被抑制的效果比使用质量大得多但尺寸接近锶的铅时要有效得多。这一发现挑战了长期以来的假设，即随机的原子质量差异是热散射的主要驱动力。

为了理解为什么尺寸如此关键，研究人员深入研究了晶体的局部环境。他们发现，原子尺寸的不匹配导致周围的氧框架发生了复杂的扭曲和变形。这种变形创造了一个随原子变化而变化的局部应力景观。研究人员开发了一种测量这种应力的方法，并发现局部应力的紊乱程度与热流的减少之间存在直接联系。局部应力的景观越混乱，携带热量的振动就会越失去相干性并发生散射。在模拟中，富钙混合物表现出的局部应力分布比富铅混合物更宽、更混乱。这种应力紊乱有效地缩短了热载体的寿命，将尖锐、有组织的波转变为弥散、散射的噪声。研究还表明，重铅原子确实引入了另一种效应，即产生了一种与热载体相互作用的软低频振动，但在阻挡热量方面，其效果不如尺寸不匹配引起的结构畸变。

这项工作为化学无序如何转化为物理特性提供了一张清晰的地图。它证明了在这些复杂晶体中，原子如何推挤和拉扯它们的邻居，从而创造出一个杂乱的局部应力景观，这比仅仅计算不同元素的数量或测量它们的质量差异更能准确地预测热性能。通过使用改进的人工智能方法来模拟这些相互作用，研究人员得以观察到热传输中那些不可见的机制。他们表明，抑制热量的关键不仅在于混合成分，还在于这些成分如何重塑晶体的局部机械环境。这一发现为设计更好的材料提供了新指南。未来的设计不再仅仅追求最大程度的化学多样性，而是可以专注于创造能最有效地散射热量的特定类型的局部应力和结构畸变。该研究证实，虽然高熵混合是一个强大的工具，但其成功取决于原子取代所带来的特定机械后果，这为开发能更高效管理热能的材料指明了路径。

技术摘要：化学紊乱通过局部应力重新分布抑制 SrTiO3 基高熵钙钛矿的热传输

问题陈述
高熵氧化物通过引入化学紊乱，为抑制热电材料中的晶格热导率提供了一条极具前景的途径。然而，高熵钛酸锶（SrTiO3 或 STO）中声子散射的原子尺度起源仍未得到解决。虽然熵工程已知可以降低热导率，但目前尚不清楚 A 位点紊乱场中的哪个特定组分——质量涨落、离子尺寸失配、局部晶格畸变或应力场——能最有效地通过 TiO6 骨架传递并抑制热传输。此外，由于机器学习力场在处理包含复杂化学空间（例如包含 La）的高熵氧化物时存在较大的拟合误差，且使用密度泛函理论（DFT）对用于热传输模拟的大型超胞和长轨迹进行计算的成本极高，构建准确的原子间势函数具有很大挑战性。

方法论
为了应对这些挑战，作者开发了一种结合机器学习势与分子动力学（MD）的两步工作流：

MACE 指导的 NEP 蒸馏：
- 教师模型： 一个预训练的 MACE（机器学习原子簇展开）模型，在涵盖 Sr/Ba/Ca/Pb/La A 位点化学空间的 615 个 DFT 标记构型上进行了微调。这种微调显著降低了预训练模型在能量、力和应力预测方面的误差。
- 蒸馏： 微调后的 MACE 模型作为“教师”，用于蒸馏出一个高效的神经进化势（NEP）。这一过程将 MACE 模型的精度转移到了一个能够进行大规模 MD 模拟的紧凑势函数中。
- 验证： 蒸馏后的 NEP 通过了针对 DFT 和直接 NEP 训练的验证。虽然直接训练的 NEP 在处理含 La 化学空间时表现不佳（力误差较高），但蒸馏后的 NEP 实现了显著降低的误差（能量误差为 2.81 meV/atom，力误差为 57.78 meV/Å），同时对于 9600 个原子的系统，其速度比 MACE 快 555 倍，从而实现了大规模模拟。
模拟与分析：
- 结构模型： 研究通过使用特殊准随机结构（SQS）和短程有序交换生成的 9600 原子超胞，考察了原始 STO 以及各种 A 位点取代组成的成分（单组分、三组分、四组分和五组分）。
- 热传输： 使用 GPUMD 中的平衡分子动力学（EMD）结合 Green-Kubo 公式计算晶格热导率。
- 描述符： 作者分析了多个描述符以与热导率建立关联：
  - A 位点质量涨落 ( $\delta_{mass}$ ) 和离子尺寸失配 ( $\delta_{size}$ )。
  - TiO6 键长畸变（均方根变化）。
  - 源自瞬时每原子维里张量的局部 von Mises 应力分布。
- 声子动力学： 计算了振动态密度（DOS）和动力学结构因子 $S(Q, E)$ ，以分析声子软化、展宽和群速度。

关键结果

热导率趋势： A 位点取代迅速降低了晶格热导率（原始 STO 值为 ~~8.93 W m⁻¹ K⁻¹）。然而，当系统接近高熵组成（四组分和五组分）时，这种降低趋于平缓，数值稳定在 ~3 W m⁻¹ K⁻¹ 附近，仍高于 Cahill-Pohl 极限（~~1.84 W m⁻¹ K⁻¹）。
掺杂特异性：
- Pb 取代： 在声子谱中引入了显著的低频类软模响应，这与重质量紊乱一致。
- Ca 取代： 尽管其质量涨落比 Pb 小，但产生了更强的全谱展宽和更大的结构异质性。
- 比较效率： Ca 引起的尺寸失配比 Pb 引起的质量紊乱能更有效地抑制热导率（例如，Sr $_{0.8}$ Ca $_{0.2}$ TiO $_3$ 的热导率低于 Sr $_{0.8}$ Pb $_{0.2}$ TiO $_3$ ）。
描述符相关性：
- 传统的基于成分的描述符（质量涨落和离子尺寸失配）与热导率仅表现出中度相关（ $r = -0.71$ 和 $-0.53$），且无法完全捕捉趋势，尤其是在区分 Pb 和 Ca 的效应方面。
- 局部应力描述符： 相对于有序 STO，A 位点 von Mises 应力分布主峰的偏移与热导率的相关性最强（ $r = -0.84$ ）。该描述符也与声子光谱相干性的丧失（ $S(Q, E)$ 中的峰值与背景比， $r = -0.87$ ）强相关。
- 机制： 结果表明，局部应力紊乱（而非单纯的质量或尺寸失配）驱动了声子光谱相干性的丧失和有效的类寿命展宽，这是热传输受抑的主要机制。
泛化性： 描述符的层级关系（应力 > 畸变 > 质量/尺寸）在含 Cd 系列（取代 La）中得到了重现，证实了这些发现并非仅与 La 相关的价态变化有关，而是 STO 中 A 位点紊乱的一个普遍特征。

意义与主张
本文提出了两个主要贡献：

方法论： 确立了 MACE 指导的 NEP 蒸馏是为难以直接拟合的复杂化学高熵氧化物生成准确、大规模原子间势函数的实用且有效途径。该方法平衡了热传输模拟所需的保真度与大规模超胞所需的计算效率。
物理见解： 确定了局部应力紊乱是低热导率的高熵钙钛矿中更优越的传输描述符。研究表明，化学紊乱不仅仅是质量和尺寸涨落的总和；通过 TiO6 网络实现的有限温度局部应力景观的重新分布，是连接 A 位点化学复杂性与非谐性声子谱及受抑热传输的关键因素。

作者得出结论，虽然高熵 A 位点紊乱强烈散射声子，但并未消除 STO 框架中所有的晶体热载流通道，因此可能需要进一步优化特定的掺杂比例（例如 Ca 对比 Pb）以接近理论最小热导率。

技术摘要：化学紊乱通过局部应力重新分布抑制 SrTiO3 基高熵钙钛矿的热传输

类似论文