Absolute value of the imaginary part of the dielectric function

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

Quxiao
Posts: 96
Joined: Fri Mar 26, 2021 11:27 am

Re: Absolute value of the imaginary part of the dielectric function

Post by Quxiao » Wed Jan 19, 2022 1:25 pm

Daniele:
I know that there are three options for BSEmod to choose from, which seem to represent three modes: resonant/retared/coupling. Regarding your previous reply, my understanding is that the difference between retared and resonant is the real part of the response function, but this is very important for the final Does the calculation of the dielectric function of , also have an impact? Because I saw the input of other people in the forum, both resonant and retarted, about the accurate application of these three modes, whether it is displayed in the learn module or other links, maybe due to my negligence, I did not find this introduction. .
Finally, if I calculate the normal BSE, should I choose resonant or retarted? Or both?
Thanks again for your selfless help!
Best wishes
Quxiao
BIT
Quxiao in BIT,calculate the exciton

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

Re: Absolute value of the imaginary part of the dielectric function

Post by Daniele Varsano » Wed Jan 19, 2022 2:33 pm

Dear Quixiao,
the definition of the time ordering of Green functions can be found in textbooks. In the retarded case the antiresonant part is also considered. When considering the coupling also the coupling between resonant and antiresonant is taken into account and in this case, the calculation is more cumbersome and usually can be neglected (Tamm Dancoff approximation). As I wrote you in the previous post you can consider the retarded mode to have an accurate real and imaginary part. Of course, you can verify that the Tamm Dancoff approximation applies by considering the coupling terms.

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/

Quxiao
Posts: 96
Joined: Fri Mar 26, 2021 11:27 am

Re: Absolute value of the imaginary part of the dielectric function

Post by Quxiao » Fri Jan 21, 2022 4:36 am

Dear Daniele:
Thanks for your help again!
I still have a question. I found that the absolute value of the imaginary part of the dielectric constant for the BSE I calculated with the yambo4.4 system is much smaller than the absolute value of the yambo5.0. Although the trend is the same, the absolute value of this The difference makes me wonder, is there a parameter I haven't noticed that controls this?
This is the result of the same system I calculated with yambo4.4 and 5.0 respectively, except for the absolute value of the imaginary part of the dielectric function, the difference is not big,please give me some help again.
Thanks for your selfness!
Best wishes
Quxiao
BIT
You do not have the required permissions to view the files attached to this post.
Quxiao in BIT,calculate the exciton

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

Re: Absolute value of the imaginary part of the dielectric function

Post by Daniele Varsano » Fri Jan 21, 2022 9:46 am

Dear Quxiao,
is this a 2D system? In this case the dielectric function is not well defined and the relevant quantity is the polarisability alpha.
When using the cutoff coulomb technique we have modified the response function consistently and this is the reason in the change of absolute value as it is used the truncated potential instead than the 1/q^2 term.

This is explained in the FAQ of the Yambo webpage:
https://www.yambo-code.org/faq/

Note that the website in these days is very slow, we are working to solve the problem.

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/

Quxiao
Posts: 96
Joined: Fri Mar 26, 2021 11:27 am

Re: Absolute value of the imaginary part of the dielectric function

Post by Quxiao » Sat Feb 12, 2022 5:23 pm

Dear Daniele:
Thanks for your reply!I have learned the link you provided,something I want to sure,
When I use the cutoff coulomb technique, the (polarisability) alpha=-X/q^2,is it?
When I use the yambo -r (cutoff coulomb technique),I could get the polarisability from the o-yambo.out.alpha_q1_diago_bse,is it?
And I also could get the imaginary part of dielectric function from the o-yambo.out.eps_q1_diago_bse ,is it?
I see that formulation:ε = 1 +4πα/d,the alpha(α)=-X/q^2?
The last question is that I have seen that there is a unit of polarizability (length unit) in some literature, so I want to know what is the unit of the polarizability output in yambo?
I am confused,please give me some help
Thanks!
Best wishes
Quxiao
Quxiao in BIT,calculate the exciton

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

Re: Absolute value of the imaginary part of the dielectric function

Post by Daniele Varsano » Sat Feb 12, 2022 5:43 pm

Dear Quxiao,
When I use the cutoff coulomb technique, the (polarisability) alpha=-X/q^2,is it?
Yes! but there is a multiplying factor that depends on the dimensionality. For 2D it is multiplied by the non periodic length "d"
And I also could get the imaginary part of dielectric function from the o-yambo.out.eps_q1_diago_bse ,is it?
Yes, it will be proportional to alpha but with very small numbers, going to zero for large volumes.
The last question is that I have seen that there is a unit of polarizability (length unit) in some literature, so I want to know what is the unit of the polarizability output in yambo?
Correct, alpha has the dimension of a length unit and it is in bohr

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/

Quxiao
Posts: 96
Joined: Fri Mar 26, 2021 11:27 am

Re: Absolute value of the imaginary part of the dielectric function

Post by Quxiao » Mon Feb 21, 2022 9:09 am

Dear Daniele:
Thanks for your help!
I learn that with your help!
Best wishes!
Quxiao
Quxiao in BIT,calculate the exciton

Quxiao
Posts: 96
Joined: Fri Mar 26, 2021 11:27 am

Re: Absolute value of the imaginary part of the dielectric function

Post by Quxiao » Tue Feb 22, 2022 4:56 pm

Dear Daniele:
I am a newcomer and have a lot of confusion.I read the yambopy code and successfully realized the exciton projection on the single electron energy band, but now there is a problem. I previously made the orbital projection of the single electron energy band and used the group theory to conclude that there is a transition between the lowest conduction band and the highest valence band. However, the exciton projection I made with yambopy is in the lowest conduction band and the second highest valence band, There is no exciton in the highest valence band, so I doubt the result.

If the projection result is correct, the formation of excitons is higher than the transition selection rule. What's more, I don't understand the physical meaning of the weight of excitons in the energy band projection, because I know that excitons can expand in the electron hole state. Is this weight the expansion coefficient, the expansion coefficient of the electron hole state with specified energy?
Thanks for your help!
Best wishes
Quxiao
Quxiao in BIT,calculate the exciton

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

Re: Absolute value of the imaginary part of the dielectric function

Post by Daniele Varsano » Wed Feb 23, 2022 8:48 am

Dear Quixiao,

the exciton weights are the square modulus of the eigenstate of the BSE which is calculated in transition space.
So,the exciton can be expressed as a linear combination of independent particle transitions:
\Psi_exc=\sum A_cv \psi_v \psi_c
the weights are |A_cv|^2

Note that the energy reported in the weight file, transition energies (Ec_Ev) are the KS energies and not the QP energies.

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/

Quxiao
Posts: 96
Joined: Fri Mar 26, 2021 11:27 am

Re: Absolute value of the imaginary part of the dielectric function

Post by Quxiao » Fri Feb 25, 2022 1:31 pm

Dear Daniele :
I understand what you said .However,I still have some questions.Please give me some help.
When I calculate the BSE of anisotropic materials, I can get the absorption spectrum in different directions by setting BlongDir.
Now I have some doubts about this BlongDir parameter. My previous understanding is the propagation direction (momentum
direction) of the incident light, but I am now suspicious, the note here is # [BSS] [cc] Electric Field,
so could this be referring to the direction of vibration of the electric field (shear wave vibration)?
Or is it the direction of propagation of the electric field (same as that of the incident light)?
This is very important to me, looking forward to your reply!
Best wishes!
Quxiao
Quxiao in BIT,calculate the exciton

Locked