CP-VASP 电势参照、电子数扫描与氢吸附电子巨正则能差

bash run.sh：仅解析和复画11份已完成的真实集群计算的数值摘录，不运行VASP、不提交作业。
reference/OUTCAR.selected.txt 和 stdout.log.selected.txt 的每行均从实际原始输出逐字选取；不是完整OUTCAR，不可用于接续VASP。inputs/ 提供完整的INCAR/POSCAR/KPOINTS。没有POTCAR、WAVECAR或原始XML。

新计算：复制 inputs/ 到自己的 work/，按各目录元素顺序准备有权使用的PAW_PBE Cu/H POTCAR；保持Cu的ZVAL=11、H=1。先完成cu-neutral，再用自己的中性参照生成目标值。使用已配置的CP-VASP CPU Slurm模板，在每个目录提交（需要CPU分区可用）。保存原始OUTCAR，并把原生标准输出保存为stdout.log。裸电极和两个吸附位点必须使用完全相同的三个TARGETMU。H2目录为独立气相分子优化，LSOL/LCEP均关闭。
全部11个计算完成后：python analyze.py --results work --output own-analysis；python plot.py --analysis own-analysis --structures work --output own-figures。

N0取POTCAR中性价电子之和，示例Cu4=44、Cu4H=45。GCE=E_sigma0-mu*(N-N0)。吸附电子巨正则能差=GCE(Cu4H)-GCE(Cu4)-0.5E(H2)，只比较同TARGETMU。
旧Cu4与Cu4H保持理想固定几何，单个H仅在上表面；H2单独优化。裸Cu4电容按两个等价表面归一，单侧吸附模型不套用双面均分。电子能量不含零点能、熵、绝对SHE/RHE标定或收敛研究，不能标成完整吸附自由能。

VASP6.3.2 CP-VASP v2 + VASPsol++专用CPU构建；当前模块依赖oneapi/2024.2（运行时）；原生程序的具体编译来源以集群软件版本说明为准。新H2在相同程序、ENCUT550eV、EDIFF1E-6、EDIFFG=-.02eV/A下优化。旧10份Cu/CuH数据复核复用，没有重计为新DFT。
