晶格常数怎么算准？从平衡定义到实验对比的完整排查

说明：本文华算科技主要介绍晶格常数在 DFT 中的平衡定义、数值误差来源、物理模型偏差和实验比较口径。

晶格矢量与零应力平衡胞的定义

晶体的周期重复由三条晶格矢量构成。立方晶系常用一个边长 a 表示，六方晶系至少要给出 a 和 c，单斜、三斜晶系还需保留晶轴夹角。讨论“晶格常数算得准不准”之前，需要先确定比较的是哪种晶胞、哪一组独立参数以及哪一个晶相。原胞和常规胞可以描述同一周期结构，却不会给出相同的边长数值；把两套晶胞直接逐列比较，误差来自表示方式，并非电子结构计算。

给定交换关联泛函、赝势、自旋态和外压后，DFT 求得的是该模型下的势能面。平衡晶胞满足能量对允许应变的一阶导数为零；零外压下各独立应力分量应接近零。原子内坐标也要随晶胞变化重新弛豫，否则得到的是“固定内部构型的能量曲线”，其中夹带了未释放的原子力。晶格参数、原子坐标和磁性自由度需要指向同一个局部能量极小值，这才是可讨论的理论平衡结构。

图1把原胞、常规胞和能量—晶格参数曲线放在同一幅图中。曲线最低点给出该计算模型的平衡参数；用 Birch–Murnaghan 等状态方程拟合一组体积点，可联合得到平衡体积 V₀ 和体积模量 B₀。对常规立方胞才有 a₀=V₀^1/3；低对称晶体不能用总体积唯一反推出每条晶轴。若能量谷很平，微小能量噪声就会造成明显的最低点位移，报告很多小数位并不会增加可信度。

可用的晶格结果需要满足三个条件：曲线或变胞弛豫找到了稳定极小值，残余应力与原子力低于研究目标允许的尺度，结构没有滑入错误的磁态或亚稳晶相。“准”对应一段可验证的误差范围，输出文件末行的单个数字还不足以界定这段范围。

数值误差怎样通过力和应力改变晶胞？

总能看似已经稳定，晶格参数仍可能缓慢移动，因为变胞弛豫直接使用应力，应力是能量对应变的导数。导数量会放大自洽迭代、布里渊区积分和基组离散带来的细小波动。金属费米面附近的占据变化、稀疏的 k 点采样或过宽的电子展宽，都可能让相邻体积点的能量曲线出现不平滑斜率。此时优化器读取到的残余压力含有数值分量，晶胞会沿着错误斜率继续伸缩。

图2比较了有限差分和解析导数路线。差分步长过小会放大 SCF 噪声，步长过大又会混入非线性响应；同一材料的误差随两者共同变化。晶格优化采用一阶应力，仍受相同数值机制支配。电子能量阈值应由残余应力和目标晶格误差反推，仅记录总能每原子的变化量无法约束晶胞误差。

平面波计算还有一个与晶胞尺寸直接相连的误差：固定截断能时，晶胞伸缩会改变落入能量球内的平面波数目，离散基组并非随体积连续变化。Pulay 应力会使能量最低点和零应力点产生偏移。截断能提高后，密度和力误差通常下降，但不同元素、赝势和化学环境下降速度不一；含半芯态或较“硬”赝势的体系常需要更高的基组分辨率。

图3中不同材料沿近似对角关系分布，说明基组造成的电子密度误差会传入原子力。对晶格优化而言，还需单独追踪应力和最终晶轴，因为总能、应力和晶轴的收敛速度未必一致。判断数值设置是否够用，可以比较加密前后 a、b、c、夹角和残余应力的变化；若研究应变引起的能带或声子变化，允许的晶格误差须小于论文中施加的应变量。

物理模型为什么会把晶格算大或算小？

数值参数收敛后，剩下的偏差主要来自哈密顿量近似。交换关联泛函改变键合能随原子间距的曲率和最低点位置。LDA 常使不少致密固体的晶格偏小，PBE 常给出偏大的平衡体积，但这只是统计趋势；离子键、共价键、金属键和弱层间作用占比不同，误差方向可能变化。赝势若遗漏会参与成键或屏蔽的半芯态，也会改变平衡键长。泛函与赝势需要作为同一个物理模型接受检验，不能把某个体系上的经验修正当成通用比例因子。

图4中，立方相的 LDA 结果普遍偏小、PBE 偏大，PBEsol 和 vdW-DF-C09 更靠近零误差；进入四方、正交和菱方相后，不同晶轴的偏差会各自变化。铁电畸变会把体积误差与极化位移耦合，一条晶轴“碰巧对上”不能替代完整晶胞检查。层状材料还需检验色散作用，强关联氧化物则要关注局域电子、轨道占据和磁序是否与目标相一致。

DFT+U 会重排局域 d 电子的占据和杂化，平衡晶格随之变化；U 施加在哪类原子上同样重要。Ni₂MnGa 的结果显示，Ni 位或 Mn 位采用相同数值的 U，晶格常数、体积模量和磁矩响应并不相同。用实验晶格单独反标 U，可能牺牲磁矩、弹性或相对相能，所得参数未必适合后续性质计算。

磁性材料可能存在多个接近的自旋解。晶格扩大或四方畸变会改变交换分裂和态密度，磁矩变化又会反过来移动能量最低点。图6沿 c/a 畸变路径并列给出相对能量和总磁矩，不同泛函与 U 位点给出的能谷形状明显不同。计算若从错误磁序出发，即使力和应力都已收敛，也可能得到另一个局部极小值。

物理模型的选择可用一组彼此相关的量约束：晶格参数检验能量最低点，体积模量检验能量曲率，相对相能检验多个极小值排序，磁矩和局域结构检验电子态是否正确。多个量均处在合理范围，比单独追求某个晶轴与实验完全重合更能排除误差抵消。

理论晶格与实验数据怎样统一口径？

常规静态 DFT 优化对应固定电子模型、静止原子核和指定外压下的零温势能最低点，室温 XRD 包含零点振动、热膨胀、缺陷、应力与样品成分的影响。两者相差千分之几时，温度口径往往已达到同一量级。实验给出的还是特定晶相的平均周期结构；若样品含相共存、非化学计量或基底约束，晶格参数不能直接当作无缺陷体相的零温基准。

需要比较有限温度时，可把振动自由能加入静态能量，写成 F(V,T)=E_DFT(V)+F_vib(V,T)，再在每个温度寻找自由能最低点。准谐近似允许声子频率随体积变化，能够描述零点膨胀和多数弱非谐固体的热膨胀；接近软模相变、强非谐区或离子扩散区时，这一近似需要额外检验。

图7中 Li₃ClO 的晶格参数随温度持续增加，300 K 附近的数值高于零温静态平衡值，曲线斜率对应线性热膨胀系数。理论值与实验值的差应分别归入数值误差、模型偏差和热力学条件差；三项混在一个百分误差里，无法判断该从基组、泛函还是温度修正入手。

实际比较时应注明晶相与空间群、原胞或常规胞、温度、压力、成分和实验不确定度，再报告各独立晶轴及体积偏差。对于近立方但存在微小四方畸变的材料，等效立方晶格会掩盖 c/a 的错误；对于薄膜，面内晶格由基底钳制，体相全弛豫计算施加的力学约束也不同。完成这些口径对齐后，晶格常数才适合用来判断后续能带、声子、弹性或缺陷计算采用了怎样的结构基准。