【求助】Dmol3中DIIS size的用法
最近优化一个构型老不收敛，试着调了下DIIS size，然后在此构型基础上又用dmol3默认的参数优化，这样可以收敛，得到能量E1，但是不知道算出来的能量对不对，因为我拿这个E1跟用smearing收敛得到的能量E2相比，E1比E2大了两个hartree，所以我有点怀疑DIIS size是不是不能随便动。
手册上给的说明是：
DIIS size: Specify the maximum size of the subspace for the DIIS procedure. If the SCF does not converge with the default number of histories, increasing this value can sometimes lead to significantly improved SCF convergence. It is not recommended to use fewer than 4 histories. Allowed values = 1 to 10.
默认值为6，我用的是8
麻烦用过的高手指点下，我这样做对不对。因为基础比较薄弱，对“the subspace for the DIIS procedure”这类概念也理解不到位，所以只要相关，您尽可以畅所欲言。
[ Last edited by freshgirl on 2009-6-20 at 12:59 ]
返回小木虫查看更多
今日热帖
京公网安备 11010802022153号
帖子沉得好快，顶上去，希望高手出现
你设的这个值是可以的。但一般是用默认值，除非你的scf收敛有问题，比如500圈scf迭代完了仍不收敛。DIIS 是“Direct Inversion of Iterative Subspace”的缩写，是用来加速收敛的，即前一次迭代产生的轨道信息会部分用于下一次迭代（类似一个restart file）,这是不同于第一次迭代的。第一迭代的原始信息是猜想出来（guess）。这也是为什么第一个scf收敛比较耗时而后面的scf收敛却越来越快的原因。
Smear的功能也是关于收敛的。 在我现用的软件中stobe，DIIS和Smear是不能同时用的。所以在dmol3里,我一般会只改smear的值。纯属个人理解和经验，不一定正确哈
帮你顶顶。呵呵。
谢谢您的详细解释。
我现在遇到的就是scf不收敛的问题，屡次调smearing值仍是不奏效，所以才考虑用DIIS Size。
刚刚我又做了个测试，就是拿比较容易收敛的体系修改DIIS Size值重新算，结果得到的能量一样，加上您这么说，我现在就敢调这个参数了
另外，就是楼上有位虫友提到的，在Dmol3中DIIS Size和smearing能不能同时用，如果不能的话，为什么？会出现什么问题
，
在Dmol3中DIIS Size和smearin可以同时用
我也来讨论一下：
我算某个体系的单点能，测试了四种情况：
1.no DIIS, no smearing
2.no DIIS, smearing=0.005Ha
3.DIIS=6, no smearing
4.DIIS=6, smearing=0.005Ha
结果发现1，2算了500步还不收敛，而3，4算了几十步就能收敛
这样看来，调DIIS可能比smearing有用
另外，3和4计算出来的能量有点差别，大约相差0.02Ha