Results of different fitting strategies | JH โ Tried using (1) all torsion terms available, (2) just parsley torsion terms. JH โ Is there any downside to only doing ONE bespoke parameter on top of Parsley? Like, if there were several torsions running through a central bond, could we bespoke-fit just one of them? DM โ Trying to think of cases where things on different sides of a central bond ARENโT as rigid as aromatic rings, and might need to know about their flexibility. We wouldnโt want to flex certain parts of end groups. JH โ Iโm assuming that theyโre relatively rigid SB โ Could we make SMARTS where itโs the entire molecule SMILES, with the 2- and 3- atoms tagged explicitly, and any atom that is in the 1- or 4- positions would be written as a wildcard. DM โ Not sure that this would work, this might break something in the construction of SMARTS.
JH โ Found some cases where we do poorly relative to GAFF. Refitting
(?starting from GAFF?) could fix these. Where would be the best place to report these? DM โ Striketeam channel. Trevor could use this info. JH โ I use a higher energy limit before downweighting, might be important to mention as well. Will also try running with different priors.
JW โ If fitting multiple torsions would produce a certain energy profile, are we sure that a single torsion would reproduce the same energy profile? Or would it be necessary to include many to encode different phases? (General) โ This isnโt mathematically proven, but it seems likely based on our experiments that fitting a single torsion is โgood enoughโ DC โ Case like dibromoethane would trigger this sort of behavior. JH โ Right, I had to split that into multiple terms.
DM โ If we come up with code to generate a single torsion SMIRKS that โfits everythingโ for a particular rotatable bond, we should share that with JM.
|