钙钛矿太阳能电池界面态计算案例 | 第一性原理计算实战

钙钛矿/TiO₂界面态的第一性原理计算

这个项目一开始并没有显得特别棘手。

体系是MAPbI₃/TiO₂的异质结界面，目标是计算界面态密度，分析陷阱态的能级分布。看起来是一个标准的DFT任务——建模、设参数、跑、出图。项目在第一步卡了将近两周。

不是卡在参数上，而是卡在界面模型本身。

MAPbI₃是典型的极性材料，TiO₂是过渡金属氧化物，两者的晶格常数相差约8.3%。界面超晶胞建起来之后，晶格失配产生的赝应力会直接污染电子结构计算结果，让界面态”凭空出现”或者”异常消失”，数据没有任何可信度可言。前期版本的计算结果在这里大幅抖动——两次收敛到完全不同的能级位置，甚至连定性结论都相互矛盾。

解决这个问题的方向，最终落在界面应变分配策略上：让两侧材料各承担一部分晶格失配，而不是强制让某一侧完全压缩或拉伸。这个方案本身不复杂，但执行起来需要对两种材料的弹性常数都有足够的了解，否则应变分配的比例缺乏物理依据，效果反而适得其反。

建模稳定之后，INCAR参数的设置是第二道关卡。

这个体系用了HSE06杂化泛函。原因直接：TiO₂的GGA带隙严重低估（约1.8 eV vs 实验值3.0 eV），这个偏差会让界面的能带对齐完全失真。HSE06可以把带隙修正到接近实验值，代价是计算量约为GGA的8-10倍。

INCAR关键设置如下：

“` ALGO = Damped TIME = 0.4 HFSCREEN = 0.2 AEXX = 0.25 ENCUT = 520 EDIFF = 1E-5 NSW = 100 IBRION = 2 ISIF = 2 ISMEAR = 0 SIGMA = 0.05 “`

ALGO选Damped而非All，原因是杂化泛函下All算法在这个体系里出现了振荡不收敛的问题，切换到Damped之后稳定了。这个经验不在官方文档里，但在实际操作中出现的频率不低。

第三个麻烦出现在界面真空层的厚度上。

初始版本的真空层设置了12 Å，结果界面两侧的静电势出现了周期性噪声，偶极修正没有生效。排查之后发现是超晶胞的c轴方向尺寸不足，导致DIPOL标签指向了错误的位置。调整之后，界面处的静电势平滑了，能带对齐的结果也稳定了下来。

最终计算结果显示，MAPbI₃/TiO₂界面在价带顶位置存在两个明显的界面态，能量分别在VBM+0.3 eV和VBM+0.7 eV处。这两个位置与实验上观测到的快速复合行为在能级上吻合，为后续界面钝化方案的设计提供了定量依据。

这个项目的整个计算周期是27天。其中建模和参数调试占了将近60%的时间，真正的计算时间反而不是瓶颈。

回过头看，界面体系的第一性原理计算和单一材料的计算在难度上不是同一量级。建模阶段的每一个决策——失配处理方式、超晶胞尺寸、真空层厚度、偶极修正——都会向下传递影响，最终体现在数据的可信度上。这些环节的每一个如果处理不到位，后面跑出来的数字，其实没有太大意义。

如需了解类似异质结界面的第一性原理计算服务，欢迎访问科研学术网（https://www.keyanxueshu.com）查看更多案例。