TS计算正常收敛，虚频正常；但是IRC曲线错误的原因

chenxin199261 · 发表于 Post on 2014-11-7 20:11:24

我用TS对一个过渡态进行优化，结果收敛，检查虚频，虚频方向正确。
我用这个过渡态结果计算IRC。但结果显然是错误的。通过看动画发现，原子向反应物方向走几步，又反返回了。
IRC输入文件输入如下

################################################
%nprocshared=24
%chk=FeONH36-TS-5-IRC-p.chk
# irc=(maxpoints=10,calcall) ub3lyp/genecp geom=connectivity scf=xqc

5 2
....

Fe 0
lanl2dz
****
C H O N 0
6-31g(d,p)
****

Fe 0
lanl2dz
###############################################

由于多重度是5，计算SCF收敛有点困难，故加上XQC

麻烦各位了，如果遇到同样问题，麻烦指导一下。

sobereva · 发表于 Post on 2014-11-7 20:54:19

你的多重度设的是2，电荷是+5，和你描述不符。

先确认过渡态结构定位得很准确，用int=ultrafine结合tight优化收敛限，并且最好calcall，然后再做IRC试试。

如果问题依旧，把IRC的步长调小点，或者尝试别的IRC生成算法。

chenxin199261 · 发表于 Post on 2014-11-7 20:59:13

sobereva 发表于 2014-11-7 20:54
你的多重度设的是2，电荷是+5，和你描述不符。

先确认过渡态结构定位得很准确，用int=ultrafine结合tigh ...

好的，我在论坛上打错了，gaussian计算时多重度跟电荷没反。

我TS优化calcall已经用上了。
接下来。我打算按您说的int=ultrafine 并结合tight重新计算。
谢谢您

chenxin199261 · 发表于 Post on 2014-11-11 12:22:36

sobereva 发表于 2014-11-7 20:54
你的多重度设的是2，电荷是+5，和你描述不符。

先确认过渡态结构定位得很准确，用int=ultrafine结合tigh ...

按您说的，我的问题解决了，谢谢你

wyy015020 · 发表于 Post on 2019-1-14 19:43:40

chenxin199261 发表于 2014-11-11 12:22
按您说的，我的问题解决了，谢谢你

您好，我现在遇到了和您一样的问题，就是用TS找过渡态，正常收敛，有一个虚频，虚频模式正确，但是跑irc的时候只能往产物方向跑，往反应物跑几步就会折返或者终止。所以想请问您最后是怎么解决的呢？重新找到的过渡态和之前的结构差别大吗？

sobereva · 发表于 Post on 2019-1-15 02:36:24

wyy015020 发表于 2019-1-14 19:43
您好，我现在遇到了和您一样的问题，就是用TS找过渡态，正常收敛，有一个虚频，虚频模式正确，但是跑irc ...

仔细看
在Gaussian中计算IRC的方法和常见问题
http://sobereva.com/400（http://bbs.keinsci.com/thread-7736-1-1.html）

wyy015020 · 发表于 Post on 2019-1-17 20:13:16

sobereva 发表于 2019-1-15 02:36
仔细看
在Gaussian中计算IRC的方法和常见问题
http://sobereva.com/400（http://bbs.keinsci.com/threa ...

谢谢sob老师

wyy015020 · 发表于 Post on 2019-1-17 20:37:29

sobereva 发表于 2019-1-15 02:36
仔细看
在Gaussian中计算IRC的方法和常见问题
http://sobereva.com/400（http://bbs.keinsci.com/threa ...

sob老师，我现在加了tight和int=ultrafine优化过渡态，然后进行irc分析两边曲线走出来还是一样的，调大步长也不行。现在用calcall对过渡态重新优化，发现过渡态的结构和之前的结构没有区别，键长键角以及能量都差不多，请问我还有别的尝试方法吗？还是这个过渡态结构压根就是不对的？

sobereva · 发表于 Post on 2019-1-18 00:18:15

wyy015020 发表于 2019-1-17 20:37
sob老师，我现在加了tight和int=ultrafine优化过渡态，然后进行irc分析两边曲线走出来还是一样的，调大步 ...

看虚频模式，足够80%确定过渡态是不是你想要的
即便过渡态不是你想要的，只要过渡态是准确的势能面的一阶鞍点，原理上IRC总能正常跑下去
走IRC相关问题该说的我在博文里都说了，没有可以补充的

wyy015020 · 发表于 Post on 2019-1-20 11:27:04

sobereva 发表于 2019-1-18 00:18
看虚频模式，足够80%确定过渡态是不是你想要的
即便过渡态不是你想要的，只要过渡态是准确的势能面的一 ...

好的，谢谢您，我再试试

wyy015020 · 发表于 Post on 2019-1-20 11:50:33

本帖最后由 wyy015020 于 2019-1-21 20:35 编辑

sob老师，不好意思再打扰您一下，我刚刚看了优化过渡态的能量是E=-1221.36780221,然后将这个过渡态结构跑IRC，对应坐标为0的能量却是-1221.37933361，请问这是什么原因呢？我是将优化出的TS结构直接写在IRC输入文件里的，其他的过渡态我都是这样操作的，都能正常的跑IRC，就这个出现了这种问题

sobereva · 发表于 Post on 2019-1-20 15:19:27

wyy015020 发表于 2019-1-20 11:50
sob老师，不好意思再打扰您一下，我刚刚看了优化过渡态的能量是E=-1227.36780221,然后将这个过渡态结构跑IR ...

8成是做优化和跑IRC用的影响能量的关键词不同

wyy015020 · 发表于 Post on 2019-1-21 09:24:57

老师，我把关键词发给您看一下，优化过渡态：# M062X/6-311G(d,p) opt(calcfc,ts,noeigen,maxstep=3,notrust) freq scrf=smd int=ultrafine
跑IRC: # M062X/6-311G(d,p) irc(calcfc,maxpoint=20,stepsize=15,LQA) int=ultrafine scrf=smd
还请您指导一下

sobereva · 发表于 Post on 2019-1-21 16:33:57

wyy015020 发表于 2019-1-21 09:24
老师，我把关键词发给您看一下，优化过渡态：# M062X/6-311G(d,p) opt(calcfc,ts,noeigen,maxstep=3,notrus ...

光从这上看不出什么来。走IRC时候加上guess=read读取优化任务的chk再试

wyy015020 · 发表于 Post on 2019-1-22 17:18:47

感谢sob老师，按您说的方法，我的问题解决了，谢谢您

		自动登录 Automatic login	找回密码 Forget password
密码 Password			注册 Register

[新手求助] TS计算正常收敛，虚频正常；但是IRC曲线错误的原因