----------------------------------------
> Date: Thu, 31 Jan 2008 15:29:50 +0100
> From: [EMAIL PROTECTED]
> To: [email protected]
> Subject: RE : FEP : separating components of dgdl
>
> 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 ?
>
No, those terms are always included.
> 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 ?
I have no clue.
But the easy way to find out is to take your system with only N->DUM
and do two runs of zero steps at two different lambda values and see which
energy term
in the log file has changed.
>
> Other question : is there some information available about the "best"
> settings of soft core parameters with amber99 ?
As long as you are working with bio-molecules, the optimal soft-core parameters
are the same for all force-fields.
Berk.
_________________________________________________________________
Express yourself instantly with MSN Messenger! Download today it's FREE!
http://messenger.msn.click-url.com/go/onm00200471ave/direct/01/_______________________________________________
gmx-users mailing list [email protected]
http://www.gromacs.org/mailman/listinfo/gmx-users
Please search the archive at http://www.gromacs.org/search before posting!
Please don't post (un)subscribe requests to the list. Use the
www interface or send it to [EMAIL PROTECTED]
Can't post? Read http://www.gromacs.org/mailing_lists/users.php
