[QE-users] 回复： "Error in routine davcio (13): error writing file ./CrI3.hub " when running pw2bgw

[半城烟沙]

用BerkeleyGW 安装包里的pw2bgw, 覆盖QE的安装文件，重新编译即可


半城烟沙
1364263610 at qq.com



        



         原始邮件
         
       
发件人：猛猛 <lumengmeng22 at mails.ucas.ac.cn>
发件时间：2025年10月30日 21:56
收件人：users at lists.quantum-espresso.org <users at lists.quantum-espresso.org>
主题：[QE-users] "Error in routine davcio (13): error writing file ./CrI3.hub " when running pw2bgw





Dear all:
I'm using QE-7.5 +BerkeleyGW4.0 to calculate the excitons in Monolayer CrI3, following the calculation in paper :"Physical origin of giant excitonic and magneto optical responses in two-dimensional ferromagnetic insulators" by Meng Wu et. al.  First, I did the pw.x calculation which has nothing wrong, and it generated some wfc files like:
CrI3.hub1, CrI3.hub2, ...
CrI3.hubnoS1, CrI3.hubnoS2,...
CrI3.wfc1, CrI3.wfc2, ...
then I did the pw2bgw.x calculation, the output showed as the title:
Error in routine davcio (13):
     error writing file "/public3/home/sc40681/qe-test/CrI3/wfn/./CrI3.hub20"
 Error in routine davcio (13):
     error writing file "/public3/home/sc40681/qe-test/CrI3/wfn/./CrI3.hub37"
 Error in routine davcio (13):
     error writing file "/public3/home/sc40681/qe-test/CrI3/wfn/./CrI3.hub39"
......
The following are the pw.x input files and the pw2bgw.x input files:




pw.in:
&CONTROL
  prefix='CrI3'
  calculation='bands',
  pseudo_dir='../../../QE/pseudo/ld1',
  outdir='./',
  restart_mode ='from_scratch'
  verbosity='high',
  tprnfor=.true.,
  tstress=.true.,
 forc_conv_thr=1.0d-4,
  nstep=100,
/
&SYSTEM
  ibrav= 0,
  nat= 8,
  ntyp= 2,
  nbnd=130
  occupations = 'smearing',
  smearing = 'gauss',
  nspin=4
  noncolin = .true.
  lspinorb = .true.
  starting_magnetization(1)=0.6
  starting_magnetization(2)=0.2
  degauss = 1.0d-2,
  ecutwfc = 50,
  ecutrho = 500,
/
&ELECTRONS
  conv_thr = 1.0d-8
  mixing_beta = 0.7d0

/
&IONS
/
&CELL
  press_conv_thr=0.1
/
ATOMIC_SPECIES
  Cr 51.9961 Cr_ONCV_PBE_fr.upf
  I 126.90447 I_ONCV_PBE_fr.upf
CELL_PARAMETERS (angstrom)
6.8670000000  0.0000000000  0.0000000000
-3.4335000000  5.9469964478  0.0000000000
 0.0000000000  0.0000000000  18.0000000000
ATOMIC_POSITIONS (crystal)
Cr 0.333333333  0.666666667  0.000377778
Cr 0.666666667  0.333333333  0.999622226
I  0.349896669  0.998803318  0.913106084
I  0.001196661  0.351093352  0.913106084
I  0.648906648  0.650103331  0.913106084
I  0.650103331  0.001196663  0.086893886
I  0.998803318  0.648906648  0.086893886
I  0.351093352  0.349896699  0.086893886

K_POINTS crystal
   36
  0.000000000  0.000000000  0.000000000   1.0
  0.000000000  0.166666667  0.000000000   1.0
  0.000000000  0.333333333  0.000000000   1.0
  0.000000000  0.500000000  0.000000000   1.0
  0.000000000  0.666666667  0.000000000   1.0
  0.000000000  0.833333333  0.000000000   1.0
  0.166666667  0.000000000  0.000000000   1.0
  0.166666667  0.166666667  0.000000000   1.0
  0.166666667  0.333333333  0.000000000   1.0
  0.166666667  0.500000000  0.000000000   1.0
  0.166666667  0.666666667  0.000000000   1.0
  0.166666667  0.833333333  0.000000000   1.0
  0.333333333  0.000000000  0.000000000   1.0
  0.333333333  0.166666667  0.000000000   1.0
  0.333333333  0.333333333  0.000000000   1.0
  0.333333333  0.500000000  0.000000000   1.0
  0.333333333  0.666666667  0.000000000   1.0
  0.333333333  0.833333333  0.000000000   1.0
  0.500000000  0.000000000  0.000000000   1.0
  0.500000000  0.166666667  0.000000000   1.0
  0.500000000  0.333333333  0.000000000   1.0
  0.500000000  0.500000000  0.000000000   1.0
  0.500000000  0.666666667  0.000000000   1.0
  0.500000000  0.833333333  0.000000000   1.0
  0.666666667  0.000000000  0.000000000   1.0
  0.666666667  0.166666667  0.000000000   1.0
  0.666666667  0.333333333  0.000000000   1.0
  0.666666667  0.500000000  0.000000000   1.0
  0.666666667  0.666666667  0.000000000   1.0
  0.666666667  0.833333333  0.000000000   1.0
  0.833333333  0.000000000  0.000000000   1.0
  0.833333333  0.166666667  0.000000000   1.0
  0.833333333  0.333333333  0.000000000   1.0
  0.833333333  0.500000000  0.000000000   1.0
  0.833333333  0.666666667  0.000000000   1.0
  0.833333333  0.833333333  0.000000000   1.0

HUBBARD (ortho-atomic)
U Cr-3d 1.5
J Cr-3d 0.5


pw2bgw.in:
&input_pw2bgw
  prefix = 'CrI3'
  real_or_complex = 2
  wfng_flag = .true.
  wfng_file = 'WFN'
  wfng_kgrid = .true.
  wfng_nk1 = 6
  wfng_nk2 = 6
  wfng_nk3 = 1
  wfng_dk1 = 0
  wfng_dk2 = 0
  wfng_dk3 = 0
/


Best Wishes,
MM
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.quantum-espresso.org/pipermail/users/attachments/20251030/aa4c6774/attachment.html>