钙钛矿dft计算

钙钛矿材料（通式 ABX₃）这两年在光伏和催化两头都很热，但我这个项目我要回答的不是”它能不能发电”，而是一个更基础也更烧钱的问题：某体系钙钛矿在合成和服役过程中到底稳不稳、带隙合不合适、缺陷会不会成为非辐射复合中心。客户拿了一组卤化物钙钛矿配方来做钙钛矿dft计算，想用计算先把实验配方筛一遍。

第一步是结构搭建与弛豫。ABX₃ 里 A 位是有机/无机阳离子、B 位是金属（Pb、Sn）、X 是卤素。我先做容差因子 t = (r_A + r_X) / √2(r_B + r_X) 的几何判据初筛，t 落在 0.8–1.0 之间才算结构稳定。然后做 full relaxation 拿到平衡晶格常数，和实验值对比误差控制在 2% 以内才算模型可靠。这里有个坑：有机 A 位阳离子（比如 MA、FA）柔性很大，初始构象设错会让弛豫跑偏，我用分子力学先预优化阳离子再塞进晶格。

第二步是稳定性判据。我算了形成能 ΔH_f 和分解焓，负值越大越稳定。更进一步，算出了容忍因子和八面体因子的关联图，发现卤素混合（I/Br 共掺）能在保持带隙可调的同时显著提升相稳定温度。针对热稳定性，我还算了 Goldschmidt 容差因子随温度的变化，部分配方在 330 K 以上容差因子越界，提示会有相分离风险——这点和客户后来实测的 85℃ 衰减曲线能对上。

第三步是电子结构。带隙是我这个项目的关键输出。标准 PBE 泛函对带隙系统性低估 30%–50%，所以我用 HSE06 杂化泛函做了修正，杂化比取 0.25，算出来带隙和实验光吸收边差不到 0.1 eV。我还画了能带和态密度，确认导带底由 Pb 的 6p 主导、价带顶由 X 的 p 轨道主导，是典型的直接带隙特征，这对光伏是有利的。缺陷形成能是另一个重点：我算了几种本征点缺陷（空位、间隙、反位）在不同费米能级下的形成能，发现碘空位 V_I 形成能最低，是主要的 p 型浅缺陷，会引入非辐射复合通道——这解释了客户样品的低荧光量子产率。

踩坑方面，第一是杂化泛函的算时。HSE 比 PBE 慢一个数量级，全体系跑 HSE 不现实，我的做法是 PBE 跑结构、只在平衡结构上做单点 HSE 修正带隙，兼顾精度和成本。第二是缺陷计算的自洽。缺陷形成能和费米能级强相关，必须用裕度修正（potential alignment）处理带电缺陷的静电修正，否则形成能会漂移几 eV。第三是 spin-orbit coupling（SOC），含 Pb 的重元素体系不考虑 SOC，带隙会系统性偏低 0.1–0.2 eV，我最后都加上了。

从我的工程经验看，钙钛矿dft计算 的价值不在于”算出一个漂亮带隙”，而在于把稳定性、带隙、缺陷活性放在一张图里权衡。很多配方带隙完美但缺陷形成能太低，一做就漏电流；有些很稳定但带隙太大吸不了光。我给客户的建议是优先选 I/Br 混合、控制 A 位尺寸在容差窗口中段的配方。更多钙钛矿计算的实战细节，见 [钙钛矿dft计算](https://www.keyanxueshu.com/category/dft/)；需要搭建这类计算流程，[钙钛矿dft计算](https://www.keyanxueshu.com/) 上也有服务说明。