Hello everyone,
Iam doing DFT calculations on ZnO wurtzite, means hexagonal structure.
It is a wide gap semiconductor with gap= 3.4 eV.
the value a=6.239 bohr, c= 9.8379 bohr
So I use celldm(1)=6.239
celldm(3)=1.6048 // c/a
So to find out celldm(1) using espresso I varied the celldm(1) from 5.1 to 7 in steps of 0.1 but I doesnt get minima anywhere, so I then vary the celldm(1) from 5 to 10 in step of1 to check the minima value. But I havnt get minima but my energy gets converge after 7bohr.
I checked my input file many times but doesnt able to find my mistake. Please check the graph which I got using Zno.scf.in and please guide me what Iam doing wrong.?
input file:
&control
calculation='scf',
restart_mode='from_scratch',
outdir='./tmp',
prefix='zno',
pseudo_dir='/home/abhi/src/espresso-5.2.1/pseudo/',
/
&system
ibrav=4,
celldm(1)=lp, celldm(3)=1.6048 // I varied lp from 5 to 10 IN STEP OF 1
nat=2,
ntyp=2,
ecutwfc = 40,
ecutrho = 240
nspin=1,
occupations='fixed',
/
&electrons
mixing_beta=0.7,
diagonalization='david'
/
ATOMIC_SPECIES
Zn 65.39 Zn.pbe-van.UPF
O 15.9994 O.pbe-van_ak.UPF
ATOMIC_POSITIONS alat
Zn 0.3333 0.6667 0.0000
O 0.3333 0.6667 0.3750
K_POINTS automatic
10 10 10 0 0 0
Thanks,
Abhinav Nag