Re: electron-phonon matrix elements

Dear Fairoja,

You can consider any state. We only show the el-ph matrix element at k=\Gamma and then changed q.

However you can look at any k-point.

In the code you have:

**Code:**

! consider only initial state k=0

ik = 1

WRITE(6,'(a,3f10.6)') 'xk ', (xk(n,ik),n=1,3)

where ik=1 is in this case the Gamma point. You can just change that value to get different k-points.

Note that you can also add k-point loop around all of that like this

**Code:**

DO ik=1, MAX_NUMBER_OF_K_POINTS

WRITE(6,'(a,3f10.6)') 'xk ', (xk(n,ik),n=1,3)

write(6,*) ' ibnd jbnd imode eig_i (eV) eig_j (eV) omega_nu (meV) |g| (meV)'

...

ENDDO

Be careful that this will generate very large txt files (lots of data).

Best wishes,

Samuel