qvasp：Si16真实DFPT声子与GaAs弹性教程

只看已计算的数据和图片：
  bash view.sh
得到 analysis/*.csv、figures/*.png。参考力常数、原始频率YAML和qvasp数据在reference中；声学支微小负值原样保留，没有ASR或裁零。

重新计算：
  bash prepare.sh
  sinfo -p cpu-8173m
  sbatch job.slurm

先确认CPU队列可用。模板申请整节点56核，实际同时运行两个条件、各8MPI和1线程，合计16物理核。不要把计算脚本直接当作普通登录节点生产任务运行。run-allocated.sh在已分配资源内运行VASP并调用post.sh。后处理可在自己的轻量工作目录运行。当前示例是CPU路径，不是GPU加速数据或CPU节点性能比较。

实际VASP6.5.1CPU、Intel oneAPI2024.2/IntelMPI2021.13，Si16来自已优化Si2原胞的2x2x2超胞；ENCUT450eV、EDIFF1E-6、IBRION8、Gamma积分2^3/3^3。NCORE1、KPAR1，每个条件8MPI。运行脚本仅为当前任务设置unlimited stack，不改系统配置。qvasp2.26调用Phonopy2.43.6，以521.471输出cm^-1。view.sh用独立Phonopy4.7.2 Rust和相同转换因子核对，不把cm^-1写成THz。

Si原始Gamma光学频率约513.353763/505.381113 cm^-1，两网格相差1.58%；声学最小值约-0.0033356 cm^-1，不删去。这是有限超胞、有限积分网格的教学例，不是完整收敛研究；qvasp完整声子接口需要实际VASP DFPT Hessian，普通有限位移单点的XML不能直接代替。

GaAs使用已有真实VASP6.6.1、IBRION6/NFREE2/ISIF3计算的总弹性张量，不计作新DFPT。C11=96.40676、C12=40.54473、C44=48.93888GPa；体模量59.16540667GPa。原生Voigt顺序XX YY ZZ XY YZ ZX；1kBar=0.1GPa。总弹性包含离子放松贡献，不能与固定离子的张量混用。这里没有把IBRION8写成可直接计算完整弹性张量。

不包含POTCAR、WAVECAR或完整原始OUTCAR/XML。POTCAR由prepare.py按本地许可势生成，重新计算会在自己的work目录自然产生原始文件。下载包只分发自己的实际数值力常数、频率、结构和输入；GaAs原始张量行明确标注为数值摘录，view.sh在analysis/GaAs-native中生成用于qvasp-elasticd的精简练习OUTCAR，并不伪装成完整计算输出。

环境：qvasp/2.26、phonopy/4.7.2；qvasp的phonond子进程使用phonopy/2.43.6。Phonopy两版为GNU13.3 C/OpenMP后端，4.7.2另含phonors0.5.0 Rust；数学库线程固定1。
