<div dir="ltr">Hi,<br><br>The same compensating homogeneous background charge is (implicitly) applied in CP2K if your system is charged. The compensating background charge is required because the electrostatic energy of a periodic charged system diverges. This is standard operation in most codes. <br><br>A direct comparison of energies calculated with different charges is therefore meaningless because the total energy of a charged system converges very slowly w.r.t. to system size. You can however apply a finite-size correction to the energy of the charged system. It should be safe to compare energy differences after the correction. Some codes apply a correction automatically but with CP2K you have to do it manually. Applying e.g. a Madelung correction to the energy should be straightforward. There should be plenty of literature regarding finite-size corrections in charged vacancy modelling and even comparisons of different methods. This <a href="https://journals.aps.org/prb/abstract/10.1103/PhysRevB.86.045112">paper</a> could be a useful starting point.<br><br><br><br>BR, <br><br>Nico<br><br>tiistai 20. kesäkuuta 2017 13.43.56 UTC+3 Cristiana Passiu kirjoitti:<blockquote class="gmail_quote" style="margin: 0;margin-left: 0.8ex;border-left: 1px #ccc solid;padding-left: 1ex;"><div dir="ltr">Dear all,<div>the linked paper, from 2010, (<a href="http://aip.scitation.org/doi/10.1063/1.3491271" target="_blank" rel="nofollow" onmousedown="this.href='http://www.google.com/url?q\x3dhttp%3A%2F%2Faip.scitation.org%2Fdoi%2F10.1063%2F1.3491271\x26sa\x3dD\x26sntz\x3d1\x26usg\x3dAFQjCNH5rSS2OPdxMym8ZULdHuaph--hvA';return true;" onclick="this.href='http://www.google.com/url?q\x3dhttp%3A%2F%2Faip.scitation.org%2Fdoi%2F10.1063%2F1.3491271\x26sa\x3dD\x26sntz\x3d1\x26usg\x3dAFQjCNH5rSS2OPdxMym8ZULdHuaph--hvA';return true;">http://aip.scitation.org/doi/<wbr>10.1063/1.3491271</a>), in which calculations have been performed with the CRYSTAL package, states in the second paragraph: "the total energies of charged systems have no physical meaning in CRYSTAL06 because of the interaction with the balancing background of charge".</div><div>Do you know if this also applies to the cp2k package? I would like to compare the stability of a neutral crystal with those of the same crystal with positive or negative charges, which I'm using to model charged oxygen vacancies. Does it make sense to compare the total energies obtained after geometry optimization?.</div><div><br></div><div>Many thanks,</div><div>Cristiana</div><div><br></div><div><br></div></div></blockquote></div>