求助NO2+进攻乌洛托品过程，柔性扫描应该怎么定位特定的C-N断裂？

freshman234

      在复现一篇文献的乌洛托品（六亚甲基四胺）硝解反应的DFT计算，在复现过程中，我使用了与文章中相同的初始结构（乌洛托品,由剑桥晶体学数据中心获得）、理论方法和基组（m062x/6-311G(d)），引入了SMD模型来模拟发烟硝酸环境。然而，经过很多尝试，我找到的过渡态虚频（-363 cm-1，文章中是-317cm-1）以及后续IRC两端优化的反应物与产物的E+ZPE差值（shermo程序计算）与文章结果存在较大偏差，由于此前是根据文献图片手搭的，所以怀疑是过渡态初猜结构没有搭好，现在想尝试柔性扫描帮忙找过渡态。
      柔性扫描思路是让硝化剂NO2+去靠近乌洛托品中的N原子，把NO2+与N原子距离从1.5扫到1.3埃，但结果是能量一直在增长（文献里过渡态时的距离是1.4埃），并且乌洛托品的三个C-N键同时拉长且没有一个是显著于其它两个拉长的，这与文献仅断裂一个C-N键不符，求助大佬们这种硝解反应应该怎么设置柔性扫描？

wzkchem5

仔细阅读http://sobereva.com/678