说明：本文主要介绍结构优化在第一性原理计算中真正处理的对象，包括总能量、原子受力、晶胞应力、收敛判据和后续性质计算中的使用边界。
基本概念：结构优化在最小化什么？
在 DFT 语境里，结构优化处理的不是一张图片中“更顺眼”的几何外形，而是在给定模型、泛函、赝势、自旋态和边界条件下，寻找 Born-Oppenheimer 势能面上的一个平衡构型。可变的量通常写成原子坐标 R 和晶胞矩阵 h，程序在每一个几何点求解电子结构，得到总能量 E(R,h)。优化目标首先是这个几何点对应的总能量函数，而不是带隙、吸附能、磁矩或某个实验峰位。
总能量给出高度，原子受力给出斜率。对第 i 个原子，力表示为 Fi = -∂E/∂Ri；允许晶胞变化时，应力张量来自能量对晶胞形变的响应。结构优化真正追踪的是能量、力、应力和位移之间的一致收敛，其中任何一个量失控，几何结果都可能停在不合适的位置。
模型自由度先限定了“能优化什么”。固定晶胞的表面吸附模型只移动部分原子，变晶胞体相计算还会改变晶格常数和夹角，外加压力时常对应焓 H = E + pV 的平衡。冻结底层 slab、固定分子构型、保持对称性或限定磁态，都会把搜索空间收窄。约束不是装饰参数，而是优化对象的一部分，它决定程序能否离开初始构型附近的某些方向。
这里的低能不是脱离模型的绝对低能。不同泛函、U 值、色散修正、赝势和磁态都会改变 E(R,h) 的形状，优化坐标只沿当前能量面移动，同一初始结构在不同能量面上可能到达不同局部平衡。
用一句话概括，结构优化是在指定自由度内寻找力和应力接近平衡的局部极小点。若自由度没有定义好，“优化结果”对应的势能面和变量集合就不清楚，例如同一个 slab 可以只放松吸附物，也可以放松表面两层原子，还可以重新比较不同终止面。
图1. 不同切面、表面结构和化学势区域共同限定可被优化的表面模型。DOI：10.1039/D3NR06468A
结构优化中残余力有何含义？
只优化原子位置时，晶胞像固定坐标系，原子在其中寻找力接近平衡的位置。初始结构可能来自实验 CIF、手工构型、结构搜索、吸附位点枚举或上一轮计算；程序不预设某个键长的目标值，只根据当前电子密度计算力，再更新坐标。残余力小表示局部能量面的一阶斜率接近零，并不自动表示构型是全局最低点。
局部几何的变化常体现为键长、键角、配位环境和吸附高度的重新分配。沸石数据集中，优化前后的 Si-O 键长和 Si-O-Si 角度分布出现明显位移，说明原子坐标弛豫会改变局域四面体连接方式。坐标优化改变的是原子间相对位置，电子结构图谱和能量差只是后续从该几何基准上计算出来的结果。
图2. 沸石结构优化前后 Si-O 键长、Si-Si 距离和 Si-O-Si 角度分布发生改变。DOI：10.1038/s41597-022-01160-5
在离子步序列里，能量下降、最大力变化和原子位移会以不同速度靠近阈值。能量可以变化很小，而局部力仍然偏大；相反，大体系中某个原子的残余力异常也可能被平均能量掩盖。最大力、均方根力、位移和能量变化共同构成坐标弛豫判据，最后一个总能值不能覆盖全部几何状态。
对于吸附、缺陷和界面模型，残余力还会指向局部配位张力。吸附物转动、表面层起伏、缺陷邻近原子外移，都会改变后续 Eads、Bader 电荷或 PDOS。结构优化没有直接优化这些派生量，它只把几何推向当前能量模型允许的局部平衡。
若一个吸附构型残余力没有降下来，后续电荷差分图中的电子积累和耗散区域就可能混入几何未平衡的影响。残余力本身不是性能描述符，它是检查几何状态是否适合作为后续电子结构计算起点的门槛。
图3. 优化轨迹中的能量差、原子梯度、相对体积和应力张量分布。DOI：10.1038/s41597-022-01160-5
晶胞优化为何需考量应力与体积？
周期材料的晶胞不是盛放原子的静态盒子。体相、二维层状材料、分子晶体和高压结构中，晶格矢量决定原子周期重复方式，也决定倒易空间采样、体积、密度和应变状态。当晶胞自由度打开时，优化变量从原子坐标扩展到晶格长度、夹角和体积，应力张量便成为与力同等级的几何响应量。
在晶胞被允许改变的计算中，固定晶胞和变晶胞对应两种平衡条件。固定晶胞下，程序只让内部坐标适应给定晶格；变晶胞下，晶格会朝目标压力或目标应力状态移动。分子有机晶体的单胞优化显示，PBE-D3 和 B97-D 下不同体系的单胞体积可以增加，也可以减小。晶胞优化回答的是给定计算模型下的平衡体积和形状，不是把实验室温结构逐点复制出来。
图4. X23 分子晶体在 PBE-D3 和 B97-D 下优化晶胞参数后出现不同幅度的单胞体积变化。DOI：10.3390/cryst9120665
应力没有收敛时，晶胞通常仍带着压缩、拉伸或剪切倾向。二维材料若面内晶格没有放松，能带和声子会带有预应变特征；层状材料若缺少合适色散描述，层间距可能偏离合理范围；磁性材料若磁态设定改变，平衡体积也可能随之改变。晶胞、电子态和磁态之间存在能量耦合，结构优化读数必须回到同一个计算条件下理解。
当截断能、k 点或赝势设置改变时，晶胞优化会把数值误差写进应力和体积。截断能不足会带来基组相关应力误差，k 点不足会扰动金属体系总能曲面，赝势价电子设置会改变平衡体积。这里的判断边界不是“参数越大越好”，而是目标能量差、残余应力和关键结构量在可接受范围内稳定。晶胞结果的物理含义由计算目标限定，例如体相形成能、二维应变、表面 slab 或分子晶体 packing 对晶格自由度的要求并不相同。
对于外加压力或应变题目，目标函数的物理含义会再次改变。零压体相优化看平衡体积，高压结构看给定压力下的焓，二维应变计算则常固定面内晶格再放松内部坐标。同样叫结构优化，优化自由度不同，所得几何含义就不同。
优化结构可直接作为计算结论吗？
优化结束后得到的是一个在指定条件下满足收敛阈值的几何基准。沸石数据集中，最后一步的最大核梯度、能量变化和最大原子位移都被统计出来，说明“收敛结构”本身仍有清楚的判据。残余力、能量变化、位移和必要的应力标准共同限定几何质量，单写“结构已经优化”很难判断几何状态。
图5. 沸石数据集中最后一个优化步的最大核梯度、能量变化和最大原子位移分布。DOI：10.1038/s41597-022-01160-5
优化结构不是全局稳定性的同义词。局部极小点可能来自初始构型、约束、对称性、磁态或覆盖度选择；声子虚频、AIMD 轨迹、形成能、凸包位置和反应自由能会继续检验不同物理层面的稳定与可用范围。结构优化给出几何平衡点，后续计算检验热力学、动力学和电子结构问题，两类结果不能互相替代。
在后续性质计算中，优化结构相当于共同坐标基准。静态总能、DOS、能带、功函数、吸附能、NEB 路径和声子谱都从某个几何状态出发；几何若没有达到相应收敛，性质差异会混入结构噪声。单一指标不能把优化质量包圆，例如晶格常数接近实验并不能保证带隙准确，残余力很小也不保证吸附构型已覆盖全部候选位点。几何收敛只说明计算起点足够清洁，后续性质仍要由对应物理量来判定。
优化结果的物理含义来自具体模型：体相对应是否变晶胞和目标压力，表面对应冻结层与真空层，吸附对应初始位点和覆盖度，缺陷对应电荷态和超胞尺寸，磁性体系对应磁序。结构优化到底优化的是给定自由度上的能量平衡；它为后续计算提供几何起点，也保留了方法、模型和约束带来的边界。