<div dir="ltr">Dear Ramanish,<div><br></div><div>Looking at the rest of your data, it appears that the problem was a value of the relative cutoff that was too small. For example, if you consider the convergence of energy with respect to the total cutoff for REL_CUTOFF = 60, the suspicious behaviour vanishes and even small values of CUTOFF yield converged energies.</div><div><br></div><div>I would suggest you stick with REL_CUTOFF = 50 or 60. These are also the values suggested in the input manual to get accurate results. For the absolute cutoff, you are better off using higher values than the default one as well. Some even suggest using values as high as 800 or 1000 Ry, if accuracy is of the highest importance. The correct choice of the value will also depend on the system (e.g. its size, or the elements present). In any case, from your tests it appears that smaller values of CUTOFF are sufficient as long as you use REL_CUTOFF = 50 or higher, though.</div><div><br></div><div>What I personally do is start with a reasonably high value for CUTOFF (e.g. 700 or 800 Ry). If I find that a few extra percent of performance are critical for my calculations, I might test a smaller value for the CUTOFF by looking at the quantity of my interest. For example, if I were to study binding energies of small molecules on surfaces, I would study the effect of a smaller value of CUTOFF on that exact quantity. If I find that the effect is negligible, I might use a smaller cutoff (e.g. 500, or even 350 Ry).</div><div><br></div><div>However, from my experience, other parameters, such as the exchange and correlation functional used, or the basis sets, have a much stronger effect on the computational resources required.</div><div><br></div><div>Hope this helps.</div><div>Yours sincerely,</div><div>Patrick Gono</div><div><br></div><div><br></div></div><br><div class="gmail_quote"><div dir="ltr" class="gmail_attr">On Thu, 27 Jun 2019 at 16:03, Ramanish Singh <<a href="mailto:sing...@umn.edu">sing...@umn.edu</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"><div dir="ltr">Hi,<div><br></div><div>I am simulating a system of 64 Iodine atoms. I performed the cutoff test to find the suitable cutoff and reliable_cutoff but the data I got is somewhat weird. I have attached the data with this post.</div><div>For example, for REL_CUTOFF=30 I got the following data: </div><div><p style="font-stretch:normal;font-size:11px;line-height:normal;font-family:Menlo;color:rgb(86,32,244)"><span style="font-variant-ligatures:no-common-ligatures"># Number of Grids: 4</span></p>
<p style="font-stretch:normal;font-size:11px;line-height:normal;font-family:Menlo;color:rgb(86,32,244)"><span style="font-variant-ligatures:no-common-ligatures"># Relative Cutoff (Ry): 30</span></p>
<p style="font-stretch:normal;font-size:11px;line-height:normal;font-family:Menlo;color:rgb(86,32,244)"><span style="font-variant-ligatures:no-common-ligatures"># Cutoff (Ry) | Total Energy (Ha) | NG on grid 1 | NG on grid 2 | NG on grid 3 | NG on grid 4</span></p>
<p style="font-stretch:normal;font-size:11px;line-height:normal;font-family:Menlo;color:rgb(0,0,0)"><span style="font-variant-ligatures:no-common-ligatures"> 200.00 -725.7833908007 0 15624 15391 4590</span></p>
<p style="font-stretch:normal;font-size:11px;line-height:normal;font-family:Menlo;color:rgb(0,0,0)"><span style="font-variant-ligatures:no-common-ligatures"> 250.00 -725.7860435954 0 12298 13049 10258</span></p>
<p style="font-stretch:normal;font-size:11px;line-height:normal;font-family:Menlo;color:rgb(0,0,0)"><span style="font-variant-ligatures:no-common-ligatures"> 300.00 -725.7849325650 0 7090 18177 10338</span></p>
<p style="font-stretch:normal;font-size:11px;line-height:normal;font-family:Menlo;color:rgb(0,0,0)"><span style="font-variant-ligatures:no-common-ligatures"> 350.00 -725.7837028953 0 3652 19891 12062</span></p>
<p style="font-stretch:normal;font-size:11px;line-height:normal;font-family:Menlo;color:rgb(0,0,0)"><span style="font-variant-ligatures:no-common-ligatures"> 400.00 -725.7839064559 0 2188 17023 16394</span></p>
<p style="font-stretch:normal;font-size:11px;line-height:normal;font-family:Menlo;color:rgb(0,0,0)"><span style="font-variant-ligatures:no-common-ligatures"> 450.00 -725.7853932286 0 1619 17592 16394</span></p>
<p style="font-stretch:normal;font-size:11px;line-height:normal;font-family:Menlo;color:rgb(0,0,0)"><span style="font-variant-ligatures:no-common-ligatures"> 500.00 -725.8084189293 0 527 16090 18988</span></p>
<p style="font-stretch:normal;font-size:11px;line-height:normal;font-family:Menlo;color:rgb(0,0,0)"><span style="font-variant-ligatures:no-common-ligatures"> 550.00 -725.7834120461 0 0 16617 18988</span></p>
<p style="font-stretch:normal;font-size:11px;line-height:normal;font-family:Menlo;color:rgb(0,0,0)"><span style="font-variant-ligatures:no-common-ligatures"> 600.00 -725.7833214512 0 0 15624 19981</span></p>
<p style="font-stretch:normal;font-size:11px;line-height:normal;font-family:Menlo;color:rgb(0,0,0)"><span style="font-variant-ligatures:no-common-ligatures"> 650.00 -725.7826948381 0 0 15624 19981</span></p>
<p style="font-stretch:normal;font-size:11px;line-height:normal;font-family:Menlo;color:rgb(0,0,0)"><span style="font-variant-ligatures:no-common-ligatures"> 700.00 -725.7826948381 0 0 15544 20061</span></p>
<p style="font-stretch:normal;font-size:11px;line-height:normal;font-family:Menlo;color:rgb(0,0,0)"><span style="font-variant-ligatures:no-common-ligatures"> 750.00 -725.7844897119 0 0 12298 23307</span></p>
<p style="font-stretch:normal;font-size:11px;line-height:normal;font-family:Menlo;color:rgb(0,0,0)"><span style="font-variant-ligatures:no-common-ligatures"> 800.00 -725.7829016340 0 0 12218 23387</span></p></div><div><br></div><div>The energy is exactly the sane for CUTOFF = 650 and 700 Ry , but it changes for CUTOFF=750 and 800 ( with error more than 1e-3).</div><div>Similar thing happens for other values of REL_CUTOFF as well. The error decreases to 1e-7 as CUTOFF increases but suddenly jumps to 1e-3 when CUTOFF is increased beyond 600 Ry.</div><div>Do I need to increase my CUTOFF beyond 800 to see if it converges? And is it worth it to keep it that high keeping in mind the increase in the simulation time?</div><div><br></div><div>Thanks,</div><div>Ramanish</div></div>
<p></p>
-- <br>
You received this message because you are subscribed to the Google Groups "cp2k" group.<br>
To unsubscribe from this group and stop receiving emails from it, send an email to <a href="mailto:cp...@googlegroups.com" target="_blank">cp...@googlegroups.com</a>.<br>
To post to this group, send email to <a href="mailto:cp...@googlegroups.com" target="_blank">cp...@googlegroups.com</a>.<br>
Visit this group at <a href="https://groups.google.com/group/cp2k" target="_blank">https://groups.google.com/group/cp2k</a>.<br>
To view this discussion on the web visit <a href="https://groups.google.com/d/msgid/cp2k/f3462318-02c6-4f95-b97b-97a567821932%40googlegroups.com?utm_medium=email&utm_source=footer" target="_blank">https://groups.google.com/d/msgid/cp2k/f3462318-02c6-4f95-b97b-97a567821932%40googlegroups.com</a>.<br>
For more options, visit <a href="https://groups.google.com/d/optout" target="_blank">https://groups.google.com/d/optout</a>.<br>
</blockquote></div>