FR38 · 热输运

界面热阻：边界上的温度跳变

把两种固体接在一起，让热流过去。为什么温度会在边界上跳一下？这个跳变什么时候占主导？ 声子走到两种固体的边界时，并不是都能过去，很多会被弹回，所以温度在接缝处突然下降。两种固体的软硬和密度差得越多、温度越低，这个跳变就越大。一层材料如果比卡皮查长度还薄（在这个模型里，失配很大的一对材料在室温下约为 1.4 微米），边界对热流的阻碍就比材料本身还大。

让热流过硅和另一层薄膜。声子在边界上被弹回，温度就在那里跳变；比较声学失配和漫散射失配两种模型。

界面热导 $90.9\,\mathrm{MW\,m^{-2}\,K^{-1}}$

热阻相同的硅层厚度（卡皮查长度） $1.41\times10^{3}\,\mathrm{nm}$

界面上的温降占比 $44.4\,\%$

温度计拟合的热导 / 模型值 $1.00 \pm 0.0089$

设置第二种固体与叠层

温度（K）

300

T

薄膜声速（km/s）

2

v_2

薄膜密度（g/cm³）

11

\rho_2

叠层厚度（nm）

1000

d

沿叠层拖动白色界面。

界面位置、噪声、声子包与种子

叠层中硅所占的比例

0.5

温度计噪声（K）

0.03

穿越统计的声子包数

1500

噪声与声子包的随机种子

3803

模型、单位与来源

$G=\frac{v_1}{4}\,\Gamma_{12}\,C_1(T;\,\omega<\omega_{D,2}),\qquad \Gamma_{12}=2\!\int_0^1\alpha_{12}(\mu)\,\mu\,d\mu,\qquad \alpha^{\mathrm{DMM}}_{12}=\frac{v_2^{-2}}{v_1^{-2}+v_2^{-2}}$

每一侧都是灰体德拜固体，三支声学支用同一个速度。高于薄膜截止频率的声子无法进入薄膜。声子包统计按同样的权重抽取入射角和频率，应在误差线内与模型一致。叠层两端的热库维持 10 K 温差；跳变用教科书里的热导，它比较的是射向界面的声子的温度，所以两侧相同时仍有辐射极限对应的跳变。

Wikipedia, Interfacial thermal resistance (revision of 2025-12-05) · Song and Chen (2021) · Ioffe NSM, Si mechanical properties

教程章节：界面热阻

模型假设

两个灰体德拜固体，每个有三支速度相同的声学支，在平整界面处相接。第 1 侧用硅的密度和薄膜实验里的灰体速度；第 2 侧是假设的固体，由速度和密度决定，原子数密度与硅相同。声学失配模型把声子当作平面波，界面完美、平整地结合；漫散射失配模型让界面忘掉声子从哪边来。两者都是教科书里的极限，不是对真实界面的预测。这里的热导比较的是射向界面的声子的温度；即使两侧完全相同，这个公式也会留下一个小跳变，教程中有解释。体热导率用 300 K 的动理学公式计算，温度变化时保持不变。温度计读数是合成数据：模型温度加带种子的高斯噪声。界面局部放大图是示意图。

叠层中的温度分布

黄色：模型。粉色：合成温度计读数（误差线为两倍噪声标准差）。绿色虚线：两侧各自的直线拟合。白线：界面。

界面热导随温度的变化

蓝色：声学失配。橙色：漫散射失配。白色虚线：辐射极限。黄色圆圈：当前模型值。粉色：温度计拟合（误差线为两倍标准误差）。

硅中声子穿过界面的概率

蓝色：声学失配。橙色：漫散射失配。绿色：所选模型的声子包统计（误差线为两倍标准误差）。

更多数值与检验方法

德拜温度：硅，薄膜: $0$
声阻抗：硅，薄膜: $0$
按热流加权的穿越概率: $0$
辐射极限：所有射来的声子都穿过: $0$
从薄膜一侧算出的同一热导（细致平衡）: $0$
体热导率：硅，薄膜: $0$
通过叠层的热流: $0$
温度计拟合给出的跳变及一倍标准误差: $0$
穿过的声子包比例及一倍标准误差: $0$
模型预测的声子包穿越比例: $0$
统计的声子包数: $0$
温度计噪声: $0$

细致平衡要求从两侧算出的热导相同，上面两行检验了这一点。声子包统计独立抽样入射角和频率，不使用计算模型热导的数值积分。

物理教程

热为什么在边界上绊一下

背景热从一种固体流进另一种固体时，温度在接缝处并不是平滑变化的，而是在边界上突然下降，好像那里藏着一层很薄的隔热片。这份额外的热阻叫界面热阻，也叫卡皮查热阻，得名于最早在液氦温度下测到它的物理学家。

为什么值得理解芯片把硅、氧化物和金属叠成几十纳米厚的层。一旦某层比它的卡皮查长度还薄，边界对热流的阻碍就比材料本身还大。

先抓重点

核心问题: 把两种固体接在一起，让热流过去。为什么温度会在边界上跳一下？这个跳变什么时候占主导？
一句话直觉: 只有一部分声子能过去。界面热导等于射来的声子的弹道热流乘上它们穿过的概率， $G=\tfrac14 v_1 C_1 \Gamma_{12}$ 。声学失配模型用声阻抗 $Z=\rho v$ 和斯涅尔定律算这个概率；漫散射失配模型只用细致平衡。

核心数学模型

声学失配的透射率

\alpha_{12}(\theta_1)=\frac{4Z_1Z_2\cos\theta_1\cos\theta_2}{\left(Z_1\cos\theta_2+Z_2\cos\theta_1\right)^2},\qquad \frac{\sin\theta_1}{v_1}=\frac{\sin\theta_2}{v_2}

超过临界角，即 $\sin\theta_1>v_1/v_2$ 时，没有声子能穿过。

卡皮查长度

\ell_K=\frac{\kappa_1}{G}

与界面热阻相同的硅层厚度。远比它薄的叠层，热阻主要来自界面。

常见难点

完美接合在这个公式里仍有跳变

容易误解两侧材料相同时，公式给出的热阻为零。

正确理解两侧相同时，所有声子都能穿过，但 $G$ 停在辐射极限 $\tfrac14 vC$ ，而不是无穷大。这个公式比较的是射向界面的声子的温度，不是温度计在各层内部读到的温度。教程推导了消掉这个剩余跳变的修正。

动手实验

01
看跳变
从“软而重的薄膜：明显的跳变”开始。
你会看到： 叠层颜色在白色界面处突变，局部放大图里大多数声子包以紫色被弹回，两条温度计拟合线在界面处相差好几开尔文。
02
找临界角
选择“硬而轻的薄膜：临界角”。
你会看到： 在透射率图中，蓝色曲线在约 32 度之后降到零，绿色的声子包统计点跟着它走。
03
降温
选择“低温界面：温降几乎都在界面”，再把温度滑块往回拖高。
你会看到： 在 10 K 时，几乎全部温降都集中在界面。低温下热导大致按温度的三次方下降，对数图上的圆点显示了这一点。