<div dir="ltr">Thanks for the suggestion, Frederick.<div><br></div><div>I exchanged the XC section with the one from `regtest-hybrid-3/H2O-hybrid-b3lyp_shortcut.inp`, so</div><div><br></div><div><div class="prettyprint" style="background-color: rgb(250, 250, 250); border-color: rgb(187, 187, 187); border-style: solid; border-width: 1px; overflow-wrap: break-word;"><code class="prettyprint"><div class="subprettyprint"><div class="subprettyprint">    &XC</div><div class="subprettyprint">      &XC_FUNCTIONAL B3LYP</div><div class="subprettyprint">      &END XC_FUNCTIONAL</div><div class="subprettyprint">      &HF</div><div class="subprettyprint">        &SCREENING</div><div class="subprettyprint">          EPS_SCHWARZ 1.0E-10</div><div class="subprettyprint">        &END</div><div class="subprettyprint">        &MEMORY</div><div class="subprettyprint">          MAX_MEMORY  5</div><div class="subprettyprint">        &END</div><div class="subprettyprint">      &END</div><div class="subprettyprint">    &END XC</div></div></code></div><br>But unfortunately it's still not converging. I attached the input & output of this run.</div><div><br></div><div>(By the way, I was previously using B3LYP which is also used in the test `regtest-libxc/H2O-hybrid-b3lyp_libxc.inp`)<br><br><br>On Monday, October 1, 2018 at 6:46:44 PM UTC+2, Frederick Stein wrote:<blockquote class="gmail_quote" style="margin: 0;margin-left: 0.8ex;border-left: 1px #ccc solid;padding-left: 1ex;"><p dir="ltr">Dear Kristjan,</p>
<p dir="ltr">You have to delete the line with the keyword FRACTION in the HF-section. And I would use the B3LYP shortcut. There is an example file in regtest-hybrids-3, I think.</p>
<p dir="ltr">Cheers.<br>
Frederick</p>
<br><div class="gmail_quote"><div dir="ltr">Kristjan Eimre <<a href="javascript:" target="_blank" gdf-obfuscated-mailto="QjUXpRk4AgAJ" rel="nofollow" onmousedown="this.href='javascript:';return true;" onclick="this.href='javascript:';return true;">kris...@gmail.com</a>> schrieb am Mo., 01.10.2018, 17:41:<br></div><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div dir="ltr">Dear CP2K community,<div><br></div><div>I am having trouble converging the SCF of a b3lyp/6-311++G** calculation for a fairly small molecule consisting of C and H.</div><div><br></div><div>Without too much trouble, I was able to converge the b3lyp/6-311++G** for benzene and also b3lyp/6-311G** for my molecule (so without the diffuse basis functions).</div><div><br></div><div>I optimized the geometries also with b3lyp/6-311G**.</div><div><br></div><div>Additionally, for my molecule, I converged blyp/6-311++G**, which I was using as the initial guess to b3lyp calculation.</div><div><br></div><div>I attached one of the inputs scripts that I tried and the corresponding output.</div><div><br></div><div>Here are the different things that I tried:</div><div>* Played with the parameters in QS (most notably EPSFIT, ALPHA0_H)</div><div>* OT with CG minimizer</div><div>* Diagonalization</div><div>* MT poisson solver</div><div>* UKS calculation<br></div><div>* smaller basis set, 6-311+G* (from EMSL)</div><div><br></div><div>But unfortunately everything failed.</div><div><br></div><div>Which knobs should I turn to possibly reach convergence? Any hints are very much appreciated.</div><div><br></div><div>best regards,<br>Kristjan</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:" rel="nofollow" target="_blank" gdf-obfuscated-mailto="QjUXpRk4AgAJ" onmousedown="this.href='javascript:';return true;" onclick="this.href='javascript:';return true;">cp2k+uns...@googlegroups.<wbr>com</a>.<br>
To post to this group, send email to <a href="javascript:" rel="nofollow" target="_blank" gdf-obfuscated-mailto="QjUXpRk4AgAJ" onmousedown="this.href='javascript:';return true;" onclick="this.href='javascript:';return true;">cp...@googlegroups.com</a>.<br>
Visit this group at <a href="https://groups.google.com/group/cp2k" rel="nofollow" target="_blank" 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>
For more options, visit <a href="https://groups.google.com/d/optout" rel="nofollow" target="_blank" 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></div>