xTB6.7.1稳定C API：电子能、梯度、电荷、偶极、键级和轨道取数

解压进入xtb-c-api，在已分配CPU作业内bash run.sh；或者sbatch job.slurm。cpu-8173m按56物理核整节点分配和计费，例子只用一个CPU线程，不使用GPU。不在登录节点做生产模拟。

脚本用现有oneAPI2024.2的icx编译C、ifx链接现有libxtb.a与顺序MKL，参数-nofor-main适用于C主程序链接Fortran静态库。没有打包程序二进制，也不需要单独安装Python-xTB绑定。原生开发头和库由现有xtb/6.7.1提供，Python绘图使用sisl/0.16.4。

水、NH4+、OH自由基和H2分别使用GFN1/GFN2；所有模型的坐标、电荷与未成对电子数在inputs/中。XYZ是Å，C API输入DAT是Bohr；本包DAT换算使用0.529177210903Å/Bohr，旧发行版CLI的转换/文本精度会带来极小舍入差，分析已逐坐标比较。实际代码读取每种方法的轨道数，不把GFN1和GFN2的基组维数写死。

当前6.7.1Intel发行包对本例恰好沿坐标轴放置的水构型，梯度与能量有限差分存在明显不一致；API和CLI彼此一致不足以排除这个问题。本包使用一般三维方向，并额外检查完整梯度、总力矩、旋转能量和偶极原点变换；换自己的体系时仍需独立核对。orientation-reference.json只保留可用于能量/电荷/偶极旋转对照的数据，不把原方向的异常梯度当作可靠力参考。

四个分子共8对API/CLI中心点；水的两个方法各9个坐标的正负0.001Bohr独立位移，另有一次对象复制后的坐标更新/还原。每次共52个原生进程、54个电子单点；不会按位移数或反复运行把它们当作54个新物理方法。

raw/是本次真实原生结果，results/含CSV和全数组，plots/为本次新图，reference-plots/为原始示例图。python -B show_results.py直接读取水GFN2结果，不启动计算。

梯度单位Eh/Bohr，原子受力为负梯度。API完整Wiberg矩阵含弱项；CLI的wbo文件只写>0.1的原子对，未列出不等于零。部分电荷具有方法依赖性，不是唯一的形式氧化态；总偶极不能简单用部分电荷乘坐标替代。带电分子的偶极依赖坐标原点，平移t后增加Q*t。孤立分子virial不在此解释成周期应力；MO图不当作晶体能带。没有重叠矩阵时不使用C.T@C检验轨道正交。

此包固定这四个教学输入和参考。另算体系时另建目录，准备自己的电荷、自旋和Bohr坐标，核对方法及API支持范围，不沿用固定参考。api-symbol-coverage.csv列出头文件中的全部入口并标记本例实际调用情况；GFN0/GFNFF、溶剂和外电荷等未在本例运行的接口仍单列，不因有符号就声称已经实算。
