24小时热门版块排行榜    

查看: 547  |  回复: 3

泽润东方

金虫 (正式写手)

[交流] 【求助】Gaussian 计算能量不能收敛

本人刚刚接触Gaussian计算。
想利用Gaussian03计算一个有机分子/金属基底的能量,总原子数少于100.但是总是不能收敛,计算选用的是B3LYP/LANL2MB,各位能否给在下解释一下不能收敛的原因,以及有几种解决方法可以解决这个问题 ?

拜谢。
回复此楼
希望就在不远处,成功就在路尽头。
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

manson1998

金虫 (著名写手)

OLED达人

★ ★
泽润东方(金币+5): 2010-09-06 19:31:28
luoqiquan(金币+2):thanks 2010-09-06 19:37:36
首先要看你优化结构的走向是否正确。如果不正确的话就要调整结构。一般根据文献或者实验测得的晶体结构入手。在输出文件中可以看到这一类不能收敛的立场都是远大于收敛值的。
如果结构基本正确的话,但是总不能收敛的话。一般可以通过缩小优化步长,使用小基组的方法来达到收敛。欢迎继续补充细节来讨论。
惺惺相惜情不自禁
2楼2010-09-06 17:51:12
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

manson1998

金虫 (著名写手)

OLED达人


dongdong3881(金币+1):多谢参与 2010-09-07 10:01:58
泽润东方(金币+5): 2010-09-09 18:21:04
泽润东方(金币+5):谢谢,小基组不适合我计算的体系里的金属元素。已经解决了。 2010-11-16 15:00:24
缩小步长可以加入关键词iop(1/8=10)默认为30貌似,小基组如果你原先用的是6-31g* 改为6-31G或者3-21G*都叫减小基组。这样的设置一般会达到收敛,但是精度不高。一般我们在小基组上收敛之后再用大基组优化,也能得到收敛的结果。再有一种方法是用QC 或者降低收敛精度。前者十分耗时,后者在个人看来没有太大意义。
惺惺相惜情不自禁
3楼2010-09-07 08:42:03
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

manson1998

金虫 (著名写手)

OLED达人

忘记补充一下 不合适的计算方法也可能导致收敛的失败
惺惺相惜情不自禁
4楼2010-09-07 08:45:43
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 泽润东方 的主题更新
普通表情 高级回复 (可上传附件)
信息提示
请填处理意见