24小时热门版块排行榜    

查看: 3452  |  回复: 31

byin

铜虫 (正式写手)


[交流] vasp离子弛豫不收敛,求方法

supercell有50个原子,并且约束只能Z向移动,即一共有49个独立的自由度。弛豫结束的判据是原子受力<=0.01 eV/A。

可是弛豫快一周了,受力却一直在0.015附近徘徊,还请各位高手支招!

INCAR:

#------------------------
PREC=Accurate
ISTART=1
ICHARG=1
ISPIN=2
MAGMOM=19*1 6*-1 20*3 5*1
#------------------------
ENCUT=600
NELM=100
NELMIN=6
#------------------------
EDIFFG=-1E-2
NSW=1000
IBRION=2
ISIF=2
POTIM=0.1
#------------------------
ISMEAR=0
SIGMA=0.1
#------------------------
IALGO=48
#------------------------
~
~
回复此楼

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

vasp的收敛

» 猜你喜欢

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

» 抢金币啦!回帖就可以得到:

查看全部散金贴

已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
回帖支持 ( 显示支持度最高的前 50 名 )

bird007

荣誉版主 (职业作家)


★ ★
byin(金币+5):谢谢参与
zzy870720z(金币+1): 谢谢建议 2011-09-20 21:59:00
你可以逐步逼近,先设置精度粗一些
再逐步提高精度
19楼2011-09-20 17:50:10
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

mengnan0501

铜虫 (小有名气)



小木虫(金币+0.5):给个红包,谢谢回帖
更改IBRION=1;ADDGRID=.TRUE.
27楼2011-09-22 20:51:18
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
普通回帖

gougou531

禁虫 (正式写手)


byin(金币+5):谢谢参与
本帖内容被屏蔽

2楼2011-09-20 11:43:37
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

ly2005306

木虫 (著名写手)


★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
bingmou(金币-49): 惩罚! 2011-09-20 12:31:22
祝福楼主,有好运!
14楼2011-09-20 12:01:45
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

ly2005306

木虫 (著名写手)


youzhizhe: 抗议无效,第一性原理版规3.3! 2011-09-20 22:05:37
你这就是自动取款机掉钱,还要判拿到钱的人坐牢。
哎,无语!
15楼2011-09-20 13:01:34
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

ly2005306

木虫 (著名写手)


youzhizhe: 第一性原理版规3.3! 2011-09-20 22:05:42
没天理啊!
16楼2011-09-20 13:02:21
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
★ ★
byin(金币+5):谢谢参与
zzy870720z(金币+1): 谢谢提示 2011-09-20 21:58:45
把IBRION=2改为IBRION=1试一试
17楼2011-09-20 13:06:03
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

byin

铜虫 (正式写手)


引用回帖:
2楼: Originally posted by gougou531 at 2011-09-20 11:43:37:
你是在模拟分子动力学模拟吗?呵呵,弛豫判据一定非得以受力为原则吗?试试以能量为判据说不定会收敛。还有就是以能量为判据和以受力为判据这两者之间的区别联系是什么?哪一个更好,或者是说他们分别适合什么样 ...

感谢回复!

从理论上讲,力收敛和能量收敛是等价的吧,F=-dE/dx。F为0对应能量极值。
数值上的区别我就不清楚了,还请赐教。

我主要是参考了相关文献,一般都是这个力收敛标准,也就借用了。
18楼2011-09-20 14:57:58
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

userhung

禁虫 (文学泰斗)


★ ★
byin(金币+5):谢谢参与
zzy870720z(金币+1): 谢谢建议,希望可以具体些 2011-09-20 21:59:39
换一换判据哦,祝福好运~~~~~~~~
20楼2011-09-20 18:24:17
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

userhung

禁虫 (文学泰斗)



byin(金币+5):谢谢参与
the best wishes to you ~~~
22楼2011-09-20 18:24:55
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

userhung

禁虫 (文学泰斗)



byin(金币+5):谢谢参与
the best wishes to you ~~~
23楼2011-09-20 18:25:09
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

userhung

禁虫 (文学泰斗)



byin(金币+5):谢谢参与
the best wishes to you ~~~
24楼2011-09-20 18:25:26
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

userhung

禁虫 (文学泰斗)



byin(金币+5):谢谢参与
the best wishes to you ~~~
25楼2011-09-20 18:25:43
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

userhung

禁虫 (文学泰斗)


★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
byin(金币+5):谢谢参与
youzhizhe(金币-20): 恶意灌水,扣回所得后惩罚!第一性原理版规3.3 2011-09-20 22:07:14
the best wishes to again ~~~
26楼2011-09-20 18:26:05
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

wujing110

禁虫 (正式写手)


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

28楼2011-09-24 19:10:48
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

byin

铜虫 (正式写手)


引用回帖:
28楼: Originally posted by wujing110 at 2011-09-24 19:10:48:
ENCUT 是不是太高了,可以降低点试试看
我最近在做一个体系也老不收敛,烦死了

截断能高反而会导致不收敛???

一个老师还建议我把Ecut提高……
29楼2011-09-24 19:57:50
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

wujing110

禁虫 (正式写手)


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

30楼2011-09-24 22:13:46
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

chgzhong

木虫 (正式写手)


★ ★
小木虫(金币+0.5):给个红包,谢谢回帖
uuv2010(金币+1): 欢迎参加交流 2011-09-26 21:20:49
我们一般是以能量为判据,决定收敛与否的。
31楼2011-09-26 13:59:15
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
简单回复
ly20053063楼
2011-09-20 11:58   回复  
byin(金币+5):谢谢参与
ly20053064楼
2011-09-20 11:58   回复  
byin(金币+5):谢谢参与
cenwanglai(金币-35): 禁止恶意求助帖里纯表 2011-09-20 21:43:36
ly20053065楼
2011-09-20 11:58   回复  
byin(金币+5):谢谢参与
ly20053066楼
2011-09-20 11:58   回复  
byin(金币+5):谢谢参与
ly20053067楼
2011-09-20 11:58   回复  
byin(金币+5):谢谢参与
ly20053068楼
2011-09-20 11:58   回复  
byin(金币+5):谢谢参与
ly20053069楼
2011-09-20 11:59   回复  
byin(金币+5):谢谢参与
ly200530610楼
2011-09-20 11:59   回复  
byin(金币+5):谢谢参与
ly200530611楼
2011-09-20 11:59   回复  
byin(金币+5):谢谢参与
ly200530612楼
2011-09-20 11:59   回复  
byin(金币+5):谢谢参与
ly200530613楼
2011-09-20 11:59   回复  
bingmou(金币-49): 扣除恶意回帖赚取金币 2011-09-20 12:31:03
userhung21楼
2011-09-20 18:24   回复  
byin(金币+5):谢谢参与
zzy870720z(金币-30): 扣回非有意义回复,谢谢合作 2011-09-20 22:00:48
2012-03-09 21:14   回复  
相关版块跳转 我要订阅楼主 byin 的主题更新
普通表情 高级回复 (可上传附件)
信息提示
请填处理意见