在上一章《第十六章：INCAR文件常见错误分析与解决方案！——2026新版VASP基础教程》中，华算科技朱老师详细介绍了VASP示例任务。VASP计算一般从结构优化开始，然后对优化后的结构进行性质计算。本章将介绍VASP结构优化计算实例，让大家更好的了解计算流程，具体包括结构优化、半导体性质案例、表面与催化案例。
VASP晶体结构优化的根本目标是寻找体系总能量（E）的全局或局部极小值，即寻找原子在晶格中的“物理平衡位置”。
这一过程基于密度泛函理论（DFT）：
给定原子位置后，VASP通过自洽迭代求解Kohn-Sham方程，计算出体系的电子能量。
基于Hellmann-Feynman定理，计算每个原子受到的力（F = -∂E/∂R）。
根据计算得到的力，调整原子位置和/或晶格参数，使得能量降低。
不断重复上述过程，直到力和能量的变化都满足预设的收敛标准（如EDIFFG）。
结构优化不仅仅是移动原子，还涉及到严格的物理量控制：
最关键的收敛指标。通常要求每个原子的残余力小于0.01 eV/Å（严格计算可能要求更小，如0.001 eV/Å）。
每一步的能量变化必须足够小，通常要求能量差小于10^-5 eV。
对于需要优化晶胞形状和体积的计算，还必须使晶胞内部的应力张量趋近于零。
VASP的结构优化并非单一过程，而是可以细分为不同的维度，分别针对不同的物理需求：
|
维度 |
功能 |
关键参数 (INCAR) |
适用场景 |
|
原子位置优化 (Ionic Relaxation) |
固定晶格参数，仅移动原子坐标。 |
IBRION=2 (共轭梯度法), NSW, EDIFFG |
分子动力学模拟、低温结构弛豫 |
|
晶胞体积优化 (Cell Volume Optimization) |
调整晶胞体积，保持形状不变。 |
ISIF=2 |
计算晶格常数、体积模量 |
|
全结构优化 (Full Relaxation) |
同时优化晶胞形状、体积和原子位置。 |
ISIF=3 |
完全弛豫、寻找稳定相 |
注意：全结构优化往往是最消耗计算资源的步骤，尤其是对于大体系（>100个原子），计算成本会呈指数级增长。
基于上述原理，标准的VASP结构优化流程通常包括以下步骤：
来源：可以是实验测得的晶体结构（如CIF文件），也可以是通过材料基因组学（Materials Project）等数据库获取的预测结构。
注意：需要确保原子之间没有不合理的重叠或过大的距离（常见的坑是初始结构离得太近导致力过大）。
POSCAR：晶胞信息和原子坐标。
POTCAR：赝势文件（PAW势）。
KPOINTS：k点网格（通常需要足够稠密，尤其是金属体系）。
INCAR：控制参数文件（最关键）：
ENCUT：平面波截断能，通常设置为POTCAR推荐值的1.3倍以上。
PREC：精度设置，通常选Accurate。
EDIFF：电子收敛精度。
EDIFFG：离子力收敛阈值。
IBRION：离子弛豫算法。
ISIF：控制优化维度。
资源配置：对于较大的晶胞，通常需要数十甚至上百个CPU核心，并行计算。
计算时间：可能从几小时到几天不等，取决于体系大小和优化维度。
OUTCAR：核心输出文件。通过检查其中的FORCES段落，确认所有原子的力是否都小于设定阈值。
CONTCAR：优化后的结构文件。若计算收敛，即可使用该文件进行后续计算。
对称性检查：VASP在优化过程中可能会破坏原有的晶体对称性（如降至P1），因此优化后通常需要使用软件（如VESTA、VaspView）重新识别对称性。
稳定性分析：对于新发现的结构，常需计算声子谱（Phonon）或弹性常数，验证是否满足Born准则。
结构优化原理：介绍了电子结构求解、力的计算，结构更新
收敛标准：介绍了力和能量收敛
结构优化流程：介绍了结构准备、输入文件准备、提交计算
下一章将正式引入本次教程的核心—MgO晶体优化。我们将从前期准备工作、电子结输入文件设置、任务提交与监控方面详细介绍VASP结构优化案例，敬请期待！
