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

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
不知道为何出错,请高手指点!谢谢
回复此楼

» 猜你喜欢

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

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

查看全部散金贴

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

小花虫

新虫 (初入文坛)


★
小木虫: 金币+0.5, 给个红包,谢谢回帖
L.Z请问你的问题解决了吗,我也遇到了同样的问题
10楼2015-06-25 21:44:39
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 11 个回答

子虚乌有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楼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  在计算过程中都要一致!
4楼2011-03-24 11:17:57
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
普通表情 龙 兔 虎 猫 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[有机交流] OpenAI又攻克数学难题了,离AI做有机合成还有多远? +4 wangyl_123 2026-10-10 8/400 2026-10-12 18:09 by wangyl_123
[硕博家园] 极限拉扯,如何抉择 +5 Isak0320 2026-10-12 7/350 2026-10-12 18:03 by 有人要害朕
[教师之家] 课题组会帮你通过考核么? +3 好运来~9527 2026-10-10 6/300 2026-10-12 14:55 by 简单化xn
[公派出国] 售SCI一区文章,我:8O5.5.1.O5.4,科目全,可伽急 +5 qR2AdCWTVtkw 2026-10-09 10/500 2026-10-11 22:26 by XQNyDCa1Tk5n
[找工作] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +5 qR2AdCWTVtkw 2026-10-09 8/400 2026-10-11 22:26 by XQNyDCa1Tk5n
[考博] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +5 qR2AdCWTVtkw 2026-10-09 9/450 2026-10-11 22:15 by XQNyDCa1Tk5n
[考博] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +3 qR2AdCWTVtkw 2026-10-09 7/350 2026-10-11 22:11 by XQNyDCa1Tk5n
[教师之家] 售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急 (EPI+-1)(金币-50) +5 qR2AdCWTVtkw 2026-10-09 10/500 2026-10-11 22:06 by XQNyDCa1Tk5n
[基金申请] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +5 qR2AdCWTVtkw 2026-10-09 9/450 2026-10-11 22:06 by XQNyDCa1Tk5n
[教师之家] 售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急 (EPI+-1)(金币-50) +5 qR2AdCWTVtkw 2026-10-09 10/500 2026-10-11 14:25 by VOMexwnDzOFZ
[教师之家] 售SCI文章,我:8O.5.5.1O.54,科目全,可十急 (EPI+-1)(金币-50) +4 qR2AdCWTVtkw 2026-10-09 7/350 2026-10-11 14:05 by VOMexwnDzOFZ
[考研] 售SCI文章,我:8O.5.5.1O.54,科目全,可十急 +5 qR2AdCWTVtkw 2026-10-09 8/400 2026-10-11 13:51 by VOMexwnDzOFZ
[硕博家园] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +5 qR2AdCWTVtkw 2026-10-09 8/400 2026-10-11 13:45 by VOMexwnDzOFZ
[找工作] 售SCI一区T0P文章,我:8O.55.1.O.5.4,科目齐全,可+急 +3 qR2AdCWTVtkw 2026-10-09 7/350 2026-10-11 10:54 by VOMexwnDzOFZ
[找工作] 售SCI一区T0P文章,我:8O.55.1.O.5.4,科目齐全,可+急 +4 qR2AdCWTVtkw 2026-10-09 7/350 2026-10-11 10:34 by VOMexwnDzOFZ
[考博] 招收博士生一名 +3 mandolin 2026-10-10 3/150 2026-10-11 10:21 by 不是山谷!
[公派出国] 售SCI一区T0P文章,我:8O.55.1.O.54,科目全,可伽急 +4 qR2AdCWTVtkw 2026-10-09 4/200 2026-10-11 03:14 by M0zrYVSAnxpE
[有机交流] 各位大佬,路线设计求助 +6 zapen 2026-10-08 8/400 2026-10-10 15:59 by zapen
[基金申请] 入职后第一次管课题组的专利缴费,对账才发现官费比别人多掏了一倍 +5 13108017953 2026-10-08 5/250 2026-10-10 11:15 by wdwd123
[有机交流] 我应该选哪条合成路线让实验成功率高一些? +5 jing816y 2026-10-08 9/450 2026-10-10 09:05 by zapen
信息提示
请填处理意见