BSE spectra considering momentum

You can find here problems arising when using old releases of Yambo (< 5.0). Issues as parallelization strategy, performance issues and other technical aspects.

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

pangrt
Posts: 26
Joined: Thu Jun 18, 2020 1:38 pm

BSE spectra considering momentum

Post by pangrt » Sun Mar 14, 2021 1:36 am

Dear all,
When calculating the BSE containing momentum, for example,the Absorption @ Q(99) iku: 0.281250030 0.453125030 0.00000000
in the output file o.eps_q99_diago_bse corresponds to the transition from which point in the VB to which point in the CB, and what rules are used to arrange it ?
Thanks!
Pangrt
RongTianPang, Inner Mongolia University, Hohhot, China :)

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

Re: BSE spectra considering momentum

Post by Daniele Varsano » Sun Mar 14, 2021 12:16 pm

Dear pangrt,
please add your full name in your signature, it is a rule of the forum.

The excitation energies from which the spectrum is build are calculated by diagonalization of the BSE equations so different transitions do participate with a certain weight (i.e. it is not an independent particle spectrum). The BSE matrix is built considering all the transition in the BZ connecting valence bands at (k-q) to conduction bands at k (v(k-q)-->c(k)).

You can inspect the weight of each transition to a given excitation using the ypp post processing:

ypp -e a

In order to do that, you need to store the eigenvector of the BSE solution, so remember you need to add the flag WRbsWF in your yambo BSE diagonalization input file.

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/

pangrt
Posts: 26
Joined: Thu Jun 18, 2020 1:38 pm

Re: BSE spectra considering momentum

Post by pangrt » Sun Mar 14, 2021 2:26 pm

Dear Daniele,
Thank you very much for your help. But when I use ypp -e a to generate the file o.exc_qpt12_weights_at_31, the content of the file is
# Band_V Band_C Kv-q ibz Symm_kv Kc ibz Symm_kc Weight Energy
#
72.0000000 73.0000000 96.0000000 1.00000000 39.0000000 3.00000000 0.135884881 1.41313267
72.0000000 73.0000000 97.0000000 1.00000000 25.0000000 3.00000000 0.105637647 1.42940605
72.0000000 73.0000000 91.0000000 1.00000000 51.0000000 3.00000000 0.990237594E-1 1.41552889.
while the smallest direct band gap is at the 102nd K point.I am confused if the transition from v(96)->c(39) is included, then why there are only 102 absorption spectra instead of 102×102.
Thanks!
Pangrt
Last edited by pangrt on Mon Mar 15, 2021 1:41 am, edited 1 time in total.
RongTianPang, Inner Mongolia University, Hohhot, China :)

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

Re: BSE spectra considering momentum

Post by Daniele Varsano » Sun Mar 14, 2021 3:33 pm

Dear RongTianPang,

Here the output shows that the excitation number 31 at qpt =12 is composed by several transitions.
The most relevant (weight=0.13) is the transition from a valence (72) band at kpt=96 to conduction (73) at kpt=39 rotated by the symmetry operation number 3. Then, also other transitions contribute.
The kpt and symmetry operation are reported in the report file.
sorry, but I cannot understand your doubts, can you reformulate in a more clear way?

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/

pangrt
Posts: 26
Joined: Thu Jun 18, 2020 1:38 pm

Re: BSE spectra considering momentum

Post by pangrt » Mon Mar 15, 2021 1:59 am

Dear Daniele,
Thank you very much, with your help, I have understood.

Thanks!
Pangrt
RongTianPang, Inner Mongolia University, Hohhot, China :)

pangrt
Posts: 26
Joined: Thu Jun 18, 2020 1:38 pm

Re: BSE spectra considering momentum

Post by pangrt » Tue Mar 16, 2021 3:46 am

Dear Daniele,
I cannot understand why the exciton oscillator is several orders of magnitude stronger when considering the transition of momentum. In the direct transition, my calculation result is displayed in the file o.exc_qpt1_E_sorted as

Code: Select all

#    Maximum Residual Value = .84226E-01
#   
#    E [ev]             Strength           Index
#
     1.76183784        0.138285998E-1      1.00000000
     1.76186168        0.133084534E-1      2.00000000
     1.79497492        0.282442314E-2      3.00000000
     1.79499960        0.291901291E-2      4.00000000
     1.8261855         0.31177094E-09      5.0000000
After considering the momentum, the oscillator strength of the exciton is displayed in the file o.exc_qpt2_E_sorted as

Code: Select all

#    Maximum Residual Value = .23171E+06
#   
#    E [ev]             Strength           Index
#
     1.79297149        0.140369602E-2      1.00000000
     1.79298830        0.298437499E-2      2.00000000
     1.82057118        0.215879246E-2      3.00000000
     1.82070005        0.312032588E-1      4.00000000
Even this intensity will reach#
# Maximum Residual Value = .54333E+09.
I don’t know if my understanding is wrong or the calculation result is a problem.
Thanks!
Pangrt
RongTianPang, Inner Mongolia University, Hohhot, China :)

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

Re: BSE spectra considering momentum

Post by Daniele Varsano » Tue Mar 16, 2021 7:22 am

Dear RongTianPang,

thanks for reporting, there could be a prefactor problem, I will investigate it.
Here the reason is that generalised dipoles are different from the q=0 and finite q case, being q.<v|r|c> for q=0 (dipole approximation), and <v|e^iq.r|c> for a finite q case. I need to check the prefactors in the finite q case. In any case it is just a multiplication factor. Let me check.

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/

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

Re: BSE spectra considering momentum

Post by Daniele Varsano » Wed Mar 17, 2021 9:16 am

Dear RongTianPang,
indeed there is a problem we will fix soon.

The q=0 max strength is correct, while for finite q the pre-factor is wrong.
Note the values reported are normalised to the maximum value so the relative intensity for each q is correct.
The absolute value instead should be modified as S=S*10^-10/|q|^2

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/

pangrt
Posts: 26
Joined: Thu Jun 18, 2020 1:38 pm

Re: BSE spectra considering momentum

Post by pangrt » Wed Mar 17, 2021 10:46 am

Dear Daniele,
Thank you for solving my doubts very quickly!

Thanks!
Pangrt
RongTianPang, Inner Mongolia University, Hohhot, China :)

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

Re: BSE spectra considering momentum

Post by Daniele Varsano » Wed Mar 17, 2021 10:49 am

Dear Pangrt,
please note I wrote the wrong sign in the exponent, now corrected.

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/

Locked