计算化学公社»论坛首页 › 理论与计算化学 (Theoretical and Computational Chemistry) › 量子化学 (Quantum Chemistry) › GAMESS诡异的SCF收敛情况

GAMESS诡异的SCF收敛情况

查看数: 13225 | 评论数: 14 | 收藏 Add to favorites 3

关灯 | 提示：支持键盘翻页<-左右->

帖子模式

airyang

发布时间: 2017-8-26 17:02

正文摘要:

本帖最后由 airyang 于 2017-8-28 03:20 编辑楼主计算一个有C2对称性分子中两个片段之间的相互作用能量能量分解（GKS-EDA），为方便我们把片段1,2叫M1，M2。补充说明下的是，执行GKS-EDA的GAMESS是问苏培峰老师 ...

回复 Reply

sobereva 发表于 Post on 2018-8-1 05:20:56

yang_lzu 发表于 2018-7-31 09:32
明白了，非常感谢老师的回答，另外老师有推荐的关于描述泛函的书籍吗？对于wb97x-d了解得很少。

没必要看什么特殊文章。
非要看，直接看原文，以及看较新的泛函测试、综述文章

yang_lzu 发表于 Post on 2018-7-31 09:32:12

sobereva 发表于 2018-7-31 04:30
SCF收敛与否有一定数值巧合性，你当前的问题没有什么必然原因

明白了，非常感谢老师的回答

，另外老师有推荐的关于描述泛函的书籍吗

？对于wb97x-d了解得很少。

sobereva 发表于 Post on 2018-7-31 04:30:46

yang_lzu 发表于 2018-7-30 13:01
敢问大神，之前我用的M06-2X完美收敛，但后面换成wb97x-d就不能收敛了（在m06-2x基础上上改了方法，加了N ...

SCF收敛与否有一定数值巧合性，你当前的问题没有什么必然原因

yang_lzu 发表于 Post on 2018-7-30 13:01:21

wei 发表于 2017-8-29 12:09
我用的是wb97xd，这也是苏老师推荐的之一。至于m06-2x，你试试 NRAD=99 nthe=36 nphi=72试试。

敢问大神，之前我用的M06-2X完美收敛，但后面换成wb97x-d就不能收敛了（在m06-2x基础上上改了方法，加了NRAD=99 nthe=36 nphi=72），这是为什么呢？

airyang 发表于 Post on 2017-8-30 02:08:47

本帖最后由 airyang 于 2017-8-30 02:18 编辑

补充说明：
出错应该与切换为极坐标格点无关，且尝试后发现解决楼主的个例问题中也不必动用 NRAD=99 nthe=36 nphi=72 这么大的格点（GAMESS说明书中称为army grade），下面仍然使用Lebedev 格点：

径向的NRAD从99设置成100，角度NLEB仍然选用meta-GGA默认的590，就可以解决收敛错误的问题；
或者径向NRAD=99不变，角度选择GAMESS提供的下一档770也可以解决SCF错误收敛的问题。
很可能是GAMESS-2012自身的问题，因为同样结构在GAMESS-2014中全部正常收敛

airyang 发表于 Post on 2017-8-29 20:56:40

wei 发表于 2017-8-29 12:09
我用的是wb97xd，这也是苏老师推荐的之一。至于m06-2x，你试试 NRAD=99 nthe=36 nphi=72试试。

多谢啊，换了你建议的格点类型与精度，果然解决了！

wei 发表于 Post on 2017-8-29 12:09:29

我用的是wb97xd，这也是苏老师推荐的之一。至于m06-2x，你试试 NRAD=99 nthe=36 nphi=72试试。

airyang 发表于 Post on 2017-8-28 01:28:38

wei 发表于 2017-8-27 19:28
你遇到的情况我也遇到过，加大格点积分精度后都正常了。现在头疼的是，大体系做EDA，内存不够报错。（我的 ...

因为用了明尼苏达泛函，整个系列的积分格点一直使用等同gaussian的ultrafine级别（NRAD=99 NLEB=590），请问下你解决的时候用多大的积分格点呢？

EDA确实太吃内存，我64G瑟瑟发抖，啥都做不动 = =！

wei 发表于 Post on 2017-8-27 19:28:35

本帖最后由 wei 于 2017-8-27 19:33 编辑

你遇到的情况我也遇到过，加大格点积分精度后都正常了。现在头疼的是，大体系做EDA，内存不够报错。（我的内存是256G）无法解决。

airyang 发表于 Post on 2017-8-26 17:59:37

sobereva 发表于 2017-8-26 17:11
改成其它基组或泛函（特别是高HF成份的）能解决不？或者索性关了SOSCF，一直用DIIS
也可以尝试写上$DFT SW ...

汇报下测试结果：
1.M06-2X-D3/def-TZVP 出错
2.B3LYP-D3/def2-TZVPP 正常
3.原始的M06-2X-D3/def2-TZVPP 出错

airyang 发表于 Post on 2017-8-26 17:41:25

因为用了M06-2X, 格点是用的 NRAD=99 NLEB=590，等同gaussian的ultrafine

airyang 发表于 Post on 2017-8-26 17:38:55

sobereva 发表于 2017-8-26 17:11
改成其它基组或泛函（特别是高HF成份的）能解决不？或者索性关了SOSCF，一直用DIIS
也可以尝试写上$DFT SW ...

$SCF DIIS=.T. SOSCF=.F. FDIFF=.F. ...，一直用的DIIS。错误的例子尝试过关了DIIS，只用SOSCF结果很不好收敛，放宽收敛限后，可以进入M1的with all basisi set 但还是出一样的错误。

因为查了公社帖子：http://bbs.keinsci.com/thread-435-1-1.html
SWITCH=1.0 切积分格点的也试了，甚至设置到了100 (因为那个density change 会大到50，所以干脆阈值给100)，也没奏效。

airyang 发表于 Post on 2017-8-26 17:24:43

sobereva 发表于 2017-8-26 17:11
改成其它基组或泛函能解决不？

基组改小很多有用def2-tzvpp换3-21g，小一点儿（def2-tzvpp换def2-tzvp）没用，介于之间的没尝试，我试试

方法M06-2X没换过，因为是一系列的比较还是想着最终结果能用，我去试试B3LYP

一共将近70个做GKS-EDA的相似结构，三个结构SCF异常。

sobereva 发表于 Post on 2017-8-26 17:11:29

改成其它基组或泛函（特别是高HF成份的）能解决不？或者索性关了SOSCF，一直用DIIS
也可以尝试写上$DFT SWITCH=1.0 $END避免迭代前期用低精度积分格点看能否奏效

		自动登录 Automatic login	找回密码 Forget password
密码 Password			注册 Register