Re: BSE spectra considering momentum
Posted: Wed May 18, 2022 9:49 am
Sorry I forgot to attach files.
Welcome to the Yambo forum! Post requests, look for help, and discuss the code with the community of users and developers.
https://www.yambo-code.eu/forum/
Code: Select all
384 ! POLARIZABILITY
385 !----------------
386 !
387 if (eval_alpha_loc) then
388 !
389 if(l_col_cut.or.trim(alpha_dim)=="volume") then
390 alpha_factor=eps_2_alpha/4._SP/pi*bare_qpg(iq,1)**2/q_norm(iq)**2
391 alpha(:,1)=Epsilon_ii(:,2)*alpha_factor
392 alpha(:,2)=Epsilon_ii(:,3)*alpha_factor
393 if (l_print_col4) alpha(:,3)=Epsilon_ii(:,4)*alpha_factor
394 else
395 alpha_factor=eps_2_alpha/4._SP/pi
396 alpha=(1._SP+eels)*alpha_factor
397 endif
398 !
I think problem in bare_qpg(1,1) and q_norm(1).! POLARIZABILITY
!----------------
!
if (eval_alpha_loc) then
!
if(l_col_cut.or.trim(alpha_dim)=="volume") then
alpha_factor=eps_2_alpha/4._SP/pi*bare_qpg(1,1)**2/q_norm(1)**2
alpha(:,1)=Epsilon_ii(:,2)*alpha_factor
alpha(:,2)=Epsilon_ii(:,3)*alpha_factor
if (l_print_col4) alpha(:,3)=Epsilon_ii(:,4)*alpha_factor
else
alpha_factor=eps_2_alpha/4._SP/pi
alpha=(1._SP+eels)*alpha_factor
endif