<div dir="ltr"><span style="color:rgb(0,0,0);font-size:12.8px">I will third the comment that not using a shifted cell might be important (might). In some cells  shifting the k-point origin with MSR1a leads to slightly unbalanced forces which are hard to converge to the "right" symmetric result. If the forces are slightly off, this is an indication that the density is also slightly off. Exactly why this occurs I do not know, I suspect very soft modes associated with numerical errors in finite arithmetic.</span><div style="color:rgb(0,0,0);font-size:12.8px"><br></div><div style="color:rgb(0,0,0);font-size:12.8px">N.B., Wien2k is quite good with these numerical errors. I've noticed that Vasp calculations that collaborators have done often have much larger symmetry breaking.</div><div class="gmail_extra"><br><div class="gmail_quote">On Sun, Nov 26, 2017 at 11:39 AM, Karel Vyborny <span dir="ltr"><<a href="mailto:vybornyk@fzu.cz" target="_blank">vybornyk@fzu.cz</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">I suppose that this does not have to do (much) with centering the mesh.<br>
My guess based on other QMO calculations is that some contributions to<br>
mat. els. of e.g. vx*vx from different parts of the BZ don't cancel<br>
(numerically) even if they actually should.<br>
<br>
It is instructive to repeat the calculation without SO and see how big the<br>
difference between sig_xx and sig_yy (for any M direction) is then.<br>
<br>
Cheers,<br>
<br>
Karel<br>
<br>
<br>
--- x ---<br>
dr. Karel Vyborny<br>
Fyzikalni ustav AV CR, v.v.i.<br>
Cukrovarnicka 10<br>
Praha 6, CZ-16253<br>
tel: +420220318459<br>
<br>
<br>
On Sun, 26 Nov 2017, Fecher, Gerhard wrote:<br>
<br>
> There was a recent discussion on magnetic anisotropy, With a remark by Peter,<br>
> Did you use a Gamma centered k mesh (and enough k points)<br>
><br>
> Ciao<br>
> Gerhard<br>
><br>
> DEEP THOUGHT in D. Adams; Hitchhikers Guide to the Galaxy:<br>
> "I think the problem, to be quite honest with you,<br>
> is that you have never actually known what the question is."<br>
><br>
> ==============================<wbr>======<br>
> Dr. Gerhard H. Fecher<br>
> Institut of Inorganic and Analytical Chemistry<br>
> Johannes Gutenberg - University<br>
> 55099 Mainz<br>
> and<br>
> Max Planck Institute for Chemical Physics of Solids<br>
> 01187 Dresden<br></blockquote></div>-- <br><div class="gmail_signature" data-smartmail="gmail_signature"><div dir="ltr"><div><div dir="ltr"><div><div dir="ltr"><div dir="ltr"><div dir="ltr"><div dir="ltr"><div dir="ltr"><span style="font-size:12.8px">Professor Laurence Marks</span><br></div><div dir="ltr"><span style="font-size:12.8px">"Research is to see what everybody else has seen, and to think what nobody else has thought", </span><span style="font-size:12.8px">Albert Szent-Gyorgi</span><br><a href="http://www.numis.northwestern.edu" target="_blank">www.numis.northwestern.edu</a> ; <span style="font-size:12.8px">Corrosion in 4D: </span><a href="http://MURI4D.numis.northwestern.edu" style="font-size:12.8px" target="_blank">MURI4D.numis.northwestern.edu</a><div><span style="font-size:12.8px">Partner of the CFW 100% program for gender equity, </span><a href="http://www.cfw.org/100-percent" style="font-size:12.8px" target="_blank">www.cfw.org/100-percent</a></div><div>Co-Editor, Acta Cryst A</div></div></div></div></div></div></div></div></div></div></div>
</div></div>