<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>