计算化学公社

标题: 求助：QE计算能带时出错不知如何解决 [打印本页]
作者
Author: wtz1013703083    时间: 2022-5-18 14:33
标题: 求助：QE计算能带时出错不知如何解决
提交bands计算后出现图片中这种报错，求教如何解决
作者
Author: 丁越    时间: 2022-5-18 15:45
不贴出输入文件没法凭着这种报错一下就看出问题
作者
Author: wtz1013703083    时间: 2022-5-18 16:25
@丁越 这是我的relax和band输入文件
作者
Author: wtz1013703083    时间: 2022-5-18 18:14
现在是计算出来了。。不过得到的bd.dat.gnu文件是空的orz
作者
Author: 丁越    时间: 2022-5-18 19:01
本帖最后由 丁越 于 2023-8-15 22:17 编辑

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


作者
Author: 丁越    时间: 2022-5-18 22:56
wtz1013703083 发表于 2022-5-18 18:14
现在是计算出来了。。不过得到的bd.dat.gnu文件是空的orz

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

谢谢指导，tot_charge=-1在我的计算里是必须的，我现在怀疑是bands.x后处理的问题,因为我用QEtoolkit用bd.dat画出了能带图，但是bd.dat.gnu没有数据
作者
Author: tianmafei    时间: 2023-8-15 08:02
丁越 发表于 2022-5-18 19:01
参考下面两个文件重新计算。

丁老师好！请您指导，我想用HSE06来计算能带，您可以指导一下简要的做法和步骤吗？先谢谢您！
作者
Author: 丁越    时间: 2023-8-17 21:08
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. /
复制代码







欢迎光临 计算化学公社 (http://bbs.keinsci.com/) Powered by Discuz! X3.3