<div dir="auto"><div dir="auto">Hello Samuel,</div><div dir="auto"><br></div><div dir="auto">unfortunately these problems are not isolated to only you. Several users including myself have had the same problem. </div><div dir="auto"><br></div><div dir="auto">I suspect most of it comes from an over eagerness to fit the points as best as possible while getting the most localised wf. </div><div dir="auto"><br></div><div dir="auto">While this sounds like what we want, only fitting the eigenvalues means that the description between energy eigenvalues, the bands, like local symmetry and derivatives are not conserved, imposing these properties would likely give a less localised but more 'true' approximant. </div><div dir="auto"><br></div><div dir="auto">Unfortunately I don't have a programmatic solution for the authors of w90 but I can tell you a 'hack' to get some of your fits looking better. The first step is to realise that when your dft package plots the bands it is actually doing a non self consistent probing of the energy landscape using your current configuration. If you complete a run with a coarse mesh, you can restart with a fine one (this is often great for large problems with lattice distortions). With this in mind, calculate your converged system with a mesh you like, then take your converged charge density, ionic positions with a large number of kpoints. </div><div dir="auto"><br></div><div dir="auto">The eigenvalues you get from this should look a lot like a regular band plot. As w90 tries to fit all of these kpoints exactly you can force it into a more realistic space. Not ideal but it works.</div><div dir="auto"><br></div><div dir="auto">Cheers,</div><div dir="auto">Jesse Vaitkus</div></div><div class="gmail_extra"><br><div class="gmail_quote">On 12/12/2016 7:39 am, "Samuel Ponce" <<a href="mailto:samuel.ponce@materials.ox.ac.uk">samuel.ponce@materials.ox.ac.uk</a>> wrote:<br type="attribution"><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">




<div>
<div style="direction:ltr;font-family:Tahoma;color:#000000;font-size:10pt">Dear Wannier team,<br>
<div>
<div>
<div>
<div>
<div>
<div>
<div>
<div>
<div>
<div>
<div><br>
</div>
I have a small question about Wannier90. <br>
<br>
</div>
Somehow I find it very difficult to get very good Wannierization of Silicon. I'm interested in the top of the VB and the bottom of the CB (between G-X). 
<br>
<br>
</div>
I can get okish one but when I try to improve them, I do not manage. <br>
<br>
</div>
Attach to this email is a comparison between the L-G-X path of Silicon with PW/bands.x (in green) VS Wannier90 in red.
<br>
<br>
</div>
As you can see, the Wannier code fail to reproduce the second branch of the CBM (on the right), which is problematic. The wannier seems to have some unphysical oscillation in the CB.
<br>
<br>
</div>
Attach to this email are all the file needed to obtain both plots and even the line needed to do the gnuplot (you have to scale them). In each of them I've place a README file with all the steps for convenience. I've also placed the psp.
<br>
</div>
<div>The calculations are very fast. <br>
</div>
<div><br>
</div>
Maybe its something very stupid. <br>
<br>
</div>
So far, I've try:<br>
</div>
- increasing the nscf grid (you can see a comparison between 8x8x8 and 10x10x10 grid. The latest 10x10x10 gives dramatic oscillations)<br>
</div>
<div>- changing the number of bands<br>
</div>
- playing with the frozen windows<br>
</div>
<div>- increasing the number of Wannier iteration<br>
</div>
<div><br>
</div>
<div>Nothing made it better than that. <br>
<br>
</div>
<div>Would you have any suggestions? <br>
<br>
</div>
<div>PS: The lattice parameter is slightly lower (10.17 Bohr) than the DFT one (10.208 Bohr) on purpose but is the same in both case. In the case of the DFT lattice 10.208 similar issues exists but its slightly better.
<br>
<br>
</div>
<div>PPS: I also noticed that the Wannierzation was breaking the symmetries. For example doubly degenerate bands along G-X were slightly lifted. As a results, effective masses are quite bad. Is there a way to impose crystal symmetries?<br>
</div>
<div><br>
</div>
<div>Thank you,<br>
<br>
</div>
Samuel</div>
</div>

<br><br>______________________________<wbr>_________________<br>
Wannier mailing list<br>
<a href="mailto:Wannier@quantum-espresso.org">Wannier@quantum-espresso.org</a><br>
<a href="http://mailman.qe-forge.org/cgi-bin/mailman/listinfo/wannier" rel="noreferrer" target="_blank">http://mailman.qe-forge.org/<wbr>cgi-bin/mailman/listinfo/<wbr>wannier</a><br>
<br></blockquote></div></div>