<div dir="ltr"><div>Dear Matthias:</div><div>thank you for the kind reply.</div><div>I did not wish to be "too specific" about the system in public because the results are under <br></div><div>a peer review.</div><div>Here I attach one of the exemplary cases: pyridazine at graphene. <br></div><div>Thank you.</div><div>With best regards,</div><div>Victor<br></div><div><br></div> </div><br><div class="gmail_quote"><div dir="ltr" class="gmail_attr">On Wed, Oct 19, 2022 at 12:12 PM Krack Matthias (PSI) <<a href="mailto:matthias.krack@psi.ch">matthias.krack@psi.ch</a>> wrote:<br></div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><div class="msg-5433807862262957151"> <div style="overflow-wrap: break-word;" lang="en-CH"> <div class="m_-5433807862262957151WordSection1"> <p class="MsoNormal"><span style="font-size:11pt" lang="DE-CH">Hi Victor<u></u><u></u></span></p> <p class="MsoNormal"><span style="font-size:11pt" lang="DE-CH"><u></u> <u></u></span></p> <p class="MsoNormal"><span style="font-size:11pt" lang="EN-US">Without the full input, it is difficult to make any suggestion. Do you have a &BSSE section defining the FRAGMENTS?<u></u><u></u></span></p> <p class="MsoNormal"><span style="font-size:11pt" lang="EN-US"><u></u> <u></u></span></p> <p class="MsoNormal"><span style="font-size:11pt" lang="EN-US">Best<u></u><u></u></span></p> <p class="MsoNormal"><span style="font-size:11pt" lang="EN-US"><u></u> <u></u></span></p> <p class="MsoNormal"><span style="font-size:11pt" lang="EN-US">Matthias<u></u><u></u></span></p> <p class="MsoNormal"><span style="font-size:11pt" lang="EN-US"><u></u> <u></u></span></p> <div style="border-color:rgb(181,196,223) currentcolor currentcolor;border-style:solid none none;border-width:1pt medium medium;padding:3pt 0cm 0cm"> <p class="MsoNormal" style="margin-right:0cm;margin-bottom:12pt;margin-left:36pt"> <b><span style="font-size:12pt;color:black">From: </span></b><span style="font-size:12pt;color:black"><a href="mailto:cp2k@googlegroups.com" target="_blank">cp2k@googlegroups.com</a> <<a href="mailto:cp2k@googlegroups.com" target="_blank">cp2k@googlegroups.com</a>> on behalf of Victor Volkov <<a href="mailto:volkovskr@gmail.com" target="_blank">volkovskr@gmail.com</a>><br> <b>Date: </b>Saturday, 15 October 2022 at 19:53<br> <b>To: </b>cp2k <<a href="mailto:cp2k@googlegroups.com" target="_blank">cp2k@googlegroups.com</a>><br> <b>Subject: </b>[CP2K:17873] understanding BSSE print out<u></u><u></u></span></p> </div> <div> <p class="MsoNormal" style="margin-left:36pt"><span style="font-size:11pt">Dear members:<u></u><u></u></span></p> </div> <div> <p class="MsoNormal" style="margin-left:36pt"><span style="font-size:11pt">good evening/morning.<u></u><u></u></span></p> </div> <div> <p class="MsoNormal" style="margin-left:36pt"><span style="font-size:11pt"><u></u> <u></u></span></p> </div> <div> <p class="MsoNormal" style="margin-left:36pt"><span style="font-size:11pt">it is not clear to me how to interpret CP2K print out for BSSE.<u></u><u></u></span></p> </div> <div> <p class="MsoNormal" style="margin-left:36pt"><span style="font-size:11pt"><u></u> <u></u></span></p> </div> <div> <p class="MsoNormal" style="margin-left:36pt"><span style="font-size:11pt">If I take a small guest (like benzene) next to a carbon layer using GWP/PBE/DZVP<u></u><u></u></span></p> </div> <div> <p class="MsoNormal" style="margin-left:36pt"><span style="font-size:11pt">Of course, since interactions are so week, DFTD3 is a must - otherwise a benzene may "fly away".<u></u><u></u></span></p> </div> <div> <p class="MsoNormal" style="margin-left:36pt"><span style="font-size:11pt"><u></u> <u></u></span></p> </div> <div> <p class="MsoNormal" style="margin-left:36pt"><span style="font-size:11pt">Then, using BSSE for S1(guest) and S2 (carbon layer)<u></u><u></u></span></p> </div> <div> <p class="MsoNormal" style="margin-left:36pt"><span style="font-size:11pt">CP2K prints<u></u><u></u></span></p> </div> <div> <p class="MsoNormal" style="margin-left:36pt"><span style="font-size:11pt">S1 + S2ghost: Total energy -455.12112693238430 Hartree<br> S1ghost + S2: Total energy -455.12112693238430 Hartree<br> S1 + S2: Total energy -455.12112693238441 Hartree<u></u><u></u></span></p> </div> <div> <p class="MsoNormal" style="margin-left:36pt"><span style="font-size:11pt"><u></u> <u></u></span></p> </div> <div> <p class="MsoNormal" style="margin-left:36pt"><span style="font-size:11pt">The differences are none or tiny.<u></u><u></u></span></p> </div> <div> <p class="MsoNormal" style="margin-left:36pt"><span style="font-size:11pt"><u></u> <u></u></span></p> </div> <div> <p class="MsoNormal" style="margin-left:36pt"><span style="font-size:11pt">However, the final print is:<u></u><u></u></span></p> </div> <p class="MsoNormal" style="margin-left:36pt"><span style="font-size:11pt">CP-corrected Total energy: 0.015729 <br> 1-body contribution: -44.96634268204556 <br> 1-body contribution: -410.13905550190344<br> 2-body contribution: 455.121127 <br> BSSE-free interaction energy: 455.121127<u></u><u></u></span></p> <div> <p class="MsoNormal" style="margin-left:36pt"><span style="font-size:11pt"><u></u> <u></u></span></p> </div> <div> <p class="MsoNormal" style="margin-left:36pt"><b><span style="font-size:11pt">First, it is not clear what is "CP".</span></b><span style="font-size:11pt"><u></u><u></u></span></p> </div> <div> <p class="MsoNormal" style="margin-left:36pt"><span style="font-size:11pt"><u></u> <u></u></span></p> </div> <div> <p class="MsoNormal" style="margin-left:36pt"><span style="font-size:11pt">Next, i<b>f this is the value of the BSSE artifact</b> (due to the wave-function truncation), <u></u><u></u></span></p> </div> <div> <p class="MsoNormal" style="margin-left:36pt"><span style="font-size:11pt">then <b>its huge</b> 41.2964895 kJ/mole.<u></u><u></u></span></p> </div> <div> <p class="MsoNormal" style="margin-left:36pt"><span style="font-size:11pt">This would be <b>almost twice larger than the binding energy</b> of benzene-like guests when at carbonaceous surfaces.<u></u><u></u></span></p> </div> <div> <p class="MsoNormal" style="margin-left:36pt"><span style="font-size:11pt">This would mean that quality evaluations of binding energies for such systems using CP2K would be very difficult. <u></u><u></u></span></p> </div> <div> <p class="MsoNormal" style="margin-left:36pt"><span style="font-size:11pt"><u></u> <u></u></span></p> </div> <div> <p class="MsoNormal" style="margin-left:36pt"><span style="font-size:11pt">Any comments to hear would be nice.<u></u><u></u></span></p> </div> <div> <p class="MsoNormal" style="margin-left:36pt"><span style="font-size:11pt"><u></u> <u></u></span></p> </div> <div> <p class="MsoNormal" style="margin-left:36pt"><span style="font-size:11pt">Finally, it is not clear why 2-body contribution and<b> BSSE-free interaction are printed as positive numbers</b>.<u></u><u></u></span></p> </div> <div> <p class="MsoNormal" style="margin-left:36pt"><span style="font-size:11pt"><u></u> <u></u></span></p> </div> <div> <p class="MsoNormal" style="margin-left:36pt"><span style="font-size:11pt">I hope the questions I ask be interesting for users and helpful for the on-line manual to update.<u></u><u></u></span></p> </div> <div> <p class="MsoNormal" style="margin-left:36pt"><span style="font-size:11pt"><u></u> <u></u></span></p> </div> <div> <p class="MsoNormal" style="margin-left:36pt"><span style="font-size:11pt">After all, I wish to note that it was by a <b>rare chance I was able to track</b> the old <u></u><u></u></span></p> </div> <div> <p class="MsoNormal" style="margin-left:36pt"><span style="font-size:11pt">notice <b>on the importance of declaring</b><u></u><u></u></span></p> </div> <div> <p class="MsoNormal" style="margin-left:36pt"><span style="font-size:11pt"> &KIND X_ghost<br> BASIS_SET DZVP-MOLOPT-SR-GTH<br> GHOST<br> &END KIND <u></u><u></u></span></p> </div> <div> <p class="MsoNormal" style="margin-left:36pt"><span style="font-size:11pt"><u></u> <u></u></span></p> </div> <div> <p class="MsoNormal" style="margin-left:36pt"><span style="font-size:11pt">Even though ... <u></u><u></u></span></p> </div> <div> <p class="MsoNormal" style="margin-left:36pt"><b><span style="font-size:11pt">if I skip such definition, the numbers would come "a bit better"</span></b><span style="font-size:11pt"><u></u><u></u></span></p> </div> <div> <p class="MsoNormal" style="margin-left:36pt"><b><span style="font-size:11pt">:-)</span></b><span style="font-size:11pt"><u></u><u></u></span></p> </div> <div> <p class="MsoNormal" style="margin-left:36pt"><span style="font-size:11pt">CP would be about 21 kJ/mole, what is about comparable to the binding energy <u></u><u></u></span></p> </div> <div> <p class="MsoNormal" style="margin-left:36pt"><span style="font-size:11pt">of a benzene-like guest to a carbonaceous material. <u></u><u></u></span></p> </div> <div> <p class="MsoNormal" style="margin-left:36pt"><span style="font-size:11pt"> <u></u><u></u></span></p> </div> <div> <p class="MsoNormal" style="margin-left:36pt"><span style="font-size:11pt">Thank you for your time and attention.<u></u><u></u></span></p> </div> <div> <p class="MsoNormal" style="margin-left:36pt"><span style="font-size:11pt"><u></u> <u></u></span></p> </div> <div> <p class="MsoNormal" style="margin-left:36pt"><span style="font-size:11pt">I wish everybody be well and happy.<u></u><u></u></span></p> </div> <div> <p class="MsoNormal" style="margin-left:36pt"><span style="font-size:11pt"><u></u> <u></u></span></p> </div> <div> <p class="MsoNormal" style="margin-left:36pt"><span style="font-size:11pt">Victor<u></u><u></u></span></p> </div> <div> <p class="MsoNormal" style="margin-left:36pt"><span style="font-size:11pt"><u></u> <u></u></span></p> </div> <div> <p class="MsoNormal" style="margin-left:36pt"><span style="font-size:11pt"> <u></u><u></u></span></p> </div> <div> <p class="MsoNormal" style="margin-left:36pt"><span style="font-size:11pt"><u></u> <u></u></span></p> </div> <div> <p class="MsoNormal" style="margin-left:36pt"><span style="font-size:11pt"><u></u> <u></u></span></p> </div> <div> <p class="MsoNormal" style="margin-left:36pt"><span style="font-size:11pt"><u></u> <u></u></span></p> </div> <div> <p class="MsoNormal" style="margin-left:36pt"><span style="font-size:11pt"><u></u> <u></u></span></p> </div> <p class="MsoNormal" style="margin-left:36pt"><span style="font-size:11pt">-- <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+unsubscribe@googlegroups.com" target="_blank">cp2k+unsubscribe@googlegroups.com</a>.<br> To view this discussion on the web visit <a href="https://groups.google.com/d/msgid/cp2k/44686fb2-400e-443b-956d-e5abefa70cebn%40googlegroups.com?utm_medium=email&utm_source=footer" target="_blank"> https://groups.google.com/d/msgid/cp2k/44686fb2-400e-443b-956d-e5abefa70cebn%40googlegroups.com</a>.<u></u><u></u></span></p> </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="mailto:cp2k+unsubscribe@googlegroups.com" target="_blank">cp2k+unsubscribe@googlegroups.com</a>.<br> To view this discussion on the web visit <a href="https://groups.google.com/d/msgid/cp2k/ZRAP278MB08278F88114F9F43DE58BEDAF42B9%40ZRAP278MB0827.CHEP278.PROD.OUTLOOK.COM?utm_medium=email&utm_source=footer" target="_blank">https://groups.google.com/d/msgid/cp2k/ZRAP278MB08278F88114F9F43DE58BEDAF42B9%40ZRAP278MB0827.CHEP278.PROD.OUTLOOK.COM</a>.<br> </div></blockquote></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="mailto:cp2k+unsubscribe@googlegroups.com">cp2k+unsubscribe@googlegroups.com</a>.<br /> To view this discussion on the web visit <a href="https://groups.google.com/d/msgid/cp2k/CAK%2B-SSgY_70DJ4kwXhs8Df3OUMjURoxm11SssC1P1b-Nj9y%2BRQ%40mail.gmail.com?utm_medium=email&utm_source=footer">https://groups.google.com/d/msgid/cp2k/CAK%2B-SSgY_70DJ4kwXhs8Df3OUMjURoxm11SssC1P1b-Nj9y%2BRQ%40mail.gmail.com</a>.<br />