结构优化

本文总结了Abinit结构优化的计算流程，以及代码的实现方法。

在ABINIT软件中，当设置计算参数 ionmov=2 时，程序采用 Broyden-Fletcher-Goldfarb-Shanno (BFGS) 的拟牛顿法来进行晶体结构和原子位置的几何优化。BFGS方法通过近似目标函数（系统总能量）的逆海森矩阵来引导原子坐标朝着能量极小值收敛，在计算效率与优化收敛性之间取得了极佳的平衡。

在ABINIT的结构优化过程中，参数 ionmov=22 对应于有限内存Broyden-Fletcher-Goldfarb-Shanno (Limited-memory BFGS, L-BFGS) 方法。该方法是一种著名的拟牛顿法（Quasi-Newton Method），不需要显式计算并存储庞大的Hessian矩阵及其逆矩阵，而是通过保存最近几步的迭代历史信息来近似逆Hessian矩阵，特别适用于拥有大量自由度（如大体系多原子结构的弛豫和晶胞优化）的问题。

结构优化（structure optimization 或 structure relaxation）是第一性原理计算中最基础、也最容易"算了很久却不一定算对"的环节。它的目标是通过调整原子坐标、晶格形状或晶胞体积，使体系到达一个能量更低、受力更小的稳定构型。对于 VASP 而言，优化的自由度主要由 ISIF 控制，优化算法主要由 IBRION 控制，而是否收敛通常由 EDIFFG 等参数判断。