24小时热门版块排行榜

返回列表

当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖

byin

铜虫 (正式写手)

1ST强帖: 1
应助: 0 (幼儿园)
金币: 245
帖子: 395
在线: 174.4小时
虫号: 948127

[交流] 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
#------------------------
~
~

回复此楼

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

查看全部散金贴

1楼 2011-09-20 11:24:08

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

byin

铜虫 (正式写手)

1ST强帖: 1
应助: 0 (幼儿园)
金币: 245
帖子: 395
在线: 174.4小时
虫号: 948127

引用回帖:

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

感谢回复！

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

我主要是参考了相关文献，一般都是这个力收敛标准，也就借用了。