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
_______________________________________________________________________________
The Quantum ESPRESSO Foundation stands in solidarity with all civilians
worldwide who are victims of terrorism, military aggression, and indiscriminate
warfare.
--------------------------------------------------------------------------------
Quantum ESPRESSO is supported by MaX (www.max-centre.eu)
users mailing list [email protected]
https://lists.quantum-espresso.org/mailman/listinfo/users
