<html><head><meta http-equiv="Content-Type" content="text/html; charset=utf-8"></head><body style="word-wrap: break-word; -webkit-nbsp-mode: space; line-break: after-white-space;" class="">Dear Linfeng, <div class=""><br class=""></div><div class="">a couple of comments. Setting EPS_SCF to ridiculous low values below 1E-8 makes no sense - you are basically „optimizing“ noise. </div><div class="">Using OT, for a geo_opt EPS_SCF 1E-7 should be fine but you should also set EPS_DEFAULT to 1E-12. Also setting all MAX_* and </div><div class="">RMS_* values to 1E-4 is a rather stringent criterion, I typically go with 1/3 of the default values, i.e. </div><div class=""><br class=""></div><div class="">MAX_DR 1E-3</div><div class="">MAX_FORCE 1.5E-4</div><div class="">RMS_DR 5E-4</div><div class="">RMS_FORCE 1E-4</div><div class=""><br class=""></div><div class="">Lastly, within CP2K LBFGS is typically faster converging than BFGS, though slightly less intuitive in its usage. </div><div class=""><br class=""></div><div class="">Cheers, </div><div class="">Thomas<br class=""><div><br class=""><blockquote type="cite" class=""><div class="">Am 18.05.2020 um 09:45 schrieb Linfeng Gan <<a href="mailto:linfe...@gmail.com" class="">linfe...@gmail.com</a>>:</div><br class="Apple-interchange-newline"><div class=""><div dir="ltr" class=""><div class="">Here is the result.</div><div class=""><br class=""></div><br class="">On Monday, May 18, 2020 at 8:41:45 AM UTC+8, Linfeng Gan wrote:<blockquote class="gmail_quote" style="margin: 0;margin-left: 0.8ex;border-left: 1px #ccc solid;padding-left: 1ex;"><div dir="ltr" class=""><br class=""><div class=""><span lang="EN-US" class=""> </span><br class="webkit-block-placeholder"></div><div class=""><span lang="EN-US" class=""> </span><br class="webkit-block-placeholder"></div><p class="MsoNormal"><span lang="EN-US" class="">Thank you for your advice, Travis.</span></p><p class="MsoNormal"><span lang="EN-US" class="">I had tested EPS_SCF from 1e-5 to 1e-14. I posted
the 1e-5 one only because it took the shortest time to get the final structure.
And I will try convergence criteria with 1e-4 as you suggest. </span></p><div class=""><span lang="EN-US" class=""> </span><br class="webkit-block-placeholder"></div><p class="MsoNormal"><span lang="EN-US" class="">Linfeng</span></p>
<div class=""><br class=""></div><div class=""><br class=""></div>On Sunday, May 17, 2020 at 12:36:33 PM UTC+8, Travis wrote:<blockquote class="gmail_quote" style="margin:0;margin-left:0.8ex;border-left:1px #ccc solid;padding-left:1ex"><div dir="ltr" class="">Hi,<div class=""><br class=""></div><div class="">I see you have EPS_SCF commented out in your input file. The default convergence threshold is 1e-5 au which is usually not tight enough for accurate forces (and certainly not for vibrational analysis, which is computed numerically - hence it will take ages over Gaussian). The value you commented out is actually better and may take less steps to converge. Your criteria for convergence are also extremely tight, 1e-4 au is typically enough. Also, just to be clear, CP2K when setup appropriately will reproduce Gaussian output for the total energies for non-periodic calculations. So you CAN use Gaussian for the gas phase reaction and CP2K for the surface reaction.</div><div class=""><br class=""></div><div class="">-T</div><div class=""><br class=""><br class="">On Saturday, May 16, 2020 at 10:47:56 PM UTC-4, Linfeng Gan wrote:<blockquote class="gmail_quote" style="margin:0;margin-left:0.8ex;border-left:1px #ccc solid;padding-left:1ex"><div dir="ltr" class=""><br class="">Thank you very much Vladimir,<br class=""><div class=""><br class=""></div><div class="">1) I know CP2K is not a choice for small molecules. But I am investigating an organic reaction which could take place in gas phase and on C3N4 surface. I have to compare the data between the gas phase and the surface. So, I can’t use gaussian or orca to that. DMol3 or CASTEP maybe is best choice for me, but I think we should give open source software a chance. So, CP2K is my choice right now. I am not very familiar with CP2K that it’s why I go to this forum and ask for some help. </div><br class="">2) Agree. <br class=""><br class="">3) Yes, it was. <br class=""><br class="">4) 10000 times is an exaggeration, I apologize for that. Form the same initial structure, it took gaussian 41 seconds with a single Xeon 2680v2 CPU to finish the geometry optimization, but it took CP2k 40 minutes with 20 CPU to do the same thing. If we include the frequencies calculation the time is 4 minutes vs 71 minutesX20CPU. A single SCF calculation of CP2K was pretty fast, but most of time was wasting in locating the minimize. <br class="">This efficiency issue only happened in this 1-butene molecular. Another molecular I tested, cyclobutene, seemed much batter than 1-butene. <br class=""><br class="">I attach all the input and the output file. So, you guy can help me to improve this. Thanks all.<br class=""><br class=""><br class="">Best,<br class="">LInfeng<br class=""><br class=""><br class="">On Friday, May 15, 2020 at 7:41:12 PM UTC+8, Vladimir Rybkin wrote:<blockquote class="gmail_quote" style="margin:0;margin-left:0.8ex;border-left:1px #ccc solid;padding-left:1ex"><div dir="ltr" class="">Dear Linfeng,<div class=""><br class=""></div><div class="">a few general remarks:</div><div class="">1) CP2K is not the most efficient tool for small molecules. For instance, it does a lot of integrations over cells in real and reciprocal space, and for a small molecule the cell is almost everywhere empty.</div><div class="">2) CP2K does not use internal coordinates for optimisation. This has to do with the fact that it's mostly meant for large periodic systems. Gaussians' Berny algorithm does use internal coordinates and is very efficient for small systems.</div><div class="">3) It's very likely that you 150 steps have not change the geometry significantly. </div><div class="">4) It's not completely clear where this 10000 times come from. It may have to do with how you compile and run your applications.</div><div class=""><br class=""></div><div class="">That said, there's no specific problems with optimisers and efficiency. It's about applicability of tools for specific purposes. If you going to work with small molecules, you'd better stick to programs made for them, i.e. Gaussian. CP2K is meant for large and/or periodic systems.</div><div class=""><br class=""></div><div class="">Yours,</div><div class=""><br class=""></div><div class="">Vladimir<br class=""><br class="">пятница, 15 мая 2020 г., 2:51:51 UTC+2 пользователь Linfeng Gan написал:<blockquote class="gmail_quote" style="margin:0;margin-left:0.8ex;border-left:1px #ccc solid;padding-left:1ex"><div dir="ltr" class=""><p class="MsoNormal"><span lang="EN-US" class="">Hi all,</span></p><p class="MsoNormal"><span lang="EN-US" class="">When I optimized 1-butene molecular at
B3lYP/6-311G** level, it took much much longer than gaussian16 did. <span class=""> </span>I had test optimizer CG, BDGS and LBFGD,
nothing better. It seems the optimizers have some efficiency issues. All three
optimizer had taken at least 150 optimization steps to get the final structure,
but the initial structure had been optimized by gaussian. <span class=""> </span><span class=""> </span>How
could I improve that?</span></p><div class=""><span lang="EN-US" class=""> </span><br class="webkit-block-placeholder"></div><p class="MsoNormal"><span lang="EN-US" class="">Best</span></p><p class="MsoNormal"><span lang="EN-US" class="">Linfeng</span></p>
</div></blockquote></div></div></blockquote></div></blockquote></div></div></blockquote></div></blockquote></div><div class=""><br class="webkit-block-placeholder"></div>
-- <br class="">
You received this message because you are subscribed to the Google Groups "cp2k" group.<br class="">
To unsubscribe from this group and stop receiving emails from it, send an email to <a href="mailto:cp...@googlegroups.com" class="">cp...@googlegroups.com</a>.<br class="">
To view this discussion on the web visit <a href="https://groups.google.com/d/msgid/cp2k/e60c3ce3-6df6-48f1-8026-460bfc1a7798%40googlegroups.com?utm_medium=email&utm_source=footer" class="">https://groups.google.com/d/msgid/cp2k/e60c3ce3-6df6-48f1-8026-460bfc1a7798%40googlegroups.com</a>.<br class="">
<span id="cid:3d41f338-9400-4d26-b4f6-2f20dcfd7972"><opt_geo.out></span><span id="cid:d5d73003-bec2-42fb-958e-0b03967f3d6d"><opt_geo.inp></span></div></blockquote></div><br class=""></div><br class=""><br class=""><div class="">
<div style="color: rgb(0, 0, 0); letter-spacing: normal; text-align: start; text-indent: 0px; text-transform: none; white-space: normal; word-spacing: 0px; -webkit-text-stroke-width: 0px; word-wrap: break-word; -webkit-nbsp-mode: space; -webkit-line-break: after-white-space;" class=""><div style="color: rgb(0, 0, 0); letter-spacing: normal; text-align: start; text-indent: 0px; text-transform: none; white-space: normal; word-spacing: 0px; -webkit-text-stroke-width: 0px; word-wrap: break-word; -webkit-nbsp-mode: space; -webkit-line-break: after-white-space;" class=""><div style="color: rgb(0, 0, 0); letter-spacing: normal; text-align: start; text-indent: 0px; text-transform: none; white-space: normal; word-spacing: 0px; -webkit-text-stroke-width: 0px; word-wrap: break-word; -webkit-nbsp-mode: space; -webkit-line-break: after-white-space;" class=""><div style="color: rgb(0, 0, 0); font-family: Helvetica; font-size: 12px; font-style: normal; font-variant-ligatures: normal; font-variant-caps: normal; font-variant-east-asian: normal; font-variant-position: normal; font-weight: normal; letter-spacing: normal; line-height: normal; text-align: start; text-indent: 0px; text-transform: none; white-space: normal; word-spacing: 0px; -webkit-text-stroke-width: 0px; orphans: 2; widows: 2;" class="">==============================</div><div style="color: rgb(0, 0, 0); font-family: Helvetica; font-size: 12px; font-style: normal; font-variant-ligatures: normal; font-variant-caps: normal; font-variant-east-asian: normal; font-variant-position: normal; font-weight: normal; letter-spacing: normal; line-height: normal; text-align: start; text-indent: 0px; text-transform: none; white-space: normal; word-spacing: 0px; -webkit-text-stroke-width: 0px; orphans: 2; widows: 2;" class="">Thomas D. Kühne</div><div style="color: rgb(0, 0, 0); font-family: Helvetica; font-size: 12px; font-style: normal; font-variant-ligatures: normal; font-variant-caps: normal; font-variant-east-asian: normal; font-variant-position: normal; font-weight: normal; letter-spacing: normal; line-height: normal; text-align: start; text-indent: 0px; text-transform: none; white-space: normal; word-spacing: 0px; -webkit-text-stroke-width: 0px; orphans: 2; widows: 2;" class="">Dynamics of Condensed Matter</div><div style="color: rgb(0, 0, 0); font-family: Helvetica; font-size: 12px; font-style: normal; font-variant-ligatures: normal; font-variant-caps: normal; font-variant-east-asian: normal; font-variant-position: normal; font-weight: normal; letter-spacing: normal; line-height: normal; text-align: start; text-indent: 0px; text-transform: none; white-space: normal; word-spacing: 0px; -webkit-text-stroke-width: 0px; orphans: 2; widows: 2;" class="">Chair of Theoretical Chemistry</div><div style="color: rgb(0, 0, 0); font-family: Helvetica; font-size: 12px; font-style: normal; font-variant-ligatures: normal; font-variant-caps: normal; font-variant-east-asian: normal; font-variant-position: normal; font-weight: normal; letter-spacing: normal; line-height: normal; text-align: start; text-indent: 0px; text-transform: none; white-space: normal; word-spacing: 0px; -webkit-text-stroke-width: 0px; orphans: 2; widows: 2;" class="">University of Paderborn</div><div style="color: rgb(0, 0, 0); font-family: Helvetica; font-size: 12px; font-style: normal; font-variant-ligatures: normal; font-variant-caps: normal; font-variant-east-asian: normal; font-variant-position: normal; font-weight: normal; letter-spacing: normal; line-height: normal; text-align: start; text-indent: 0px; text-transform: none; white-space: normal; word-spacing: 0px; -webkit-text-stroke-width: 0px; orphans: 2; widows: 2;" class="">Warburger Str. 100</div><div style="color: rgb(0, 0, 0); font-family: Helvetica; font-size: 12px; font-style: normal; font-variant-ligatures: normal; font-variant-caps: normal; font-variant-east-asian: normal; font-variant-position: normal; font-weight: normal; letter-spacing: normal; line-height: normal; text-align: start; text-indent: 0px; text-transform: none; white-space: normal; word-spacing: 0px; -webkit-text-stroke-width: 0px; orphans: 2; widows: 2;" class="">D-33098 Paderborn</div><div style="color: rgb(0, 0, 0); font-family: Helvetica; font-size: 12px; font-style: normal; font-variant-ligatures: normal; font-variant-caps: normal; font-variant-east-asian: normal; font-variant-position: normal; font-weight: normal; letter-spacing: normal; line-height: normal; text-align: start; text-indent: 0px; text-transform: none; white-space: normal; word-spacing: 0px; -webkit-text-stroke-width: 0px; orphans: 2; widows: 2;" class="">Germany</div><div style="color: rgb(0, 0, 0); font-family: Helvetica; font-size: 12px; font-style: normal; font-variant-ligatures: normal; font-variant-caps: normal; font-variant-east-asian: normal; font-variant-position: normal; font-weight: normal; letter-spacing: normal; line-height: normal; text-align: start; text-indent: 0px; text-transform: none; white-space: normal; word-spacing: 0px; -webkit-text-stroke-width: 0px; orphans: 2; widows: 2;" class=""><a href="mailto:tdku...@mail.upb.de" class="">tdku...@mail.upb.de</a></div><div style="color: rgb(0, 0, 0); font-family: Helvetica; font-size: 12px; font-style: normal; font-variant-ligatures: normal; font-variant-caps: normal; font-variant-east-asian: normal; font-variant-position: normal; font-weight: normal; letter-spacing: normal; line-height: normal; text-align: start; text-indent: 0px; text-transform: none; white-space: normal; word-spacing: 0px; -webkit-text-stroke-width: 0px; orphans: 2; widows: 2;" class="">+49/(0)5251/60-5726</div></div></div></div>
</div>
<br class=""></body></html>