结构优化收敛等于稳定吗？从收敛含义到势能曲率、竞争相、能垒与有限温运动的综合判断

说明：本文华算科技主要介绍结构优化的收敛含义，以及势能曲率、竞争相、转变能垒和有限温运动对结构稳定判断的影响。

结构优化“收敛”的基本含义是什么？

结构优化是在选定的电子结构近似下，调整原子坐标与允许变化的晶胞，使能量函数达到指定的停止条件。对每组核坐标，DFT 求得电子态、总能量与原子力；核坐标改变，能量随之改变。收敛描述一次数值迭代的结束条件，结构稳定则涉及微小扰动、竞争相以及温度下的原子运动。

原子力写成 F_i = −∂E/∂R_i。采用力条件时，程序检查允许移动方向上的残余力是否小于阈值；采用能量条件时，检查相邻几何步的能量差。平坦能量面上，能量变化可以很小，原子位置仍未充分确定。达到最大迭代次数只表示计算结束，收敛条件应另行核实。

在 SrTiO₃ 的静态能量扫描中，Ti 的极性位移与氧八面体转动共同改变能量。图1d中，P = 0 附近的斜率为零，却沿极性畸变方向降能。各条曲线分别以相同转角下的非极性构型为零点，曲线起点重合不能用于比较不同转角结构的总能。该扫描描述静态畸变能，实验中的量子顺电状态还涉及核的量子涨落。

固定晶胞的计算只调整内坐标；外延薄膜可固定面内晶格而释放面外尺寸；表面模型常冻结底部原子。冻结方向可以保留约束反力，面内应力也可以由衬底承担。这样的优化结果对应指定夹持条件；自由块体则允许额外的剪切、伸缩或原子位移。

Ti–Zr–Hf–Ta 合金的计算给出了具体差别：图2上图把原子固定在理想晶格点，下图允许原子弛豫，hcp、ω 与 bcc 的相对能量及交叉组成随之改变。比较结构时，释放了哪些内坐标必须一致。不同磁序、局域轨道占据或交换关联近似，也可能得到不同的平衡构型；它们各自的力收敛后，能量排序仍需在同一计算定义下比较。

残余力很小，结构为什么仍会沿某些位移降能？

在驻点附近，微小位移 u 引起的能量变化可写成 ΔE ≈ ½u^THu，其中 H 是能量对坐标的二阶导数矩阵，即 Hessian。力消失约束一阶导数，局部极小值要求位移方向没有负曲率。若 H 存在负本征值，沿相应集体位移就能降低能量。整体平移等零模应与结构畸变区分；真正的零曲率软方向还可能受高阶项控制。

高对称结构可以因对称性而恰好具有零力，优化过程中若始终保留该对称性，破坏对称性的位移就不会出现。bcc 向 ω 的转变涉及不同原子面按特定相位协同移动；图3用黑色和红色位置标出两种排列，并区分理想转变、带缺陷排列及更大位移后的孪晶 bcc。随机扰动可以产生候选结构，振动本征矢则能指定待检验的位移方向。

周期晶体中，动力学矩阵 D(q) 的本征值给出 ω²。ω² ，常以负频率画出。Γ 点只检验与所用晶胞同周期的振动；更长周期的畸变需要相容超胞或非零 q 点。只沿少数高对称路径取点，还可能漏掉路径之外的软模。

Γ 点附近很小的负值可能受残余应力、力常数截断和声学求和规则影响，二维弯曲声学支还对旋转不变性敏感。确认物理失稳时，可检验负模随超胞和力精度是否保持，并沿振型畸变后重新弛豫。均匀剪切或伸缩则由弛豫离子弹性矩阵检验；有外部应力时，应采用相应载荷下的稳定条件。声子与弹性响应检查的扰动并不完全相同。

局部极小值为什么仍可能分解或相变？

一个构型抵抗邻近的小位移，其他原子排列仍可能具有更低能量。同组成多晶型可以直接比较匹配条件下的每原子能量；存在多种组成的竞争相时，最低能量组合必须保持总元素比例。在常见的 0 K 静态筛选中，候选相高于该组合的能量差记为 ΔE_hull，按同一原子数归一。

Al–Fe 的例子将相竞争画成了凸包。图4左图采用 1250 K 的热力学描述：金属间化合物使用 DFT 形成焓，AlFe 固溶体计入有限温混合自由能。某相相对于元素的形成能为负时，仍可高于其他化合物的组合。右图换成逆凸包能量坐标，稳定相可位于负侧；它与通常非负的 ΔE_hull 使用不同定义。

凸包比较描述分解的热力学驱动，真正发生相变还涉及断键、配位重排和新相成核。石墨到金刚石的候选路径中，形成两相界面与界面继续推进具有不同能垒。图5a比较 10 GPa 下的计算路径，b、c记录两个阶段的压力依赖，d、e给出含界面中间晶体的原子重排；终态能量较低的路径仍可具有较高峰顶。

局部优化沿当前能量梯度更新坐标，通常不会跨越这样的有限转变能垒去寻找另一个极小值。图5使用可变晶胞 NEB 在指定压力下计算路径，能量按每原子归一；实际成核速率还受临界核尺寸、界面能和缺陷影响。零温定压相变比较焓，有限温速率所用的活化自由能还包含振动和熵贡献。

有限温度下，优化构型能保持多久？

有限温度下比较相稳定，需要在相应条件下计算自由能。定温定压时，G 可包含静态电子能、振动自由能、电子热贡献与 pV，化学无序还会引入构型熵。非谐振动能够改变平均结构及有效回复力，高温平均位置接近高对称晶格时，瞬时原子仍可能偏离那些位置。

TiZrHf 的 54 原子 AIMD 模型中，图6的位置分布和平均位置给出了两种记录：300、900 K 的平均坐标沿特定方向偏移，1500 K 时接近 bcc 晶格点。这个尺寸的结果用于识别温度下的结构变化，定量转变温度还受超胞采样限制；单帧原子排列与时间平均晶格应分别描述。

AIMD 在轨迹各步计算电子结构产生的力，所观察的重排受温度、晶胞尺寸、系综和时长限制。温度和总能量的有界涨落可以与结构重排并存：恒温器持续调节动能，原子仍可能扩散、失去长程有序或改变配位。键长分布、平均位置、结构因子和配位随时间的变化，记录了能量曲线未直接标出的结构演化。

锡烯的大尺度模拟提供了时长差异的例子。图7采用经 DFT 数据训练的机器学习势，并非逐步求解电子态的 AIMD：200 K 轨迹从 0.3 ns 时的起伏单层，逐渐转为含 bcc 配位的三维结构；500 K 的万原子轨迹则在数百皮秒内形成液滴。对应温度和转变时长属于该势、尺寸及模拟条件。

短轨迹未出现重排时，稀有事件可能仍未被采样。用单一热激活事件作估算，平均等待时间 τ ≈ ν⁻¹exp[ΔG^‡/(k_BT)]。若假设尝试频率 ν = 10¹³ s⁻¹、活化自由能 0.8 eV、温度 300 K，τ 约为 3 s；在泊松事件模型下，一条 10 ps 轨迹采到该事件的概率仅约 4 × 10⁻¹²。