pymatgen2026.9.24：实际Si五点DFT E(V)数据+各向同性准谐Debye近似。原表每原子数据按8原子晶胞换算，Poisson比假定0.25，零压100至500K。不是逐体积声子QHA或已收敛Si热物性。
bash results.sh
bash plot.sh
自行分析：bash run.sh，使用新解压目录生成analysis。
随包原图在reference-figures，重新绘图写到figures。
