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

24小时热门版块排行榜    

Znn3bq.jpeg
小木虫论坛-学术科研互动平台 » 计算模拟区 » 第一性原理 » Vasp&MedeA » 【求助】GW0-Band structure erro
111/1返回列表
查看: 2900  |  回复: 10
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看

B.C.Wang

金虫 (小有名气)

[交流] 【求助】GW0-Band structure erro

我现在正在用GW计算能带，计算过程中出现了下面的错误：
trying to continue reading WAVECAR, but it might fail
POSCAR, INCAR and KPOINTS ok, starting setup
FFT: planning ...(           1 )
reading WAVECAR
the WAVECAR file was read sucessfully
calculate exact exchange contribution

-----------------------------------------------------------------------------
|                                                                             |
|     EEEEEEE  RRRRRR   RRRRRR   OOOOOOO  RRRRRR      ###     ###     ###     |
|     E        R     R  R     R  O     O  R     R     ###     ###     ###     |
|     E        R     R  R     R  O     O  R     R     ###     ###     ###     |
|     EEEEE    RRRRRR   RRRRRR   O     O  RRRRRR       #       #       #      |
|     E        R   R    R   R    O     O  R   R                               |
|     E        R    R   R    R   O     O  R    R      ###     ###     ###     |
|     EEEEEEE  R     R  R     R  OOOOOOO  R     R     ###     ###     ###     |
|                                                                             |
|      The derivative of the wavefunctions with respect to k (WAVEDER)        |
|      can not be found. You must redo the groundstate calculations           |
|      using LOPTICS=.TRUE. in order to write the WAVEDER file.               |
|                                                                             |
-----------------------------------------------------------------------------

the WAVEDER file was not read
energies w=
   0.0   0.0   0.4   0.0   0.7   0.0   1.1   0.0   1.5   0.0
   1.8   0.0   2.2   0.0   2.5   0.0   2.9   0.0   3.2   0.0
   3.6   0.0   3.9   0.0   4.3   0.0   4.6   0.0   4.9   0.0
   5.3   0.0   5.6   0.0   5.9   0.0   6.3   0.0   6.6   0.0
   6.9   0.0   7.2   0.0   7.5   0.0   7.8   0.0   8.1   0.0
   8.4   0.0   8.7   0.0   9.0   0.0   9.3   0.0   9.6   0.0
   9.9   0.0  10.3   0.0  10.6   0.0  10.9   0.0  11.2   0.0
  11.5   0.0  11.8   0.0  12.1   0.0  12.4   0.0  12.8   0.0
  13.1   0.0  13.4   0.0  13.8   0.0  14.1   0.0  14.5   0.0
  14.8   0.0  15.2   0.0  15.6   0.0  16.0   0.0  16.4   0.0
  16.8   0.0  17.2   0.0  17.7   0.0  18.1   0.0  18.6   0.0
  19.1   0.0  19.7   0.0  20.2   0.0  20.8   0.0  21.5   0.0
  22.1   0.0  22.8   0.0  23.6   0.0  24.4   0.0  25.2   0.0
  26.1   0.0  27.1   0.0  28.2   0.0  29.4   0.0  30.7   0.0
  32.1   0.0  33.7   0.0  35.4   0.0  37.4   0.0  39.6   0.0
  42.1   0.0  45.0   0.0  48.4   0.0  52.3   0.0  57.0   0.0
internal ERROR in RE_READ_KPOINTS: the total number of non zero k-points in the
  full Brillouine zone has changed         399         798
internal ERROR in RE_READ_KPOINTS: the total number of non zero k-points in the
  full Brillouine zone has changed         399         798
internal ERROR in RE_READ_KPOINTS: the total number of non zero k-points in the
  full Brillouine zone has changed         399         798
internal ERROR in RE_READ_KPOINTS: the total number of non zero k-points in the
  full Brillouine zone has changed         399         798
internal ERROR in RE_READ_KPOINTS: the total number of non zero k-points in the
  full Brillouine zone has changed         399         798
internal ERROR in RE_READ_KPOINTS: the total number of non zero k-points in the
  full Brillouine zone has changed         399         798
internal ERROR in RE_READ_KPOINTS: the total number of non zero k-points in the
  full Brillouine zone has changed         399         798
我的INCAR文件为：
Start parameters:
  ISPIN   = 1             non-spin polarized calculation?
  PREC    = Accurate      other possibilities: Normal - Low|Medium|High
  LREAL   = .FALSE.
  ISTART  = 1
  ICHARG  = 11
#  ROPT=1E-3 1E-3 1E-3
Electronic relaxation:
#  ALGO   = Fast        algorithm Normal-Blocked Davidson Very_Fast-RMM-DIIS Fast-Mixture
  EDIFF  = 1E-05       stopping criterion
  EDIFFG = -0.01
  IALGO=48
Ionic relaxation:
  NSW    = 0          number of ionic relaxation steps
  IBRION = -1           update method -1-no 0-MD 1-quasi Newton 2-CG 3-damped MD 5-Hessian
#  POTIM  = 0.1         step size | timestep | scale optimization
#  ISIF   = 3           relax 2-only atoms 3-everything 4-atoms,cellshape 5-cellshapE
  ENCUT  = 600
DOS related options:
  ISMEAR = 0          smearing 0-Gaussian 1-Methfessel-Paxton -5-tetrahedron-Blöchl
  SIGMA  = 0.05         width
  NEDOS  = 1000
  NBANDS = 48
  NPAR   = 16
Analysis:
  LORBIT = 11      orbital analysis False-no 10-DOSCAR & PROOUT 11-lm decomp 12-phases
#  LOPTICS = .TRUE.
Hydrid-Functional
#  LHFCALC   = .TRUE.
#   ALGO     = All
#  TIME      = 0.4
#  HFSCREEN  = 0.2
#  ENCUTFOCK = 0
#  NKRED     = 2
#  PRECFOCK  = A
  NELM      = 1
  ALGO      = GW0
  NOMEGA    = 80
  ENCUTGW   = 150
不知道为何出错，请高手指点！谢谢
回复此楼

» 猜你喜欢

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

» 抢金币啦！回帖就可以得到:

查看全部散金贴
1楼 2011-03-23 00:46:02
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

子虚乌有5388

木虫 (正式写手)
★ ★
B.C.Wang(金币+1):谢谢参与
bingmou(金币+1): 感谢分享 2011-03-23 16:44:05
B.C.Wang(金币+1): 2011-03-23 22:51:23
LOPTICS=.TRUE.才行，你的注释掉了。
赞一下(2人) 回复此楼
2楼2011-03-23 09:46:13
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

B.C.Wang

金虫 (小有名气)
引用回帖:
Originally posted by 子虚乌有5388 at 2011-03-23 09:46:13:
LOPTICS=.TRUE.才行，你的注释掉了。

我试过LOPTICS=.TRUE.，也行不通，出现同样的错误！你觉得还有哪些原因导致计算出错？
赞一下 回复此楼
3楼2011-03-23 22:51:06
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

luckstar05

木虫 (小有名气)
★ ★
B.C.Wang(金币+1):谢谢参与
gzqdyouxia(金币+1): 谢谢提示 2011-03-25 10:53:26
Kpoint  在计算过程中都要一致！
赞一下(2人) 回复此楼
4楼2011-03-24 11:17:57
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

liddrr

铜虫 (小有名气)
★
B.C.Wang(金币+1):谢谢参与
你在计算时，KPOINT文件里的是怎么写的？计算后能画出能带结构图吗？
赞一下(1人) 回复此楼
5楼2011-10-11 16:46:33
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

liliangfang

荣誉版主 (著名写手)
★
B.C.Wang(金币+1):谢谢参与
k点设置有问题
回复此楼
6楼2011-10-12 08:33:37
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

may-sucess

金虫 (正式写手)
★
小木虫(金币+0.5):给个红包，谢谢回帖
大家好！我也是出现同样的问题，请 多多帮助！我的INCAR是：
SYSTEM=WZ-AlN(1x1x1)
   3 #Precsion of calculations
  4 #ICHARG=11
  5 ENCUT =400
  6 #PREC = Accurate
   NBANDS=64
  9 #type of run
10 #ISIF=3
11 ISMEAR = 0
12 SIGMA = 0.05
13 #ISPIN=1
14 GGA = PE
15 IBRION=-1
16 NSW=0
17 #EDIFFG=-0.02
18 EDIFF=1E-8
19 NELM=1
20 #NSIM=4
21 #LMAXMIX=4
23 #HSE06
24 ALGO=GW0
25 LHFCALC = .TRUE.
26 HFSCREEN = 0.2
27 TIME = 0.4
28 #ENCUTFOCK = 0
29 #NKRED = 2 #########################
30 AEXX = 0.25
31 ALDAC =0.0
32 AGGAc =1.0
33 PRECFOCK = M
35 #LDOS
36 LORBIT=10
37 NEDOS=2000
38 EMIN=-30.0
39 EMAX=30.0
40 #OPTIC
41 LOPTICS = .TRUE.
42 NOMEGA = 50
出现问题是：
random initialization beyond band           65
38  the WAVECAR file was read sucessfully
39  initial charge from wavefunction
40  calculate exact exchange contribution
41
42  -----------------------------------------------------------------------------
43 |                                                                             |
44 |     EEEEEEE  RRRRRR   RRRRRR   OOOOOOO  RRRRRR      ###     ###     ###     |
45 |     E        R     R  R     R  O     O  R     R     ###     ###     ###     |
46 |     E        R     R  R     R  O     O  R     R     ###     ###     ###     |
47 |     EEEEE    RRRRRR   RRRRRR   O     O  RRRRRR       #       #       #      |
48 |     E        R   R    R   R    O     O  R   R                               |
49 |     E        R    R   R    R   O     O  R    R      ###     ###     ###     |
50 |     EEEEEEE  R     R  R     R  OOOOOOO  R     R     ###     ###     ###     |
51 |                                                                             |
52 |      The derivative of the wavefunctions with respect to k (WAVEDER)        |
53 |      can not be found. You must redo the groundstate calculations           |
54 |      using LOPTICS=.TRUE. in order to write the WAVEDER file.               |
55 |                                                                             |
56  -----------------------------------------------------------------------------
57
58  the WAVEDER file was not read
59 energies w=
60    0.0   0.0   0.6   0.0   1.1   0.0   1.7   0.0   2.3   0.0
61    2.8   0.0   3.4   0.0   3.9   0.0   4.4   0.0   4.9   0.0
62    5.4   0.0   5.9   0.0   6.4   0.0   6.9   0.0   7.4   0.0
63    7.9   0.0   8.3   0.0   8.8   0.0   9.3   0.0   9.8   0.0
64   10.3   0.0  10.8   0.0  11.3   0.0  11.8   0.0  12.3   0.0
65   12.9   0.0  13.5   0.0  14.1   0.0  14.8   0.0  15.5   0.0
。。。。。  
和
internal ERROR in RE_READ_KPOINTS: the total number of non zero k-points in the
  full Brillouine zone has changed          76         384
赞一下(1人) 回复此楼
7楼2011-11-16 06:34:34
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

may-sucess

金虫 (正式写手)
★
小木虫(金币+0.5):给个红包，谢谢回帖
我的KPOINTS文件是：计算OPTICS步骤的K点加上高对称点，如下：
0.50000000000000    0.00000000000000    0.50000000000000             3
39 0.12500000000000    0.12500000000000    0.50000000000000             6
40 0.25000000000000    0.12500000000000    0.50000000000000            12
41 0.37500000000000    0.12500000000000    0.50000000000000            12
42 0.25000000000000    0.25000000000000    0.50000000000000             6
43 0.37500000000000    0.25000000000000    0.50000000000000             6
44  0.0000000E+00  0.0000000E+00  0.5555556E-01  0
45  0.0000000E+00  0.0000000E+00  0.1111111E+00  0
46  0.0000000E+00  0.0000000E+00  0.2222222E+00  0
47  0.0000000E+00  0.0000000E+00  0.2777778E+00  0
48  0.0000000E+00  0.0000000E+00  0.3888889E+00  0
49  0.0000000E+00  0.0000000E+00  0.4444444E+00  0
50 -0.3703333E-01  0.7407444E-01  0.5000000E+00  0
51 -0.7406667E-01  0.1481489E+00  0.5000000E+00  0
52 -0.1111000E+00  0.2222233E+00  0.5000000E+00  0
53 -0.1481333E+00  0.2962978E+00  0.5000000E+00  0
54 -0.1851667E+00  0.3703722E+00  0.5000000E+00  0
。。。。。。。。。。。。。。。。
其中没有重复的K点，望多多帮助哈！
赞一下(1人) 回复此楼
8楼2011-11-16 06:38:58
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

gloomy2004

金虫 (著名写手)
★
小木虫(金币+0.5):给个红包，谢谢回帖
请问Lz, 你的是什么体系？最后怎么解决的？我和你的报同样的错误。
赞一下(1人) 回复此楼
9楼2012-01-31 21:55:33
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

小花虫

新虫 (初入文坛)
★
小木虫: 金币+0.5, 给个红包，谢谢回帖
L.Z请问你的问题解决了吗，我也遇到了同样的问题
赞一下 回复此楼
10楼2015-06-25 21:44:39
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

单调主义者

新虫 (初入文坛)
★
小木虫: 金币+0.5, 给个红包，谢谢回帖
报错原因就是文件夹当中没有WAVEDER文件吧。
将这个文件从上一部计算的文件夹当中拷贝过来。
https://jingyan.baidu.com/article/363872ecce02276e4ba16fe8.html
我是按照这个教程计算的，从第二部计算当中将WAVEDER文件拷贝过来
赞一下 回复此楼
11楼2018-09-16 13:13:01
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 B.C.Wang 的主题更新
111/1返回列表
普通表情 龙 兔 虎 猫 高级回复 (可上传附件)
信息提示
关闭
请填处理意见
关闭