Getting Bandgap from DOS calculated by GW method

Queries about input and output files, running specific calculations, etc.


Moderators: Global Moderator, Moderator

Post Reply
Message
Author
felixvasp

Getting Bandgap from DOS calculated by GW method

#1 Post by felixvasp » Mon May 16, 2011 9:16 am

Hi everyone,
I am newbie to VASP,especially DOS&BAND structure calculations. I ran the Si diamond bulk DOS calculation using GW method by two steps:

1. ran DFT calculation on Si bulk to get WAVEDER.
INCAR:
System = diamond Si
NBANDS=160
ENCUT=250
ISMEAR = 0; SIGMA = 0.05;
LOPTICS=.TRUE.

2. make another directory, copy POTCAR, KPOINT, POSCAR to the diretory, and ran GW calculations using WAVEDER from step 1 and using INCAR as follows:
INCAR
system = diamond Si
NBANDS=160
ENCUT=250
ISMEAR = 0; SIGMA = 0.05;
LOPTICS=.TRUE.
ALGO=GW
LWAVE=.TRUE.

3. plot DOS from DOSCAR.
But the calculated DOS graph doesn't give the BANDGAP of ~1.2 eV reported by literature but gives a BANDGAP of ~0.5 eV. Is there anything I didn't wrong? but someone reproduced the bandgap of GW by bandstructure. Can DOS produced by GW calculation gives right bandgap or if I did something wrong during producing DOS for GW, how to get the right DOS from GW calculation. Thanks a lot for help.






<span class='smallblacktext'>[ Edited ]</span>
Last edited by felixvasp on Mon May 16, 2011 9:16 am, edited 1 time in total.

felixvasp

Getting Bandgap from DOS calculated by GW method

#2 Post by felixvasp » Tue May 17, 2011 9:59 pm

Anyone has the answer or I need to provide more specific inf.?
Last edited by felixvasp on Tue May 17, 2011 9:59 pm, edited 1 time in total.

Post Reply