钙钛矿光伏材料能带结构计算中杂化泛函的选择

接到MAPbI₃的电子结构计算需求时，需求方给的第一个要求就是”带隙要准，后面还要接着算光吸收谱”。用标准PBE跑了一轮，带隙1.04eV。而同一批样品在实验室用UV-Vis吸收边测量的带隙是1.55eV——差了0.5eV，将近50%的低估。

这在钙钛矿体系里不是意外。PBE对带隙的系统性低估是DFT本身的问题——Kohn-Sham方程里交换关联势的导数不连续性没被考虑，所有半导体的带隙都会被压缩。但钙钛矿的情况比传统半导体（如Si、GaAs）更严重，因为Pb的6s孤对电子和I的5p电子之间还有额外的自相互作用误差，进一步压窄了带隙。接下来我们需要看两个主流杂化泛函能把这0.5eV的差距拉回来多少，以及在不同精度需求下怎么选。

GGA的缺陷在钙钛矿体系中被放大了

PBE带隙偏低不是新鲜事，但MAPbI₃体系的偏低幅度（~0.5eV）比传统半导体大得多。Si的PBE带隙约0.65eV，实验值1.17eV，偏低约45%。GaAs的PBE带隙约0.2eV，实验值1.42eV，偏低近85%。数字上MAPbI₃的45%似乎还好，但绝对值上0.5eV的差已经足以把1.55eV的直接带隙材料误判为近乎半金属——后续算出来的光吸收边会大幅红移。

问题的根源在Pb-6s和I-5p的杂化态。这两组轨道在价带顶和导带底都有贡献，而且空间分布弥散，自相互作用误差比局域的d电子还难搞定。GGA+U可以部分修正d电子的自相互作用，但对s-p杂化几乎没有作用——U值作用在局域的d或f轨道上，s和p根本套不进Hubbard模型的框架。

结论很清楚：钙钛矿的带隙问题，GGA+U搞不定，只能上杂化泛函。

HSE06 vs PBE0：两个方案的实测对比

两个体系——MAPbI₃（四方相，I4/mcm）和CsPbBr₃（立方相，Pm3̄m）——各用三种泛函跑了能带计算。截断能统一用520eV（Pb的PAW赝势ENMAX约322eV，520eV覆盖），k点在Γ-centered网格上取4×4×4。SOC（自旋轨道耦合）在所有计算中都打开——Pb的SOC效应在钙钛矿里不可忽略，能带会进一步被SOC压窄0.2-0.3eV。

两组数据摆在一起能看出几个规律。

HSE06对MAPbI₃的带隙修复几乎是完美的，1.48eV离实验值只差0.07eV。这个精度在DFT框架内已经是顶格表现了——剩下的0.07eV可能是有限温度晶格振动导致的带隙重整化效应，这是电子结构计算本身算不出来的。

但HSE06到CsPbBr₃上开始出现0.06eV的低估。虽然差距仍然远小于PBE，但趋势值得注意：随着带隙变大，HSE06的低估有增大的迹象。这与HSE06的屏蔽参数μ=0.2Å⁻¹有关——这个值在中等带隙体系里是经验最优的，但到了宽带隙材料，屏蔽长度可能不够匹配。

PBE0在两个体系里都偏高：MAPbI₃高0.14eV，CsPbBr₃高0.16eV。PBE0混合25%的精确交换，没有短程/长程拆分，交换部分直接全局替换。对于Pb基钙钛矿这种电子屏蔽相对较强的体系，25%的精确交换偏多了。

方案选择取决于后续计算任务

两个方案的选择不是”哪个更准”，是”你的下一步要干什么”。

如果下一步算光吸收谱和激子结合能，带隙精度优先——HSE06。PBE0的0.14eV正向偏差虽然不大，但到了GW-BSE那一步，起始点的带隙偏高会让激子结合能被高估。

如果下一步算缺陷形成能和掺杂，需要更准确的带边位置而非只关心带隙绝对值——PBE0对带边的修正更均匀，不像HSE06对导带的修正偏强。而且掺杂体系涉及电荷局域化，精确交换比例高一点对缺陷态的定域性描述有帮助。

如果只是快速筛选材料，时间成本是瓶颈——HSE06。几十个原子的小晶胞用HSE06大概2小时，PBE0要3.5小时左右（VASP 6.3，32核）。批次计算里，每个结构多1.5小时，20个结构就多出一天。

在这个项目里，客户后续要做激子效应计算，我们选了HSE06。如果他要做的是缺陷工程筛选，我们会推荐PBE0。杂化泛函没有”默认最优选项”——你的下一步决定你的起跑线。