热导率、电导率、热扩散系数与比热：物理内涵、测量关系与材料实证

说明：本文华算科技主要介绍热导率、电导率、热扩散系数和比热分别对应的物理过程、常用单位与测量关系，并以岩石、玻璃和半导体材料的实验结果解释四个量为何不能互相替代。

热导率、电导率、热扩散系数和比热

一块材料两端存在温差时，热导率κ描述给定温度梯度下形成多大的传导热流密度。傅里叶定律写作 q＝−κ∇T，负号表示热流指向低温处；κ 的单位是 W·m⁻¹·K⁻¹。同样的温差、厚度与截面积下，κ 较大的均匀试样通常传出更多热量。

受热后的温度分布还会随时间改变。热扩散系数 α表征温度不均匀消退的快慢，单位为 m²·s⁻¹。榴辉岩实验把 α 与 κ 分别列为两组曲线：升温时两者都降低，但下降幅度和单位不同，不能把任意一条曲线改名为另一个量。

图1. 榴辉岩的热扩散系数与热导率随温度、压力变化的实验曲线。DOI：10.1038/s41467-022-32484-w

热导率还容易与整块试样的导热量混淆。长度为 L、截面积为 A 的均匀棒在一维稳态下，热导 G_th＝κA/L，单位为 W·K⁻¹。κ 是材料沿指定方向的参数，G_th则随棒的粗细和长度改变；比较两根棒传热多少，不能只看材料名。

温度为什么会升高，取决于输入热量相对于材料储热量有多大。比热容 c_p通常指单位质量材料在定压条件下升高 1 K 所需的热量，单位为 J·kg⁻¹·K⁻¹；论文中的“比热”常是这个质量比热容。材料总质量增加，所需热量增加，但材料的 c_p不会因试样切大而翻倍。

若改问“这整块试样升高 1 K 需要多少热量”，量名是热容 C＝mc_p，单位为 J·K⁻¹。把 1 g 样品换成 2 g，同一温区的质量比热不因质量翻倍而改变，但整块试样的热容翻倍；论文表格中同样写着“heat capacity”时要看其单位是否已经按质量归一。

给材料施加电场时，电导率σ则描述电场产生电流密度的能力。在线性响应区，J＝σE，σ 的 SI 单位为 S·m⁻¹。传热涉及能量，导电涉及电荷；金属电子可以同时参与两者，绝缘体却仍能靠晶格振动导热。

四个量因此回答不同问题：κ 与 σ 分别对应热流和电流，c_p对应升温所需的能量，α 对应温度分布演化的时间尺度。若长度尺度为 L，在均匀材料的扩散近似下，温差消退的特征时间与L²/α同数量级；换一块厚度加倍的试样，时间尺度约增至四倍。

比热和密度如何换算热扩散系数？

给 1 kg 材料输入相同热量，c_p较大的材料温升较小。在玻璃态 SiO₂中，低温时部分振动模式尚未充分参与储热，随着温度升高，比热逐渐增加；原论文同时给出了计算的定容比热曲线和实验的定压比热点。两种热容口径在这组数据的温区内接近，但计算曲线仍不能当成对 c_p的直接测量。

图2. 玻璃态 SiO₂的计算定容比热与实验定压比热随温度变化。DOI：10.1038/s41524-023-01033-4

定压比热与定容比热分别对应外界压力固定、体积固定的升温过程。固体在许多常见温区两者数值接近，接近不代表定义相同；温度升高带来的膨胀、相变附近的响应，都会让两种实验条件的差异变得重要。玻璃态 SiO₂的图上，空心圆是定压实验点，实线主要是定容计算值。

单位质量的 c_p还要乘以密度 ρ，才能得到每单位体积升高 1 K 需要的热量，即体积热容 ρc_p，单位为 J·m⁻³·K⁻¹。相同的 c_p若对应不同孔隙率和体积密度，单位体积材料吸收同样热量后的温升也可能不同。

多孔体里的 ρ 应与热流经过的试样体积对应。把致密晶粒的真密度代入多孔块体的有效热扩散系数，会把孔洞占据的体积当成实体，从而高估换算出的热导率。体积口径一致比保留更多有效数字更关键，烧结体与粉体压片尤其如此。

在局部热平衡、均匀各向同性材料的常用模型中，热传导方程给出α＝κ/(ρc_p)。分子κ 使热量沿温度梯度传开，分母 ρc_p决定局部升温需要多少热量。若两块试样 κ 相同、ρ 相同，而其中一块 c_p加倍，它的 α 就减半。

这项关系也解释了“导热强”与“升温快”为何不是同一判断。榴辉岩中绿辉石比例与含水量变化时，实验分别记录 α 和 κ；两者随含水量下降的趋势相近，但比较温度传播的快慢还须计入样品的体积热容。

从局部能量守恒看，热流进入一个小体积后，一部分能量用于提高该处温度；ρc_p越大，相同净入热流造成的温升越小。若试样内部还有相变潜热、化学反应或强烈温度依赖，单个常数 α 已不足以描述整段升温过程，需要把随温度变化的 κ 和 c_p代入热传导方程。

图3. 榴辉岩中绿辉石比例与含水量对热扩散系数、热导率的影响。DOI：10.1038/s41467-022-32484-w

实际计算常采用κ＝αρc_p：热扩散系数来自瞬态温度响应，c_p来自量热测量，ρ 来自密度测量。若 α＝2 mm²·s⁻¹、ρ＝5000 kg·m⁻³、c_p＝500 J·kg⁻¹·K⁻¹，先将 α 换成 2×10⁻⁶ m²·s⁻¹，得到 κ＝5 W·m⁻¹·K⁻¹。

传热与储热如何决定热扩散系数？

单一电子型载流子且处于低场线性区时，σ≈neμ，n 为参与传输的电子浓度，e 为元电荷，μ 为迁移率。增加 n 可提高电导率，但若杂质散射同时降低 μ，σ 的变化要由两个因子共同决定；多载流子材料还需要把各通道的电导贡献相加。

该式也解释了为何两块样品的电导率可以一样：一块具有较高 n 与较低 μ，另一块则相反。σ是宏观电流响应，单凭它无法反推出载流子浓度或迁移率。半导体掺杂比较中，霍尔测量常用来补充 n 的信息，但多带传输时单一霍尔斜率还会受到不同通道权重影响。

Bi 掺杂SnSe 单晶的霍尔响应反映了载流子类型与浓度的变化。原论文比较未掺杂样品和几种 Bi 含量，并追踪 n 随温度升高的变化；这些数据是讨论 σ 的输入信息，霍尔斜率本身并非热导率的读数。

图4. 未掺杂与 Bi 掺杂 SnSe 单晶的霍尔响应和载流子浓度。DOI：10.1038/ncomms13713

未掺杂SnSe 与 Bi 掺杂 SnSe 的霍尔曲线斜率方向相反，分别对应空穴型和电子型的主要载流子。载流子符号改变不会让热导率自动反号：电导率和热导率都是非负的纵向输运系数，正负号主要出现在霍尔系数或 Seebeck 系数这类响应中。

金属的传热可包含电子贡献 κ_e，晶格振动则提供声子贡献 κ_ph。在适用的金属近似中，Wiedemann–Franz 关系写作 κ_e≈LσT，其中 L 为洛伦兹数；它联系的是电子热导率与电导率，不能直接把总热导率 κ 全部换成 σT。

常见的金属参考值 L₀≈2.44×10⁻⁸ W·Ω·K⁻²，来自简并自由电子近似。实际洛伦兹数会受能带和散射机制影响，L 不必等于 L₀；半导体的载流子统计与声子贡献不同，直接用 L₀σT 估计总 κ 尤其容易失准。

SnSe 的结果给出另一种情况：Bi 掺杂显著改变沿 c 轴的 σ，几组样品的总热导率 κ 接近。论文指出该半导体的热流主要由晶格振动承担，电子热流占比较小；所以电流增强并不要求 κ 按相同比例增加。

图5. 不同 Bi 掺杂 SnSe 单晶的电导率、热导率及热电参数随温度变化。DOI：10.1038/ncomms13713

即使材料是金属，温度、散射类型也会改变电子携热与携电之间的比例；含有显著声子热流时，κ_e与总 κ 更需分开。比较两种材料的 σ 和 κ 时，应交代温度、载流子浓度和声子贡献，“导电好就一定导热好”并非跨材料的普遍定律。

测量方向与方式如何影响数值读数？

层状晶体的原子连接与能带色散沿晶轴不同，热流和电流也会表现出方向依赖。Bi 掺杂 SnSe 单晶沿 a、b、c 轴测得的 σ 与 κ 均不相同；其中沿 a 轴的电导率和热导率低于另两个方向。写一个不标方向的单晶数值，容易把不同晶轴的输运混在一起。

图6. Bi掺杂 SnSe 单晶沿不同晶轴测得的电导率、热导率及热电参数。DOI：10.1038/ncomms13713

同一SnSe 样品沿 a 轴与 b 轴的数据差别，来自晶体内不同方向的成键和输运通路；把它磨成多晶块体后，晶粒取向分布还会改变有效数值。单晶轴向值与多晶平均值必须区分，织构明显的压片尤其不能简单视为各向同性。

各向异性材料的 κ 和 σ 严格说可用张量表示，日常实验常报告某个主轴或试样方向的分量。ρ 与 c_p是标量，但由它们换算 α 时，κ 与 α 应取同一方向；SnSe 晶片若垂直 a 轴受热，所得扩散系数不能与沿 b 轴的热导率直接相乘。

激光闪射法让脉冲加热片状试样的一面，记录另一面的温度随时间升高，再由厚度与时间响应拟合α。差示扫描量热可测 c_p，密度测量给出 ρ，三者合起来计算 κ；因此一张表中的 κ 可能是组合结果，而非仪器直接读出的原始信号。

薄膜还可以用瞬态热光栅测面内热扩散：交叉光束制造周期性温度分布，探测信号随时间衰减。硅薄膜原论文利用这类衰减测热传输，并比较光激发载流子前后的响应；该光学信号不是四探针电压，尽管自由载流子也会影响热传输过程。

图7. 硅薄膜的瞬态热光栅测量装置及光激发前后探测信号的衰减曲线。DOI：10.1038/s41467-020-19938-9

激光闪射通常测片材厚度方向的瞬态温升，瞬态热光栅可选择面内温度起伏的尺度；两种测量几何不能在各向异性薄膜上直接合并。电学四探针则沿电极排列方向测量电压降，电极间距、厚度及接触状态影响从电阻到 σ 的几何换算。

电导率通常从电压、电流及试样几何尺寸换算，四探针测量还可减小接触电阻的影响。比较这些物性时，应标明温度、试样方向和孔隙率；交流电学数据还需注明频率。多孔陶瓷的有效 κ 与 σ还包含孔隙和界面造成的传输阻力。