<div dir="auto"><div dir="auto">Dear Tamas Karpati,</div><div dir="auto"><br></div><div dir="auto">    Yes definitely sorry for the unclear original question. The E(a+b) - [E(a) +E(b)] quantity is the one of interest here. My concern though is that the reference for total dft energies here matters as it crops up in the finals result for the formula above, even if all references are the same. It seems from trying several things that this is ends up being problematic in the case of cells with different charge specified.</div><div dir="auto"><br></div><div dir="auto">   I know the QE faqs say reference levels are not well defined when using PBCs, but is there any way to try to extrapolate to vacuum (all cells are same size with a lot of vacuum in all directions) or am I maybe missing something in the approach and something else is the problem? Thanks again for the information!</div><div dir="auto"><br></div><div dir="auto">Regards,</div><div dir="auto">Robert Stanton</div><div dir="auto">Graduate Student</div><div dir="auto">Clarkson University </div></div><div><br><div class="gmail_quote"><div dir="ltr" class="gmail_attr">On Tue, Dec 29, 2020 at 3:16 AM Tamas Karpati <<a href="mailto:tkarpati@gmail.com" target="_blank" rel="noreferrer">tkarpati@gmail.com</a>> wrote:<br></div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left-width:1px;border-left-style:solid;padding-left:1ex;border-left-color:rgb(204,204,204)">Dear Robert,<br>
<br>
Thanks for clarification. If I understand correctly,<br>
you intend to model A + B interaction energies using<br>
E(A) and E(B) values obtained for the separate systems.<br>
<br>
Separate means non-interacting so you lost your<br>
interaction energy in such simulations.<br>
At the atomic scale the situation is more complex<br>
than it is for eg. point charges (think of the distributed<br>
electron and shielded core charges at small scales).<br>
>From this point of view, your current problem is also<br>
different to the separation of vdW parameters into<br>
atomic (~ isolated) parameters where they model<br>
the interaction of point-like objects. Yours is not one.<br>
<br>
Most probably you need the association energy<br>
     E(A+B) - [ E(A) + E(B) ]<br>
evaluated for carefully selected cells and params.<br>
<br>
(Sorry for the inconvenience: last time E meant a cell<br>
or material, here it denotes total energy.)<br>
<br>
Bests,<br>
  t<br>
<br>
On Tue, Dec 29, 2020 at 1:33 AM Robert Stanton <<a href="mailto:stantor@clarkson.edu" target="_blank" rel="noreferrer">stantor@clarkson.edu</a>> wrote:<br>
><br>
> Dear Tamas Karpati,<br>
><br>
>      Definitely, the question is basically what association energies are for things which are bound ionically. It seems I require calculation of the lone ions with compensating background charge, as well as the neutral compound. I am just wondering how I can partially correct/change what these absolute DFT energies are with respect to for valid comparison between calculations. Is there a way to do this or do you think even if achieved, this still might be a problematic approach? Thanks in advance!<br>
><br>
> Regards,<br>
> Robert Stanton<br>
><br>
> On Mon, Dec 28, 2020 at 12:20 PM Tamas Karpati <<a href="mailto:tkarpati@gmail.com" target="_blank" rel="noreferrer">tkarpati@gmail.com</a>> wrote:<br>
>><br>
>> Dear Robert,<br>
>><br>
>> Since each member of the series E(2-) ~ E(1-) ~ E ~ E(1+) ~ E(2+)<br>
>> etc... dwell on a different<br>
>> potential energy hypersurface, there is no direct physical meaning of<br>
>> their energy differences.<br>
>> Can you give the concrete question you try to answer?<br>
>><br>
>> On the other hand, it can be quite useful to compare eg. the<br>
>> stabilization energies (or other<br>
>> derived energy-like quantities, even band positions). In case of the<br>
>> stab. ens. it is common to<br>
>> make hypothetical reactions and evaluate their energies which can then<br>
>> directly be compared.<br>
>><br>
>> Hope this helps,<br>
>>   t<br>
>><br>
>><br>
>> On Mon, Dec 28, 2020 at 7:23 PM Robert Stanton <<a href="mailto:stantor@clarkson.edu" target="_blank" rel="noreferrer">stantor@clarkson.edu</a>> wrote:<br>
>> ><br>
>> > Dear all,<br>
>> ><br>
>> >     I am just wondering if there is, in general, an approach to comparing energies between two cells with different tot charge specifications? Is there a way for example to extrapolate the energy values to be w.r.t. vacuum to have a common comparison point?<br>
>> ><br>
>> >     I have seen somewhat similar questions, generally these are specifically for 2D system, where it's stated there are examples addressing this, but have not been able to find those examples in the QE install. Any information, or a pointer to this specific example (if it applies to more than just 2D) would be greatly appreciated! Thanks in advance for any help.<br>
>> ><br>
>> > Regards,<br>
>> > Robert Stanton<br>
>> > Graduate Student<br>
>> > Clarkson University<br>
>> > _______________________________________________<br>
>> > Quantum ESPRESSO is supported by MaX (<a href="http://www.max-centre.eu" rel="noreferrer noreferrer" target="_blank">www.max-centre.eu</a>)<br>
>> > users mailing list <a href="mailto:users@lists.quantum-espresso.org" target="_blank" rel="noreferrer">users@lists.quantum-espresso.org</a><br>
>> > <a href="https://lists.quantum-espresso.org/mailman/listinfo/users" rel="noreferrer noreferrer" target="_blank">https://lists.quantum-espresso.org/mailman/listinfo/users</a><br>
>> _______________________________________________<br>
>> Quantum ESPRESSO is supported by MaX (<a href="http://www.max-centre.eu" rel="noreferrer noreferrer" target="_blank">www.max-centre.eu</a>)<br>
>> users mailing list <a href="mailto:users@lists.quantum-espresso.org" target="_blank" rel="noreferrer">users@lists.quantum-espresso.org</a><br>
>> <a href="https://lists.quantum-espresso.org/mailman/listinfo/users" rel="noreferrer noreferrer" target="_blank">https://lists.quantum-espresso.org/mailman/listinfo/users</a><br>
><br>
> _______________________________________________<br>
> Quantum ESPRESSO is supported by MaX (<a href="http://www.max-centre.eu" rel="noreferrer noreferrer" target="_blank">www.max-centre.eu</a>)<br>
> users mailing list <a href="mailto:users@lists.quantum-espresso.org" target="_blank" rel="noreferrer">users@lists.quantum-espresso.org</a><br>
> <a href="https://lists.quantum-espresso.org/mailman/listinfo/users" rel="noreferrer noreferrer" target="_blank">https://lists.quantum-espresso.org/mailman/listinfo/users</a><br>
_______________________________________________<br>
Quantum ESPRESSO is supported by MaX (<a href="http://www.max-centre.eu" rel="noreferrer noreferrer" target="_blank">www.max-centre.eu</a>)<br>
users mailing list <a href="mailto:users@lists.quantum-espresso.org" target="_blank" rel="noreferrer">users@lists.quantum-espresso.org</a><br>
<a href="https://lists.quantum-espresso.org/mailman/listinfo/users" rel="noreferrer noreferrer" target="_blank">https://lists.quantum-espresso.org/mailman/listinfo/users</a><br>
</blockquote></div></div>