[gmx-users] non-bonded parameters

nk yx nkyx at hotmail.com
Thu Nov 17 08:40:11 CET 2005


hi all,
    I am confused about the calculation of non-bonded LJ parameters in the 
ffG43a1nb.itp. According to the manual, the combination rule of C6 and C12 
is:
Cij(6)=(Cii(6)*Cjj(6))^0.5
Cij(12)=(Cii(12)*Cjj(12))^0.5
However,in the ffG43a1nb.itp file, the values seem not to be in agreement 
with the rule. For example,
O  C6=0.0022619536 C12=7.4149321e-07
OM C6=0.0022619536 C12=7.4149321e-07
OA C6=0.0022619536 C12=1.505529e-06

According to the rule, 
O OA   C6=0.0022619536 
       C12=(7.4149321e-07 * 1.505529e-06)^0.5 = 1.0565697e-06 
but the value is 1.380375e-06 in the ffG43a1nb.itp.
In addition, C12 value between O and OM, and tween O and OA should be the 
same, but in the file, C12 is 1.38380375e-06 and 2.258907e-06,respectively.
Could anybody tell me why? Thanks in advance.

Regards,
nkyx

_________________________________________________________________
与联机的朋友进行交流,请使用 MSN Messenger:  http://messenger.msn.com/cn  




More information about the gromacs.org_gmx-users mailing list