qvasp 2.26：输入准备与真实 VASP 结果后处理

本包使用鸢算科技集群已经完成的 Si、Cu、Fe 计算，以及本次实际生成的 Cu/Fe/石墨烯结构。没有把文件准备或数据提取当成新的 VASP 计算。

1. 使用有软件使用权限的普通账号，解压后进入目录。
2. bash run.sh
3. 查看 work/cases 各例、work/analysis-core、work/analysis-extra 和 work/figures。

环境：module load ase/3.29.0 qvasp/2.26；Python、NumPy、Spglib、Matplotlib 来自 ASE 环境。qvasp 为官方 Linux 预编译 2.26；不推测其编译器版本。脚本加载模块并使用单线程。后处理不启动 VASP；实际生产 VASP 计算仍应提交 Slurm。

POTCAR、WAVECAR、完整原始 OUTCAR/XML 不在下载包内。POTCAR 在集群按所选元素现场生成，用户仍需有相应软件使用权限。reference/TOOL-PRACTICE-METADATA.json 只含真实计算中选取的参数、最后数值行和最终力。run.py 用它生成明确标注为 DERIVED TOOL-PRACTICE 的精简 OUTCAR，供工具练习；这不是完整 VASP 输出，也不能据此证明任意其他计算收敛。处理自己的任务时请使用自己的完整原始文件。

60 个目录包括参数对照和格式检查，不代表 60 种独立科学方法。15 类逻辑工具操作包括模板、势、网格、线模式、HSE路径、声子路径准备、格式、超胞、固定原子、原胞/对称、能量、能带、DOS/LDOS、Bader和功函数。

重要文件约定：
- qvasp -e 在本例读取 energy(sigma->0)，不是 TOTEN。五个应变、固定/放松原子两组分别比较，不把不同计算设置混排。
- c2p/c2p2 依赖 CIF 列顺序。包内提供实际通过元素、原子数、晶胞检查的 CIF。原生 p2c 在本例写出无效空间群 T；分析脚本显式改为 P 1 后按完整原子表读入，原文件保留。真实空间群另外用 Spglib 判断，不把 P1 当作材料对称性结论。
- 原生 -sw 在当前版本实际调用 CIF 导出，不进行轴交换；swap_axes.py 是明确的 ASE 替代操作。
- 原生 -3k/-k 有取整与自动维数差异；保留 .03 对照，Cu/Fe 的 .02 正例满足本例 .03 的最大间隔目标。二维 graphene-explicit2D 明确只把真空 z 改为1，面内16x16。
- khse-si-explicit11 重新生成 11x11x11 的不可约网格，整数权重之和1331，并加200个零权重路径点。它不会自动继承原 IBZKPT，不证明现有 WAVECAR 可继续使用。
- 连续且不重复端点的101点 Si 路径，原生 bandd 少输出3点。bands-Si-endpoint-format 只把已计算的三个端点重复一次，形成104条格式输入，未增加采样或插值；原生输出保留全部101个原始样点。SOURCE-POINT-MAPPING.json 给出逐行映射。此适配只对应本包这条路径，其他结构不要照抄索引。
- 半导体 band.dat 以 VBM=0 输出；分析 CSV 转换到实际 SCF 费米能5.8214eV。DOS 则以各自 OUTCAR 费米能为零。
- DOS 自旋向下列为镜像负值，不是负态密度。LDOS 最后一列是全部显示列的代数和，磁性体系不要把它当成物理的上下自旋总态密度。
- vline.dat 保存未减费米能的平面平均势，原生坐标从一个网格步长起。功函数本例为最大平面平均势减费米能；正式表面研究仍应检查真空平台、层数、真空厚度和偶极设置。

七张图来自真实结构/数据。小网格的 Bader 电荷、PBE Si 带隙和 Cu 薄层功函数用于学习流程，不是收敛材料基准。声子/弹性完整链和所有剩余命令不属于本包，后续教程另补。

参考：https://qvasp.com/ 与 https://qvasp.com/Tutorials.html （在线手册标记2.25，以当前模块实际2.26帮助为准）。
