It is currently Fri Dec 06, 2019 5:34 pm

Tell a friend!


Post new topic Reply to topic  [ 2 posts ] 
phonon assisted absorption 
Author Message

Joined: Tue Jun 11, 2019 1:53 am
Posts: 2
University: The University of To
Post phonon assisted absorption
Hello,

For the phonon assisted absorption calculations, in the PRL paper (Phys. Rev Lett. 108, 167402 (2012)), the generalized optical matrix elements are:

S1(k,q) = Sum_m v_im(k)*g_mj,nu(k,q) / (e_mk - e_ik - hbar*omega + i*eta)
S2(k,q) = Sum_m g_im,nu(k,q)*v_mj(k+q) / (e_mk+q - e_ik +/- hbar*omega_nu,q + i*eta)

the plus and minus signs in the above equation are for phonon emission and absorption.

As of the implementation, in the source indabs.f90 file, the matrix elements are:

s1a(:) = s1a(:) + epf(mbnd, jbnd,imode) * 0.5 * vkk(:,ibnd, mbnd) / ( ekmk - ekq + wq(imode) + ci * eta(m) )
s1e(:) = s1e(:) + epf(mbnd, jbnd,imode) * 0.5 * vkk(:,ibnd, mbnd) / ( ekmk - ekq - wq(imode) + ci * eta(m) )
s2a(:) = s2a(:) + epf(ibnd, mbnd,imode) * 0.5 * vkq(:,mbnd, jbnd) / ( ekmq - ekk - wq(imode)+ ci * eta(m))
s2e(:) = s2e(:) + epf(ibnd, mbnd,imode) * 0.5 * vkq(:,mbnd, jbnd) / ( ekmq - ekk + wq(imode)+ ci * eta(m))

As we can see, the S2 agrees with each other, however, S1a and S1e are not the S1 in the PRL paper. Discrepancy can also be found for the carrier and phonon statistic function P = (n_nu,q + 1/2 +/- 1/2)(f_ik - f_jk+q) in the paper. For the implementation, it is however:

pfac = nqv(imode) * wgkk *(one- wgkq ) - (nqv(imode)+one)*(one-wgkk) * wgkq
pface = (nqv(imode)+one) * wgkk *(one- wgkq ) - nqv(imode) *(one-wgkk) * wgkq

for phonon absorption and emission, respectively. I really appreciate if anyone could give some comments and suggestions.

Best,
Xiaoming Wang


Sat Jun 15, 2019 2:35 pm
Profile E-mail

Joined: Tue Jun 11, 2019 1:53 am
Posts: 2
University: The University of To
Post Re: phonon assisted absorption
Hi,

Btw, anyone know what the '0.5' comes from in the s1a s2a s1e and s2e equations?

Best,
Xiaoming


Tue Jun 18, 2019 6:16 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:  
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