Energy gap from optical absorption

[DavidPolito93]

Dear developers,

I am treating a 1D system oriented along the z direction with the ws cutoff. I have computed the band structure within the G0W0 approximation and I obtain an energy gap of ~ 2 eV.
Then I computed the optical spectra (without BSE): Yambo gives me the alpha file (by the way, why does it have four columns? I would have expected only two, one for the imaginary part and one for the real part. Instead I have two imaginary parts and two real parts).

My problem is that I would have expected a peak around 2 eV for the imaginary part, but I get a peak at exactly half of the value of the energy gap (i.e. ~ 1 eV).
Is there some factor of 2 that have to be taken into account? Or am I doing something wrong?

I carefully checked the convergence for all of the parameters. Maybe I missed something?

I am attaching the input:

optics # [R OPT] Optics
rim_cut # [R RIM CUT] Coulomb potential
chi # [R CHI] Dyson equation for Chi.
dipoles # [R ] Compute the dipoles
DIP_Threads=0 # [OPENMP/X] Number of threads for dipoles
X_Threads=0 # [OPENMP/X] Number of threads for response functions
RandQpts=0 # [RIM] Number of random q-points in the BZ
RandGvec= 1 RL # [RIM] Coulomb interaction RS components
CUTGeo= "ws z" # [CUT] Coulomb Cutoff geometry: box/cylinder/sphere/ws X/Y/Z/XY..
CUTwsGvec= 1.100000 # [CUT] WS cutoff: number of G to be modified
Chimod= "HARTREE" # [X] IP/Hartree/ALDA/LRC/PF/BSfxc
NGsBlkXd= 3 Ry # [Xd] Response block size
% QpntsRXd
1 | 1 | # [Xd] Transferred momenta
%
% BndsRnXd
1 | 30 | # [Xd] Polarization function bands
%
% EnRngeXd
0.00000 | 40.00000 | eV # [Xd] Energy range
%
% DmRngeXd
0.20000 | 0.20000 | eV # [Xd] Damping range
%
ETStpsXd= 4001 # [Xd] Total Energy steps
% LongDrXd
0.000000 | 0.000000 | 1.000000 | # [Xd] [cc] Electric Field
%

Sincerely,

Davide Romanin
-----------------------------------------------------
PhD student in Physics XXXIII cycle
Representative of the PhD students in Physics
Applied Science and Technology department (DiSAT)
Politecnico di Torino
Corso Duca degli Abruzzi, 24
10129 Torino ITALY
------------------------------------------------------

[Palummo Maurizia]

Dear Davide

you are obtaining 4 columns because you are doing an RPA + local-field calculation. the column 4 and 5 are the RPA without LF (real IP calculation) while 2 and 3
are with LF included.
Regarding the peak at 1 eV are you sure you are reading the GW corrections?

Best
Maurizia

[DavidPolito93]

Dear Maurizia,

Thank you for your rapid reply!

Effectively, I checked and it is the PBE gap (the one without the GW corrections). I thought I gave the correct SAVE folder. I will try to redo the computation and I'll let you know

Thanks a lot,

Davide
-----------------------------------------------------
PhD student in Physics XXXIII cycle
Representative of the PhD students in Physics
Applied Science and Technology department (DiSAT)
Politecnico di Torino
Corso Duca degli Abruzzi, 24
10129 Torino ITALY
------------------------------------------------------

[Palummo Maurizia]

Dear Davide

in any case you have to read the ndb.QP file or add the correction as scissor + streaching in the input file
if you want to do a calculation of epsilon2 at RPA+LF starting from the corrected bands.

Maurizia

[Daniele Varsano]

Dear Davide,
in order to include your GW quasiparticle you need to add the following keywords in the input file:

Code: Select all

XfnQPdb= "E < ./dir/ndb.QP"

where dir is the directory containing the QP database (SAVE if you have not redirected in other directory).

Alternatively you can assign a scissor/streching correction by using:

Code: Select all

% XfnQP_E
 delta-E | slople_conduction | slope_valence |        # [EXTQP Xd] E parameters  (c/v) eV|adim|adim
%

These variables are activated adding the qp verbosity un the input file generation:

yambo -o c -V qp

Best,
Daniele