Any way to transform the GW result by Abinit to db.QP?

Concerns issues with computing quasiparticle corrections to the DFT eigenvalues - i.e., the self-energy within the GW approximation (-g n), or considering the Hartree-Fock exchange only (-x)

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

Post Reply
Guangfu Luo
Posts: 17
Joined: Mon Dec 14, 2009 10:46 am

Any way to transform the GW result by Abinit to db.QP?

Post by Guangfu Luo » Tue Mar 08, 2011 8:22 am

Dear friends,

To do a BSE calculation, one can input the GW results either using the approximate scissor-and-stretching method (with "KfnQP_E") or directly importing the original results db.QP (with "KfnQPdb"). Sometimes, the GW corrections are complex and using db.QP is the only good choice.

The problem is that if I do the GW calculation in Abinit, how can I transform the Abinit results to the binary db.QP? What I hope to know is the file structure of db.QP. Can you enlighten me?

Sincerely,

PS. If KfnQPdb can also accept ASCII file, that whould be great.
Guangfu Luo
Institute for Molecular Science
National Institutes of Natural Sciences, Japan

User avatar
claudio
Posts: 458
Joined: Tue Mar 31, 2009 11:33 pm
Location: Marseille
Contact:

Re: Any way to transform the GW result by Abinit to db.QP?

Post by claudio » Tue Mar 08, 2011 2:44 pm

Dear Guangfu Luo

it is not possible to import directly the ABINIT results in Yambo.

But as fas as I know other people modified yambo in order to do it.

Try to contact Silvana Botti

http://www.tddft.org/bmg/people.php?oid=5040305

Sincerely
Claudio Attaccalite
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

Guangfu Luo
Posts: 17
Joined: Mon Dec 14, 2009 10:46 am

Re: Any way to transform the GW result by Abinit to db.QP?

Post by Guangfu Luo » Wed Mar 09, 2011 4:10 am

Dear Claudio,

Thanks very much for your infomation. I will contact Silvana Botti.

Sincerely,
Guangfu Luo
Institute for Molecular Science
National Institutes of Natural Sciences, Japan

User avatar
claudio
Posts: 458
Joined: Tue Mar 31, 2009 11:33 pm
Location: Marseille
Contact:

Re: Any way to transform the GW result by Abinit to db.QP?

Post by claudio » Mon May 23, 2011 4:46 pm

Dear Guangfu Luo

in the last version of yambo 3,2,4,

http://www.yambo-code.org/download.php

we added a tool to write quasi-particle correction in the yambo database, ypp -q

regards
Claudio Attaccalite
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

gfluo
Posts: 8
Joined: Sun May 15, 2011 5:10 pm

Re: Any way to transform the GW result by Abinit to db.QP?

Post by gfluo » Fri Jun 17, 2011 7:50 am

Dear Claudio,

This is good news and thanks for your work. I execute "ypp -q" but found some part unclear for me. Specifically, I have two questions: (1) what is meaning of "Img(E)"? (2) Can I use more than one line to specify the QP correction one state by one state?

---------------parameters generated by "ypp -q"---------------------

qpdb # [R] Quasi-particle database manager
%QPkrange # [QPDB]QPDB energy correction( kp1| kp2| bnd1| bnd2| E-Eo[eV]| Img(E)[eV]| Z|)
-1| -1| -1| -1| 0.0| 0.0| 0.0|
%
---------------------------------------------------------------------------

Sincerely,
Guangfu Luo
Institute for Molecular Science
National Institutes of Natural Sciences
38 Nishigo-Naka, Myodaiji, Okazaki 444-8585, Japan

User avatar
claudio
Posts: 458
Joined: Tue Mar 31, 2009 11:33 pm
Location: Marseille
Contact:

Re: Any way to transform the GW result by Abinit to db.QP?

Post by claudio » Fri Jun 17, 2011 3:29 pm

Dear Guangfu Luo

ypp -q works like this, image you have two bands to whith
you want apply the QP correction the number 3 and 4 and the kpoints 1 2 3 4

1| 4| 3| 3| -0.5| 0.0| 0.0|
1| 4| 4| 4| +1.0| 0.0| 0.0|

the band 3 will be shifted down by 0.5 eV for the kpoints 1,2,3,4
while the band 4 will be shifted up by 1.0 eV for the kpoints 1,2,3,4

Img(E) is the immaginary part of the Self-Energy in case you have it,
but I guess if you have done normal self-consistent GW, you have only the real-part.
Z is the renomalization factor.

Last thing, please download the yambo version on SVN (revision 18)
because we corrected minor bugs

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

gfluo
Posts: 8
Joined: Sun May 15, 2011 5:10 pm

Re: Any way to transform the GW result by Abinit to db.QP?

Post by gfluo » Mon Jun 20, 2011 2:37 pm

Dear Claudio,

Thanks for your explanation and kind reminder.

Sincerely,
Guangfu Luo
Institute for Molecular Science
National Institutes of Natural Sciences
38 Nishigo-Naka, Myodaiji, Okazaki 444-8585, Japan

Post Reply