It is currently Mon Oct 14, 2019 6:22 am

Tell a friend!


Post new topic Reply to topic  [ 2 posts ] 
Negative phonon linewidth !! 
Author Message

Joined: Mon Jul 11, 2016 11:45 am
Posts: 30
University: Sungkyunkwan Univers
Post Negative phonon linewidth !!
Dear all,

I am curious to know about the flag parameter degaussw!! Is the the same value as in electron-phonon calculations (i.e. degauss in phonon dispersion) or we need to cross check in EPW calculations ??

Also, I am quite surprised and puzzled with negative phonon linewidths in all acoustic modes as well as highest optical mode (occurring at the middle of Gamma to M (0.5,0.0) direction). My wannier spreads are between 1.8 to 1.9 Angstroms for all atoms and my wannier projected band structures in line with standard DFT. Indeed my electron phonon dispersion showing completely positive modes at acoustic level itself. I also checked with different eptemp, and found negative linewidths with different energies. Below I am attaching the part my epw.in for calculating the phonon line widths.

Kindly suggest me for further

Quote:

elecselfen = .false.
phonselfen = .true.
a2f = .false.

parallel_k = .true.
parallel_q = .false.

fsthick = 2.0 ! eV
eptemp = 10 ! K (same as PRB 76, 165108)
degaussw = 0.25 ! eV

dvscf_dir = '../phonons/save'
filukk = './$prefix.ukk'
filqf = '/home/saint/Work/PWSCF//EPW/New/meshes/path.dat'
nkf1 = 16
nkf2 = 16
nkf3 = 1

nk1 = 8
nk2 = 8
nk3 = 1

nq1 = 8
nq2 = 8
nq3 = 1



with regards
S. Appalakondaiah


Mon Oct 03, 2016 5:33 am
Profile E-mail
User avatar

Joined: Thu Jan 14, 2016 10:52 am
Posts: 151
University: University of Vienna
Post Re: Negative phonon linewidth !!
Dear S. Appalakondaiah,

You should test your results with respect to degaussw, as it corresponds to the broadening of the delta functions in the calculation of the phonon linewidth.
Regarding the negative phonon linewidths you are getting, I can comment that the formula implemented is not positive definite so it is possible to find negative numbers. This cannot happen if you use the less accurate 'delta-delta approximation' (input delta_approx=.true.).
You could try incresing the fine k-mesh and decreasing degaussw in your calculations and check if the results change.

Best
Carla


Tue Oct 04, 2016 6:19 pm
Profile E-mail
Display posts from previous:  Sort by  
Post new topic Reply to topic  [ 2 posts ] 


Who is online

Users browsing this forum: No registered users and 1 guest


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:  
cron
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