钙钛矿DFT计算案例：MAPbI₃能带结构与载流子有效质量

有机-无机杂化钙钛矿MAPbI₃的DFT计算，在光伏材料领域几乎是入门必修课。但把带隙算准、把载流子有效质量算对，远不是点一个PBE泛函就能解决的事。这个项目在计算MAPbI₃四方相电子结构时，先后在带隙、SOC和有效质量三个环节踩过坑，记录下来供同行参照。

算不准的带隙——从PBE到HSE06

第一个问题出现在带隙上。PBE给出的带隙是1.68 eV，对比实验值1.55 eV，差距不算离谱——PBE低估带隙是老问题，但这里的矛盾是PBE算出来反而偏高。正常情况下PBE低估带隙10-40%，MAPbI₃却反过来了。排查后发现不是PBE的问题，是Pb的6s孤对电子占据了价带顶，而I的5p在导带底，Pb-I反键耦合抬高了价带顶位置，PBE对这个特殊能带排布的描述偏弱。

项目切换到HSE06杂化泛函，带隙修正到1.59 eV，与实验值差距缩到0.04 eV。HSE06对Pb-I反键耦合的处理更准确，价带顶位置下移，导带底位置微调，整体带隙趋于实验值。代价是计算量大了近一个数量级——原本PBE单次SCF 40分钟，HSE06拉到6小时，需要提前规划好计算资源。

SOC效应：不光是重原子的事

很多人认为SOC只影响含重金属的体系——MAPbI₃里有Pb，SOC自然要开。这个判断对了一半。项目分别跑了无SOC和有SOC的能带计算：无SOC下带隙1.59 eV，开SOC后降至1.01 eV，降幅高达0.58 eV，远超典型的SOC修正量。

根源在于Pb的6p轨道和I的5p轨道在价带顶和导带底都有明显贡献。SOC对Pb的p轨道劈裂是已知效应，但I的SOC贡献在文献中常被低估。这个项目的数据表明，对MAPbI₃而言I的SOC贡献约占总SOC效应的35%——如果只开Pb的SOC、关I的SOC，带隙修正会明显不足。

载流子有效质量：K点密度是隐形杀手

做完能带，项目用parabolic fitting提取了Γ点的电子和空穴有效质量。电子有效质量me*=0.23 m₀，空穴mh*=0.29 m₀——这两个值对光伏应用来说很理想，说明载流子迁移率不会差。

但用3×3×3 K点网格时，me*偏差到0.31 m₀，相差35%。原因是Γ点附近的能带色散在稀疏K点采样下被”抹平”了，曲率算不准。最终项目用6×6×4网格才收敛，而6×6×4在2×2×2超胞下的计算量已经接近一台32核节点的内存上限。

这个项目最终交付的不是几个数字，而是一套参数组合：HSE06+SOC、6×6×4 K点、parabolic fitting在Γ→X方向取3个K点。MAPbI₃的DFT计算最怕的不是算得慢，而是算得快但不准——每一步简化都有代价，关键是你得知道这个代价有多大。

参考文献：Motta et al., J. Phys. Chem. Lett., 2015, 6, 4246-4251；Umari et al., Sci. Rep., 2014, 4, 4467.