[gmx-users] RE : FEP : separating components of dgdl

[BON Michael]

Hi,

You were right for the bonds, I forgot that I was using constraints = all-bonds.

Could the missing term be the soft core then ? Or is it included in VdW ?

Anyway, I noticed that this problem of hidden extra term happens when I modify something else than a hydrogen. Does it make sense ?
For example, I did N->DUM with all the parameters kept the same : every components of dVdl are 0 (VdW, Coulomb, angle, proper dihedral, Ryck.,  LJ 1-4, position restraint, disp. corr. ), but dVdl/dlambda is -6. Any idea ?

Other question : is there some information available about the "best" settings of soft core parameters with amber99 ?
-------------- next part --------------
A non-text attachment was scrubbed...
Name: winmail.dat
Type: application/ms-tnef
Size: 2821 bytes
Desc: not available
URL: <http://maillist.sys.kth.se/pipermail/gromacs.org_gmx-users/attachments/20080131/69b649a2/attachment.bin>