Energy band values with p2y/a2y

Come here to discuss issues with using the p2y converter from PWscf to Yambo, how to best run PWscf , link with the IOTK library, and so on.

Moderators: andrea.ferretti, Conor Hogan

Forum rules
Before posting, make sure to first read carefully all the detailed steps and suggestions given in the main documentation:
https://www.yambo-code.eu/learn/#learn-forum
Post Reply
Laura Caputo
Posts: 88
Joined: Sun Apr 11, 2021 3:02 pm

Energy band values with p2y/a2y

Post by Laura Caputo » Fri Jun 04, 2021 4:44 pm

Good evening,

I am doing GW band gap calculations using Quantum Espresso. I noticed that the energies of the bands that I consider for calculating the GW gap are different when using Yambo+Abinit and Yambo+Quantum Espresso.
Before moving to Yambo, I checked the value of the gap at the DFT level which was the same for both Abinit and QE. When moving to Yambo, I also checked the gap at the DFT level which again is the same. So, the only thing that changes are the band energies. Is this due to some differences in p2y and a2y? Or withing the DFT codes?

Sincerely,
Laura
Laura Caputo
Ph.D. Student
Université Catholique de Louvain
https://uclouvain.be/fr/repertoires/laura.caputo

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

Re: Energy band values with p2y/a2y

Post by Daniele Varsano » Fri Jun 04, 2021 4:55 pm

Dear Laura,
So, the only thing that changes are the band energies. Is this due to some differences in p2y and a2y? Or withing the DFT codes?
What do you mean exactly by band energies?
p2y and a2y are just meant to read the ground state electronic structure for both the codes and rewrite in yambo format. Of course, if the overall electronic structure is different, apart from the gap, you can find different results. Anyway it seems to me quite odd that you find different electronic structure with the two codes if you use the same set up, e.g. pseudo potentials and functionals.

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/

andrea.ferretti
Posts: 206
Joined: Fri Jan 31, 2014 11:13 am

Re: Energy band values with p2y/a2y

Post by andrea.ferretti » Fri Jun 04, 2021 4:57 pm

Dear Laura,
Is (1) a rigid shift of the bands compatible with what you observe ?
or is (2) the shape of the bands actually modified ?

(1) is expected, since yambo shifts the DFT eigenvalues such that the KS-HOMO (or the Fermi level) is set to zero.
(2) instead would be a bug.

take care
Andrea
Andrea Ferretti, PhD
CNR-NANO-S3 and MaX Centre
via Campi 213/A, 41125, Modena, Italy
Tel: +39 059 2055322; Skype: andrea_ferretti
URL: http://www.nano.cnr.it

Laura Caputo
Posts: 88
Joined: Sun Apr 11, 2021 3:02 pm

Re: Energy band values with p2y/a2y

Post by Laura Caputo » Fri Jun 04, 2021 5:00 pm

Dear Daniele,

I mean like this. With Quantum Espresso:

# K-point Band Eo [eV]
#
1 36 -0.754478
1 37 0.790687

With Abinit:
1 36 0.000000
1 37 1.551236

The gap is quite the same (I used a higher cut off in one of the two codes, that's why one is only slightly higher). However, the values of the energies of the bands are different, as a shift has happened.
Last edited by Laura Caputo on Fri Jun 04, 2021 5:23 pm, edited 1 time in total.
Laura Caputo
Ph.D. Student
Université Catholique de Louvain
https://uclouvain.be/fr/repertoires/laura.caputo

Laura Caputo
Posts: 88
Joined: Sun Apr 11, 2021 3:02 pm

Re: Energy band values with p2y/a2y

Post by Laura Caputo » Fri Jun 04, 2021 5:02 pm

Dear Andrea,

Sorry maybe I was answering in the meantime that you were replying. As I said, yes it is like a shift and the HOMO value is set to zero only within the Abinit framework.

Regarding the shape of the bands, you mean by checking the Density of States/band structure in DFT framework?
You do not have the required permissions to view the files attached to this post.
Laura Caputo
Ph.D. Student
Université Catholique de Louvain
https://uclouvain.be/fr/repertoires/laura.caputo

andrea.ferretti
Posts: 206
Joined: Fri Jan 31, 2014 11:13 am

Re: Energy band values with p2y/a2y

Post by andrea.ferretti » Fri Jun 04, 2021 5:09 pm

Hi Laura,

no worries then,
In the case of QE (at least) yambo applies a rigid shift as discussed above.
It could well be that in the case of Abinit this does not happen, or happens differently.

A
Andrea Ferretti, PhD
CNR-NANO-S3 and MaX Centre
via Campi 213/A, 41125, Modena, Italy
Tel: +39 059 2055322; Skype: andrea_ferretti
URL: http://www.nano.cnr.it

Laura Caputo
Posts: 88
Joined: Sun Apr 11, 2021 3:02 pm

Re: Energy band values with p2y/a2y

Post by Laura Caputo » Fri Jun 04, 2021 5:10 pm

Dear Andrea,

In my case, it seems like it is the inverse (QE does not apply it and Abinit does). But however, as you said, it is a rigid shift so, if I understand, it should be fine.

Thanks for your support
Laura Caputo
Ph.D. Student
Université Catholique de Louvain
https://uclouvain.be/fr/repertoires/laura.caputo

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

Re: Energy band values with p2y/a2y

Post by Daniele Varsano » Fri Jun 04, 2021 9:52 pm

Dear Laura,
the VBM in Yambo is set to zero.
Are you using the same k point grid (number and shifts) in the two codes?
Is it possible that the k-point index 1 does not correspond to the same point?
you can check in the report the corresponding kpoint and also see at what k point the direct gap is found.

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