极化激元纳米线中的辐射热传导

论文信息： L. C. McCormack and M. Francoeur, Radiation conduction in polaritonic nanowires, arXiv:2605.24030 [physics.comp-ph] (2026).

论文链接：https://doi.org/10.48550/arXiv.2605.24030

研究背景

在微纳电子器件中，结构尺寸不断缩小会增强边界散射，从而显著降低传统声子热导率。近年来，表面声子极化激元被认为可能提供一种额外的热输运通道：它们由电磁波与极性介质中的光学声子耦合形成，可以沿材料界面传播，并在纳米尺度热管理中发挥作用。已有实验报道表明，在某些极性纳米结构中，极化激元可能对热导率产生可观贡献。然而，目前对这些实验结果的解释仍存在争议，核心原因是以往大量理论分析依赖动理学理论和色散关系，而缺少基于涨落电动力学的全波模型验证。本文正是针对这一问题，重新评估 SiO₂ 极化激元纳米线中的辐射导热能力。

研究内容

本文研究的是 SiO₂ 纳米线中的辐射热导率，也可以理解为由热电磁涨落和声子极化激元共同介导的轴向热输运。研究对象是直径为 66 nm 和 132 nm 的 SiO₂ 纳米线，工作温度接近室温至 400 K。为了从严格电磁理论出发计算热输运，本文采用基于涨落电动力学的离散系统格林函数方法，即 DSGF 方法。该方法把纳米线离散成许多小体元，通过计算体元之间的系统格林函数和谱透射系数，得到沿纳米线方向的辐射热导率。

图 1：沿 SiO₂ 纳米线在扩散输运区中的热电磁输运。该纳米线的直径为 D ，长度为 L 。示意图中，纳米线直径为 66 nm，当长度达到 200 nm 时进入扩散区。纳米线被离散为 32,184 个均匀立方子体元。其中一半纳米线保持在温度 T + δ T ，另一半保持在温度 T 。

计算模型中，纳米线被人为分成两个等体积区域，一半温度为 (T+δ T)，另一半温度为 (T)，在 δT 趋于 0 的极限下提取扩散区间的辐射热导率。对于 66 nm 和 132 nm 直径的 SiO₂ 纳米线，长度取 200 nm 时已经达到扩散极限，继续增加长度对结果影响很小。DSGF 结果显示，在 400 K 时，66 nm 纳米线的辐射热导率为 0.0264 W m⁻¹ K⁻¹，132 nm 纳米线为 0.0165 W m⁻¹ K⁻¹。也就是说，辐射导热存在，但并不强。

图 2：采用 DSGF 方法、传统动力学理论（KT）以及引入有效平均自由程的动力学理论计算得到的扩散区中 SiO₂ 纳米线轴向总辐射热导率随温度的变化关系：（A） D = 66 nm ，（B） D = 132 nm 。采用 DSGF 方法预测的 SPhPs 和 BPhPs 对总辐射热导率的贡献随温度的变化关系：（C） D = 66 nm ，（D） D = 132 nm 。

随后，本文区分了两类极化激元贡献：表面声子极化激元和体声子极化激元。SiO₂ 在两个 Reststrahlen 波段内介电函数实部为负，支持沿表面传播的 SPhPs；在这些波段之外，介电函数实部为正，热电磁输运则可由纳米线体内传播的 BPhPs 介导。DSGF 计算表明，SiO₂ 纳米线中的辐射热导率主要由 SPhPs 贡献，而不是由体内传播的 BPhPs 主导。这一点对于理解纳米线中“极化激元导热”的物理本质很关键，因为它说明界面模式仍然是主要通道。

图 3：无限长度 SiO₂ 纳米线在直径为 66 nm 和 132 nm 时的色散关系：（A）角频率 ω 随轴向波矢实部 R e ( k z ) 的变化关系；（B）角频率 ω 随轴向波矢虚部 I m ( k z ) 的变化关系。在（A）和（B）中均绘出了 SiO₂ 中的光线。（C）当 D = 66 nm 和 132 nm 时的谱平均自由程。在所有图中，Reststrahlen 带均以灰色阴影区域标出。

本文进一步将 DSGF 结果与传统动理学理论进行对比。传统方法先由 Maxwell 方程推导无限长圆柱纳米线的色散关系，再由复轴向波矢得到辐射平均自由程，并代入类似 Boltzmann 输运方程的表达式。结果发现，传统动力学理论显著高估了辐射热导率。例如，对于 66 nm SiO₂ 纳米线，在 400 K 时，动理学理论预测值约为 1.18 W m⁻¹ K⁻¹，约为 DSGF 结果的 45 倍，甚至接近 SiO₂ 的体声子热导率。这种高估主要来自 BPhPs 的长平均自由程，但该理论没有充分考虑纳米线有限横向尺寸导致的边界散射效应。为改善这一问题，本文在动理学理论中引入了有效平均自由程，用 Matthiessen 规则修正 Reststrahlen 波段外 BPhPs 的尺寸效应。修正后，66 nm 纳米线在 400 K 的辐射热导率变为 0.0227 W m⁻¹ K⁻¹，与 DSGF 结果只相差约 14%；132 nm 纳米线的结果也有所改善，但仍低于 DSGF 结果的一半。更重要的是，即使总热导率可以通过有效平均自由程得到较好修正，谱分布仍然无法准确再现。换言之，动理学理论经过修正后可以用于估算总量，但不能可靠描述不同频率、不同极化激元模式对辐射导热的具体贡献。

图 4：400 K 下，采用 DSGF 方法、传统动力学理论（KT）以及引入有效平均自由程的动力学理论计算得到的扩散区中 SiO₂ 纳米线轴向谱辐射热导率：（A） D = 66 nm ，（B） D = 132 nm 。在两幅图中，Reststrahlen 带均以灰色阴影区域标出。

结论与展望

这篇论文的关键价值在于，它用全波涨落电动力学方法重新审视了极化激元纳米线中的辐射导热问题。过去基于动理学理论的研究往往预测声子极化激元能够显著补偿纳米结构中声子热导率的下降，甚至给出接近常规声子热导率量级的辐射热导率。然而，本文的 DSGF 计算表明，对于 66 nm 和 132 nm 的 SiO₂ 纳米线，400 K 下的辐射热导率远低于传统动理学理论预测。造成差异的核心在于，传统理论高估了体声子极化激元在纳米线体内的传播平均自由程，而全波模型自然包含了有限尺寸、边界和相干电磁输运效应。引入有效平均自由程后，动理学理论对总热导率的预测可以明显改善，但仍无法准确给出谱分布。这说明，在设计基于声子极化激元的微纳热管理器件时，仅依赖色散关系和 Boltzmann 图像是不够的，必须结合涨落电动力学全波模型来判断极化激元到底能贡献多少实际热输运。

下面是会议通知：