24小时热门版块排行榜    

查看: 1293  |  回复: 7
当前主题已经存档。
当前只显示满足指定条件的回帖,点击这里查看本话题的所有回帖

herosmzd

银虫 (正式写手)

[交流] 【求助】关于提高精度就不收敛的问题

算一个优化结构时,先用hf/3-21g,正常的,再用b3lyp/3-21g,就有2个不收敛
Maximum Force            0.000002     0.000450     YES
RMS     Force            0.000000     0.000300     YES
Maximum Displacement     0.011336     0.001800     NO
RMS     Displacement     0.001408     0.001200     NO
这是什么原因?我应该怎么做呢?谢谢高手指点我
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

hairan

木虫 (著名写手)

瞌睡木虫

★ ★
lei0736(金币+2,VIP+0):多谢详细解答
引用回帖:
Originally posted by herosmzd at 2008-3-3 09:02:
算一个优化结构时,先用hf/3-21g,正常的,再用b3lyp/3-21g,就有2个不收敛
Maximum Force            0.000002     0.000450     YES
RMS     Force            0.000000     0.000300     YES
Maximum Displa ...

如果按照Gaussian默认的收敛标准,你的结果应该收敛了,因为当force小于0.000004, RMS force小于0.000003的时候,程序会认为force已经可以忽略不计,已经找到了stationary point。你应该是使用了scf=tight一类的关键词
我的计算经验是,遇到此种情况,先把计算停了,然后将结构调出来,重新创建输入文件,开新的计算,输入文件的关键词中不要加入guess=read或者geom=allcheck等从chk文件读内容的选项,有时后两项会很快收敛的,但并不是任何时候都有效。
据权威报纸调查,84%的青少年犯罪者有玩电子游戏的经历,故应该禁止经营电子游戏。据我们所知,100%的青少年犯罪者有穿鞋的经历,故制鞋厂应该停业。
8楼2008-03-05 15:46:45
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 8 个回答

codyypj

银虫 (小有名气)

★ ★ ★
herosmzd(金币+2,VIP+0):多谢帮助
zzgyb(金币+1,VIP+0):谢谢你的帮助源,欢迎继续关注计算模拟版!
一般如果结构合理后面的两个关系不大。不收敛和没关系,你也可以试试scf=tight
实在还没收敛就:SCF=VERYTIGHT
2楼2008-03-03 09:39:13
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

好好生活

木虫 (著名写手)

★ ★ ★ ★
herosmzd(金币+2,VIP+0):谢谢您的意见
zzgyb(金币+2,VIP+0):谢谢你的帮助源,欢迎继续关注计算模拟版!
用 scf=xqc 更好,我试过很多,一般用 scf=xqc比用scf=tight 成功几率更大

[ Last edited by 好好生活 on 2008-3-3 at 10:32 ]
3楼2008-03-03 10:22:45
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

K372

铜虫 (小有名气)

★ ★ ★
lei0736(金币+2,VIP+0):多谢详细解答
herosmzd(金币+1,VIP+0):多谢指教,能帮我解答下对结果的影响吗?
你的FORCE项已经很小了,但是DISPLACEMENT项还很大,即势能梯度很小、优化位移很大,说明势能面非常平坦。这种情况很可能是由于DFT对势能面的描述定性不正确造成的。
1楼的朋友说的是不对的,SCF=TIGHT/VERYTIGHT只能提高HF-SCF的收敛精度,对DFT的优化过程特别是优化的收敛性没有直接的帮助。
由于DFT方法和传统的从头计算法出发点不同,当DFT方法描述的势能面出现问题的情况下,建议在计算条件许可时考虑尝试用Post-HF方法尝试进行优化,如原子数不超过50个,可以尝试采用MP2方法,不超过20个,可以尝试采用CCSD方法。由于这些方法对计算机性能要求较高,望LZ慎重考虑!
4楼2008-03-03 13:00:22
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
普通表情 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 关于豆爷回答的JTJC与%2F数量 +5 yang182083 2026-08-06 7/350 2026-08-08 18:29 by zhanghaozhu
[基金申请] 2026国自然放榜时间 +3 布布和一二 2026-08-08 3/150 2026-08-08 18:15 by Lightingo
[基金申请] 国基金的申报应该改成非等额制,评价高的钱多评价低的钱少,但是增加资助率 +7 a089 2026-08-07 7/350 2026-08-08 18:05 by gltch
[基金申请] 好奇怪的filecode +3 布布和一二 2026-08-08 4/200 2026-08-08 17:46 by zhanghaozhu
[考研] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +3 qTvzQBRAHjCi 2026-08-07 5/250 2026-08-08 16:47 by oEVWOejN9taj
[论文投稿] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +3 qTvzQBRAHjCi 2026-08-07 4/200 2026-08-08 16:47 by oEVWOejN9taj
[考博] 售SCI一区T0P文章,我:8O.55.1.O.54,科目全,可伽急 +3 Vi50GxzrFcSG 2026-08-07 4/200 2026-08-08 16:27 by oEVWOejN9taj
[基金申请] 售SCI文章,我:8O.5.5.1O.54,科目全,可十急 +3 Vi50GxzrFcSG 2026-08-07 5/250 2026-08-08 16:27 by oEVWOejN9taj
[基金申请] 这样的filecode谁见过 +5 布布和一二 2026-08-08 8/400 2026-08-08 14:54 by archvillain
[教师之家] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +4 HEQlVqMTIA7d 2026-08-07 5/250 2026-08-08 14:27 by oEVWOejN9taj
[论文投稿] 售SCI一区T0P文章,我:8O.55.1.O.5.4,科目齐全,可+急 +3 HEQlVqMTIA7d 2026-08-07 5/250 2026-08-08 14:22 by oEVWOejN9taj
[基金申请] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +4 KXLV3nuBVBY7 2026-08-07 5/250 2026-08-08 14:07 by oEVWOejN9taj
[教师之家] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +3 KXLV3nuBVBY7 2026-08-07 5/250 2026-08-08 14:07 by oEVWOejN9taj
[论文投稿] 售SCI文章,我:8O.5.5.1O.54,科目全,可十急 +3 KXLV3nuBVBY7 2026-08-07 6/300 2026-08-08 14:01 by oEVWOejN9taj
[基金申请] 关于代码变化问题,想知道的进来 +11 且听虎啸 2026-08-07 14/700 2026-08-08 12:50 by 布布和一二
[硕博家园] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +3 Vi50GxzrFcSG 2026-08-07 6/300 2026-08-08 10:19 by 3OOjAIS77qg2
[基金申请] 求各位大神看下 100+6 hpkpkpkp 2026-08-05 33/1650 2026-08-06 14:49 by zhiyanjiang
[基金申请] 8月时间戳变的,举个手。玩一下,释放压力 +9 archvillain 2026-08-04 11/550 2026-08-05 20:06 by wlwhappy
[基金申请] 面上提前没消息,有中的吗 +14 archvillain 2026-08-02 18/900 2026-08-04 14:42 by archvillain
[基金申请] 什么时候能放榜呀? +3 Jacob678 2026-08-03 3/150 2026-08-03 16:14 by gltch
信息提示
请填处理意见