计算化学公社

标题: 求助：轴手性分子（C/H/O/N/S/I）IRC只走两步自动停止，是否为TS不对 [打印本页]

作者
Author: TIANTTT 时间: 2025-4-4 20:32
标题: 求助：轴手性分子（C/H/O/N/S/I）IRC只走两步自动停止，是否为TS不对
本帖最后由 TIANTTT 于 2025-4-4 21:23 编辑

各位老师好！初学者冒昧前来求助。我目前在算一个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
如果有仗义侠士出手相救,我不甚感激！

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

作者
Author: TIANTTT 时间: 2025-4-4 21:17

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

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

作者
Author: KAIMISITERUI 时间: 2025-4-5 00:25
你分子的结构式里确定是有个C-S双键的吗？

作者
Author: TIANTTT 时间: 2025-4-5 14:01

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

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

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

作者
Author: TIANTTT 时间: 2025-4-5 16:57

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

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

作者
Author: sai77 时间: 2025-4-6 16:02
你找的这个过渡态可能高估了旋转能垒，因为从扫描曲线看另一个高点对应的结构能量更低一点

欢迎光临计算化学公社 (http://bbs.keinsci.com/)