DFT计算功函数：MoS2/MoS2异质结界面工程的电荷转移分析

功函数这个概念在DFT计算里经常出现，但真正做项目的时候，很多人对它的理解停留在”真空能级减费米能级”这个公式上。一旦遇到异质结界面、表面偶极修正、电荷重分布这些实际因素，功函数的计算就远不是一个简单差值能解决的问题。

项目背景：异质结能带对齐的决定性参数

团队去年处理了一个二维材料异质结项目，研究对象是MoS2/WSe2范德华异质结。客户是做光电探测器的，需要知道两层材料堆叠之后，界面处的能带排列方式——是Type-I还是Type-II对齐，功函数差异能驱动多大的内建电场。这个信息直接决定器件的工作机制和量子效率。

在开始建模之前，团队先分别计算了单层MoS2和单层WSe2的功函数。MoS2的计算结果为5.12 eV（PBE泛函，500 eV截断能，12×12×1 k点），WSe2为4.57 eV。两个值之间的差异0.55 eV，看起来不大，但放在二维材料的尺度上，这个功函数差对应的内建电场足以在界面区域产生可观的电荷重分布。

真空层的选取：功函数计算中最容易被忽视的参数

DFT计算功函数需要一个真空层来定义真空能级。对于二维材料，真空层厚度对功函数的影响比体相材料大得多——因为表面偶极的长程尾部衰减较慢，真空层不够厚时，两侧表面的偶极会通过真空层相互作用，导致真空能级被低估。

团队做了真空层厚度从12 Å到25 Å的收敛测试。结果显示，12 Å时功函数偏低约0.08 eV，15 Å时偏差缩小到0.03 eV，18 Å以上基本收敛。考虑到单层MoS2本身的厚度约6.1 Å，18 Å的真空层意味着周期性单元的总高度约为30 Å。这个值看起来浪费，但功函数对真空层厚度的敏感度确实需要这个余量。

在异质结计算中，真空层的选取更棘手——两层材料的组合高度更大，需要的真空层也要相应增加。团队最终在异质结模型中使用了20 Å的真空层，周期高度约48 Å，k点密度降至9×9×1以控制计算量。

界面偶极与功函数偏移

单层功函数给出了材料的本征值，但异质结形成之后，功函数会发生变化——这个变化来源于界面处的电荷转移。MoS2的电子亲和力比WSe2强，所以电子会从WSe2向MoS2方向转移，在界面区域形成一个内建偶极矩。

根据团队的DFT计算，异质结形成后MoS2侧的功函数从5.12 eV下降到4.89 eV，WSe2侧从4.57 eV上升到4.71 eV。功函数差的收敛过程从0.55 eV缩小到0.18 eV——这个”功函数对齐”是界面电荷重分布的自然结果，也是异质结能带对齐的物理根源。

为了量化界面电荷转移，团队用平面平均的差分电荷密度做了分析。电荷转移主要集中在界面两侧约5 Å的范围内，积分结果为每单位面积0.018 e/Ų。这个量级在范德华异质结中属于中等偏强的水平，说明两层之间的耦合不仅仅是物理吸附，还包含了显著的电子相互作用。

计算精度与实际器件的对应

需要说明的是，PBE泛函的带隙低估问题会影响功函数的绝对值。MoS2的PBE带隙约1.67 eV，实验单层值约1.88 eV；WSe2的PBE带隙约1.56 eV，实验约1.65 eV。这些带隙差异会传导到功函数的计算中，使得PBE给出的功函数系统性偏低约0.1-0.15 eV。

如果项目需要功函数的绝对值精度在0.1 eV以内，建议使用HSE06杂化泛函——但计算量会增加到PBE的10倍以上。在这个项目中，团队选择PBE的原因是：客户关注的核心是功函数差（即内建电场强度），而PBE对功函数差的描述与HSE06的差异通常在0.03 eV以内，远小于绝对值本身的系统误差。

最终交付的分析报告中，功函数的计算结果为客户判断异质结类型（Type-II对齐）提供了定量依据。结合有效质量和载流子迁移率的计算，团队还评估了该异质结作为光电探测器活性层的光响应范围和响应速度，形成了完整的器件性能预测框架。功函数看起来只是一个数，但在这个项目中，它连接了材料的电子结构、界面的物理化学、以及器件的工作原理。

如需了解更多的科研计算案例与项目经验，可以参考 https://www.keyanxueshu.com/ 上的技术文章。