Moderators: Daniele Varsano, andrea.ferretti, andrea marini, Conor Hogan, myrta gruning
593 if (BS_blk(i_block)%mode=="R") BS_blk(i_block)%mat(i_Tk,i_Tp)= sqrt_f_itk*(H_x*Co-H_c)*sqrt_f_itp
594 if (BS_blk(i_block)%mode=="A") BS_blk(i_block)%mat(i_Tk,i_Tp)=cI*sqrt_f_itk*(H_x*Co-H_c)*sqrt_f_itp
write(unit,*) i_k_bz, i_v_k , i_c_k
write(unit,*) i_p_bz, i_v_p , i_c_p
write(unit,*) real(H_x*Co),aimag(H_x*Co),real(H_c),aimag(H_c)
1. Could you please point out which command (or calculation) will call the module for outputting the Columb matrix?
2. If I want the screened electron-hole interaction matrix, H_c would be enough, right? Don't need exchange part, Co?
3. Is the Columb matrix k-point dependent? i.e., the code will output the Columb matrix between KS orbitals at every k-point, right? Actually, I just need the Columb matrix at Gamma point right now.
what's the unit of the elements of Columb matrix, eV, or Ry, and or Ha?
H_x, Co, H_c, i_p_bz, i_v_p , i_c_p?
could you please point out the right way where I need to make the changes, and if possible, how to make the changes.
I do not know what you mean at gamma point, you have two k indexes (k,k'), if you want only the diagonal part k=0 and k'=0 just insert and if condition for i_k_bz=0 and i_p_bz=0
But I just need the Columb matrix at Gamma k-point from YAMBO output, we can do this right?
Just not sure if we could extract the Columb matrix between the KS orbitals at the interested k-point.
if (i_k_bz==1.and.i_p_bz==1) then
write(unit,*) i_v_k , i_c_k, i_v_p , i_c_p, real(H_c),aimag(H_c)
endif
Users browsing this forum: No registered users and 1 guest