ms功函数计算

我做 ms功函数计算 最初是给一个阴极材料客户算发射性能。文献给的块体功函数是 4.3 eV，我切了个 (100) 面算出来 3.9 eV，客户质疑”面怎么可能差这么多”。后来确认是表面弛豫让最外层原子内移，电子云外推变了——那次让我记住，ms功函数计算本质上算的是”表面”，不是”块体”。

一、为什么功函数对表面这么敏感

功函数 W = E_vac − E_F，是电子从费米能级逃逸到真空所需的能量。真空能级由表面外推的静电势决定，而表面原子只要一弛豫、一吸附，电势尾巴立刻变。我算过 Au(111) 和 Au(100)，差 0.3 eV 以上；同一面有氧吸附，功函数能降 1 eV（向真空方向推电子云）。所以 ms功函数计算 第一步永远是把表面模型建对，而不是急着看数。

二、核心原理：静电势平面平均

ms功函数计算 在 CASTEP 里靠 “Electrostatic potential” 任务输出三维势，再用 Analysis 做沿表面法向的平面平均。真空区势趋于常数就是 E_vac，费米能级从能级本征值取。我做ms功函数计算时，真空层至少 15 Å，否则两侧表面势互相干扰，平面平均在中段平不了。曾经真空只有 8 Å，算出的 W 震荡严重，加厚到 18 Å 才收敛。

三、关键技术要点：slab 对称性与两端

非对称 slab（只暴露一面）要固定底部几层模拟体相，否则两端都弛豫功函数无定义。我一般固定底部两层、弛豫上面 3~4 层。对称 slab（两面相同）则全弛豫但注意它是双面功函数。另一个坑是 k 点：功函数对法向不设 k（为 1），但面内 k 要够密，否则 E_vac 抖动。我面内取 0.02 Å⁻¹ 间距较稳。

四、实操流程：从表面到 W

我的步骤：Build 表面（设厚度与真空）→ 几何优化（固定底层）→ 算 Electrostatic potential → 平面平均取真空势 → 取费米能级 → W = E_vac − E_F。做吸附体系时，我先算洁净面 W 作基准，再算吸附后变化，这样相对值比绝对值更可靠。最容易被省的是底部固定，我宁可多两层也不让体相漂移。

实算一个 Cs 吸附在 W(100) 的体系：洁净面 W=4.6 eV，吸附 Cs 后掉到 2.0 eV 以下，这正是碱金属降低逸出功、用于热阴极的原理。但如果不固定底部、真空层又只有 8 Å，算出来只掉到 3.8 eV，差出近一倍，会让你误判掺杂量。所以做这类计算时，那两个”笨功夫”——厚真空、固底层——省掉任何一个，相对变化都不可信，更别提拿去解释界面的电子发射行为。

五、常见踩坑：自旋与掺杂忘了带

磁性表面做 ms功函数计算 忘了开自旋，费米能级位置和电势都错。掺杂体系也类似，表面态被改变。我算过一个 N 掺杂石墨烯，不开自旋 W 算高 0.4 eV。还有用 PBE 算出的 W 普遍比实验浅约 0.1~0.2 eV（缺乏精确交换），定性排序没问题，定量对标要谨慎或上 HSE。

六、复盘：功函数是”表面活性的尺子”

回过头看，ms功函数计算 最实用的不是报一个块体值，而是看”表面怎么改、吸附怎么调”的相对变化，它直接关联催化、发射、接触电阻。我接项目时默认先要表面/吸附场景，再决定建模精度。把真空层、底部固定、自旋三件事做对，W 才敢拿去解释界面行为。被证明有用的，是愿意为 2 Å 真空层反复加厚的执着。下次有人问”为什么你的 W 和文献差 0.2″，我会先查真空层和自旋，而不是换泛函。