求助：轴手性分子（C/H/O/N/S/I）IRC只走两步自动停止，是否为TS不对

TIANTTT

各位老师好！初学者冒昧前来求助。我目前在算一个55个原子的轴手性分子（C/H/O/N/S/I），已经以手性轴所在二面角进行了柔性扫描，然后我以扫描结构的相对能量最高点对应化合物 作为初猜结构进行TS和IRC验证，但总是出现对应虚频并非我要的运动方向的情况.更严重的是，我将看上去比较合理的过渡态进行IRC验证也总出现只跑两步自动结束但不报错的情况。
我的输入文件关键词如下#p  irc=(forward,reverse,Maxpoints=50,stepsize=30,calcall,LQA) b3lyp/genecp em=gd3bj。
柔扫log文件我附上链接 https://pan.baidu.com/s/1bKEFh6u9SZSE1OD3kan4IA?pwd=4mgf 提取码: 4mgf
如果有仗义侠士出手相救,我不甚感激！

Stardust0831

给出TS任务的log文件和反应方程式。
网盘里文件上传前最好也压缩一下。
请不要同时在多处重复提问。

TIANTTT

Stardust0831 发表于 2025-4-4 20:48
给出TS任务的log文件和反应方程式。
网盘里文件上传前最好也压缩一下。
请不要同时在多处重复提问。

非常感谢您的提醒！这是补充、压缩后的文件： https://pan.baidu.com/s/1bKEFh6u9SZSE1OD3kan4IA?pwd=4mgf  TS反应物生成物过渡态均为同一物质，仅由于手性轴存在，分子会发生消旋，为探求这一过程故求旋转势垒。

KAIMISITERUI

你分子的结构式里确定是有个C-S双键的吗？

TIANTTT

KAIMISITERUI 发表于 2025-4-5 00:25
你分子的结构式里确定是有个C-S双键的吗？

是的，但opt之后自动变成C-S=O这样的结构，这个会产生较大影响吗？

Stardust0831

帮你试了一下。stepsize=40的时候就能跑出来。默认步长不能跑是因为TS附近的势能面太平缓了。相关知识可以看：《在Gaussian中计算IRC的方法和常见问题》
在ORCA里用NEB可以在找到TS的同时直接获得MEP（MEP和IRC类似，不过是定义在非质权坐标下的），此时可以避免反应路径由于势能面过于平缓而跑不出来的问题。

TIANTTT

Stardust0831 发表于 2025-4-5 15:41
帮你试了一下。stepsize=40的时候就能跑出来。默认步长不能跑是因为TS附近的势能面太平缓了。相关知识可以 ...

非常感谢您的帮助！我会认真阅读的！

sai77

你找的这个过渡态可能高估了旋转能垒，因为从扫描曲线看另一个高点对应的结构能量更低一点