VASP Si有限带数RPA与EXX总能
新目录解压，bash preview.sh仅算参考数据的能量组合并重画2图；bash run.sh重新分析完整RPA原生数值与EXX分量并绘图。参考source为真实Si2 PBE/Exact/RPA计算，新增定轨道EXX一次，未重跑已有自洽/Exact/RPA。参考图figures，新提取results，新图plots；复跑换新目录。
严格区分：EXHF一行不等于完整EXX总能；完整EXX包括动能/Hartree/离子/PAW等。RPA输出的Hartree contr.toMP2也不是完整MP2。有限总能=完整EXX+RPA相关能+占据修正，本例占据修正0。
本例固定两原子Si原胞，550eV/EDIFF1e-6、Gamma3^3、64带、ENCUTGW150、NOMEGA16。64是有限虚轨道空间，不能视为所有平面波或收敛的材料能量。原生8个内部响应截断点与其converged value只属于内部外推，不能代替带数/波函数截断/k网格/频率收敛。
完整参考输入inputs。要从头计算，在包根目录准备合法匹配Si的POTCAR，再sbatch run.slurm。脚本做PBE→Exact→两个共享相同Exact波函数的EXX与RPA分支，并读取结果。不会自动取得POTCAR，已有work停止保护。原始波函数、完整电子输出和程序不随包分发。
从自己的已完成目录提取，module load pymatgen/2026.9.24，然后python read_rpa.py work/si-rpa --output my-rpa，python read_exx.py work/si-exx-fixed-PBE --output my-exx，cp my-exx/EXX-DATA.json my-rpa/，python analyze.py --data my-rpa --output my-results，python plot.py --data my-results --output my-plots。
该analyze脚本针对当前2原子/64带/8内部截断和4项q输出示例，换模型要同步修改检查范围。不同算法GWR/RPAR与不同版本输出不能直接假设同格式。新DFT提交和只重画参考数据的命令分开使用。
