[toc]

在光滑的界面上“刻”一层无序：非晶表面如何让二维材料散热快三倍

导语

二维材料与三维衬底之间的界面热传导，一直是纳米电子器件散热的核心瓶颈。这篇发表在《International Journal of Heat and Mass Transfer》上的分子动力学模拟研究，发现了一个反直觉的规律：在晶体SiO₂衬底表面引入一层原子级薄的非晶层，可以让h-BN/SiO₂界面的热导提升约2.8倍——从12 MW m⁻² K⁻¹增至35 MW m⁻² K⁻¹。更关键的是，作者揭示了背后的声子机制：非晶表面的无序拓扑“唤醒”了原本沉默的中间频率和局域声子模式，为热量从二维材料流入三维衬底开辟了新的通道。

背景：2D/3D界面的“热堵车”

石墨烯、六方氮化硼（h-BN）等二维材料与SiO₂、GaN、SiC等三维衬底结合形成的异质结构，是新一代纳米电子器件的核心架构。但一个根本性问题始终困扰着这类器件：二维材料与三维衬底之间靠弱的范德华力（van der Waals force）结合，界面热导（interfacial thermal conductance, ITC）极低，导致热量在界面处“堵车”，形成巨大的温度跳跃。这会直接降低器件效率、缩短寿命、损害可靠性。

此前的研究尝试了多种策略来改善界面热传导：施加压力或应变、引入空位、修改表面化学性质、插入缓冲层等。其中，在衬底表面引入非晶层是一个有趣的方向——有研究报道，在化学键合的3D/3D界面（如Si/Ge）中，非晶层通常会降低界面热导，因为结构无序会强烈散射热载流子；但在2D/3D范德华界面中，情况似乎相反。然而，非晶层增强ITC的微观声子机制尚不清晰，尤其是非晶表面如何影响声子模式耦合、局域化行为以及频率分布，缺乏系统研究。

研究方法：给SiO₂表面“刻”上一层无序

作者采用非平衡分子动力学（NEMD）模拟，构建了h-BN单层支撑在SiO₂衬底上的异质结构模型。核心设计是：在晶体SiO₂（c-SiO₂）衬底顶部，通过“加热-淬火”方法生成一层约14 Å厚的非晶SiO₂（a-SiO₂）薄层，同时保持下方衬底仍为晶体结构，从而在不牺牲衬底本征高热导率的前提下，单独研究非晶表面对界面热传输的影响。

通过改变淬火条件，作者制备了不同rms粗糙度的非晶表面，系统考察了粗糙度、界面耦合强度、空位浓度及空位距界面距离对ITC的影响。

核心发现与逐图解读

图1：模型长什么样

图1（原文Fig. 1） 展示了模拟体系的示意图。子图a是h-BN/SiO₂范德华异质结构的整体构型，h-BN单层“半悬浮”在SiO₂衬底上，分为与衬底接触的重叠区和悬空区。子图b和c分别对应晶体（c-SiO₂）和非晶（a-SiO₂）表面。子图d–f展示了通过淬火退火技术生成的、具有不同rms粗糙度的非晶表面。

这张图的核心信息是：非晶层只存在于衬底表面约14 Å的薄层内，下方仍是晶体SiO₂。这种“三明治”式设计既利用了非晶表面对声子的调制作用，又保留了晶体衬底的高热导率，避免了全非晶衬底本身热导极低（0.1–1 W/mK）的弊端。

图2：非晶表面让ITC提升2.8倍

图2（原文Fig. 2） 是全文最核心的定量结果。子图a是c-SiO₂和a-SiO₂的径向分布函数（RDF），验证了非晶结构的合理性——c-SiO₂的Si-O、O-O、Si-Si第一峰分别位于1.62 Å、2.79 Å和3.12 Å，与文献和实验数据一致。

子图b给出了ITC和PDOS匹配因子S随表面粗糙度的变化。晶体表面（rms=0.12 Å）的ITC约为12 MW/m²K，而非晶表面在rms=2.1 Å时达到峰值约35 MW/m²K，提升约2.8倍。值得注意的是，这一数值甚至高于MoS₂/h-BN异质结构在硅衬底上的ITC（约32–33 MW/m²K）。S因子与ITC同步增长，说明PDOS重叠度的增加是ITC提升的关键驱动因素。

子图c展示了弛豫过程中h-BN/SiO₂异质结构的温度演化，验证了系统达到稳态。子图d对比了不同衬底厚度下的ITC——ITC几乎不随衬底厚度变化，而h-BN的本征热导率则表现出强烈的长度依赖性。这一对比清晰地表明：界面热导由界面处的声子透射决定，而非衬底内部的声子传输。

通俗地说：在光滑的晶体表面上“刻”出一层纳米级的无序，反而让热量更容易从二维材料流入三维衬底——但粗糙度存在一个最优值，超过rms≈2.1 Å后，过大的粗糙度会产生界面空洞，减少有效接触面积，ITC反而下降。

图3：非晶表面如何改变声子“频谱”

图3（原文Fig. 3） 从声子态密度（PDOS）角度揭示了ITC增强的微观机制。子图a对比了悬空h-BN、c-SiO₂支撑h-BN和a-SiO₂支撑h-BN的PDOS。最关键的差异出现在面外（ZA）声子分支：在悬空h-BN中低于3 THz的声子频率，在支撑结构中向高频移动。而在a-SiO₂支撑的h-BN中，约4 THz处出现了一个新峰，来源于面内振动模式向面外模式的转化。这种面内-面外模式的增强耦合对界面能量传递至关重要，因为低频声子才是界面热传输的主要贡献者。

子图b展示了c-SiO₂和a-SiO₂的PDOS分布。c-SiO₂在37.3 THz处有一个尖锐的声子谱峰，来源于晶格的平移对称性，而a-SiO₂由于缺乏长程周期性，该峰完全消失。子图c和d分别展示了c-SiO₂和a-SiO₂界面区域的层分辨PDOS。在c-SiO₂界面，前三个原子层在约2 THz处出现一个新的高度局域化声子模式，但迅速衰减；而在a-SiO₂中，5–25 THz中间频率区域的PDOS显著增加，高频和低频区域则有所下降，说明非晶层将高频声子散射并重新分布到了中间频率。

图4：热量主要走低频通道

图4（原文Fig. 4） 给出了频率分辨的ITC和累积ITC。子图a显示，无论是c-SiO₂还是a-SiO₂界面，约12 THz以下的低频声子贡献了约95%的总ITC，14 THz以上的声子贡献极为有限。这解释了为什么高频声子对ITC贡献小——它们不仅群速度低（图6a中高频分支平坦），而且更容易被散射（图6b中弛豫时间短）。

随着rms增加，约3 THz处的低频模式增强，9–13 THz范围出现蓝移且峰高增加。在2–12 THz范围内，非晶界面贡献了约92%的净谱热流增强，而12–23 THz范围仅贡献约8%。子图b和c对比了c-SiO₂和a-SiO₂支撑h-BN在重叠区域附近的层分辨面外PDOS。a-SiO₂支撑的h-BN在约1.5 THz处表现出更强的局域声子峰，且衰减更慢，这些局域声子与扩展声子耦合，将能量重新分配到面外模式，进而传入衬底。

图5：非晶表面让声子“既局域又离域”

图5（原文Fig. 5） 用声子参与率（PPR）刻画了声子的局域化程度。子图a对比了悬空h-BN、c-SiO₂支撑和a-SiO₂支撑h-BN的PPR。在低频范围，a-SiO₂支撑的h-BN出现PPR偏移，对应PDOS中的面外振动峰，说明局域声子与传播声子耦合，导致声子重新分布到不同频率模式。

子图b和c分别展示了c-SiO₂和a-SiO₂界面区域的层分辨PPR。两者在低频范围都出现蓝移，表明存在局域表面声子模式，但非晶表面的声子局域化显著更强，源于随机原子分布增强了声子散射。更有趣的是，a-SiO₂界面在低频范围表现局域化，而在5–25 THz中间频率区域则表现出更强的离域化——这些中间频率振动模式更容易通过声子散射穿过界面，从而促进界面声子传输。这与图3b中a-SiO₂的PDOS结果一致。

图6：声子色散与寿命

图6（原文Fig. 6） 展示了声子色散曲线和弛豫时间。子图a中，支撑h-BN在低频区中心附近出现显著蓝移，并出现了悬空h-BN中不存在的光学声子模式。这种蓝移源于两个因素：支撑h-BN的反射对称性破缺限制了面外振动模式的声子-声子散射相空间；弱的范德华相互作用改变了h-BN层内的原子间力常数。在Γ点附近，弯曲声子带变平（群速度趋近于零），这是面外方向平移对称性破缺的结果。

子图b显示，支撑体系中低频声学声子的弛豫时间显著降低，而a-SiO₂支撑的h-BN进一步降低，因为非晶表面的静态声子散射增加了额外的散射通道。声子弛豫时间的降低看似不利于热传导，但在界面处，更强的声子散射反而促进了声子模式的重新分布和转换，将高频声子转化为更容易穿过界面的低频声子——这与h-BN本征热导率从悬空的538 W/mK降至a-SiO₂支撑的330 W/mK形成对比：衬底支撑牺牲了h-BN的本征热导，却换来了界面热导的大幅提升。

图7：耦合强度激活更多声子

图7（原文Fig. 7） 考察了界面耦合强度χ（1、2、3）对ITC的影响。子图a和b显示，随着χ增加，c-SiO₂和a-SiO₂体系的ITC均单调增加，谱热导峰显著展宽并蓝移，扩展到13 THz。对于非晶表面，14–22 THz范围的高频声子也被激活，意味着更宽频率范围的声子参与界面热传输。子图c对比了NEMD和谱热流方法计算的ITC，两者趋势一致。

图8：空位效应取决于“离界面有多远”

图8（原文Fig. 8） 系统考察了c-SiO₂衬底中空位浓度和空位距界面距离对ITC的影响。结果显示，空位的增强效应高度非均匀：在距界面最近的Layer1中，1.5%空位浓度使ITC提升60%，但超过此浓度后ITC下降；在Layer2中，2%空位浓度使ITC提升84%；在更远的Layer3中，2%空位仅提升25%。PDOS分析表明，空位将35–38 THz的高频声子散射到26–30 THz和2–3 THz的低频模式，增强了h-BN面外声子与SiO₂的耦合。空位的影响随距界面距离增加而迅速衰减，说明界面附近的结构调控最为关键。

机理、意义与局限

这项工作的核心机理可以这样理解：在晶体衬底表面引入一层原子级无序，相当于在声子高速公路上增设了一个“换乘站”。原本，h-BN中携带热量的高频声子无法与SiO₂的声子模式匹配，就像两种不同轨距的列车无法直接换乘；而非晶表面的无序拓扑打破了这种“轨距壁垒”，通过声子散射将高频声子“降频”为低频声子，同时激活了中间频率的局域声子模式，让热量得以高效地从二维材料“换乘”到三维衬底。

作者直接计算的结果包括：ITC从12增至35 MW/m²K（NEMD），PDOS重叠因子S随粗糙度增加，空位浓度与ITC的非单调关系。作者对机制的解释——声子重新分布、面内-面外模式转换、局域-扩展声子耦合——则是基于PDOS、PPR、SED和谱热流分析的合理推断，但需要实验验证。

需要指出的是，这是一项纯分子动力学模拟研究，实验上直接测量非晶表面对ITC的影响极具挑战，因为很难将表面动力学与体材料动力学区分开。模拟中使用的L-J势函数和Tersoff势虽然经过验证，但实际界面可能存在氧化、污染、化学键合等复杂因素。此外，非晶层的热稳定性、长期可靠性以及在真实器件制造中的可控性，都是从模拟走向应用需要回答的问题。

写在后面

这篇论文真正想解决的问题是：2D/3D范德华异质结构中界面热导过低，限制了纳米电子器件的散热能力，而传统改善策略（压力、应变、化学修饰）各有局限，需要一种实验上可行、不牺牲衬底本征热导的界面工程方案。

它声称的贡献是：通过MD模拟证明，在晶体SiO₂表面引入原子级薄的非晶层，可将h-BN/SiO₂界面热导提升约2.8倍，并揭示了“非晶表面激活中间和局域声子模式、促进面内-面外声子耦合、将高频声子重新分布为低频热载流子”的微观机制。

最容易被reviewer攻击的地方在于：第一，纯模拟研究缺乏实验验证，L-J势函数对界面耦合强度的描述是否准确会直接影响ITC绝对值；第二，非晶层的生成方法（5000 K淬火）在实验上是否可行、能否精确控制rms粗糙度至2.1 Å，存在工程不确定性；第三，论文声称ITC提升2.8倍，但这一数值高度依赖于所选势函数和模型尺寸，与其他文献报道的绝对值存在差异（如实验值25.6 MW/m²K vs 模拟值12 MW/m²K），需要谨慎解读。

如果我是同方向的科研工作者，值得精读的两部分是：

第一部分：图3–图4的声子机制分析。 这是全文最具洞察力的部分。重点理解：非晶表面如何在PDOS中引入4 THz新峰（面内-面外模式转换）；c-SiO₂界面2 THz局域声子模式如何快速衰减，而a-SiO₂如何将高频声子重新分布到5–25 THz中间频率；频率分辨ITC如何证明12 THz以下声子贡献95%的热传输。这套“PDOS + 谱热流 + 累积ITC”的分析框架，可以直接迁移到其他2D/3D界面体系。

第二部分：图8关于空位距离效应的系统分析。 重点理解：空位增强ITC存在最优浓度（Layer1为1.5%，Layer2为2%），超过后因有效接触面积减少而下降；空位将高频声子散射为低频模式的具体频率通道（35–38 THz → 26–30 THz和2–3 THz）；空位影响随距界面距离增加而迅速衰减的定量规律（Layer1提升60%，Layer2提升84%，Layer3仅提升25%）。这部分为界面缺陷工程提供了可操作的设计准则。

可以先跳过的部分是：图2a的RDF验证和图6a的声子色散基础分析——前者是标准的结构合理性验证，后者关于蓝移和弯曲声子带变平的讨论在声子物理文献中已有充分覆盖，对于熟悉该领域的读者信息增量有限。

论文信息

标题：Surface amorphization enables phonon mode activation and enhanced thermal transport at 2D/3D van der Waals heterostructures

期刊：International Journal of Heat and Mass Transfer

年份：2026

DOI：10.1016/j.ijheatmasstransfer.2026.129663

本文为学术传播与知识分享目的撰写，不用于商业用途。如本文任何内容无意中侵犯了您的版权，请联系我们将及时核实并删除相关内容。