"Error in routine davcio (13): error writing file ./CrI3.hub " when running pw2bgw

16 views

Skip to first unread message

MM L

unread,

Oct 30, 2025, 9:55:02 AMOct 30

to he...@berkeleygw.org

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
/