[Pw_forum] Not getting right Band structure
Neel Singh
neelphysics at yahoo.in
Wed Jul 29 10:35:54 CEST 2009
Hi sir,
Sir i have calculated the band structure of YVO4 and it is coming as near about 0.5 ev .. sir as suggested by you in earlier reply that
>means band gap at gamma = 0.3058-(-8.9833)= 8.2891
> which is very large ??
Not necessarily, the offset of band in a periodic system is arbitrary. You
need to know where you Fermi energy is, or which are the HOMO and LUMO of
your system to find the gap.
I have calculated Fermi energy and it is coming as 6.108. The band gap of YVO4 is 3.9 eV experimental and DFT calculates the band gap as nearly half i..e. it should come as 1.8 eV approx. Please help.. Thanks In advance.My Input files are listed below...
SCF INPUT
===============================================
&control
calculation='scf',
restart_mode='from_scratch',
prefix='yvo'
pseudo_dir = '/home/physics/pseudo/',
outdir='/home/physics/work/yvo/temp'
/
&system
ibrav =1, a=7.118300514, b= 7.118300514, c=6.289300454
, nat =24, ntyp= 3, ecutwfc = 40.0,ecutrho=300,
occupations='smearing', smearing='mp', degauss=0.01,
/
&electrons
mixing_beta = 0.7
diagonalization = 'david'
conv_thr = 1.0e-6
/
&ions
/
ATOMIC_SPECIES
Y 88.90585 Y.pbe-nsp-van.UPF
V 50.9415 V.pbe-sp-van.UPF
O 15.9994 O.pbe-van_ak.UPF
ATOMIC_POSITIONS {angstrom}
Y 0.000000000 5.338725356 0.786162553
/
/
K_POINTS {automatic}
4 4 4 0 0 0
=======================================================
Bands Input
======================================================
&control
calculation='bands'
pseudo_dir = '/home/physics/pseudo/',
outdir='/home/physics/work/yvo/temp',
prefix='yvo',
/
&system
ibrav =1, a=7.118300514, b= 7.118300514, c=6.289300454, nat=24, ntyp=3,
ecutwfc = 40.0,ecutrho=300, occupations='smearing', smearing='mp', degauss=0.01, nbnd= 115,
/
&electrons
mixing_beta=0.7
conv_thr=1.0e-6
diagonalization = 'david'
/
&ions
/
ATOMIC_SPECIES
Y 88.90585 Y.pbe-nsp-van.UPF
V 50.9415 V.pbe-sp-van.UPF
O 15.9994 O.pbe-van_ak.UPF
ATOMIC_POSITIONS {angstrom}
Y 0.000000000 5.338725356 0.786162553
\
\
K_POINTS
26
0.000 0.000 0.000 1
0.000 0.000 0.125 1
0.000 0.000 0.250 1
0.000 0.000 0.375 1
0.000 0.000 0.500 1
0.125 0.125 0.500 1
0.250 0.250 0.500 1
0.375 0.375 0.500 1
0.500 0.500 0.500 1
0.500 0.500 0.333 1
0.500 0.500 0.166 1
0.500 0.500 0.000 1
0.375 0.375 0.000 1
0.250 0.250 0.000 1
0.125 0.125 0.000 1
0.000 0.000 0.000 1
0.000 0.000 0.166 1
0.000 0.000 0.333 1
0.000 0.166 0.500 1
0.000 0.333 0.500 1
0.000 0.500 0.500 1
0.000 0.500 0.333 1
0.000 0.500 0.166 1
0.000 0.500 0.000 1
0.000 0.333 0.000 1
0.000 0.166 0.000 1
============================================================
===========================================================
&inputpp
prefix = 'yvo'
outdir = '/home/physics/work/yvo/temp'
filband = 'yvobands.dat'
/
========================================
Plotband INPUTS
=======================================
yvobands.dat
-10 12
yvobands.xmgr
yvobands.ps
6.1088
1.0 0.00
================================================
============================================
With regards,
Neel Singh
Univeristy of Delhi
--- On Tue, 28/7/09, Lorenzo Paulatto <paulatto at sissa.it> wrote:
From: Lorenzo Paulatto <paulatto at sissa.it>
Subject: Re: [Pw_forum] Not getting right Band structure
To: "PWSCF Forum" <pw_forum at pwscf.org>
Date: Tuesday, 28 July, 2009, 10:07 AM
In data 28 luglio 2009 alle ore 11:54:14, Neel Singh
<neelphysics at yahoo.in> ha scritto:
> means band gap at gamma = 0.3058-(-8.9833)= 8.2891
> which is very large ??
Not necessarily, the offset of band in a periodic system is arbitrary. You
need to know where you Fermi energy is, or which are the HOMO and LUMO of
your system to find the gap.
> 2) In the output of SCF and Bands , the output file has 4 warnings
This a very frequently asked question, please see the corresponing entry
in the manual:
http://www.pwscf.org/guide/2.1.2/html-node/node56.html
http://www.pwscf.org/guide/2.1.2/html-node/node57.html
> for this I have already increase the following quantities in the
> increments
> kinetic-energy cutoff = 40.0000 Ry to 60
> charge density cutoff = 250.0000 Ry to 400
> But still the warnings are there.
It is not a big increment of the cutoff that will solve this problem, but
a small increment or nr1, nr2 nr3
cheers, and good work
--
Lorenzo Paulatto
SISSA & DEMOCRITOS (Trieste)
phone: +39 040 3787 511
skype: paulatz
www: http://people.sissa.it/~paulatto/
*** save italian brains ***
http://saveitalianbrains.wordpress.com/
_______________________________________________
Pw_forum mailing list
Pw_forum at pwscf.org
http://www.democritos.it/mailman/listinfo/pw_forum
See the Web's breaking stories, chosen by people like you. Check out Yahoo! Buzz. http://in.buzz.yahoo.com/
-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://www.democritos.it/pipermail/pw_forum/attachments/20090729/40d20325/attachment-0001.htm
More information about the Pw_forum
mailing list