Page 1 of 1
Coulomn Interaction
Posted: Thu Mar 22, 2012 9:47 am
by WEI Wei
Dear All,
I am going to do a GW calculation about a 1D ribbon structure along y direction, how can I set the variable of "CutGeo"? "CutGeo=none" or "CutGeo=box xz" or "CutGeo=box"?
Many thanks in advance!!!
Re: Coulomn Interaction
Posted: Thu Mar 22, 2012 10:11 am
by Daniele Varsano
Dea Wei,
Find here an example of input file.
Kepp in mind that the given Lx, the interacting box goes
from -Lx/2 to Lx/2, so Lx/2 should be large enough to include
all your system and the size of supercell big enough that two replica
do not interact in a distance of Lx/2. (see the picture below).
Code: Select all
xxxxxx )( xxxxxxxxx )
<--------|--------->
-Lx/2 Lx/2
Moreover in order to build up the box Coulomb potential, you need to use the RIM technique
(Random integration method (see the Yambo paper (Comp. Phys. Comm. 180, 1392 (2009)),
as integral of Coulomb potential over the Brillouin zones are needed (the number in this input files
are totally indicative, and may be you need much less RandGvec).
Code: Select all
RandQpts= 1000000 # [RIM] Number of random q-points in the BZ
RandGvec= 123 RL # [RIM] Coulomb interaction RS components
CUTGeo= "box xz" # [CUT] Coulomb Cutoff geometry: box/cylinder/sphere
% CUTBox
Lx| 0.0000 | Lz | # [CUT] [au] Box sides
%
Cheers,
Daniele
Re: Coulomn Interaction
Posted: Thu Mar 22, 2012 10:18 am
by WEI Wei
Dear Daniele,
Great! You always help me a lot!! Thank you so much!!!