<div dir="ltr">Hi Pietro,<div><br></div><div>Thanks for your detailed answer, also for the unfold.x code and the recent updates, I've used it several times now and have found your examples very helpful. I'll keep your description in mind when reading the PRB, it often helps (me) to have a rough answer in mind when reading through the details. </div><div><br></div><div>On BandUP, I looked into this some months ago and found that the code could not work with versions of QE past 6.3 (I think - maybe there was some change to the output of pw.x after this version?). I've just looked at the github page again (<a href="https://github.com/band-unfolding/bandup">https://github.com/band-unfolding/bandup</a>) and it has not been updated yet, so still not compatible.</div><div><br></div><div>Thanks for your offer, I may contact you in the future,</div><div><br></div><div>Cheers</div><div><br></div><div><div>Will Hewett</div><div>Post Doctoral Researcher</div><div>Victoria University of Wellington</div><div>New Zealand</div></div><div><br></div><div><br></div></div><br><div class="gmail_quote"><div dir="ltr" class="gmail_attr">On Tue, Jan 5, 2021 at 10:34 AM Pietro Bonfa' <<a href="mailto:pietro.bonfa@unipr.it">pietro.bonfa@unipr.it</a>> wrote:<br></div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex">Dear Will,<br>
<br>
the quantities computed by the code are presented in Phys. Rev. B 85, <br>
085201 (2012).<br>
<br>
More precisely, you can collect the quantity in Eq. 8 (detailed in <br>
Appendix A) for the k points specified in input. Roughly speaking, that <br>
is the amount unfolded states with equivalent k vector contribute to the <br>
folded state at a given energy. Admittedly this is far from clear but <br>
the article mentioned above describes very accurately both the problem <br>
and the proposed solution. Still I hope that my awkward description will <br>
be sufficient to understand that, in a perfect supercell, this number <br>
depends on the degeneracy of the eigenvalue in the base cell. If the <br>
symmetry is low enough (and forgetting about spin), this number is <br>
either 0 or 1.<br>
<br>
The picture that you are sharing is actually Eq. 9, which is the <br>
quantity described above times a Dirac delta function of the energy, <br>
which is approximated with a Gaussian.<br>
<br>
As a consequence, the 'weight' in a perfect supercell depends on the <br>
degeneracy (or almost degeneracy) of a state, the discretization of the <br>
energy interval and the width of the Gaussian approximating the Dirac delta.<br>
I hope this partially answers your question.<br>
<br>
That being said, as you may have read, I've been always advocating for <br>
bandUP, since unfold.x was just my exercise to learn Fortran and the <br>
internals of QE.<br>
However, lately I found a little bit of time to fix some problems, add <br>
some tests and parallel execution, so I'm a little more confident than <br>
before on its correctness. Still let me remind you to carefully check <br>
your results and feel free to contact me for further details.<br>
<br>
Best regards,<br>
Pietro<br>
<br>
--<br>
Pietro Bonfà<br>
Department of Mathematical, Physical and Computer Sciences<br>
University of Parma<br>
<br>
<br>
<br>
<br>
<br>
On 1/3/21 5:37 AM, William Hewett wrote:<br>
> Hi all,<br>
> <br>
> I'm using QE to calculate the band structure of rare-earth nitride <br>
> materials, currently looking at the effect of nitrogen vacancies and the <br>
> resulting states created. I'm running calculations on 3x3x3 (primitive) <br>
> supercells, then using unfold.x to process the results.<br>
> <br>
> Unfold.x gives an output where each point (k,energy) also has a 'weight' <br>
> clearly this is zero where no states are present and non-zero where they <br>
> are present. My question is, what exactly is this weight? Some sort of <br>
> DOS for a single k-point?<br>
> <br>
> i.e. in the image below most points in the VB (near 6 on the y-axis) are <br>
> quite dark, while some points in the CB (flat 4f bands near -1) are <br>
> lighter in color.<br>
> <br>
> image.png<br>
> <br>
> Kind regards,<br>
> <br>
> Will Hewett<br>
> Post Doctoral Researcher<br>
> Victoria University of Wellington<br>
> New Zealand<br>
> <br>
> <br>
> _______________________________________________<br>
> Quantum ESPRESSO is supported by MaX (<a href="http://www.max-centre.eu" rel="noreferrer" target="_blank">www.max-centre.eu</a>)<br>
> users mailing list <a href="mailto:users@lists.quantum-espresso.org" target="_blank">users@lists.quantum-espresso.org</a><br>
> <a href="https://lists.quantum-espresso.org/mailman/listinfo/users" rel="noreferrer" target="_blank">https://lists.quantum-espresso.org/mailman/listinfo/users</a><br>
> <br>
_______________________________________________<br>
Quantum ESPRESSO is supported by MaX (<a href="http://www.max-centre.eu" rel="noreferrer" target="_blank">www.max-centre.eu</a>)<br>
users mailing list <a href="mailto:users@lists.quantum-espresso.org" target="_blank">users@lists.quantum-espresso.org</a><br>
<a href="https://lists.quantum-espresso.org/mailman/listinfo/users" rel="noreferrer" target="_blank">https://lists.quantum-espresso.org/mailman/listinfo/users</a></blockquote></div>