鸢算科技 VASP性质教程第四组，2026-10-25

六篇：分子电离与电子亲和能、Cu(111)吸附位点、H吸附自由能、DOS与ELF、DFT+U、色散修正。
NiO的U扫描保持同一结构和相反初始Ni磁矩；双层石墨烯是固定堆垛刚性距离扫描，各方法以自身4.4Å能量为参照，不是无限分离结合能。
先读官网教程，再选择对应目录，不要把所有任务无区别一起提交。
准备自己的合法POTCAR，元素顺序按POSCAR；执行sbatch run.sh。

1. 分子：中性水参考结构已固定，两个电荷态使用完全相同内部几何。改变盒子时不要缩放分子键长。水NELECT=8/7；F=7/8。Γ点分子用vasp_gam。若F18使用restart目录，需先运行precondition目录获得CHGCAR；不要遗漏这个依赖。
2. 表面：先运行四个*-preconditioned目录，固定底部两层，放松自由原子。将各自CONTCAR复制为最终偶极修正目录的POSCAR、CHGCAR复制为CHGCAR，再做静态。参考静态POSCAR已附，是本次收敛结构。clean/fcc用dipole-centered，hcp/ontop用dipole目录。不同位点保持相同晶胞、k点、参数和能量口径。
3. H2：hydrogen-reference得到优化结构和能量；h2-vibrations以其CONTCAR为POSCAR做振动。fcc六个位移从同一收敛fcc结构对H施加±0.015 Å；Cu坐标固定。本包中这六份POSCAR已完成位移，重新优化后须基于新结构重新生成。
4. 自由能：加载module load ase/3.29.0或使用含ASE/Numpy的Python，运行python extract.py。脚本只读取各自计算输出，不需要复制私有result.json。H-only部分Hessian忽略基底振动变化；气相H2用298.15K、1bar理想气体。
5. DOS与ELF：si-electronic产生DOSCAR/CHGCAR/ELFCAR，cu-pdos产生轨道投影。DOS-DATA.json是本次参考曲线；DENSITY.json是本次Si最近邻键的线段数据。铜d带中心包含全部五个d轨道，窗口[-10,0]eV相对Ef。
6. reference图和JSON是本次结果，不会随自己的VASP作业自动变化。extract.py用于分子/吸附/热化学的重算；DOS/ELF可按教程用VASPKIT或其他读取工具重新提取。不要把固定几何的离子态大原子力当成垂直IP计算失败，也不要把垂直IP称为绝热IP。

不含POTCAR、WAVECAR、CHGCAR、原始私有日志。输入含相应的科学近似和条件，不能作为任意体系通用收敛参数。
