计算化学公社

标题: NSCF/PROJWFC 计算时报错找不到 wfc1 文件 [打印本页]

作者
Author:
indec    时间: 2021-3-7 04:32
标题: NSCF/PROJWFC 计算时报错找不到 wfc1 文件
本帖最后由 indec 于 2021-3-8 01:20 编辑

在 QE 里跑完 SCF 计算(能顺利跑完)后,跑 NSCF 任务时,out文件里报错说打不开 xxxx.wfc1 文件,但是在 xxxx.save 文件夹里有 wfc1.dat ,想知道尝试过把 .dat 文件改名成 .wfc1 也还是报错找不到文件。
scf 文件输入如下:
  1. &CONTROL
  2.     calculation = "scf"
  3.     max_seconds =  8.64000e+06
  4.         restart_mode='from_scratch',
  5.         prefix='xxxx',
  6.         pseudo_dir='/pseudo/sssp_eff',
  7. /

  8. &SYSTEM
  9.     a                         =  1.69752e+01
  10.     b                         =  1.89553e+01
  11.     c                         =  2.11052e+01
  12.     cosab                     =  3.68527e-01
  13.     cosac                     =  2.17276e-02
  14.     cosbc                     =  3.54418e-01
  15.     degauss                   =  1.00000e-02
  16.     ecutrho                   =  480
  17.     ecutwfc                   =  60
  18.     ibrav                     = 14
  19.     nat                       = 380
  20.     ntyp                      = 7
  21.     nbnd                      = 1400
  22.     occupations               = "smearing"
  23. /

  24. &ELECTRONS
  25.     conv_thr         =  1.00000e-06
  26.     electron_maxstep = 200
  27.     mixing_beta      =  0.7
  28. /

  29. K_POINTS {gamma}

  30. ATOMIC_SPECIES
  31. xxxxxx

  32. ATOMIC_POSITIONS {angstrom}
  33. xxxxxx
复制代码


NSCF out文件报错如下:
  1.      negative rho (up, down):  5.911E-02 0.000E+00

  2. %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
  3.      Error in routine diropn (10):
  4.      error opening ./xxxx.wfc1
  5. %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

  6.      stopping ...
  7. ~                  
复制代码

作者
Author:
kyuu    时间: 2021-3-7 14:04
哥,你可知你算的是scf
作者
Author:
indec    时间: 2021-3-8 01:19
kyuu 发表于 2021-3-7 14:04
哥,你可知你算的是scf

知道。跑完SCF后跑NSCF也是一样的报错,打不开 wfc1 文件。
作者
Author:
卡开发发    时间: 2021-3-9 09:10
indec 发表于 2021-3-8 01:19
知道。跑完SCF后跑NSCF也是一样的报错,打不开 wfc1 文件。

按道理SCF做NSCF应该不需要*.wfc1而是*.save,但PROJWFC需要*.wfc1。你可以在SCF过程指定outdir到工作目录试试,如果还不行看看是否是版本问题或者编译问题。
作者
Author:
SOC_LYC    时间: 2021-11-2 10:29
我也遇到了同样的问题(6.3版本),对我的体系有用的解决方法是降低ecutwfc、ecutrho值使输出的wfcx.dat别太大




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