三维材料载流子迁移率计算

结合案例图分析，我把这个项目的完整过程、参数设置、结果解读和踩坑经验整理如下。

项目过程：三维材料载流子迁移率计算

这个项目的客户是半导体器件方向，要估算硅（Si）在不同温度和掺杂浓度下的电子/空穴迁移率，用于器件模拟参数校准。实验数据有，但客户希望从第一性原理给出一条” clean”的理论曲线，和实验做对比。

我采用的是形变势理论（Deformation Potential Theory）+ 玻尔兹曼输运方程（Boltzmann transport）的方案，软件用 VASP 算能带和形变势，再用 BoltzTraP 或自己写的脚本算迁移率。之所以没直接用 MD 算散射，是因为 MD 需要超大体系（>1000 原子）和超长轨迹，成本太高；形变势理论对三维块体半导体的声学声子散射已经能给出合理趋势。

参数设置：PBE，ENCUT=520 eV，K 点 12×12×12，计算了 Γ 点附近的能带曲率以得到有效质量；形变势通过 ±1% 单轴应变下的 VBM/CBM 偏移获得；温度范围 100–600 K，掺杂浓度 1×10¹⁴–1×10²⁰ cm⁻³。

案例分析：从数据到结论

第一张案例图是迁移率随温度的变化。纵轴是迁移率（cm²·V⁻¹·s⁻¹），横轴是温度（K），双对数坐标。红色曲线代表电子 μ_n，蓝色曲线代表空穴 μ_p。室温（RT，灰色阴影区）附近，电子迁移率约 1.4×10³ cm²·V⁻¹·s⁻¹，空穴约 5×10² cm²·V⁻¹·s⁻¹（与硅实验值 1350 / 480 cm²·V⁻¹·s⁻¹ 量级一致）。高温区电子迁移率近似按 T⁻²·⁴ 下降，这接近声学声子散射主导的 T⁻¹·⁵ 到 T⁻²·⁵ 范围，说明我们的模型捕获了主要散射机制。

第二张案例图是迁移率随掺杂浓度的变化。低掺杂区（<10¹⁵ cm⁻³）迁移率基本饱和；当掺杂浓度超过 10¹⁶ cm⁻³ 后，电子和空穴迁移率都开始明显下降，电子下降更剧烈，整体近似按 N⁻⁰·⁷ 变化。这与电离杂质散射的理论预期（Brooks-Herring 模型）一致。

总结教训：复盘必须警惕的三件事

第一，有效质量对迁移率影响极大。我用 12×12×12 K 点算 Γ 点附近曲率，得到电子有效质量约 0.19 m₀，和文献 0.19–0.26 m₀ 范围吻合；如果 K 点只用 6×6×6，曲率误差会让迁移率低估约 15%。

第二，掺杂浓度计算要考虑简并效应。浓度高于 10¹⁹ cm⁻³ 时，费米能级进入导带，简单的非简并近似会失效，需要改用费米-狄拉克积分。

第三，温度依赖的指数不是固定值。我得到电子是 T⁻²·⁴，空穴更接近 T⁻²·⁰，这与有效质量和形变势的差异有关，不能硬套 T⁻¹·⁵。

项目经验：下次再做我会怎么做

以后再做块体半导体迁移率，我会按这个顺序：1) 高精度能带（dense K 点）→ 2) 有效质量和形变势 → 3) 声子散射模型 → 4) 与实验曲线对标。如果客户要的是定量准确的数值，必须提醒 PBE 带隙低估会带来 10–30% 的系统误差，必要时用 HSE06 或 GW 修正带隙后再算形变势。

三维材料载流子迁移率计算

计算工程师团队：500+ ｜ 成功计算案例：5W+ ｜ 1W+ 助力顶刊发布经验 需要模拟代算服务可以联系我们，可免费评估项目的计算费用和免费出计算方案。