24小时热门版块排行榜    

Znn3bq.jpeg
查看: 2056  |  回复: 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 的主题更新
普通表情 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 同样的基金本子,换个专家直接从C变A! +3 国自然国社科中 2026-05-19 3/150 2026-05-19 08:50 by Equinoxhua
[考博] 找博士生导师 +5 小代想上岸 2026-05-15 6/300 2026-05-19 08:28 by zhyzzh
[教师之家] 上海大学实验技术岗位非升即走 +10 嘻嘻哈哈乐呵呵 2026-05-15 10/500 2026-05-19 08:05 by 六两废铜
[考博] 博士申请 +5 星…… 2026-05-18 6/300 2026-05-18 23:49 by 糊糊涂涂好
[基金申请] 国自然上会要求 +5 无名者登山 2026-05-18 9/450 2026-05-18 17:50 by BlakeReary
[基金申请] 重磅!青年科学基金项目(C类)资助增幅预计超过50% +7 水和泥不是水泥 2026-05-13 10/500 2026-05-18 07:50 by 水和泥不是水泥
[硕博家园] 我在等一个没有答案的答案 +3 Love_MH 2026-05-17 3/150 2026-05-18 02:22 by 竹林孤影
[文学芳草园] 半夜喝咖啡 +3 myrtle 2026-05-15 5/250 2026-05-18 01:03 by 小沈2018
[找工作] 售SCI一区T0P文章,我:8.O.5.5.1.O.5.4,科目齐全,可+急 +4 ky2p12rrjj 2026-05-15 4/200 2026-05-17 19:47 by Equinoxhua
[考研] 售SCI一区T0P文章,我:8.O.5.5.1.O.5.4,科目齐全,可+急 +4 xx7gd5zq4e 2026-05-15 6/300 2026-05-17 19:36 by Equinoxhua
[公派出国] 售SCI一区T0P文章,我:8.O.5.5.1.O.5.4,科目齐全,可+急 +4 v9tggjlwd0 2026-05-15 5/250 2026-05-17 18:39 by Equinoxhua
[考博] 光量子物理方向 博士招生 1人(2026.09) +3 sandyworld 2026-05-15 4/200 2026-05-17 14:38 by sandyworld
[硕博家园] 售SCI一区T0P文章,我:8.O.5.5.1.O.5.4,科目齐全,可+急 +4 v9tggjlwd0 2026-05-15 4/200 2026-05-17 08:06 by 11n4dfd8yn
[考研] 售SCI一区T0P文章,我:8.O.5.5.1.O.5.4,科目齐全,可+急 +4 l7k6xnh0yc 2026-05-14 8/400 2026-05-17 07:26 by 11n4dfd8yn
[找工作] 售SCI一区T0P文章,我:8.O.5.5.1.O.5.4,科目齐全,可+急 +3 k37jurhrau 2026-05-16 3/150 2026-05-17 01:37 by ue3ir18jc3
[博后之家] 售SCI一区T0P文章,我:8.O.5.5.1.O.5.4,科目齐全,可+急 +4 k37jurhrau 2026-05-16 4/200 2026-05-17 01:35 by ue3ir18jc3
[考研] 售SCI一区T0P文章,我:8.O.5.5.1.O.5.4,科目齐全,可+急 +3 ky2p12rrjj 2026-05-15 4/200 2026-05-17 00:50 by ue3ir18jc3
[公派出国] 售SCI一区T0P文章,我:8.O.5.5.1.O.5.4,科目齐全,可+急 +3 x0mp7owy2b 2026-05-15 4/200 2026-05-17 00:35 by ue3ir18jc3
[有机交流] 如何实现卤原子转化 +3 BT20230424 2026-05-15 5/250 2026-05-16 16:20 by czyzsu
[教师之家] 教学课件你会给同学吗 +8 硕士研究生吗 2026-05-13 8/400 2026-05-14 22:23 by 常规沥青
信息提示
请填处理意见