24小时热门版块排行榜

返回列表

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

xiaoluoj

金虫 (正式写手)

应助: 12 (小学生)
金币: 692.1
散金: 5
红花: 17
帖子: 487
在线: 504.1小时
虫号: 1526088
注册: 2011-12-06
专业: 通信理论与系统

[求助] ENCUT测试

体系为280个原子，选用PAW_PBE势。2种原子A,B的 ENMIN  ，ENMAX分别为
                  A                B
                           184:245    300:400
为了测试ENCUT，之前我首先在ENCUT=700时测试了KPOINTS，INCAR如下
SYSTERM
ISTART=0
#ICHARG=
LWAVE=.False.
PREC=high
ALGO=Fast
LREAL=Auto
LELF=.TRUE.
#LCHARG=.FALSE.
################electronic relaxation############
NELM=80
ENCUT=700
EDIFF=1E-05
#ISPIN=2
###################ion relaxtion################
#NSW=400
#IBRION=1
#ISIF=2
#EDIFFG=-0.001
ISMEAR=-5
分别测试了1x1x1，2x2x2，3x3x3，4x4x4的K点，KPOINTS格式如下：
Automatic
0
Gamma
x x x
0 0 0
其能量分别为-2290.2356 -2290.7178 -2290.7170  -2290.7176，所以测试ENCUT时我采用了2x2x2的KPOINTS。
  分别测试ENCUT=400 500 600 700 800 900 1000，能量如下：
400 -2289.7561
500 -2290.1378
600 -2290.1329
700 -2290.7178
800 -2290.7593
900 -2290.8174
1000 -2290.8880
我的问题是为什么都达到1000了，能量还没有收敛。哪里出现了问题吗？

回复此楼

» 猜你喜欢

各位大神，目前国内有哪些比较好用的逆合成软件？已经有14人回复
国社科系统bug了，是不是要放榜了？已经有10人回复
现代”学阀”该如何界定已经有9人回复
课题组招2027级博士上海工程技术大学激光智能制造方向已经有4人回复
上海工程技术大学激光智能制造课题组 2027级博士研究生招生已经有3人回复
上海工程技术大学激光智能制造课题组｜2027级博士研究生招生公告已经有11人回复
申博发邮件已经有11人回复
我的奶奶已经有3人回复

1楼 2013-01-09 11:21:55

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

xiaoluoj

金虫 (正式写手)

应助: 12 (小学生)
金币: 692.1
散金: 5
红花: 17
帖子: 487
在线: 504.1小时
虫号: 1526088
注册: 2011-12-06
专业: 通信理论与系统

引用回帖:

8楼: Originally posted by emilyoyang at 2013-01-09 17:30:34
可能是你自洽计算不收敛虽然你自洽计算结束了但计算的能量可能不对
你计算的体系有没有磁性？自洽计算磁矩是否收敛到了优化时的磁矩？（我猜可能没有）所以造成你的能量计算本身是不准确的进而造成不同cut ...

PS：但是500和700时的总能量如1楼所示，差了大概580meV，如果选用500的话，如何说服审稿人