求助Ni金属捕获自由基的过渡态寻找

wjclear

  各位计算大佬，您们好，我最近在计算Ni捕获自由基的过渡态寻找，先通过固定键长扫描获得最高势能面的结构，然后用这个优化结构再进行了一步固定优化获得仅一个虚频（-101.66
）的优化结构，然后取消固定用opt=(calcfc,ts,noeigen)进行计算，可计算结果出错，还请各位大佬指点一二。附件有固定优化的输出文件以及最后寻找过渡态的文件。

sai77

取消限制性优化要读取第一步的波函数

wjclear

sai77 发表于 2026-4-24 07:40
取消限制性优化要读取第一步的波函数

好的 谢谢。我是使用第一步优化的最高势能面的结构作为初猜，进行键长固定优化的。我先按照您说的方式提交运行看看

sai77

wjclear 发表于 2026-4-24 10:59
好的 谢谢。我是使用第一步优化的最高势能面的结构作为初猜，进行键长固定优化的。我先按照您说的方式提 ...

这没问题，固定键长优化后，确认虚频震动方向；取消固定再找过渡态时要读取前面优化好的波函数。你这体系是自旋极化单重态吧，还要注意验证波函数稳定性

wjclear

sai77 发表于 2026-4-24 11:52
这没问题，固定键长优化后，确认虚频震动方向；取消固定再找过渡态时要读取前面优化好的波函数。你这体系 ...

谢谢您，按照您给体的提示，最后确实找到一早期的过渡态中间体，因为虚频只有一个，振动方向也符合，IRC也基本符合，就是成键两原子的距离到了3.2埃，也有在类似计算的文献里找到有这般远距离的早期过渡态的先例。