[gmx-users] Result difference between Gromacs and other pakages, is it a bugs??(oplsaa)

PeiQuan Chen gromacs at 163.com
Wed Jun 11 13:05:01 CEST 2003


Dear gmx-users:
>  Have you done an exact comparison of parameters?
>(share/top/ffoplaabon.itp)
  Yes. I also use the gmxdump -s topol.tpr to check the parametes.
The following is part of the file:

            functype[0]=LJ, c6= 2.03186437e-03, c12= 3.73510579e-06
            functype[1]=LJ, c6= 4.99228307e-04, c12= 3.34443996e-07
            functype[2]=LJ, c6= 4.99228307e-04, c12= 3.34443996e-07
            functype[3]=LJ, c6= 1.22660160e-04, c12= 2.99463281e-08
            functype[4]=BONDS, b0A= 1.09000e-01, cbA= 2.84702e+05, b0B= 1.09000e-01, cbB= 2.84702e+05
            functype[5]=ANGLES, thA= 1.07800e+02, ctA= 2.76329e+02, thB= 1.07800e+02, ctB= 2.76329e+02

         Bond:
            nr: 12
            multinr[division over processors]: 12
            iatoms:
               0 type=4 (BONDS) 0 1
               1 type=4 (BONDS) 0 2
               2 type=4 (BONDS) 0 3
               3 type=4 (BONDS) 0 4
  The parametes is just as the same I have used in tinker expect for the unit. And I have give the final 
conformation after em. the C-Cl bonds is not equal as I have got using tinker.
  The four C-Cl bond distances is 1.090000 1.081665 1.091833  1.091833, which tinker result is four 1.09000.
Do you think whether is a precision problem?


Sincerely yours
---- 
PeiQuan Chen
________________________________________________________________________
PeiQuan Chen 	Department of Chemistry, Nankai University
WeiJin Road 94, TianJin, 300071, P.R.China
phone:	(86)2223506728		Mobile Phone: (+86)13920675030
gromacs at 163.com	
++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++





More information about the gromacs.org_gmx-users mailing list