xTB6.7.1 C接口补充：GFN0/GFN-FF、外部点电荷梯度、GBSA有限网格

在已分配的CPU作业内bash run.sh，或者sbatch job.slurm。CPU分区按56物理核整节点分配与计费；例子使用一个线程，不使用GPU。不在登录节点运行生产计算。
使用已有xtb/6.7.1头文件和静态库、oneapi/2024.2的icx编译C和ifx链接；Python绘图来自sisl/0.16.4。不安装新共享软件，不打包二进制库。
流程：编译助手→GFN0与GFN-FF气相→GFN1/GFN2外部电荷→GFN2水GBSA5810网格→全坐标能量差分→230/770网格对照→输出CSV、全数组和两张图。
共121原生计算进程、126能量计算，其中107个电子单点和19个GFN-FF力场单点；不是126种物理方法。方法选择器3代表GFN-FF，不是GFN3。
inputs/water.dat坐标Bohr，首行原子数/电荷/未成对电子数；inputs/pointcharges.dat首行电荷数，随后原子序数/电荷e/x/y/z(Bohr)。点电荷固定，不是自洽极化MM，也没有MM内部相互作用能。
能量Eh，梯度Eh/Bohr，受力是负梯度。外电荷的6个分量和分子的9个分量分别做中心差分；总平移/转动关系必须包含点电荷。卸载外电荷/溶剂后会重算并回到对应气相结果。
此6.7.1 C接口xtb_setSolvent调用GBSA/Still，不等同于新版本文档或CLI的ALPB。三网格能量和力矩由实际输出获得。5810角点仍有约1.11e-6 Eh转动力矩残差，高于原1e-7目标；不称高精度旋转不变性已收敛。梯度差分通过不等于网格已收敛。
results/VALIDATION.json只表示列出的接口操作、差分和数据复跑检查通过；results/GBSA-ROTATIONAL-LIMIT.json单独保留未通过的原始转动力矩门槛。不得用前一个passed=true宣称所有物理精度通过。
raw/是本次实际原生计算；results/为本次CSV/全数组，plots/是本次新图，reference-plots/为原示例图。第一次完整运行生成raw后，再次查看只运行python -B verify.py --reference和python -B plot.py，不会重复单点。修改体系请另建目录并重设参考，不沿用本例的固定数值。
