在上一章《第二十七章VASP计算实战：电压曲线计算流程讲解！｜2026新版VASP基础教程》中，华算科技朱老师详细介绍了VASP示例任务。VASP计算一般从结构优化开始，然后对优化后的结构进行性质计算。本章将介绍VASP电池材料计算实例，让大家更好的了解计算流程，具体包括结构优化、半导体性质案例、表面与催化案例。
离子嵌入能（Insertion Energy）和脱出能（Extraction Energy）本质上是一个热力学过程的能量差。对于一个通用的离子 M 嵌入到宿主晶格A
中形成固溶体A-M，其平均嵌入能
的定义为：
其中：
为包含n个嵌入离子的体系总能量（已经几何优化的）；
为纯宿主晶格的总能量；
为离子M 在参考态（通常是金属相或气相）的化学势（通常对应于其体相或分子相的总能量）；
n为插入的离子数量。
如果计算的是离子的脱出能，则其物理意义与嵌入能相同，只是符号相反（取决于定义方式），通常关注的是正值（脱出困难）或负值（易于脱出）。
VASP 计算离子嵌入脱出能的核心流程可以概括为体系构建 → 参数设定 → 几何优化 → 总能量提取。为了确保结果的可比性，必须保证所有计算（含有插入离子、纯宿主以及参考态）的计算条件（K点网格、平面波截断能、赝势等）一致。
以下是该流程的具体步骤：
宿主晶格 (Host Lattice)
首先获取宿主晶格的晶体结构（如 LiFePO
、石墨烯、金属氧化物等）。可以通过 Materials Project、ICSD 等数据库下载 POSCAR 文件，或使用 Materials Studio / VESTA 构建。
如果研究的是离子在晶格间隙位的嵌入，需要先对原始晶格进行体积扩展或超胞构建（如 2x2x2），以容纳插入离子并避免周期性镜像之间的相互作用。
插入离子 (Inserted Ion)
在宿主晶格中选择合适的嵌入位置（如八面体空位、四面体空位或层间空位）。将离子原子添加到 POSCAR 文件中相应的坐标位置。
如果是脱出能的计算，只需要保留宿主晶格结构，不需要进行任何改变。
参考态 (Reference State)
计算
时，需要构建离子M的参考态体系。例如，对于金属 Li，需要构建金属 Li 的体相晶格（通常为 BCC 结构），对于气体
，则需要构建 分子的超胞（如15×15×15 Å 的盒子），并进行单点能计算。
VASP 的计算完全由输入文件驱动。针对离子嵌入能计算，通常需要准备以下文件：
POSCAR
包含原子种类、晶格矢量和原子坐标（已插入离子的坐标）。
POTCAR
赝势文件。对于体系中涉及的每一种元素，都需要对应的 PAW 赝势。需要注意的是，离子
的赝势必须与参考态计算时使用的赝势一致。
KPOINTS
K-点网格。对于能量差分计算（如形成能），通常采用高精度的K-点网格。例如，对于 2x2x2 超胞，可能需要使用 
或更密集的网格，以确保总能量收敛。
INCAR
控制参数文件。对于离子嵌入能计算，关键参数包括：
PREC = Accurate：高精度计算。
ENCUT = 520 eV（或更高）：平面波截断能。
ISMEAR = 0 (Gaussian) 和 SIGMA = 0.05：对于绝缘体或半导体。
ISIF = 2：仅弛豫离子位置，不优化晶格（如果需要考虑晶格膨胀效应，则设为 3）。
IBRION = 2 (Conjugate Gradient) 或 IBRION = 1 (RMM-DIIS)：离子弛豫算法。
NSW = 200：最大离子步数。
EDIFF = 1e-5：电子收敛标准。
EDIFFG = -0.02：离子力收敛标准（0.02 eV/Å）。
ISPIN = 2：如果涉及磁性离子（如 Fe、Co 等），需要设置自旋极化计算。
提交脚本 (Submission Script)
根据 HPC 环境（如 SLURM、PBS）编写提交脚本，调用 VASP 可执行文件（如 vasp_std）。
几何优化 (Geometry Optimization)
计算的第一步是对含有插入离子的体系进行离子弛豫（ISIF=2），使体系达到力学平衡。此过程涉及 VASP 内部的自洽循环 (SCF) 和离子步进循环。
SCF 循环：读取 POSCAR、POTCAR 和 KPOINTS，猜测初始电荷密度，解 Kohn-Sham 方程，计算总能量 E，并判断电子密度是否收敛
。
离子步进：如果电子密度收敛，计算每个离子上的力F；如果F 未满足
，则根据F 更新离子坐标，生成新的 POSCAR（即 CONTCAR），进入下一轮 SCF 循环。
循环直至所有离子力低于阈值，记录最终的总能量
。
纯宿主体系计算 (Host Lattice Calculation)
对纯宿主体系（未插入离子）进行相同的几何优化和能量计算，获得
。
参考态计算 (Reference State Calculation)
对离子M 的参考态（如金属相或气相）进行能量计算，获得
。如果是气相分子（如
），需要在 ISIF=2（固定盒子）下进行单点能计算，确保盒子足够大以避免相互作用。
计算完成后，VASP 会生成一系列输出文件：
OUTCAR：详细的计算日志（包括能量、力、收敛信息）。
OSZICAR：简洁的能量收敛记录。
CONTCAR：最终的结构文件。
CHGCAR：电荷密度文件。
WAVECAR：波函数文件。
提取总能量：
通常使用 grep 或脚本提取 OUTCAR 中的 free energy 或 energy(sigma->0)。
确保提取的能量是静态能量（SCF 收敛后），而非弛豫过程中的瞬时能量。
计算嵌入能：
将提取的
、 和
和 代入公式，计算出
。
理论基础：介绍嵌入能与脱出能定义
计算流程：介绍了构建模型、准备输入文件、数据提取
下一章将正式引入本次教程的核心—NEB方法简介。我们将从NEB方法物理背景、NEB方法原理、VASP计算流程方面详细介绍VASP 过渡态与反应路径，敬请期待！
