24小时热门版块排行榜    

查看: 2102  |  回复: 11

490105179

木虫 (正式写手)

[求助] 石墨烯能带问题求助 已有1人参与

我计算3×3石墨烯超胞的能带,但是画出来的能带图很诡异,如下图。
INCAR:
SYSTEM = C               
ENCUT = 500              
                        
PREC     = Accurate      
                        
Ionic Relaxation 1      
IBRION   =  1            
NSW      =200            
EDIFFG   =-0.02         
ISIF = 2                 
POTIM=0.1               
                        
Electronici Relaxation 1
ISMEAR = 0               
SIGMA=0.2               
NELM= 100               
NELMIN= 4               
ALGO=Fast               
EDIFF= 1E-5              
                        
LREAL= Auto              
GGA     = PE            
ISTART=1                 
ICHARG = 11              
LORBIT=11               
NPAR=4                  
LUSE_VDW=.TRUE.         
Zab_vdW=-1.8867         
AGGAC=0.0000            
K点
k-points along high symmetry lines  
40                                 
Line-mode                           
  rec                              
0.000  0.000  0.000 ! G            
0.000  0.500  0.000 ! M            
                                    
0.000  0.500  0.000 ! M            
0.333  0.667  0.000 ! K            
                                    
0.333  0.667  0.000 ! K            
0.000  0.000  0.000 ! G            
                                    
哪位大神帮帮忙,急求。金币不是问题,不够可以追加
回复此楼

» 猜你喜欢

» 本主题相关商家推荐: (我也要在这里推广)

» 本主题相关价值贴推荐,对您同样有帮助:

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

490105179

木虫 (正式写手)

为什么不能上传图片了?
2楼2016-08-08 16:13:17
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

490105179

木虫 (正式写手)

http://pan.baidu.com/s/1o7WRTrW 无法直接上传图片,只能放在网盘里了
3楼2016-08-08 16:16:06
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

c19940408

银虫 (初入文坛)

试试把能带计算的标准输出帖出来吧。。以及可以把能带步的INCAR单独贴出来吗?1L的有点混乱
4楼2016-08-08 16:46:00
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

490105179

木虫 (正式写手)

引用回帖:
4楼: Originally posted by c19940408 at 2016-08-08 16:46:00
试试把能带计算的标准输出帖出来吧。。以及可以把能带步的INCAR单独贴出来吗?1L的有点混乱

你好,我的账号今天不能上传图片跟文件。1L的贴出的就是能带步的INCAR,不知道您说的能带计算的标准是什么意思?
5楼2016-08-08 16:52:02
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

490105179

木虫 (正式写手)

引用回帖:
4楼: Originally posted by c19940408 at 2016-08-08 16:46:00
试试把能带计算的标准输出帖出来吧。。以及可以把能带步的INCAR单独贴出来吗?1L的有点混乱

或者您能加我QQ聊吗?小木虫账号就是我QQ号
6楼2016-08-08 16:52:39
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

c19940408

银虫 (初入文坛)

★ ★
漫天飘雪: 金币+2, 谢谢交流 2016-08-08 18:39:03
引用回帖:
5楼: Originally posted by 490105179 at 2016-08-08 16:52:02
你好,我的账号今天不能上传图片跟文件。1L的贴出的就是能带步的INCAR,不知道您说的能带计算的标准是什么意思?...

为啥你的能带计算是和几何结构优化合在一起的。。。不是应该先优化几何结构,再做能带计算吗。。标准输出就是vasp的输出结果。
7楼2016-08-08 16:57:26
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

490105179

木虫 (正式写手)

引用回帖:
7楼: Originally posted by c19940408 at 2016-08-08 16:57:26
为啥你的能带计算是和几何结构优化合在一起的。。。不是应该先优化几何结构,再做能带计算吗。。标准输出就是vasp的输出结果。...

我是在静态计算基础上计算能带的,您的意思是不是说在计算能带时要把INCAR中结构优化的参数去掉?下面是asp输出结果
running on   48 total cores
distrk:  each k-point on   48 cores,    1 groups
distr:  one band on   12 cores,    4 groups
using from now: INCAR     
vasp.5.3.5 31Mar14 (build Sep 28 2014 11:30:43) complex                        
  
POSCAR found type information on POSCAR  C
POSCAR found :  1 types and      18 ions
scaLAPACK will be used

-----------------------------------------------------------------------------
|                                                                             |
|  ADVICE TO THIS USER RUNNING 'VASP/VAMP'   (HEAR YOUR MASTER'S VOICE ...):  |
|                                                                             |
|      You have switched on vdW-DFT.                                          |
|      This routine was written and supplied by Jiri Klimes.                  |
|      We recommed that you carefully read and cite the following             |
|      publication                                                            |
|      J. Klimes, D.R. Bowler, A. Michelides                                  |
|        J. Phys.: Cond Matt. 22 022201 (2010)                                |
|      J. Klimes, D.R. Bowler, A. Michelides                                  |
|        Phys. Rev. B. 83, 195131 (2011)                                      |
|      and references therein.                                                |
|                                                                             |
-----------------------------------------------------------------------------

LDA part: xc-table for Pade appr. of Perdew
found WAVECAR, reading the header
  number of k-points has changed, file:     6 present:   120
  trying to continue reading WAVECAR, but it might fail
WARNING: stress and forces are not correct
POSCAR, INCAR and KPOINTS ok, starting setup
WARNING: small aliasing (wrap around) errors must be expected
FFT: planning ...
reading WAVECAR
the WAVECAR file was read successfully
charge-density read from file: C                                    
entering main loop
       N       E                     dE             d eps       ncg     rms          rms(c)
Total vdW correction in eV:     11.1729193
DAV:   1    -0.165193018116E+03   -0.16519E+03   -0.37785E+02 13096   0.538E+01
RMM:   2    -0.168452125373E+03   -0.32591E+01   -0.32591E+01 12719   0.147E+01
RMM:   3    -0.168358523362E+03    0.93602E-01   -0.10100E+00 14772   0.302E+00
RMM:   4    -0.168366502218E+03   -0.79789E-02   -0.78854E-02 15034   0.740E-01
RMM:   5    -0.168366950658E+03   -0.44844E-03   -0.44416E-03 13419   0.155E-01
RMM:   6    -0.168366987254E+03   -0.36597E-04   -0.36009E-04 10346   0.401E-02
RMM:   7    -0.168366989162E+03   -0.19074E-05   -0.18382E-05  7823   0.105E-02
Total vdW correction in eV:     11.1729193
   1 F= -.16836699E+03 E0= -.16835751E+03  d E =-.189501E-01
writing wavefunctions
8楼2016-08-08 17:02:02
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

c19940408

银虫 (初入文坛)

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
感谢参与,应助指数 +1
490105179: 金币+50, ★★★★★最佳答案 2016-08-08 17:17:04
引用回帖:
8楼: Originally posted by 490105179 at 2016-08-08 17:02:02
我是在静态计算基础上计算能带的,您的意思是不是说在计算能带时要把INCAR中结构优化的参数去掉?下面是asp输出结果
running on   48 total cores
distrk:  each k-point on   48 cores,    1 groups
distr: ...

可以试一下把几何优化参数去掉,另外设置istart=0,以及把lreal去掉,重新拷贝静态计算的电荷,并且增加k点。看看结果如何

[ 发自手机版 http://muchong.com/3g ]
9楼2016-08-08 17:13:46
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

490105179

木虫 (正式写手)

引用回帖:
9楼: Originally posted by c19940408 at 2016-08-08 17:13:46
可以试一下把几何优化参数去掉,另外设置istart=0,以及把lreal去掉,重新拷贝静态计算的电荷,并且增加k点。看看结果如何
...

好的,我试试,谢谢您
10楼2016-08-08 17:16:43
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 490105179 的主题更新
信息提示
请填处理意见