VASP RPA实操：相关能与总能怎样读

完成四阶段输入与提交后，最容易混淆的是“OUTCAR里哪一个数才是我要的能量”。本篇用同一个两原子Si原胞，把提取、相加、输出和图的含义逐项说明。

1. 先重做附件中的数据分析#

下载Si RPA与精确交换实算包：

python3 -m zipfile -e vasp-Si-RPA-EXX-wave305-*.zip .
cd vasp-Si-RPA-EXX-wave305
bash run.sh
ls results plots

results/DATA.json和CSV保留全部8个内部截断点、4项q贡献、实际晶胞/坐标及精确交换能量组成。不是从截图读数，也没有重新拟合修改参考能量。

2. 从自己的原生输出读取#

下面在资料包根目录执行，work/是上一教程的实际新计算目录：

module purge
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
cat my-exx/EXX-DATA.json
cat my-results/DATA.json

read_rpa.py读实际RPA OUTCAR和结构，read_exx.py读完整HF能量和占据修正，analyze.py组合两者。该分析器的检查范围对应本例2个原子、64条带、8个内部截断和4项q贡献；换成其他体系要同步调整范围，不能把不同输出强塞成相同形状。

3. 八行截断表不是八次独立计算#

本次一次ACFDT输出中的完整表如下，能量均为eV/两原子原胞：

这些行是同一次计算内部生成的不同响应截断，不是八个新的独立SCF。原生输出的内部外推相关能为−11.7029328497 eV/原胞。它反映当前有限带数与k网格下的内部响应基组外推，不能代替NBANDS、ENCUT、k点或频率收敛。详见VASP ENCUTGW说明。

右图保留四项原生q贡献，实部分别为−0.999352、−3.301890、−2.454368、−4.745574 eV，虚部均为0。图按输出出现次序编号，不给它们编造q坐标或均等权重；求和与150 eV相关能在原生打印精度内一致。

4. 怎样组合完整RPA总能#

本例采用相同固定PBE轨道，按统一的能量约定相加：

E(RPA) = E(EXX，完整HF总能) + Ec(RPA) + 占据修正
        = −18.52251733 − 11.5011832588 + 0
        = −30.0237005888 eV/两原子原胞

本例完整HF自由能与零展宽能量相同，占据修正为0；其他体系尤其金属必须检查自己的输出，不能把这个零当成通用规则。这里的EXX在固定PBE轨道上评价，不是另一套自洽HF参考。

OUTCAR中的EXHF=+58.03164388 eV是能量分解中的一项，不是完整的−18.52251733 eV HF总能。 完整总能还包括动能、Hartree、离子、PAW等组成。只拿EXHF那一行与RPA相关能相加会得到错误含义的数字。

若仅把相关能换成这次内部外推值，同一参考下的组合为−30.2254501797 eV/原胞。它仍然继承有限64带、3³网格和16频率点的限制；数值小数位多不表示材料精度达到那些位数。

5. Hartree-only MP2不是完整MP2#

RPA表里标成Hartree contr. to MP2的数字只有相应直接/Hartree部分，不包含完整MP2所需的交换项。不能将表中的−26.2535419734 eV写成“本次完整MP2相关能”。要学习直接项、交换项和总相关能的提取，使用LiH MP2三阶段教程。

6. 做相对能量时怎样保持参照#

比较两种结构或吸附前后时，每个状态都需要自己的相容PBE轨道、完整EXX和RPA相关能；保持赝势、带数策略、k采样和截断精度一致，并对最终的能量差检查收敛。不要用不同体系中一个单独相关能的绝对高低判断稳定性。

本例没有计算孤立原子参照，因此不输出Si结合能；固定晶胞也不构成RPA平衡晶格常数或结构优化。需要精度结论时，先增加虚轨道、响应与波函数截断、频率和k网格，再比较目标性质的变化。