[Wien] Possible bugs in Mixer

Peter Blaha peter.blaha at tuwien.ac.at
Sat Dec 20 10:01:06 CET 2025


Clearly, :DIS seems to converge. So "mixer" cannot be the problem, since 
it produces a density which is identical to the new density in the next 
iteration (within 0.0000001).
In any case, on machine2 you still have a non-zero :dis, while on 
machine 1 it is zero to all digits. So scf quality is not EXACTLY the 
same (although for ALL practical purposes it is the same).

Sure, you should check the convergence of the total energy (and other 
quantities like  :FGL, :EFG   :VCOUL :FER ...) as well.

> Eigenenergies are
> 
>      57    0.555640059910   0.55363321799308D-03
>      58    0.555640061603   0.55363321799308D-03
>      59    0.740271026751   0.00000000000000D+00
>      60    0.740271028515   0.00000000000000D+00
>      61    0.748205005635   0.00000000000000D+00
>      62    0.748205006555   0.00000000000000D+00
> -------------------------------------------------------------------
>      57    0.555640058339   0.55363321799308D-03
>      58    0.555640061273   0.55363321799308D-03
>      59    0.740271025639   0.00000000000000D+00
>      60    0.740271027711   0.00000000000000D+00
>      61    0.748205004954   0.00000000000000D+00
>      62    0.748205005071   0.00000000000000D+00

I do not understand these numbers. What should they be ???

Check convergence of eigenvalues with :EIGxxxx


Momentum matrix elements: They are determined only up to a phase (same 
with the eigenvectors). So you cannot compare plain numbers between 2 
machines for these quantitites.

Regards
-- 
-----------------------------------------------------------------------
Peter Blaha,  Inst. f. Materials Chemistry, TU Vienna, A-1060 Vienna
Email: peter.blaha at tuwien.ac.at
WWW:   http://www.imc.tuwien.ac.at      WIEN2k: http://www.wien2k.at
-------------------------------------------------------------------------



More information about the Wien mailing list