结构优化不收敛的十种原因与解决办法！

说明：本文华算科技主要介绍结构优化不收敛的十类来源，以及针对原子接触、优化路径、电子态、数值误差和晶胞条件的解决办法。

结构优化的收敛条件是什么？

键长、键角、分子取向和晶胞形状改变时，材料的电子能量随之改变。结构优化沿这些几何自由度更新模型，每到一个几何点，DFT 都要求出该结构的电子态，计算原子力和应力。优化路径使用的是当前几何上的能量梯度，梯度受初始构型、电子态及数值近似影响；轨迹走得缓慢、来回振荡和持续重构，可能有不同来源。

原子受力写成 F_i = −∂E/∂R_i，对应能量对核坐标的负梯度。允许移动方向上的残余力应低于目标阈值；晶胞参与弛豫时，还需比较内部应力与指定外部应力。冻结原子可承受约束反力，受限方向的力没有必要全部消失。原子力通常用 eV Å⁻¹ 表示，应力则用 GPa 等单位表示。

孤立甲磺酸分子的结构图给出了优化后的原子位置及 S=O 键位置。平衡几何由所采用的电子结构模型确定，它的检验对象是残余力和局部能量变化。分子增加水合接触、改换电荷态或处在外部约束下，原先几何未必仍满足同一组平衡条件。

能量变化很小仍可能伴随较大的力，例如步长过小，原子几乎没有移动。沿低曲率方向，即使力已很小，到达平衡点所需的位移仍可较大。判断一段轨迹应比较最大残余力、位移方向、电子自洽误差和适用的应力分量。终止在最大步数处只说明迭代预算耗尽，轨迹的下降趋势与误差水平决定延长计算是否有效。

初始结构与优化路径为何反复振荡？

重复原子、周期拼接后的重叠、吸附物过低或错误的分数坐标，会产生很大的短程排斥力。少数原子的大幅移动可改变配位乃至电子态，后续梯度随之突变。异常大力的位置应对应到具体原子对；对每个原子计算它与周期复制原子的最短距离，跨晶胞的接触也属于待优化模型。

甲磺酸–水团簇图展示的是平衡水合构型，短接触随氢键网络和异构体改变。重建初始坐标时，按原子身份、键型和邻近配位区分真实氢键与重叠，修正重复位置和不合理吸附高度。适用的低成本模型可用于预弛豫，再用目标 DFT 模型优化；预弛豫后的成键与配位仍须核对，尤其是可能发生质子转移的体系。

同一根键交替拉长和缩短，力的方向随步数反复反转，常见于更新幅度过大。对局部二次势 E = ½kx²，简单梯度下降给出 x_new = (1−αk)x；α 与局部曲率 k 的乘积决定更新稳定程度。在 k > 0 时，稳定更新要求 0

减小最大原子位移或更新尺度，使试探几何保留在当前能量近似有效的范围内。线搜索或信赖域可依据实际能量与力调整步长；初始结构离平衡很远时，采用能容忍大初始力的下降算法。若大步导致配位重排，旧梯度历史对新几何的曲率近似已可能失效，应重新建立相应历史。

共价键伸缩较硬，分子转动、层间滑移或长波形变较软，势能面沿不同方向的曲率可相差很大。统一步长为避免硬方向发散而受限，软方向却每次移动很少。Hessian 的条件数反映曲率差异，刚体平移、转动等零模应从该比较中排除。

预条件矩阵或合适的内部坐标按各方向的曲率缩放更新，使硬键与软运动获得不同尺度。PM6 分子测试图横轴是条件数比值的平方根，纵轴是函数/梯度调用次数之比，说明改善条件数可以减少优化调用；这是模型势的算法测试。应用到 DFT 时，残余力和电子计算误差仍决定实际可达到的阈值。

表面重构需要位移的原子若全部固定，或低对称畸变所需位移若受母相对称操作限制，优化将始终停留在受限空间。被限制的畸变方向无法通过增加步数产生。另一种情况是受限空间已达到平衡，但检查时把约束反力也列入自由结构的收敛要求。

按研究对象核对冻结层、选择性位移和对称操作。有物理依据的约束应保留，例如外延夹持或表面深层近似；无意限制的自由度应释放。对待检验的低对称结构施加小而明确的初始扰动，再观察其是否回到原相或形成新的平衡构型，同时区分自由方向的残余力与约束反力。

电子态与数值误差怎样干扰受力？

固定核坐标后，电子密度仍在变化，离子更新就会使用尚未稳定的梯度。能量变化达到较松阈值时，原子力误差仍可能超过目标力阈值，准牛顿法还会把这些误差写入曲率历史。大胞金属的长波电荷振荡和表面–真空的非均匀屏蔽，也可使电子残差下降困难。

对同一几何收紧电子自洽精度，比较各原子的力是否稳定，使电子求解引起的力变化小于需要分辨的残余力。密度振荡时调整混合幅度与合适的预条件，残差仍持续下降而触及电子步上限时才延长循环。恢复稳定电子力后再更新结构，受到噪声污染的曲率历史应重新建立。

金属能带随原子位置移动，费米能级附近的占据也随之变化。粗 k 点网格配合过窄展宽，可能使力对微小几何变化表现出不平滑跳动，即使每个几何点都已电子自洽。k 点离散误差与展宽引入的系统误差分别影响受力，需要在同一模型上分辨。

1T-MoS₂ 图中的原子有固定小位移，参考力本来就非零。增大展宽可降低粗网格的采样误差，但过大展宽又使力偏离低展宽参照；下方柱图比较了计算耗时与系统误差。联合核验 k 点网格和展宽下的力，选取能支持目标力精度的组合；绝缘体则采用适合有带隙占据的处理。

过渡金属氧化物和缺陷体系可能有多个自洽电子解。局域磁矩、轨道占据突然改变时，原子力也可能突变，优化器积累的历史来自不同电子态。每步电子自洽与沿同一电子态平滑优化是两项不同要求；总磁矩相近的两组局域占据也可能具有不同能量与受力。

NiO 图中的纵轴是距目标占据矩阵的距离，释放占据约束后，一些轨迹偏离初始目标，却仍达到了自洽状态。可对有依据的初始磁矩和轨道占据分别开展弛豫，追踪局域磁矩及占据的变化。约束电子态有明确物理用途时才使用；解除约束后复查受力和能量，避免把人为固定的电子态当作自由体系基态。

截断能、实空间网格或局域基组不足时，电子自洽可以完成，几何梯度却仍带有离散误差。粗网格可能使原子相对网格移动时的能量和力产生非物理起伏。晶胞变化还会改变有限基组的表示，Pulay 应力可偏移体积平衡；这类误差与电子残差来源不同。

受压 Si 图横轴为截断能，纵轴为压缩应力 −σ₃；两种应力求值随基组改善接近各自稳定值。固定一组代表性几何，增加基组与网格精度，比较力、应力及微小位移能量差的变化。力阈值应高于已确认的数值误差水平，实际采用的精度由目标体系验证，图中 Si 的截断能无需推广到其他元素。

晶胞条件与真实失稳为什么使优化持续进行？

体相零压平衡、外延面内夹持和固定真空的表面模型允许不同的晶胞变化。让带真空的 slab 任意收缩面外尺寸，或用平均压力替代规定的应力张量，会改变优化目标。晶胞平衡由允许应变与外部载荷共同限定，外加静水压时通常最小化 H = E + p_extV；一般非静水载荷还涉及具体机械条件。

Si 的加载图在不同侧向应力下呈现不同变形路径，F₁、F₃ 是形变梯度分量，和原子力的符号含义不同。这些曲线由受约束平衡结构构造，并非离子迭代记录。修正办法是按模型选择可变晶胞分量和目标应力：slab 保留所需真空，外延模型保留面内夹持，体相则检查允许应变方向上的应力残差。

持续的八面体转动、层间滑移或剪切重构，可能对应能量沿某个位移方向继续降低。高对称位置甚至可具有零力，却沿某个模式呈负曲率。Si 的失稳应力图说明，局部稳定范围随完整载荷条件改变；超出稳定范围后，原来的晶相不再是该条件下可维持的平衡分支。

在排查电子和数值误差后，计算相关振动模式或沿观察到的畸变扫描能量，用位移模式定位下降方向。按真实不稳定模构造扰动并重新弛豫，必要时扩大超胞以容纳非 Γ 点畸变；同时核验近零虚频的力精度和有限尺寸误差。若原相只在有限温度下稳定，零温静态优化得到的畸变结构对应另一组物理条件，温度效应需由相应自由能或有限温计算描述。