24小时热门版块排行榜    

查看: 2074  |  回复: 30

ellemiu

铁杆木虫 (正式写手)

[求助] IRC 结果分析

TS寻找过渡态,完了以后用IRC扫描势能面,顺利收敛,,在分析结果的时候却意外的发现,irc的势能面扫描过程中,开始点和结束点的几何构型几乎一样,反应产物,感觉只是键在来回震荡。

好吧,问题出来了
1、这样的结果说明什么? 如果这是一个正常的结果,是不是意味着我的反应不能发生

2、如果不是正常的结果,,改如何避免这样的情况呢?
拜请各位大神不吝赐教。

本人计算的是一个是H原子的迁移过程,涉及自由基。
IRC输入文件如下:
%chk=h3tsirc2.chk
%nprocshared=4
%mem=4GB
# b3lyp/6-31g scfcyc=300 irc=(calcfc,maxpoints=20)

TS-Ph-H-OC-Ph

0 2

第一张,irc扫描完以后的结果

irc扫描完以后的结果
回复此楼
业精于勤荒于嬉,成于思而毁于随
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
回帖支持 ( 显示支持度最高的前 50 名 )

necturnal

铜虫 (初入文坛)

引用回帖:
10楼: Originally posted by ellemiu at 2012-05-22 08:46:16:
昨天去掉了scfcyc=300,结果出现了9999的错误。看样子还是这个不能少。

我找找IOp的帖子,谢谢

那你去掉scfcyc=300的同时加了iop(5/13=1)吗?
11楼2012-05-22 11:34:20
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
普通回帖

meteoric30

木虫 (正式写手)

【答案】应助回帖

★ ★
感谢参与,应助指数 +1
ellemiu: 金币+2, 有帮助, irc直接拿的TS的结果来继续的,不太明白irc的产物是指什么,你指的是最后一步的状态进行优化吗》?步长设大一点怎样设呢? 2012-05-18 16:19:23
步长设大一点,irc的产物优化了吗?优化的结果看看是不是和ts一样
2楼2012-05-18 15:35:18
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

ellemiu

铁杆木虫 (正式写手)

引用回帖:
2楼: Originally posted by meteoric30 at 2012-05-18 15:35:18:
步长设大一点,irc的产物优化了吗?优化的结果看看是不是和ts一样

irc直接拿的TS的结果来继续的,不太明白irc的产物是指什么,你指的是最后一步的状态进行优化吗?

还有步长设大一点怎样设呢?
业精于勤荒于嬉,成于思而毁于随
3楼2012-05-18 16:20:22
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

ldzh0531

木虫 (小有名气)

【答案】应助回帖


感谢参与,应助指数 +1
ellemiu: 金币+1, 有帮助, 我算出的虚频为-1415,其实我只是想扫势能面看看TS是不是正确的 2012-05-21 10:14:58
IRC只是根据频率计算单点,并不优化。想要得到准确构型要提取出终端点优化。另外,通常氢转移的虚频要大于1000
4楼2012-05-18 17:58:22
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

necturnal

铜虫 (初入文坛)

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
感谢参与,应助指数 +1
gmy1990: 金币+2 2012-05-19 08:59:34
ellemiu: 金币+4, ★★★很有帮助, 分子体系是比较小,我先按照你给的建议来试试。对于IOP的设定,我一直弄的不是很懂,能否做一下解释? 2012-05-21 10:21:57
ellemiu: 金币+4, ★★★很有帮助, 楼主很热心,追加金币 2012-05-22 14:46:41
ellemiu: 金币+5, ★★★很有帮助, 很热心的网友。 2012-05-26 10:21:29
按照你的说法,这个IRC应该不是一个正确的IRC,不过这不代表你的TS就不对。我觉得可以尝试 去掉scfcyc=300,因为从能量上看,你的分子不是很大,小分子迭代次数过多,容易结构不合理,也就是可能跑偏了,改用iop(5/13=1)。还有就是可以尝试减小步长(stepsize=10是default)。
5楼2012-05-19 06:55:35
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

ellemiu

铁杆木虫 (正式写手)

引用回帖:
4楼: Originally posted by ldzh0531 at 2012-05-18 17:58:22:
IRC只是根据频率计算单点,并不优化。想要得到准确构型要提取出终端点优化。另外,通常氢转移的虚频要大于1000

我算出的虚频为-1415,其实我只是想扫势能面看看TS是不是正确的
业精于勤荒于嬉,成于思而毁于随
6楼2012-05-21 10:15:58
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
7楼2012-05-21 10:22:17
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

ellemiu

铁杆木虫 (正式写手)

引用回帖:
5楼: Originally posted by necturnal at 2012-05-19 06:55:35:
按照你的说法,这个IRC应该不是一个正确的IRC,不过这不代表你的TS就不对。我觉得可以尝试 去掉scfcyc=300,因为从能量上看,你的分子不是很大,小分子迭代次数过多,容易结构不合理,也就是可能跑偏了,改用iop( ...

分子体系是比较小,我先按照你给的建议来试试。对于IOP的设定,我一直弄的不是很懂,能否做一下解释?

其实还有一个问题,,什么样的情况下才知道我的TS是不是合理的呢
业精于勤荒于嬉,成于思而毁于随
8楼2012-05-21 10:22:21
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

necturnal

铜虫 (初入文坛)

引用回帖:
8楼: Originally posted by ellemiu at 2012-05-21 10:22:21:
分子体系是比较小,我先按照你给的建议来试试。对于IOP的设定,我一直弄的不是很懂,能否做一下解释?

其实还有一个问题,,什么样的情况下才知道我的TS是不是合理的呢

iop(5/13=1)的意思忽略不收敛的failure,继续算下去。
高斯官网有个IOp reference,
http://www.gaussian.com/g_tech/g_iops/iops_toc.htm
大部分IOp的意思可以在里面查到,但是还是不全。Sobereva大神有几个帖子专门介绍IOp的,感兴趣的可以搜一下看看。
一般情况下,判断TS是否合理,可以看虚频,如果虚频的指向是你所需要的成键和断键的方向,一般情况可以初步判断合理,但是还是需要通过IRC才算是比较准确了。
9楼2012-05-22 08:35:41
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

ellemiu

铁杆木虫 (正式写手)

引用回帖:
9楼: Originally posted by necturnal at 2012-05-22 08:35:41:
iop(5/13=1)的意思忽略不收敛的failure,继续算下去。
高斯官网有个IOp reference,
http://www.gaussian.com/g_tech/g_iops/iops_toc.htm
大部分IOp的意思可以在里面查到,但是还是不全。Sobereva大神有几个

昨天去掉了scfcyc=300,结果出现了9999的错误。看样子还是这个不能少。

我找找IOp的帖子,谢谢
业精于勤荒于嬉,成于思而毁于随
10楼2012-05-22 08:46:16
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 ellemiu 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 一志愿西安交通大学材料工程专业 282分求调剂 +10 枫桥ZL 2026-03-18 12/600 2026-03-21 22:02 by peike
[考研] 354求调剂 +6 Tyoumou 2026-03-18 9/450 2026-03-21 20:47 by lbsjt
[考研] 333求调剂 +5 87639 2026-03-21 7/350 2026-03-21 19:31 by ColorlessPI
[考研] 【考研调剂】化学专业 281分,一志愿四川大学,诚心求调剂 +11 吃吃吃才有意义 2026-03-19 11/550 2026-03-21 18:23 by 学员8dgXkO
[考研] 306求0703调剂一志愿华中师范 +5 纸鱼ly 2026-03-21 5/250 2026-03-21 17:11 by 学员8dgXkO
[考研] 268求调剂 +9 简单点0 2026-03-17 9/450 2026-03-21 15:37 by lature00
[考研] 279求调剂 +5 红衣隐官 2026-03-21 5/250 2026-03-21 14:59 by lature00
[考研] 279分求调剂 一志愿211 +14 chaojifeixia 2026-03-19 15/750 2026-03-21 13:24 by zhukairuo
[考研] 一志愿华中科技大学,080502,354分求调剂 +5 守候夕阳CF 2026-03-18 5/250 2026-03-21 01:06 by JourneyLucky
[考研] 304求调剂 +7 司空. 2026-03-18 7/350 2026-03-20 23:08 by JourneyLucky
[考研] 287求调剂 +7 晨昏线与星海 2026-03-19 8/400 2026-03-20 22:19 by JourneyLucky
[考研] 中南大学化学学硕337求调剂 +3 niko- 2026-03-19 6/300 2026-03-20 21:58 by luoyongfeng
[考研] 0817 化学工程 299分求调剂 有科研经历 有二区文章 +22 rare12345 2026-03-18 22/1100 2026-03-20 20:39 by zhukairuo
[考研] 求调剂 +3 eation27 2026-03-20 3/150 2026-03-20 19:32 by JourneyLucky
[考研] 环境工程调剂 +9 大可digkids 2026-03-16 9/450 2026-03-20 17:38 by 醉在风里
[考研] 344求调剂 +6 knight344 2026-03-16 7/350 2026-03-18 20:13 by walc
[考研] 0703化学调剂 +3 妮妮ninicgb 2026-03-17 3/150 2026-03-18 10:29 by macy2011
[考研] 一志愿南京大学,080500材料科学与工程,调剂 +4 Jy? 2026-03-16 4/200 2026-03-17 11:02 by gaoqiong
[考研] 070303 总分349求调剂 +3 LJY9966 2026-03-15 5/250 2026-03-16 14:24 by xwxstudy
[考研] 085601材料工程315分求调剂 +3 yang_0104 2026-03-15 3/150 2026-03-15 10:58 by peike
信息提示
请填处理意见