<div dir="ltr"><div>Hello All,</div><div><br></div><div>I'm trying to calculate the binding energy of an ion inside a condensed phase unit cell.</div><div><br></div><div>The energy difference equals the energy of the full (neutral) liquid system minus the energy of two oppositely charged systems: vacuum-phase molecule (in this case, negatively charged) and liquid/hole combo (positively charged). </div><div><br></div><div>Simply setting the CHARGE setting in my input file leads to spurious results, especially large dependence on vacuum-phase ion's energy on the unit cell size. </div><div><br></div><div>I've been playing around with turning the liquid unit cell into a cluster, and trying to calculate the energy without periodicity, i.e. using WAVELET poisson with ANALYTIC solver. But I'm not sure whether this plays nice with the CHARGE keyword. </div><div><br></div><div>Maybe a better way would be to stay within periodic boundary conditions, but turn on MULTIPOLE poisson method...again, not sure how this plays with having a charged unit cell. </div><div><br></div><div>I can't find in the documentation what exactly setting a nonzero CHARGE does... does it do nothing, or does it automatically turn on a background jellium?</div><div><br></div><div>Any suggestions would be greatly appreciated. </div><div><br>Thanks,</div><div><br></div><div>L.K.</div><div><br></div><div><br></div><div><br></div></div>