I am new to quantum espresso and yambo, but has some experience with abinit and DP codes. Recently, I find something that is wired to me:
I start from LDA+U, giving the k point grid as 4 4 4 0 0 0 (Gamma centered 4x4x4 grid) in QE that should give 64 k points in the full BZ (12 k points in the Irreducible part BZ), however, yambo finds only 58 k points in the FBZ. When I changed to LDA calculation of Si, giving the same k grid, yambo finds the correct number 64 for the FBZ....
I attached " r_setup" file from the LDA+U calculation to make my question clear.
Can someone tell me the reason? Maybe it is better to tell me how does yambo find the kpoints in FBZ from the k points in IBZ using the symmetry, e.g., tell me the routine where this is done.
Many thanks!
Cheers,
Sky
the total number of k point in the full BZ
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:
http://www.yambocode.org/doc/databases/index.php
Before posting, make sure to first read carefully all the detailed steps and suggestions given in the main documentation:
http://www.yambocode.org/doc/databases/index.php

 Posts: 8
 Joined: Mon Feb 19, 2018 12:32 pm
the total number of k point in the full BZ
Jianqiang Sky ZHOU
Institut des NanoSciences de Paris (INSP)
European Theoretical Spectroscopy Facility
Institut des NanoSciences de Paris (INSP)
European Theoretical Spectroscopy Facility
 Daniele Varsano
 Posts: 2326
 Joined: Tue Mar 17, 2009 2:23 pm
 Contact:
Re: the total number of k point in the full BZ
Dear Sky,
please note the r_setup is not attached.
yambo at the moment does not support DFT+U calculation as starting point: we are working on that but at the moment it is not supported
yambo reads the K points in the IBZ and symmetries from QE, from that reconstruct the points in the BZ. Please note that nonsymmorphic symmetries are not supported.
You can compare if all the symmetries of PW are correctly reported in the yambo report.
Best,
Daniele
please note the r_setup is not attached.
yambo at the moment does not support DFT+U calculation as starting point: we are working on that but at the moment it is not supported
yambo reads the K points in the IBZ and symmetries from QE, from that reconstruct the points in the BZ. Please note that nonsymmorphic symmetries are not supported.
You can compare if all the symmetries of PW are correctly reported in the yambo report.
Best,
Daniele
Dr. Daniele Varsano
S3CNR Institute of Nanoscience and MaX Center, Italy
MaX  Materials design at the Exascale
http://www.nano.cnr.it
http://www.maxcentre.eu/
S3CNR Institute of Nanoscience and MaX Center, Italy
MaX  Materials design at the Exascale
http://www.nano.cnr.it
http://www.maxcentre.eu/

 Posts: 8
 Joined: Mon Feb 19, 2018 12:32 pm
Re: the total number of k point in the full BZ
Dear Daniele,
1. I am sorry for the missing file, but now I realize that the problem is in QE instead of yambo... So it does not matter anymore...
2. I understand your concern about DFT+U but as far as I know, it only matters at q=0 limit because yamboo needs the commutator to the Hamiltonian which contains U. All the rest q values should be safe, right? I hope you can confirm that so I will be able to calculate some chi at q!= 0 .
3. The problem now is that the symmetries report by QE is not complete.... I am working on that these days, and find that QE report 12 symmetries but what it actually uses to get the inequivalent k poitns are 23 symmetries....maybe because of +k and k problem.... I do not know the reason yet but trying to figure out this point.
4. I put force_symmorphic in QE calculation so that the nonsymmorphic symmetries are not taken into account in QE.
Many thanks!
Sky
1. I am sorry for the missing file, but now I realize that the problem is in QE instead of yambo... So it does not matter anymore...
2. I understand your concern about DFT+U but as far as I know, it only matters at q=0 limit because yamboo needs the commutator to the Hamiltonian which contains U. All the rest q values should be safe, right? I hope you can confirm that so I will be able to calculate some chi at q!= 0 .
3. The problem now is that the symmetries report by QE is not complete.... I am working on that these days, and find that QE report 12 symmetries but what it actually uses to get the inequivalent k poitns are 23 symmetries....maybe because of +k and k problem.... I do not know the reason yet but trying to figure out this point.
4. I put force_symmorphic in QE calculation so that the nonsymmorphic symmetries are not taken into account in QE.
Many thanks!
Sky
Daniele Varsano wrote:Dear Sky,
please note the r_setup is not attached.
yambo at the moment does not support DFT+U calculation as starting point: we are working on that but at the moment it is not supported
yambo reads the K points in the IBZ and symmetries from QE, from that reconstruct the points in the BZ. Please note that nonsymmorphic symmetries are not supported.
You can compare if all the symmetries of PW are correctly reported in the yambo report.
Best,
Daniele
Jianqiang Sky ZHOU
Institut des NanoSciences de Paris (INSP)
European Theoretical Spectroscopy Facility
Institut des NanoSciences de Paris (INSP)
European Theoretical Spectroscopy Facility
 Daniele Varsano
 Posts: 2326
 Joined: Tue Mar 17, 2009 2:23 pm
 Contact:
Re: the total number of k point in the full BZ
Dear Sky,
It depends on what you need to calculate. If you are interested in linear response it should be sage, but if you need to do a GW calculation then you need to subctract from the self energy the U contribution to each KS state and this part needs coding we are planning to do in the next days/months.
About symmetry problem you can check in the r_setup if yambo recognize that there is the time reversal symmetry.
Best,
Daniele
Code: Select all
2. I understand your concern about DFT+U but as far as I know, it only matters at q=0 limit because yamboo needs the commutator to the Hamiltonian which contains U. All the rest q values should be safe, right? I hope you can confirm that so I will be able to calculate some chi at q!= 0 .
About symmetry problem you can check in the r_setup if yambo recognize that there is the time reversal symmetry.
Best,
Daniele
Dr. Daniele Varsano
S3CNR Institute of Nanoscience and MaX Center, Italy
MaX  Materials design at the Exascale
http://www.nano.cnr.it
http://www.maxcentre.eu/
S3CNR Institute of Nanoscience and MaX Center, Italy
MaX  Materials design at the Exascale
http://www.nano.cnr.it
http://www.maxcentre.eu/

 Posts: 8
 Joined: Mon Feb 19, 2018 12:32 pm
Re: the total number of k point in the full BZ
Dear Daniele,
Thank you very much for your explanation. At this moment I think I will only do a linear response RPA calculation. If I want to do GW, I will start from LDA but not LDA+U...
I have checked carefully the time reversal symmetry in QE code, it does not make any difference on this problem.... Actually, I will use shifted grid or double grid such that symmetry will not be used, so at this moment, it is fine to run calculations. But I will also try to figure out what is happening when a nonshifted grid is used.
Cheers,
Sky
Thank you very much for your explanation. At this moment I think I will only do a linear response RPA calculation. If I want to do GW, I will start from LDA but not LDA+U...
I have checked carefully the time reversal symmetry in QE code, it does not make any difference on this problem.... Actually, I will use shifted grid or double grid such that symmetry will not be used, so at this moment, it is fine to run calculations. But I will also try to figure out what is happening when a nonshifted grid is used.
Cheers,
Sky
Daniele Varsano wrote:Dear Sky,
It depends on what you need to calculate. If you are interested in linear response it should be sage, but if you need to do a GW calculation then you need to subctract from the self energy the U contribution to each KS state and this part needs coding we are planning to do in the next days/months.Code: Select all
2. I understand your concern about DFT+U but as far as I know, it only matters at q=0 limit because yamboo needs the commutator to the Hamiltonian which contains U. All the rest q values should be safe, right? I hope you can confirm that so I will be able to calculate some chi at q!= 0 .
About symmetry problem you can check in the r_setup if yambo recognize that there is the time reversal symmetry.
Best,
Daniele
Jianqiang Sky ZHOU
Institut des NanoSciences de Paris (INSP)
European Theoretical Spectroscopy Facility
Institut des NanoSciences de Paris (INSP)
European Theoretical Spectroscopy Facility