Note: This post will be moved to Force Field / FF releases / Parsley minor release(s) once it is ready to be shared.
Fitting Data and Results
Fitting targets: 2nd generation training sets (link for the details of the training set generation scheme: http://doi.org/10.5281/zenodo.3777278)
...
Initial X2 | Final X2 | |
|---|---|---|
v1.2.0-RC1(https://openforcefield.atlassian.net/wiki/x/AQC4F ) | 3.619E+04 | 6.877E+03 |
v1.2.0-RC2 | 8.710E+03 | 6.843E+03 03 |
Benchmark Results
Benchmark data
...
| objective value (X2) |
|---|---|
initial guess | 29,469 |
v1.1.0 | 20,097 |
v1.2.0-preliminary (link: http://doi.org/10.5281/zenodo.3781313 ) | 16,939 |
v1.2.0-RC1 | 16,713 |
v1.2.0-RC2 | 16,910 |
(2) v1.2.0-RC1 vs. v1.2.0-RC2: Parameter comparison
Direct parameter comparison has found no significant differences in equilibrium bond lengths and equilibrium angles while showing notable differences in k values of some angle/ torsion parameters.
Angle terms with significant different final k values between RC1 and RC2:
a6([#1:1]-[*;r3:2]~;!@[*:3], k value in SMIRNOFF99Frosst: 100 kcal/mol/radian2)a3([*;r3:1]1~;@[*;r3:2]~;@[*;r3:3]1, k value in SMIRNOFF99Frosst: 200 kcal/mol/radian2)a15([#8X1:1]~[#6X3:2]~[#8:3], k value in SMIRNOFF99Frosst: 126 kcal/mol/radian2)
: Based solely on my intuition without no strong evidence, RC2 angle k values for the angle terms(doesn’t change much from 1.1.0 during the optimzation) seem physical; 400kcal/mol/radian2 for angle seems too large compare to other angle k values. Final gradients for the angle k values are also higher in RC1 (a6: 5.580e+00 , a15 : 4.750e+00) compared to the gradients in RC2 (a6: 1.599e+00 , a15 : 1.391e+00)
: bond/ angle coverage check? (Wed or later today)
...
Torsion terms with significant different final k values between RC1 and RC2
...
running now
...
:
t146,t147: 6 periodicitiest15([*:1]-[#6X4;r3:2]-@[#6X4;r3:3]-[*:4]),t16([#6X4;r3:1]-[#6X4;r3:2]-[#6X4;r3:3]-[*:4]): in-ring rotation
...