<div dir="ltr"><div>Hi,</div><div><br></div><div>Just to add, the 40 referenced in "Now assuming this is the highest among all elements present, take this number and multiply it by 40" should be set to whatever REL_CUTOFF is. I typically use a value of 50 for static cell calculations and 80 for variable cell ones. These figures are based on my experience and from literature on condensed phase simulations. Using the basis set listed in Patrick's post, this would give cutoffs of 520 Ry and 850 Ry with some rounding up to 'nicer' numbers. The point of this is to ensure that the tightest orbitals make it to the finest grid. Importantly, this assumes you don't use the smoothing function in https://manual.cp2k.org/trunk/CP2K_INPUT/ATOM/METHOD/XC/XC_GRID.html.<br></div><div><br></div><div>-T</div><div><br></div><div><br></div><br>On Wednesday, December 4, 2019 at 6:27:18 AM UTC-4, 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"><div>Dear Ramanasish,</div><div><br></div><div>You are attempting to do ab-initio molecular dynamics. At each step, the forces are calculated by performing a single point energy calculation of your system. If your CUTOFF is too small, the single point calculation leads to the wrong electron density, and hence your forces will also be wrong.</div><div><br></div><div>Choosing CUTOFF, you only need to make sure that your energy in a single point calculation is converged. To be completely clear: the latter of the two options you propose.</div><div><br></div><div>In general, look at the highest exponent of the basis you use. For example, if your system contains oxygen, and you use the double-zeta MOLOPT basis set:</div><div> O DZVP-MOLOPT-SR-GTH DZVP-MOLOPT-SR-GTH-q6<br> 1<br> 2 0 2 5 2 2 1<br> 10.389228018317 0.126240722900 0.069215797900 -0.061302037200 -0.026862701100 0.029845227500<br> 3.849621072005 0.139933704300 0.115634538900 -0.190087511700 -0.006283021000 0.060939733900<br> 1.388401188741 -0.434348231700 -0.322839719400 -0.377726982800 -0.224839187800 0.732321580100<br> 0.496955043655 -0.852791790900 -0.095944016600 -0.454266086000 0.380324658600 0.893564918400<br> 0.162491615040 -0.242351537800 1.102830348700 -0.257388983000 1.054102919900 0.152954188700</div><div><br></div><div>The highest exponent is 10.389... in the first row of values, on the very left.</div><div><br></div><div>Now assuming this is the highest among all elements present, take this number and multiply it by 40. That should give you some idea of the appropriate value of CUTOFF, in this case CUTOFF = 400. Run a few ENERGY calculations with values, say, between 350 and 600, and check that the energy is, indeed, converged.</div><div><br></div><div>Yours sincerely,</div><div>Patrick Gono<br></div></div><br><div class="gmail_quote"><div dir="ltr">On Tue, Dec 3, 2019 at 10:42 PM Ramanish Singh <<a href="javascript:" target="_blank" gdf-obfuscated-mailto="bewy8DYXAQAJ" 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 have a question regarding the plane wave CUTOFF in the MGRID section. I am doing NVT MD simulations with a liquid box.</div><div><br></div><div>For deciding the CUTOFF value do I need to run the MD simulations at different values of CUTOFF and then decide the appropriate CUTOFF based on evolution of energy with time for different CUTOFF values</div><div><br></div><div>OR</div><div><br></div><div>a single point energy of the system (<span style="color:rgb(0,0,0)">RUN_TYPE ENERGY) is fine? (as given here: </span><a href="https://www.cp2k.org/howto:converging_cutoff" target="_blank" rel="nofollow" onmousedown="this.href='https://www.google.com/url?q\x3dhttps%3A%2F%2Fwww.cp2k.org%2Fhowto%3Aconverging_cutoff\x26sa\x3dD\x26sntz\x3d1\x26usg\x3dAFQjCNFHeI16YFpXApR5fGmB8YOPSHi7nQ';return true;" onclick="this.href='https://www.google.com/url?q\x3dhttps%3A%2F%2Fwww.cp2k.org%2Fhowto%3Aconverging_cutoff\x26sa\x3dD\x26sntz\x3d1\x26usg\x3dAFQjCNFHeI16YFpXApR5fGmB8YOPSHi7nQ';return true;">https://www.cp2k.org/<wbr>howto:converging_cutoff</a> )</div><div><br></div><div>Please also give any relevant references for the same if available.<br></div><div><br></div><div><br></div><div>Thanks,</div><div>Ramanish Singh</div><div><br></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="bewy8DYXAQAJ" rel="nofollow" onmousedown="this.href='javascript:';return true;" onclick="this.href='javascript:';return true;">cp...@googlegroups.<wbr>com</a>.<br> To view this discussion on the web visit <a href="https://groups.google.com/d/msgid/cp2k/722d97d6-2566-4562-a571-c38dec41ed81%40googlegroups.com?utm_medium=email&utm_source=footer" target="_blank" rel="nofollow" onmousedown="this.href='https://groups.google.com/d/msgid/cp2k/722d97d6-2566-4562-a571-c38dec41ed81%40googlegroups.com?utm_medium\x3demail\x26utm_source\x3dfooter';return true;" onclick="this.href='https://groups.google.com/d/msgid/cp2k/722d97d6-2566-4562-a571-c38dec41ed81%40googlegroups.com?utm_medium\x3demail\x26utm_source\x3dfooter';return true;">https://groups.google.com/d/<wbr>msgid/cp2k/722d97d6-2566-4562-<wbr>a571-c38dec41ed81%<wbr>40googlegroups.com</a>.<br> </blockquote></div> </blockquote></div>