I use the yambo 4.5.3 (sorry for no upgrade because I have no authority in our cluster). As an example, I run a spin polarization calculation on two-dimensional CN system using the PW program. As shown in PW output (see attached file yb_pwscf.out), the system has 9 electrons with 5 up and 4 down. Then I run the BSE calculation and modify the code of ypp (excitons_driver.F) to print out the BSS_eh_table (see attached file YBBSEEHT.txt).
Code: Select all
do I = 1, BSS_n_eig
if(n_sp_pol==1) write(200,'(1X,3I4," 1 ")') BSS_eh_table(I,:)
if(n_sp_pol==2) write(200,'(1X,4I4)') BSS_eh_table(I,:)
enddo
Code: Select all
k v c spin
1 4 5 1
1 4 5 2
1 4 6 1
1 4 6 2
1 4 7 1
1 4 7 2
Code: Select all
% BSEBands
4 | 7 | # [BSK] Bands range
Could you give me some hints or references for constructing the BSE matrix in the spin polarization calculation?
Best regards.
Youzhao Lan