功函数是什么？
DFT 中的功函数定义为电子从材料费米能级移动到表面外真空区所需的能量，常写作 Φ = Evac – EF。
这里的 Evac 并非代码自动给出的绝对零点，它来自 slab 模型外侧真空层中的静电势平台；EF 是同一个自洽计算得到的电子化学势。只有这两个量取自同一套 slab 计算，差值才对应该表面的功函数。
slab 模型把晶体沿某个晶面切开，在表面法向加入真空层，再对 平面平均静电势 V(z) 做读取。靠近原子层的位置，V(z) 会随着离子核和电子密度强烈起伏；进入足够宽的真空区后，曲线应出现近似平坦的区域，这个平坦值就是 Evac 的来源。
如果真空层中还保持明显斜率，常见原因是真空太窄、slab 上下表面不等价、吸附层带来净偶极，或电荷修正没有处理好。
图1. 高通量表面功函数数据库的构建流程，包含体相结构筛选、slab 生成、终止层识别与 DFT 计算。DOI：10.48550/arXiv.2011.10905
把这套流程放到功函数问题中，读者能看到体相结构、晶面、表面终止、slab、DFT 自洽电势和最终 Φ 之间的次序。
材料名无法固定 Φ。Ni、Pt、MoS2 或 TiO2 这类名称只告诉我们化学组成，无法替代表面取向和终止层。真正进入公式的是某一个具体 slab 的真空电势与 EF。
对称 slab 的上下表面等价时，两个真空平台通常相同或很接近；非对称 slab 会给出 两个外侧真空平台，对应两侧不同表面。
此时只写一个 Φ 会遮住表面来源，较好的写法是分别给出 top surface 与 bottom surface 的值，或者构造双面对称模型来避免净偶极。
若研究对象本来就是极性表面、吸附分子或带外电场的界面，净偶极是物理问题的一部分，文中应说明电势平台取自哪一侧。
图2. 对称与非对称 slab 中表面终止数的判别示意，体现上下表面是否等价对 Φ 的影响。DOI：10.48550/arXiv.2011.10905
上下表面等价时，slab 中心相当于镜面对称参考；上下终止不同，真空两端会感受到不同表面偶极。
原子层外侧的电子尾部、离子层排列、配位不饱和位点、吸附物取向都会改变表面偶极，Evac 跟着移动，而 EF 未必同步移动。功函数的数值变化常在这一段空间里发生。
电势图谱若没有可识别的真空平台，公式本身也给不出可复查的 Φ。没有真空区的体相周期模型没有外部 Evac，直接谈 Φ 会失去参照。
slab、薄膜、二维材料、带外表面的界面模型才自然拥有真空参照；纯体相导体更适合讨论 EF、DOS、载流子浓度或接触势，表面功函数不适合作为主结果。
金属表面：费米能级清楚时，Φ 可以直接比较
金属和准金属最适合直接计算功函数，因为它们在 EF 附近有可占据电子态，电子化学势在自洽计算中定义清楚。
研究金属电极、场发射、热电子发射或金属/二维半导体接触时，Φ 的物理含义也比较直接：数值越大，电子从该表面进入真空所需的能量差越大。
同一类表面条件指的是同一晶面、同一终止和同一吸附状态。
金属的 Φ 会随晶面变化，(111)、(100)、(110) 的表面原子密度、配位不饱和程度和电子溢出不同；同一晶面的 O、S、H、卤素或有机分子吸附，也会改变表面偶极。计算里常见的偏差来自 slab 层数未收敛、真空平台不平、表面弛豫不足，以及非对称 slab 中两侧表面混用。
图3. DFT slab 数据中的功函数与表面能分布，Φ 的主体区间集中在 3 到 5 eV 附近。DOI：10.48550/arXiv.2011.10905
大量金属和准金属表面的 Φ 在高通量数据中主要分布在几 eV 范围内，但两端数值可以相差很大。晶面、终止、表面能和电子结构共同决定这个分布，元素周期表上的单一趋势解释不了全部结果。
低表面能的晶面在热力学上更容易出现，低 Φ 或高 Φ 的终止层则未必就是最稳定终止。
金属表面比较功函数时，文中至少要写出 材料、晶面、终止层、弛豫状态、真空层和偶极修正处理。若研究吸附引起的 Φ 变化，应把 ΔΦ 写成同一表面在吸附前后的差值，不要拿两个不同 slab 的绝对值硬比。
吸附物电负性较强时，电子密度往表面外侧或吸附原子附近重排，真空侧电势台阶会改变；碱金属覆盖时，电子外溢和表面偶极方向又会把 Φ 推向另一侧。
图4. 顶层和次表层元素类型对应的功函数趋势，表面化学成分会改变真空能级台阶。DOI：10.48550/arXiv.2011.10905
把 顶层元素和次表层元素 分开统计后，功函数对最外层原子的敏感程度会变得很直观。顶层出现碱金属、碱土金属时，低 Φ 表面比例增加；顶层富含 O、卤素、C、N、S、Se 时，高 Φ 表面更常见。
单个元素类别无法决定全部结果，最外层电荷分布、局部偶极和金属态屏蔽 共同改变真空平台。
在金属催化或电化学界面中，Φ 还常和电极电势、表面电荷、双电层联系起来。真空 slab 的 Φ 反映干净或吸附表面的电子逸出能；含水、离子、电场的模型中，外侧参照已经不再是理想真空，电势降分布在金属、吸附层和溶剂层之间。
此时直接把真空功函数替换成实验电极电势，会把参照体系改掉。较清晰的写法是分别报告真空 slab 的 Φ、界面电势差和所用电位标尺。
半导体表面：Φ 能算，但费米能级位置要写明
在理想本征半导体的 slab 里，EF 位于带隙区域，那里没有像金属那样连续的可占据电子态。
计算程序给出的 EF 常由占据数、展宽方式和带隙中点约定影响。若表面有缺陷、吸附物、掺杂原子或表面态，EF 可能受这些态钉扎，Φ 才对应一个特定表面的电子化学势。
半导体表面更常和两个量一起读：电离能 IP = Evac – EVBM，电子亲和能 EA = Evac – ECBM。
IP 描述价带顶到真空能级的距离，EA 描述导带底到真空能级的距离。对光催化、半导体接触、电子抽取层和空穴传输层来说，EVBM 与 ECBM 往往比单个 Φ 更能说明载流子注入方向。
图5. 少层与体相 TMD 的带边能量以真空能级为零点排列，VBM 和 CBM 可分别对应电离能与电子亲和能。DOI：10.1103/PhysRevB.103.085404
少层 TMD 的带边可以直接放在真空能级下面读取。MoS2、MoSe2、MoTe2、WS2、WSe2、WTe2 的 VBM 和 CBM 随层数移动，IP 与 EA 随之改变。若只给出一个 Φ，读者无法判断变化来自价带顶、导带底、费米能级位置，还是表面偶极。
PBE、HSE 和 GW 给出的带边位置并不相同，PBE 低估带隙时，EVBM、ECBM 相对真空能级的位置会随之偏移；如果文章讨论光吸收、能带对齐或 Schottky 势垒，带边误差会传到 IP、EA 和势垒高度。
真空平台的读取仍来自静电势，但带边位置会受泛函和准粒子修正影响。写作时应把“PBE slab 功函数”和“HSE/GW 带边相对真空能级”分清。
图6. Mo、W 二硫族化物少层到体相的功函数和带隙变化，厚度会改变带边与 Φ。DOI：10.1103/PhysRevB.103.085404
TMD 功函数随层数变化并不剧烈，但单层和多层仍有差别；带隙变化更明显。二维半导体的屏蔽弱，真空层、层数、偶极修正和层间距都会影响电势平台。
对于单层材料，表面就是整张薄片；对于多层薄片，两侧表面可能相同，也可能因基底、吸附层或外电场而分裂成两个不同参照。
半导体表面报告 Φ 时，可以分成两种写法。研究 有明确掺杂、缺陷浓度、表面态或电极接触的实际表面，可以报告该模型的 Φ，并说明 EF 由哪些态决定。研究本征半导体的能带对齐，IP、EA、EVBM、ECBM 与真空能级的相对位置更透明，Φ 只能作为附加信息。
异质结：接触后的势垒来自界面电势重排
拿 两个孤立材料的功函数相减，再直接写成接触后的电荷转移量或内建电场强度，是异质结讨论里常见的误用。
孤立表面 Φ 确实能提供初始判断，特别是金属/半导体接触或弱耦合 vdW 结；但两种材料真正接触后，费米能级会重新排列，界面处发生电荷重排，局部偶极、电势降和带弯曲都会改变原来的真空参照。
弱相互作用的二维异质结常用 真空能级对齐 作为初始模型：分别计算两个单层材料的 EVBM、ECBM 和 Evac，再比较带边相对位置。这种处理适合筛选 Type I、Type II、Type III 带边排列，尤其在层间杂化较弱、界面距离较大、没有强共价键重构时较有用。
图7. 两种半导体形成结时的三类带边排列，真空能级可作为弱耦合模型的外部参照。DOI：10.1103/PhysRevB.103.085404
三类带边排列来自真空能级下的 VBM 和 CBM。
Type I、Type II、Type III 分别对应包覆型、交错型和断裂型带隙：
Type I 中一个材料的带隙包含另一个材料的带边，电子和空穴倾向留在同一侧；
Type II 中 VBM 和 CBM 分属两种材料，更利于空间分离；
Type III 常与隧穿器件和强载流子转移联系起来。这里讨论的是带边位置，单个 Φ 的大小排序无法替代它。
若异质结已经建成 显式界面超胞，计算对象就变了。周期界面带真空层时，体系只有一个共同 EF，但两侧外表面可以有不同 Evac，于是出现侧向功函数或表观功函数。若模型没有真空层，例如体相 A/B 埋入式界面，外部真空参照消失，文章应改用宏观平均静电势、核心能级对齐、局域 DOS 和带边偏移来讨论势垒。
图8. Ti3C2T2/graphene 界面差分电荷密度，表面终止改变界面偶极和电荷重排。DOI：10.48550/arXiv.1902.00776
MXene/graphene 的例子中，Ti3C2T2 的 O、OH、F 终止改变界面附近的电子积累和耗尽区域，也改变 graphene 的狄拉克点相对 EF 的位置。功函数差可以提示电子流动方向，但界面形成后，真正进入势垒的是 电荷密度重排、界面偶极和能带杂化 后的电势台阶。
对金属/半导体接触，孤立金属 Φ 与半导体 EA 可给出理想 Schottky-Mott 图像；但真实界面常有金属诱导能隙态、化学键、缺陷态和界面氧化层，势垒会偏离简单差值。
对半导体/半导体结，孤立带边对齐能筛选候选组合；强耦合界面则要读取接触后局域势和投影能带。对催化异质界面，Φ 适合描述电子供受倾向，反应速率仍要回到吸附构型、反应能垒、d 带中心、轨道杂化和中间体稳定程度。
一个可复查的功函数结果，写法应具体到“哪一个表面、哪一侧真空、哪一种 EF、哪一套带边”。金属表面可直接比较 Φ；半导体表面应补上 IP、EA 或带边位置；异质结则要区分 接触前的真空对齐 和 接触后的界面电势重排。这些参照写清后，功函数就从孤立数字变成连接表面电子结构、界面偶极和能带排列的可计算量。