有机-无机钙钛矿表面建模的难点，不在于把 CIF 切成 slab，而在于同一化学式下存在晶相、A 位阳离子取向、PbX₆ 八面体倾斜、表面终止和偶极方向等多个自由度。钙钛矿表面模型首先是结构身份问题，后续吸附能、缺陷态和带边对齐都建立在这个身份之上。
本文以 MAPbI₃、MAPbBr₃ 和 FAPbBr₃ 这类杂化钙钛矿为对象，讨论 VASP 建模时怎样确定晶体母体、怎样生成表面终止、怎样处理 MA⁺/FA⁺ 取向，以及怎样用 AIMD 和电子结构检查优化后的模型。正文不讨论小分子吸附，重点放在表面模型本身。
晶相和取向怎样选？
钙钛矿 APbX₃ 的 A 位有机阳离子、PbX₆ 八面体和卤素骨架共同决定晶体对称性。室温下常见四方或近似立方结构，低温可能出现更低对称晶相。建模前需要确认目标实验或论文讨论的是 bulk、薄膜取向还是特定表面；图1 展示 MAPbBr₃ 与 FAPbBr₃ 薄膜取向差异，说明晶面偏好会直接影响应选择 (100)、(110) 还是其他表面。
图1 MAPbBr₃ 与 FAPbBr₃ 薄膜取向的 GIWAXS 结果。该图用于说明晶面取向会影响表面模型选择。 DOI: 10.1021/acs.chemmater.3c00075。
如果研究薄膜表面，应优先让 slab 晶面和实验取向对应；如果研究本征材料性质，则先优化 bulk 晶胞，再从低指数面建立 slab。晶面选择应来自物理问题，而不是来自切面工具默认值。同一篇计算中，不宜把 MAPbI₃ 的四方相、MAPbBr₃ 的立方相和 FAPbBr₃ 的不同取向直接放在同一个表面能排序里。
bulk 优化后先核对晶格常数、Pb-X-Pb 键角、A 位取向和总偶极方向。PBE、PBEsol、SCAN 或含 vdW 修正的泛函会给出不同晶格和氢键强度，杂化泛函可用于电子结构校验，但不一定适合作为所有大 slab 的几何优化起点。计算软件可使用 VASP，结构准备可由 VESTA、ASE 或 pymatgen 辅助完成。
A 位阳离子能简化吗？
MA⁺ 或 FA⁺ 不能只当作背景电荷。图3 给出 MAPbI₃ 中 MA⁺ 的八种取向定义，说明 C-N 轴方向、倾斜角和周围 PbI₆ 八面体之间存在明确几何关系。若把有机阳离子随意旋转，优化可能收敛到局部构型，表面偶极和氢键网络也会改变。
图2 MAPbI₃ 中 MA⁺ 有机阳离子的八种取向定义。该图用于说明有机阳离子不是可以任意摆放的填充物。 DOI: 10.1038/srep21687。
常用简化是先选取几组代表取向，例如面内、面外和交替偶极排列，并在同一超胞下优化比较总能。A 位取向应作为模型变量显式比较，而不是藏在 POSCAR 初始坐标里。对于表面 slab，顶层和底层的 MA⁺ 取向还会改变两侧真空平台，必要时应做对称 slab 或启用偶极修正。
如果计算目标是表面吸附或缺陷，A 位取向还会影响局域氢键和带边态。以 MAPbI₃ 表面为例，表层 MA⁺ 的 N-H 指向卤素空位附近时，局域静电势可能不同于 N-H 背向空位的模型；若这两种结构能量接近，应分别继续做吸附或缺陷计算。
PbX₂ 和 AX 终止怎样比较？
杂化卤化物钙钛矿的同一晶面可能有PbX₂ 终止，也可能有 AX 终止。图2 给出 MAPbBr₃ 与 FAPbBr₃ 代表性 (100)/(110) 表面的 DFT 表面能比较，说明终止面不是形式标签，而是会改变表面稳定性和后续电子结构。表面终止面必须通过表面能或目标环境来选择。
图3 MAPbBr₃ 与 FAPbBr₃ 代表性 (100)/(110) 表面的 DFT 表面能和终止结构。该图用于说明 slab 终止面必须单独比较。 DOI: 10.1021/acs.chemmater.3c00075。
表面能可写作 γ = (E(slab) – nE(bulk) – ΣN_i μ_i)/A。若 slab 两面相同，分母可按 2A；若只有一面为研究表面，另一面经钝化或不同终止处理，表面能公式要和该 slab 的几何身份对应。对非化学计量终止，Pb、卤素和有机盐化学势会决定稳定性窗口。
γ = [E(slab) – nE(bulk) – ΣN_i μ_i] / A 式(1)
slab 厚度通常从 4-6 个无机层起步测试，真空层从 15-20 Å 起步，并检查 LOCPOT 平面平均静电势是否有平台。含有机阳离子的 slab 常有内建偶极；若上下表面不等价，`LDIPOL = .TRUE.` 和正确的 `IDIPOL` 方向应用于静态检查。
固定层怎样设置？
底层固定可以减少表面整体漂移，但过多固定会把 PbX₆ 八面体重构压住。通常可固定底部 1-2 层无机骨架，释放表面层和附近 A 位阳离子；对于对称 slab，可以释放两侧表面并保持中心层接近 bulk。固定层设置应跟研究对象对应：若研究表面重构，表层必须释放；若只研究带边对齐，中心层要保持 bulk-like 区域。
优化后要检查 Pb-X 键长、Pb-X-Pb 角、表层 A 位阳离子是否伸入真空、卤素是否异常脱离骨架，以及两侧真空平台是否可读。若表面在优化中发生大幅重构，应把它作为新的重构表面命名，不再沿用初始终止名称。
对含缺陷表面，还要把缺陷放在不同层位测试。表层卤素空位、次表层空位和体相空位对应不同物理问题；若只选择一个缺陷层位，会把表面终止效应和缺陷层位效应混在一起。
氢键模式怎样影响能量？
图4 展示 MA⁺ 与 PbI₆ 八面体之间两类氢键作用模式。N-H···I 的方向会影响局部八面体倾斜和 A 位分子稳定性，也会改变表面附近的偶极。有机阳离子取向通过氢键和静电势影响表面结构，不是只改变图片外观。
图4 MA⁺ 与 PbI₆ 八面体之间两类氢键作用模式。该图用于说明有机阳离子取向会改变局域结构稳定性。 DOI: 10.1038/srep21687。
实际计算中，可以先在 bulk 中比较几种低能 MA⁺/FA⁺ 排列，再把它们切成相同表面。也可以先构建表面终止，再在表层手动设置几种代表取向。两条路线都可以，但同一组比较中应保持晶胞、层数、真空和固定层一致。
若只做一个取向，最好选择与实验取向或低能 bulk 排列相符的结构，并在方法部分说明。若研究表面电势、功函数或带边对齐，至少要比较面内偶极和面外偶极两类模型，因为它们对真空平台和局域电场很敏感。
取向会不会改电子结构？
图5 比较了不同氢键模式下 MAPbI₃ 的能带和 PDOS。这个图说明，取向差异会影响带边态分布、局域轨道贡献和带隙附近态密度。对于吸附或缺陷研究，取向引起的电子结构差异可能与吸附物或缺陷本身的影响同量级。
图5 不同氢键模式下 MAPbI₃ 的能带和 PDOS。该图用于说明结构取向变化会传导到电子结构。 DOI: 10.1038/srep21687。
因此表面模型完成后，不应只看总能。需要同时检查DOS、带边态 PARCHG、LOCPOT 平面平均和 Bader 或差分电荷。结构取向和电子结构必须成对对应，否则后续解释吸附增强、缺陷钝化或带边移动时会缺少模型来源。
如果两个取向能量差很小，但带边态分布不同，后续吸附计算应让两者都参与对比。相反，若某个取向能量明显偏高且优化后会转向低能构型，它不应作为主要表面模型继续计算。
AIMD 检查哪些对象？
杂化钙钛矿柔性较强，静态 0 K 优化只能说明一个局域极小点。AIMD 可用于检查表面 MA⁺/FA⁺ 是否旋转、PbX₆ 八面体是否保持骨架、表面卤素是否迁移，以及 slab 是否出现明显层间塌缩。温度、时间步和模拟时长应按研究问题设置，例如 300 K 下 1 fs 时间步作为起点，并跟踪能量、温度和关键键长随时间变化。
AIMD 不是替代结构优化。更合适的流程是：bulk 优化、slab 优化、代表取向比较、短程 AIMD 检查、抽帧再优化、电子结构计算。AIMD 用来检验模型在热扰动下是否保持同一结构身份。若 AIMD 后表面终止或有机阳离子取向完全改变，后续静态吸附应从抽帧优化后的结构重新开始。
模型可靠性怎样收尾？
一个钙钛矿表面模型能不能继续用于吸附、缺陷或能级对齐，要同时看晶相、晶面、终止面、slab 层数、真空层、固定层、A 位取向、偶极修正、表面能、AIMD 后结构变化和电子结构。对每一个表面模型，表面能要能对应到 POSCAR，LOCPOT 要能给出真空平台，DOS/PARCHG 要能解释带边态来源。
如果后续要做小分子吸附，吸附物应放在已经通过表面身份检查的模型上。PbX₂ 终止、AX 终止、表面卤素空位和不同 MA⁺ 取向不应混在同一列里比较吸附能。钙钛矿表面吸附能的可比性来自同一表面母体。
形成结论前，模型选择本身就应接受检验，而不是只把一个 POSCAR 当作结构来源。钙钛矿表面越柔性，越需要确认晶相、取向、终止和热稳定性；这样后续的吸附、缺陷和能级对齐才有稳固的结构基础。
要点：有机-无机钙钛矿表面模型要同时定义晶相、表面终止、A 位阳离子取向和偶极方向。
要点：PbX₂ 与 AX 终止需要通过表面能或目标环境选择，不能只按切面名称决定。
要点：MA⁺/FA⁺ 取向会通过氢键和局域静电势改变结构与电子态，应作为模型变量标明。
要点：AIMD 可检查柔性表面在热扰动下是否保持同一结构身份，抽帧后仍要回到静态优化和电子结构分析。
