说明：本文华算科技主要介绍功函数的定义、DFT slab 中真空能级的读取方式，以及表面终止、空位和氧化层怎样改变功函数。
表面如何决定材料功函数
功函数 Φ 描述的是电子从材料表面移动到真空静止态所需要跨过的能量差。对金属和许多导电表面来说，常用写法是 Φ = Evac – EF，其中 EF 是体系的费米能级，Evac 是表面外侧真空区的电子能量参照。这个定义里必须盯住两个字：表面。它不是单独由化学式决定的体相常数，同一种 TiN、MoS2 或石墨烯，只要暴露晶面、吸附层、缺陷浓度、覆盖度发生变化，表面偶极和真空能级台阶就可能跟着改变。
把电子移到真空，并不是把电子拿到无限远的抽象空间里。计算里通常要求真空区足够厚，使表面外侧的静电势进入平台；实验里 UPS、Kelvin probe 或 KPFM 也都要说明参考探针、样品表面状态和环境。功函数读的是某个表面的电子逸出能量，而不是“材料整体电子越多或越少”的直接刻度。金属表面可以用 EF 作为清楚的电子化学势；半导体表面还要写清费米能级位于带隙内的哪一侧、有没有表面态钉扎、模型里是否加入外加电荷或掺杂。
图 1 中的 TiN(100)、TiN(110) 和 TiN(111) 看起来只是不同切面，功函数却会给出不同数值。原因在于电子云在表面结束的位置并不相同：有的表面外层 Ti/N 排列更对称，有的表面带有更强的层状极性，电子从固体一侧泄出到真空一侧后形成的表面偶极也不同。同一个体相材料，换一个晶面就等于换了一套 Φ 的表面模型条件。
DFT slab 里怎样读取真空能级
DFT 中计算 Φ，一般先做 slab 模型：中间是若干层表面原子，上下或一侧留真空层，周期方向仍按超胞重复。自洽计算后，体系给出 EF；再沿表面法向取平面平均静电势，找到真空区的平坦平台，把这个平台当作 Evac。若 slab 两侧表面等价，两个真空平台应接近一致；若上下表面不等价，平面平均势会带有斜率或两侧平台差，偶极修正和读取哪一侧真空就会进入结果说明。
这里的关键不是把输出文件里某个数字抄出来，而是核对 势能曲线是否真的出现平台。真空层太薄时，两个周期镜像表面仍可能相互作用；slab 太薄时，中间层还没有恢复接近体相的电子结构；非对称吸附或极性表面没有处理好时，真空区会残留人为电场。这样的 Φ 即使有三位小数，也只是在不合适模型上的数值。
图2. O:TiN(100) 和 O:TiN(111) 的 DOS 以及以 vacuum 为零点的费米能级位置，右侧图中 WF 对应 vacuum 与 EF 的能量差。DOI：10.48550/arXiv.2102.08646
图 2 右侧把 vacuum 放在 0 eV，红线和蓝线分别标出氧化前后 TiN 表面的 EF。从图上看，功函数对应的就是 vacuum 到 EF 的竖直距离；左侧 DOS 则提醒读者，EF 附近的电子态也会因表面氧化而改变。不要把 DOS 横轴的 0 eV 和真空能级 0 eV 混成同一个参照：很多 DOS 图默认把 EF 设为 0，而功函数图常把 Evac 设为 0，两个零点服务的判断对象不同。
不同程序输出的静电势可能带有不同零点，跨体系比较时要统一能量参照。同一篇文章里，如果一个图以 EF 为 0 eV，另一个图以真空能级为 0 eV，直接把峰位、带边和 Φ 数值放在同一坐标上会出错。操作上应先说明零点，再说明差值：Φ 来自 Evac 与 EF，带边来自 VBM/CBM 相对同一真空参照或同一平均势参照的位置。
表面终止、厚度和空位怎样改写Φ
功函数变化常常来自两类量的移动：一类是 Evac 被表面偶极推高或拉低，另一类是 EF 因成分、缺陷或载流子浓度变化而移动。实际材料里两者经常一起变，所以只说“功函数升高”还不够，需要进一步问：是真空平台移动了，还是费米能级移动了，还是二者都有贡献。
图 3 中，TiN(111) 的 Φ 接近 4.7 eV，TiN(100) 和 TiN(110) 分别在约 3.0 eV 和 3.2 eV 附近；同一取向从 1 nm 增厚到 4 nm 后变化很小。这个案例给出一个实用判断：厚度收敛和表面终止要分开看。如果 slab 厚度改变后 Φ 仍大幅跳动，模型可能还没有收敛；如果厚度已经稳定，不同晶面的差别就更可能来自外层原子排列和表面偶极。
判断空位对 Φ 的影响时，缺陷浓度、所在原子层和表面取向要同时写清。图 4 先给出 VN 位于不同原子层时的模型和形成能，说明缺陷放在表层、次表层或内部层，热力学代价并不相同。对金属 TiN 这类体系，单个 VN 的局部扰动可能被高电子密度屏蔽，表现为结构和 Φ 的变化都不剧烈；当缺陷浓度提高，成分从 TiN 向 TiNx 偏离，电子计数和表面电荷分布才会明显推着 Φ 走。
图 5 对应的是另一组缺陷浓度结果：单个 VN 位于不同层时，TiN(100) 与 TiN(111) 的 Φ 没有明显跃迁；多氮空位浓度升高后，TiN(100) 的 Φ 逐渐上移，TiN(111) 则有下降趋势。这里不能写成“空位一定提高功函数”或“缺陷越多 Φ 越低”。空位改变的是具体表面上的电子分布和成分比例，方向取决于晶面、空位浓度、缺陷所在层以及金属态屏蔽。
把功函数用于材料判断时要限定什么
功函数常被拿来讨论电子发射、金属/半导体接触、界面能级对齐、催化表面电子供给能力和 KPFM 表面电势差。它能回答的是“电子从这个表面到真空需要多大能量”以及“表面偶极或 EF 有没有改变”。它不能单独给出反应速率、吸附强弱、载流子迁移率或界面电子转移量；这些还需结合吸附能、PDOS、差分电荷密度、Bader 电荷、能垒、带边位置和实际电势条件。
图 6 中的 O:TiN(100) 和 O:TiN(111) 已经不是洁净 TiN 表面，外层出现 Ti-O、Ti-O-N 与 Ti-N 区域。这样的表面再去读 Φ，数值包含氧化层偶极、表面成分和 EF 改变的共同结果。若把它和洁净 TiN 的 Φ 直接放在一起比较，必须写明模型差异：是否同一晶面、是否同一厚度、真空层是否收敛、氧覆盖度是否一致、是否两侧对称、有没有偶极修正。
在论文图谱里，功函数最好和具体证据配套使用。若讨论电子逸出或场发射，低 Φ 只是起点，表面稳定、重构和吸附污染都要进入同一组核对项；若讨论金属接触半导体，Φ 可给 Schottky-Mott 估计，但界面态和化学键会改写势垒；若讨论催化，Φ 可以提示表面给电子倾向，不能替代吸附中间体自由能和反应能垒。把 Φ 写成有条件的表面能级指标，比把它当作万能活性描述符更接近材料计算结果本身。
功函数和电极电势也要分开写。Φ 使用真空能级作外部参照，电极电势通常相对 SHE、RHE 或实验设定的参比电极；在固液界面中，双电层、电荷补偿、溶剂取向和恒电位条件均会改变表面电势。把真空 slab 的 Φ 直接换算成反应电势，会漏掉界面环境和参比标尺。
