过渡态往反应物方向只跑了一步，但貌似并不属于sob老师博文中的情况

dzdhp

如题，irc只往反应物方向跑了一步，生成物方向正常，详细拜读了sob老师的博文“在Gaussian中计算IRC的方法和常见问题”论坛也看了好多帖子，但未找到解决方法，甚至还觉得我这个irc应该是正常的，但是因为是新手，只能在论坛学习，仍不确定自己的想法，来请教一下大神们的意见。
首先是过渡态，虚频为330.39，且方向连接产物和反应物，看起来的比较合理的，唯一的问题是与反应物挨的太近，这个应该是导致反应物方向只走一步的原因。过渡态搜索用的是混合基组，关键词如下
opt=(calcfc,ts,noeigen) freq m062x/gen nosymm em=gd3


1-5 8-93
6-31g(d,p)
****
6 7
6-31+g(d)
****



6和7号原子是脱氢羧基上的两个氧，为了节省时间只给他们加了弥散。irc用的也是混合基组，与ts一致，的关键词为：
IRC(calcfc,stepsize=3,maxpoint=60) nosymm em=gd3 m062x/gen


因为只往反应物方向跑了一步，参考sob老师的博文使用calcall，得到的过渡态与用calcfc一样，且虚频大小方向一致，所以应该不存在过渡态定位不准确的情况，继续跑过渡态，其他关键词无改动，反应物方向仍为一步。尝试加大步长，一旦超过5就会报错：Maximum number of corrector steps exceded。已经检查过IRC任务第一次输出的SCF Done能量和找过渡态最后一步的SCF Done是一致的。
下面是我跑的结构，图中6,43,53分别是O,H,N。过渡态O-H和N-H距离分别为1.16028，1.35170。反应物O-H和N-H距离分别为1.14710，1.36773。产物O-H和N-H距离分别为1.45117，1.11785。irc路径中过渡态那一点也是能量最高的点。之后单独取两段结构进行优化，都在20步左右就完成，结构与初始结构基本无差异，且无虚频。因为学生是第一次写这样的论文，有点拿不准是否正常，所以请教一下各位老师的意见，敬谢。

sobereva

IRC的stepsize越小越容易出这个问题

博文明确说了，较大IRC步长遇到Maximum number of corrector steps exceded可以靠用LQA、GS2等其它算法，或者recalc/calcall解决

dzdhp

sobereva 发表于 2023-12-28 11:13
IRC的stepsize越小越容易出这个问题

博文明确说了，较大IRC步长遇到Maximum number of corrector steps  ...

谢谢sob老师解答，发帖的时候忘记说了。试过LQA，步长试过10，20,30，均只往反应物方向跑一步。之后换GS2，会报错GS2结构优化失败。calcall也试过了，还是只往反应物方向跑一步，然后就去走另一端了。

ionexchangeC

尝试使用IOp(1/108=-1)

sobereva

dzdhp 发表于 2023-12-28 17:20
谢谢sob老师解答，发帖的时候忘记说了。试过LQA，步长试过10，20,30，均只往反应物方向跑一步。之后换GS2 ...

还可以尝试更大步长

dzdhp

ionexchangeC 发表于 2023-12-28 22:00
尝试使用IOp(1/108=-1)

谢谢建议，这个在论坛看到了，加上之后还是没用，不过现在已经用别的法子解决了

dzdhp

sobereva 发表于 2023-12-29 04:17
还可以尝试更大步长

谢谢sob老师关于步长的建议，之后尝试更大步长，结果只跑三个点，有点怀疑步长太大，遂尝试缩小，现在LAQ加上calcall，步长6，8都试了一下，已经成功跑出来了反应物端，两种步长分别跑出来3,5个点，反应物端曲线都及其平缓，接近直线。

超限制抱怨

dzdhp 发表于 2023-12-29 08:53
谢谢sob老师关于步长的建议，之后尝试更大步长，结果只跑三个点，有点怀疑步长太大，遂尝试缩小，现在LAQ ...

佬，看看你irc的输入文件呗