qvasp提取优化能量与批量对照
用Si的五组真实应变结果，提取零展宽能量、检查优化结束并整理表格。
阅读本篇目录
本篇使用鸢算科技集群实际生成的结构或已经完成的VASP计算数据。环境为官方Linux预编译qvasp2.26（模块 qvasp/2.26），辅助环境为ASE3.29.0、Spglib2.7.0、NumPy和Matplotlib；qvasp原生二进制的编译器版本未公开确认，不将它写成本地重新编译版本。
下载输入、真实数值数据与绘图脚本。解压进入目录，以普通账号运行：
bash run.sh
运行后查看 work/cases、work/analysis-core、work/analysis-extra 和 work/figures。这一命令执行文件准备和后处理，不启动VASP。包内不含POTCAR、WAVECAR或完整原始OUTCAR/XML；精简OUTCAR由明确标注的实际数值摘录生成，仅供工具练习。处理自己的任务时应使用自己的完整原始文件。
1 在一个计算目录读取#
cd work/cases/energy-relaxed-+0.000
qvasp -e
本例输出−43.39961206 eV及收敛提示。这里读取的是 energy(sigma->0)，不是电子自由能 TOTEN。比较不同结构时，能量定义、原子数、势、泛函、ENCUT、k点和展宽处理应一致。
2 在父目录一次读取多个条件#
cd work/cases/energy-batch-relaxed
qvasp -e
程序扫描本例五个 strain0至 strain4子目录中的OUTCAR，给出能量和收敛列。固定原子组在 energy-batch-clamped；放松原子组在 energy-batch-relaxed。
| xx单轴应变 | 放松原子后的E(sigma→0) / eV | 晶胞体积 / ų |
|---|---|---|
| -1.0% | -43.39196341 | 161.935220 |
| -0.5% | -43.39776058 | 162.753074 |
| +0.0% | -43.39961206 | 163.570929 |
| +0.5% | -43.39755051 | 164.388784 |
| +1.0% | -43.39161905 | 165.206638 |
图中两组分别以各自零应变能量为零，避免用零点差掩盖曲线形状。五组放松计算的最终最大力为0.0145–0.0151 eV/Å，低于0.02阈值。模型为Si8，xx方向施加固定应变，放松组只放松内部原子；这不是各向同性体积扫描，也没有由这五点宣称完整弹性张量。
3 回到原始输出检查结束状态#
处理自己的完整OUTCAR时可用：
grep 'energy(sigma->0)' OUTCAR | tail -n 1
grep 'reached required accuracy' OUTCAR
grep 'because EDIFF is reached' OUTCAR | tail -n 1
tail -n 20 OUTCAR
qvasp -e的YES按相应结束标记判断，不能代替核查所有物理设置。静态自洽和离子优化的停止条件也不同。没有完整结束、电子未收敛或力仍过大的任务，应先检查输入及输出，再做能量比较。含外压时工具可能切换为焓和PV输出，本篇零外压数据表不将两者混合。
4 导出CSV和图片#
资料包的 work/analysis-core/Si-strain-energies.csv包含实际能量和相对值，前五行为固定原子组，后五行为放松组。所有图由 plot.py读取数值生成，可用普通文本工具查看：
head work/analysis-core/Si-strain-energies.csv
参考：qvasp官方教程与模块实际 qvasp -help。在线手册标记2.25，个别行为与当前2.26不同，本文按本例实际输入输出说明。
继续：完整声子与弹性操作#
已经准备好基本输入后，可接着看qvasp完整DFPT声子流程与GaAs弹性张量提取。前者包含Si16实际DFPT、两种k网格、匹配DIM、Phonopy提取与画图；后者展示实际总弹性张量、单位换算和离子放松含义。都提供终端操作、完整输入和真实数值资料，不将模板生成等同于计算完成。