Skip to Content.
Sympa Menu

forum - Re: Re: [abinit-forum] a bug in GW code (cchi0q0.F90)?

forum@abinit.org

Subject: The ABINIT Users Mailing List ( CLOSED )

List archive

Re: Re: [abinit-forum] a bug in GW code (cchi0q0.F90)?


Chronological Thread 
  • From: deyulu@yahoo.com
  • To: forum@abinit.org
  • Subject: Re: Re: [abinit-forum] a bug in GW code (cchi0q0.F90)?
  • Date: Fri, 23 Mar 2007 18:48:42 +0100

Fabien:
Thank you for your explanation, but I don't agree with your argument.
In doing scissor+G0W0, self-consistent GW on the energies-only, one should
use
H', i.e., H_{KS}+scissor or H_{KS}+Delta sigma, to evaluate quasi-particle
energies.
It applies also to the commutator trick. If we takes the 1st order app.
on H' to assume it is diagonal, one has:
< v | r | c > = < v | [ H' , r ] | c > / (E'_v - E'_c),
where E'_v and E'_c are quasi-particle energies instead of KS energies.
For a wide range of materials, it has been shown in Fiorentitni and
Balderesch's paper (PRB, 51:17916, 1995) that the effect of Delta sigma can
be largely attributed to a static screened exchange term, which commutes with
r. Then it leads to < v | r | c > = < v | [ H_{KS} , r ] | c > / (E'_v -
E'_c).
In my opinion, this expression is consistent with the formula for G!=0
terms.

Best

Deyu***************************************************************************
Deyu Lu (Ph.D)

190 Chemistry Building
University of California, Davis
One Shields Avenue
Davis, CA 95616
Office phone: (530) 754-9663
Group Webpage: http://angstrom.ucdavis.edu/

***************************************************************************



Archive powered by MHonArc 2.6.16.

Top of Page