It is currently Sun Jan 20, 2019 10:47 am

Tell a friend!

Post new topic Reply to topic  [ 3 posts ] 
graphene obtains a "band gap" in the spectral function 
Author Message
User avatar

Joined: Thu Jul 14, 2016 6:09 am
Posts: 31
University: POSTECH (Korea)
Post graphene obtains a "band gap" in the spectral function
Dear all,

I am facing a weird problem when calculating the spectral function for graphene.

I have first converged the line-width of the electrons which works fine at an nscf-mesh of 33x33x1 and a q-point mesh of 11x11x1. I interpolate at a fine q mesh of 500x500x1 which gives very nice smooth linewidth plots.

however when I calculate the spectral function the material suddenly has a sizable "gap" at the K point of ~0.250 eV. I am curious where this might arise?

I used a fine mesh in the "path.dat" file close to K

        0.3318285135    0.3318285135    0.0000000000  0.00
        0.3323290090    0.3323290090    0.0000000000  0.00
        0.3328295045    0.3328295045    0.0000000000  0.00
        0.3333300000    0.3333300000    0.0000000000  0.00
        0.3333333333    0.3333333333    0.0000000000  0.00
        0.3338305105    0.3323290090    0.0000000000  0.00
        0.3343310210    0.3313280180    0.0000000000  0.00
        0.3348315315    0.3303270270    0.0000000000  0.00
        0.3353320420    0.3293260360    0.0000000000  0.00

the interpolated bands and phonons look excellent (I use the common trick with two carbon atoms)

  proj(1)     = 'C1:sp2;pz'
  proj(2)     = 'C2:pz'

and I checked "by hand" that the Fermi level is properly located at the "Dirac point". Increasing the k and q grid did not change the situation so I think I am stuck at something else - does anyone know what it might be?

Thanks in advance for your help!


Mon Dec 10, 2018 1:32 pm
Profile E-mail
Site Admin
User avatar

Joined: Wed Jan 13, 2016 7:25 pm
Posts: 534
University: Oxford
Post Re: graphene obtains a "band gap" in the spectral function

Well the underlying perturbation theory is by construction ill defined for metals/semi-metals.

The real part of the el-ph self energy corresponds to the quasi-particle shift with respect to KS. This is responsible for the zero-point motion.
By definition the spectral function contains the real-part of the self-energy (not only the imaginary part).

As a result this will open a gap.

See for example PRB 92, 085137 (2015) for more info on spectral functions.

An add-hoc "solution" is just to shift back the spectral function so that both end touch. This is off course a big approximation and one need to be careful.

Hope this helps.

Dr. Samuel Poncé
Department of Materials
University of Oxford
Parks Road
Oxford OX1 3PH, UK

Thu Dec 13, 2018 6:32 pm
Profile E-mail
User avatar

Joined: Thu Jul 14, 2016 6:09 am
Posts: 31
University: POSTECH (Korea)
Post Re: graphene obtains a "band gap" in the spectral function
Hi Samuel,

thanks a lot for your reply! Interestingly it did eventually work out quite well by doubling the coarse k-mesh - the resulting spectral function has a nice Dirac-cone like shape. It also helped to determine the fermi level more precisely using tetrahedron method in a separate calculation and read it later in EPW. I am not sure if that is a "cancellation of error" kind of situation but the results look good and are comparable to literature!

Hope that helps the next guy who stumbles across graphene :)


Mon Jan 07, 2019 1:33 pm
Profile E-mail
Display posts from previous:  Sort by  
Post new topic Reply to topic  [ 3 posts ] 

Who is online

Users browsing this forum: No registered users and 2 guests

You cannot post new topics in this forum
You cannot reply to topics in this forum
You cannot edit your posts in this forum
You cannot delete your posts in this forum

Search for:
Jump to:  
Powered by phpBB © phpBB Group.
Designed by Vjacheslav Trushkin

(All content on this board is governed by and is the sole responsibility of the board administrator.)

Gratis forum Free forum hosting| gratis phpbb3 forum | phpbb3 styles