鸢算科技：热位移带隙与一次有效声子拟合。
25个si-thermal目录按POSCAR顺序准备合法Si POTCAR，再逐目录sbatch run.sh。NSW=0保留快照位移，不进行优化。
完成后在附件根目录module load ase/3.29.0，python extract_thermal_gap.py；需要有效力常数时module load phonopy/4.8.3，python extract_effective_fc.py。后者使用已有symfc，比较6快照模型的保留样本力误差和8快照最终模型。
原始二阶力常数si-phonon-v100-result.yaml来自两个旧有限位移计算目录（本包也含其输入），可按官网晶体声子教程重算。THERMAL-GAP-PLAN.json含温度、配对编号与位移RMS。
0K包含量子零点位移；每温度4个独立样本及其负位移，不把8个快照当成8个独立样本。固定16原子超胞体积，未做热膨胀/完整AHC/自洽声子/高温相稳定性。RESULTS.json和图仅为参考，不随新任务自动改变。
附件不含POTCAR、波函数或私有原始输出。
