钙钛矿分子动力学模拟

我做过一个有机-无机杂化钙钛矿的分子动力学模拟项目，客户关注的是 MA⁺ 有机阳离子在晶格里的转动和 I⁻、Pb²⁺ 的扩散行为，因为这两个因素直接决定钙钛矿的离子迁移率和器件稳定性。

左侧图是 300 K 下的钙钛矿 MD 快照。可以看到 Pb²⁺ 和 I⁻ 构成八面体骨架，MA⁺ 位于 A 位 cage 里。橙色箭头表示热振动方向，绿色虚线圈标出了 MA⁺ 的取向旋转。钙钛矿的高温立方相里，MA⁺ 可以自由旋转；低温四方相和正交相里，MA⁺ 的取向会被有序化。我模拟的温度是 300 K，接近室温，体系处于四方-立方过渡区，所以 MA⁺ 既有转动也有部分有序。

右侧图是不同离子的均方位移 MSD 随时间变化。绿色实线是 MA⁺，红色是 I⁻，蓝色虚线是 Pb²⁺。200 ps 内，MA⁺ 的 MSD 最大，接近 1.3 nm²；I⁻ 次之，约 0.45 nm²；Pb²⁺ 几乎不动，约 0.28 nm²。从斜率可以估算扩散系数 D = MSD/(6t)（三维），MA⁺ 的 D 大约在 1×10⁻¹¹ m²/s 量级，I⁻ 约 3×10⁻¹² m²/s。这个量级和文献报道的杂化钙钛矿离子迁移率一致。

模拟设置上，我用 LAMMPS 跑全原子 MD，力场选的是 GAFF 对有机物部分加 Buckingham 对无机骨架，参数从文献中来。超胞取 4×4×4（约 3200 原子），NPT 平衡后 NVT 采样，温度 300 K，总时长 1 ns，前 200 ps 用于平衡。MSD 分析取后 800 ps，分不同离子种类统计。

踩坑经验：第一个是力场选择。钙钛矿力场很多，有的只拟合晶格常数，不考虑离子迁移；我选力场时特别看重对 Pb-I 距离和 MA⁺ 转动势的描述。第二个是模拟时长：200 ps 对看 MA⁺ 转动够了，但要看 I⁻ 长程扩散，最好到 5-10 ns。第三个是温度控制：Nosé-Hoover 链对钙钛矿这种软晶格比较合适，Berendsen 弱耦合容易低估温度涨落。

工程经验上，我建议客户把离子迁移和器件退化联系起来：MA⁺ 和 I⁻ 迁移会在电极附近积累空间电荷，导致滞后和效率衰减。抑制迁移的策略包括用 larger A-site cation（如 FA⁺、Cs⁺ 混合）降低 cage 尺寸，或者用二维钙钛矿包覆三维晶粒。我这个项目最后也给了一个不同温度下的 MSD 趋势预测，供客户做加速老化实验参考。

更多 LAMMPS/GROMACS 钙钛矿和离子迁移模拟的案例，我整理在 [钙钛矿分子动力学模拟](https://www.keyanxueshu.com/category/md/) 栏目里。[钙钛矿分子动力学模拟](https://www.keyanxueshu.com/)

右侧 MSD 曲线还有一个细节值得注意：MA⁺ 的 MSD 呈明显直线，说明它是在做近似 Brownian 扩散；而 I⁻ 的 MSD 在 50 ps 内有轻微弯曲，可能是因为 I⁻ 在八面体角共享位置之间有跳跃式迁移，不是连续扩散。这种跳跃机制在无机 CsPbI3 中更明显，可以用 van Hove 关联函数进一步确认。另外，MA⁺ 的转动可以用取向序参量 <P2(cosθ)> 来量化，从有序到无序相变时这个值会从接近 1 掉到 0。我们 300 K 的模拟里，序参量大约在 0.2-0.4 之间，说明已经部分无序，和四方相的特征一致。如果客户研究的是高温立方相，建议把温度升到 350 K 以上再跑。