超越Peierls-Boltzmann粒子图像：MOF-5与HKUST-1中声子寿命和带间相干耦合主导的热输运机制

《Small》：Unraveling the Mechanism for Phonon Based Heat Transport in Metal-Organic Frameworks

　　

热输运是金属-有机框架（metal-organic frameworks, MOFs）许多应用中的关键过程，包括气体吸附、催化和生物封装。设计具有最优热导率的MOFs的前提是深入理解其底层原子机制。这是该出版物的核心，研究人员首先表明，将晶态绝缘体中的热传导描述为声子（phonons；作为振动能量量子）的粒子状传播这一常用方法是不充分的。这是由于在此考虑的高度复杂材料（MOF-5和HKUST-1）中密集排列的声子带之间存在强带间耦合（inter-band couplings）。当通过适当的物理模型结合高精度机器学习势（machine-learned potentials）对此加以考虑时，与实验单晶数据实现定量一致成为可能。有趣的是，主导粒子状输运的声子主要是极低频率（高达～1.5 THz）的（声学）模式，而带间跃迁也发生在高得多的频率。这可归因于声子寿命（phonon lifetimes）随频率急剧下降，从而削弱粒子状输运但促进带间耦合。研究人员的分析不仅合理化了MOF-5与HKUST-1之间热导率的差异，而且为建立基于声子的MOFs热输运过程理解奠定了基础。

该论文发表于《Small》，围绕金属-有机框架（metal-organic frameworks, MOFs）中的热输运机制展开。MOFs由无机节点和有机连接体组装成高比表面积晶态多孔网络，广泛用于气体存储与分离、催化、传感及生物封装等。其吸附/脱附、催化周转和生物封装中的温度敏感性使热管理成为性能关键。然而，MOFs热导率（thermal conductivity, κ）的微观机制仍部分不明，单晶实验数据稀缺，MOF-5和HKUST-1的室温单晶κ分别约0.32和0.69 W m^?1 K^?1。过去模拟多基于分子动力学（molecular dynamics, MD）和Green-Kubo方法，能给出实空间瓶颈，但难以揭示特定振动本征模作用；基于声子（phonons）的非谐晶格动力学（anharmonic lattice dynamics, ALD）和Peierls-Boltzmann输运方程（BTE）在MOFs中应用极少。在复杂MOFs中，大晶胞导致密集振动谱、强非谐展宽和带重叠，准粒子图像可能不足，因此需要更一般的Wigner输运方程（Wigner transport equation, WTE）。研究人员以完全活化、无缺陷、三维周期的MOF-5和HKUST-1为原型，结合BTE/WTE和系统特异性机器学习矩张量势（moment tensor potentials, MTPs），旨在定量解释其热导率差异并建立基于声子的理解。研究得出，仅粒子状声子传播不足以描述热导率，必须包含带间相干（coherences）输运；WTE结合高精度MTPs可在室温与单晶实验数据定量一致；主导粒子状输运的是极低频率声学模式，带间跃迁延伸到更高频率，声子寿命随频率急降是关键。这为MOFs热输运的声子图像奠定基础，并提示热导率调控需同时考虑粒子与相干通道。

关键技术方法：研究人员采用色散校正密度泛函理论（density functional theory, DFT；VASP，PBE-D3(BJ)，PAW）获得参考能量、力与应力，并在主动学习框架下进行NpT/NVT采样；训练系统特异性MTPs（MLIP-2）；利用Phonopy/Phono3py有限位移计算二阶与三阶原子间力常数；在弛豫时间近似下求解BTE（BTE-RTA）和WTE，将晶格热导率分解为粒子状κ_P和相干κ_C；用DynaPhoPy从MD轨迹计算光谱能量密度（spectral energy density, SED），评估热膨胀和声子频率重整化；计算模式Grüneisen参数和散射相空间。实验单晶数据来自Huang et al.（MOF-5）和Babaei et al.（HKUST-1）等，用于验证。

2.1 基于声子的热输运描述：研究人员在非谐晶格动力学框架中，以声子频率、群速度、模式热容和线宽描述热输运；BTE-RTA给出粒子状κ_P，WTE总热导率为κ_P与κ_C之和，κ_C源于寿命展宽带间的Wigner速度耦合。结论是复杂MOFs需超越Peierls-Boltzmann图像。

2.2 MOF-5和HKUST-1的超低热导率：单晶实验显示二者κ极低；MTP计算表明HKUST-1更致密、更硬，但两者平均声速相近；简单声子气体模型无法解释低κ及二者差异。

2.3 光谱组成与振动状态：声子色散和态密度显示MOFs单位体积声学态远少于致密材料，光学声子密集；MOF-5在约0.6–0.8 THz出现谱拥挤，HKUST-1在约1 THz有平坦光学流形，提示寿命和带间作用重要。

2.4 BTE极限—population与coherence输运机制相互作用：BTE-RTA预测κ_P随温度近似T^?1下降，300 K低估实验约22%（MOF-5）和45%（HKUST-1）；引入WTE后κ_C随温度增加，κ_tot与实验差约10%，MOF-5为0.36 W m^?1 K^?1、HKUST-1为0.60 W m^?1 K^?1，且温度依赖减弱。缺陷或边界散射可改善MOF-5，但长度尺度未知。

2.5 热输运通道的光谱分析：κ_P主要来自1 THz以下声学声子；κ_C延伸到约10 THz且几乎由光学声子贡献；光学声子对κ_P约37%、对κ_C约90%，总贡献约55%。声子寿命随频率急降导致粒子通道窄、相干通道宽。

2.6 声子寿命趋势理解：模式Grüneisen参数显示MOF-5中约0.6 THz和1.6 THz以上平坦光学带强非谐，HKUST-1中约1.0 THz平坦带强非谐；散射相空间在HKUST-1低于0.4 THz较小，导致其低声学/光学声子寿命较大，κ_P较高。

2.7 光谱热导率的温度依赖：比较100和400 K，低温下声子寿命整体增大，κ_P升高、κ_C降低；κ_P频率依赖几乎不变，κ_C在低温下低频贡献增强、高频减弱，符合模式热容变化。

2.8 超越常规晶格动力学：热膨胀和声子频率重整化：HKUST-1负热膨胀很小，对声子谱和热导率影响可忽略；MOF-5负热膨胀约?16×10^?6 K^?1，远大于HKUST-1约?4×10^?6 K^?1；单独考虑300 K晶格参数会使κ_P坍缩、总κ降至0.23 W m^?1 K^?1，但SED表明非谐频率重整化与负热膨胀效应大部分补偿，因此基本晶格动力学近似仍可用于讨论主导机制。

讨论部分总结：该研究表明，基于声子可定量描述原型MOFs热导率，但必须包含WTE中的带间隧穿/跳跃相干输运；否则室温κ被低估。粒子状κ_P随温度下降，相干κ_C随温度上升，约450 K两者相当，导致总κ温度依赖弱。低频声子寿命急降使κ_P主要来自极低声学声子，而κ_C由更高频光学声子主导，平坦强非谐光学带贡献显著。HKUST-1室温κ约为MOF-5两倍，源于κ_P和κ_C均更大：前者与HKUST-1低频散射相空间较小、声子寿命较大有关，后者与其2–10 THz间125条声子带多于MOF-5的63条有关。结论翻译：总体而言，上述考虑表明，MOFs中的热输运是粒子状声子迁移与带间声子隧穿/跳跃之间的精细相互作用。对于前者，陡峭能带和高声子寿命有利；后者更受益于邻近能带的寿命展宽，这些能带理想情况下彼此大致平行。这意味着通常对一种输运有利的因素对另一种不利，因此不同声子类型通常主导两种输运通道。这使得先验地提炼出调控MOFs热导率的清晰设计准则极为困难。事实上，这需要更广泛的数据基础，需要对系统调控刚度、剪切模量等的多个材料系列进行详尽筛选。在此，当前研究提供了如何进行此类筛选以及需要考虑哪些方面的工具和指南。