Post slides here and recording here
(slide 6) MG – One thing I’ve noticed is that the approach to the minimization strongly affects the final structure. Did you consider using restraints?
BS – Proton transfer cases made it through benchmark workflow -
(slide 17) AG – Seems like this could be the effect of angle terms being a little too stiff.
DM – CBayly fretted a lot about angle values in rings, to ensure we get the right stiffness.
AG – Are angle terms different for atoms in rings?
DM – I don’t recall. But I do remember that we wondered about different parameters for different sized rings.
AG – I’d hope that the FF could figure that out. If you assume a ring, the ring is enforces by connectivity+bond length.
DM –
AG – I’m thinking that the out of ring angles may be too stiff
LW – For the molecule shown on this slide, it wouldn’t get any ring-specific angle parameters.
DM – So maybe for flexible terms we should differentiate that.
PB – I think we have specific angle parameters for 3- and 5-membered rings.
(mol 00971)
(mol 00553)
00946
00066
00660
DM – I was expecting to see more sulfonamides
BS – These were to act as a showcase - We have lots of examples of sulfonamides as well. We also had targets that were all-orphan which might arguably be worse.
MG – It’d be good to do an analysis of which parameters appear in these bad structures.
DM – I can think of two broad approaches to this - 1) parameter splitting/chemical perception based refinement, and 2) Training/testing dataset generation. Those will feed in well to our FF creation/fitting pipeline.
BS – You’re pointing out that there are point problems that could be fixed - I’m pointing out that these metrics should be included in training.
MG – I know that we’d been talking about this quite a lot - Basically in things like torsion drives, how do we measure energy/structure deviation.
BS – I wonder how these metrics are factoring into training.
PB – When we use the torsion profile target, we we try to match to the minimized profile. JH pointer out that we could match….
DM – What about opt geo?
PB – We can imporve that if we include the internal coordinate differences, it kinda improves the
DM – I don’t think to ddEs at all
PB – Right, we fit to relative energies from torsion profiles, not optimized geometries. For optgeo we train on geometry but not energies.
AG – Is that an opportunity to add, or would that make it more complicated?
PB – We could add that.
AG – Would also need to make a database of test cases.
DM – We have lots of data already, but this may help us construct a good, well-defined test set.
AG – I’m not sure if a “well defined” set is possible - This is where the FF has to consider “diversity” - And you could use the entire benchmarking set for that.
PB – I think we could test the idea to fit to rel energies among optgeo.
BS – So one thing would be to fit to rel energies in matched conformers. And to
PB – Right.
AG – I think this can bring in a few things that are not represented in torsion profiles, like through-space interactions and ring conformations.
DM – Could we get the list of orphans as well? I’d want to check for parameter enrichment for those.
LW – I’d made a method that tries to figure out which energy errors come can be attributed to which parameters, so I’d be interested to try running on the chemistries that make lots of orphans.
DM – We generally try to avoid structures with a lot of steric congestion since those make it hard to figure out the sources of error.
AG – Could we do a sensitivity analysis to see which parameters are related to high ddEs - Like to find which situations have instances of parameter application pulling the parameter values in different directions.
PB + LW – We can work with BSwope to support this analysis.