应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 第一性原理 » Vasp&MedeA » ENCUT测试
131/2返回列表12下一页
查看: 2328  |  回复: 12
只看楼主@他人 存档 新回复提醒 (忽略) 收藏    在APP中查看

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了，能量还没有收敛。哪里出现了问题吗？
回复此楼

» 猜你喜欢
1楼2013-01-09 11:21:55
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

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的回帖

colinchen2

金虫 (小有名气)

【答案】应助回帖

感谢参与，应助指数 +1
....楼主这得跑多久......

小点的晶胞做就行了吧。
赞一下 回复此楼
5楼2013-01-09 16:00:43
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

xiaoluoj

金虫 (正式写手)
引用回帖:
4楼: Originally posted by emilyoyang at 2013-01-09 15:54:10
计算量太大了啊   一般都是用小的胞测试 然后扩展到supercell上的...

我怕单胞扩展到超胞会出现状况，所以用超胞测试。撇开计算量不谈，因为结果已经算出来了，现在的问题是ENCUT不收敛
赞一下 回复此楼
6楼2013-01-09 16:03:03
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

xiaoluoj

金虫 (正式写手)
引用回帖:
5楼: Originally posted by colinchen2 at 2013-01-09 16:00:43
....楼主这得跑多久......

小点的晶胞做就行了吧。

24核，自洽计算1,2个小时搞定
    现在问题是为什么不收敛，1000应该已经太大了，文献上是400，我想用500（这也挺大的），但是500和后面的值差距有点大。如果用700的话，我看了一下时间大概是前者的2倍有点耗时间
赞一下 回复此楼
7楼2013-01-09 16:40:14
 已阅   回复此楼   关注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的回帖

xiaoluoj

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

running on   11 nodes
distr:  one band on    1 nodes,   11 groups
vasp.5.2.2 15Apr09 complex
POSCAR found type information on POSCAR  N  Si
POSCAR found :  2 types and  280 ions
scaLAPACK will be used
LDA part: xc-table for Pade appr. of Perdew
POSCAR, INCAR and KPOINTS ok, starting setup
FFT: planning ...(           1 )
reading WAVECAR
entering main loop
       N       E                     dE             d eps       ncg     rms          rms(c)
DAV:   1     0.648341706010E+04    0.64834E+04   -0.62749E+05  6292   0.148E+03
DAV:   2    -0.198633704846E+04   -0.84698E+04   -0.82035E+04  7590   0.316E+02
DAV:   3    -0.260881746868E+04   -0.62248E+03   -0.61848E+03  7282   0.898E+01
DAV:   4    -0.262567178229E+04   -0.16854E+02   -0.16800E+02  7876   0.148E+01
DAV:   5    -0.262604758145E+04   -0.37580E+00   -0.37522E+00  7381   0.226E+00    0.103E+02
RMM:   6    -0.229243796515E+04    0.33361E+03   -0.13557E+03  6700   0.332E+01    0.288E+01
RMM:   7    -0.229130484895E+04    0.11331E+01   -0.13174E+02  7342   0.996E+00    0.187E+01
RMM:   8    -0.229118137113E+04    0.12348E+00   -0.16175E+01  7178   0.440E+00    0.578E+00
RMM:   9    -0.229007908022E+04    0.11023E+01   -0.38694E+00  7377   0.183E+00    0.131E+00
RMM:  10    -0.229033387579E+04   -0.25480E+00   -0.28771E+00  7375   0.147E+00    0.247E+00
RMM:  11    -0.229013411261E+04    0.19976E+00   -0.59105E-01  7374   0.741E-01    0.318E-01
RMM:  12    -0.229013715473E+04   -0.30421E-02   -0.43355E-02  7378   0.204E-01    0.451E-02
RMM:  13    -0.229013781928E+04   -0.66455E-03   -0.17970E-03  7418   0.521E-02    0.856E-02
RMM:  14    -0.229013780946E+04    0.98257E-05   -0.24151E-03  7377   0.437E-02    0.276E-02
RMM:  15    -0.229013780307E+04    0.63870E-05   -0.86009E-05  6487   0.860E-03
   1 F= -.22901378E+04 E0= -.22901378E+04  d E =0.000000E+00
-----------------------------------------------
   这是ENCUT=500时，自洽计算的输出，INCAR，KPOINTS（2x2x2）如1楼所示。体系没有磁性，“虽然你自洽计算结束了 但计算的能量可能不对”什么意思，为什么自洽计算结束了，能量会不正确？
赞一下 回复此楼
9楼2013-01-09 20:52:11
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

xiaoluoj

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

PS：但是500和700时的总能量如1楼所示，差了大概580meV，如果选用500的话，如何说服审稿人
赞一下 回复此楼
10楼2013-01-09 20:56:06
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 xiaoluoj 的主题更新
131/2返回列表12下一页
普通表情 龙 兔 虎 猫
信息提示
关闭
请填处理意见
关闭