MS计算费米能级：从能带结构到功函数的完整工作流

费米能级在Materials Studio里的计算流程看起来很短——优化结构、算单点能、看能带图就完事了。但真正在学术写作中引用费米能级数值的时候，反复验证的次数往往比表面上看起来多得多。

费米能级的第一个变量来自结构优化。表面体系的费米能级对表层原子的弛豫极度敏感。做过一个TiO₂（110）面的项目，表面第一、二层原子放开弛豫后，费米能级相对于真空的位移超过了0.2eV。不弛豫表层的计算结果在这个项目里完全不具备参考价值。经验是至少弛豫表面三层原子，确认层间距的变化趋于收敛，再锁定结构去做能带。

做表面弛豫时必须注意底层原子的固定策略。MS默认只弛豫原子坐标而不优化超胞的晶格参数——这对于表面slab模型是正确的，因为面内晶格常数应该保持与体相计算结果一致。但如果弛豫表面时没有固定底层几排原子，底层原子会跟着上层一起漂移，整个slab发生整体性位移，导致费米能级参照的真空能级也跟着偏移。正确做法是弛豫上面三层，底部两到三层原子的坐标固定。

能带计算阶段的K点路径怎么选，直接影响费米能级附近的色散关系的采样质量。CASTEP里默认的高对称点路径对常规晶体够用，但涉及到低对称性表面和异质结构时，必须手动编辑K点路径文件。一个实用的经验：沿表面法线方向的K点采样取Γ点就行，面上取至少两个方向的高对称点路径。如果沿某个方向价带顶和导带底各自的能级没有在K空间里被完整捕捉到，你得到的费米能级位置就是一个在Brillouin区里做了不完全积分的结果。

费米能级数值本身不足以描述电子性质。真正有用的信息是费米能级与价带顶、导带底的相对位置。一个常见的陷阱是：CASTEP输出的能带数据默认把费米能级归零，也就是直接把价带最高点对齐到零。这个做法对于带隙本身的报告是方便的，但对于功函数计算就是致命的。解决方法很简单：在单点能计算时勾选输出静电势，在.castep文件里找到真空能级的绝对数值，再把费米能级转回绝对坐标。

如果把费米能级加真空能级得到功函数是一个合理的目标，那这个计算在DFT框架内通常会比实验值偏低0.3~0.5eV。这和泛函的带隙低估不同——功函数的系统偏差主要来自表面偶极矩在DFT中的不完全描述。对自己的计算结果要有一个清醒的认知：数值上的绝对偏差是方法本身带来的，关键在于能正确描述不同表面终端的相对差异。
掺杂对费米能级的影响是一个在实际项目中经常遇到的变数。以一个n型掺杂的ZnO体系为例，Al掺杂浓度从1%提到3%，费米能级从导带底下方0.1eV的位置一路推进导带——从非简并半导体变成了简并半导体。这个变化不只是数值上的位移，而是材料导电机制的根本转变。在做掺杂体系的费米能级计算时，超胞里掺杂原子的实际浓度和空间分布要尽量接近实验条件——一个2×2×2的ZnO超胞（含32个Zn位）里放一个Al原子等效掺杂浓度约为3%，这个浓度下得到的费米能级位移和实验对比才有物理意义。更低的掺杂浓度需要更大的超胞，但计算量不是线性增长——超胞原子数翻倍，CASTEP的单点能计算时间可能翻三到四倍。做掺杂计算前先想清楚：你关心的掺杂浓度范围内，用多大规模的超胞去模拟，计算结果在统计上才有代表性。

https://www.keyanxueshu.com/category/dft/