From: AnkiReddy katha (kathaankireddy_at_gmail.com)
Date: Thu May 07 2015 - 02:12:03 CDT
Dear NAMD users
I am trying to use OPLS force field in NAMD.
vdwGeometricSigma yes
1-4 scaling 0.5
The above are changes we need to make for non-bonded interactions.
Bonded interactions (bonds term and angle term) needs no change as both
OPLS and CHARMM type used same expressions.
Dihedral term expression also is similar, so one can easily adapt OPLS
dihedral parameters into the CHARMM type.
There seems to be an issue with Improper style.
OPLS force field uses periodic improper dihedral, CHARMM type uses
harmonic improper dihedral
One idea I had was that we can have have impropers and propers under the single
head "dihedrals" in the parameter file, like OPLS style as shown below.
74 21 5 9 8 OPLS_IMPROPER 0.000000 0.000000 2.000000
0.000000 0.000000 0.000000 0.000000
75 15 11 10 5 OPLS_IMPROPER 0.000000 0.000000 8.000000
0.000000 0.000000 0.000000 0.000000
If we do that, then NAMD may not calculate "improper" energy term
explicitly, however it includes in the "dihedral".
Also it may consider 1-4 interactions between the atoms which are 1
and 4 atoms, which are actually not separated by three
bonds.
If some of you have used OPLS in NAMD, could you please clarify how
you handled OPLS_impropers.
Thank you very much
Best Regards
Anki
This archive was generated by hypermail 2.1.6 : Thu Dec 31 2015 - 23:21:52 CST