<div dir="ltr">Hi Chris,<div><br></div><div>This is not a CP2K-specific question. It is about how you would like to define the BSSE-corrected interaction energy between two "monomers". In the original counterpoise scheme by Boys and Bernardi, the geometries of the two "monomers" were fully optimised with the "dimer" basis. In CP2K, you will have to set up some GHOST atoms manually in order to relax the monomer geometries with "dimer" basis. One approximation you can make is that you calculate the interaction energy at fixed monomer geometries in the "dimer", i.e. you neglect the geometry relaxations of both "monomers" due to additional basis sets on GHOST atoms, and this is exactly what CP2K is doing. In some other quantum chemistry codes (e.g. Gaussian), the geometry relaxations of monomers may be automatically taken into account if you switch on the counterpoise correction scheme.</div><div><br></div><div>BTW: You didn't break the symmetry in the electronic configuration of Cu ions, and currently you have a closed-shell singlet solution. Are you sure that's the electronic ground state?</div><div><br></div><div>SL</div><div><br></div><div><br></div></div><div class="gmail_extra"><br><div class="gmail_quote">On 16 May 2016 at 13:15, Chris Campbell <span dir="ltr"><<a href="mailto:chris...@gmail.com" target="_blank">chris...@gmail.com</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div dir="ltr">Afternoon SL,<br><br>Apologies for the hassle but wondered if I could quickly ask a follow up question.<br><br>Initially for my system using a PBE functional I optimised my Adsorbent and Adsorbate independently and then together, when i found the difference I got an interaction energy of ~25kJ/mol.<br><br>Following this BSSE fragment Calculation (using the fully optimised combination of Adsorbent and Adsorbate together), I get a value in the output of BSSE-free interaction energy: ~-20kJ/mol, which seems very reasonable (to me).<br><br>This energy, however , calculated for the Adsorbate and adsorbent is not from the independent optimisation positions so can I use this energy? or do I have to work backwards to calculate the BSSE correction and then apply it to the original -25kJ/mol?<br><br>Apologies if this is not worded very clearly.<br><br>Kindest Regards,<br>Chris<br></div><div class="HOEnZb"><div class="h5">

<p></p>

-- <br>
You received this message because you are subscribed to the Google Groups "cp2k" group.<br>
To unsubscribe from this group and stop receiving emails from it, send an email to <a href="mailto:cp2k+uns...@googlegroups.com" target="_blank">cp2k+uns...@googlegroups.com</a>.<br>
To post to this group, send email to <a href="mailto:cp...@googlegroups.com" target="_blank">cp...@googlegroups.com</a>.<br>
Visit this group at <a href="https://groups.google.com/group/cp2k" target="_blank">https://groups.google.com/group/cp2k</a>.<br>
For more options, visit <a href="https://groups.google.com/d/optout" target="_blank">https://groups.google.com/d/optout</a>.<br>
</div></div></blockquote></div><br></div>