<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>