24小时热门版块排行榜    

查看: 918  |  回复: 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的回帖

meteoric30

木虫 (正式写手)

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

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的回帖
相关版块跳转 我要订阅楼主 wurongyan 的主题更新
普通表情 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 0817 化学工程 299分求调剂 有科研经历 有二区文章 +14 rare12345 2026-03-18 14/700 2026-03-19 19:45 by maocaozhuxi
[考研] 一志愿南京理工大学085701资源与环境302分求调剂 +3 葵梓卫队 2026-03-18 5/250 2026-03-19 19:35 by 给你你注意休息
[考研] 材料专硕274一志愿陕西师范大学求调剂 +8 薛云鹏 2026-03-13 8/400 2026-03-19 15:36 by haoshis
[考研] 0703化学调剂,求各位老师收留 +10 秋有木北 2026-03-14 10/500 2026-03-19 05:52 by anny19840123
[考研] 344求调剂 +6 knight344 2026-03-16 7/350 2026-03-18 20:13 by walc
[考研] 26调剂/材料/英一数二/总分289/已过A区线 +7 步川酷紫123 2026-03-13 7/350 2026-03-18 17:12 by 尽舜尧1
[考研] 材料,纺织,生物(0856、0710),化学招生啦 +3 Eember. 2026-03-17 9/450 2026-03-18 10:28 by Eember.
[考研] 材料与化工专硕调剂 +5 heming3743 2026-03-16 5/250 2026-03-17 14:03 by 勇敢太监王公公
[考研] 302求调剂 +4 小贾同学123 2026-03-15 8/400 2026-03-17 10:33 by 小贾同学123
[考研] 一志愿,福州大学材料专硕339分求调剂 +3 木子momo青争 2026-03-15 3/150 2026-03-17 07:52 by laoshidan
[考研] 274求调剂 +5 时间点 2026-03-13 5/250 2026-03-17 07:34 by 热情沙漠
[考研] 11408 一志愿西电,277分求调剂 +3 zhouzhen654 2026-03-16 3/150 2026-03-17 07:03 by laoshidan
[考研] 药学383 求调剂 +3 药学chy 2026-03-15 4/200 2026-03-16 20:51 by 元子^0^
[考研] 22408总分284求调剂 +3 InAspic 2026-03-13 3/150 2026-03-15 11:10 by zhq0425
[考研] 288求调剂 +4 奇点0314 2026-03-14 4/200 2026-03-14 23:04 by JourneyLucky
[基金申请] 现在如何回避去年的某一个专家,不知道名字 +3 zk200107 2026-03-12 6/300 2026-03-14 17:13 by zk200107
[考研] 297求调剂 +4 学海漂泊 2026-03-13 4/200 2026-03-14 11:51 by 热情沙漠
[考研] 26调剂/材料科学与工程/总分295/求收留 +9 2026调剂侠 2026-03-12 9/450 2026-03-13 20:46 by 18595523086
[考研] 308求调剂 +3 是Lupa啊 2026-03-12 3/150 2026-03-13 14:30 by 求调剂zz
[考研] 290求调剂 +3 ADT 2026-03-13 3/150 2026-03-13 10:19 by peike
信息提示
请填处理意见