ms功函数计算

结合案例图分析，本文以这个金属电极选型项目为例，讲功函数对表面有多敏感、模型设置有哪些雷。这是为一家电子器件客户做的金属电极功函数评估项目。客户在选栅极金属，关心不同晶面发射特性和接触势垒，需要一组可信的功函数数据。第一张案例图是某金属表面沿法线方向 z 的静电势分布，真空区电势平稳对应真空能级 E_vac，材料内部因电子气振荡在 2–4 eV 波动，费米能级 E_F 与 E_vac 的差就是功函数 Φ=3.10 eV。第二张案例图是 Cu、Al、Ag 不同晶面的功函数对比，Cu(111) 高达 4.98 eV，明显高于 Cu(110) 的 4.48 eV。

项目过程：金属电极选型项目的执行链路

方法上我选了 Materials Studio 的 CASTEP 做平面波 DFT，而不是 VASP，原因是客户要的是多晶面对比且课题组熟悉 MS 的 Workfunction 模块。但核心难点不在引擎，在 slab 模型本身。我的设置：真空层 20 Å、slab 7 层且表层弛豫、k 点 8×8×1。真空层 20 Å 是我收敛测试定的——10 Å 时 E_vac 在真空区还明显倾斜，15 Å 基本平但边缘仍有 0.03 eV 波动，加到 20 Å 后 E_vac 平台起伏 < 0.01 eV，再厚纯浪费机时。 客户最初只要”某个晶面的功函数一个数”，我主动把 Cu 的 (111)/(110) 和 Al、Ag 都算全了，因为栅极接触势垒对晶面取向极敏感，单报一个会误导器件设计。

案例分析：从静电势到分晶面功函数

图里读出的数要能落地。第一步读平面平均静电势的真空平台值 E_vac，再减去费米能级 E_F 得 Φ——但平面平均内部有 Friedel 振荡，必须用 macroscopic average 抹平短周期振荡才能稳取 E_vac。我对真空区取多个 z 切片算 E_vac 标准差，确认 < 0.01 eV 才认为平台可靠；否则功函数会随切片位置跳 0.05 eV。 第二步把 Cu(111)/Cu(110)/Al/Ag 的 Φ 列成表，Cu(111) 4.98±0.05 eV 最高，Cu(110) 4.48±0.05 eV 最低，各向异性达 0.5 eV。误差量级上，PBE 和 LDA 对同一 slab 的功函数差通常只有 0.1 eV，说明误差九成来自 slab 设置而非泛函选择。

总结教训：真空层与偶极修正的坑

做的过程中我踩过实打实的坑。第一坑是真空层：10 Å 时 Φ 比收敛值低 0.18 eV 且 E_vac 倾斜，15 Å 偏差收窄到 0.04 eV，20 Å 后稳定在 3.10±0.02 eV——这条收敛曲线让客户直观看到”为什么必须 20 Å”。 第二坑是不对称 slab 的偶极修正：忘开 dipole correction 时，真空两侧电势平台差了 0.25 eV，功函数直接被系统偶极矩带偏；加上后两侧平台对齐，偏差降到 0.02 eV 以内。第三坑是 slab 层数，5 层时表层弛豫未充分传播、功函数比 7 层偏低 0.06 eV，故取 7 层。

项目经验：功函数交付的可复用清单

作为可复用的项目经验，第一，真空层、slab 层数、表面弛豫、偶极修正四项必须全查并公示，缺一不可。第二，给客户一定要分晶面给数据并标注清洁表面，不同晶面差 0.5 eV，只报平均值他会设计出错。第三，涉及异质结还要和能带对齐一起给，误差带 ±0.05 eV 必须标注。客户最终选了高功函数的 Cu(111) 做栅极接触层。说到底，功函数项目交付的不是一个数，而是”分晶面 + 清洁标注 + 能带对齐”的一组可比参数，任何一项缺失都会让客户的器件设计踩坑。

ms功函数计算

计算工程师团队：500+ ｜ 成功计算案例：5W+ ｜ 1W+ 助力顶刊发布经验

需要模拟代算服务可以联系我们，可免费评估项目的计算费用和免费出计算方案。