I have done the real time simulation by solving the time dependent BSE as implemented in Yambo_rt. By applying a small field and using the DELTA field I could obtained a optical spectrum which looks like the spectrum obtained by solving the GW+BSE method. I want to applied a large intensity field (non-linear case) and see the effect of it in the absorption spectrum. The field related parameter which is used is given as follows.
Code: Select all
% Field1_Freq
2.230000 | 2.230000 | eV # [RT Field1] Frequency
%
Field1_DFreq= 0.000000 eV # [RT Field1] Frequency step
Field1_Int=100000000.000000 kWLm2 # [RT Field1] Intensity
Field1_Width=0.000000 fs # [RT Field1] Width
Field1_FWHM=100.000000 fs # [RT Field1] Full Width at Half Maximum (overwrite width if set)
Field1_kind= "QSSIN" # [RT Field1] Kind(SIN|RES|ANTIRES|GAUSS|DELTA|QSSIN)
Field1_pol= "linear" # [RT Field1] Pol(linear|circular)
% Field1_Dir
0.000000 | 1.000000 | 0.000000 | # [RT Field1] Versor
%
% Field1_Dir_circ
0.000000 | 1.000000 | 0.000000 | # [RT Field1] Versor_circ
Code: Select all
# Time[fs] DirGap [meV] IndGap [meV] DeltaDirGap [meV] DeltaIndGap [meV]
#
0.00000000 0.00000000 0.00000000 0.00000000 0.00000000
0.100000001 0.00000000 0.00000000 0.00000000 0.00000000
0.200000003 0.00000000 0.00000000 0.00000000 0.00000000
0.300000012 0.00000000 0.00000000 0.00000000 0.00000000
i) Is this file(o.BGR) related to band gap renormalization?
ii) At each time step why the gap(direct or indirect) value is not updated and gives zero value for all time start to end. What went wrong?
Regards,