vasp功函数怎么计算：从静电势提取到真空能级对齐的完整路径

vasp功函数怎么计算，表面物理和器件界面的同行问得最多。功函数的定义很简单：Φ = VAC − EF，即真空能级减去费米能级。但真要在VASP里算准，难点全在”真空能级怎么取”——静电势在表面法向上有振荡，直接读最大值会系统性偏高。我做过Au(111)和Si(100)两个体系，把每一步的坑都踩了一遍。

slab构建时真空层不够功函数直接废

功函数依赖表面的静电势渐近线，而静电势要”跑到真空”才平稳。我最初给Au(111) slab的z方向真空层只留了10Å，静电势在真空区还没拉平就开始周期性重复，平均电势提取出来的”真空能级”其实还带着表面镜像的尾巴，功函数算出来比实验低0.3eV。

把真空层加到18Å后，z方向静电势中段明显出现一段平坦平台，平台值才可信。可复用要点：功函数计算的真空层至少15–20Å，比普通吸附计算（12–15Å）更厚，因为吸附计算关心的是吸附质能量，功函数关心的是体相到真空的势垒，对真空区的平稳度要求更高。

真空层之外，slab厚度对功函数也有要求，但逻辑和吸附不同：功函数关心的是表面偶极，而表面偶极来自最外层原子弛豫和电荷重新分布，一般4-6层就收敛。我对比过Au(111)用4层和6层，功函数差0.02eV，6层已稳；但少于4层时，两面原子的弛豫会穿透整个slab，上下表面互相干扰，功函数能偏0.1eV。所以功函数slab的判据是”弛豫不穿透”，而非”体相长好”——这和吸附能要厚slab体相电子结构收敛的侧重点不一样，值得新手分清。

偶极修正（IDIPOL）开不开差别巨大

在vasp功函数怎么计算里，表面slab天然带偶极矩（体相截断在表面产生电荷非补偿），不加偶极修正时，周期性边界条件下静电势会线性倾斜，整个势能面被”拧”歪。我第一次算Si(100)功函数，IDIPOL没开，真空区电势斜率明显，提取的平台值虚高0.4eV。

开IDIPOL=3（z方向偶极修正）后，静电势变成水平渐近线，平台值稳定。这条铁律：做表面功函数、表面能、任何依赖真空能级的计算，IDIPOL必须开。顺便，沿z方向的k点只需1个，但面内k点要密（Au(111)用9×9×1），否则费米能级EF（从DOS读）会因k点稀疏抖0.05eV，传导到功函数。

从LOCPOT提取静电势的两步法

VASP跑完静态计算（LAECHG=.TRUE.或写LOCPOT），静电势存在LOCPOT里。我用VTST的vaspkit或自己写的脚本，沿z方向把每个xy平面上的电势做平均，得到V(z)一维曲线。关键在”平均化”——单点读数受原子平面振荡干扰，平面平均后取真空区平台值作为VAC。

然后EF从OUTCAR的E-fermi读（静态计算输出的费米能级）。Φ = VAC − EF。Au(111)实验功函数约5.31eV，我算得5.28eV，偏差0.03eV，在PBE误差带内。Si(100)实验约4.6eV（n型表面态影响），PBE算得4.55eV，吻合。数字站出来说话：平台取法对了，偏差就在方法误差内。

平面平均这一步有个细节能决定VAC取没取对：LOCPOT存的是原始三维网格电势，直接按z索引平均会受网格数（NGX NY NZ）影响，必须先把电势重采样到均匀间隔再平均。我用的脚本先把LOCPOT读成三维数组，沿z做每一层的xy平均，再对z画曲线；取平台时不是取最后一个点，而是取平台段（通常真空区中段1/3）的平均值减掉微小线性漂移。曾有一次图省事取了z最大那一点的电势，结果该点正好在周期性边界的振荡处，VAC虚高0.15eV。平台取中段平均，是功函数算准的最后一公里。

半导体功函数比金属多一层麻烦

金属费米能级在能带里明确，EF直接减VAC即可。半导体麻烦在：块体半导体费米能级位置依赖掺杂，VASP默认算的是本征半导体，EF落在带隙中间，算出的”功函数”其实是电子亲和势+带隙一半的混合物，不是实验测的功函数（实验功函数受表面态和Pinning强烈影响）。

我在Si(100)上专门做了对照：本征计算EF在带隙中，Φ=4.55eV；若强行把EF钉在价带顶（模拟p型重掺杂），Φ变成4.9eV接近实验。结论：半导体功函数计算必须说明掺杂条件，否则数字没有对标意义。这一点很多教程不提，导致新手算出来和实验对不上就慌。

再补一个半导体场景的坑：表面重构。Si(100)实际表面是2×1二聚体重构，若不建重构直接算理想(1×1)终止，表面态密度完全不同，功函数能差0.2eV以上。我在交付Si功函数时，先弛豫出2×1重构，再算功函数，才和实验的4.6eV对齐。所以半导体功函数的另一半诚实边界是：表面终止和重构必须和实验制备条件一致，否则算的是”理想清洁面”而非”实际表面”。这条在氧化物（如TiO2(110)的桥氧缺失）上同样致命。

三个最容易被忽略的坑

第一，slab两面的不对称性。如果slab上下表面终止不同（如一面重构、一面未重构），偶极修正只能校正净偶极，两面功函数会不同——这时要分别取两面的真空平台，得到两个功函数。第二，结构没充分弛豫。表面弛豫改变最外层原子间距，直接改表面偶极，功函数可偏0.1eV，必须ISIF=3充分优化。第三，真空层不够厚时偶极修正也救不回来，平台取不到真值。

回过头看，功函数计算的诚实边界

回过头看，vasp功函数怎么计算这件事，技术门槛不在公式（Φ=VAC−EF谁都懂），而在让VAC可信的那串前置设置：厚真空、开偶极修正、面内密k点、充分弛豫、平面平均取平台。任何一道松手，VAC就带假，功函数必偏。

vasp功函数怎么计算的可复用要点：真空≥18Å、IDIPOL=3、面内k点≥9×9、结构ISIF=3优化、LOCPOT平面平均取平台。这条链路在金属、氧化物、二维材料上都稳。它的边界在半导体/绝缘体——那时功函数受表面态钉扎，计算值要和具体表面终止、掺杂、温度下的实验条件对齐才有意义，单纯体相计算只能给一个参考值。

功函数算准后还有个实用延伸：用功函数反推表面能。表面能可以从slab两侧真空平台的差值和体相能量推出，是薄膜生长、催化位点稳定性的关键量。我常在功函数任务里顺手算表面能，给客户多一个维度判断表面稳定性。一次算出金红石和锐钛矿的(101)表面能差，正好解释客户实验中锐钛矿更易暴露(101)的现象。

功函数这件事，门槛不在公式而在前置设置。厚真空、偶极修正、密k点、充分弛豫，四道关过完，VAC才可信。我把这套前置做成固定清单，新人照着勾就不会出大错——功函数计算的可靠性，本质是流程可靠性。