钙钛矿表面小分子吸附容易被写成普通表面吸附问题：切一个 slab，放上 H2O、O2、CO2 或 NH3，优化后比较吸附能。这个流程本身没有错，但对卤化物钙钛矿和有机-无机杂化钙钛矿来说，表面模型的身份比吸附物本身更容易改变结果。有机阳离子取向、卤素空位和表面终止面会同时改变局域静电场、未饱和配位和带边态位置。 吸附能排序需要在这些变量被控制后再解释。
本文的重点是：当小分子靠近卤化物钙钛矿表面时，模型里哪些自由度必须作为变量，哪些变量应该固定成对照。只有这样，吸附能、差分电荷密度和 DOS 才能说明同一个界面事件。
终止面先改变什么？
ABX3 钙钛矿切表面时，常见终止包括 AX 终止和 BX2 终止；在卤化物钙钛矿中，也常见 PbI2 类终止和有机阳离子相关终止。不同终止的表面电荷分布、未饱和 Pb 位、卤素外露程度不同，小分子看到的配位环境就不同。
同一个 H2O 分子，在富 Pb 的表面可能通过 O 端与 Pb 配位；在富卤素表面可能主要形成氢键或弱静电吸附。若只算一个终止面，很容易把终止面效应误认为分子种类效应。小分子吸附能的第一层对照，应是同一分子在不同终止面上的能量和构型排序。 之后才适合比较不同分子的强弱。
图1 混合卤化物钙钛矿晶体结构、带边位置和气体响应测试示例。结构中 B 位金属、卤素框架和表面分子共同决定后续吸附模型的边界条件。DOI: 10.1002/smll.202404430。
VASP 建模时，slab 厚度要能让中间层接近体相配位，真空层通常不低于 15 Å，并对非对称 slab 使用 LDIPOL 和 IDIPOL=3。若研究的是极性表面，还要特别检查 slab 两侧是否等价；上下表面不等价时，吸附能里会混入宏观偶极和表面重构贡献。
有机阳离子怎样进入吸附判断？
在 MAPbI3、FAPbI3 或混合阳离子体系中，有机阳离子的取向会改变局域偶极。分子吸附在表面时，靠近有机端基的一侧和靠近无机 Pb-I 框架的一侧，静电势可以明显不同。若每个吸附构型都让有机阳离子自由翻转，吸附能变化可能来自阳离子重排，而不是小分子-表面作用。
可以先固定几种代表性阳离子取向，例如朝外、平行表面、朝内；每种取向下再放置同一小分子，比较最终吸附能和阳离子弛豫幅度。若不同取向下同一吸附位能差超过热能量级，就需要把阳离子取向作为模型变量，而不是只报最低能构型。有机阳离子取向改变的是界面静电势和氢键网络，不应被简单并入吸附物取向误差。 这个比较能把分子取向和表面取向拆开。
图2 不同缺陷和掺杂模型的表面结构对照。表面模型改变后，小分子吸附构型和电子转移通道也会随之改变。DOI: 10.1002/smll.202404430。
计算设置上，PBE 或 PBEsol 可用于几何优化，卤化物体系通常需要加入 D3 或 Tkatchenko-Scheffler 这类色散修正；若讨论带边、缺陷态或电荷转移，可在关键构型上补做 HSE06 或 SOC 对照。几何筛选和电子结构分析可以分层处理，但不能把低精度筛选得到的趋势直接当成精确带边结论。
卤素空位改变哪两类量？
卤素空位 VX 会留下未饱和金属配位，并可能在带隙中引入缺陷态。小分子靠近空位时，吸附增强可能来自配位空位，也可能来自分子轨道与缺陷态耦合。两种机制在吸附能上都可能表现为更负，但差分电荷密度、PDOS 和 Bader 电荷的指向不同。
因此，缺陷表面的小分子吸附至少要有三个能量：完美表面吸附能、缺陷形成能、缺陷表面吸附能。若缺陷形成代价很高，单看缺陷表面吸附很强并不能说明实际表面一定由该缺陷主导。缺陷吸附模型要同时回答缺陷是否可形成，以及小分子是否优先占据缺陷位。 这能避免把低概率缺陷位当作主表面。
Eads = E(slab + molecule) – E(slab) – E(molecule) 式(1)
Ef(VX) = E(defective slab) – E(perfect slab) + μX + q(EF + EVBM) + Ecorr 式(2)
式(1)比较小分子进入同一表面后的能量变化，式(2)说明卤素空位本身是否容易形成。Eads 很负只能说明吸附态稳定；若 Ef(VX) 在目标化学势下过高，缺陷位就不一定能代表主要表面。
图3 不同表面模型对 O2 吸附能的影响示例。无缺陷、卤素空位、金属空位和掺杂模型之间的吸附能差异，体现出缺陷位和终止面不能混作同一变量。DOI: 10.1002/smll.202404430。
O2 在完美表面弱吸附、却在卤素空位附近出现明显电荷转移时，应检查 O-O 键长、O2 分子态和缺陷态是否发生耦合。H2O 在缺陷附近吸附增强，但优化后伴随卤素迁移或表面分解时，吸附能不再代表完整分子吸附，而应按解离吸附重新定义反应式。
分子取向怎样覆盖关键端？
小分子不是只有一个放置方式。H2O 可以 O 端朝向未饱和 Pb，也可以 H 端靠近卤素；CO2 可以线性平躺、弯曲靠近缺陷位，也可以通过 O 端弱配位；NH3 通常要比较 N 端配位和氢键取向。若初始取向覆盖不足，最低能吸附构型可能根本没有被搜索到。
取向筛选可分两轮：第一轮用较低精度快速优化多个初猜，保留几何合理且能量靠前的构型；第二轮统一精度复算，包括相同 ENCUT、KPOINTS、vdW、偶极修正和收敛阈值。吸附取向比较必须在同一 slab 身份下完成，否则构型排序会混入表面终止和缺陷浓度差异。 统一精度复算是能量排序的最低条件。
图4 O2 吸附后的差分电荷密度和电子转移示例。吸附能增强后，需要用电荷密度和轨道贡献判断增强来自配位作用、静电作用还是缺陷态耦合。DOI: 10.1002/smll.202404430。
差分电荷密度的参考体系要与优化后的吸附构型对齐：吸附体系、去掉分子的 slab、孤立分子三者使用同一晶胞和同一 FFT 网格。若 slab 或分子参考重新优化，差分图会混入结构弛豫造成的密度差，难以对应实际成键区域。
吸附能能单独解释稳定性吗？
钙钛矿表面吸附常伴随结构软化、卤素迁移和局域相变。一个分子吸附能很负，可能代表真实强配位，也可能代表表面已经发生不可逆重构。判读时要同时给出关键键长、表面层位移、缺陷附近配位数和吸附后带边态变化。
某个缺陷表面对 CO2 给出很强吸附，但 CO2 优化后明显弯曲并形成碳酸盐样结构时，这已经不是弱吸附传感模型；另一个完美表面对 O2 吸附能较弱，但带边态几乎不变，反而更接近可逆吸附。强吸附只有在分子身份、表面身份和电子结构证据同时成立时，才适合解释界面相互作用。 判断时要同步查看构型是否保留分子吸附身份。
图5 不同表面模型吸附 O2 前后的态密度对照。吸附引起的新态、带边偏移和缺陷态填充，是解释气体-钙钛矿相互作用的重要证据。DOI: 10.1002/smll.202404430。
若文章目标是催化或传感，还要把吸附能和电荷转移拆开报告。催化更关心中间体是否适度稳定、后续断键或质子电子转移是否可达；传感更关心吸附前后电导、载流子浓度或表面态是否发生可逆变化。同一套 DFT 数据可以服务两个问题，但判据不同。
计算矩阵怎样压住变量？
可执行流程可以这样定：第一步，选定体相结构和表面终止，收敛 slab 厚度、真空层和固定层；第二步，枚举 2–3 个阳离子取向并筛掉不稳定表面；第三步，在同一表面上建立完美、VX、可能的金属空位或掺杂对照；第四步，对每个表面放置同一组分子取向并统一优化。
后处理按同一矩阵输出：吸附能、关键键长、Bader 电荷、差分电荷密度、PDOS、吸附前后功函数或带边位置。若涉及带电缺陷，再补充电荷修正、有限尺寸效应和费米能级范围。少变量、多对照的矩阵，比一次性堆很多最低能结构更能说明有机阳离子、卤素空位和终止面的分工。 这种矩阵也便于后续补算。
最终写结论时，可以把三类变量分开：终止面定义外露配位环境，有机阳离子调节界面静电和氢键，卤素空位提供局域缺陷态和未饱和金属位。小分子吸附强弱来自这三者的组合，而不是某一个吸附能数值本身。
要点：钙钛矿表面小分子吸附不能只比较一个最低吸附能。
要点：终止面定义外露配位，有机阳离子改变界面静电和氢键。
要点：卤素空位会同时改变未饱和配位和缺陷态，需要和缺陷形成能一起判断。
要点：分子取向、表面终止和缺陷浓度应放在少变量、多对照矩阵中比较。
