Low and Anisotropic Thermal Conductivity in Mixed-Valent Sn$_2$S$_3$

这项研究表明，混合价态 Sn$_2$S$_3$ 沿 c 轴表现出低且各向异性的热导率，这是由孤对电子库仑相互作用诱发的弱键合 Sn(II) 原子的非谐性“碰撞”振动所驱动的。

要点

想象一种像热量“交通堵塞”一样，能阻止热量轻易流过的材料。这就是研究人员通过研究 Sn₂S₃（硫化锡）这一化合物来揭示其为何如此擅长阻挡热量的故事。

以下是他们研究结果的拆解，使用了简单的类比：

1. 材料：热量的一单行道

不要把 Sn₂S₃ 仅仅看作一个实心块，而要把它看作一捆绑在一起的吸管或面条。

• 强方向（面条）： 如果你尝试沿着面条的长度（b轴）推动热量，它移动得非常快。这里的原子紧密相连，就像一条润滑良好的高速公路。
• 弱方向（间隙）： 如果你尝试横向（a轴和c轴）推动热量，它就会卡住。面条之间存在间隙，就像碗里面条之间的空隙一样。热量很难跳过这些间隙。
• 结果： 这种材料具有高度的“各向异性”，这意味着它对热量的处理方式取决于你试图发送热量的方向。它就像一条单行道，交通在某个方向流动迅速，但在其他方向却陷入了瘫痪。

2. “摇摆”原子：松动的螺丝

在这个结构内部，有两种类型的锡原子：Sn(IV) 和 Sn(II)。

• Sn(IV) 就像一颗紧紧拧进墙里的螺丝。它纹丝不动。
• Sn(II) 就像一颗头部摇晃的松动螺丝。它拥有“孤对”电子（可以把它们想象成隐形的、具有排斥力的气球），这些气球会推挤它的邻居。
• 摇摆现象： 由于这些排斥性的气球的存在，Sn(II) 原子并不能牢固地固定在原位。它们在自己的小笼子里摇晃、剧烈且混乱地振动。研究人员称这些原子为**“摇摆者”（rattlers）**。

3. 摇摆如何阻挡热量

通常，热量在固体中的传播就像波浪在体育场人群中传递（人们依次起立又坐下）。这被称为“声学声子”。

• 干扰： 当这些“松动的螺丝”（Sn(II)）开始摇摆时，它们就像体育场里的人突然开始随机地跳来跳去。这种混乱状态会散射有序的热波，将其打碎并阻止其流动。
• 惊喜之处： 研究人员发现，这些摇摆原子创造了非常缓慢、平坦的振动（低频光学声子）。通常，科学家认为只有快速、有序的波才能携带热量。但在这种材料中，这些混乱的、摇摆的振动实际上携带了惊人的热量（沿快速方向约占 63%），这是一个罕见且有趣的发现。

4. 温度的转折

通常，随着物体变热，热量的流动方式也会发生变化。

• 论文的发现： 在大多数材料中，热量流动会随着温度升高而呈现可预测的下降。但在 Sn₂S₃ 中，无论温度如何变化，热量流动都保持着令人惊讶的稳定且处于低水平。这是因为“摇摆”机制在散射热量方面极其有效，以至于无论你增加多少能量，交通堵塞依然存在。

总结

论文得出结论，Sn₂S₃ 是一种“混合价态”材料（意味着它具有两种不同状态的原子），其中 Sn(II) 原子就像硬盒子里摇晃的弹珠。由于电子排斥作用，这些弹珠会在盒子里剧烈摇摆，创造出一个散射热波的混乱环境。这使得该材料非常擅长阻挡热量，尤其是在特定方向上，为寻找能够保持冷却或高效管理热量的材料提供了一个新的蓝图。

技术摘要

技术摘要：混合价态 $\text{Sn}_2\text{S}_3$ 的低热导率与各向异性

问题陈述
寻找具有本质低晶格热导率 ($\kappa_l$) 的材料对于推进热电能量转换和热管理至关重要。虽然原子“摇摆”（rattling）和孤对电子的概念已被已知可降低笼状结构中的热导率，但在混合价态化合物中支配这些现象的具体机制仍未得到完全理解。本研究重点关注硫化锡 ($\text{Sn}_2\text{S}_3$)，这是一种具有 $\text{Sn(II)}$ 和 $\text{Sn(IV)}$ 两种氧化态的混合价态化合物。尽管已有关于其低热导率的实验测量结果，但其背后的摇摆机制、其准一维结构对热各向异性的影响，以及光学声子对热输运的具体贡献，仍需要更深入的理论阐明。

方法论
作者采用全面的第一性原理计算框架来研究 $\text{Sn}_2\text{S}_3$ 的晶格动力学和声子输运：

• 电子结构： 使用基于 PBE-GGA 泛函和 Grimme DFT-D1 校正（以考虑准一维结构中的范德华相互作用）的 VASP（Vienna Ab initio Simulation Package）进行密度泛函理论 (DFT) 计算。
• 晶格动力学： 利用基于从头算分子动力学 (AIMD) 模拟（NVT 系综，200–700 K）的温度依赖有效势 (TDEP) 方法提取相互作用力常数 (IFCs)。
• 热输运建模： 通过两种互补的方法计算晶格热导率：
  1. 通过 ShengBTE 和 Fourphonon 软件包进行玻尔兹曼输运方程 (BTE) 计算，考虑了三声子和四声子散射。
  2. 使用维格纳输运方程 (WTE)，该方程统一了粒子性（对角项）和波动性（非对角项）对热导率的贡献。
• 成键与非谐性分析： 研究利用投影晶体轨道哈密顿密度 (pCOHP) 分析成键相互作用，通过计算均方位移 (MSD) 和原子轨迹来识别摇摆行为，并通过评估格吕奈森参数 (Grüneisen parameters) 和冻结声子势来量化非谐性。

主要贡献与结果

1. 确定 $\text{Sn(II)}$ 为摇摆原子：
   研究证实，以活性孤对电子为特征的 $\text{Sn(II)}$ 原子比 $\text{Sn(IV)}$ 和 S 原子表现出显著更大且各向异性的原子位移。势能垒分析表明，$\text{Sn(II)}$ 在沿 c 轴方向的位移具有最低能量势垒，表明其存在弱成键和由相邻孤对电子静电斥力驱动的“摇摆”行为。这导致了低频、平坦的光学声子支。

2. 各向异性热导率：
   由于其准一维结构，$\text{Sn}_2\text{S}_3$ 表现出显著的热各向异性。在室温下，沿 b 轴（成键方向，6.55 $\text{W m}^{-1} \text{K}^{-1}$）的晶格热导率显著高于沿 a 轴和 c 轴（范德华方向，约 1.6–1.7 $\text{W m}^{-1} \text{K}^{-1}$），其各向异性比为 3.86。

3. 光学声子的主导作用：
   与声子是主要热载体的传统观点相反，分析表明，沿 b 轴的光学声子对热导率的贡献约为 63%。这种异常归因于沿 b 轴的特定低频光学模式（5–32 meV）具有较高的群速度，并与声学支发生了耦合。

4. 波动性输运与温度依赖性：
   利用维格纳输运模型，研究发现总晶格热导率表现出微弱的温度依赖性（$T^{-0.68}$），偏离了粒子型声子气体模型所预测的标准 $T^{-1}$ 趋势。虽然在 300 K 时粒子性贡献占主导地位，但波动性贡献（源于声子支之间的隧道效应）也非常显著，特别是涉及光学声子的部分。

5. 非谐性机制：
   低热导率是由与 $\text{Sn(II)}$ 原子相关的强非谐性驱动的。这可以通过低频光学和声学支的大格吕奈森参数，以及冻结声子势与谐振拟合的显著偏差得到证实。

意义
本文声称，这些发现为深入理解混合价态化合物中的热输运提供了依据。通过将 $\text{Sn(II)}$ 特有的空间活性孤对电子与摇摆行为、低频光学模式以及强非谐性联系起来，该工作为发现具有本质低热导率的材料提供了路径。特别是，光学声子作为主要热载体的发现，挑战了在各向异性、弱结合系统中关于声子的传统假设，并为热电材料的设计提供了新见解。该研究并未提出具体的应用方案，而是强调了混合价态化合物作为热管理和热电优化领域一类材料的潜力。