I am running calculations to obtain binding energy of a molecule on the graphene sheet. After obtaining the energy value of this system (bonded), I calculate the energy values for lone molecule and lone graphene atoms (I tried with the same coordinates and varying them). Thus, I can roughly find the excess (adsorption) energy as follows:<div><br /></div><div>E (excess)=E(bonded)-E(lone molecule)-E(lone graphene)</div><div><br /></div><div>After this step I want to convert this energy to kJ/mol (*2625.5). It worked pretty well for benzene on graphene (obtained 40 kJ/mol) and other aromatic structures but for ammonia and triethylamine I can obtain something like 700 kJ/mol or even more which seems not physically corrected values.</div><div><br /></div><div>I would like to know where I did wrong and what can I do to improve results.</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 visit <a href="https://groups.google.com/d/msgid/cp2k/15f7ade8-1915-40d9-899d-fdb81fdd7d6fn%40googlegroups.com?utm_medium=email&utm_source=footer">https://groups.google.com/d/msgid/cp2k/15f7ade8-1915-40d9-899d-fdb81fdd7d6fn%40googlegroups.com</a>.<br />