24小时热门版块排行榜    

Znn3bq.jpeg
查看: 3344  |  回复: 10

cenwanglai

荣誉版主 (知名作家)

老和山猥琐派九段

优秀版主

[交流] 【求助结束】结构优化的时候如何提高力的收敛精度? 已有9人参与

INCAR中:
..............................................
IBRION=2
EDIFF=1E-4
EDIFFG=-0.01
POTIM=0.2
LREAL=Auto
..............................................


电子步收敛效果很好,甚至到了E-8,但是力一直不收敛。请问:
1)EDIFFG为正的时候表示离子步总能收敛精度,为负的时候表示受力收敛,二者在数量级上的考虑有什么不同?有的文献中设为正,有的设为负,不知道应不应该二者都考虑。或者先考虑正,在考虑负的。
2)我的结构优化EDIFF能很好很快的满足,但是EDIFFG(负)一直不能达到有求。请问怎么提高力的收敛精度?提高截断能?增加k点密度?增加FFT网格密度?

[ Last edited by cenwanglai on 2010-5-10 at 21:00 ]
回复此楼

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

量子化学基础

» 猜你喜欢

» 本主题相关价值贴推荐,对您同样有帮助:

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

gavinliu7390

木虫 (著名写手)

叶落鹰飞

★ ★
hometownlove(金币+2):谢谢交流 2010-04-06 15:54
cenwanglai(金币+2):非常感谢! 2010-04-06 19:56
1) 应该是随便满足一个条件就会收敛。一般来说应该是把EDIFFG调成EDIFF的十分之一。
2) 你可以选择换换IBRION或者是POTIM.因为不同的算法,对收敛影响的还是很大的。也有可能你的结构离平衡太远了。
真理是一点点接近的!
2楼2010-04-06 11:39:58
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

dazao

金虫 (正式写手)

★ ★
hometownlove(金币+2):谢谢交流 2010-04-06 15:54
cenwanglai(金币+2): 2010-04-06 19:54
cenwanglai(金币+1):非常感谢你的答复,我去试试. 2010-04-06 20:00
可以尝试现设EDIFFG=-0.5 收敛了换成-0.03 最后-0.01
逐级逼近收敛 计算速度也会提高
3楼2010-04-06 15:28:17
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

锐利的碎片

木虫 (正式写手)

star watcher

★ ★
小木虫(金币+0.5):给个红包,谢谢回帖交流
aylayl08(金币+1):谢谢交流 2010-04-09 08:28
ediffg取正值是表示体系自由能最小。这个通常做md时采用。
4楼2010-04-09 00:33:12
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

alick1102

铜虫 (小有名气)

★ ★ ★ ★
小木虫(金币+0.5):给个红包,谢谢回帖交流
gavinliu7390(金币+3):谢谢回答! 2010-04-09 10:45
引用回帖:
Originally posted by gavinliu7390 at 2010-04-06 11:39:58:
1) 应该是随便满足一个条件就会收敛。一般来说应该是把EDIFFG调成EDIFF的十分之一。
2) 你可以选择换换IBRION或者是POTIM.因为不同的算法,对收敛影响的还是很大的。也有可能你的结构离平衡太远了。

一版推荐去负值。因为用能量做判标的话不一定能找到真正的稳定点。
到了博士的时候才知道,科学也不过是为了吃饭的。
5楼2010-04-09 09:57:04
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

lono75

铁杆木虫 (著名写手)

★ ★ ★ ★
小木虫(金币+0.5):给个红包,谢谢回帖交流
gavinliu7390(金币+3):谢谢解答! 2010-04-09 10:45
vi INCAR.coarse
$system =  coarse geom optimisation
......
EDIFF  = 1E-2
EDIFFG = -1E-2
......

vi INCAR.fine
$system =  fine geom optimisation
......
EDIFF  = 1E-6
EDIFFG = -1E-4
......
由天命,成人事!
6楼2010-04-09 10:41:25
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

cenwanglai

荣誉版主 (知名作家)

老和山猥琐派九段

优秀版主

★ ★
zzy870720z(金币+2):多谢建议 2010-05-10 08:46:09
综合大家的意见:力的收敛速度,首先要看初始结构的合理性.尽量建立合理的物理模型.然后是逐步提高收敛精度要求.

PNAS上有篇文章说
total energy was coverged to at least the desired accuracy of ~10meV per cell.这个好象精度不高呢.一个CeO2的cell包括4个CeO2 formula. 对于3*3*2的slab,有18/4=4.5个cell.总能精度就是4.5*10E-3=0.045eV.
vasp默认1E-3(0.001).
文献中一般是EDIFF=1E-4,EDIFFG=-1/-2E-3.这样对总能收敛精度就没有要求了.不知道隐含的总能收敛精度是多少.

PNAS:
Optimization fo ionic conductivity in doped ceria (2006, David A. Anderrson,etc)

[ Last edited by cenwanglai on 2010-5-10 at 08:49 ]
7楼2010-05-10 08:35:26
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

gougou531

禁虫 (正式写手)


小木虫(金币+0.5):给个红包,谢谢回帖
本帖内容被屏蔽

8楼2012-01-07 16:18:17
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

fengangelo

禁虫 (小有名气)

★ ★
小木虫(金币+0.5):给个红包,谢谢回帖
liliangfang(金币+1): 多谢交流 2012-01-27 20:05:29
本帖内容被屏蔽

9楼2012-01-27 19:53:57
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

fyuewen

金虫 (正式写手)


小木虫: 金币+0.5, 给个红包,谢谢回帖
引用回帖:
19548398楼: Originally posted by 锐利的碎片 at 2010-04-09 00:33:12
ediffg取正值是表示体系自由能最小。这个通常做md时采用。

VASP手册中说EDIFFG是不适用MD的

VASP 5.2 第58页 可以检索EDIFFG-tag得到。
Physicsloveukulele!Aloha!
10楼2012-07-23 19:29:26
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 cenwanglai 的主题更新
普通表情 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 材料与化工调剂 10+5 下一站上岸@ 2026-04-10 14/700 2026-04-10 10:45 by l_paradox
[考研] 求调剂 +6 雪逢冬 2026-04-10 6/300 2026-04-10 10:15 by Delta2012
[考研] 085400电子信息类(川大控制工程)求调剂可跨专业 求老师联系 +3 626776879 2026-04-08 3/150 2026-04-09 16:05 by 猪会飞
[考研] 化学308分求调剂 +21 你好明天你好 2026-04-07 23/1150 2026-04-08 22:32 by 凯凯要变帅
[考研] 一志愿211,0703化学305分求调剂 +20 严西西戏 2026-04-06 26/1300 2026-04-08 22:12 by 凯凯要变帅
[考研] 考研求调剂 +4 雯??? 2026-04-08 4/200 2026-04-08 21:44 by 土木硕士招生
[考研] 一志愿郑州大学085600求调剂 +21 吃的不少 2026-04-05 24/1200 2026-04-08 16:47 by sunhuadong
[考研] 求调剂,现在还能填的 +3 上岸小莹加油 2026-04-08 3/150 2026-04-08 14:30 by zhq0425
[考研] 323求调剂 +3 林zlu 2026-04-07 4/200 2026-04-07 23:21 by lbsjt
[考研] 277求调剂 数一104分 +9 瓶子PZ 2026-04-05 14/700 2026-04-07 17:52 by 蓝云思雨
[考研] 080500求调剂 +12 黄宇博 2026-04-06 12/600 2026-04-07 12:41 by upczlm1989
[考研] 工科 22408 267求推荐 +4 wanwan00 2026-04-05 5/250 2026-04-06 22:47 by chenzhimin
[考研] 331求调剂 +8 于征yz 2026-04-05 8/400 2026-04-06 00:54 by fmesaito
[考研] 材料调剂 +9 革微桂 2026-04-04 9/450 2026-04-05 08:27 by 544594351
[考研] 344材料与化工调剂 +9 调剂上岸玘 2026-04-03 9/450 2026-04-04 23:10 by happyddm
[考研] 316求调剂 +9 墨辰_Orion926 2026-04-04 9/450 2026-04-04 21:35 by lbsjt
[考研] 求生物学专业调剂-332分 +5 云朵遛弯指南 2026-04-04 5/250 2026-04-04 10:05 by rzh123456
[考研] 322求调剂 +6 FZAC123 2026-04-03 6/300 2026-04-03 22:23 by 科研小专家
[考研] 285求调剂 +5 AZMK 2026-04-03 8/400 2026-04-03 18:17 by AZMK
[考研] 081200-11408-276学硕求调剂 +5 崔wj 2026-04-03 5/250 2026-04-03 15:06 by arrow8852
信息提示
请填处理意见