功函数vaspkit计算

结合案例图分析，本文以这个异质结能带对齐项目为例，讲 VASPKIT 算功函数的便利和 slab 模型的坑。这是一个半导体异质结能带对齐的支持项目。客户在做光电探测器，需要知道几种半导体功函数以判断界面能带偏移。第一张案例图是 VASPKIT task 426/427 算出的表面静电势曲线，蓝色 planar average 在内部有 Friedel 振荡，红色 macroscopic average 把短周期振荡抹平。E_vac≈7.1 eV、E_F≈4.1 eV，功函数 Φ=3.0 eV。第二张案例图是 Si、Ge、GaAs、InP、ZnO 五种的功函数对比。

项目过程：异质结能带对齐项目的执行链路

方法上 VASPKIT 的 426/427 自动读 LOCPOT 做平面和宏观平均，比手写 Python 处理 LOCPOT 省事太多，但前提是 slab 模型到位。我所有体系统一了泛函（PBE）、赝势、截断能（500 eV）和真空层（20 Å），否则比较毫无意义——不同参数下功函数可差 0.2 eV 以上。 关键设置：真空层 20 Å 是收敛测试定的，12 Å 时两侧 E_vac 平台差 0.15 eV，15 Å 收窄到 0.05，20 Å 后 < 0.02 eV；表面在 z 居中避免镜像相互作用；极性表面开 LDIPOL。客户要一组可比数据，参数统一是底线，我据此把五种半导体逐一算完。

案例分析：从静电势到可比功函数数据

图里读出的数要能落地。第一步读 macroscopic average 的真空平台 E_vac 减 E_F 得 Φ，planar average 内部的 Friedel 振荡必须靠宏观平均抹平才能稳取。我对真空区多个 z 切片算 E_vac 标准差确认 < 0.02 eV 才采信，否则功函数随切片跳 0.05 eV。 第二步把五种半导体 Φ 列成表：GaAs 5.22±0.05 eV 最高、ZnO 4.45±0.05 最低，差值 0.77 eV 对异质结能带偏移是真有物理意义的。误差量级上，宏观平均窗口取一个晶面间距（如 Si 约 5.4 Å），窗口偏差会让 E_vac 读值偏 0.03 eV，必须固定。

总结教训：宏观平均与偶极修正的坑

做的过程中我踩过实打实的坑。第一坑是真空层：12 Å 时 Φ 比收敛值低 0.15 eV 且两侧平台不平，15 Å 偏差 0.05 eV，20 Å 稳定在 3.0±0.02 eV——这条收敛曲线让客户看到”为什么 20 Å 是必须的”。 第二坑是 macroscopic average 窗口宽度，通常取一个晶面间距，窗口太大把真实电势变化也抹掉、E_vac 读值偏差 0.03 eV，太小则平滑不彻底仍有振荡。第三坑是极性表面忘开 LDIPOL，真空两侧平台差 0.10 eV，功函数被带偏；加上后两侧对齐。

项目经验：VASPKIT 对比交付的可复用清单

作为可复用的项目经验，第一，所有体系必须同一套参数，否则比较无意义。第二，半导体功函数要和能带对齐、电子亲和能一起交付，单给 Φ 客户用不了。第三，误差带 ±0.05 eV 必须标注。客户最终选了高功函数的 GaAs 做界面接触层。说到底，VASPKIT 把”算功函数”从体力活变成点几下，但 slab 模型不过关，点出来的也是漂亮的错误。

功函数vaspkit计算

计算工程师团队：500+ ｜ 成功计算案例：5W+ ｜ 1W+ 助力顶刊发布经验

需要模拟代算服务可以联系我们，可免费评估项目的计算费用和免费出计算方案。