24小时热门版块排行榜    

查看: 1931  |  回复: 9

meizhixin

新虫 (小有名气)

[交流] vasp优化非晶体系一直不收敛怎么办?

我用vasp优化100个原子的非晶合金体系,精度设置为EDIFF=1E-5,EDIFFG=-0.01,一直达不到精度该怎么办?而且后来我发现即使EDIFF和EDIFFG都不设置,低精度下计算60步也不收敛,,,,,肿么办??

发自小木虫Android客户端
回复此楼

» 猜你喜欢

已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

ra2ghgzh

木虫 (正式写手)

学痴

★ ★ ★
小木虫: 金币+0.5, 给个红包,谢谢回帖
meizhixin(漫天飘雪代发): 金币+2, 谢谢交流 2017-07-07 10:02:57
好像原子数多了就是会这样,我正在算一个200多原子的任务,有时电子步要走900多才收敛,离子步走100还不收敛。但是只要时间足够最终还是会收敛的。
我永远当不了老板
2楼2017-07-07 09:28:44
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

meizhixin

新虫 (小有名气)

引用回帖:
2楼: Originally posted by ra2ghgzh at 2017-07-07 09:28:44
好像原子数多了就是会这样,我正在算一个200多原子的任务,有时电子步要走900多才收敛,离子步走100还不收敛。但是只要时间足够最终还是会收敛的。

我之前帮别人算一个128个原子的掺杂体系,是晶态的,先低精度下优化,计算14步就达到精度了,提高精度设置EDIFF=1E-5, EDIFFG = -0.01计算8步就收敛了,,,可能也与我的结构有关系吧

发自小木虫Android客户端
3楼2017-07-07 10:04:00
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

meizhixin

新虫 (小有名气)

引用回帖:
2楼: Originally posted by ra2ghgzh at 2017-07-07 09:28:44
好像原子数多了就是会这样,我正在算一个200多原子的任务,有时电子步要走900多才收敛,离子步走100还不收敛。但是只要时间足够最终还是会收敛的。

那你200个原子体系达到精度,计算了多长时间???

发自小木虫Android客户端
4楼2017-07-07 10:07:05
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

ra2ghgzh

木虫 (正式写手)

学痴


小木虫: 金币+0.5, 给个红包,谢谢回帖
引用回帖:
4楼: Originally posted by meizhixin at 2017-07-07 10:07:05
那你200个原子体系达到精度,计算了多长时间???
...

快的用了两天,慢的已经有两周了
我永远当不了老板
5楼2017-07-07 10:49:01
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

西风月123

木虫 (著名写手)


小木虫: 金币+0.5, 给个红包,谢谢回帖
优化不收敛的话,除了降低门限,如EDIFF,EDIFFG之外,还可以采用换IBRION优化方式,以及换算法ALGO等等,这些都有利于收敛。

此外,如果体系很大,计算资源也有限的话,推荐精度设置可以设置为EDIFF=1E-4,EDIFFG=-0.05

还有就是在60步内收敛属于正常现象,60步内不收敛也没必要难过,不管多少步关键还是要看能不能达到收敛标,以及优化出来的结构理不理想。

我们这边也做非晶,你是做非晶是体材料,还是界面呢?
6楼2017-07-07 15:55:39
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

meizhixin

新虫 (小有名气)

引用回帖:
6楼: Originally posted by 西风月123 at 2017-07-07 15:55:39
优化不收敛的话,除了降低门限,如EDIFF,EDIFFG之外,还可以采用换IBRION优化方式,以及换算法ALGO等等,这些都有利于收敛。

此外,如果体系很大,计算资源也有限的话,推荐精度设置可以设置为EDIFF=1E-4,EDIFFG ...

我做的是体材料,我是想用能量应变法计算下弹性常数,但施加应变后的结构进行优化很难收敛,估计结构很不理想,原子间距离过小,你说的那些我也考虑了,还是不行,你有计算过非晶的弹性常数吗?

发自小木虫Android客户端
7楼2017-07-07 22:29:08
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

西风月123

木虫 (著名写手)


小木虫: 金币+0.5, 给个红包,谢谢回帖
引用回帖:
7楼: Originally posted by meizhixin at 2017-07-07 22:29:08
我做的是体材料,我是想用能量应变法计算下弹性常数,但施加应变后的结构进行优化很难收敛,估计结构很不理想,原子间距离过小,你说的那些我也考虑了,还是不行,你有计算过非晶的弹性常数吗?
...

我做的是非晶二氧化硅,内容有电子顺磁共振参数(ESR),以及做CI-NEB模拟等。非晶的弹性常数没算过。
但是,好像记得在施加应用时,要渐进地增加,这样有利于收敛~~~
8楼2017-07-08 11:19:13
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

meizhixin

新虫 (小有名气)

引用回帖:
8楼: Originally posted by 西风月123 at 2017-07-08 11:19:13
我做的是非晶二氧化硅,内容有电子顺磁共振参数(ESR),以及做CI-NEB模拟等。非晶的弹性常数没算过。
但是,好像记得在施加应用时,要渐进地增加,这样有利于收敛~~~...

谢谢你的回答,我觉得我这可能就是加应变后结构不理想,所以需要很长时间才能收敛,我用的是候博士的脚本来获得应变后的结构的,打算再改下参数试试,看多久能收敛

发自小木虫Android客户端
9楼2017-07-08 18:26:30
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
10楼2017-07-08 19:05:26
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 meizhixin 的主题更新
普通表情 高级回复 (可上传附件)
信息提示
请填处理意见