questions about the results of BSE calculations?

Deals with issues related to computation of optical spectra in reciprocal space: RPA, TDDFT, local field effects.

Moderators: Davide Sangalli, andrea.ferretti, myrta gruning, andrea marini, Daniele Varsano, Conor Hogan

pyadav
Posts: 86
Joined: Thu Nov 26, 2020 2:56 pm
Contact:

Re: questions about the results of BSE calculations?

Post by pyadav » Sat Feb 04, 2023 9:09 am

Dear team Yambo,

Can someone please comment on the following two points?

1. How can I convert the RL numbers (wavevectors for energy ) into Rydberg. Any formula or any rough estimation method?
2. The dipole socillator strength reported in file 'o.exc_qpt1_E_sorted' is absolute or relative to the largest dipole oscillator strength?

Thank you,
Pushpendra
Pushpendra Yadav
Ph.D. Research Scholar
Quantum Transport and Theory Group
Department of Physics
Indian Instittute of Technology Kanpur, India.

https://sites.google.com/site/amitkag1/

User avatar
Daniele Varsano
Posts: 3773
Joined: Tue Mar 17, 2009 2:23 pm
Contact:

Re: questions about the results of BSE calculations?

Post by Daniele Varsano » Mon Feb 06, 2023 10:28 am

Dear Pushendra,

1) the relation is Number of RL such that G^2/2 < Ecut. I do not have a rough estimation. Anyway, when you assign an energy cutoff you have the number of RL in the report. Please note that instead of Ha, or Ry, you can also use mHa (mRy).

2) They are normalized and you have the normalization factor at the top of the output file:

Code: Select all

Maximum Residual Value = 
If you want to change it you can modify the source in ./ypp/excitons/excitons_sort_and_report.F
setting:

Code: Select all

R_normalize="no"
Best,
Daniele
Dr. Daniele Varsano
S3-CNR Institute of Nanoscience and MaX Center, Italy
MaX - Materials design at the Exascale
http://www.nano.cnr.it
http://www.max-centre.eu/

pyadav
Posts: 86
Joined: Thu Nov 26, 2020 2:56 pm
Contact:

Re: questions about the results of BSE calculations?

Post by pyadav » Mon Feb 06, 2023 1:36 pm

Dear Daniele,

Thank you so much for the reply and clarifications. They are really helpful.

Wishes,
Pushpendra.
Pushpendra Yadav
Ph.D. Research Scholar
Quantum Transport and Theory Group
Department of Physics
Indian Instittute of Technology Kanpur, India.

https://sites.google.com/site/amitkag1/

Dean
Posts: 98
Joined: Thu Oct 10, 2019 7:03 am

Re: questions about the results of BSE calculations?

Post by Dean » Thu Aug 17, 2023 3:28 am

Dear Daniele,
After setting "Weight_treshold=0.0 eV", I get all the transitions in the file xxx.exc_qpt1_weights_at_1.
But I find that the weights for one K point in the first part of outfile do not equal the sum of weights in the second part of outfile, such as Gamma point (0.144) vs (0.02+......).
Please see the attached file.
Best wishes,
You do not have the required permissions to view the files attached to this post.
Dr. Yimin Ding
Soochow University, China.

User avatar
Daniele Varsano
Posts: 3773
Joined: Tue Mar 17, 2009 2:23 pm
Contact:

Re: questions about the results of BSE calculations?

Post by Daniele Varsano » Fri Aug 25, 2023 9:35 am

Dear Yimin,

this is just because the weight reported summed up on transitions for each K point and the individual contributions are normalized differently, i.e. the K point contribitions are normalized to its maximum:
see line 154 in
ypp/excitons/excitons_amplitudes.F

Code: Select all

148    K_weight=0._SP
149    do neh = 1,BS_H_dim
150      ikbz  = BSS_eh_table(S_indx(neh),1)
151      ikibz = Xk%sstar(ikbz,1)
152      K_weight(ikibz)=K_weight(ikibz)+A_weight(neh)
153    enddo
154    K_weight=K_weight/maxval(K_weight)
155    titles(1:4)=(/'             ','K-point [iku]','             ','Weight       '/)
156    call msg('o weight','#',titles(:4),INDENT=0,USE_TABS=.true.)
157    do ikibz=1,Xk%nibz
158      if (K_weight(ikibz)>EXCITON_weight_treshold) then
159        call msg('o weight','# ',(/Xk%pt(ikibz,:),K_weight(ikibz)/),INDENT=0,USE_TABS=.true.)
160      endif
161    enddo
162    call msg('o weight','#','',INDENT=0,USE_TABS=.true.)
Indeed, they do not sum-up to 1, but they indicate that:
# : 0.375000000 0.714285672E-1 0.00000000 1.00000000
has the maximum contribution.

Best,
Daniele
Dr. Daniele Varsano
S3-CNR Institute of Nanoscience and MaX Center, Italy
MaX - Materials design at the Exascale
http://www.nano.cnr.it
http://www.max-centre.eu/

Dean
Posts: 98
Joined: Thu Oct 10, 2019 7:03 am

Re: questions about the results of BSE calculations?

Post by Dean » Mon Aug 28, 2023 10:00 am

Dear Daniele,
Thanks for your reply.
Best,
Dr. Yimin Ding
Soochow University, China.

Post Reply