24小时热门版块排行榜    

查看: 2318  |  回复: 12
当前只显示满足指定条件的回帖,点击这里查看本话题的所有回帖

xiaoluoj

金虫 (正式写手)

[求助] 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了,能量还没有收敛。哪里出现了问题吗?
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

xiaoluoj

金虫 (正式写手)

例如500,600时总能变化小于0.28eV(见1楼),但是继续增大ENCUT时,总能变化就大于0.28eV了,这怎么理解?
12楼2013-01-09 21:24:37
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 13 个回答

emilyoyang

木虫 (正式写手)

哪有采用supercell测试encut的啊~~~没见过
2楼2013-01-09 13:43:16
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

xiaoluoj

金虫 (正式写手)

引用回帖:
2楼: Originally posted by emilyoyang at 2013-01-09 13:43:16
哪有采用supercell测试encut的啊~~~没见过

为什么不能用超胞?
3楼2013-01-09 14:54:23
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

emilyoyang

木虫 (正式写手)

【答案】应助回帖

感谢参与,应助指数 +1
引用回帖:
3楼: Originally posted by xiaoluoj at 2013-01-09 14:54:23
为什么不能用超胞?...

计算量太大了啊   一般都是用小的胞测试 然后扩展到supercell上的
4楼2013-01-09 15:54:10
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
信息提示
请填处理意见