<div>Dear members:</div><div>good evening/morning.</div><div><br></div><div>it is not clear to me how to interpret CP2K print out for BSSE.</div><div><br></div><div>If I take a small guest (like benzene) next to a carbon layer using GWP/PBE/DZVP</div><div>Of course, since interactions are so week, DFTD3 is a must - otherwise a benzene may "fly away".</div><div><br></div><div>Then, using BSSE  for S1(guest) and S2 (carbon layer)</div><div>CP2K prints<br></div><div>S1 + S2ghost: Total energy -455.12112693238430 Hartree<br>S1ghost + S2: Total energy -455.12112693238430 Hartree<br>S1           + S2: Total energy -455.12112693238441 Hartree<br></div><div><br></div><div></div><div>The differences are none or tiny.<br></div><div></div><div><br></div><div>However, the final print is:</div>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<div><br></div><div><b>First, it is not clear what is "CP".</b></div><div></div><div><br></div><div>Next, i<b>f this is the value of the BSSE artifact</b> (due to the wave-function truncation), <br></div><div>then <b>its huge</b> 41.2964895 kJ/mole.</div><div>This would be <b>almost twice larger than the binding energy</b> of benzene-like guests when at carbonaceous surfaces.</div><div></div><div></div><div>This would mean that quality evaluations of binding energies for such systems using CP2K would be very difficult. <br></div><div><br></div><div>Any comments to hear would be nice.<br></div><div><br></div><div>Finally, it is not clear why 
2-body contribution and<b> BSSE-free interaction are printed as positive numbers</b>.</div><div><br></div><div>I hope the questions I ask be interesting for users and helpful for the on-line manual to update.</div><div><br></div><div>After all, I wish to note that it was by a <b>rare chance I was able to track</b> the old <br></div><div>notice <b>on the importance of declaring</b><br></div><div>    &KIND X_ghost<br>     BASIS_SET DZVP-MOLOPT-SR-GTH<br>     GHOST<br>    &END KIND </div><div><br></div><div>Even though ... <br></div><div><b>if I skip such definition, the numbers would come "a bit better"</b></div><div><b>:-)<br></b></div><div>CP would be about 21 kJ/mole, what is about comparable to the binding energy <br></div><div>of a benzene-like guest to a carbonaceous material. <br></div><div> <br></div><div>Thank you for your time and attention.</div><div><br></div><div>I wish everybody be well and happy.</div><div><br></div><div>Victor<br></div><div><br></div><div> <br></div><div><br></div><div><br></div><div><br></div><div><br></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/44686fb2-400e-443b-956d-e5abefa70cebn%40googlegroups.com?utm_medium=email&utm_source=footer">https://groups.google.com/d/msgid/cp2k/44686fb2-400e-443b-956d-e5abefa70cebn%40googlegroups.com</a>.<br />