<html><head><meta http-equiv="content-type" content="text/html; charset=us-ascii"></head><body style="overflow-wrap: break-word; -webkit-nbsp-mode: space; line-break: after-white-space;">Xavier,<div><br></div><div>Thank you for the reference and the comments.  This is very helpful indeed.</div><div><br></div><div>Best,</div><div>Eric<br><div><br><blockquote type="cite"><div>On Nov 4, 2022, at 2:05 PM, Xavier Bidault <jazzquark@gmail.com> wrote:</div><br class="Apple-interchange-newline"><div><div dir="ltr">Hi Eric,<div><br></div><div>I would recommend this paper: <a href="https://doi.org/10.1063/5.0041511">https://doi.org/10.1063/5.0041511</a></div><div>where we used CP2K/PHONOPY for phonon calculations of molecular crystals.</div><div>We used a large cutoff of 1500 Ry for CP2K, but you can use 600 Ry, which is already tightly converged. But DO NOT use lower than 600 Ry for CHON systems. You can do your own convergence tests, of course.</div><div>REL_CUTOFF is not that important (60 Ry or default is good).</div><div>But NGRIDS 5 is important! Faster calculations.<br></div><div>The driving convergence parameter is MAX_FORCE 1E-5. You can keep the others at 1E-3. </div><div>Yes, EPS_SCF 1E-8.</div><div>Yes, NLCC.</div><div>I think SR is not bad, but do some comparison (optimization and phonon calculations) with non-SR to see if SR is accurate enough.</div><div>PBE-D2 or PBE-D3(BJ)? Look at my concerns in the paper. We ended up with D2 with the pseudo potentials, but D3(BJ) with all-electron calculations.<br></div><div>GAPW is for all-electron calculations.</div><div><br></div><div>Best,</div><div>Xavier</div></div><br><div class="gmail_quote"><div dir="ltr" class="gmail_attr">On Fri, Nov 4, 2022 at 7:29 AM Eric Patterson <<a href="mailto:eric.v.patterson@gmail.com">eric.v.patterson@gmail.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">Dear all,<div><br></div><div>I am soon to embark on a series of optimizations and harmonic frequency calculations for a set of molecular solids (mostly C, H, N, O organics).  I have been combing this forum for advice, and would like to check that the following are general best practices for these calculations.  Any advice and corrections greatly appreciated!</div><div><br></div><div>I plan to use PBE-D3(BJ)/DZVP-MOLOPT-SR-GTH.</div><div><br></div><div>I believe I understand how to select CUTOFF and REL_CUTOFF, and plan to use EPS_SCF 1E-8.</div><div><br></div><div>It seems best to use NLCC pseudopotentials with GPW, correct?</div><div>Also use smoothing?  If so, which is recommended?</div><div>Use USE_FINER_GRID, or not?</div><div><br></div><div>If using GAPW, still use NLCC?</div><div>Same questions regarding smoothing and USE_FINER_GRID.</div><div><br></div><div>I have read that extremely tight convergence on the geometry is recommended.  Do the following seem appropriate?</div><div><br></div><div>    MAX_FORCE 4.5E-6<br>    RMS_FORCE 3.0E-6<br>    MAX_DR 3.0E-5<br>    RMS_DR 1.5E-5<br></div><div><br></div><div>Anything else experienced users might wish to share?  I appreciate any advice as it could save me a ton of time benchmarking various options.</div><div><br></div><div>Thank you,</div><div>Eric</div><div><br class="webkit-block-placeholder"></div>

-- <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:cp2k+unsubscribe@googlegroups.com" target="_blank">cp2k+unsubscribe@googlegroups.com</a>.<br>
To view this discussion on the web visit <a href="https://groups.google.com/d/msgid/cp2k/d527e09a-28eb-41d3-b751-2e3ae555b0e3n%40googlegroups.com?utm_medium=email&utm_source=footer" target="_blank">https://groups.google.com/d/msgid/cp2k/d527e09a-28eb-41d3-b751-2e3ae555b0e3n%40googlegroups.com</a>.<br>
</blockquote></div><div><br class="webkit-block-placeholder"></div>

-- <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:cp2k+unsubscribe@googlegroups.com">cp2k+unsubscribe@googlegroups.com</a>.<br>
To view this discussion on the web visit <a href="https://groups.google.com/d/msgid/cp2k/CAD0N%2BNVDhq8WCzP9qg418ZSGRFiJwiNsrHrRC8_6gr%3DW67Cz2Q%40mail.gmail.com?utm_medium=email&utm_source=footer">https://groups.google.com/d/msgid/cp2k/CAD0N%2BNVDhq8WCzP9qg418ZSGRFiJwiNsrHrRC8_6gr%3DW67Cz2Q%40mail.gmail.com</a>.<br>
</div></blockquote></div><br></div></body></html>

<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:cp2k+unsubscribe@googlegroups.com">cp2k+unsubscribe@googlegroups.com</a>.<br />
To view this discussion on the web visit <a href="https://groups.google.com/d/msgid/cp2k/0C2705DD-488E-4CF4-8D0A-7A35C422E8BB%40gmail.com?utm_medium=email&utm_source=footer">https://groups.google.com/d/msgid/cp2k/0C2705DD-488E-4CF4-8D0A-7A35C422E8BB%40gmail.com</a>.<br />