CREST 3.0.2 水与硫化氢：优化、数值Hessian、振动频率与热化学

解压后进入crest-harmonic，使用sbatch job.slurm提交；或者在已分配的CPU作业内执行bash run.sh。CPU分区按56核整节点分配和计费；资料包同时运行最多4个单线程小计算。不要把本教学脚本用作登录节点生产模拟。

两个孤立小分子的原生ohess任务分别完成结构优化和数值Hessian。随后38个独立原生单点读取完整梯度：每个分子的中心点和9个笛卡尔方向的正负0.0025Bohr位移。原Hessian、独立梯度矩阵、质量加权频率、零点能与热化学公式相互核对。结果在results/，图在plots/，所有新原生输出在cases/。

GFN2/tblite、电子温度300K、accuracy0.0001、extreme优化；热化学使用原生默认298.15K及其温度表。numhess1是Eh/Bohr²的矩阵；非线性三原子分子共有3个内部振动和6个刚体自由度。3.0.2的vibspectrum1内部模强度1000和g98.out的一些99字段为程序占位值，不是真实IR强度、约化质量或Raman强度。图只展示实际频率与热化学修正，不把占位值画成光谱。

G(RRHO) total是加入电子能的热化学修正，已包含零点能，不能重复加ZPVE。孤立分子气相统计不等于溶液、液态水或反应自由能；反应比较还需一致的组分配平、方法、标准态与环境。

reference.json仅是这两个固定教学分子的参考；换体系另建目录准备自己的结构和TOML，不能沿用固定参考比较。
