这套思路的起始点，估计是Daniel A. Singleton在2008年发的一篇文章，你找一下DOI号为10.1021/ja802577v里“dynamic matching”的内容。真要抠细节的话，去看proggenHP文件里的代码。

不过你首先应该想的不是重复文献中的做法，而是你到底要解释什么实验现象。是post-transition state bifurcation，还是某个intermediate尚未充分relax的时候反应就进入到了下一步，还是某个intermediate在很平缓的势能面上瞎晃悠。