[QE-users] error in band calculation Error in routine read_wfc
Pietro Delugas
pdelugas at sissa.it
Sat Sep 8 15:20:46 CEST 2018
hello
in the bands.in you are specifying startingwfc = 'file', which means
taht you are telling the program to get the starting wave functions from
a file generate in a previous calculation. But your previous calculation
was, if I have understood correctly, the nscf.in where the kpoint were
different in number and coordinates. The program is not able to read a
wavefunction file for each of the kpoint and your current calculation
crashes.
If you delete the startingwfc ='file' line from the input the program
will most likely work.
regards - pietro
On 08/09/2018 14:59, venkatesh chandragiri wrote:
> 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 <http://nscf.in> (with
> pw.x), band.in <http://band.in> (with pw.x) and bandx.in
> <http://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 <http://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 <http://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 <http://band.in> ..?
>
> FYI: the nscf.in <http://nscf.in> file almost similar with band.in
> <http://band.in> expect the K_POINTS {automatic} and calculation = 'nscf'
>
> thanks in advance for your replies.
>
> best wishes,
> venkatesh
>
>
> _______________________________________________
> users mailing list
> users at lists.quantum-espresso.org
> https://lists.quantum-espresso.org/mailman/listinfo/users
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.quantum-espresso.org/pipermail/users/attachments/20180908/53f4a05b/attachment.html>
More information about the users
mailing list