MS计算功函数：CASTEP中表面模型与静电势参考点选取的实战决策

MS计算功函数的操作步骤在CASTEP中不算复杂——建一个表面slab模型，加足够厚的真空层，做自洽计算后从LOCPOT文件提取静电势沿z方向的分布，功函数就是真空能级与费米能级之差。但”足够厚”的真空层是多少？静电势的”真空能级”取哪个位置？表面模型的终止面该怎么选？这些看似细节的问题，每一个都可能让功函数的计算结果偏移0.1-0.3 eV。

项目背景：二维材料功函数对异质结能带对齐的影响

团队处理过一个MoS₂/WSe₂范德华异质结项目，目标是判断两层材料堆叠后的能带对齐方式。客户做光电探测器，需要知道异质结是Type-I还是Type-II对齐——这直接决定器件的工作机制和量子效率。而功函数差异是判断能带对齐类型的决定性参数。

DFT计算功函数的第一步是分别计算单层MoS₂和单层WSe₂的功函数。团队用CASTEP的PBE泛函，截断能500 eV，k点网格12×12×1。MoS₂的功函数计算结果为5.12 eV，WSe₂为4.57 eV，差异0.55 eV。两个值看起来差距不大，但放在二维材料尺度上，0.55 eV的功函数差对应的内建电场足以在界面区域产生可观的电荷重分布。

真空层厚度：功函数计算中最容易被忽视的参数

功函数计算需要一个真空层来定义真空能级。对于二维材料，真空层厚度对功函数的影响比体相材料大得多——因为表面偶极的长程尾部衰减较慢，真空层不够厚时，两侧表面的偶极会通过真空层相互作用，导致真空能级被低估。

团队做了真空层厚度从12 Å到25 Å的收敛测试。12 Å时功函数偏低约0.08 eV，15 Å时偏差缩小到0.03 eV，18 Å以上基本收敛。考虑到单层MoS₂本身厚度约6.1 Å，18 Å的真空层意味着周期性单元总高度约30 Å。这个值看起来浪费，但功函数对真空层厚度的敏感度确实需要这个余量。

在异质结计算中，真空层的选取更棘手——两层材料的组合高度更大，需要的真空层也要相应增加。团队最终在异质结模型中使用了20 Å的真空层，周期高度约48 Å，k点密度降至9×9×1以控制计算量。

表面模型的终止面：不是随意切一刀

MS计算功函数中表面slab的构建需要选择终止面。对于MoS₂这类层状材料，终止面只有一种（S-Mo-S三层结构），不存在选择问题。但对于像TiO₂这样的体相材料，(110)、(101)、(001)三个低Miller指数面的功函数可能差异0.2-0.5 eV——因为不同面的表面偶极矩方向和大小不同。

团队在一个TiO₂功函数评估的辅助计算中验证了这一点：锐钛矿(101)面功函数4.83 eV，(001)面5.02 eV——差异0.19 eV来自(001)面更高密度的表面氧悬挂键导致的偶极矩翻转。如果客户需要的是特定晶面的功函数（比如薄膜生长的主导面），就必须建对应面的slab模型，不能随意选一个面算完当作”TiO₂的功函数”。

静电势参考点的选取：真空区中央不是唯一选择

从LOCPOT提取的静电势沿z方向的分布曲线，在真空区域呈现一个平台——这就是真空能级。功函数等于这个平台值减去费米能级。但”真空区中央”这个参考点选取方式，在偶极修正开启的情况下需要调整。

CASTEP中功函数计算通常需要开启偶极修正（Dipole Correction）——在真空层中央引入一个额外的势能跳跃，消除slab不对称导致的真空层两侧静电势差异。不开偶极修正时，如果slab上下表面不对称（比如一面是S另一面也是S但不同构型），真空层两侧的静电势平台值可能差异0.1-0.2 eV——这个差异直接导致功函数定义的模糊化。

开启偶极修正后，真空层两侧的静电势平台值趋于一致，功函数的提取变得更干净。但偶极修正的参数（修正方向和参考点）需要手动指定——CASTEP默认的修正方向是z轴正方向，如果slab的极化方向与默认不一致，修正反而会引入额外误差。

PBE泛函对功函数的系统偏差

PBE泛函的功函数计算结果通常比实验偏低0.1-0.2 eV，这是GGA自相互作用误差对静电势分布的影响。MoS₂的PBE功函数5.12 eV，文献实验值约5.2-5.3 eV；WSe₂ PBE值4.57 eV，实验约4.6-4.7 eV。绝对值偏差在0.1-0.2 eV量级，但功函数差（0.55 eV）的偏差通常小于0.05 eV——因为两个材料的系统性偏差方向一致，差值几乎不受影响。

这个项目最终认定PBE的功函数差值可以直接用于异质结能带对齐判断，不需要额外做HSE06修正——客户关注的核心是功函数差而非绝对值，PBE在这个维度上足够可靠。

回过头看，MS计算功函数的可靠性取决于三个前置决策：真空层厚度是否收敛、表面终止面是否对应目标晶面、偶极修正是否正确设置。这三个决策都有物理依据可循而非随意猜测。方法论，可在中找到系统整理的案例，也欢迎回到https://www.keyanxueshu.com/浏览全站内容架构。