应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 量子化学 » Gaussian » IRC 结果分析
311/4返回列表1234下一页
查看: 2495  |  回复: 30
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看

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扫描完以后的结果
回复此楼

» 收录本帖的淘帖专辑推荐

过渡态

» 猜你喜欢

» 本主题相关价值贴推荐，对您同样有帮助:
业精于勤荒于嬉，成于思而毁于随
1楼 2012-05-18 14:59:37
已阅   回复此楼   关注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）吗？
赞一下(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的回帖

Single_love
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 的主题更新
311/4返回列表1234下一页
普通表情 龙 兔 虎 猫
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[博后之家] 售SCI一区文章，我:8O5.5.1.O5.4,科目全,可伽急 +4 h4CP7TrQR8Lg 2026-08-22 4/200 2026-08-23 20:46 by 0p9y1isJp2jM
[论文投稿] 售SCI文章，我:8O.5.5.1O.54,科目全,可十急 +3 DpSrDtM079iu 2026-08-22 3/150 2026-08-23 15:33 by KM5EcsNQRBPn
[论文投稿] 售SCI一区T0P文章，我:8.O.55.1.O54,科目全,可伽急 +3 DpSrDtM079iu 2026-08-22 4/200 2026-08-23 15:30 by KM5EcsNQRBPn
[考博] 售SCI一区T0P文章，我:8.O.55.1.O.54,科目齐全,可+急 +3 h4CP7TrQR8Lg 2026-08-22 5/250 2026-08-23 12:09 by LR9qGULyN2ew
[基金申请] 什么时候开奖？ +9 CrisMessi 2026-08-18 10/500 2026-08-23 12:03 by 丶昵称占用
[硕博家园] 售SCI一区T0P文章，我:8.O.55.1.O54,科目全,可伽急 +3 2JOx3r2CYEgw 2026-08-22 5/250 2026-08-23 07:41 by OEbVnUOu01ol
[博后之家] 售SCI一区T0P文章，我:8.O.55.1.O54,科目全,可伽急 +3 DpSrDtM079iu 2026-08-22 5/250 2026-08-23 04:07 by OEbVnUOu01ol
[论文投稿] 售SCI一区文章，我:8O5.5.1.O5.4,科目全,可伽急 +3 2JOx3r2CYEgw 2026-08-22 4/200 2026-08-23 04:04 by OEbVnUOu01ol
[硕博家园] 售SCI一区T0P文章，我:8.O.55.1.O.5.4,科目全,可+急 +3 2JOx3r2CYEgw 2026-08-22 7/350 2026-08-23 03:55 by OEbVnUOu01ol
[硕博家园] 售SCI一区文章，我:8O5.5.1.O5.4,科目全,可伽急 +3 2JOx3r2CYEgw 2026-08-22 7/350 2026-08-23 03:43 by OEbVnUOu01ol
[硕博家园] 售SCI一区T0P文章，我:8O.55.1.O.54,科目全,可伽急 +4 QTy3jDtz1uLt 2026-08-21 11/550 2026-08-23 02:55 by OEbVnUOu01ol
[基金申请] 93BebMhtakh前后11位开头都是大写 +7 且听虎啸 2026-08-17 8/400 2026-08-22 21:55 by 医学老男孩
[基金申请] filecode，4个jtjc了 +13 ziyangfang 2026-08-19 16/800 2026-08-22 17:08 by WH3796
[教师之家] 售一区SCI文章T0P，我:8O.551.O54,科目全,可十急 +4 2JOx3r2CYEgw 2026-08-21 4/200 2026-08-22 16:52 by sunzitan
[基金申请] 人气不行了 +8 fansofjerry 2026-08-21 8/400 2026-08-22 16:30 by zyqchem
[基金申请] 今日不放榜？网传国自然预计 8 月 27 日可查结果 +16 医学老男孩 2026-08-20 20/1000 2026-08-21 21:13 by Ldrop2023
[基金申请] 科研孤儿太难了 +17 我4大白菜 2026-08-20 18/900 2026-08-21 20:57 by zhangev
[基金申请] 看来今天不会放榜了？ +8 chengyan1220 2026-08-21 11/550 2026-08-21 17:52 by dcqxinyang
[论文投稿] 投稿咨询 +5 wwm09 2026-08-17 7/350 2026-08-21 10:11 by 期刊论文帮手
[基金申请] 估计是周四 +3 archvillain 2026-08-18 3/150 2026-08-21 01:48 by jnhyjjm
信息提示
关闭
请填处理意见
关闭