ip rpa eps and eps0

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
zhukov_vladlen
Posts: 1
Joined: Wed Aug 06, 2014 12:39 pm

ip rpa eps and eps0

Post by zhukov_vladlen » Fri Aug 15, 2014 7:05 am

Dear all,
I have just compilated yambo 3.4.1 and done IP and RPA calculations for optics in metallic titanium and rutile. The commands used are usual:
yambo -o c -k hartree,
with Chimod = IP or Hartree in the yambo.in.

What is strange I find in the calculated o.eps...... files of RPA: here the eps-results are almost identical to eps0, with differences no more that 0.00001. I find this for both Re and Im parts, for all the Q-vectors. Moreover, the results of IP-calculos are almost identical to those of RPA. Is this true?

Vladlen Zhukov.

User avatar
Daniele Varsano
Posts: 3829
Joined: Tue Mar 17, 2009 2:23 pm
Contact:

Re: ip rpa eps and eps0

Post by Daniele Varsano » Fri Aug 15, 2014 8:32 am

Dear Vladlen,
please fill your signature with your affiliation, this is a rule of the forum.
Indeed your results are quite strange, in order to look what is going wrong is always useful to post the output and report files.
In any case I wonder if you raised the variable NGsBlk until convergences, this variable controls the G vector block of the response function taking account of the local field effects. If you leave the NGsBlk=1 no local field are included and the Hartree approximation is equivalent to the IP approximations.
In the tutorial section (fantastic dimension) there are some examples on RPA calculations:

http://www.yambo-code.org/theory/docs/doc_Xd.php#RPA
http://www.yambo-code.org/tutorials/fan ... ns/LiF.php

Best,

Daniele
Dr. Daniele Varsano
S3-CNR Institute of Nanoscience and MaX Center, Italy
MaX - Materials design at the Exascale
http://www.nano.cnr.it
http://www.max-centre.eu/

Post Reply