<div dir="ltr">Dear Fabian<div><br></div><div>Thank you for your suggestions, I'm running new tests based on them and hopefully the charge inaccuracy can be solved...</div><div><br></div><div>As for the basis set, I agree that they should not be directly compared since different basis set of course yield different total energy, however I expect somewhat "smaller" difference, not few hundred hartree... With this much of difference I can only assume one of the calclation is completely nonsense...</div><div><br></div><div>Best regards</div><div><br></div><div>Liu Jing</div><div><br>在 2020年4月27日星期一 UTC+2下午5:05:03，Fabian Ducry写道：<blockquote class="gmail_quote" style="margin: 0;margin-left: 0.8ex;border-left: 1px #ccc solid;padding-left: 1ex;"><div dir="ltr"><div><font size="2">Dear Jing,</font></div><div><font size="2"><br></font></div><div><font size="2">Some pointers to improve your simulations:<br></font></div><div><font size="2"> - The charge you find is an artefact of the numerical inaccuracy in your system. EPS_DEFAULT should be tighter, just leave it at the default value of 1.0E-10<br></font><font size="2"> - How did you choose CUTOFF and REL_CUTOFF, did you follow <a href="https://www.cp2k.org/howto:converging_cutoff" target="_blank" rel="nofollow" onmousedown="this.href='https://www.google.com/url?q\x3dhttps%3A%2F%2Fwww.cp2k.org%2Fhowto%3Aconverging_cutoff\x26sa\x3dD\x26sntz\x3d1\x26usg\x3dAFQjCNFHeI16YFpXApR5fGmB8YOPSHi7nQ';return true;" onclick="this.href='https://www.google.com/url?q\x3dhttps%3A%2F%2Fwww.cp2k.org%2Fhowto%3Aconverging_cutoff\x26sa\x3dD\x26sntz\x3d1\x26usg\x3dAFQjCNFHeI16YFpXApR5fGmB8YOPSHi7nQ';return true;">https://www.cp2k.org/howto:<wbr>converging_cutoff</a>? For CELL_OPT you might need to increase these numbers</font></div><div><font size="2"></font></div><div><font size="2"> - SZV is a very small basis set, useful for testing but not for anything else. DZVP often produces okish results but not always. You should check the convergence of the basis set.<br></font></div><div><font size="2"><br></font></div><div><font size="2">The total energy is very much dependent on the basis set that you use. So you cannot compare the energy of your "reference" with the one from cp2k and expect them to be the same.</font></div><div><font size="2"><br></font></div><div><font size="2">Best,</font></div><div><font size="2">Fabian<br></font></div></div></blockquote></div></div>