michellab / BioSimSpace_examples

GNU General Public License v2.0
0 stars 0 forks source link

GAFF2 and parmchk #1

Open lohedges opened 6 years ago

lohedges commented 6 years ago

Since we want to support both GAFF and GAFF2 we will need to use parmchk2 rather than parmchk. See here or here for details.

To use different force fields with antechamber

# GAFF (default)
antechamber -at gaff ...
# or
antechamber ...

# GAFF2
antechamber -at gaff2 ...

To use different force fields with parmchk2

# GAFF (default)
parmchk2 -s 1 ...
# or
parmchk2 ...

# GAFF2
parmchk2 -s 2 ...

Note that parmchk2 does not produce the same frcmod output as parmchk when run in GAFF mode, i.e. running parmchk or parmchk2 without the -s 2 flag and comparing the output. I've tested with several example Mol2 files that were generated by antechamber using GAFF and there are differences in some improper terms. (I assume that these are due to changes in the methodology/algorithms used within parmchk2?)

This means that we'll need to re-run all of our existing parmchk examples with parmchk2, regardless of which force field is used. It would probably be a good idea to try GAFF and GAFF2 for the same systems. Since GAFF2 is still in beta we might want to work out cases where it is prone to failure.