<div>Hello, <br /></div><div>Could you please send the .inp files so that I can have a look ? It looks like maybe something went wrong with the geometry optimization of either your lone components or your adduct.</div><div>Best,</div><div>Quentin</div><br /><div class="gmail_quote"><div dir="auto" class="gmail_attr">Le mercredi 26 février 2025 à 13:22:13 UTC+1, Oksana Grinevich a écrit :<br/></div><blockquote class="gmail_quote" style="margin: 0 0 0 0.8ex; border-left: 1px solid rgb(204, 204, 204); padding-left: 1ex;">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></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 visit <a href="https://groups.google.com/d/msgid/cp2k/d5cefdbe-b5c2-4935-a389-e0b9f3c1ff93n%40googlegroups.com?utm_medium=email&utm_source=footer">https://groups.google.com/d/msgid/cp2k/d5cefdbe-b5c2-4935-a389-e0b9f3c1ff93n%40googlegroups.com</a>.<br />