writing \chi_0(q, G=G'=0) and rho(q, G=0)

Deals with issues related to computation of optical spectra in reciprocal space: RPA, TDDFT, local field effects.

Moderators: Davide Sangalli, andrea.ferretti, myrta gruning, andrea marini, Daniele Varsano, Conor Hogan

Post Reply
sky-zhou
Posts: 8
Joined: Mon Feb 19, 2018 12:32 pm

writing \chi_0(q, G=G'=0) and rho(q, G=0)

Post by sky-zhou » Thu Nov 22, 2018 3:55 pm

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
Jianqiang Sky ZHOU
Institut des NanoSciences de Paris (INSP)
European Theoretical Spectroscopy Facility

User avatar
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)

Post by claudio » Fri Nov 23, 2018 9:34 am

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
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

Post Reply