[gmx-users] dihedral potential structure

[Филипп Орехов]

Hello gromacs users,

I want to ask a simple thing, which nevertheless is not really clear to me,

but it's sometimes too important. For instance, I have in my molecule
torsinal

angle around a double bond: -C-C(H)=C(H)-C-. The dihedral potential contains

following four parts:

-C-C=C-C-, a U1 from topology

-H-C=C-C-, a U2=U3 from topology

-C-C=C-H-, a U3=U2 from topology

-H-C=C-H-, a U4 from topology

Attention please, here the question is:  how  does the final potential look
like?  Is it just only a sum

Utot=(U1+U2+U3+U4) or gromacs divides it on the number of elements of the
dihedral, I mean

Utot=0.25*(U1+U2+U3+U4)?

A lot of thanks in advance.
-- 
With best regards/С наилучшими пожеланиями,
Orekhov Ph.S./Орехов Ф.С.
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://maillist.sys.kth.se/pipermail/gromacs.org_gmx-users/attachments/20090404/fbc63c19/attachment.html>