鸢算科技：SrVO3 频率依赖 cRPA 有效库仑相互作用
先看参考图：bash preview.sh。完整读取、原生报告、绘图：bash run.sh。
输出results/data、tools-output、plots/figures。使用新解压目录，输出目录受保护。
后处理：py4vasp0.11.3/Python3.11；NumPy2.3.5/Matplotlib3.11.2，已分配资源内1CPU即可，不运行VASP。
实际原生计算：VASP6.6.1、Intel oneAPI2024.2/IntelMPI、HDF51.14.6、Wannier903.1.0接口，MPI8、每进程1线程。此算例需要同时支持HDF5和Wannier90的可执行程序；不能用仅有HDF5的程序直接替代。
SrVO3立方5原子、a3.84652A、550eV、Gamma4^3、144能带、ENCUTGW200。复用自己的已收敛WAVECAR/WAVEDER/CHGCAR和真实WANPROJ，三个V t2g-like MLWF；不在响应步骤重新优化Wannier函数。
inputs保存四个外部实频0/5/10/20eV的完整INCAR/POSCAR/KPOINTS；自行准备匹配PAW、波函数、响应导数和Wannier投影文件。它们不在本包内。现有MLWF实际占用孤立能带21-23，不将本例说成纠缠带分数占据优势展示。
carriers为自有全分辨率选定张量NPZ，运行时创建明确标记的最小HDF读取载体，不是完整原始VASP输出。未附源代码、可执行文件、PAW或波函数/投影系数。
每次原生响应只输出一个外部频率，NOMEGA=8不能当成八个已导出的频率点。四个点来自四次真实响应计算，未宣称网格、空带、截断或频率收敛。
py4vasp单频effective_coulomb.read()未返回频率字段，脚本独立读取native comega实部/虚部。这里是实轴推迟响应，不是Matsubara频率。
关键版本约定：此VASP6.6.1驱动写cijkl=V_high-V_low；py4vasp0.11.3名为bare_low_cutoff的返回值在这里是截断修正量。真正低截断裸相互作用=bare_high_cutoff-bare_low_cutoff。屏蔽U已含修正，不要重复加。
API张量索引(i,l,k,j)对应原生文本I,J,K,L，转为native四指标使用api[field][0].transpose(0,3,2,1)。脚本与四个频率全部648个原生VIJKL/UIJKL打印条目逐一对照。
U=平均T[i,i,i,i]，u=平均T[i,i,j,j]且i!=j，J=平均T[i,j,j,i]且i!=j。静态平均U约3.277326eV，裸高截断约16.270731eV，修正约.306126eV。
10eV处Re U为负是本次动态响应原值；不能把它当静态U、删去负值或据此断言材料静态不稳定。虚部采用原生推迟响应符号。
绘图连接线仅引导视线，不是计算了中间频率。projector-metrics是原投影系数导出的Gram误差/迹/本征值等标量，未包含原投影系数。
