鸢算科技 VASP性质教程第二组（2026-10-25）

对应官网八篇：晶体声子、准谐热力学、非谐热导、二维力学、中性空位、AIMD、NEB、电荷与COHP。
calculations 下有实际INCAR/POSCAR/KPOINTS及参考结果；每目录自行准备合法POTCAR后 sbatch run.sh。
VASP 6.6.1/PBE PAW64/550eV/EDIFF=1E-6；优化力0.02eV/A。默认P100单卡，按账号分区修改模板。
phonon 目录是位移结构，不得再优化。三阶位移编号保留被截断项的空缺，不能自行重编号。
生成全部vasprun.xml后，用配套Phonopy4.8.3、Phono3py4.8.2与ASE运行两个分析脚本。分析将写新JSON；先备份本包参考JSON。较重热导后处理也应提交Slurm。
NEB需进入neb，提供H的POTCAR，加载VASP及ASE Python环境后提交run.sh；workflow.py从端点开始优化。最终图像文件与JSON是本次参考结果。不是把五张图像独立静态计算就等于重跑NEB。
LOBSTER使用随包lobsterin，先完成VASP生成WAVECAR，再加载LOBSTER6.0.0运行lobster；COHP文件为参考数据。
Bader需重算AECCAR0/AECCAR2/CHGCAR，再按官网chgsum.pl与bader命令处理；reference-ACF.dat是本次参考值。
POTCAR、WAVECAR和密度文件不分发；原始私有作业日志不在包内。
热导网格尚未收敛；空位为16位点小超胞；MD只有0.25ps，不用于报告可靠扩散系数。详见对应教程。

分析脚本附带本集群ASE模块依赖路径，不安装任何新软件。完成所有对应VASP计算后可在解压根目录 sbatch post-phonons.sh。
