<html>
<head>
<meta http-equiv="Content-Type" content="text/html; charset=utf-8">
</head>
<body>
<div class="BodyFragment"><font size="2"><span style="font-size:10pt;">
<div class="PlainText">Dear Elias,<br>
I don't know it is useful in your case, one should probably look deeper into the band structure, but you may want to give a look to the attached file: it is a new feature (contributed by Gabriele Sclauzero, ETH Zurich) for disentangling only within spherical
regions of k space (typically, near specific k-points).<br>
<br>
The attached PDF is the documentation of the new feature, that is still only in the private development branch and not yet released. If you think it could be useful for you, we can send you the source code.<br>
<br>
Best,<br>
Giovanni<br>
<br>
</div>
</span></font></div>
<div class="BodyFragment"><font size="2"><span style="font-size:10pt;">
<div class="PlainText"><br>
-- <br>
Giovanni Pizzi<br>
Post-doctoral Research Scientist<br>
EPFL STI IMX THEOS<br>
MXC 340 (Bâtiment MXC)<br>
Station 12<br>
CH-1015 Lausanne (Switzerland)<br>
Phone: +41 21 69 31124<br>
<br>
<br>
<br>
<br>
On 19 Feb 2015, at 18:12, Elias Assmann wrote:<br>
<br>
> Hi,<br>
> <br>
> I have a calculation where I find that disentanglement does not<br>
> converge when I use a frozen window. Without the frozen window,<br>
> disentanglement converges quickly:<br>
> <br>
> +---------------------------------------------------------------------+<-- DIS<br>
> | Iter Omega_I(i-1) Omega_I(i) Delta (frac.) Time |<-- DIS<br>
> +---------------------------------------------------------------------+<-- DIS<br>
> 1 15.18420341 13.87841728 9.409E-02 1.42 <-- DIS<br>
> 2 14.48064415 13.51066622 7.179E-02 1.55 <-- DIS<br>
> 3 13.96026696 13.24672414 5.387E-02 1.69 <-- DIS<br>
> 4 13.57678297 13.04599162 4.069E-02 1.82 <-- DIS<br>
> …<br>
> 116 12.23570367 12.23570367 8.045E-11 16.87 <-- DIS<br>
> 117 12.23570367 12.23570367 6.777E-11 17.01 <-- DIS<br>
> 118 12.23570367 12.23570367 5.708E-11 17.14 <-- DIS<br>
> 119 12.23570367 12.23570367 4.809E-11 17.27 <-- DIS<br>
> <br>
> (The complete ‘wout’ file is attached.) But when I turn the frozen<br>
> window on, the procedure seems to get stuck in a loop until it runs<br>
> out of iterations:<br>
> <br>
> +---------------------------------------------------------------------+<-- DIS<br>
> | Iter Omega_I(i-1) Omega_I(i) Delta (frac.) Time |<-- DIS<br>
> +---------------------------------------------------------------------+<-- DIS<br>
> 1 12.67123899 12.48798995 1.467E-02 1.37 <-- DIS<br>
> 2 12.54971488 12.44729442 8.228E-03 1.47 <-- DIS<br>
> 3 12.49131486 12.42411347 5.409E-03 1.57 <-- DIS<br>
> 4 12.45357352 12.41102843 3.428E-03 1.68 <-- DIS<br>
> …<br>
> 997 12.39231118 12.39299287 -5.501E-05 106.04 <-- DIS<br>
> 998 12.39231119 12.39299288 -5.501E-05 106.15 <-- DIS<br>
> 999 12.39231118 12.39299287 -5.501E-05 106.26 <-- DIS<br>
> 1000 12.39231119 12.39299288 -5.501E-05 106.36 <-- DIS<br>
> <br>
> I tried making the frozen window smaller, but the end result is the<br>
> same.<br>
> <br>
> Also, I have a little trouble understanding the output. I would have<br>
> expected “Omega_I(i)” in one row to be equal to “Omega_I(i-1)” in the<br>
> next row. And if the spread oscillates, why is Delta always negative?<br>
> <br>
> I have actually posted about this calculation before:<br>
> <<a href="http://mailman.qe-forge.org/pipermail/wannier/2014-June/000806.html">http://mailman.qe-forge.org/pipermail/wannier/2014-June/000806.html</a>>.<br>
> The problem was that I cannot get a good Wannier projection on an<br>
> N×N×N k-mesh with N even.<br>
> <br>
> I am revisiting it now because for my current application I need the<br>
> even-N (specifically, power-of-2) k-mesh. I noticed that an extrinsic<br>
> band touches my target bands just at the R-point of the BZ (0.5, 0.5,<br>
> 0.5); and that R is included in even-N but not odd-N k-meshes. So I<br>
> tried to use disentanglement with that one extra band.<br>
> <br>
> But disentanglement does not help. The Wannierization seems to have<br>
> the same problems in all the cases.<br>
> <br>
> Finally, I tried including another band in the disentanglement. In<br>
> this case, disentanglement does converge nicely, but the<br>
> Wannierization is not satisfactory either.<br>
> <br>
> <br>
> Elias<br>
> <br>
> <br>
> PS: The immediate motive for this post was that I was the puzzling<br>
> behavior of the disentanglement. But obviously I would also be<br>
> grateful for any other hints on how to converge the even-N<br>
> calculation.<br>
> <br>
> -- <br>
> Elias Assmann<br>
> Institute of Solid State Physics<br>
> Vienna University of Technology<br>
> <<a href="http://www.ifp.tuwien.ac.at/cms/">http://www.ifp.tuwien.ac.at/cms/</a>><br>
> <W14_k8.wout>_______________________________________________<br>
> Wannier mailing list<br>
> Wannier@quantum-espresso.org<br>
> <a href="http://mailman.qe-forge.org/cgi-bin/mailman/listinfo/wannier">http://mailman.qe-forge.org/cgi-bin/mailman/listinfo/wannier</a><br>
<br>
</div>
</span></font></div>
</body>
</html>