Phonopy4.8.3 位移准备、真实力数据整理与QHA操作

在分配CPU内用普通账号新解压目录 bash run.sh，1物理CPU，无GPU、无新DFT。prepare.py实际运行4次phonopy-init和3次phonopy-strain-cells；thermal.py从真实数值HDF5数据运行mesh密度长度24/48两次原生QHA并独立求二次自由能面的极小值，plot.py输出真实图。输出已有则另解压。

inputs包含真实既有Si五体积96/98/100/102/104%原胞体积的数据：10次VASP6.6.1 GPU有限位移力、5次静态能量（原研究已完成），ENCUT550eV、EDIFF1e-6、位移超胞4x4x4电子k网格、2Si静态原胞9x9x9电子k网格。source实际原生phonopy-collect及phonopy-anisotropic-qha-dataset已从完整原始XML执行并逐项核对，完整RAW不公开。本包从数值FORCE_SETS/带力YAML/数值HDF5开始，绝不把数值参考说成新的电子计算。

DFT-inputs保留15条件的INCAR/POSCAR/KPOINTS；势文件需用户自行按授权准备，不包含POTCAR/WAVE/完整OUTCAR/XML。重做DFT须按各自Slurm资源配置提交，run.sh不会提交作业。普通用户已有真实输出时，phonopy-collect phonopy_disp.yaml force1.xml force2.xml写FORCE_SETS；加--sp写包含力的phonopy_params.yaml。两类输入和不同网格点需要严格对应。

QHA数据参考常规胞8Si、原胞2Si和同一个16Si超胞。常规胞的非对角超胞矩阵[0 1 1;1 0 1;1 1 0]与原来2x2x2原胞为同一超胞，原始力按周期位置精确重排，无新力或插值。静态能量与声子自由能都按2Si原胞。原生type1分支使用传统有限差分求FC然后symfc投影对称化；--fc-calculator symfc只对type2密集/随机位移选择拟合器。不能把本例说成无ASR的原始力常数比较。

本例无电子温度自由能、二次F(a)拟合、无晶格平滑，热统计去掉Gamma声学模。mesh标量24/48是密度长度参数，实际声子网格为8^3/15^3，不是24^3/48^3。输入温度0至500K、100K一步，中央差分输出0至400K；最大两网格a差约9.98e-6A。小体系/小网格教学结果，不作为收敛材料预测。更一般各向异性需要更多独立晶格采样；此立方例仅一个自由晶格参数。
