| 查看: 2515 | 回复: 10 | |||
[交流]
【求助】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 不知道为何出错，请高手指点！谢谢 |
回复此楼» 猜你喜欢
寻求一种能扛住强氧化性腐蚀性的容器密封件
已经有5人回复
-
真诚求助：手里的省社科项目结项要求主持人一篇中文核心，有什么渠道能发核心吗
已经有7人回复
-
论文投稿，期刊推荐
已经有6人回复
-
请问哪里可以有青B申请的本子可以借鉴一下。
已经有4人回复
-
孩子确诊有中度注意力缺陷
已经有14人回复
-
请问下大家为什么这个铃木偶联几乎不反应呢
已经有5人回复
-
请问有评职称，把科研教学业绩算分排序的高校吗
已经有5人回复
-
2025冷门绝学什么时候出结果
已经有3人回复
-
天津工业大学郑柳春团队欢迎化学化工、高分子化学或有机合成方向的博士生和硕士生加入
已经有4人回复
-
康复大学泰山学者周祺惠团队招收博士研究生
已经有6人回复
» 本主题相关价值贴推荐，对您同样有帮助:
谁能告诉我无水硫酸镁除水大概能到多少0.1%能行吗？
已经有8人回复- TECPLOT10.0
已经有30人回复
- Accurate DOS and Band-structure calculations里的几个脚本
已经有4人回复
- 晶体结构分析入门 Crystal Structure Analysis: A Primer 国际晶体学会的晶体学读本
已经有822人回复
- 求助化学工艺设计软件 RSTT-LC3.0及其注册码
已经有1人回复
- 有关CASTEP Band Structure问题
已经有8人回复
- Journal of Molecular Structure 投稿状态疑问
已经有12人回复
- 【求助】用Materiaol 的 CASTEP计算的band structure中费米能都是0吗?
已经有13人回复
- 重磅推荐经典书Crystal Structure Refinement
已经有18人回复
» 抢金币啦！回帖就可以得到:
-
中国地质大学（武汉）—国家级青年人才杨明教授组-招收博士-新能源材料化学及催化材料
+1/88
-
加拿大/英属哥伦比亚大学曹彦凯课题组招收全奖博士/博后 [机器学习/优化/控制方向]
+1/85
-
硫化物全固态电池的产业化破局：手套箱如何实现全线稳定制造
+1/84
-
华东师范大学 程义云 课题组招2026年博士研究生 - 有机化学、材料化学、高分子合成等
+1/82
-
湘潭大学化学学院陈华杰教授课题组招收有机/高分子方向的博士研究生
+1/81
-
中国科学技术大学 精准智能化学重点实验室 武建昌课题组招聘博士后
+1/79
-
Call for papers，征稿
+1/57
-
国家青年人才叶立群教授课题组招收2026级博士研究生
+1/33
-
中国矿业大学博士招生
+1/33
-
华北电力大学（北京）（第一性原理计算）博士招生——学博，专博各1人
+2/14
-
澳门科技大学2026年数学博士招生——计算物理与数学课题组: 相场与计算流体动力学
+1/14
-
招若干有分子生物，细胞培养，动物实验背景的人员(中山大学)
+1/14
-
南京大学能源与资源学院蔡亮课题组诚招2026年申请-考核制博士生2-3名
+1/10
-
招收26年秋季入学博士生（北科大高精尖学院 力学超材料/机器学习/增材制造相关方向）
+1/10
-
伦敦大学学院（University College London）机械工程系博士招生/CSC联培招生
+1/7
-
招收26年秋季入学博士生（北科大高精尖学院 力学超材料/机器学习/增材制造相关方向）
+1/6
-
北理工柔性电子国家杰青团队招【博士后】【博士】【科研助理】
+1/4
-
2025君科院博士研究生招生欢迎报考
+2/4
-
中科院新疆理化技术研究所人工智能生物医学团队现招聘博士后
+1/3
-
南京航空航天大学核科学与技术方向招收博士生
+1/1
2楼2011-03-23 09:46:13
3楼2011-03-23 22:51:06
4楼2011-03-24 11:17:57
5楼2011-10-11 16:46:33
6楼2011-10-12 08:33:37
★
小木虫(金币+0.5):给个红包，谢谢回帖
小木虫(金币+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 |
7楼2011-11-16 06:34:34
★
小木虫(金币+0.5):给个红包，谢谢回帖
小木虫(金币+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点，望多多帮助哈！ |
8楼2011-11-16 06:38:58
9楼2012-01-31 21:55:33
10楼2015-06-25 21:44:39
11楼2018-09-16 13:13:01