管理热量如何在固体材料中传递是现代技术面临的一项基本挑战，从防止计算机芯片过热到将废热转化为电能，皆是如此。在大多数固体材料中，热量的传递主要通过构成晶体结构的原子振动来实现，而非通过电子的流动。科学家们长期以来一直寻求控制这种运动的方法，通常是通过改变材料的形状或添加杂质来散射这些振动，从而减缓热量的传递。然而，一个更困难的目标是寻找一种能够同时控制电荷运动和热量运动的方法，因为这两项特性往往以复杂的方式相互关联。如果研究人员能够设计出一种材料，使其能在像金属一样导电与像半导体一样阻断电流之间进行切换，同时提高其导热性能，那将为先进的热管理系统开启新的大门。
一支研究团队现在已经证明，对特定晶体进行极微小的化学成分改变，可以实现这一目标。他们从一种由钡、硼和碳组成的立方晶体 BaB3C3 开始。在其天然状态下，这种材料表现得像金属，允许电子自由流动。研究人员随后进行了一次精确的化学替代，在每个晶体重复单元中，仅用一个碳原子替换了一个硼原子。这一个简单的替换改变了结构中的总电子数，并迫使原子重新排列成一种新的、略有不同的形状。结果引发了剧烈的转变：该材料不再表现得像金属，而是变成了一种半导体——一种可以控制电流流动的物质。与此同时，这一变化使该材料的导热性能显著提升。
这一转变的关键在于原子之间如何进行键合。在原始的类金属晶体中，原子的连接方式留有一些空的强键合机会。当研究人员用碳替换硼时，多出的电子填补了这些空位，使得碳原子能够彼此形成新的、非常强且短的键。这种原子网络的重构就像是加紧了原子之间的弹簧。用物理术语来说，这种现象“硬化”了原子的振动，将最高振动频率推向更高处，并使整个结构更加刚硬。由于原子现在被结合得更加紧密，携带热量的振动可以移动得更快，且传播距离更长，而不会被散射或阻断。
通过使用强大的计算机模拟来模拟这些原子的行为，研究人员发现，新材料（现称为 BaB2C4）的导热能力大约是原始材料的两倍半。在室温下，原始立方晶体的导热速率为 7.6 瓦特每米每开尔文。经过替代后，新的四方晶体导热速率达到了 18.7 瓦特每米每开尔文。此外，热流变得具有方向性；它比垂直穿过晶体时更容易沿着晶体的平面方向移动，这种各向异性属性源于晶体结构新的、较低的对称性。研究人员通过将计算机模型与标准且高精度的计算进行对比验证了他们的发现，证实了他们所使用的机器学习势函数能够可靠地预测这些复杂的物理行为。
这项工作表明，定向化学替代是工程化材料的一种强大工具。通过简单地将一种原子类型改为另一种，是有可能重塑电子结构以创造半导体的，同时增强原子键以提升热导率。该研究并未声称已经解决了所有的热管理问题，但它为调节硼-碳笼状结构的这些特性提供了一种清晰且经证实的策略。研究人员提出，这种方法可以扩展到同家族的其他成分，从而可能揭示出一系列从金属向半导体过渡的材料谱系。此类材料对于未来需要精确控制电能和热能的技术至关重要，能将简单的化学微调转化为一种复杂的工程解决方案。
技术摘要：通过组分驱动实现 B→C 取代克拉特结构的金属-半导体转变及增强的声子传输
问题陈述
通过化学设计同时控制电子结构和热传输仍然是一个重大挑战，这对于热管理和能源材料至关重要。虽然外部刺激（应变、压力、场）和结构转变可以调节热导率，但化学取代提供了一种直接途径，用以改变电子计数并重构共价框架。然而，实现一种能够同时诱导金属-半导体转变并增强晶格热导率 (κ) 的策略并非易事，因为这些性质通常是相互竞争的。本研究旨在建立一种化学设计规则，使电子和热响应在单一固态材料内实现耦合，特别关注 B-C 克拉特系统。
方法论
本研究采用结合密度泛函理论 (DFT) 和机器学习原子间势 (MLIP) 的多尺度计算方法。
- DFT 计算： 使用带有投影缀加平面波 (PAW) 方法和 HSE06 杂化泛函的 Vienna Ab initio Simulation Package (VASP) 进行结构优化和电子结构分析，以确保精确的电子相关性。使用 LOBSTER 进行化学键合分析，计算投影晶体轨道哈密顿布居 (COHP) 和积分 COHP (iCOHP)。
- 热稳定性： 在 300 K 下进行了有限温度从头算分子动力学 (AIMD) 模拟，以验证结构完整性并生成热可及构型。
- 机器学习： 从头开始训练了一个消息传递原子簇展开 (MACE) 模型，该模型基于 DFT 生成的数据集（能量、力和坐标）。模型的准确性通过与 DFT 基准值的对比进行了验证，其能量和力的均方根误差 (RMSE) 均较低。
- 声子传输： 使用训练好的 MACE 势函数，通过 Phonopy 和 Phono3py 计算二阶和三阶原子间力常数 (IFC)。在 100–500 K 的温度范围内，通过求解完整的线性化玻尔兹曼输运方程 (LBTE) 来确定晶格热导率 (κ)。
核心贡献与结果
本研究调查了通过单个公式单位内的 B→C 取代，实现从立方相 BaB3C3 到四方相 BaB2C4 的转变。
电子结构转变：
- BaB3C3： 表现出金属特性，能带穿过费米能级 (EF)。
- BaB2C4： 取代过程捐赠了一个价电子，填充了此前未占据的成键态并使 EF 上移。这触发了金属-半导体转变，开启了一个 0.33 eV 的间接带隙。
- 键合重构： LOBSTER 分析揭示了共价网络显著的重新分布。虽然 B–B 和 B–Ba 相互作用适度增强，但最显著的变化是 BaB2C4 中强 C–C 键的形成（键长从 ~2.55 Å 缩减至 ~1.95 Å；iCOHP 量值从 ~-0.29 eV 增加至 ~-3.9 eV）。
晶格动力学与热传输：
- 声子硬化： BaB2C4 中强 C–C 键的形成使声子谱发生硬化，将最大频率从立方相的 ~800 cm−1 移动到四方相的 ~900 cm−1。
- 增强的导热性： 在 300 K 时，立方相 BaB3C3 的晶格热导率是各向同性的（κx=κy=κz），数值为 7.6 Wm−1K−1。相比之下，BaB2C4 表现出各向异性行为（κx=κy>κz），其面内 κ 显著增强至 18.7 Wm−1K−1，实现了约 2.5 倍的提升。
- 机制： 这种增强归因于两个因素：(1) 由于谱线硬化导致的声子群速度增加；(2) 由于在整个频率范围内非简谐散射事件减少而导致的声子寿命延长。
方法论验证：
- MACE-MLIP 成功高保真地重现了 BaB3C3 的 DFT 声子色散和热导率数值，验证了其在处理更复杂的 BaB2C4 系统时，能以极低的计算成本实现高精度模拟。
意义
本文证明，微小的化学取代 (B→C) 是在 B-C 克拉特系统中协同工程化电子和热学性质的一种可行策略。通过重构键合网络，本研究实现了双重目标：将材料从金属转变为半导体，同时增强了晶格热传输。这挑战了传统的权衡关系，即半导体通常由于不同的散射机制而表现出比金属更低的热导率；相反，研究表明特定的键合修改可以同时改善两者。这项工作强调了基于 MACE 的 MLIP 在准确捕捉复杂热传输现象中谐波和非简谐晶格动力学方面的效用，为研究结构和化学性质相似的化合物提供了一个可靠的框架。
每周获取最佳 applied physics 论文。
受到斯坦福、剑桥和法国科学院研究人员的信赖。
请查收邮箱确认订阅。
出了点问题，再试一次？
无垃圾邮件，随时退订。