Yambo k points from quantum espresso

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
jmullen
Posts: 29
Joined: Wed Apr 01, 2009 6:29 pm

Yambo k points from quantum espresso

Post by jmullen » Thu Oct 22, 2009 7:20 pm

Greetings all,

Just an hour ago, I computed the GW corrections for a graphene sheet and noticed that the kpoints generated in the NSCF Quantum Espresso output file and the QP's in Yambo are not the same. Specifically, the Dirac point is not the same XYZ as in the NSCF file. The index is the same in both output files. Under what conditions does Yambo alter the kpoints instead of evaluating/importing the exact kpoints as in the espresso output?

Regards
Jeff Mullen
NCSU Physics

jmullen
Posts: 29
Joined: Wed Apr 01, 2009 6:29 pm

Re: Yambo k points from quantum espresso

Post by jmullen » Thu Oct 22, 2009 7:44 pm

On closer look, it seems to be just a different coordinate system.
Jeff Mullen
NCSU Physics

User avatar
andrea marini
Posts: 325
Joined: Mon Mar 16, 2009 4:27 pm
Contact:

Re: Yambo k points from quantum espresso

Post by andrea marini » Thu Oct 22, 2009 8:30 pm

jmullen wrote:Under what conditions does Yambo alter the kpoints instead of evaluating/importing the exact kpoints as in the espresso output?
Dear Jeff, of course yambo does not alter the PW k-points list ;) As you said it is only a matter of reference system. In yambo k-points coordinates are written in IKU=Internal K Units. Check any report file and search for

Code: Select all

 Unit cell scaling factors [a.u.]:  4.71637   4.08450  12.17584


The IKU is simply obtained by applying the tranformation

(K1,K2,K3)_IKU=(Kx*2*pi/4.72,Kx*2*pi/4.08,Ky*2*pi/12.18)

It is simply a generalization of the 2*pi/alat used by PW. Note that the unit cell scaling factors are completely arbitrtary as they are internal to the code. They are simply used to have a more user friendly notation for the vectors in the Reciprocal lattice. So these factors can differ from the PW lattice constant. This is way, even in a cubic system, PW and yambo can print appraently different k-points.
Andrea MARINI
Istituto di Struttura della Materia, CNR, (Italy)

Post Reply