[gmx-users] How to convert C6&C12 from GROMOS54a7 to sigma&epsilon for OPLSAA

Justin Lemkul jalemkul at vt.edu
Fri Aug 18 19:03:53 CEST 2017



On 8/18/17 3:43 AM, Ming Tang wrote:
> Dear all,
>
> I want to calculate sigma&epsilon for OPLSAA.ff from C6&C12 developed based on GROMOS54a7.ff.
>
>  From the manual, I got the equations:
>
>
>
> C12 = Sigma^(6)*C6
> C6 = 4*sigma^(6)*epsilon
> However, I found sigma and epsilon are different in OPLSAA and AMBER, and that C6&C12 from GROMOS and sigma&epsilon from OPLSAA do not match these equations. Could anybody give me some guidance?

Mathematically speaking, you can certainly transform one to another, but 
you shouldn't do it.  GROMOS parameters aren't guaranteed to mix with 
any other force field, so convention wisdom is to not waste time trying 
to do it.  Force fields are self-consistent entities. Mixing and 
matching yields a physical model of unknown (likely junk) quality.

-Justin

-- 
==================================================

Justin A. Lemkul, Ph.D.
Assistant Professor
Virginia Tech Department of Biochemistry

303 Engel Hall
340 West Campus Dr.
Blacksburg, VA 24061

jalemkul at vt.edu | (540) 231-3129
http://www.biochem.vt.edu/people/faculty/JustinLemkul.html

==================================================



More information about the gromacs.org_gmx-users mailing list