<div dir="ltr">Correction: You should calculate the BSSE-corrected binding energy with both "monomers" in the dimer geometry (this is what &BSSE does in CP2K), not from independently optimised "monomer" geometry (your first calculation), especially if the host-guest interaction is strong. See this document for more discussion, http://vergil.chemistry.gatech.edu/notes/cp.pdf<div><div><br></div><div>As I pointed out, the electronic configuration of Cu2+ ions will affect the host-guest interaction in your case, so you will need to set up the electronic configurations correctly.<div><br></div><div>SL</div><div><br><br>On Monday, May 16, 2016 at 5:06:34 PM UTC+1, S Ling wrote:<blockquote class="gmail_quote" style="margin: 0;margin-left: 0.8ex;border-left: 1px #ccc solid;padding-left: 1ex;"><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><br><div class="gmail_quote">On 16 May 2016 at 13:15, Chris Campbell <span dir="ltr"><<a href="javascript:" target="_blank" gdf-obfuscated-mailto="Aw3kUCoZBgAJ" rel="nofollow" onmousedown="this.href='javascript:';return true;" onclick="this.href='javascript:';return true;">chri...@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><div>

<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="javascript:" target="_blank" gdf-obfuscated-mailto="Aw3kUCoZBgAJ" rel="nofollow" onmousedown="this.href='javascript:';return true;" onclick="this.href='javascript:';return true;">cp2k+uns...@googlegroups.<wbr>com</a>.<br>
To post to this group, send email to <a href="javascript:" target="_blank" gdf-obfuscated-mailto="Aw3kUCoZBgAJ" rel="nofollow" onmousedown="this.href='javascript:';return true;" onclick="this.href='javascript:';return true;">cp...@googlegroups.com</a>.<br>
Visit this group at <a href="https://groups.google.com/group/cp2k" target="_blank" rel="nofollow" onmousedown="this.href='https://groups.google.com/group/cp2k';return true;" onclick="this.href='https://groups.google.com/group/cp2k';return true;">https://groups.google.com/<wbr>group/cp2k</a>.<br>
For more options, visit <a href="https://groups.google.com/d/optout" target="_blank" rel="nofollow" onmousedown="this.href='https://groups.google.com/d/optout';return true;" onclick="this.href='https://groups.google.com/d/optout';return true;">https://groups.google.com/d/<wbr>optout</a>.<br>
</div></div></blockquote></div><br></div>
</blockquote></div></div></div></div>