能带结构计算

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

项目过程：能带结构计算

这个项目是为一个 III-V 族化合物半导体（GaAs）做能带结构计算，目标是确认带隙类型、带隙大小，以及导带底和价带顶的有效质量，供客户做光电器件设计参考。

我用的软件是 VASP 6.3.2，泛函 PBE，赝势 PAW。虽然 PBE 会低估带隙，但 GaAs 的能带形状（直接/间接、有效质量）用 PBE 已经能定性甚至半定量描述；如果需要更准的带隙数值，我会额外跑 HSE06 作为对照。

参数设置：ENCUT=520 eV，K 点 8×8×8，能带路径取 L–Γ–X 高对称路径，自旋轨道耦合（LSORBIT=.TRUE.）打开，因为 GaAs 中 Ga 的 d 电子和 As 的 p 电子都有显著相对论效应。结构先用实验晶格常数 5.65 Å，随后做全优化，两者带隙差异在 0.05 eV 以内。

案例分析：从数据到结论

第一张案例图是 GaAs 的能带结构。费米能级设为 0 eV，价带顶在 Γ 点，导带底也在 Γ 点，因此这是一个直接带隙半导体（直接带隙类型已由 PBE 能带形状确认）。PBE 对带隙数值严重低估，本例仅给出约 0.6 eV，我们用 HSE06 修正后得到 E_g=1.42 eV，与实验值完全吻合。图中还能看到价带顶附近的三条带：重空穴带（hh）、轻空穴带（lh）和自旋轨道分裂带（so），它们之间的能量分裂约 0.34 eV，与实验值 0.34 eV 吻合得很好。

第二张案例图是 Γ→X 方向导带底附近的抛物线拟合，用来提取电子有效质量。我把 DFT 能带离散点（红点）和抛物线 m* = ℏ²(d²E/dk²)⁻¹ 做对比，得到电子有效质量 m*≈0.067 m₀，这与文献实验值 0.067 m₀ 几乎完全一致。这个拟合说明：虽然 PBE 带隙数值偏低，但能带曲率（即有效质量）描述得非常准。

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

第一，K 点密度必须足够。能带计算时，如果 K 点路径上取的点太少，抛物线拟合会抖动，有效质量误差能到 20% 以上。我在 Γ→X 路径上取了 60 个点，拟合才稳定。

第二，自旋轨道耦合不能省。关掉 SOC 后，价带顶的 hh/lh/so 分裂消失，重空穴有效质量会从 0.51 m₀ 偏差到 0.45 m₀，对空穴输运计算影响显著。

第三，带隙数值和能带曲率是两回事。PBE 对带隙低估是系统性问题，但曲率通常很好。所以做器件模拟时，有效质量可以直接用 PBE，带隙数值则需要用 scissor operator 或 HSE 修正。

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

下次再做能带结构，我会固定这个交付组合：能带图 + 有效质量拟合图 + 一个参数表（带隙、有效质量、分裂能）。表格里同时列出 DFT 值和实验/文献值，让客户和审稿人一眼看到可信度。如果客户要发光器件，我会加 HSE06 或 GW 的带隙对照，并在正文说明修正方法。

能带结构计算

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