[gmx-users] Wrong bond and Lennard-Jones potentials in cyclohexane

Mark Abraham Mark.Abraham at anu.edu.au
Tue Jan 5 12:22:21 CET 2010


Henri Ervasti wrote:
> Dear all,
> 
> Thank you for your help and comments on abovementioned problem. In the
> end though it turned out that the large potentials were caused by
> overlapping rings (two rings were going through inside each other) which
> genbox unfortunately somehow had generated. I hadn't though of that
> myself until a colleague made a remark that this might be the case. I
> used genconf to do another simulation box, where there is no chance of
> this happening, and my results turned out rather good.

Unlucky! Still, it underscores the old story that one needs to visualize 
results to check they follow one's preconceptions. In extremis, with 
ring structures, use genconf with a dummy space-filling atom that you 
edit out later.

Mark



More information about the gromacs.org_gmx-users mailing list