It is currently Mon Dec 16, 2019 6:59 am

Tell a friend!

Post new topic Reply to topic  [ 3 posts ] 
Negative Scattering Rates 
Author Message

Joined: Tue Sep 19, 2017 12:31 am
Posts: 7
University: Uni. at Buffalo
Post Negative Scattering Rates
I am currently calculating the shortrange scattering rate for alpha-Ga2O3, but scattering rates at some energy values are coming out to be negative. I donot understand what could be the possible reason for this result.
Here is the input:
prefix = 'aGa2O3'
amass(1) = 69.723
amass(2) = 15.999
outdir = './scratch/'

iverbosity = 1

elph = .true.
epbwrite = .true.
epbread = .false.

etf_mem = 0

epwwrite = .true.
epwread = .false.

nbndsub = 1
nbndskip = 44

wannierize = .true.
num_iter = 500
iprint = 4
proj(1) = 'f=0.6435000,0.6435000,0.6435000:s'
dis_froz_min= 12.5287
dis_win_min = 12.5287
wdata(1) = 'num_bands = 45'

efermi_read = .true.
fermi_energy= 12.1866

asr_typ = 'crystal'

lpolar = .true.
shortrange = .true.

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

scattering = .true.
scattering_serta = .true.

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

dvscf_dir = './save/'

filkf = './kpoints.dat'

nqf1 = 20
nqf2 = 20
nqf3 = 20

nk1 = 6
nk2 = 6
nk3 = 6

nq1 = 6
nq2 = 6
nq3 = 6
32 cartesian
0.000000000000000E+00 0.000000000000000E+00 0.000000000000000E+00
-0.177208793294881E+00 -0.102311544511585E+00 0.661588117531448E-01
-0.354417586589762E+00 -0.204623089023169E+00 0.132317623506290E+00
0.531626379884643E+00 0.306934633534754E+00 -0.198476435259434E+00
-0.177208793294881E+00 0.102311544511585E+00 0.132317623506298E+00
-0.354417586589762E+00 0.000000000000000E+00 0.198476435259443E+00
0.531626379884643E+00 0.511557722557924E+00 -0.132317623506281E+00
0.354417586589762E+00 0.409246178046339E+00 -0.661588117531359E-01
0.177208793294881E+00 0.306934633534754E+00 0.888178419700125E-14
-0.354417586589762E+00 0.204623089023169E+00 0.264635247012597E+00
0.531626379884643E+00 0.716180811581093E+00 -0.661588117531270E-01
0.354417586589762E+00 0.613869267069508E+00 0.177080572427712E-13
0.531626379884643E+00 -0.306934633534754E+00 -0.396952870518895E+00
0.000000000000000E+00 0.000000000000000E+00 0.198476435259443E+00
-0.177208793294881E+00 -0.102311544511585E+00 0.264635247012588E+00
0.708835173179524E+00 0.409246178046339E+00 -0.661588117531359E-01
0.531626379884643E+00 0.306934633534754E+00 0.885402862138562E-14
0.354417586589762E+00 0.204623089023169E+00 0.661588117531536E-01
-0.177208793294881E+00 0.102311544511585E+00 0.330794058765742E+00
0.708835173179524E+00 0.613869267069508E+00 0.177358128183869E-13
0.531626379884643E+00 0.511557722557924E+00 0.661588117531625E-01
0.354417586589762E+00 0.409246178046339E+00 0.132317623506307E+00
0.708835173179524E+00 -0.409246178046339E+00 -0.330794058765750E+00
0.531626379884643E+00 -0.511557722557924E+00 -0.264635247012606E+00
0.354417586589762E+00 -0.613869267069508E+00 -0.198476435259461E+00
0.531626379884643E+00 -0.306934633534754E+00 -0.198476435259452E+00
0.000000000000000E+00 0.000000000000000E+00 0.396952870518886E+00
0.886043966474405E+00 0.511557722557924E+00 0.661588117531625E-01
0.708835173179524E+00 0.409246178046339E+00 0.132317623506307E+00
0.886043966474405E+00 -0.511557722557924E+00 -0.264635247012606E+00
0.708835173179524E+00 -0.613869267069508E+00 -0.198476435259461E+00
0.000000000000000E+00 0.000000000000000E+00 -0.595429305778329E+00

The total number of valence bands are 44 and I am just interested in the 1st conduction band for the scattering rate calculation hence nbndskip equals 44.
Any help in this regard will be appreciated.

Thank You,
Ankit Sharma
University at Buffalo

Mon Oct 28, 2019 6:25 pm
Profile E-mail

Joined: Tue Sep 19, 2017 12:31 am
Posts: 7
University: Uni. at Buffalo
Post Re: Negative Scattering Rates
The fermi energy is calculated separately using the DOS calculation from Quantum espresso and for a doping density of 1e17 cm-3(assuming complete ionization). The top of the valence band in around 10 eV and the bottom of the conduction band is around 12.5 eV obtained under the LDA.
For any more information, please let me know.

Thank You,
Ankit Sharma

Tue Oct 29, 2019 4:38 pm
Profile E-mail

Joined: Tue Sep 19, 2017 12:31 am
Posts: 7
University: Uni. at Buffalo
Post Re: Negative Scattering Rates
In the transport.f90 module, the matrix element are read as only the real part of the square of the matrix element in case of shortrange and not Gamma point. I am not able to figure out why this is the case, because in my case some of the shortrange matrix elements are pure imaginary and hence the square is a negative number. I think this might be the reason for the negative shortrange scattering rates.
Any clarification in this regard will be greatly appreciated. Below is the code snippet present in the transport.f90 and selfen_elec.f90 files.

IF ( shortrange .AND. ( abs(xqf (1, iq))> eps8 .OR. abs(xqf (2, iq))> eps8 &
.OR. abs(xqf (3, iq))> eps8 )) THEN
! SP: The abs has to be removed. Indeed the epf17 can be a pure imaginary
! number, in which case its square will be a negative number.
g2 = REAL( (epf17 (jbnd, ibnd, imode, ik)**two)*inv_wq*g2_tmp )
g2 = (ABS(epf17 (jbnd, ibnd, imode, ik))**two)*inv_wq*g2_tmp

Thank You,

Tue Dec 03, 2019 5:18 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