| 查看: 3115 | 回复: 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 不知道为何出错,请高手指点!谢谢 |
» 猜你喜欢
FileCode能看出啥?
已经有23人回复
应该是93bebmhtak前后十一个字符比较关键
已经有36人回复
关于Filecode分析方法
已经有9人回复
2019年青年基金涵评意见,大家看看几个A,几个B?
已经有10人回复
不应该看fileCode
已经有6人回复
好奇怪的filecode
已经有6人回复
长年满屏的广告,版主太不责任了。基金也等的急
已经有4人回复
filecode
已经有15人回复
职称评审,求友友推荐见刊最快的期刊
已经有4人回复
奇怪,两个人的filecode固定段从头到尾一模一样
已经有10人回复
» 本主题相关价值贴推荐,对您同样有帮助:
谁能告诉我无水硫酸镁除水大概能到多少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人回复
» 抢金币啦!回帖就可以得到:
今天看到省基金挂了,希望国基面上能中吧
+5/775
北京理工大学集成电路学院-柔性电子器件与智造研究所-2027级博士招生
+1/468
北京信息科技大学机电工程学院硕士研究生招生线上宣讲会-8月16日
+5/460
期待、焦虑与不确定:国自然发榜前心理调查
+5/370
诚征女友,坐标江西
+2/252
武汉女征男
+2/212
北京信息科技大学机电工程学院硕士研究生招生线上宣讲会——8月16日
+1/96
毛细管流变仪,DSC, 力学性能,注塑、挤出、平板硫化机等成型与测试
+1/90
中国科学院空天院精密干涉测量团队 精密机械结构工程师招聘启事
+1/85
招收本科/研究生实习生(2026年9月开始可以长期实习)
+1/81
香港理工大学PhD/RA/PostDoc in OR/AI/OM
+1/79
物理学报录用,散币感谢合作者和审稿人
+2/44
北京理工大学光电学院拟神经计算器件方向,2027级博士/硕士招生
+1/30
哈工大(深圳)何自开团队诚招2027级化学工程与技术专业博士生
+1/26
南京航空航天大学招收27年力学硕士和博士生研究生
+1/11
北京理工大学-集成电路与电子学院杰青团队-招博士后
+1/7
同济大学环境学院肖倩课题组2026年补充招收推免生
+1/6
海南大学化学院—功能分子器件团队2027博士招生+博士后招聘
+1/5
中山大学李国强课题组招收博士生(柔性电子方向)
+1/4
悉尼科技大学(QS 第87)数据科学/AI招收2027年秋季入学全奖及CSC博士生
+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










回复此楼