[QE-users] error in band calculation Error in routine read_wfc
venkatesh chandragiri
venkyphysicsiitm at gmail.com
Sat Sep 8 14:59:40 CEST 2018
Dear espresso users,
I am very new to this code. Recently, I have done scf calculation for
Mn2VGa and successfully completed it. Before doing scf, I have thoroughly
checked the ecut_rho, ecut_wfc, k_mesh for the PAW potentials that have
been used. I have used 9x9x9 k-mesh for 16 atom unit cell which give number
of k-points 125. Everything works fine up to scf. Now, I try to calculate
the band structure for a given k-path, by creating the input files such as
nscf.in (with pw.x), band.in (with pw.x) and bandx.in (with bands.x). For
simplicity, I have used same k_mesh (9x9x9) in the nscf calculation. Now I
have some error in my band.in (with pw.x) calculation.
================= error is =====
The potential is recalculated from file :
/home2/phd/ph13d050/scratch/qe_sim/qe_temp/out/MVG.save/charge-density.dat
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
Error in routine read_wfc (29):
cannot open restart file
/home2/phd/ph13d050/scratch/qe_sim/qe_temp/out/MVG.save/wfcup126 for reading
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
stopping …
================================
My band.in file is
=================
&control
calculation = 'bands'
restart_mode='from_scratch'
prefix='MVG'
tprnfor = .true.
tstress = .true.
verbosity='high'
outdir = '/home2/phd/ph13d050/scratch/qe_sim/qe_temp/out/'
pseudo_dir = '/home2/phd/ph13d050/soft/qe-6.3/pseudo/upf_files'
/
&system
ibrav= 0,
celldm(1) = 10.96954,
!celldm(3) = 0.6842,
nat = 16,
ntyp= 4,
ecutwfc = 80,
ecutrho = 400,
starting_magnetization(2) = 0.7
starting_magnetization(3) = 0.7
starting_magnetization(4) = 0.4
occupations = 'tetrahedra'
degauss = 0.003,
smearing = 'marzari-vanderbilt'
nspin = 2
/
&electrons
mixing_beta = 0.2,
conv_thr = 1.0d-8,
mixing_mode='local-TF'
electron_maxstep=200
diagonalization='cg'
* startingwfc='file'*
/
&ions
ion_dynamics = 'bfgs'
/
&cell
cell_dynamics = 'bfgs'
/
CELL_PARAMETERS {alat}
1.003862527 0.000000000 0.000000000
0.000000000 1.003862363 -0.000000044
0.000000000 -0.000000044 1.003862363
ATOMIC_SPECIES
Ga 69.72300 Ga.pbe-dn-kjpaw_psl.0.2.UPF
Mn1 54.938000 Mn.pbe-spn-kjpaw_psl.0.3.1.UPF
Mn2 54.938000 Mn.pbe-spn-kjpaw_psl.0.3.1.UPF
V 50.94150 V.pbe-spnl-kjpaw_psl.1.0.0.UPF
ATOMIC_POSITIONS (crystal)
Ga 0.500000000 0.000000000 0.500000000
Ga 0.500000000 0.500000000 0.000000000
Ga 0.000000000 0.000000000 0.000000000
Ga 0.000000000 0.500000000 0.500000000
Mn1 0.749988367 0.249989736 0.750010264
Mn1 0.250011633 0.249989736 0.750010264
Mn1 0.749988367 0.750010264 0.249989736
Mn1 0.250011633 0.750010264 0.249989736
Mn2 0.250011523 0.249989721 0.249989721
Mn2 0.749988477 0.249989721 0.249989721
Mn2 0.250011523 0.750010279 0.750010279
Mn2 0.749988477 0.750010279 0.750010279
V 0.000000000 0.000000000 0.500000000
V -0.000000000 0.500000000 0.000000000
V 0.500000000 0.000000000 0.000000000
V 0.500000000 0.500000000 0.500000000
K_POINTS crystal_b
7
0.5000000000 0.2500000000 0.7500000000 100
0.5000000000 0.5000000000 0.5000000000 100
0.0000000000 0.0000000000 0.0000000000 100
0.5000000000 0.0000000000 0.5000000000 100
0.5000000000 0.2500000000 0.7500000000 100
0.3750000000 0.3750000000 0.7500000000 100
0.0000000000 0.0000000000 0.0000000000 100
=======================
Is there any error that is connected
to the input parameter *startingwfc='file'* in band.in ..?
FYI: the nscf.in file almost similar with band.in expect the K_POINTS
{automatic} and calculation = 'nscf'
thanks in advance for your replies.
best wishes,
venkatesh
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.quantum-espresso.org/pipermail/users/attachments/20180908/d7f0dfa0/attachment.html>
More information about the users
mailing list