openff-forcefields icon indicating copy to clipboard operation
openff-forcefields copied to clipboard

Distorted sulfonamide geometries in 2.1.0-rc.1

Open j-wags opened this issue 2 years ago • 3 comments

          @j-wags Hi Jeff Thank you for the help.

I have retested this with the new 2.1.0 but I don't think this has been fixed yet? Thank you. image

Originally posted by @xiki-tempula in https://github.com/openforcefield/openff-toolkit/issues/1522#issuecomment-1516634852

j-wags avatar Apr 20 '23 16:04 j-wags

@j-wags I can take a look.

pavankum avatar Apr 20 '23 17:04 pavankum

Thank you.

xiki-tempula avatar Apr 20 '23 18:04 xiki-tempula

It's worth seeing what 2.2.0 does with this molecule

mattwthompson avatar Apr 22 '24 14:04 mattwthompson

@mattwthompson using your code with openff-2.2.0.offxml gives me a much nicer minimized structure.

image

jaketanderson avatar Sep 07 '24 21:09 jaketanderson

The origin of the issue was that the k's for the sulfamide angles (a31 and a32) were optimizing to a very low value during force field optimization. We think this low k led to the angle term being overwhelmed by nonbonded interactions (just a hypothesis). We’ve temporarily fixed it in Sage 2.2 by very heavily regularizing (all) the parameters during optimization, preventing them from deviating too much from the initial modified seminario method value, but that’s not a permanent solution as the regularization did cause some minor issues with other parameters. We're looking into more lasting ways to fix it for future releases, but from the user side, this should be fixed now.

amcisaac avatar Sep 09 '24 20:09 amcisaac