计算化学公社

 找回密码 Forget password
 注册 Register

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

查看数: 17699 | 评论数: 15 | 收藏 Add to favorites 3
关灯 | 提示:支持键盘翻页<-左 右->
    组图打开中,请稍候......
发布时间: 2014-11-7 20:11

正文摘要:

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

回复 Reply

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

请问就你说的这种情况你最后怎么解决的安?
wyy015020 发表于 Post on 2019-1-22 17:18:47
感谢sob老师,按您说的方法,我的问题解决了,谢谢您
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-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-20 15:19:27
wyy015020 发表于 2019-1-20 11:50
sob老师,不好意思再打扰您一下,我刚刚看了优化过渡态的能量是E=-1227.36780221,然后将这个过渡态结构跑IR ...

8成是做优化和跑IRC用的影响能量的关键词不同
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,就这个出现了这种问题
wyy015020 发表于 Post on 2019-1-20 11:27:04
sobereva 发表于 2019-1-18 00:18
看虚频模式,足够80%确定过渡态是不是你想要的
即便过渡态不是你想要的,只要过渡态是准确的势能面的一 ...

好的,谢谢您,我再试试
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-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对过渡态重新优化,发现过渡态的结构和之前的结构没有区别,键长键角以及能量都差不多,请问我还有别的尝试方法吗?还是这个过渡态结构压根就是不对的?
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老师
sobereva 发表于 Post on 2019-1-15 02:36:24
wyy015020 发表于 2019-1-14 19:43
您好,我现在遇到了和您一样的问题,就是用TS找过渡态,正常收敛,有一个虚频,虚频模式正确,但是跑irc ...

仔细看
在Gaussian中计算IRC的方法和常见问题
http://sobereva.com/400http://bbs.keinsci.com/thread-7736-1-1.html
wyy015020 发表于 Post on 2019-1-14 19:43:40
chenxin199261 发表于 2014-11-11 12:22
按您说的,我的问题解决了,谢谢你

您好,我现在遇到了和您一样的问题,就是用TS找过渡态,正常收敛,有一个虚频,虚频模式正确,但是跑irc的时候只能往产物方向跑,往反应物跑几步就会折返或者终止。所以想请问您最后是怎么解决的呢?重新找到的过渡态和之前的结构差别大吗?
chenxin199261 发表于 Post on 2014-11-11 12:22:36
sobereva 发表于 2014-11-7 20:54
你的多重度设的是2,电荷是+5,和你描述不符。

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

按您说的,我的问题解决了,谢谢你
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重新计算。
谢谢您

手机版 Mobile version|北京科音自然科学研究中心 Beijing Kein Research Center for Natural Sciences|京公网安备 11010502035419号|计算化学公社 — 北京科音旗下高水平计算化学交流论坛 ( 京ICP备14038949号-1 )|网站地图

GMT+8, 2024-11-24 04:46 , Processed in 0.176207 second(s), 26 queries , Gzip On.

快速回复 返回顶部 返回列表 Return to list