能带理论计算：从Bloch定理到平面波方法的物理基础

能带理论计算这个词，刚进计算材料领域的时候很容易被当成”在软件里点几下就出图”的操作。但真正在项目中用它来解释实验现象时，会发现自己缺的不是操作步骤，而是对底层物理逻辑的理解——为什么周期性体系要用平面波展开？赝势到底在”近似”什么？这些问题不搞清楚，参数调错了方向可能自己都发现不了。

Bloch定理：周期性体系中电子波函数的必然形式

能带理论计算的逻辑起点，是Bloch定理。一个在周期性势场中运动的电子，其波函数必然具有平面波包络乘上周期函数的形式。这不是近似假设，而是平移对称性的数学推论——只要晶格有平移不变性，单电子波函数就自然满足这个结构。

反过来想，为什么自由电子气模型里只需要平面波就够了？因为没有周期势的调制，包络中的周期函数退化为常数。而固体里的电子既不是完全自由的，也不是完全束缚的——Bloch波恰好处于中间地带。

但Bloch定理只是给出了波函数的数学形式，实际问题变成了：那个周期函数要怎么展开？这就是后面所有方法分歧的根源。

平面波基组：自然的周期性基矢，代价是陡峭区域吃不消

用平面波来展开Bloch波函数的周期部分，好处是天然的完备正交性，而且跟Fourier变换有直接数学映射——实空间的精细变化对应倒空间的高动量分量，这是做k点采样和动能截断的理论依据。

但平面波有一个致命弱点：靠近原子核的区域波函数振荡剧烈，需要极高频的平面波分量才能描述，计算量随截断能呈~EKBY³ 增长。把这个区域直接从波函数里”模糊掉”，同时对价电子做精确描述——这就是赝势的核心思路。

一个具体的项目例子：ZnO体系中，用300 eV的截断能配合超软赝势，总能收敛到~1 meV/atom，但用模守恒赝势需要500 eV以上才能达到同样的精度。差距来自赝势对应的芯区赝波函数光滑程度不同，这个选择直接影响整个计算周期。

k点采样：布里渊区的离散化如何在积分精度和计算成本之间找平衡

物理可观测量在周期性体系中是k空间整个BZ的积分，实际计算中只能用有限个离散点近似。k点网格的密度直接决定计算精度——对电子态密度敏感的物理量（如带隙、光学性质），k点网格的影响远大于截断能。

不同的k点网格生成方案（Monkhorst-Pack、Γ-centered）对应不同的对称性利用策略。Γ-centered方案在六方晶系中对带边能量的描述往往更准确，因为六方晶系的高对称k路径经过Γ点的比例更高，Monkhorst-Pack可能遗漏对称性敏感点。以六方GaN体系为例，k点从4×4×4换到6×6×6，直接带隙数值从1.72 eV收敛到1.74 eV，变化虽然只有20 meV，但对于能带对齐分析来说已经足够影响界面态的能量排布判断。

赝势与全电子方法：精度-效率的取舍不是非此即彼

赝势把内层电子和原子核的效应等效成一个平滑的有效势，换来的是平面波数量大幅减少。PAW（投影缀加平面波）方法是目前主流的折中方案——在赝势框架内通过投影算符还原全电子波函数的信息，既保住了计算效率，又拿回了靠近原子核区域的物理精度。

但PAW并不等价于全电子方法。在计算EFG（电场梯度）、超精细结构这类对近核区域敏感的物理量时，PAW与全电子LAPW方法之间仍存在~5-10%的系统偏差。这个差距在轻元素体系（如Li、B）里更明显，因为这些元素的价电子原本就分布在较靠近核的区域。

回溯到这个项目里最关键的结论：能带理论计算的真正门槛不在于软件操作，而在于理解每个近似在被研究的体系中引入的误差量级。当计算结果与实验出现偏离时，知道该往哪个方向排查，比盲调参数有效得多。

更多内容请访问 https://www.keyanxueshu.com/