<div dir="ltr"><div dir="ltr"><div dir="ltr"><div class="gmail_default" style="font-family:verdana,sans-serif;color:rgb(0,0,0)">Sorry Fabien, I am dubious that k-points is the issue, I suspect a bad model.</div><div class="gmail_default" style="font-family:verdana,sans-serif;color:rgb(0,0,0)"><br></div><div class="gmail_default" style="font-family:verdana,sans-serif;color:rgb(0,0,0)">Let me explain, repeating some things I have said before or are in the mixer notes albeit perhaps buried.</div><div class="gmail_default" style="font-family:verdana,sans-serif;color:rgb(0,0,0)"><br></div><div class="gmail_default" style="font-family:verdana,sans-serif;color:rgb(0,0,0)">An "easy" assumption that many people make is that dft problems are well posed, so have "simple" fixed-point solutions. In fact this is not true. In reality dft problems have stacks of numerical issues which leads to "noise". You can think about this in terms of finding the minimum of "F(x) = x*x + Constant*Noise(x)" where Noise(x) is some randomly generated noise. The larger Constant is, the larger the deviation is from x=0 of the apparent minimum. Averaging over more k-points may reduce Constant, but often in my experience does relatively little.</div><div class="gmail_default" style="font-family:verdana,sans-serif;color:rgb(0,0,0)"><br></div><div class="gmail_default" style="font-family:verdana,sans-serif;color:rgb(0,0,0)">A symptom of noise is poor convergence -- which is clearly what is occurring here.</div><div class="gmail_default" style="font-family:verdana,sans-serif;color:rgb(0,0,0)"><br></div><div class="gmail_default" style="font-family:verdana,sans-serif;color:rgb(0,0,0)">If you improve the convergence rate you will reduce the noise. Some suggestions:</div><div class="gmail_default" style="font-family:verdana,sans-serif;color:rgb(0,0,0)"><br></div><div class="gmail_default" style="font-family:verdana,sans-serif;color:rgb(0,0,0)">1) Use HDLO and reduce RMT</div><div class="gmail_default" style="font-family:verdana,sans-serif;color:rgb(0,0,0)">2) Reduce RMT and use HDLO</div><div class="gmail_default" style="font-family:verdana,sans-serif;color:rgb(0,0,0)">3) Larger LM (e.g. 8)</div><div class="gmail_default" style="font-family:verdana,sans-serif;color:rgb(0,0,0)">4) Larger gmax in case.in2 (e.g. 16)</div><div class="gmail_default" style="font-family:verdana,sans-serif;color:rgb(0,0,0)">5) Use TEMPS 0.0018 (room temperature)</div><div class="gmail_default" style="font-family:verdana,sans-serif;color:rgb(0,0,0)">6) Peter's suggestions</div><div class="gmail_default" style="font-family:verdana,sans-serif;color:rgb(0,0,0)"><br></div></div><div class="gmail_quote"><div dir="ltr" class="gmail_attr">On Sun, Jan 20, 2019 at 10:27 AM <a href="mailto:tran@theochem.tuwien.ac.at">tran@theochem.tuwien.ac.at</a> <<a href="mailto:tran@theochem.tuwien.ac.at">tran@theochem.tuwien.ac.at</a>> wrote:<br></div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex">Hi,<br>
<br>
If not too expensive, you could try more k-points (20000 and 25000)<br>
to see if the fluctuations persist.<br>
<br>
Questions: how many atoms in the unit cell? Are you optimizing the<br>
position of atoms?<br>
<br>
F. Tran<br>
<br>
On Sunday 2019-01-20 17:17, Victor Zenou wrote:<br>
<br>
>Date: Sun, 20 Jan 2019 17:17:55<br>
>From: Victor Zenou <<a href="mailto:zanov@post.bgu.ac.il" target="_blank">zanov@post.bgu.ac.il</a>><br>
>Reply-To: A Mailing list for WIEN2k users <<a href="mailto:wien@zeus.theochem.tuwien.ac.at" target="_blank">wien@zeus.theochem.tuwien.ac.at</a>><br>
>To: A Mailing list for WIEN2k users <<a href="mailto:wien@zeus.theochem.tuwien.ac.at" target="_blank">wien@zeus.theochem.tuwien.ac.at</a>><br>
>Subject: Re: [Wien] proper k-points for Nd<br>
><br>
>Dear Peter, LDM and wien2k users<br>
>I did some more calculation with energy and CHARGE convergence criteria<br>
>equal to 10^-4 Ryd and 10^-3, respectively. The energy usually converged to<br>
>1*10-5 or even 1*10^-6. Still I get energy fluctuation up to 4*10-2 Ryd<br>
>(please see below)<br>
>I didn't check AFM yet, but this is not the point. The point is to get<br>
>"stable" calculation for optimal k-points.<br>
>I have to say that for each calculation I totally aresed all file except<br>
>case.struc.<br>
>Also I used all defaults as well as gmax=12, also rmt(min)*kmax=8.5, for all<br>
>calculations.<br>
>I can repeat all calculation with smaller RMT, such as 2.5, or to check <br>
>HDLO.<br>
>I think I will start with smaller RMT. Any comments / suggestions?<br>
>Victor<br>
><br>
>K-pointsIBZ E [Ryd], 10-4 CHARGE convergence CHARGE DISTANCE  MMTOT<br>
>3000 208 -77040.50516761 -.000504 0.000825   16.00318<br>
>5000 280 -77040.48748499 -.000544 0.000758   16.68523<br>
>7500 455 -77040.51668342 -.000862 0.000183   16.62547<br>
>10000 540 -77040.48747210 -.000817 0.000499   16.67544<br>
>15000 840 -77040.46976815 -.000216 0.001362   16.64536<br>
><br>
><br>
></blockquote></div><br clear="all"><div><br></div>-- <br><div dir="ltr" class="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>