Dear David,
> > The reason can be in order of atoms in IMPROPER angles for NH2 and COO- groups.
>
> This will only really be fixed by a change in the way the impropers are
> defined, I think. Someone will have to take this on as a project; it will
> be easier once we implement the CHARMM force field equations (but don't hold
> your breath.)
It can be done also by reordering atom positions in topology file
without any code or energy formula change. I have written program in
python which reads original topology file and AMBER libraries for
amino acids and nucleic to create new topology file with correct order
of atoms for IMPROPER's. With new topology file any permutations of
equivalent atoms give the same results - potential is symmetric.
I can share this program, please let me know if somebody is interested.
With best regards,
Edyta Malolepsza
-----------------------------------------------------------------------
The AMBER Mail Reflector
To post, send mail to amber.scripps.edu
To unsubscribe, send "unsubscribe amber" to majordomo.scripps.edu
Received on Fri Apr 18 2008 - 21:12:19 PDT
