Dear developers:
When I read paper PRL 105, 026401 (2010), in it reference 15, there is a setence "G0W0 calculations employ the plasmon pole approximation
and include local-field effects in the dielectric function and self-energy sigma with an energy cutoff of 6 Ry."
I am confused the self-energy with cutoff of 6 Ry. In the input file of G0W0 , where can I set the self-energy?
In general, My procedure is:
>yambo -i
>yambo -x
>yambo -g n -p p
S D Wang
self-energy sigma
Moderators: Davide Sangalli, andrea.ferretti, myrta gruning, andrea marini, Daniele Varsano, Conor Hogan
-
- Posts: 287
- Joined: Fri Apr 09, 2010 12:30 pm
- myrta gruning
- Posts: 240
- Joined: Tue Mar 17, 2009 11:38 am
- Contact:
Re: self-energy sigma
The authors refer to local-field effects. In the dielectric function those are controlled by the response block size, thus in Yambo by NGsBlkXp.sdwang wrote: "G0W0 calculations employ the plasmon pole approximation
and include local-field effects in the dielectric function and self-energy sigma with an energy cutoff of 6 Ry."
I am confused the self-energy with cutoff of 6 Ry. In the input file of G0W0 , where can I set the self-energy?
The dielectric function is then used for W in the self-energy, and I guess this is what the authors of that paper mean when they mention the cutoff for local field effects in the self-energy. In Yambo (yambo -g n -p p) you need just to specify NGsBlkXp (at least this is my interpretation).
Best,
m
Dr Myrta Grüning
School of Mathematics and Physics
Queen's University Belfast - Northern Ireland
http://www.researcherid.com/rid/B-1515-2009
School of Mathematics and Physics
Queen's University Belfast - Northern Ireland
http://www.researcherid.com/rid/B-1515-2009
-
- Posts: 287
- Joined: Fri Apr 09, 2010 12:30 pm
- Daniele Varsano
- Posts: 3816
- Joined: Tue Mar 17, 2009 2:23 pm
- Contact:
Re: self-energy sigma
Hi Shudong,
you can set the cutoff in Rydberg, i.e. you can include all the Gvectors, such that
G^2/2 < E_cut(Ry). In order to do that you can replace the RL unit in the yambo input
with mHa (i.e. milli Hartree). 1Hartree=2Ry
Cheers,
Daniele
you can set the cutoff in Rydberg, i.e. you can include all the Gvectors, such that
G^2/2 < E_cut(Ry). In order to do that you can replace the RL unit in the yambo input
with mHa (i.e. milli Hartree). 1Hartree=2Ry
Cheers,
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/
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/