24小时热门版块排行榜    

查看: 2075  |  回复: 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 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 考研调剂 +3 呼呼?~+123456 2026-03-21 3/150 2026-03-21 20:04 by 无际的草原
[考研] 一志愿东华大学控制学硕320求调剂 +3 Grand777 2026-03-21 3/150 2026-03-21 19:23 by 简之-
[考研] 311求调剂 +3 勇敢的小吴 2026-03-20 3/150 2026-03-21 17:40 by ColorlessPI
[考研] 296求调剂 +4 www_q 2026-03-20 4/200 2026-03-21 17:26 by 学员8dgXkO
[考研] 一志愿重庆大学085700资源与环境总分308求调剂 +7 墨墨漠 2026-03-20 7/350 2026-03-21 16:36 by barlinike
[考研] 330求调剂0854 +3 assdll 2026-03-21 3/150 2026-03-21 13:01 by 搏击518
[考研] 070300化学319求调剂 +7 锦鲤0909 2026-03-17 7/350 2026-03-21 03:46 by JourneyLucky
[考研] 材料 336 求调剂 +3 An@. 2026-03-18 4/200 2026-03-21 01:39 by JourneyLucky
[考研] 一志愿武理材料305分求调剂 +6 想上岸的鲤鱼 2026-03-18 7/350 2026-03-21 01:03 by JourneyLucky
[考研] 一志愿华南师大 070300(化学)304分求调剂 +3 0703武芊慧雪304 2026-03-18 3/150 2026-03-21 00:48 by JourneyLucky
[考研] 294求调剂材料与化工专硕 +15 陌の森林 2026-03-18 15/750 2026-03-20 23:28 by JourneyLucky
[考研] 一志愿 西北大学 ,070300化学学硕,总分287,双非一本,求调剂。 +4 晨昏线与星海 2026-03-19 4/200 2026-03-20 22:15 by JourneyLucky
[考研] 材料与化工 322求调剂 +4 然11 2026-03-19 4/200 2026-03-20 22:12 by luoyongfeng
[考研] 329求调剂 +9 想上学吖吖 2026-03-19 9/450 2026-03-20 22:01 by luoyongfeng
[考研] 289求调剂 +6 怀瑾握瑜l 2026-03-20 6/300 2026-03-20 20:30 by 学员8dgXkO
[考研] 261求B区调剂,科研经历丰富 +3 牛奶很忙 2026-03-20 4/200 2026-03-20 19:34 by JourneyLucky
[考研] 081700化工学硕调剂 +3 【1】 2026-03-16 3/150 2026-03-19 23:40 by edmund7
[考研] 344求调剂 +6 knight344 2026-03-16 7/350 2026-03-18 20:13 by walc
[考研] 一志愿,福州大学材料专硕339分求调剂 +3 木子momo青争 2026-03-15 3/150 2026-03-17 07:52 by laoshidan
[考研] 321求调剂 +5 大米饭! 2026-03-15 5/250 2026-03-16 16:33 by houyaoxu
信息提示
请填处理意见