[gmx-developers] Possible bug in grompp?

[David van der Spoel]

On Mon, 2005-04-18 at 16:43 -0600, Walter Ash wrote:
> I've been working on a new molecule topology and I've noticed what I 
> believe is a bug in grompp...  I am using the ffopls forcefield with 
> genpairs=yes, fudgeQQ and fudgeLJ = 0.5,  and combination rule 3.  This 
> is with the 3.2.1 distribution version.
> 
> In the molecular topology, under [ pairs ], it appears I have to specify 
> LJ14 parameters in C6 and C12 format despite the use of sigma and 
> epsilon everywhere else in this forcefield (comb rule 3).  If I use 
> sigma and epsilon in a [ pairs ] definition, my LJ14 energies are orders 
> of magnitude too high as grompp intereprets these parameters as C6 and 
> C12 values instead.  Putting C6 and C12 parameters here gives the 
> expected behavior.
> 
> This is contrary to what is stated in the manual and what I would expect 
> the behavior to be.
> 
It seems that this was fixed recently for pair_types but I'm not sure
whether it was fixed for [ pairs ] as well. You could try the latest
CVS.

> 
> _______________________________________________
> gmx-developers mailing list
> gmx-developers at gromacs.org
> http://www.gromacs.org/mailman/listinfo/gmx-developers
> Please don't post (un)subscribe requests to the list. Use the 
> www interface or send it to gmx-developers-request at gromacs.org.
-- 
David.
________________________________________________________________________
David van der Spoel, PhD, Assoc. Prof., Molecular Biophysics group,
Dept. of Cell and Molecular Biology, Uppsala University.
Husargatan 3, Box 596,          75124 Uppsala, Sweden
phone:  46 18 471 4205          fax: 46 18 511 755
spoel at xray.bmc.uu.se    spoel at gromacs.org   http://xray.bmc.uu.se/~spoel
++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++