[gmx-users] Computing the "energy of interaction"
golshan.hejazi at yahoo.com
Tue Dec 10 18:26:54 CET 2013
I have a paracetamol crystal and I would like to estimate the paracetamol/paracetamol energy of interaction per unitcell.
To do this I am doing the following:
Energy of interaction of paracetamols per unitcell=
Total energy of a unitcell - (Energy of a molecule in the unitcell* number of the molecules in the unitcell)
Energy of interaction of paracetamols per unitcell= Term1 + Term2
How to compute Term1:
Took a crystal slab. Compute the total energy(using steepest decent) and divide it by the number of unitcells. The results is: -2443.29 kj/mol
How to compute Term2:
number of the molecules in the unitcell is 4. Energy of a single molecule in vacuum is 490.20 kj/mol
Energy of interaction of paracetamols per unitcell= -2443.29 kj/mol-(-490.20*4) kj/mol= -482.378 kj/mol
Which is more than a covalent bond!
Could you please help me to figure out what is going wrong?
More information about the gromacs.org_gmx-users