[QE-users] error in band calculation Error in routine read_wfc
venkatesh chandragiri
venkyphysicsiitm at gmail.com
Sun Sep 9 08:27:56 CEST 2018
Dear Dr. Pietro Delugas
<https://www.mail-archive.com/search?l=users@lists.quantum-espresso.org&q=from:%22Pietro+Delugas%22>
thanks for you reply. Now pw.x bands.in works fine without error. However,
I have not got any data in plotbands.xmgr.* and a pdf file generated from
ps2pdf <plotband.ps> , do not contain any data to be plot. please check my
bandx.in (run with bandx.x) and plotbands.in (run with plotbands.x) as
given below.
====================bandx.in===========
&BANDS
outdir='/home2/phd/ph13d050/scratch/qe_sim/qe_temp/out/'
prefix='MVG'
filband="bandx.dat"
/
===============
=====================plotbands.in=============
bandx.dat
-6,3
plotbands.xmgr
plotbands.ps
18.1295
1.0,0
===================
thanks
venkatesh
On Sat, Sep 8, 2018 at 6:29 PM venkatesh chandragiri <
venkyphysicsiitm at gmail.com> 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 (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/20180909/4fdf6828/attachment.html>
More information about the users
mailing list