鸢算科技：冻结位移非谐势能与单坐标量子自由能。
calculations共9组Si2单点输入，准备合法Si POTCAR后各目录sbatch run.sh。全部完成后在附件根目录module load ase/3.29.0，python extract_anharmonic_mode.py。
正负u对应两个Si沿x反向移动，其他自由度固定；不做结构优化。能量用与力一致的F；实际力仅检验能量拟合。
一维Schrodinger有限区间模型使用有效质量2mSi；不等于完整四阶力常数/四声子散射/体相非谐自由能/相变温度。RESULTS与figures仅作参考。无POTCAR或原始VASP输出。
