---
title: "声子导热：从平均自由程到第二声、界面热阻与热电转换"
description: "看懂傅里叶导热定律何时失效：声子自由程、薄膜、热光栅、第二声、界面跳变和热电效率。"
canonical_url: "https://sandboxphysics.com/zh/stories/phonon-heat-transport/"
language: "zh-CN"
page_type: "story"
alternate_languages:
  en: "https://sandboxphysics.com/stories/phonon-heat-transport/"
last_modified: "2026-10-31"
---

[声子导热](https://sandboxphysics.com/zh/stories/phonon-heat-transport/)

FR38 · 连续互动课程

# 声子导热：从平均自由程到第二声、界面热阻与热电转换

Sandbox Physics · 2026-10-30

在绝缘体和半导体里，热主要由声子携带，也就是晶格的量子化振动。傅里叶定律把声子的运动当成扩散，前提是所关心的每段距离上都有很多次碰撞。这门课讨论前提不成立时会怎样：比平均自由程还薄的膜，衰减太慢甚至翻负的热光栅，两种材料交界处的温度跳变，以及声子导热怎样限制热电发电。每个公式都算出数来，并用可以动手的模拟检验。

[从平均自由程开始](https://sandboxphysics.com/zh/stories/phonon-heat-transport/heat-carriers-and-mean-free-path/) 先动手试试

声子在两次碰撞之间大约飞行一个平均自由程 $\Lambda$。膜比 $\Lambda$ 厚得多时（左），热靠许多短步穿过，傅里叶定律成立。克努森数 $\mathrm{Kn}=\Lambda/L$ 接近 1 时（右），许多声子不经散射，从热壁直接飞到冷壁。原创示意图，未按比例。

## 从一次声子飞行，到热波、界面跳变和热电

第一章用三个手册数值估算声子在硅里飞多远，第二章写出追踪这些飞行的玻尔兹曼方程。第三章精确求解两壁之间的薄膜，说明傅里叶定律在哪里给出过高的热流。第四章用两束激光画出热光栅，从衰减读出扩散率；第五章加入守恒动量的碰撞，衰减变成阻尼的热波，叫作第二声。第六章跟着声子穿过两种材料的边界，第七章给一条热电臂记能量账。最后一章读 Huberman 等人 2019 年关于石墨第二声的论文。

需要导数、指数函数和一点气体动理论。玻尔兹曼方程、热光栅的余弦条纹以及塞贝克效应和珀尔帖效应，在第一次用到时再介绍。 [FR38](https://github.com/yangzichao/sandboxphysics/issues/68)

1. 01 · L1

   ## [晶体里的热由谁携带，声子能飞多远](https://sandboxphysics.com/zh/stories/phonon-heat-transport/heat-carriers-and-mean-free-path/)

   从硅的比热、声速和热导率反推声子平均自由程，约 38 nm；再用克努森数判断一段距离上的导热是扩散还是飞行。

2. 02 · L2

   ## [声子玻尔兹曼方程：从飞行和碰撞推出傅里叶定律](https://sandboxphysics.com/zh/stories/phonon-heat-transport/phonon-boltzmann-equation/)

   在弛豫时间近似下写出声子玻尔兹曼方程，推出热导率公式和弹道极限，并算出 40 nm 硅膜里傅里叶定律为什么要得太多。

3. 03 · L3

   ## [薄膜导热：在弹道与扩散之间](https://sandboxphysics.com/zh/stories/phonon-heat-transport/ballistic-and-diffusive-films/)

   精确求解两块黑体壁之间的灰体薄膜，算出 40 nm 硅膜的热流只有傅里叶定律的 43%，并看懂壁面温度跳变和串联热阻近似。

4. 04 · L2

   ## [瞬态热光栅：用干涉条纹测量热扩散](https://sandboxphysics.com/zh/stories/phonon-heat-transport/transient-thermal-grating/)

   两束激光在样品表面画出冷热条纹，条纹衰减得多快就给出热扩散率；周期变短时衰减变慢，说明长自由程声子不再扩散。

5. 05 · L3

   ## [正常散射与第二声：热为什么会像波一样来回振荡](https://sandboxphysics.com/zh/stories/phonon-heat-transport/normal-scattering-and-second-sound/)

   区分守恒动量的正常散射和不守恒的电阻性散射，推出热的阻尼波方程，算出第二声速度、频率窗口和热光栅翻负的时刻。

6. 06 · L3

   ## [界面热阻：热流在两种材料的交界处为什么会跳变](https://sandboxphysics.com/zh/stories/phonon-heat-transport/interface-thermal-resistance/)

   用声学失配和漫散射失配两个模型估算界面热导，算出硅与软膜之间约 1.4 µm 的卡皮查长度，并修正一个教科书公式。

7. 07 · L2

   ## [热电转换的能量账：声子热导为什么决定效率](https://sandboxphysics.com/zh/stories/phonon-heat-transport/thermoelectric-energy-ledger/)

   给一条热电臂记一笔能量账：导热、珀尔帖热和焦耳热各占多少，最大效率公式从哪里来，以及为什么降低声子热导能提高效率。

8. 08 · L4

   ## [读 Huberman 等人 2019 年的石墨第二声论文](https://sandboxphysics.com/zh/stories/phonon-heat-transport/reading-huberman-2019/)

   逐图读石墨第二声论文：300 K 的指数衰减、85 K 的信号翻负、温度窗口和第一性原理模拟，复算速度和极小值时刻，并列出读这类论文的清单。

## 沿途的实验

[01 **让飞行与扩散赛跑** 加热类硅薄膜的一侧，看声子包穿过它。把膜越做越薄，直到傅里叶定律给出的热流过高，弹道极限接管。](https://sandboxphysics.com/zh/labs/phonon-ballistic-diffusive-film/) [02 **让热振荡起来** 用两束泵浦脉冲画出冷热条纹，看它们淡去、振铃或翻转符号；再在留出的周期上检验扩散和阻尼波两个模型。](https://sandboxphysics.com/zh/labs/transient-thermal-grating/) [03 **找出边界上的跳变** 让热流过硅和另一层薄膜，看声子在边界上穿过或弹回，再从合成温度计读数里读出跳变。](https://sandboxphysics.com/zh/labs/interface-thermal-resistance/)

三个实验运行的都是模型，不是测量数据。薄膜实验在两块黑体壁之间求解灰体玻尔兹曼方程，并与独立的积分方程解和声子包蒙特卡罗核对。热光栅实验让带正常散射和电阻性散射的灰体声子分布随时间演化，并与另一条拉普拉斯域的解核对。界面实验把声子包射向边界，统计穿过的比例，与声学失配和漫散射失配公式比较。每个实验都把检验结果放在模型输出旁边，都不能代替真实样品的测量。

## 来源与复算范围

1. [Ioffe Institute, NSM Archive: Thermal properties of Silicon (Si) and Mechanical properties, elastic constants, lattice vibrations of Silicon (Si)](http://www.ioffe.ru/SVA/NSM/Semicond/Si/thermal.html)

   两页都通读过（力学页在同一目录的 mechanic.html）。用到的 300 K 数值：比热 0.7 J/(g K)、热导率 1.3 W/(cm K)、热扩散率 0.8 cm²/s、密度 2.329 g/cm³，以及沿 [100] 方向的声速 8.43 km/s（纵波）和 5.84 km/s（横波）。课程把两个声速合成一个德拜平均声速，这一步是我们自己做的，不是页面给出的数值。

2. [A. A. Maznev, J. A. Johnson and K. A. Nelson, Onset of nondiffusive phonon transport in transient thermal grating decay, Physical Review B 84, 195206 (2011), arXiv:1108.3770v2](https://arxiv.org/abs/1108.3770)

   通读全文（arXiv 第 2 版）。引言部分：动理学理论给出的几十纳米平均自由程、把灰体近似称为过度简化、理论上常用的两块黑体壁之间的平板、由泵浦光交叉角决定的光栅周期，以及硅在 1 到 10 µm 周期下约 0.3 到 30 ns 的光栅衰减时间。摘要指出，硅的有效热导率在周期大到 10 µm 时仍有明显下降。课程没有引用这篇论文的公式编号。

3. [S. Lee, D. Broido, K. Esfarjani and G. Chen, Hydrodynamic phonon transport in suspended graphene, Nature Communications 6, 6290 (2015), doi:10.1038/ncomms7290](https://doi.org/10.1038/ncomms7290)

   读过正文和方法部分，没有打开补充材料。这里用到：正常散射守恒晶体动量，倒逆散射等电阻性过程不守恒；第二声需要一个介于电阻性散射率和正常散射率之间的频率窗口；三维德拜固体中第二声速度约为声速除以根号三；以及预测的石墨烯第二声速度约 2000 到 3000 m/s。课程没有引用这篇论文的公式编号。

4. [S. Huberman, R. A. Duncan, K. Chen, B. Song, V. Chiloyan, Z. Ding, A. A. Maznev, G. Chen and K. A. Nelson, Observation of second sound in graphite at temperatures above 100 K, Science (2019), doi:10.1126/science.aav3548, arXiv:1901.09160v1](https://arxiv.org/abs/1901.09160)

   据 arXiv 第 1 版通读全文，包括补充材料。这里用到：515 nm 泵浦、532 nm 探测光外差检测的瞬态光栅装置；扩散时间常数等于周期平方除以四倍 π 平方再除以扩散率；300 K 下 37.5 µm 周期测得 11 cm²/s，以及短周期时比扩散更慢的衰减（图 1）；85 K 时信号翻负和由第一个极小值得到的约 3200 m/s 速度（图 2）；104 K 和 125 K 仍有振荡、150 K 和 50 K 时消失（图 3）；预测的 50 到 250 K 窗口（图 4）；第一性原理给出的 3650 m/s；以及补充材料中 7.5 µm 曲线在约 1.2 ns 处的极小值和振荡期间约 0.2 K 的光栅振幅。没有把《科学》正式发表版与 arXiv 第 1 版逐行对照。

5. [Wikipedia, Interfacial thermal resistance, revision 1325816272 (5 December 2025)](https://en.wikipedia.org/w/index.php?title=Interfacial_thermal_resistance&oldid=1325816272)

   通读这一修订版。这里用到：温度跳变的定义、声学失配和漫散射失配两个模型及其假设、细致平衡、“声学模型通常给出热阻上限、漫散射模型给出下限”的说法，以及大多数声子主导界面 80 到 300 MW/(m² K) 的热导范围。维基百科是二手来源，课程只用它的定义和典型范围。

6. [Q. Song and G. Chen, Evaluation of diffuse mismatch model for phonon scattering at disordered interfaces, Physical Review B 104, 085310 (2021), arXiv:2106.04745v2](https://arxiv.org/abs/2106.04745)

   部分阅读：摘要、引言，以及理论部分开头讲漫散射失配模型假设与细致平衡的小节。这里用到：该模型假定声子在漫散射后忘记来处；论文发现单个无序界面并不能抹掉这种记忆，因此模型的适用范围有限。原子尺度格林函数的结果和后面各节没有读，也没有引用。

7. [H. S. Kim, W. Liu, G. Chen, C.-W. Chu and Z. Ren, Relationship between thermoelectric figure of merit and energy conversion efficiency, PNAS 112(27), 8205 (2015), doi:10.1073/pnas.1510231112](https://pmc.ncbi.nlm.nih.gov/articles/PMC4500231/)

   在 PubMed Central 读过正文；推导所在的 SI 附录没有打开。这里用到：恒定物性模型给出的常用最大效率公式，即卡诺效率乘以一个取决于平均温度下优值的因子；吸热量等于导热加珀尔帖热再减去一半焦耳热；负载比及其最优值；以及常用公式把 SnSe 在热端 700 °C 时的效率高估两倍以上，原因是它的物性随温度变化很大。没有引用公式编号。

8. [Wikipedia, Thermoelectric materials, revision 1375991214 (21 September 2026)](https://en.wikipedia.org/w/index.php?title=Thermoelectric_materials&oldid=1375991214)

   部分阅读：导言，以及优值、器件效率和功率因子几节。这里用到：zT 等于电导率乘塞贝克系数平方再乘温度、除以热导率；最大效率公式以卡诺效率为第一个因子，只在物性不随温度变化时严格成立；商用热电制冷器的性能系数为 0.3 到 0.6。材料各节没有读。维基百科是二手来源。