[Wannier] Convergence problem of wannierisation

Arash Mostofi a.mostofi at imperial.ac.uk
Sat Nov 21 21:18:59 CET 2009


Dear Hongbin,

Others who have studied this particular system may be able to give more 
specific advice.

The first things I would check are that you are including a sufficient 
number of conduction bands for the disentanglement, using an appropriate 
inner (or "frozen") wiindow and appropriate initial guesses 
("projections") for your Wannier functions.

Best wishes,

Arash


:------------------------------------------------------------:
: Dr. Arash A. Mostofi           :: a.mostofi at imperial.ac.uk :
: Lecturer and RCUK Fellow       ::                          :
: Depts. of Materials & Physics  ::                          :
: Imperial College London        :: T  +44 (0)207 594 8154   :
: London SW7 2AZ, United Kingdom :: F  +44 (0)207 594 6757   :
:------- http://www.cmth.ph.ic.ac.uk/people/a.mostofi -------:

Zhang, Hongbin wrote:
> Dear all,
> 
> At present I am trying to play with Wannier 90 to calculate some properties like anomalous Hall effect. However, I found that sometimes it is hard to get the convergence of minimazation of both the Omega_I and Omega, that is, disentanglement and wannierization. This is true especially for the latter. I have tried on hcp Co, following the paper by E. Roman (Phys. Rev. Lett. 103, 097203). However, after 5000 iterations, the wannierisation converges to only 1E-7 (enough already? )  like below (Disentanglement has been achieved to 1E-10):
> 
>    4984    -0.148E-06     0.0003041075      106.4767765619   90254.25  <-- CONV
>    4985     0.875E-05     0.0001849761      106.4767853125   90269.00  <-- CONV
>    4986    -0.736E-05     0.0023025636      106.4767779520   90288.01  <-- CONV
>    4987    -0.155E-05     0.0010750558      106.4767764030   90301.06  <-- CONV
>    4988    -0.660E-06     0.0005507459      106.4767757427   90312.06  <-- CONV
>    4989    -0.856E-07     0.0002900457      106.4767756571   90329.45  <-- CONV
>    4990    -0.695E-07     0.0001142533      106.4767755876   90345.29  <-- CONV
>    4991     0.000E+00     0.0001822559      106.4767755865   90350.65  <-- CONV
>    4992    -0.436E-07     0.0001817570      106.4767755429   90369.36  <-- CONV
>    4993    -0.488E-07     0.0001090447      106.4767754941   90383.58  <-- CONV
>    4994     0.000E+00     0.0001206043      106.4767754934   90387.75  <-- CONV
>    4995    -0.697E-07     0.0001203373      106.4767754237   90403.39  <-- CONV
>    4996    -0.312E-06     0.0005823860      106.4767751115   90416.78  <-- CONV
>    4997    -0.471E-07     0.0001267525      106.4767750644   90433.15  <-- CONV
>    4998    -0.410E-07     0.0001274398      106.4767750233   90451.88  <-- CONV
>    4999     0.000E+00     0.0001121615      106.4767750226   90459.01  <-- CONV
>    5000    -0.797E-07     0.0001117099      106.4767749430   90477.35  <-- CONV
> 
> I have tried also different initial projections, but it does not help so much. So could you please tell me how to speed up the convergence process? And how to achieve better convergence?
> 
> I do not see a lot parameters to play with from the input file. Any suggestion is appreciated.
> 
> Best wishes,
> 
> Hongbin Zhang
> 
> ------------------------------------------------------------------------------------------------
> ------------------------------------------------------------------------------------------------
> Forschungszentrum Juelich GmbH
> 52425 Juelich
> Sitz der Gesellschaft: Juelich
> Eingetragen im Handelsregister des Amtsgerichts Dueren Nr. HR B 3498
> Vorsitzende des Aufsichtsrats: MinDir'in Baerbel Brumme-Bothe
> Geschaeftsfuehrung: Prof. Dr. Achim Bachem (Vorsitzender),
> Dr. Ulrich Krafft (stellv. Vorsitzender), Prof. Dr.-Ing. Harald Bolt,
> Prof. Dr. Sebastian M. Schmidt
> ------------------------------------------------------------------------------------------------
> ------------------------------------------------------------------------------------------------
> _______________________________________________
> Wannier mailing list
> Wannier at quantum-espresso.org
> http://www.democritos.it/mailman/listinfo/wannier



More information about the Wannier mailing list