Dear developers：
Firstly, I would like to express my sincere gratitude for your continuous assistance. Through your help, I have gained a deeper understanding of Yambo.
My recent concern is that when solving BSE using the scissor operator, two modifications are needed. They are% KfnQP_S and% XfnQP_S, respectively. My understanding is that the first parameter is equivalent to a direct
translation on the single electron energy level, which will affect the single electron energy gap when constructing BSEkernel: \ epsilonc  \ epsilonv. The second parameter should affect W. Here, my simple understanding is
that the Coulomb interaction W is approximately related to the single electron energy gap.
But I want to know how Yambo applies the parameter% XfnQP_S to the BSE solving process, is it just a value of W, or is there any corresponding formula? Because the values of% KfnQP_S and% XfnQP_S are both set to
the same translation scissor, I speculate that it may be due to the invariance of BSEkernel's diagonal elements? Keep the value of \ epsilonckepsilonvk ' Wkk' unchanged?
I would greatly appreciate it if you could explain and let me know if this operation is in a certain Fortran file.
I hope the Yambo team continues to do better and better！
Best wishes!
Quxiao
Correction of dielectric shielding for scissor operator
Moderators: Davide Sangalli, andrea.ferretti, myrta gruning, andrea marini, Daniele Varsano

 Posts: 99
 Joined: Fri Mar 26, 2021 11:27 am
Correction of dielectric shielding for scissor operator
Quxiao in BIT,calculate the exciton
 Daniele Varsano
 Posts: 3936
 Joined: Tue Mar 17, 2009 2:23 pm
 Contact:
Re: Correction of dielectric shielding for scissor operator
Dear Quxiao,
the scissor operator is applied by setting KfnQP_E and XfnQP_E in the BSE diagonal part or in the screening respectively.
The scissor operator is applied by shifting the all the conduction bands of the given quantity with respect to the maximum valence band.
When applied to W, the independent particle response Xo is calculated using shifted energies and the Dyson equation is solved: this will give a smaller screening with respect to the one calculated with Kohn Sham energies.
Anyway, the usual procedure is to apply the scissor operator to the diagonal part of the kernel to mimic quasiparticle correction and leave the screening calculated with KS energies as it somehow mimic excitation energies, and usually gives more accurate results, unless you want to apply some sort of eigenvalue selfconsistency in GW calculation to both G and W.
Best,
Daniele
the scissor operator is applied by setting KfnQP_E and XfnQP_E in the BSE diagonal part or in the screening respectively.
The scissor operator is applied by shifting the all the conduction bands of the given quantity with respect to the maximum valence band.
When applied to W, the independent particle response Xo is calculated using shifted energies and the Dyson equation is solved: this will give a smaller screening with respect to the one calculated with Kohn Sham energies.
Anyway, the usual procedure is to apply the scissor operator to the diagonal part of the kernel to mimic quasiparticle correction and leave the screening calculated with KS energies as it somehow mimic excitation energies, and usually gives more accurate results, unless you want to apply some sort of eigenvalue selfconsistency in GW calculation to both G and W.
Best,
Daniele
Dr. Daniele Varsano
S3CNR Institute of Nanoscience and MaX Center, Italy
MaX  Materials design at the Exascale
http://www.nano.cnr.it
http://www.maxcentre.eu/
S3CNR Institute of Nanoscience and MaX Center, Italy
MaX  Materials design at the Exascale
http://www.nano.cnr.it
http://www.maxcentre.eu/

 Posts: 99
 Joined: Fri Mar 26, 2021 11:27 am
Re: Correction of dielectric shielding for scissor operator
Dear Daniele:
Thank you very much for your reply. I have a better understanding of the operation process of Yambo, but I am still very interested in the specific physical form of the scissor operator. How to construct the scissor operator using KfnQP_S and XfnQP_S, and what is the specific physical form, such as scissor operator=f (KfnQP_S, XfnQP_S), what is the specific form of f? This will facilitate my deeper understanding of physical images.
Thanks again!
Best wishes!
Quxiao
Thank you very much for your reply. I have a better understanding of the operation process of Yambo, but I am still very interested in the specific physical form of the scissor operator. How to construct the scissor operator using KfnQP_S and XfnQP_S, and what is the specific physical form, such as scissor operator=f (KfnQP_S, XfnQP_S), what is the specific form of f? This will facilitate my deeper understanding of physical images.
Thanks again!
Best wishes!
Quxiao
Quxiao in BIT,calculate the exciton
 Daniele Varsano
 Posts: 3936
 Joined: Tue Mar 17, 2009 2:23 pm
 Contact:
Re: Correction of dielectric shielding for scissor operator
Dear Quxiao,
please note, the name of the variable is KfnQP_E and not KfnQP_S.
For example:
The form of the operator is a constant scissor (open of the gap) given by the first number (here 2.0eV) and a linear stretching of the conduction and valence band given by the slope in input (second and third entry).
The same for the XfnQP_E but applied to the response function X.
Best,
Daniele
please note, the name of the variable is KfnQP_E and not KfnQP_S.
For example:
Code: Select all
% KfnQP_E
2.000000  1.100000  1.100000  # E parameters (c/v) eVadimadim
%
The same for the XfnQP_E but applied to the response function X.
Best,
Daniele
Dr. Daniele Varsano
S3CNR Institute of Nanoscience and MaX Center, Italy
MaX  Materials design at the Exascale
http://www.nano.cnr.it
http://www.maxcentre.eu/
S3CNR Institute of Nanoscience and MaX Center, Italy
MaX  Materials design at the Exascale
http://www.nano.cnr.it
http://www.maxcentre.eu/

 Posts: 99
 Joined: Fri Mar 26, 2021 11:27 am
Re: Correction of dielectric shielding for scissor operator
Dear Daniele:
I'm very sorry, I made a mistake by typing the wrong variable. I know for sure that it changed KfnQP_E and XfnQP_E,But I would rather know how XfnQP_E affects the response function X and ultimately affects the W matrix, which is more urgent for me to clarify. If you have a clearer physical formula, please help me understand. Thank you!
Best Wishes!
Quxiao
I'm very sorry, I made a mistake by typing the wrong variable. I know for sure that it changed KfnQP_E and XfnQP_E,But I would rather know how XfnQP_E affects the response function X and ultimately affects the W matrix, which is more urgent for me to clarify. If you have a clearer physical formula, please help me understand. Thank you!
Best Wishes!
Quxiao
Quxiao in BIT,calculate the exciton
 Daniele Varsano
 Posts: 3936
 Joined: Tue Mar 17, 2009 2:23 pm
 Contact:
Re: Correction of dielectric shielding for scissor operator
Dear Quixiao,
you can have a look at the calculation of the independent particle polarizability:
see. Eq. 3 of the Yambo cheatsheet: https://www.yambocode.eu/wiki/images/c ... et5.0.pdf
XfnQP_E applies a scissor to the energy difference entering in Xo, so that you will have Ec  Ev+Delta.
The effect is that the larger the scissor, the larger the gap and consequently the screening will be lowered and W enhanced.
Best,
Daniele
you can have a look at the calculation of the independent particle polarizability:
see. Eq. 3 of the Yambo cheatsheet: https://www.yambocode.eu/wiki/images/c ... et5.0.pdf
XfnQP_E applies a scissor to the energy difference entering in Xo, so that you will have Ec  Ev+Delta.
The effect is that the larger the scissor, the larger the gap and consequently the screening will be lowered and W enhanced.
Best,
Daniele
Dr. Daniele Varsano
S3CNR Institute of Nanoscience and MaX Center, Italy
MaX  Materials design at the Exascale
http://www.nano.cnr.it
http://www.maxcentre.eu/
S3CNR Institute of Nanoscience and MaX Center, Italy
MaX  Materials design at the Exascale
http://www.nano.cnr.it
http://www.maxcentre.eu/