DFT周期性计算如何选定合理超胞尺寸？收敛测试与多体系建模方案

说明：本文华算科技主要介绍超胞大小在周期性 DFT 中的物理含义、尺寸收敛方法，以及缺陷、吸附、声子和分子动力学任务中的模型选择。

超胞大小的基础概念是什么

周期性 DFT 里，超胞由三条晶格矢量定义。把原胞沿 a、b、c 方向扩成 n_a × n_b × n_c 后，计算程序处理的是这个单元在三维空间中的无限重复。超胞大小同时改变真实空间距离和倒易空间尺度，缺陷、吸附物、层间位移和原子振动都带着相邻周期副本的影响。

尺寸判断从目标物理状态开始。孤立缺陷对应低浓度极限，表面吸附对应某个覆盖度，声子计算对应力常数衰减长度，AIMD 对应给定温度下的相关长度和采样体积。同一个 4 × 4 超胞可以代表稀释缺陷，也可以代表有意设定的 1/16 ML 吸附覆盖度，两种语境下的物理含义并不相同。

在低对称晶体中，超胞方向与层状、链状或极性结构绑定；若缺陷扰动主要沿层间方向传播，n_a = n_b = n_c 不自动代表三向镜像距离相当。层状氧化物、二维材料、极性 slab 和链状晶体常要分别比较面内距离、层间距离、真空厚度和电势平台。超胞矢量的长度、夹角和材料介电响应会共同改变镜像相互作用。

倒易空间也要随超胞同步折算。实空间 L_a 加倍时，对应倒易基矢约缩小一半；原胞 8 × 8 × 8 的 k 点间距，在 2 × 2 × 2 超胞中常对应 4 × 4 × 4 的采样密度。比较不同超胞时，记录 k 点间距比记录网格数字本身更可复查，否则总能差会混入布里渊区采样误差。

超胞尺寸如何影响缺陷与吸附

点缺陷在超胞中代表周期阵列里的一个缺陷副本。中性缺陷会通过弹性弛豫、局域电荷密度和缺陷态波函数与副本相互作用；带电缺陷还含有背景电荷、长程库仑项和势能对齐项。缺陷形成能的尺寸误差常随缺陷浓度、最近镜像距离和超胞形状一起变化，只报原子数会掩盖不同形状之间的差别。

可操作的尺寸测试通常把目标量画成外推变量的函数，例如 1/L、缺陷密度、Madelung 势或经过介电张量修正的势项。若相邻两到三个尺寸得到的 ΔE 小于研究容差，外推截距变化也小，超胞尺寸才有数值停止点。带电缺陷还应标出电荷态、介电常数、势能对齐口径和背景电荷处理，避免把修正公式差异写成材料差异。

带电缺陷的停止点除了形成能平台，还可检查缺陷态是否保持局域、局域电势是否在远离缺陷处恢复体相参照、结构弛豫是否被当前模型容纳。若缺陷态跨过周期边界形成能带色散，转变能级和占据状态都会带上镜像耦合；这时增加横向或体相方向距离，比单纯套用电荷修正更接近问题来源。

表面吸附的横向超胞直接定义覆盖度。以同一晶面和一个吸附物为例，2 × 2、3 × 3、4 × 4 超胞对应不同吸附物面密度，相邻吸附物之间的排斥、偶极耦合、氢键网络和表面电荷重排都会改变 E_ads。若研究孤立吸附物，横向距离、关键键长和 E_ads 的尺寸变化应共同收敛；若研究高覆盖度反应态，小超胞就是目标模型的一部分。

目标量收敛怎样给出尺寸停止点

没有通用于所有材料的原子数阈值。筛选相稳定可把能量差容差定在几十 meV/atom；缺陷转变能级常用 0.05 eV 量级检查；吸附和反应路径要看相邻构型 ΔE；力常数、弹性常数和声子频率对残余力更敏感。超胞测试必须固定泛函、赝势、截断能、展宽、自旋和 k 点间距口径，单次改变多个设置会让尺寸误差失去来源。

尺寸收敛的实际含义，是让模型误差低于研究对象本身的能量尺度。若两个吸附位点的 E_ads 只差 0.03 eV，而 3 × 3 与 4 × 4 超胞之间仍差 0.06 eV，位点排序还带着尺寸误差；若两个相结构能量差达到 0.5 eV，同样的尺寸残差对定性排序影响较小。容差应由待比较的能量差、频率差或结构位移给出，经验原子数不适合反推容差。

局域弛豫提供另一类停止信号。缺陷或间隙原子引起的位移如果延伸到超胞边界，边界原子已经不再像体相环境；若外层原子位移接近零，局域键长和配位恢复到体相参照，结构扰动才被当前超胞容纳。对缺陷几何、磁矩和局域态而言，外层位移往往比总原子数更接近真实风险，这类检查尤其适合低对称晶体和强局域弛豫缺陷。

形状测试承担方向性误差排查。对于层状材料，2 × 2 × 1 与 3 × 3 × 1 主要检查面内镜像；对于体相缺陷，2 × 2 × 2、3 × 3 × 2、2 × 3 × 3 可分辨方向性误差；对于 slab，横向覆盖度、slab 厚度和真空层是三条不同轴。尺寸停止点应由目标量平台、边界原子响应和模型物理状态三者共同限定，并在图源记录中保留所用尺寸序列。

不同计算任务怎样选用超胞尺寸

缺陷计算以缺陷浓度、最近镜像距离、形成能外推和电荷修正为主；表面计算以覆盖度、slab 厚度、真空电势平台和表面层弛豫为主；能带和态密度还会受到带折叠影响，超胞越大，本征能带投影和展开口径越重要。任务不同，尺寸测试的目标量就不同，同一套 3 × 3 × 1 模型不应直接迁移到缺陷、吸附、功函数和声子四类结论中。

表面 slab 的尺寸测试通常拆成横向、法向和真空三个方向。横向尺寸改变覆盖度和吸附物侧向相互作用；slab 厚度改变中间层是否恢复体相配位；真空层改变法向镜像、电势平台和偶极修正残差。增大真空层并不能替代横向覆盖度收敛，加厚 slab 也不会消除吸附物之间的面内排斥。

有限位移声子计算把超胞尺寸转化为可采样的力常数范围。若力常数尚未衰减到足够小，周期边界会截断长程相互作用或把远距离原子误当成近邻副本，声学支软化、虚频和简并破缺都可能出现。声子超胞的尺寸由力常数衰减长度和目标 q 点共同限定；原胞原子数只描述计算成本，不描述长程力常数范围。

力常数随距离的衰减曲线能直接暴露超胞是否截短相互作用。Sb₂S₃ 案例中，1 × 4 × 1 与 2 × 2 × 2 超胞在较远距离处仍有明显差别，2 × 4 × 2 才能覆盖更长的相互作用范围。软晶格、极性材料、低维材料和含界面的体系应额外检查长程力常数，否则小超胞中的虚频可能来自数值模型。

AIMD 和扩散计算还要把时间尺度写入尺寸判断。离子迁移、溶剂层重排和界面波动在小超胞中会反复遇到自己的周期副本，MSD 斜率、RDF 峰形和吸附物取向分布都可能偏向单一构型。可复核的超胞选择应留下六项记录：模型尺寸、最近镜像距离、浓度或覆盖度、k 点间距、目标量收敛曲线和残余图谱信号。这些记录齐全时，超胞结果才带有可检查的物理尺度和误差来源。