Dear Yamboo developer,
I want to write the independent-particle response function at G=G'=0, i.e., \chi_0(q, G=G'=0), and the oscillator strength rho defined as rho(q)=<k1,b1|e^{- iqr}| k2, b2>. Could you please how can I find them? I am a new Yambo user. I found in routine "vi src/wf_and_fft/scatter_Bamp.F" rho is difined as "rhotw(G)=<ib ik_is i_sp| e^[I(RsG-G0).r] | ob ok_os o_sp>", thus I would like to print ib, ik_is, , ob, ok_os, rhotw(G=0).
Thank you very much in advance!
Cheers,
Sky
writing \chi_0(q, G=G'=0) and rho(q, G=0)
Moderators: Davide Sangalli, andrea.ferretti, myrta gruning, andrea marini, Daniele Varsano, Conor Hogan
-
- Posts: 8
- Joined: Mon Feb 19, 2018 12:32 pm
writing \chi_0(q, G=G'=0) and rho(q, G=0)
Jianqiang Sky ZHOU
Institut des NanoSciences de Paris (INSP)
European Theoretical Spectroscopy Facility
Institut des NanoSciences de Paris (INSP)
European Theoretical Spectroscopy Facility
- claudio
- Posts: 526
- Joined: Tue Mar 31, 2009 11:33 pm
- Location: Marseille
- Contact:
Re: writing \chi_0(q, G=G'=0) and rho(q, G=0)
Dear Sky
dipole matrix elements are calculated in
src/pol_function/DIPOLE_driver.F
while the dielectric constant is constructed in
src/pol_function/X_irredux.F
best
claudio
dipole matrix elements are calculated in
src/pol_function/DIPOLE_driver.F
while the dielectric constant is constructed in
src/pol_function/X_irredux.F
best
claudio
Claudio Attaccalite
[CNRS/ Aix-Marseille Université/ CINaM laborarory / TSN department
Campus de Luminy – Case 913
13288 MARSEILLE Cedex 09
web site: http://www.attaccalite.com
[CNRS/ Aix-Marseille Université/ CINaM laborarory / TSN department
Campus de Luminy – Case 913
13288 MARSEILLE Cedex 09
web site: http://www.attaccalite.com