xTB 近似反应路径：起点、终点与输出检查

下载完整输入、参考轨迹、CSV和画图程序，保存为example.zip：

python3 -m zipfile -e example.zip .
cd xtb-dynamics-path-examples
bash results.sh
bash plot.sh

环境为xtb/6.7.1，官方Intel版，固定提交edcfbbe。例子使用单线程CPU，配套job.slurm负责加载模块和运行；下文核心命令用于说明每一步实际执行的内容。原始结果在reference-results，重算结果在work。短演示用于学会操作，生产计算需要另外选择时长、精度和采样方案。

1. 先准备两端结构#

cat inputs/ammonia-inversion-path/input.xyz
cat inputs/ammonia-inversion-path/product.xyz

两端是同样4个原子、相同顺序的氨分子镜像结构。反应路径要求原子一一对应，不能把产物XYZ重新排序后直接使用；总电荷和自旋也需要保持一致。

2. 控制文件#

$path
 nrun=1
 npoint=12
 anopt=5
 kpush=.003
 kpull=-.015
 ppull=.05
 alp=1.2
$end

npoint设置路径图像数，anopt限制本次近似路径优化轮次。推拉参数影响探索方式，不能视为通用材料常数。这里故意保留短例，展示12图像文件和后处理。

sbatch job.slurm ammonia-inversion-path

实际命令：

xtb input.xyz --gfn 2 --path product.xyz --input settings.inp --parallel 1

3. 检查路径与结束原因#

grep -E 'path trials|barrier|cycles reached' reference-results/ammonia-inversion-path/stdout.log
head -n 12 reference-results/ammonia-inversion-path/xtbpath.xyz

本例达到设定的5轮上限，得到了12帧有限坐标和能量，属于初始近似路径。程序正常退出并不表示已收敛。输出中的末端路径方向重叠还可能显示NaN，不要用这一项计算结论；应检查实际图像、能量以及中间结构是否连续。

xtbpath.xyz注释中的energy在这个入口是相对第一图像的kcal/mol，与MD的轨迹能量Eh、模式扫描的相对能量Eh都不同。xtbpath_ts.xyz中的能量又采用程序该文件的原生定义，不能仅凭同样的energy标签混在一张图上。

4. 可视化与下一步#

head -n 6 ammonia-approximate-path.csv
bash plot.sh

图中的能量曲线不作为氨反转势垒推荐值。本例还包含端点结构松弛和路径搜索影响，因此与沿单个模式的扫描曲线不能直接当作相同方法比较。

真正研究反应时，应先优化两端，检查路径图像，增加必要的搜索与优化，再把合理的候选高点交给适当的过渡态或NEB方法。收敛后检查唯一相关虚频及连接方向，最后才讨论能垒和热修正。参考官方路径说明。