<div dir="rtl"><div dir="ltr">Yes, I can try 20,000 and 25,000 or even 30,000, but I never had to go for that large number.</div><div dir="ltr">There are 4 atoms with fixed positions : 2a (0,0,0), 2c (1/3,2/3,1/4),  so relaxation of atoms can't be done.</div><div dir="ltr">I have to mention its not at its optimal volume, which will be done later.</div><div dir="ltr">Victor</div></div><br><div class="gmail_quote"><div dir="rtl" class="gmail_attr">‫בתאריך יום א׳, 20 בינו׳ 2019 ב-18:27 מאת <‪<a href="mailto:tran@theochem.tuwien.ac.at">tran@theochem.tuwien.ac.at</a>‬‏>:‬<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>
>_______________________________________________<br>
Wien mailing list<br>
<a href="mailto:Wien@zeus.theochem.tuwien.ac.at" target="_blank">Wien@zeus.theochem.tuwien.ac.at</a><br>
<a href="http://zeus.theochem.tuwien.ac.at/mailman/listinfo/wien" rel="noreferrer" target="_blank">http://zeus.theochem.tuwien.ac.at/mailman/listinfo/wien</a><br>
SEARCH the MAILING-LIST at:  <a href="http://www.mail-archive.com/wien@zeus.theochem.tuwien.ac.at/index.html" rel="noreferrer" target="_blank">http://www.mail-archive.com/wien@zeus.theochem.tuwien.ac.at/index.html</a><br>
</blockquote></div>