应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 第一性原理 » 【求助】Dmol3中DIIS size的用法
111/2返回列表12下一页
查看: 2427  |  回复: 10
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看
当前主题已经存档。

xiaowandouer

木虫 (正式写手)

[交流] 【求助】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 ]
回复此楼

» 收录本帖的淘帖专辑推荐

graphene

» 猜你喜欢
1楼 2009-05-26 14:44:18
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

xiaowandouer

木虫 (正式写手)
帖子沉得好快，顶上去，希望高手出现
赞一下 回复此楼
2楼2009-05-26 18:07:10
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

rongxue6808

新虫 (正式写手)
★ ★ ★ ★ ★
xiaowandouer(金币+3,VIP+0):谢谢您的热心应助，好运连连。。。。 5-26 20:19
qasd(金币+2,VIP+0):谢谢~ 5-26 20:44
你设的这个值是可以的。但一般是用默认值，除非你的scf收敛有问题，比如500圈scf迭代完了仍不收敛。DIIS 是“Direct Inversion of Iterative Subspace”的缩写，是用来加速收敛的，即前一次迭代产生的轨道信息会部分用于下一次迭代（类似一个restart file）,这是不同于第一次迭代的。第一迭代的原始信息是猜想出来（guess）。这也是为什么第一个scf收敛比较耗时而后面的scf收敛却越来越快的原因。
赞一下(12人) 回复此楼
3楼2009-05-26 19:03:50
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

rongxue6808

新虫 (正式写手)
★
xiaowandouer(金币+1,VIP+0):谢谢，你的提示对我很有帮助！ 5-26 20:20
Smear的功能也是关于收敛的。 在我现用的软件中stobe，DIIS和Smear是不能同时用的。所以在dmol3里,我一般会只改smear的值。纯属个人理解和经验，不一定正确哈
赞一下(1人) 回复此楼
4楼2009-05-26 19:13:53
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

cenwanglai

荣誉版主 (知名作家)

老和山猥琐派九段

优秀版主
★
xiaowandouer(金币+1,VIP+0):谢谢哈！ 5-26 20:20
帮你顶顶。呵呵。
回复此楼
5楼2009-05-26 19:53:19
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

xiaowandouer

木虫 (正式写手)
引用回帖:
Originally posted by rongxue6808 at 2009-5-26 19:03:
你设的这个值是可以的。但一般是用默认值，除非你的scf收敛有问题，比如500圈scf迭代完了仍不收敛。DIIS 是“Direct Inversion of Iterative Subspace”的缩写，是用来加速收敛的，即前一次迭代产生的轨道信息会部 ...

谢谢您的详细解释。

我现在遇到的就是scf不收敛的问题，屡次调smearing值仍是不奏效，所以才考虑用DIIS Size。
刚刚我又做了个测试，就是拿比较容易收敛的体系修改DIIS Size值重新算，结果得到的能量一样，加上您这么说，我现在就敢调这个参数了

另外，就是楼上有位虫友提到的，在Dmol3中DIIS Size和smearing能不能同时用，如果不能的话，为什么？会出现什么问题？
赞一下 回复此楼
6楼2009-05-26 20:17:39
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

rongxue6808

新虫 (正式写手)
★ ★
xiaowandouer(金币+1,VIP+0):谢谢，最后一分送给您！ 5-27 08:29
qasd(金币+1,VIP+0):谢谢~ 5-28 12:38
在Dmol3中DIIS Size和smearin可以同时用
赞一下(10人) 回复此楼
7楼2009-05-26 23:37:52
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

shrek826

金虫 (小有名气)
★
小木虫(金币+0.5):给个红包，谢谢回帖交流
我也来讨论一下：
我算某个体系的单点能，测试了四种情况：
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
赞一下(1人) 回复此楼
8楼2009-07-01 15:28:50
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

aaq2800

铁杆木虫 (知名作家)
★
小木虫(金币+0.5):给个红包，谢谢回帖交流
引用回帖:
Originally posted by shrek826 at 2009-7-1 15:28:
我也来讨论一下：
我算某个体系的单点能，测试了四种情况：
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， ...

你这个3和4的能量主要是由 smearing引起的
赞一下(1人) 回复此楼
9楼2009-07-01 20:50:56
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

117779747

金虫 (小有名气)
★
小木虫(金币+0.5):给个红包，谢谢回帖交流
引用回帖:
Originally posted by shrek826 at 2009-7-1 15:28:
我也来讨论一下：
我算某个体系的单点能，测试了四种情况：
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， ...

我现在算的话，一般先是拿DIIS和smearing同时用，如果收敛，则去掉smearing继续优化，这时得到的能量做最终计算用。
赞一下(1人) 回复此楼
10楼2009-07-02 08:13:58
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 xiaowandouer 的主题更新
111/2返回列表12下一页
普通表情 龙 兔 虎 猫 高级回复 (可上传附件)
信息提示
关闭
请填处理意见
关闭