NSCF/PROJWFC 计算时报错找不到 wfc1 文件

indec

在 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. 
   
9. &SYSTEM
   
10.     a                         =  1.69752e+01
    
11.     b                         =  1.89553e+01
    
12.     c                         =  2.11052e+01
    
13.     cosab                     =  3.68527e-01
    
14.     cosac                     =  2.17276e-02
    
15.     cosbc                     =  3.54418e-01
    
16.     degauss                   =  1.00000e-02
    
17.     ecutrho                   =  480
    
18.     ecutwfc                   =  60
    
19.     ibrav                     = 14
    
20.     nat                       = 380
    
21.     ntyp                      = 7
    
22.     nbnd                      = 1400
    
23.     occupations               = "smearing"
    
24. /
    
25. 
    
26. &ELECTRONS
    
27.     conv_thr         =  1.00000e-06
    
28.     electron_maxstep = 200
    
29.     mixing_beta      =  0.7
    
30. /
    
31. 
    
32. K_POINTS {gamma}
    
33. 
    
34. ATOMIC_SPECIES
    
35. xxxxxx
    
36. 
    
37. ATOMIC_POSITIONS {angstrom}
    
38. xxxxxx

复制代码

NSCF out文件报错如下：

1.      negative rho (up, down):  5.911E-02 0.000E+00
   
2. 
   
3. %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
   
4.      Error in routine diropn (10):
   
5.      error opening ./xxxx.wfc1
   
6. %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
   
7. 
   
8.      stopping ...
   
9. ~                  

复制代码

SOC_LYC

我也遇到了同样的问题（6.3版本），对我的体系有用的解决方法是降低ecutwfc、ecutrho值使输出的wfcx.dat别太大

卡开发发

indec 发表于 2021-3-8 01:19
知道。跑完SCF后跑NSCF也是一样的报错，打不开 wfc1 文件。

按道理SCF做NSCF应该不需要*.wfc1而是*.save，但PROJWFC需要*.wfc1。你可以在SCF过程指定outdir到工作目录试试，如果还不行看看是否是版本问题或者编译问题。

indec

kyuu 发表于 2021-3-7 14:04
哥，你可知你算的是scf

知道。跑完SCF后跑NSCF也是一样的报错，打不开 wfc1 文件。

kyuu

哥，你可知你算的是scf