计算化学公社

 找回密码 Forget password
 注册 Register

IRC问题求助

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

正文摘要:

大家好,我在做C+CO2=CO+O时,得到过渡态,然后做IRC,得到的IRC图线是对称的,并且指向的反应物和生成物是一样的,这是怎么回事呢? 这是我的输入文件: opt=(calcfc,ts,cartesian,maxcycle=200,noeigen) freq b3 ...

回复 Reply

alan 发表于 Post on 2015-12-30 14:24:48
root 发表于 2015-12-30 01:07
是的,我觉得是没有过渡态的。

对于IRC,先不看这个曲线的形状,从两端的构型上看,就是空间取向改变 ...

太谢谢你了~
root 发表于 Post on 2015-12-30 01:07:56
alan 发表于 2015-12-29 23:30
所以可以在文章里说明,这个断键通过relaxed scan后,不考虑过渡态,即不考虑此断键反应路径了,是吧?谢 ...

是的,我觉得是没有过渡态的。

对于IRC,先不看这个曲线的形状,从两端的构型上看,就是空间取向改变了一下,能量略有下降。也就是说这个“过渡态”与解离反应无关。
再说形状,从IRC=0到附近的极大值点,能量变化其实很小,只有2个波数,是因为优化过渡态的时候收敛标准不严格,以及这个构型本来就比较软,在一般的优化步长下难以准确走到驻点。所以这个IRC也不算奇怪。
alan 发表于 Post on 2015-12-29 23:30:30
root 发表于 2015-12-29 23:15
能量几乎是一直升高,到倒数第二个构型的时候,S与C的距离都已经超过4 $\AA$了,十有八九这个解离过程也 ...

所以可以在文章里说明,这个断键通过relaxed scan后,不考虑过渡态,即不考虑此断键反应路径了,是吧?谢谢~~~(这个过渡态算了好久。。)

话说,顺便提一个奇怪的现象,我求的过渡态做IRC(两个方向同时做),得到的轨迹如下图。根本不是身为过渡态的一阶鞍点图形,这个IRC路径是不是也可以从侧面反应出此过渡态的不存在?谢谢!

12121.jpg (182.34 KB, 下载次数 Times of downloads: 37)

12121.jpg

irc.rar

5.26 MB, 下载次数 Times of downloads: 3

root 发表于 Post on 2015-12-29 23:15:37
alan 发表于 2015-12-29 22:47
你好,我最近做了relaxed scan然后求过渡态,虽然求出了一个虚频不大的过渡态,但是做IRC没有得到理想的 ...

能量几乎是一直升高,到倒数第二个构型的时候,S与C的距离都已经超过4 $\AA$了,十有八九这个解离过程也是没有过渡态的。

最后一个点与倒数第二个构型相比在空间上的取向更稳定一些,你用倒数第二个构型做初始猜想找出来的过渡态,得到的也只会是连接了分子在空间相对取向上的两个局域极小值。

4 Angstrom已经可以算是解离产物了,再远算的意义不大。至于最后错误的结束,是因为最后一步的优化没结束,本质原因是两个碎片距离太远了,之间的相对作用力太小
alan 发表于 Post on 2015-12-29 22:47:40
root 发表于 2015-12-29 21:44
从NO+O’出发,随着N与O’距离靠近,势能一路下降就能到达NO2的势阱里,
也就是说,反应路径上不会出现 ...

你好,我最近做了relaxed scan然后求过渡态,虽然求出了一个虚频不大的过渡态,但是做IRC没有得到理想的结果,我怀疑我scan找的鞍点不对(by the way,我的柔性扫描没有正常结束,但看到一个弱弱的鞍点,就拿来求过渡态了)。可以麻烦你看看我的scan为什么错误结束吗?以及,图中倒数第二点是否可以认为是鞍点?

1111.png (8.09 KB, 下载次数 Times of downloads: 32)

1111.png

scan.gjf

3.02 KB, 下载次数 Times of downloads: 1

scan.log

8.99 MB, 下载次数 Times of downloads: 2

root 发表于 Post on 2015-12-29 21:53:45
alan 发表于 2015-12-29 15:32
请问你这个柔性扫描得出的结果,是怎么判断没有过渡态的呢?

有一个很有意思的例子是O(1D)+CH4,它有四个无过渡态的反应通道和4个有过渡态通道,感兴趣的话可以去看看,http://dx.doi.org/10.1039/C5CP04278J





root 发表于 Post on 2015-12-29 21:44:33
alan 发表于 2015-12-29 15:32
请问你这个柔性扫描得出的结果,是怎么判断没有过渡态的呢?

从NO+O’出发,随着N与O’距离靠近,势能一路下降就能到达NO2的势阱里,
也就是说,反应路径上不会出现垒,当然也就不存在过渡态了

再比如H2O、NH3分子解离掉一个H原子,也没有过渡态,类似的例子很多很多
alan 发表于 Post on 2015-12-29 15:32:20
dy19930508 发表于 2015-11-21 19:27
sobereva老师,我做了已经优化过的NO2其中一个N-O键的距离柔性扫描,得到一个抛物线型的曲线,这样是不是 ...

请问你这个柔性扫描得出的结果,是怎么判断没有过渡态的呢?
dy19930508 发表于 Post on 2015-11-21 22:51:33
sobereva 发表于 2015-11-21 21:25
这就是说明没过渡态
键角没必要扫了,因为你柔性扫描键长的时候键角已经自动优化到最佳了。

十分感谢老师的指导!
sobereva 发表于 Post on 2015-11-21 21:25:41
dy19930508 发表于 2015-11-21 19:27
sobereva老师,我做了已经优化过的NO2其中一个N-O键的距离柔性扫描,得到一个抛物线型的曲线,这样是不是 ...

这就是说明没过渡态
键角没必要扫了,因为你柔性扫描键长的时候键角已经自动优化到最佳了。
dy19930508 发表于 Post on 2015-11-21 08:38:39
liyuanhe211 发表于 2015-11-21 04:11
C 2s2, 2p2原子是个双自由基,这么搞多重度好像不大和谐

嗯,谢谢您的指点,十分感谢!
dy19930508 发表于 Post on 2015-11-21 08:37:15
sobereva 发表于 2015-11-21 02:40
你找到的都是过渡态,这种IRC也是合理的,但是并不对应于你说的化学反应的过渡态。
这两个反应是否有过 ...

恩,我做一下柔性扫描,十分感谢sobereva老师!
liyuanhe211 发表于 Post on 2015-11-21 04:11:12
dy19930508 发表于 2015-11-20 17:57
是C原子,是不是我的过渡态找错了。

C 2s2, 2p2原子是个双自由基,这么搞多重度好像不大和谐

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

GMT+8, 2024-11-27 20:28 , Processed in 2.879835 second(s), 26 queries , Gzip On.

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