24小时热门版块排行榜    

查看: 521  |  回复: 2

798879443

金虫 (正式写手)

[求助] 优化不收敛

虫友们好,
     我对Co3(dpa)4(ncs)2进行优化,使用的是b3lyp方法,二重态和四重态都可以算出来,但是六重态总是收敛失败,我后来降低了收敛精度,即使用conver=6,请问这样处理可以保证其收敛吗,还有与二、四重态还有可比性吗,因为我要找出能量最低的多重度,进行加电场处理,请大神指教一下。
回复此楼
蓦然回首,那人却在灯火阑珊处。
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

maxwell妖

新虫 (初入文坛)

是第三项NO吗?改一下manual displacement
2楼2013-04-02 20:24:09
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

zhouwohua

银虫 (小有名气)

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★
798879443: 金币+10, ★★★★★最佳答案, 辛苦了,谢啦 2013-04-19 19:41:10
计算外加电场对Co3(dpa)4(NCS)2结构影响时,本来Co3(dpa)4(NCS)2从晶体数据库取出来的结构是对称的,但NCS与Co3不在一条直线上,,我用bp86方法计算时加了nosymm关键词,在加小电场是,Co3(dpa)4(NCS)2变成对称的直线型了,随着电场加大,Co3(dpa)4(NCS)2变成了非对称性,请问是不是bp86方法得到的结果不好,随后我又用b3lyp计算了,不过没加nosymm关键词,请问如果比较方法的可行性,这两种处理有可比性,也就是nosymm关键词对能量等有影响吗,请大牛指点,本人新人一枚。
两者方法的关键词为:
#p opt nosymm  bp86 genecp scf=maxcycles=500
#p opt b3lyp/genecp iop(2/15=1)

这里的注意一下关键字:IOp(2/15)

Symmetry control.
-1               Turns on symmetry; same as 0 for molecules but turns on assignment of space group operations for PBC.
0                         Leave symmetry in whatever state it is presently in.
1                         Unconditionally turn symmetry off. Note that symmetry is still called, and will determine the framework group. However,                    the molecule is not oriented.
2                         Bring the molecule to a symmetry orientation, but then disable further use of symmetry.
3                         Don't even call Symm.
4                         Call Symm once with loose cutoffs, symmetrize the resulting coordinates then confirm symmetry with tight cutoffs.
5                                Recover the previous symmetry operations from the RWF, and confirm that the new structure has the same symmetry.
6                         Same as 5, but get symmetry info from the checkpoint.
00                       Default (10).
10                       Do re-orientation for PBC.
20                       Suppress re-orientation for PBC.
100                     Turn on symmetry operations for PBC.

如果你是想设置原子在同一直线上应该使用冗余坐标修改关键字。opt=modredundant

收敛指标修改了当然就和前面没有对比性了。起码只有在10的-6的级别上才有可靠的对比性。可以试下用小基组算,然后再用大基组算。如果你优化发现能量变化只是在10的-6精度内,就可以降低收敛指标,等收敛后再用之前相同的指标计算。
学无止境;天天学习;努力向上。
3楼2013-04-19 19:37:20
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 798879443 的主题更新
信息提示
请填处理意见