[gmx-users] Gromos parameters for atom type OA

Samuli Ollila Samuli.Ollila at fkem1.lu.se
Fri Dec 9 16:12:57 CET 2011


Dear all,

I was comparing Gromos 45a3 and 53a6 force fields and I found something which I cannot understand. The atom type OA is described in  ffG45a3nb.itp as:
;name  at.num   mass      charge  ptype      c6           c12
            OA    8     0.000      0.000     A  0.0022619536  1.505529e-06
 and in ffG53a6nb.itp:
            OA    8     0.000      0.000     A  0.0022619536  1.505529e-06
which are same parameters.

However, according to Gromos developers the atom type OA has been changed between Gromos 45a3 and 53a6. The parameters given by Gromos developers for 45a3 are
#       IAC     TYPE    SQRT(C6)   SQRT(C12(1))  SQRT(C12(2)) SQRT(C12(3))
              3       OA      0.04756 1.1250E-3       1.227E-3        0.0
and for 53a6
              3      OA     0.04756     1.100E-3     1.227E-3     0.0

The actual numbers are different for parameters used in Gromacs since the way to write parameters is different.

My problem is that there is difference in OA type between 45a3 and 53a6 in parameters given by Gromos developers but not in the files distributed in Gromacs package. Can someone explain this?

BR,
Samuli Ollila
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://maillist.sys.kth.se/pipermail/gromacs.org_gmx-users/attachments/20111209/c4a61e05/attachment.html>


More information about the gromacs.org_gmx-users mailing list