<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" class="gmail_attr">On Tue, Dec 3, 2019 at 10:42 PM 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 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">https://www.cp2k.org/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="mailto:cp...@googlegroups.com" target="_blank">cp...@googlegroups.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">https://groups.google.com/d/msgid/cp2k/722d97d6-2566-4562-a571-c38dec41ed81%40googlegroups.com</a>.<br>
</blockquote></div>