<div dir="ltr">Dear Patrick,<div><br></div><div>Thank you very much for your suggestion. </div><div>I still have a doubt. I think the peculiar behavior is also seen for REL_CUTOFF=60. The convergence is 1e-7 for CUTOFF=550 but after that it suddenly decreases to 1e-4 and 1e-3 at CUTOFF=750. So shall I go with 500 or with 800 for CUTOFF. </div><div><br></div><div>Thanks,</div><div>Ramanish Singh<br><br>On Saturday, June 29, 2019 at 8:53:20 AM UTC-5, Patrick Gono wrote:<blockquote class="gmail_quote" style="margin: 0;margin-left: 0.8ex;border-left: 1px #ccc solid;padding-left: 1ex;"><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">On Thu, 27 Jun 2019 at 16:03, Ramanish Singh <<a href="javascript:" target="_blank" gdf-obfuscated-mailto="ch5panYOAwAJ" rel="nofollow" onmousedown="this.href='javascript:';return true;" onclick="this.href='javascript:';return true;">si...@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># 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># 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># 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> 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> 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> 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> 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> 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> 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> 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> 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> 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> 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> 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> 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> 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="javascript:" target="_blank" gdf-obfuscated-mailto="ch5panYOAwAJ" rel="nofollow" onmousedown="this.href='javascript:';return true;" onclick="this.href='javascript:';return true;">cp...@googlegroups.<wbr>com</a>.<br>
To post to this group, send email to <a href="javascript:" target="_blank" gdf-obfuscated-mailto="ch5panYOAwAJ" rel="nofollow" onmousedown="this.href='javascript:';return true;" onclick="this.href='javascript:';return true;">c...@googlegroups.com</a>.<br>
Visit this group at <a href="https://groups.google.com/group/cp2k" target="_blank" rel="nofollow" onmousedown="this.href='https://groups.google.com/group/cp2k';return true;" onclick="this.href='https://groups.google.com/group/cp2k';return true;">https://groups.google.com/<wbr>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" rel="nofollow" onmousedown="this.href='https://groups.google.com/d/msgid/cp2k/f3462318-02c6-4f95-b97b-97a567821932%40googlegroups.com?utm_medium\x3demail\x26utm_source\x3dfooter';return true;" onclick="this.href='https://groups.google.com/d/msgid/cp2k/f3462318-02c6-4f95-b97b-97a567821932%40googlegroups.com?utm_medium\x3demail\x26utm_source\x3dfooter';return true;">https://groups.google.com/d/<wbr>msgid/cp2k/f3462318-02c6-4f95-<wbr>b97b-97a567821932%<wbr>40googlegroups.com</a>.<br>
For more options, visit <a href="https://groups.google.com/d/optout" target="_blank" rel="nofollow" onmousedown="this.href='https://groups.google.com/d/optout';return true;" onclick="this.href='https://groups.google.com/d/optout';return true;">https://groups.google.com/d/<wbr>optout</a>.<br>
</blockquote></div>
</blockquote></div></div>