Versions Compared

Key

  • This line was added.
  • This line was removed.
  • Formatting was changed.

...

  • ['TIG0', 'TIG1a', 'TIG1b'] - General torsions

  • ['TIG1c', 'TIG1d', 'TIG2', 'TIG3', 'TIG4', 'TIG5a', 'TIG5b', 'TIG6', 'TIG7', 'TIG8'] - interpolated

View file
namefit4.offxml

Fit4.1: For each of the interpolated parameter a general torsion parameter is created where the central bond can be a single, aromatic or double bond (denoted by letters p,q,r at the end of parameter id). Due to lack of enough training data that match those patterns only a subset of those are trained and here are the parameters optimized.

...

  • ['TIG0', ‘TIG1a', ‘TIG1b',
    ‘TIG3p’, ‘TIG3r', ‘TIG4p', ‘TIG5ap’, ‘TIG5bp’, ‘TIG1cp’, ‘TIG6p’, ‘TIG7p’, ‘TIG8p’, ‘TIG2p’, 'TIG2r’, 'TIG1dp’ ] - General torsions

View file
namefit4.1_trained_params_only.offxml

So, fit4.1 has two parameters extra when compared to fit4, and all are general torsion parameters. The total number of torsion parameters in [openff-1.3.0, fit4, fit4.1] are [167, 170, 172] respectively. And the parameters that are replaced with interpolated parameters are [t43, t44, t45, t48, t69, t69a, t70d, t76, t77, t78] based on smarts patterns.

...

For clarity adding more plots with subset of data:

...

Simpler Fits

Simple_fit1

Replacing t43, t44, t45 with an interpolated torsion parameter for the smarts pattern: “[*:1]~[#6X3:2]~[#6X3:3]~[*:4]” and fitting this FF to 170 targets from 14 datasets (listed below under fit0) that have a dihedral that matches this pattern. The objective function value is compared to the zeroth iteration of various other FFs and fits:

FF

X2 (obj. fn value)

simple_fit1

165.65

openff_unconstrained-1.3.0

197.78

openff_unconstrained-1.2.0

213.94

fit4

171.60

fit4.1

183.29

Performance on Benchmarking set of molecules:

...

To check whether the BM set is a good dataset to test these interpolated fits checked the performance on 'OpenFF Primary Benchmark 1 Torsion Set' from QCA:

...