24小时热门版块排行榜    

CyRhmU.jpeg
查看: 1674  |  回复: 11
当前主题已经存档。

jghe

金虫 (著名写手)

[交流] 【求助】Dmol3不收敛策略

Dmol3结构优化和计算频率时不收敛有什么好办法?
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
回帖支持 ( 显示支持度最高的前 50 名 )

lwg020810

金虫 (小有名气)

★ ★ ★ ★ ★ ★
wuchenwf(金币+5,VIP+0):谢谢回帖交流 12-3 22:06
jghe(金币+1,VIP+0):thanks 12-4 18:22
做结构优化的时候(1)从粗糙做到精细,分步进行要容易收敛;(2)从小基组到大基组,当然如果你设粗糙或精细的时候基组自动选择了;(3)选择拖尾,DIIS,将SCF收敛标准调大,迭代次数增大,但不建议设到1000,MS培训班上说的是通常50,200足够,1000没意义。
DMol3算频率是很慢的,它是对体系的每一个原子进行频率计算,所以耗的时间比较长(这是我对频率的一点体会,我没有完整的做过频率,我几次对纳米管算频率算好几天都只算几个原子逼的我直接放弃)。
一点体会,请各高手指正
2楼2009-12-03 19:49:52
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
普通回帖

archhector

引用回帖:
Originally posted by lwg020810 at 2009-12-3 19:49:
做结构优化的时候(1)从粗糙做到精细,分步进行要容易收敛;(2)从小基组到大基组,当然如果你设粗糙或精细的时候基组自动选择了;(3)选择拖尾,DIIS,将SCF收敛标准调大,迭代次数增大,但不建议设到1000,MS培训 ...

结构优化的时候还能看到收敛趋势,问题是有的scf始终无法收敛,3 4百步都不行,这咋整啊。。。
3楼2009-12-03 21:02:35
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

lwg020810

金虫 (小有名气)

★ ★ ★
wuchenwf(金币+2,VIP+0):谢谢回帖交流 12-3 22:06
jghe(金币+1,VIP+0):thanks 12-4 18:23
那说明你建立的模型和理想模型相差太远,比如键长,如果键长相差太大的话会影响到收敛,把键长设为与文献值差不多大小就要容易收敛
4楼2009-12-03 21:40:59
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

archhector

引用回帖:
Originally posted by lwg020810 at 2009-12-3 21:40:
那说明你建立的模型和理想模型相差太远,比如键长,如果键长相差太大的话会影响到收敛,把键长设为与文献值差不多大小就要容易收敛

我就用的它自带的氧化铁模型,建了一个9层左右的slab,开放表面一层驰豫,铁的spin全设成high,而且是精度中等,scf就是收敛不了,唉。莫非要先用它原始的氧化铁跑castep优化原胞结构,然后再来dmol3驰豫?
5楼2009-12-03 21:57:13
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

lwg020810

金虫 (小有名气)

引用回帖:
Originally posted by archhector at 2009-12-3 21:57:

我就用的它自带的氧化铁模型,建了一个9层左右的slab,开放表面一层驰豫,铁的spin全设成high,而且是精度中等,scf就是收敛不了,唉。莫非要先用它原始的氧化铁跑castep优化原胞结构,然后再来dmol3驰豫?

我一直用dmol3对castep不了解啊,也许你尝试下采用coast,没有忘记用赝势吧,这个还是要高手来解答,我旁观
6楼2009-12-03 22:16:03
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

archhector

引用回帖:
Originally posted by lwg020810 at 2009-12-3 22:16:

我一直用dmol3对castep不了解啊,也许你尝试下采用coast,没有忘记用赝势吧,这个还是要高手来解答,我旁观

恩。。。我就是在dmol下跑这个模型scf始终收敛不了,残念了,只能harris近似了
7楼2009-12-03 22:23:24
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

lwg020810

金虫 (小有名气)


aylayl08(金币+1,VIP+0):谢谢参与 12-5 10:00
引用回帖:
Originally posted by archhector at 2009-12-3 22:23:

恩。。。我就是在dmol下跑这个模型scf始终收敛不了,残念了,只能harris近似了

你是不是同时计算频率或其他properties了?如果是的话建议先取消,等优化做最后一步的时候再计算性质
8楼2009-12-03 22:31:34
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

archhector

引用回帖:
Originally posted by lwg020810 at 2009-12-3 22:31:

你是不是同时计算频率或其他properties了?如果是的话建议先取消,等优化做最后一步的时候再计算性质

你没有明白我的意思吧,scf是结构优化前面要做的那个自洽,自洽过关了才开始优化结构,然后再计算性质啥的。我肯定没选计算其它性质啊。自洽都过不了,结构优化不能开始,到了scf的步数上限它会给你报错的。
9楼2009-12-03 22:44:40
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

acridine

木虫 (正式写手)

★ ★ ★
jghe(金币+3,VIP+0):thanks 12-4 18:24
增加smearing有可能加速收敛。不过以前我也遇到过计算这种含过渡族金属原子表面的计算,当时只是帮别人计算表面能,而且处理的也不是磁性体系。但是测试了好久终究也没有收敛,最后用castep搞定的。感觉dmol算这种体系确实效率不高。

[ Last edited by acridine on 2009-12-3 at 23:23 ]
10楼2009-12-03 22:58:21
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 jghe 的主题更新
普通表情 高级回复(可上传附件)
信息提示
请填处理意见