求助：QE计算能带时出错不知如何解决

wtz1013703083

提交bands计算后出现图片中这种报错，求教如何解决

丁越

tianmafei 发表于 2023-8-15 08:02
丁老师好！请您指导，我想用HSE06来计算能带，您可以指导一下简要的做法和步骤吗？先谢谢您！

我这里就以面心立方Si大概简述一下计算流程，不过杂化泛函计算计算量很大，大一点体系有可能算不动，你自己试试看。
1. 做杂化泛函之前先要做一次纯泛函计算，这里就用PBE。
2. HSE06杂化泛函读取PBE收敛的波函数，用startingpot和startingwfc改为file，并且恰当设置HSE06的计算参数

1. &CONTROL
   
2.    calculation     = 'scf'
   
3.    restart_mode    = 'from_scratch'
   
4.    outdir          = './tmp'
   
5.    pseudo_dir      = '/home/dy/qe_pbe_GBRV'
   
6.    prefix          = 'Si'
   
7.    verbosity       = 'low'
   
8. /
   
9. &SYSTEM
   
10.    ibrav           = 0
    
11.    nat             = 8
    
12.    ntyp            = 1
    
13.    ecutwfc         = 44.0
    
14.    ecutrho         = 300.0
    
15.    occupations     = 'smearing'
    
16.    degauss         = 0.005
    
17.    smearing        = 'gaussian'
    
18.    input_dft       = 'hse'
    
19.    exxdiv_treatment = 'gygi-baldereschi'
    
20.    x_gamma_extrapolation = .true.
    
21.    nqx1            = 3
    
22.    nqx2            = 3
    
23.    nqx3            = 3
    
24. /
    
25. &ELECTRONS
    
26.    electron_maxstep = 128
    
27.    conv_thr        = 1.D-6
    
28.    mixing_mode     = 'plain'
    
29.    mixing_beta     = 0.7
    
30.    mixing_ndim     = 8
    
31.    diagonalization = 'david'
    
32.    startingpot     = 'file'
    
33.    startingwfc     = 'file'
    
34. /
    

复制代码

3. 继续用pw.x模块计算能带，直接复制第二步文件，然后把任务类型改为bands，后面 K_POINTS crystal字段中设置k-path，可以将晶体结构文件拖到seek-path网站上去，下面就有适合QE的高对称点路径。
4. 用band.x模块后处理能带，输入文件如下：

1. &BANDS
   
2.    prefix          = 'Si'
   
3.    outdir          = './tmp'
   
4.    filband         = 'bands'
   
5. /
   

复制代码

tianmafei

丁越 发表于 2022-5-18 19:01
参考下面两个文件重新计算。

丁老师好！请您指导，我想用HSE06来计算能带，您可以指导一下简要的做法和步骤吗？先谢谢您！

wtz1013703083

丁越 发表于 2022-5-18 22:56
还有一个问题刚才没注意，为何要tot_charge = -1,这个不是高斯算那种随便就能上个电荷。

谢谢指导，tot_charge=-1在我的计算里是必须的，我现在怀疑是bands.x后处理的问题,因为我用QEtoolkit用bd.dat画出了能带图，但是bd.dat.gnu没有数据

丁越

wtz1013703083 发表于 2022-5-18 18:14
现在是计算出来了。。不过得到的bd.dat.gnu文件是空的orz

还有一个问题刚才没注意，为何要tot_charge = -1,这个不是高斯算那种随便就能上个电荷。

丁越

参考下面两个文件重新计算。

wtz1013703083

现在是计算出来了。。不过得到的bd.dat.gnu文件是空的orz

wtz1013703083

@丁越 这是我的relax和band输入文件

丁越

不贴出输入文件没法凭着这种报错一下就看出问题