[gmx-users] Solvation Free energy Calculation

[David Mobley]

Paolo,

> As for the LJ decoupling,I cannot obtain a smooth curve and dG/dlambda
> has large deviations even if I use sc_alpha = 0.5 and sc_power = 1 as
> suggested by Shirts.

This doesn't typically give a smooth curve in the sense of
featureless, like the electrostatics curve. However, the particular
choice of parameters from Shirts typically leads to a smoother curve
than you can get with any other choice of parameters, hence it's the
way to go. I typically end up using a  lot more lambda values for the
LJ curve than the electrostatics for the same reason (16 versus 5), as
you obviously need to put lambda values closer together anytime
dG/dlambda changes rapidly, at least for TI.

Hope that helps.
David

> I use OPLS-AA forcefield
>
> I think that there is probably something that I am missing but i don't
> know what.
>
> Paolo
>
>
> _______________________________________________
> gmx-users mailing list    gmx-users at gromacs.org
> http://www.gromacs.org/mailman/listinfo/gmx-users
> Please don't post (un)subscribe requests to the list. Use the
> www interface or send it to gmx-users-request at gromacs.org.
> Can't post? Read http://www.gromacs.org/mailing_lists/users.php
>