说明：本文华算科技主要介绍 DFT 结构优化不收敛的常见原因，包括初始结构、电子步、离子步、磁性、约束、真空层和数值参数。文章给出判断顺序，帮助区分物理不稳定和计算设置问题。结构优化不收敛至少包含两类问题：电子自洽不收敛和离子结构不收敛。前者表现为 SCF 能量振荡、磁矩跳变或电荷密度难以稳定；后者表现为能量下降很慢、力始终较大、原子来回摆动或晶胞畸变异常。混淆两类问题会导致错误处理。
优化本质是在势能面上寻找局域极小点。若初始结构距离合理构型太远，算法可能沿着陡峭方向反复试探；若势能面很平，例如层状滑移、分子吸附或软声子模式，力收敛会天然较慢。此时问题不一定是程序错误，而可能是模型物理特征。
Fi=−∂E/∂Ri
式中 Fi 是第 i 个原子受到的力，E 是体系总能，Ri 是原子坐标。常见收敛标准为最大力低于 0.01–0.05 eV Å−1，高精度声子或弹性计算通常需要更严格阈值。
图1：MOF 结构建模与势能评估流程体现结构弛豫、力场近似和优化稳定性的关系。DOI：10.1038/s41524-025-01872-3。
导入晶体结构时还应检查对称性和占位。实验 CIF 中可能含有部分占位、无序溶剂或平均结构，直接转成 DFT 模型会产生非整数原子数或异常短接触。对这类结构，先构建有序超胞并合理分配占位，比直接按原始 CIF 优化更稳妥。
对于表面和界面模型，真空层不足会造成周期镜像相互作用，通常需要至少 15–20 Å 真空；二维材料若同时放开 c 轴晶胞，程序可能把真空当成可压缩空间处理。分子吸附体系还应避免偶极矩方向上的虚假电场，必要时加入偶极修正。
图2：MOF 结构建模与势能评估流程体现结构弛豫、力场近似和优化稳定性的关系。DOI：10.1038/s41524-025-01872-3。
约束设置也会制造假性不收敛。固定层太少会让表面整体漂移，固定层太多又会锁死真实弛豫；选择性动力学标签写错时，某些原子看似无法降低受力。检查优化轨迹比只看最终报错更可靠。
展宽方法也会影响离子优化路径。金属体系常用 Methfessel−Paxton 或 Fermi 展宽帮助 SCF 稳定，但最终能量比较最好回到更合适的四面体方法或较小展宽。若展宽过大，表面吸附构型可能被电子熵项牵引到不真实的局域极小点。
金属、窄带隙半导体、强关联氧化物和磁性体系容易出现 SCF 振荡。原因是电荷密度更新过猛，导致每一步都在不同电子构型之间跳动。降低混合参数、增加展宽、使用更稳健的初猜电荷，常比单纯增加最大电子步更有效。
磁性设置错误也会拖垮优化。过渡金属氧化物若未设初始磁矩，可能收敛到错误的非磁态；反铁磁体系若超胞和磁序不匹配，电子态会在多个近简并解之间振荡。此时应先固定离子做单点测试，比较不同磁构型的能量和磁矩。
Eads=Eslab+X−Eslab−EX
吸附优化中 Eads 的可靠性依赖三个结构都充分收敛。若只有复合体系放松，而清洁表面或分子未达到同等标准，吸附能误差可达到 0.1 eV 以上，足以改变催化活性排序。
图3：MOF 结构建模与势能评估流程体现结构弛豫、力场近似和优化稳定性的关系。DOI：10.1038/s41524-025-01872-3。
当优化失败出现在晶胞自由度上，应单独检查压力张量。层状材料、分子晶体和多孔材料对范德华作用很敏感，未加色散校正时层间距可能持续膨胀；压力收敛阈值设得过严，也会让晶胞在微小应力下反复振荡。
实用排查应从低成本检查开始：可视化初始和轨迹结构，确认原子距离与约束；随后做较松标准预优化，再逐步提高截断能、k 点和力阈值；若电子步仍振荡，再处理混合、展宽、磁矩和占据方式。跳过结构检查直接改高级参数，常会浪费大量算力。
当优化持续不收敛但能量已经接近平坦时，需要判断是否存在真实不稳定。可计算声子虚频、检查弹性常数或沿软模方向扰动再优化。若体系本来处在相变路径、扩散过渡态或吸附解离过程中，强行寻找局域极小点可能不符合研究目标。
图4：MOF 结构建模与势能评估流程体现结构弛豫、力场近似和优化稳定性的关系。DOI：10.1038/s41524-025-01872-3。
若研究对象是反应路径，还需要区分极小点优化与过渡态搜索。NEB 初始链中的中间构型常常力较大，不能按普通结构优化标准逐个要求收敛到极小点；真正需要检查的是路径连续性、最高能点附近切向力以及最终过渡态虚频。
图5：MOF 结构建模与势能评估流程体现结构弛豫、力场近似和优化稳定性的关系。DOI：10.1038/s41524-025-01872-3。
