Dear Jan,<div><br></div><div>I see that in your input file that you use exactly 606 bands for 606 Wannier Functions.</div><div>It would be good to maybe include more bands in your calculation and use the disentanglement</div>
<div>procedure (look at the Wannier User Guide for that).</div><div>Also, if you cell is quite big, using random projections may not be optimal. So if you have an</div><div>idea of what those Wannier Functions should be, there is definitely advantage to input some</div>
<div>carefully chosen initial projections. In the random guess case, the code may get stuck in a local</div><div>minimum that does not correspond to the proper physical minimum you seek.</div><div>As a final note, I would set the convergence threshold to 1.0d-7 at least to both disentanglement</div>
<div>and wannierisation and use a lot of steps just in case (num_iter > 2000 and dis_num_iter > 2000).</div><div><br></div><div>good luck,</div><div><br></div><div>Nicolas</div><div><br><div class="gmail_quote">On Fri, Apr 9, 2010 at 12:02 PM, Jan Felix Binder <span dir="ltr"><<a href="mailto:janfelix.binder@epfl.ch">janfelix.binder@epfl.ch</a>></span> wrote:<br>
<blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex;">Hello,<br>
I used the W90 code with PWscf (QE 4.1.2) to calculate MLWF for a<br>
amorphous structure (gamma point only, 234 atoms, 606 wf).<br>
All WF convered to reasonable values. Most of the centers are at the<br>
right positions, but approximately 10 centers hang somewhere around.<br>
What can be wrong?<br>
<br>
I add the imported lines of the input file:<br>
 gamma_only = .true.<br>
num_iter         =  100<br>
dis_num_iter     =  200<br>
conv_window = 10<br>
conv_noise_amp = 1<br>
conv_noise_num = 10<br>
dis_mix_ratio = 0.3<br>
guiding_centres = .true.<br>
<br>
iprint           =    2<br>
num_dump_cycles  =   10<br>
num_print_cycles =   10<br>
<br>
length_unit      =  Ang<br>
<br>
num_wann = 606<br>
num_bands=606<br>
<br>
exclude_bands : 607-614<br>
<br>
Begin Projections<br>
random<br>
End Projections<br>
<br>
translate_home_cell=.true.<br>
write_xyz=.true.<br>
<br>
begin unit_cell_cart<br>
 bohr<br>
47.65304  0.00000  0.00000<br>
 0.00000 23.82652  0.00000<br>
 0.00000  0.00000 23.82652<br>
end unit_cell_cart<br>
<br>
Thank you for help!<br>
Felix Binder<br>
<br>
EPFL Lausanne<br>
_______________________________________________<br>
Wannier mailing list<br>
<a href="mailto:Wannier@quantum-espresso.org">Wannier@quantum-espresso.org</a><br>
<a href="http://www.democritos.it/mailman/listinfo/wannier" target="_blank">http://www.democritos.it/mailman/listinfo/wannier</a><br>
</blockquote></div><br><br clear="all"><br>-- <br>POILVERT Nicolas<br>PhD candidate,<br>Dpt of Materials Science and Engineering<br>Massachusetts Institute of Technology<br>77, Massachusetts avenue<br>Cambridge, MA 02139<br>
USA<br>work: (617) 452-4212<br><a href="mailto:nicolas.poilvert@gmail.com">nicolas.poilvert@gmail.com</a><br>
</div>