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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 量子化学 » Gaussian » 【求助】请教IRC计算出现的一个问题
北京石油化工学院2026年研究生招生接收调剂公告
101/1返回列表
查看: 977  |  回复: 9
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看
当前主题已经存档。

dongda986

超级版主

优秀！！有木有！！！优秀！！有木有！！！优秀！！有木有！！！优秀！！有木有！！！

[交流] 【求助】请教IRC计算出现的一个问题

大家好，小弟刚接触高斯不久，昨天做IRC计算时出现这样的结果，请问是怎么回事，该怎么处理呢？
输入文件
#b3lyp/6-311g** irc=(calcfc,maxpoints=15,maxcycle=200)   
                                                         
irc                                                      
                                                         
0  1                                                      
6                    1.480515   -0.269866   -0.000058
6                   -0.430759   -0.007681   -0.000244
8                   -0.507553    1.248671    0.000012
8                   -1.053900   -0.999675    0.000106
1                    1.600926   -0.865519    0.901495
1                    2.271501    0.503300    0.000376
1                    1.601541   -0.865907   -0.901295
1                    0.719116    0.901435    0.000291
得出结果的最后部分：
Convergence criterion not met.
SCF Done:  E(RB+HF-LYP) =  -211.543169048     A.U. after  129 cycles
             Convg  =    0.7611D-02             -V/T =  1.9483
             S**2   =   0.0000
Convergence failure -- run terminated.
Error termination via Lnk1e in C:\G03W\l502.exe at Fri Nov 06 21:39:04 2009.
Job cpu time:  0 days  0 hours 49 minutes 39.0 seconds.
File lengths (MBytes):  RWF=     21 Int=      0 D2E=      0 Chk=     13 Scr=      1

[ Last edited by dongda986 on 2009-11-7 at 12:43 ]
回复此楼

» 猜你喜欢
1楼 2009-11-07 12:15:15
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

erylingjet

实习版主

优秀！！有木有！！！优秀！！有木有！！！优秀！！有木有！！！优秀！！有木有！！！
不收敛,就是结构不好呗
赞一下 回复此楼
2楼2009-11-07 15:49:45
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

erylingjet

专家顾问

优秀！！有木有！！！优秀！！有木有！！！优秀！！有木有！！！优秀！！有木有！！！
过渡态找到了吗?
赞一下 回复此楼
3楼2009-11-07 15:49:59
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

dongda986

实习版主

优秀！！有木有！！！优秀！！有木有！！！优秀！！有木有！！！优秀！！有木有！！！
引用回帖:
Originally posted by erylingjet at 2009-11-7 15:49:
过渡态找到了吗?

过渡态找到了，请问该在哪里做调整啊。我后来又算了一遍，把 iop（5/13=1）加到输入文件里去，这次算出来30种物质，感觉已经差不多了，但还是没有算成功，没能得到反应坐标对应的能量值，请问我该怎么处理一下。
赞一下 回复此楼
4楼2009-11-07 18:59:41
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

3867826

版主

优秀！！有木有！！！优秀！！有木有！！！优秀！！有木有！！！优秀！！有木有！！！
★
dongda986(金币+1,VIP+0):谢谢 11-8 18:34
分别向前  向后扫描 看情况出在哪里
赞一下 回复此楼
拒稿，终于体会到了！乱写接受，认真写拒稿，这就是2和4的区别吗？
5楼2009-11-07 19:56:56
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

dongda986

管理员

优秀！！有木有！！！优秀！！有木有！！！优秀！！有木有！！！优秀！！有木有！！！
引用回帖:
Originally posted by 3867826 at 2009-11-7 19:56:
分别向前  向后扫描 看情况出在哪里

请问你说的是用scan扫描势能面吗，我刚接触高斯，请说具体点，最好帮我简单的写一下例子，谢谢
赞一下 回复此楼
6楼2009-11-07 20:45:15
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

majun04

主管区长

majun

优秀！！有木有！！！优秀！！有木有！！！优秀！！有木有！！！优秀！！有木有！！！
★ ★ ★ ★ ★ ★
yjcmwgk(金币+1,VIP+0): 11-8 15:51
dongda986(金币+5,VIP+0):谢谢，我这就试试看 11-8 18:34
沿着反应路径向前找：
# IRC=(calcfc,forward) rb3lyp/6-31+g(d) geom=allcheck guess=read
也就是得到ＩＲＣ路径图上坐标为正的点。
沿着反应路径向后找：
# IRC=(calcfc,reverse) rb3lyp/6-31+g(d) geom=allcheck guess=read
也就是得到ＩＲＣ路径图上坐标为负的点。

这两个例子需要读取你计算过渡态频率得到的check文件。

看看《手册》109页。

[ Last edited by majun04 on 2009-11-8 at 10:17 ]
赞一下(11人) 回复此楼
乐观、执着、勤奋
7楼2009-11-08 10:15:36
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

dongda986

兑换贵宾

优秀！！有木有！！！优秀！！有木有！！！优秀！！有木有！！！优秀！！有木有！！！
引用回帖:
Originally posted by majun04 at 2009-11-8 10:15:
沿着反应路径向前找：
# IRC=(calcfc,forward) rb3lyp/6-31+g(d) geom=allcheck guess=read
也就是得到ＩＲＣ路径图上坐标为正的点。
沿着反应路径向后找：
# IRC=(calcfc,reverse) rb3lyp/6-31+g(d) geom=al ...

你好，我想再请教一下，我刚把向前和向后都算了一下，都没有发生错误，但是两次反应坐标零点对应的能量值相差较大，分别是 -229.02786和-229.07551，请问这是怎么回事，这样的结果说明我的计算在哪里出了问题呢？
赞一下 回复此楼
8楼2009-11-08 21:08:37
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

bylin

版主

优秀！！有木有！！！优秀！！有木有！！！优秀！！有木有！！！优秀！！有木有！！！
★
小木虫(金币+0.5):给个红包，谢谢回帖交流
两次计算所用的check文件都是原始的过渡态的check文件吗？如果第二次计算用的是第一次算完之后的check文件，那就不对了。
如果没有上面的问题，两次计算所用基组、方法也完全相同，两个零点坐标的能量应该是一样的。
赞一下(1人) 回复此楼
9楼2009-11-09 09:03:03
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

dongda986

实习版主

优秀！！有木有！！！优秀！！有木有！！！优秀！！有木有！！！优秀！！有木有！！！
引用回帖:
Originally posted by bylin at 2009-11-9 09:03:
两次计算所用的check文件都是原始的过渡态的check文件吗？如果第二次计算用的是第一次算完之后的check文件，那就不对了。
如果没有上面的问题，两次计算所用基组、方法也完全相同，两个零点坐标的能量应该是一样 ...

哦，我的确是没注意那个细节，多谢啊
赞一下 回复此楼
10楼2009-11-09 09:54:05
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 dongda986 的主题更新
101/1返回列表
普通表情 龙 兔 虎 猫 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 0856，材料与化工321分求调剂 +13 大馋小子 2026-03-27 14/700 2026-04-01 15:58 by yanflower7133
[考研] 317求调剂 +14 蛋黄咸肉粽 2026-03-26 15/750 2026-04-01 15:57 by BruceLiu320
[考研] 08工科275求调剂，可跨考。 +5 AaAa7420 2026-03-31 5/250 2026-04-01 15:21 by 159357hjz
[考研] 318求调剂 +8 七忆77 2026-04-01 8/400 2026-04-01 10:37 by Jaylen.
[考研] 环境工程 085701，267求调剂 +15 minht 2026-03-29 16/800 2026-04-01 10:13 by li_sujuan99
[考研] 土木304求调剂 +3 兔突突突， 2026-03-31 3/150 2026-04-01 09:42 by JourneyLucky
[考研] 一志愿西安交大材料学硕（英一数二）347，求调剂到高分子/材料相关专业 +6 zju51 2026-03-31 8/400 2026-04-01 09:29 by JourneyLucky
[考研] 考研调剂 +9 小蜡新笔 2026-03-29 10/500 2026-03-31 19:52 by Dyhoer
[考研] 085600 295分求调剂 +13 W55j 2026-03-30 15/750 2026-03-31 13:29 by 王亮_大连医科大
[考研] 276求调剂 +3 赵久华 2026-03-29 3/150 2026-03-31 10:06 by cal0306
[考研] 085600材料与化工329分求调剂 +5 这是有太阳哇 2026-03-25 5/250 2026-03-31 09:59 by luoyongfeng
[考研] 285求调剂 +6 AZMK 2026-03-29 9/450 2026-03-30 21:02 by dophin1985
[考研] 310求调剂 +10 争取九点睡 2026-03-30 10/500 2026-03-30 16:45 by ztnimte
[考研] 考研调剂 +5 Sanmu-124 2026-03-26 5/250 2026-03-30 13:36 by chemdavid
[考研] 085600，材料与化工321分求调剂 +10 大馋小子 2026-03-28 10/500 2026-03-29 23:35 by 飞行日记西
[考研] 332求92调剂 +8 蕉蕉123 2026-03-28 8/400 2026-03-29 10:46 by 周梓丹
[考研] 316求调剂 +7 江辞666 2026-03-26 7/350 2026-03-28 21:28 by sanrepian
[考研] 材料与化工（0856）304求B区调剂 +8 邱gl 2026-03-27 8/400 2026-03-28 12:42 by 唐沐儿
[考研] 考研调剂 +9 小蜡新笔 2026-03-26 9/450 2026-03-27 11:10 by 不吃魚的貓
[考研] 081200-11408-276学硕求调剂 +4 崔wj 2026-03-26 4/200 2026-03-27 08:04 by chemisry
信息提示
关闭
请填处理意见
关闭