[gmx-users] OPLS force field issue
jalemkul at vt.edu
Thu Apr 30 19:33:49 CEST 2015
On 4/30/15 11:28 AM, Mary Less wrote:
> Dear Justin,
> Thank you for your reply.
> I am trying to write an opls unied atom forcefield for Tetrahydrofuran (THF). I
> thought the common way to write a forcefield
> is to consider the atomtypes in the atomtypes.atp and choose those types which
> are suitable for the problem.
> (Considering things like the chemical surroundings as well as the charge to be zero)
> If I consider opls_071 for example I found the description "CH2 (SP3) ALKANES"
> which describes the chemical behaviour for
> my problem very well, as I think.
> I choosed opls_110 because of the description "United-atom Ether CH2 (-O)"
> because of the
> fact that the CH2 is bonded to an oxygen atom there as well.
> So do I just have to choose other atomtypes to build the united atom THF forcefield?
> What points do I have to take into account when choosing suitable atomtypes?
> Nevertheless I am wondering, that there are in general no bonding, angle or
> dihedral parameters for C2 given in the ffbonded.itp.
> (For example C2 and C3 or C2 with any arbitrary other atomtype)
> Does this mean that the type C2 is never bonded to any other atomtype?
More likely it's just a set of parameters that was never ported over. The
OPLS-UA force field is quite old, likely outdated at this point, but perhaps
people are using it for something. The OPLS-AA is more modern and more
frequently used, so it is fully supported. If there are missing parameters,
either (1) those are the wrong types to use (sorry, can't help you there, and
the documentation is rather cryptic) or (2) it's just something missing in the
force field implementation that should be somewhere in the literature.
Justin A. Lemkul, Ph.D.
Ruth L. Kirschstein NRSA Postdoctoral Fellow
Department of Pharmaceutical Sciences
School of Pharmacy
Health Sciences Facility II, Room 629
University of Maryland, Baltimore
20 Penn St.
Baltimore, MD 21201
jalemkul at outerbanks.umaryland.edu | (410) 706-7441
More information about the gromacs.org_gmx-users