--- Berk Hess <[EMAIL PROTECTED]> wrote:
> >Is there something I'm missing? A bug/feature in the code? A
> >limitation?
>
> If you are using PME, the mesh contribution is not excluded when
> using energygrp_excl, grompp should give you a warning about this.
> Adding all the pairs as exclusions solves this problem.
Thanks. Would that account for such a large difference, considering
it's a neutral molecule (and at lambda=1, with no partial charges)?
Well... since it's the derivative, and that means the free energy
change from lambda=0 to lambda=1... yes that could be it. Is this
mentioned in the manual? 'Cause I read about the [exclusions] section
and it gives the advice to use energygrp_excl, but doesn't warn about
this.
Would normal Ewald or reaction field work as (I) expected? I guess I
will keep using PME anyway, if it's as easy as adding the [exclusions].
> But why do you want to freeze the methonal?
> Things would be much easier with a free methanol.
At the moment this is only a test. I need the method to work for
arbitrary molecules for which there are no force field parameters, only
LJ. The atomic charges and intramolecular contributions I get them from
other QM/MM calculations. I also need to compare with other rigid-body
simulations. I'm not saying this is the best way to calculate free
energies, but it's the way I need now :)
Thanks a lot for the reply!
PS. The manual mentions using "no" for "comm_mode", but grompp prefers
"none" instead.
___________________________________________________________
All New Yahoo! Mail Tired of [EMAIL PROTECTED]@! come-ons? Let our SpamGuard
protect you. http://uk.docs.yahoo.com/nowyoucan.html
_______________________________________________
gmx-users mailing list gmx-users@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 [EMAIL PROTECTED]
Can't post? Read http://www.gromacs.org/mailing_lists/users.php
