<div dir="auto">I suggest using the "Check" command, which shows more useful information and post the output. Without that it is hard to guess.<br><br><div data-smartmail="gmail_signature">_____<br>Professor Laurence Marks<br>"Research is to see what everybody else has seen, and to think what nobody else has thought", Albert Szent-Gyorgi<br><a href="http://www.numis.northwestern.edu">www.numis.northwestern.edu</a></div></div><div class="gmail_extra"><br><div class="gmail_quote">On Dec 7, 2017 4:20 AM, "Pavel Ondračka" <<a href="mailto:pavel.ondracka@email.cz">pavel.ondracka@email.cz</a>> wrote:<br type="attribution"><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">Dear Wien2k mailing list,<br>
<br>
I'm seeing some strange energy oscillations in one calculation near the<br>
end when everything seems almost converged.<br>
<br>
:ENERGY convergence:  0 0.0001 .0039009700000000<br>
:ENERGY convergence:  0 0.0001 .0168906450000000<br>
:ENERGY convergence:  0 0.0001 .0096598000000000<br>
:ENERGY convergence:  0 0.0001 .0162053350000000<br>
:ENERGY convergence:  0 0.0001 .0089683200000000<br>
:ENERGY convergence:  0 0.0001 55.6398533450000000<br>
:ENERGY convergence:  0 0.0001 27.8210445050000000<br>
:ENERGY convergence:  0 0.0001 27.8235252600000000<br>
<br>
the charge is converged already:<br>
<br>
:CHARGE convergence:  1 0.001 -.000920<br>
:CHARGE convergence:  1 0.001 -.000977<br>
:CHARGE convergence:  1 0.001 -.000933<br>
:CHARGE convergence:  1 0.001 -.000972<br>
:CHARGE convergence:  1 0.001 -.000984<br>
:CHARGE convergence:  1 0.001 -.000988<br>
:CHARGE convergence:  1 0.001 -.000997<br>
:CHARGE convergence:  1 0.001 -.000998<br>
<br>
This is quite large calculation of TiSiO amorphous-like cell (96 atoms,<br>
 struct file attached, compressed in order to pass the moderation<br>
queue). The numerical parameters are quite standard,<br>
2x2x2 shifted k-grid (can also reproduce this with 3x3x3), Rkmax 7.5,<br>
PBE, rest of stuff at default. I have done a lot of calculation with<br>
similar cells and such parameters without any problems (albeit mostly<br>
some time ago with older Wien2k versions)...<br>
<br>
I'm trying to calculate the binding energy of the O 1s for the first<br>
oxygen using the Slaters transition state approach. Surprisingly when I<br>
use almost touching spheres it converges (O:1.48, Ti: 1.72, Si:1.48),<br>
however when I reduce them by 3% (as in the attached struct file) this<br>
problems surfaces. However I need this reduced spheres for comparison<br>
with another calculation. BTW the spheres around 1.5 for oxygen<br>
should be reasonable anyway?<br>
<br>
No errors or warnings anywhere, except for the expected<br>
WARN  :        CHARGED CELL with  -0.500<br>
<br>
This is with Wien2k 17.1, compiled with MKL 2017.0.098 and ifort.<br>
<br>
Any ideas what might be causing this problem?<br>
Best regards<br>
Pavel</blockquote></div></div>