<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">https://www.cp2k.org/howto: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 dependant 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>