VASPsol++ 实际小体系教程：输入、数值与9张图

先看现成结果：bash preview.sh
图在 preview/figures；原数值在 reference/，完整输入在 examples/。
本包只含自有结构/输入/数值，不含POTCAR、WAVECAR或VASP程序。

实际计算：选择一个examples内的算例，复制到自己的工作目录；按POSCAR顺序准备有授权的POTCAR，检查INCAR的ISTART/ICHARG是否需要自己已有的匹配WAVECAR。若没有续算波函数，使用ISTART=0、ICHARG=2开始新的计算。
把submit.sh复制到工作目录，确认CPU分区可用后sbatch submit.sh。
模板为CPU 56 MPI进程的正式资源入口；本教程不包含该节点的性能测试。计算大任务使用Slurm，勿在登录节点运行生产DFT。
同一例的气相/溶剂、clean/ads、H2参照应保持版本、势文件与能量设置一致。

已有全部原生输出按examples结构组织时：
bash prepare.sh ./raw
bash extract.sh
提取读取自己的OUTCAR、OSZICAR、vasprun.xml、vasp.log和相关场文件，只建立链接，不重算电子结构。
结果 results/data/*.csv 与 results/figures/*.png。

H2O/H3O+为14A周期盒、450eV、EDIFF1E-6；气相优化力0.02eV/A，溶剂比较沿用同一气相几何，不能视作溶剂中重新优化的最低能量。
Cu8两层表面及顶位H是有限教学模型。两化学势没有标定到SHE。H2参照的巨势能差不含ZPE、熵和pH修正，不是完整反应Gibbs自由能。
ISOL1与2默认空腔参数不同。ISOL1的PHI在此构建为溶剂势修正，ISOL2为总势；请勿按相同定义直接比较。
LSOL_SCF=false且不写势场时只得到气相SCF；需要LVHAR=true触发固定密度溶剂后处理。此时原生E0仍是气相能，单列A_corr而不冒称自洽溶液总能。
q_sol取价电子密度零频分量减Z，在有限FFT网格与名义NELECT-Z有小差别；本例以RHOION积分独立复核。
