2023-08-31 Force Field Release Meeting notes

 Date

Aug 31, 2023

 Participants

  • @Lily Wang

  • @David Mobley

  • @Chapin Cavender

  • @Alexandra McIsaac

  • @Pavan Behara

  • @Brent Westbrook

  • @Christopher Bayly

  • @Jeffry Setiadi

  • @Jeffrey Wagner

  • @Matt Thompson

  • @Trevor Gokey

  • @Anika Friedman

  • @Michael Shirts

Meeting recording: https://drive.google.com/file/d/1yStrjLU4AMy3y77F491-kTb6pHhZCIH-/view?usp=drive_link

 Discussion topics

Item

Presenter

Notes

Item

Presenter

Notes

Short GNN update

LW

 

  • CB – For nagl charges element limitations, people generally do selenomethionine by giving Se the same charges as S…. but Se ins’t covered by AM1BCC… so I guess it makes sense that nagl won’t handle it

  • CB – Re: Whisker plots - If AM1BCC charges would perform poorly on some unusual chemistries - then you’d expect nagl to also do poorly. So maybe it’d be good to do … (turns out to be misunderstanding, as the whisker plot was comparing to AM1BCC ESPs)

    • LW – Yeah, looking forward to training to better level of theory. But we’re starting with AM1BCC as a low-risk trial to see how it performs and how community takes it.

  • CB – Re limitations - Perchlorates and perbromates may be useful

  • DM – Reiterating that especially unusual chemistries aren’t something we should focus on.

    • CB – Agree.

  • JW – So this is now ready for the public API, but not to automatically be used in system creation?

    • LW – Yes. This is a release candidate, but it’s ready for the public API.

    • CB – Wouldn’t it be better to automatically switch over for large molecules?

    • JW + MS – It’ll be prominently featured in examples and docs, but it won’t be automatically used.

    • CB – It’d be good for there to be a setting where it switches regimes to make it automatically use nagl.

    • DM –

    • CC – One use case that might fail would be polymers with different lengths - The charge method could switch in the middle of a study.

    • JW – Could have …

    • (discussion)

    • CBy – People would want to have some flag where everything that could be handled using “old faithful” OE or AT would be, but big things would be handled by nagl.

  • MS – Would be good to document all this data+methods

    • LW – Yes, writing up docs

  •  

  • LW – Ideas for validation?

    • MS – TB will be doing solvation and mixture calcs of small molecules. Would that be helpful?

      • LW – Yes, that’d be great. I can help set those up.

      • MS – Timing/details of leave?

      • LW – Out sept 9 - early October

    • JW – I think that the FE calcs may be so noisy that nagl charges wouldn’t make a significant difference

    • CB – ESP comparison is the major thing to look at. That’s what people are gonna look at to decide when choosing whether to use it.

    • LW – I previously looked at solvation free energies on bad outliers - That showed that the GNN charges do lead to a significantly worse prediction. So I think that PL FE calcs would be a useful source of information.

    • DM – I think, if the FE calc results look comparable, this should be good to go.

    • CC – Would we need to do a valence refit? It’s possible that we’d see worse performance. Would be useful to run the geometry/ddE benchmarking.

    • CB – Yeah, things like TFDs and dimer energies would be handy to check.

    • LW – You’ve convinced me that dimer energies should be part of our standard benchmark set.

    • MT – LW, what are your leanings?

    • LW – I see lots of the flaws with NN charges. I’m happy with peptide benchmarks but some small mols are worse. But the ones where we do poorly on are fairly exotic and I don’t think they’re realistically in scope. So PL benchmarks would make me feel a lot better. And valence benchmarks are much lower hanging fruit. But I’m so close to it that outside perspectives are helpful.

    • CB – People are going to say that this is being provided as a substitute for AT/OE AM1BCC. You already have HFEs calculated, you could plot GNN charge HFEs vs. OE charge HFEs. You’d hope to see a straight diagonal line, and that would convince users.

    • LW – We have done this for molecules with different charges - We see off-diagonals. The average RMSE for this plot is <1 kJ/mol. In this analysis it does do slightly worse than AT.

    • CB – I’d love to see that in the docs - And highlighting the molecules that did poorly.

    •  

    •  

    •  

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 Action items

 Decisions