VASP五版本Libxc实算与终端操作

在新目录解压，bash preview.sh重画3张参考图；bash run.sh重新提取80条件全数值并画图。只进行1核后处理，不提交DFT。新结果results，图plots；reference/fields是完整能量、坐标、力、应力、晶胞与元素数组；inputs按版本/条件放完整INCAR/POSCAR/KPOINTS。重新运行换新目录。

原生计算为2026-10-20真实CPU8173M，VASP6.6.1/6.5.1/6.4.3/6.4.2/6.3.2，Libxc7.1.2，oneAPI2024.2/IntelMPI2021.13/MKL2024.2/HDF51.14.6。后处理pymatgen/2026.9.24。6.4.2/6.3.2采用已部署的非正密度保护分支。普通与Libxc对照使用同版本的同一个Libxc构建，分别使用内置泛函标签或GGA/METAGGA=LIBXC标签。不要把它们当GPU实现。

每版本12个主条件：Si PBE/SCAN/r2SCAN/PBE0各内置与Libxc，H2O PBE一对，Fe自旋PBE一对；另4个水分子高截断/位移对照，总80个真实静态结果。结构固定并非优化终态。使用EDIFF1e-10是为辨别微小实现差；一般参考计算可从1e-6开始并检查目标性质。水的绝对能量差保留，不做拟合平移，也不把绝对能量跨版本混合。

要计算自己的输入，复制一个目录成my-Si，准备合法匹配的POTCAR，进入该目录执行yskit，选择对应版本CPU Libxc分支生成Slurm脚本并提交。在可用的相应CPU资源内运行；选择不同GPU分区不代表这份CPU程序能GPU加速。
结束后回包根目录，module load pymatgen/2026.9.24，然后python read_vasp_result.py my-Si --output my-result。输出原生XML的完整坐标/力/应力/自由能；不把附件参考值覆盖成自己的新结果。

原始完整XML、OUTCAR、HDF、POTCAR和波函数均不随包分发；新计算由自己生成。电子自由能为e_fr_energy/TOTEN，比较时保持同一物理定义、几何、版本和泛函实现。
