钙钛矿晶体结构计算：Goldschmidt容差因子与DFT结构优化的协同验证

钙钛矿晶体结构计算在光伏和催化领域玩的是”结构决定功能”的逻辑——ABX₃框架里A位阳离子的大小决定八面体网络的倾斜程度，B位过渡金属决定电子结构，X位阴离子（通常是O、F、Cl）调控化学键极性。稳定性判断是第一步，也是很多人直接套Goldschmidt容差因子草草了事的那一步。问题的根源在于：容差因子只能做粗筛，假阳性和假阴性都不少。

容差因子的粗筛边界

Goldschmidt容差因子t = (r_A + r_X) / [√2(r_B + r_X)]，理论上t在0.9-1.0之间形成理想立方钙钛矿，0.71-0.9之间形成畸变钙钛矿（正交、菱方），低于0.71就偏离钙钛矿结构了。但经过大量文献调研和团队内部验证，这个数字规则漏洞很大。

以FASnI₃（甲脒锡碘）为例，A位FA⁺离子有效半径约2.53 Å，Sn²⁺半径1.18 Å（考虑6配位），I⁻半径2.20 Å——算出来t≈0.97，落在理想区间。按理说应该是稳定的立方结构，但实际上纯FASnI₃在室温下发生复杂相变，从高温α相（立方）→β相（四方）→γ相（正交），晶体对称性逐步降低。容差因子给出的稳定性判断过于乐观——它只考虑了离子半径的几何匹配，完全没有纳入A位分子动力学和八面体倾转的能量偏好。

DFT结构优化的四步收敛策略

钙钛矿晶体结构计算在DFT层面需要严格执行四步收敛，跳过任何一步都可能让结构落在错误的局部极小值上。团队处理过CsPbI₃（铯铅碘）体系，做的第一件事就是晶格常数的全面优化——ISIF=3，允许原子位置和晶格形状同时弛豫。

PBE泛函给出的CsPbI₃立方相优化晶格常数是6.39 Å，实验值6.29 Å——高估了约1.6%，这是PBE对钙钛矿的典型行为。修正方案是在优化中加范德华修正（optB88-vdW），晶格常数降到6.32 Å，偏差缩小到0.5%。步二是在优化后的晶格上做离子弛豫——固定晶格，放开原子位置（ISIF=2），关注Pb-I键长的分布是否均匀。步三是用HSE06杂化泛函在PBE优化好的结构上做单点计算，修正带隙。步四是声子谱验证动力学稳定性。

这里的关键判断发生在步三：PBE给出的带隙是1.44 eV，HSE06给出的带隙是1.73 eV，实验光学带隙约1.73 eV。如果只看PBE结果，带隙被低估近0.3 eV，对光伏应用的预测会产生偏差——用PBE带隙做能带对齐，给体和受体之间的LUMO-HOMO匹配会偏离真实值。

八面体倾转的对称性搜索

容差因子判据失效的核心原因之一，是它不捕获八面体倾转（octahedral tilting）。钙钛矿A位阳离子偏小时，B-X八面体会发生集体旋转来优化A-X配位。GdFeO₃型正交畸变（a⁻a⁻c⁺倾转模式，按Glazer记法）在t<0.9的钙钛矿中非常普遍。

团队在CsPbI₃的不同相之间做过比较：立方相（无倾转）→四方相（a⁰a⁰c⁻倾转）→正交相（a⁻a⁻c⁺倾转），总能依次降低约0.05 eV/f.u.和0.08 eV/f.u.——这个能量差虽然不大，但在室温下足以驱动相变。DFT结构优化如果不从低对称性结构开始，从立方相直接弛豫往往会”卡在”立方相的鞍点上，弛豫后的力收敛了，但全局结构不是能量最低点。

流程上的教训很明确：做钙钛矿晶体结构计算，不能用一个高对称性初始结构直接ISIF=3。应该先生成可能的Glazer倾转模式（a⁰a⁰a⁰、a⁰a⁰c⁻、a⁻a⁻c⁺等），每个模式跑一轮离子弛豫，选能量最低的结构做ISIF=3全局优化。

声子谱验证

结构优化收敛了，带隙也和实验对上了，最后一把钥匙是声子谱。Phonopy做有限位移法计算，2×2×2超胞（160个原子），每个原子偏移±0.01 Å。CsPbI₃正交相的计算结果：全正频率，最低声子频率约0.3 THz——动力学稳定。但如果只在PBE优化过的结构上算声子，没有做vdW修正，在Γ点附近会出现虚频（~0.2i THz），对应八面体旋转的不稳定模式。

这个虚频不是物理上的，而是PBE对晶格参数的高估（>1%）造成的。晶格偏大→八面体间空隙偏大→倾转回复力偏弱→声子模软化→虚频出现。修正晶格常数之后虚频消失——这是一个结构参数的连锁影响链条，从晶格参数到声子稳定性，偏差在每一步放大。

钙钛矿晶体结构计算是一套多工具协同的工作：容差因子粗筛→对称性搜索→DFT分级优化→声子验证，四步缺一不可。更多具体参数和运算策略，在站上同类项目复盘中有详细展开。

更多内容请访问 https://www.keyanxueshu.com/