Effective Mass Model for Thermoelectrics — A Comprehensive DerivationHeo M, Lee KH, Song J, et al. | Energy Environ. Sci., 2026, 19, 3195–3202DOI: 10.1039/d5ee06413a热电材料可将热能直接转化为电能是废热回收和固态制冷的关键技术。其性能由无量纲热电优值 zT S²σT/κ 衡量涉及塞贝克系数 S、电导率 σ 和热导率 κ 三个相互耦合的输运参数。传统上预测 zT 最大值和最优掺杂浓度需要数值求解费米积分并依赖霍尔效应测量载流子浓度——这些技术门槛限制了大量实验研究者对热电性能的准确评估。本文提出了一种基于 |S|热电势的绝对值的有效质量模型简化分析方法。其核心思想是在重掺杂简并半导体中|S| 是费米能级 EF 的单调函数因此可以作为掺杂水平的直接描述符。通过将传统 EM 模型中依赖于约化费米能级 η EF/kBT 的所有输运系数S, σ, κe重写为仅依赖于 |S| 的解析函数最终得到 zT zT(|S|; B) 的闭式表达式其中 B 是唯一的材料参数。该方法仅需最常规的热电测量数据S, σ, κ无需费米积分数值求解也无需霍尔效应测量。热电优值 zT 的定义Eq. (1) 热电优值定义zT 的三个构成参数 S、σ 和 κ 随载流子浓度 n 的变化呈现相互制约的三角关系S 和 σ 随 n 反向变化S 随 n 增加而减小σ 随 n 增加而增大而总热导率 κ κe κL 中电子贡献 κe 随 n 增加而增大通过 Wiedemann-Franz 定律 κe LσT晶格贡献 κL 近似与 n 无关。这种三重制约使得 zT 在某一最优掺杂水平达到峰值偏离该最优值任一方向都会导致 zT 下降。单抛物带模型的基本假设有效质量模型建立在单抛物带Single Parabolic Band, SPB模型之上其核心假设包括(1) 能带为抛物型E(k) ħ²k²/2m*其中 m* 为态密度DOS有效质量。(2) 载流子服从费米-狄拉克统计输运性质由费米-狄拉克积分描述。(3) 声学声子散射为主导散射机制弛豫时间 τ ∝ E^{-1/2}迁移率前置因子 μ₀ 为常数。(4) 材料为重掺杂简并/中等简并半导体|S| 可作为约化费米能级 η 的单调代理变量。尽管真实材料可能偏离这些理想假设如多带贡献、带非抛物性、电离杂质散射等但 EM 模型从大量实验数据中提炼出的经验关系使其在广泛的材料体系中仍具有出色的预测能力。此时提取的参数应理解为有效输运参数而非严格的能带常数。约化费米能级 ηη 是 EM 模型中最核心的中间变量。所有输运系数S, σ, κe, n均可表示为 η 的函数。η 的正负和大小直接反映了体系的简并程度η ≫ 0η 3对应简并半导体类金属行为η ≪ 0η -3对应非简并半导体类绝缘体η ≈ 0 对应中等简并。热电材料的最优性能通常出现在 η ≈ -2 到 2 的范围内对应于 |S| ≈ 100-300 μV/K。费米-狄拉克积分 Fj(η)费米积分是 EM 模型数学结构的基石。对于 SPB 体系所有输运系数均可通过不同阶数的费米积分及其比值表达。两个最重要的极限在简并极限η ≫ 0F_j(η) → η^{j1}/(j1)在非简并极限η ≪ 0F_j(η) → Γ(j1)·e^η。费米积分的比值 F_j(η)/F_{j-1}(η) 出现在 Seebeck 系数、Lorenz 数等关键量的表达式中是连接 η 和可观测量的桥梁。载流子浓度 n(η)在抛物带假设下载流子浓度 n 由态密度 g(E) ∝ E^{1/2} 在费米-狄拉克分布 f(E) 下的积分给出。h 为普朗克常数m* 为态密度有效质量。在简并极限η ≫ 0n ∝ η^{3/2} ∝ EF^{3/2}恢复自由电子气的经典结果。在非简并极限η ≪ 0n ∝ exp(η) ∝ exp(EF/kBT)恢复玻尔兹曼统计。塞贝克系数 S(η)Eq. (11) 塞贝克系数塞贝克系数 S 是费米面附近态密度不对称性的度量由 Mott 公式导出。正负号分别对应 p 型和 n 型半导体。S(η) 是单调递减函数在 η → −∞ 时 |S| → ∞在 η → ∞ 时 |S| → 0。这一单调性是该方法的基石——给定 |S|可以唯一地反解出 η从而将所有依赖于 η 的输运系数重写为 |S| 的函数。在简并极限下S 简化为著名的 Mott 公式简并极限下的 Mott 公式该式表明 S ∝ m*T/n^{2/3}即高有效质量、低载流子浓度和高温有利于产生大热电势。电导率 σ(η) 与加权迁移率 μw中 μ₀ 为与能带结构无关的迁移率前置因子在声学声子散射假设下为常数。为分离材料本征参数和载流子浓度依赖定义两个关键量μw 将迁移率 μ₀ 和态密度有效质量 m* 合并为一个参数。μw 的物理意义在于它直接决定了给定 |S|即给定 η下的电导率——μw 越大相同 |S| 下的 σ 越大。μw 是比 σ 或功率因子 PF 更可靠的电子输运品质因子因为它剥离了依赖于掺杂水平η的部分仅保留反映材料本征电子输运品质的部分。σE₀ 仅依赖于温度和基本常数不依赖于材料参数。在 300 K 时σE₀ 约为 3.02 × 10⁴ Ω⁻¹m⁻¹电子热导率 κe 与洛伦兹因子 LWiedemann-Franz 定律洛伦兹因子 L 在 SPB 模型中由费米积分比值给出洛伦兹因子的一般表达式在简并极限金属L → L₀ (π²/3)(kB/e)² ≈ 2.44 × 10⁻⁸ WΩK⁻²。在非简并极限L → 2(kB/e)² ≈ 1.48 × 10⁻⁸ WΩK⁻²。对于实际热电材料L 在两者之间过渡可由 |S| 通过经验公式精确估算该经验公式在 |S| 20 μV/K 范围内精度极高使得 κL 可以从总 κ 中准确分离κL κ − LσT。品质因子 B热电品质因子 B 是决定材料最大 zT 的唯一材料参数。其理论定义为Eq. (1) B 的定义B 的物理意义是加权电子输运能力与晶格热输运之比。在 300 K 的实用形式为Eq. (2) B 的实用表达式B 越大zTmax 越高。要实现 zTmax ≈ 1需要 B ≈ 0.4zTmax ≈ 2 需要 B ≈ 1.5zTmax ≈ 3 需要 B ≈ 3.5。在 300 K 和 κL ≈ 1 W/mK 的典型条件下B 0.4 对应 μw ≈ 600 cm²/Vs。关键步骤变量代换 η → |S|以上推导表明σ(η) 和 κe(η) 均通过费米积分依赖于 η。而 S(η) 是 η 的单调函数。因此可以通过以下逻辑链将 zT 重写为仅依赖于 |S| 的函数步骤 1给定 |S|通过 S(η) 的反函数得到 η S^{-1}(|S|)。由于 S(η) 严格单调反函数唯一存在。步骤 2将 η 代入 σ(η) 和 κe(η)得到 σ(|S|) 和 κe(|S|)。步骤 3代入 zT S²σT/κ得到 zT zT(|S|; B)。在实用层面作者通过数值拟合将费米积分比值直接近似为 |S| 的解析函数绕过了显式求解 η 的步骤。定义公共分式 f(|S|)Eq. (5) 加权迁移率 μw 的 |S| 表达式Eq. (6) 电导率 σ 的 |S| 表达式Eq. (5) 和 Eq. (6) 互为逆关系。Eq. (5) 从 σ 和 |S| 计算 μw剥离掺杂依赖Eq. (6) 从 μw 和 |S| 预测 σ评估 μw 改善对 zT 的影响。zT 的 |S| 表达式 — 核心结果Eq. (7)zT 的热电势表达式分子 S² ∝ η² 代表热电驱动力分母第一项 (1/B)×f(|S|) 代表电子输运的电阻分母第二项 1.34×[1.5exp(−|S|/116)] 代表电子热导的损耗其中 1.5exp(−|S|/116) 即 L 的经验表达式。在扩散热电势主导且双极效应可忽略的条件下最为可靠。对于窄带隙材料在高温下的应用计算得到的 zT 应视为近似描述符。有效质量 m* 的提取Eq. (8) 有效质量的提取虽然 EM 模型的核心分析zT 预测和优化不需要显式知道 m*但若需要区分电子结构效应和散射效应对输运的影响可以利用霍尔效应测量RH 为霍尔系数提取 m*。在理想 SPB 体系中m*/me 是一个与掺杂浓度和温度无关的材料常数。但在实际材料中m*/me 随掺杂、温度和合金化而变化反映了能带结构的变化如带收敛效应或主导散射机制的变化。基于 EM 模型的数据分析流程标准分析步骤EM 模型提供了一套仅需 S、σ 和 κ 即可完成的热电性能分析流程步骤 1 — 计算 κL由 Eq. (4) 从 |S| 估算 L再由 κL κ − LσT 分离晶格热导率。步骤 2 — 计算 μw由 Eq. (5) 从 σ 和 |S| 计算加权迁移率。步骤 3 — 计算 B由 Eq. (2) 从 μw 和 κL 计算品质因子。步骤 4 — 预测 zTmax利用 Eq. (7) 绘制 zT(|S|) 曲线读取峰值。步骤 5 — 确定最优 |S|zT 峰值对应的 |S| 即最优热电势。若当前 |S| 偏离此值则需调整掺杂水平。热电势 |S||S| 的绝对值大小直接指示材料的掺杂状态|S| 300 μV/K材料处于欠掺杂状态应增加载流子浓度以提升 σ。|S| 200 μV/K材料处于过掺杂状态应降低载流子浓度以提升 |S|。|S| ∝ T金属/简并半导体行为表明 EM 模型参数可靠。|S| 随 T 升高而降低表明双极传导开始EM 模型预测的可靠性下降。|S| 出现峰值|S|max 在 Tmax可估算热学带隙 Eg ≈ 2|S|maxTmax。加权迁移率 μwμw 的大小和温度依赖性提供了丰富的输运机制信息μw 的绝对值μw ≈ 600 cm²/Vs 是室温下实现 zTmax ≈ 1 的典型门槛κL ≈ 1 W/mK。μw ∝ T^{-3/2}声子散射主导符合 SPB 模型的基本假设。μw 随 T 升高而升高非声子散射晶界势垒、电离杂质主导。例如Mg₃(Sb,Bi)₂ 中 1 μm 晶粒尺寸的样品表现出随 T 升高的 μw归因于晶界势垒散射。μw 在固定温度下随不同样品变化 20%差异主要来自掺杂水平而非本征输运性质的变化。μw 在固定温度下随不同样品变化 50%差异来自能带结构或散射机制的根本变化。功率因子优化功率因子 PF S²σ 在 |S| ≈ 167 μV/K 处达到峰值与材料无关。然而zT 的峰值出现在 |S| ≈ 200-250 μV/K取决于 B 值。这意味着仅基于 PF 的优化会系统性导致过掺杂使 zT 低于理论最大值达 25%对于 B ≈ 0.4 的材料。这是因为 PF 优化忽略了 κe 的贡献——在 |S| 较低过掺杂时σ 虽然大但 κe 也大导致 zT 被抑制。模型验证与应用案例图1 EM 模型框架(a) 不同 B 值下的 zT(|S|) 曲线族(b) B-zTmax 对应关系(c) 标准分析流程示意图1 综合展示了 EM 模型的核心框架。(a) zT 作为 |S| 的函数不同 B 值对应不同的 zT 曲线族。每条曲线的峰值橙色圆点给出该 B 值下的理论 zTmax 和最优 |S|。(b) 将 (a) 中的峰值点映射为 B-zTmax 曲线建立了品质因子与最大优值之间的直接对应关系。(c) 展示了从 S/σ/κ 三参数出发依次计算 L → κL → μw → B → zTmax 的完整分析流程。图2 (a) 功率因子 PF 与 |S| 的关系(b) B-zTmax 曲线与实验数据对比图2(a) 展示了功率因子 PF 作为 |S| 的函数。PF 在 |S| ≈ 167 μV/K 处达到峰值与材料无关——这解释了为什么仅基于 PF 的优化往往导致过掺杂。图2(b) 将多种热电材料的实验数据与理论 B-zTmax 曲线对比验证了 EM 模型在 Bi₂Te₃ 基、PbTe 基、Mg₃Sb₂ 基、half-Heusler 等多种材料体系中的普适性。图3 μw 的分析应用(a) (Bi,Sb)₂Te₃ 中带收敛对 μw 的增强(b) Mg₃(Sb,Bi)₂ 中晶粒尺寸对 μw 的影响图3 展示了 μw 在识别输运机制变化中的诊断能力。(a) (Bi₀.₂₅Sb₀.₇₅)₂Te₃ 的 μw 比 Bi₂Te₃ 高 50% 以上归因于 Bi₂Te₃-Sb₂Te₃ 合金化引起的带收敛效应——多个能带汇聚增大了 md*而 μw ∝ (md*/me)^{3/2}。(b) Mg₃(Sb,Bi)₂ 中晶粒尺寸从 1 μm 增大到 30 μm 降低了晶界密度使 μw 升高Nb 掺杂进一步抑制了晶界势垒μw 甚至高于大晶粒样品。图4 基于 EM 模型的 zT 预测(a) n 型 Mg₃(Sb,Bi)₂(b) p 型 (Bi,Sb)₂Te₃(c) n 型 PbTe(d) p 型 SnSe图4 展示了 EM 模型在多种代表性热电材料中对 zT 的预测精度。仅需 S、σ 和 κ 的实验数据无需霍尔效应测量Eq. (7) 即可准确再现 zT 随 |S| 的变化趋势并预测最大可达 zT 和最优掺杂水平。在 Mg₃(Sb,Bi)₂、PbTe、SnSe 等体系中预测值与实测值高度一致验证了 EM 模型的广泛适用性。讨论适用条件(1) 扩散热电势主导假设热电输运由载流子扩散驱动声子曳引效应可忽略。在 T 100 K 时声子曳引可能显著EM 模型参数的物理意义变得模糊。(2) 双极效应可忽略在窄带隙材料的高温区少数载流子的激发导致 S 被部分抵消、κ 增加双极热导此时提取的 κL 和 B 被低估。(3) |S| 20 μV/Kμw 的解析表达式在该范围内精度最高。(4) 重掺杂简并/中等简并半导体|S| 作为 η 的单调代理变量最为可靠。有效参数的解释当体系偏离理想 SPB 行为如多带贡献、带非抛物性、非声学声子散射时EM 模型提取的 μw、m* 和 B 应理解为有效输运参数而非严格的能带常数。这些有效参数仍然具有实用价值——它们能够正确预测 zT 的优化方向和最大可达值但不应过度解读其微观物理含义。与霍尔效应方法的对比传统 EM 模型需要霍尔效应测量来确定 n 和 m*。本文的 |S| 基方法消除了这一要求但霍尔效应测量在以下场景中仍有价值(1) 区分电子结构效应m* 变化和散射效应μ₀ 变化对 μw 的贡献(2) 验证单带假设的合理性通过霍尔因子 rH 的温度依赖性(3) 在多相或非均匀材料中提供独立的载流子浓度信息。核心结论1. 本文建立了基于 |S| 的有效质量模型分析方法将传统上需要数值求解费米积分和霍尔效应测量的热电性能评估转化为仅需常规热电测量S, σ, κ即可完成的简洁流程。2. 模型的核心数学技巧在于利用 S(η) 的单调性将 zT 表达式中所有依赖于 η 的项通过数值拟合转换为 |S| 的解析函数最终得到 zT zT(|S|; B)。3. 加权迁移率 μw 是比 σ 或 PF 更可靠的电子输运品质因子其温度依赖性可识别主导散射机制声子 vs 晶界 vs 电离杂质。4. 品质因子 B 直接决定了材料的理论最大 zTB 0.4 → zTmax ≈ 1B 1.5 → zTmax ≈ 2B 3.5 → zTmax ≈ 3。5. 热电势 |S| 本身即是最优掺杂水平的诊断工具|S| 300 μV/K 指示欠掺杂|S| 200 μV/K 指示过掺杂。功率因子优化会系统性导致过掺杂使 zT 低于理论最大值达 25%。6. 该框架已在多种热电材料体系Bi₂Te₃ 基、PbTe 基、Mg₃Sb₂ 基、SnSe、half-Heusler、方钴矿等中得到验证适用于从实验合成到第一性原理计算再到数据驱动材料筛选的完整研究链条。