<div dir="ltr">Dear Kee Chon Wee,<br><br>Thank you so much for your help.<br><br>On Tuesday, June 11, 2019 at 1:54:09 PM UTC+4:30, Kee Choon Wee wrote:<blockquote class="gmail_quote" style="margin: 0;margin-left: 0.8ex;border-left: 1px #ccc solid;padding-left: 1ex;"><div dir="ltr">Hi Shakiba,<div><br></div><div>The purpose of a frequency calculation is to verify that the stationary point you found is not a saddle point (1 imaginary frequency, i.e a TS) or higher-order saddle point with (>1 imaginary frequency).</div><div>I don't think it has anything to do with getting struck in a local minima, In fact, if you have no imaginary frequency the structure you found is a minima (could be a global/local one).</div><div>You can try a short MD to see if there is a large transition within say a picosecond, this might give you a hint as to whether it is a minima or transition state.</div><div><br></div></div><br><div class="gmail_quote"><div dir="ltr">On Mon, Jun 10, 2019 at 8:32 PM Mohammad Shakiba <<a href="javascript:" target="_blank" gdf-obfuscated-mailto="0lh_KtYKAQAJ" rel="nofollow" onmousedown="this.href='javascript:';return true;" onclick="this.href='javascript:';return true;">mshakiba....@gmail.<wbr>com</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"><div>Hello everyone,</div><div><br></div><div>When optimizing with CP2K, we set a force limit and a maximum step size (for example less than 10meV and 0.002Bohr) and when the calculations reached those limits it is done. Before using CP2K I used Gaussian 09. When the geometry optimization was fully performed, I performed frequency calculation in order to check that there is no imaginary frequency so the optimization is not stuck in the local minima. However, when it came to big molecules, the rwf file tended to be hundreds of GB, so the calculations was not possible. Now, my question is:</div><div><br></div><div><br></div><div>Is it necessary to perform frequency calculations (like phonon calculations using phonopy) for the optimized geometry in CP2K? I mean the BFGS algorithm for geometry optimization is guaranteed to not stuck in the local minima? For example my system contains a quantum dot with many organic ligands (near 300 atoms) and the optimization is fully performed but I don't know if it would be acceptable to allege that this is the optimum geometry!  To sum up: when can we say that this is optimum geometry for this specific geometry configuration?</div><div><br></div><div>Thanks in advance. <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="0lh_KtYKAQAJ" 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="0lh_KtYKAQAJ" 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/9ff9bed8-3d0a-4ca5-b833-5435e3b9e54c%40googlegroups.com?utm_medium=email&utm_source=footer" target="_blank" rel="nofollow" onmousedown="this.href='https://groups.google.com/d/msgid/cp2k/9ff9bed8-3d0a-4ca5-b833-5435e3b9e54c%40googlegroups.com?utm_medium\x3demail\x26utm_source\x3dfooter';return true;" onclick="this.href='https://groups.google.com/d/msgid/cp2k/9ff9bed8-3d0a-4ca5-b833-5435e3b9e54c%40googlegroups.com?utm_medium\x3demail\x26utm_source\x3dfooter';return true;">https://groups.google.com/d/<wbr>msgid/cp2k/9ff9bed8-3d0a-4ca5-<wbr>b833-5435e3b9e54c%<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>