24小时热门版块排行榜    

查看: 2324  |  回复: 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的回帖

emilyoyang

木虫 (正式写手)

【答案】应助回帖

★ ★
xiaoluoj: 金币+2, 有帮助 2013-01-09 22:27:46
引用回帖:
7楼: Originally posted by xiaoluoj at 2013-01-09 16:40:14
24核,自洽计算1,2个小时搞定
    现在问题是为什么不收敛,1000应该已经太大了,文献上是400,我想用500(这也挺大的),但是500和后面的值差距有点大。如果用700的话,我看了一下时间大概是前者的2倍有点耗时间...

可能是你自洽计算不收敛  虽然你自洽计算结束了 但计算的能量可能不对
你计算的体系有没有磁性? 自洽计算磁矩是否收敛到了优化时的磁矩?(我猜可能没有) 所以造成你的能量计算本身是不准确的  进而造成不同cutoff自洽能不收敛  
另外 那么大的体系 文献上的400是合理的 你用500肯定没问题 但计算量可能大了点
8楼2013-01-09 17:30:34
已阅   回复此楼   关注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的回帖
信息提示
请填处理意见