钙钛矿晶体结构计算：从容忍因子到带边对齐的DFT实战

做钙钛矿材料计算的人，第一次算带隙的时候大概率会被震惊到——用PBE算出来的MAPbI₃带隙只有0.5 eV左右，实验值是1.6 eV，差了整整1 eV。

这不是你的计算设置错了，是半局域泛函（GGA/PBE）在描述钙钛矿这类强关联体系时，系统性的低估了带隙。

钙钛矿的结构自由度

ABX₃型钙钛矿的结构，比大多数人以为的复杂。A位离子（MA⁺, FA⁺, Cs⁺）可以旋转、倾斜，BX₆八面体可以扭曲、旋转，这些自由度耦合在一起，决定了最终的晶体结构和稳定性。

容忍因子（Tolerance Factor, t）是判断钙钛矿能否形成的关键参数：

t = (r_A + r_X) / [√2 (r_B + r_X)]

t在0.8-1.0之间的体系，倾向于形成立方或四方钙钛矿结构。但容忍因子只是几何判据，没考虑A位阳离子的氢键作用和取向有序，实际稳定性需要用DFT直接计算形成能。

VASP计算钙钛矿的全流程

第一步：结构建模。MAPbI₃的原始结构可以从Materials Project或者COD数据库下载。注意MAPbI₃有α相（立方，高温）、β相（四方，室温）、γ相（正交，低温）三种，建模时要选对相。

MA⁺（CH₃NH₃⁺）的取向在钙钛矿里可以自由旋转，DFT计算时需要做构象采样——把MA⁺摆成几个不同的取向，分别优化，取能量最低的构象。这一步很多人跳过，结果是结构没弛豫到位。

第二步：结构弛豫。INCAR里关键设置：

ISIF = 3   （同时弛豫晶格矢量和离子位置）
IBRION = 2  （CG算法）
EDIFF = 1E-5
EDIFFG = -0.01  （力收敛标准，eV/Å）

钙钛矿的晶格弛豫量通常有几个percent，ISIF=3不能省。

第三步：电子结构计算。用弛豫后的CONTCAR做静态计算（ISIF=2, IBRION=-1），设：

ICHARG = 11  （读入CHGCAR，不做电荷混合）
LORBIT = 11  （输出态密度投影）

交换关联泛函的选择

实际工作里的策略：结构弛豫用PBE或SCAN（便宜），电子结构和带隙用HSE06（在PBE弛豫结构上做单点计算）。

容忍因子 vs DFT形成能

容忍因子能快速筛选可能的钙钛矿组分，但有几个已知的失效案例：

• CsPbI₃：容忍因子0.81，理论上是稳定钙钛矿，但室温下会自发转变为非钙钛矿相（δ相）。DFT计算形成能才能捕捉到这个相变。
• 双钙钛矿（A₂BB’X₆）：容忍因子公式不适用，必须用DFT计算形成能判断稳定性。

DFT计算形成能的方法：E_form = E_total(ABX₃) – [E(AX) + E(BX₂)]，或者用convex hull距离（在Materials Project的Phase Diagram工具里可以做）。

带边对齐（Band Edge Alignment）

光伏应用里，钙钛矿的价带顶（VBM）和导带底（CBM）需要与电子传输层（ETL）和空穴传输层（HTL）匹配。

VASP里计算带边位置的步骤：

1. 用HSE06计算钙钛矿的能级（E_Fermi相对于真空能级）。
2. 计算材料的功函数（Work Function）：需要做slab模型，驰豫表面，计算E_vac – E_Fermi。
3. 把VBM和CBM对齐到真空能级，再和ETL/HTL的带边比较。

注意：平面平均电势（Planar-averaged potential）需要做真空层修正，消除周期性边界条件下的偶极修正。

常见坑点

坑一：用PBE的带隙指导材料设计。PBE系统性低估带隙，用它筛选出来的材料，实验上带隙可能完全不对。至少用HSE06验证。

坑二：忽略自旋轨道耦合（SOC）。铅基钙钛矿里，Pb的6p轨道有强自旋轨道耦合，打开带隙~0.1-0.3 eV。HSE06+SOC是目前最可靠的组合。

坑三：MA⁺取向没优化。MAPbI₃里MA⁺的取向直接影响带隙（~0.1 eV变化），必须用不同的MA⁺取向做弛豫，取能量最低的结构做电子结构计算。

钙钛矿计算的实际价值

DFT计算在钙钛矿研究里的核心价值，不是算出和实验完全一致的带隙（这几乎不可能），而是：

1. 筛选可能的组分（容忍因子 + 形成能计算）
2. 理解带宽变化的趋势（比如混合卤素后带隙怎么变）
3. 分析缺陷能级和离子迁移路径（这两者直接决定器件的稳定性和寿命）

这三件事，实验做起来很贵而且慢，DFT能在几天内给出定性的正确方向。