<div dir="ltr">Dear Nicola, <div><br></div><div>Thank you for your reply, I have tried before to project the wannier function without spin-orbit, I can get smooth bands. But in that case, I used disentanglement method instead of excluding bands. Simply excluding the bands may cause the problem. Aron also suggests me to use disentanglement method and set sp3 hybrid as initial projection. </div><div><br></div><div>Thank you for your suggestions.</div><div><br></div><div>Best regards,</div><div><br></div><div>Yu</div><div class="gmail_extra"><br><div class="gmail_quote">On Wed, Feb 11, 2015 at 10:29 PM, Nicola Marzari <span dir="ltr"><<a href="mailto:nicola.marzari@epfl.ch" target="_blank">nicola.marzari@epfl.ch</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><br>
Dear Yu,<br>
<br>
<br>
my first guess is that if you want to disentangle the 5d bands,<br>
and use 12 wannier functions, you need to include the<br>
bottom of the parabolic "s" band - i.e. bands 1-2,<br>
as in Fig. 4a of the Souza et al 2002 PRB.<br>
<br>
Why don't you try this first in gold without spin-orbit, and see<br>
if you can reproduce the results above (those where for copper,<br>
but the same considerations apply).<br>
<br>
                        nicola<div><div><br>
<br>
<br>
On 11/02/2015 20:16, Yu Zhang wrote:<br>
</div></div><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div><div>
Dear wannier90 users and developers,<br>
<br>
I am using wannier90-1.2 and vasp5.3 to obtain the wannier wave function<br>
of gold with spin-orbital coupling. I used d5s orbital as target.<br>
However, it ends up with some wiggly bands as shown by the attached<br>
figure. I have tried different projection (d5), and k-mesh, but I still<br>
failed to get the smooth bands. Does anyone know what causes the wiggly<br>
bands and how to get rid of them?<br>
<br>
<br>
My wannier90 input file is:<br>
==============================<u></u>=============<br>
num_wann  = 12  ! set to NBANDS by VASP<br>
num_bands = 12<br>
<br>
exclude_bands : 1-2, 15-16<br>
<br>
spinors = .true.<br>
Begin Projections<br>
Au: l=0;l=2<br>
End Projections<br>
<br>
begin unit_cell_cart<br>
      0.0000000     2.0340000     2.0340000<br>
      2.0340000     0.0000000     2.0340000<br>
      2.0340000     2.0340000     0.0000000<br>
end unit_cell_cart<br>
<br>
begin atoms_cart<br>
Au       0.0000000     0.0000000     0.0000000<br>
end atoms_cart<br>
<br>
mp_grid =    21    21    21<br>
<br>
begin kpoints<br>
       0.000000000000      0.000000000000      0.000000000000<br>
       0.047619047619      0.000000000000      0.000000000000<br>
       0.095238095238      0.000000000000      0.000000000000<br>
       0.142857142857      0.000000000000      0.000000000000<br>
       0.190476190476      0.000000000000      0.000000000000<br>
.......<br>
==============================<u></u>====================<br>
<br>
Thank you very much in advance!<br>
<br>
Best regards,<br>
<br>
Yu Zhang<br>
<br>
<br></div></div>
______________________________<u></u>_________________<br>
Wannier mailing list<br>
<a href="mailto:Wannier@quantum-espresso.org" target="_blank">Wannier@quantum-espresso.org</a><br>
<a href="http://mailman.qe-forge.org/cgi-bin/mailman/listinfo/wannier" target="_blank">http://mailman.qe-forge.org/<u></u>cgi-bin/mailman/listinfo/<u></u>wannier</a><br>
<br><span><font color="#888888">
</font></span></blockquote><span><font color="#888888">
<br>
-- <br>
<br>
------------------------------<u></u>------------------------------<u></u>----------<br>
Prof Nicola Marzari, Chair of Theory and Simulation of Materials, EPFL<br>
</font></span></blockquote></div><br></div></div>