铁电第一性原理计算：极化强度与自发位移的那些事

铁电材料做第一性原理计算，和普通绝缘体最大的区别在于——你必须处理一个没有对称中心的极性结构。这种结构在传统能带计算里并不特殊，但一旦涉及 Berry Phase 极化计算、压电系数、或者相变路径模拟，事情就变得微妙起来了。

铁电性的本质

铁电材料的核心特征是存在自发极化：在某个温度以下，晶格会发生微小畸变，正负电荷中心不再重合，从而产生宏观电偶极矩。钛酸钡(BaTiO₃)和 PbTiO₃ 是最经典的研究对象，它们的铁电畸变表现为 Ti 原子沿 c 轴偏离氧八面体中心。

从 DFT 的角度看，铁电相对应的是一种能量极小态——但这个极小态有多重性(铁电相与反铁电相之间存在势垒)，所以初始结构的选择会直接影响计算收敛到哪一个相。

计算铁电性的几个关键技术点

Berry Phase 方法是计算极化强度的标准手段。VASP 里通过 LCALCPOL = .TRUE. 激活Berry Phase，计算结果会输出电子极化、离子极化和总极化。这个方法本身是健壮的，但有前提条件：计算必须覆盖完整的绝热演化路径。如果你在做铁电材料的能带计算，需要先做自洽再在倒空间取密集 k 点做非自洽计算，极化信息不会自动传递，需要手动保存 CHGCAR 再做处理。

结构优化是另一个容易翻车的地方。 铁电材料的能量曲面非常平坦——铁电相和顺电相的能量差往往只有几个 meV/f.u.，如果 EDIFFG 设得太宽松，优化过程中原子坐标可能卡在某个中间状态。个人的经验是：铁电体系用IBRION=2(共轭梯度法)配合 EDIFFG = -1E-3.比默认的 CG 收敛行为更稳定。

关于赝势的选择也值得专门说。 PAW-PBE 赝势对大多数氧化物体系足够用，但做极化计算时，Ti 的赝势需要检查是否正确生成了 3d 电子的投影。PAW_PBE Ti_pv(包含 3s²3p⁶3d²4s²，半芯态)比 Ti(包含 3s²3p⁶3d³4s¹)赝势在处理铁电位移时收敛行为更可预测。

相变模拟的思路

如果你的目标是模拟铁电-顺电相变，典型的做法是沿铁电畸变方向做能量曲线(PEAD 方法)或者做晶格常数的体积优化曲线。不同体积下的铁电畸变幅度不同，联系到 Landau 理论，可以提取出 Ginzburg-Landau 系数。

对于温度效应的模拟，第一性原理只能做 0K 计算。如果需要考虑热力学贡献，通常的做法是结合准谐近似(PHONON 或 Phonopy)来估算自由能曲线，进而预测相变温度。

结语

铁电第一性原理计算的核心难点在于：体系能量差极小、收敛行为对参数敏感、极化计算需要额外设置。开始之前，先确认赝势、初始磁矩和自旋极化设置是否正确，再选择合适的优化算法和收敛阈值。如果在 Berry Phase 极化计算中得到异常值，先检查波函数是否收敛到足够精度，可以到科研学术网一起排查具体的问题。