FR38 · 连续互动课程
声子导热：从平均自由程到第二声、界面热阻与热电转换
在绝缘体和半导体里，热主要由声子携带，也就是晶格的量子化振动。傅里叶定律把声子的运动当成扩散，前提是所关心的每段距离上都有很多次碰撞。这门课讨论前提不成立时会怎样：比平均自由程还薄的膜，衰减太慢甚至翻负的热光栅，两种材料交界处的温度跳变，以及声子导热怎样限制热电发电。每个公式都算出数来，并用可以动手的模拟检验。
从一次声子飞行，到热波、界面跳变和热电
第一章用三个手册数值估算声子在硅里飞多远，第二章写出追踪这些飞行的玻尔兹曼方程。第三章精确求解两壁之间的薄膜，说明傅里叶定律在哪里给出过高的热流。第四章用两束激光画出热光栅，从衰减读出扩散率；第五章加入守恒动量的碰撞，衰减变成阻尼的热波，叫作第二声。第六章跟着声子穿过两种材料的边界，第七章给一条热电臂记能量账。最后一章读 Huberman 等人 2019 年关于石墨第二声的论文。
需要导数、指数函数和一点气体动理论。玻尔兹曼方程、热光栅的余弦条纹以及塞贝克效应和珀尔帖效应，在第一次用到时再介绍。 FR38
01 · L1
晶体里的热由谁携带，声子能飞多远
从硅的比热、声速和热导率反推声子平均自由程，约 38 nm；再用克努森数判断一段距离上的导热是扩散还是飞行。
02 · L2
声子玻尔兹曼方程：从飞行和碰撞推出傅里叶定律
在弛豫时间近似下写出声子玻尔兹曼方程，推出热导率公式和弹道极限，并算出 40 nm 硅膜里傅里叶定律为什么要得太多。
03 · L3
薄膜导热：在弹道与扩散之间
精确求解两块黑体壁之间的灰体薄膜，算出 40 nm 硅膜的热流只有傅里叶定律的 43%，并看懂壁面温度跳变和串联热阻近似。
04 · L2
瞬态热光栅：用干涉条纹测量热扩散
两束激光在样品表面画出冷热条纹，条纹衰减得多快就给出热扩散率；周期变短时衰减变慢，说明长自由程声子不再扩散。
05 · L3
正常散射与第二声：热为什么会像波一样来回振荡
区分守恒动量的正常散射和不守恒的电阻性散射，推出热的阻尼波方程，算出第二声速度、频率窗口和热光栅翻负的时刻。
06 · L3
界面热阻：热流在两种材料的交界处为什么会跳变
用声学失配和漫散射失配两个模型估算界面热导，算出硅与软膜之间约 1.4 µm 的卡皮查长度，并修正一个教科书公式。
07 · L2
热电转换的能量账：声子热导为什么决定效率
给一条热电臂记一笔能量账：导热、珀尔帖热和焦耳热各占多少，最大效率公式从哪里来，以及为什么降低声子热导能提高效率。
08 · L4
读 Huberman 等人 2019 年的石墨第二声论文
逐图读石墨第二声论文：300 K 的指数衰减、85 K 的信号翻负、温度窗口和第一性原理模拟，复算速度和极小值时刻，并列出读这类论文的清单。
沿途的实验
三个实验运行的都是模型，不是测量数据。薄膜实验在两块黑体壁之间求解灰体玻尔兹曼方程，并与独立的积分方程解和声子包蒙特卡罗核对。热光栅实验让带正常散射和电阻性散射的灰体声子分布随时间演化，并与另一条拉普拉斯域的解核对。界面实验把声子包射向边界，统计穿过的比例，与声学失配和漫散射失配公式比较。每个实验都把检验结果放在模型输出旁边，都不能代替真实样品的测量。
来源与复算范围
- Ioffe Institute, NSM Archive: Thermal properties of Silicon (Si) and Mechanical properties, elastic constants, lattice vibrations of Silicon (Si)
两页都通读过（力学页在同一目录的 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（横波）。课程把两个声速合成一个德拜平均声速，这一步是我们自己做的，不是页面给出的数值。
- 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
通读全文（arXiv 第 2 版）。引言部分：动理学理论给出的几十纳米平均自由程、把灰体近似称为过度简化、理论上常用的两块黑体壁之间的平板、由泵浦光交叉角决定的光栅周期，以及硅在 1 到 10 µm 周期下约 0.3 到 30 ns 的光栅衰减时间。摘要指出，硅的有效热导率在周期大到 10 µm 时仍有明显下降。课程没有引用这篇论文的公式编号。
- S. Lee, D. Broido, K. Esfarjani and G. Chen, Hydrodynamic phonon transport in suspended graphene, Nature Communications 6, 6290 (2015), doi:10.1038/ncomms7290
读过正文和方法部分，没有打开补充材料。这里用到：正常散射守恒晶体动量，倒逆散射等电阻性过程不守恒；第二声需要一个介于电阻性散射率和正常散射率之间的频率窗口；三维德拜固体中第二声速度约为声速除以根号三；以及预测的石墨烯第二声速度约 2000 到 3000 m/s。课程没有引用这篇论文的公式编号。
- 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
据 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 版逐行对照。
- Wikipedia, Interfacial thermal resistance, revision 1325816272 (5 December 2025)
通读这一修订版。这里用到：温度跳变的定义、声学失配和漫散射失配两个模型及其假设、细致平衡、“声学模型通常给出热阻上限、漫散射模型给出下限”的说法，以及大多数声子主导界面 80 到 300 MW/(m² K) 的热导范围。维基百科是二手来源，课程只用它的定义和典型范围。
- 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
部分阅读：摘要、引言，以及理论部分开头讲漫散射失配模型假设与细致平衡的小节。这里用到：该模型假定声子在漫散射后忘记来处；论文发现单个无序界面并不能抹掉这种记忆，因此模型的适用范围有限。原子尺度格林函数的结果和后面各节没有读，也没有引用。
- 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
在 PubMed Central 读过正文；推导所在的 SI 附录没有打开。这里用到：恒定物性模型给出的常用最大效率公式，即卡诺效率乘以一个取决于平均温度下优值的因子；吸热量等于导热加珀尔帖热再减去一半焦耳热；负载比及其最优值；以及常用公式把 SnSe 在热端 700 °C 时的效率高估两倍以上，原因是它的物性随温度变化很大。没有引用公式编号。
- Wikipedia, Thermoelectric materials, revision 1375991214 (21 September 2026)
部分阅读：导言，以及优值、器件效率和功率因子几节。这里用到：zT 等于电导率乘塞贝克系数平方再乘温度、除以热导率；最大效率公式以卡诺效率为第一个因子，只在物性不随温度变化时严格成立；商用热电制冷器的性能系数为 0.3 到 0.6。材料各节没有读。维基百科是二手来源。