<div dir="ltr">Dear all,<div>the linked paper, from 2010, (http://aip.scitation.org/doi/10.1063/1.3491271), 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><style type="text/css">
p.p1 {margin: 0.0px 0.0px 0.0px 0.0px; font: 10.0px Times}
span.s1 {font: 8.0px Times}
</style></div>