[gmx-users] Parameterization of OPLS-AA force field for MoS2 simulation in GROMACS.
bcworley at utexas.edu
Wed May 18 23:12:59 CEST 2016
I'm trying to modify the OPLS-AA force field parameters for Mo and S atoms
to simulate monolayer MoS2. I'm following the general procedure suggested
but they use existing OPLS-AA values for atomtypes, bondtypes, etc. To my
knowledge, there is no Mo atom in OPLS-AA standard force fields.
I've read the Varshney paper "MD simulations of molybdenum disulphide
(MoS2): Force-field parameterization and thermal transport behavior," and
would like to use the parameters listed there in my force field. However,
I'm having difficulty translating their parameters for set 9 into relevant
GROMACS .itp parameters.
For example, how would one translate bond coefficients K and r0 to
bondtypes values b0 and kb? Or am I even on the right track?
Any help on this matter would be greatly appreciated!
Dept. of Chemistry
The University of Texas at Austin
More information about the gromacs.org_gmx-users