[gmx-users] LJ potential definitions?

Anton Feenstra feenstra at chem.vu.nl
Thu Jan 23 11:08:45 CET 2003


Hi All,


I am trying to figure out some strange desagreement between LJ potential
definitions.

 From the Gromacs manual (p50/51) I have:

	E(r) = 4? ( (?/r)^12 - (?/r)^6 )

but from Cornell et al. (& Kollman), JACS 117 (1995), eq. 1, I have:

	E(r) = A/(r^12) - B/(r^6)
where (from Table 13):	A = ? ?^12 and B = 2? ?^6
(they use R for r and R* for ?)
Rewriting this, gives the following:

	E(r) = ? ( (?/r)^12 - 2 (?/r)^6 )

So the question is, compared to the 'Gromacs' formula, where did the
factor 4 in 4? go, and second, where does the factor 2 for the sixth
power term come from?


-- 
Groetjes,

Anton
  ________ ___________________________________________________________
|        | Anton Feenstra                                            |
| .      | Dept. of Pharmacochemistry - Vrije Universiteit Amsterdam |
| |----  | De Boelelaan 1083 - 1081 HV Amsterdam - The Netherlands   |
| |----  | Tel: +31 20 44 47608 - Fax: +31 20 44 47610               |
| ' __   | Feenstra at chem.vu.nl - http://www.chem.vu.nl/afdelingen/FAR|
|  /  \  |-----------------------------------------------------------|
| (    ) | Dept. of Biophysical Chemistry - University of Groningen  |
|  \__/  | Nijenborgh 4 - 9747 AG Groningen - The Netherlands        |
|   __   | Tel +31 50 363 4327 - Fax +31 50 363 4800                 |
|  /  \  | K.A.Feenstra at chem.rug.nl - http://md.chem.rug.nl/~anton   |
| (    ) |-----------------------------------------------------------|
|  \__/  | "If You See Me Getting High, Knock Me Down"               |
|        | (Red Hot Chili Peppers)                                   |
|________|___________________________________________________________|




More information about the gromacs.org_gmx-users mailing list