24小时热门版块排行榜    

查看: 849  |  回复: 4
当前只显示满足指定条件的回帖,点击这里查看本话题的所有回帖

wurongyan

木虫 (小有名气)

[交流] 【求助】请大家帮我分析下程序 已有4人参与

我用ub3lyp/6-31g(d) 计算CBrF3+正离子能量,可结构优化总是不收敛,怎么办呢?我尝试了加大循环次数(maxcycle=200),也尝试了用相同结构的分子优化所得结构作为初始结构,结果还是没有收敛。请各位高手给予指点,谢谢!!

输入文件:
# opt freq ub3lyp/6-31g(d) geom=connectivity

Title Card Required

1 2
C
F                  1              B1
F                  1              B2    2              A1
F                  1              B3    3              A2    2              D1
Br                 1              B4    3              A3    2              D2

   B1             1.35000000
   B2             1.35000000
   B3             1.35000000
   B4             1.91000000
   A1           109.47121829
   A2           109.47121829
   A3           109.47123134
   D1          -120.00003407
   D2           119.99998297

1 2 1.0 3 1.0 4 1.0 5 1.0
2
3
4
5

输出文件:
最后显示:
LOGIC IS A SYSTEMATIC METHOD OF COMING TO THE
WRONG CONCLUSION WITH CONFIDENCE.
Error termination request processed by link 9999.
Error termination via Lnk1e in C:\G03W\l9999.exe at Tue Sep 14 12:32:52 2010.
Job cpu time:  0 days  0 hours 13 minutes 39.0 seconds.
File lengths (MBytes):  RWF=     22 Int=      0 D2E=      0 Chk=      9 Scr=      1
输出文件:
http://g.zhubajie.com/urllink.php?id=9936328ad7kmgt1qyq4qcaf

[ Last edited by wurongyan on 2010-9-18 at 10:14 ]
回复此楼

» 猜你喜欢

» 本主题相关商家推荐: (我也要在这里推广)

» 本主题相关价值贴推荐,对您同样有帮助:

怎么样才能将知识变成智慧呢?
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

manson1998

金虫 (著名写手)

OLED达人

★ ★
小木虫(金币+0.5):给个红包,谢谢回帖交流
heyo_123(金币+1):鼓励交流~!~ 2010-09-19 09:36:49
这么小的体系,可以用QC。或者取出optmization中的倒数第二个点计算一下力常数,然后再优化。最重要的还是看一下结构有没有问题。
惺惺相惜情不自禁
3楼2010-09-19 09:22:18
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 5 个回答

meteoric30

木虫 (正式写手)

★ ★ ★
小木虫(金币+0.5):给个红包,谢谢回帖交流
小红豆(金币+2):感谢交流 2010-09-18 17:42:23
l9999的错。。。。不算错吧,继续算呢?要不就是能垒跳跃了。
2楼2010-09-18 15:50:29
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

zhangmt

至尊木虫 (著名写手)

我叫MT


小木虫(金币+0.5):给个红包,谢谢回帖交流
1、下载了你的输出文件,但用gv5打开失败,无法读取中间过程,因为手边没有gv4,所以简单的用文本软件看了看结果,关键出错信息为:
Optimization stopped.
    -- Number of steps exceeded,  NStep=  20
就是说超过20个循环而没有收敛,优化失败,终止。
把所有的SCF Done的能量作图,显示能量是显著下降趋势的,所以感觉如果加上opt(maxcycle=100)之类的,也许能收敛的。
2、用你的结构用g09跑,7步之内优化结束,只用了不到5分钟。但结果有一个-10的小虚频,所以按虚频方向微调结构,结果6步之内收敛。
3、但从振动虚频发现,分子的原有对称性已经被打破,仔细看g09的log文件发现,除了第一次输入结构系统判定为C3v,随后的优化一律自动改为C1——也就是忽略了对称性——这在g03中似乎是很罕见的——g03默认维持原有对称性,回去再看楼主的log文件,果然g03始终保持C3v的对称性——也许这也是03和09的区别之一吧。
4、虽然09的结果变成了C1,也就是说三个C-F键长已经不一样了,但相差细微,所以我认为这不过是基组的误差引起的而已。
以上供楼主参考。浪费我将近一个小时的时间。呵呵。
一群自以为正义凛然的年轻人将一切不能以科学解释的事情定性为封建迷信并大刀阔斧地进行消灭,其实这是修养不足学识浅薄的一种体现,也是可恶的偏执和愚蠢的自以
4楼2010-09-19 21:27:16
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

wurongyan

木虫 (小有名气)

谢谢!

我改变了输入文件的初始结构,已经搞定了!
谢谢大家的帮助!祝各位中秋快乐,天天有个好心情!!
怎么样才能将知识变成智慧呢?
5楼2010-09-20 10:09:17
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
普通表情 高级回复 (可上传附件)
信息提示
请填处理意见