2025-04-09 FF fitting meeting

2025-04-09 FF fitting meeting

Participants

  • @Lily Wang

  • @Matt Thompson

  • @Barbara Morales

  • @Chapin Cavender

  • Chris Fennell: Oklahoma State U - Chemistry

  • @Jennifer Clark

 

Recording: https://us06web.zoom.us/rec/share/80QI43MvPyxmJRoCx6aZ9Eq_4SPebduB7MVAh8Ppbwq24PThIXm-F84tAZetKyg.6bddilcQhYhNK0bp

Passcode: .^D8&=cf

(skip a few min in for technical difficulties)

Discussion topics

Item

Presenter

Notes

Item

Presenter

Notes

vdW terms, ML water models

Chris Fennell

  • LW: The water model you talked about on the last slide: is that still in progress?

    • CF: Yes, I didn’t actually show the best 3-point model, that was a very recent (last week) development that involved using our existing hyper-surface ML model but changing a new factor, the bond length (maybe other things too).

  • CC: Have you considered changing the combining rules for the FF?

    • CF: We are using LB, I don’t like it for the size parameter, but the effect of that change has not been studied and would / could be a large change.

  • LW: Was that 50 or 15 ns?

    • CF: 50 ns to converge the dielectric constant of methanol, water needs ~100 ns.

  • LW: How much do values bounce around during in-serial parameterisation?

    • CF: It depends, you need to choose how to perturb it, is it a 10% change? It’ll diverge. We then make 2% changes and then linearly interpolate the target. Sometimes if you don’t run the simulation long enough the dielectric constant won’t converge. It usually takes 2 to 4 loops.

  • LW: Did you look at non-alcohols e.g. carboxylic acids?

    • CF: Not yet, we are interested in moving on to amines and carboxylate groups

  • LW: Have you looked at the hydrocarbon parameters? When we refit GAFF the hydrocarbon size decreased but that’s the opposite of what you’re seeing.

    • CF: Yes that is the opposite, I would recommend looking at the diffusivity, given I haven’t looked at hydration free energies. I’m sure your forcefield is great though.

  • JAC: It seems like the traditional GAFF FF trends toward experiment but does it become slower than experiment as the hydrocarbon chain continues to increase?

    • CF: looking at this now.

  • JAC: Have you considered testing the viscosity of these last 7 hydroxyl models? There is usually a trade-off with modeling diffusivity and viscosity.

    • CF: We are working on it, but not yet.

NAGL vdW refit

LW

POSTPONED

Cresset torsion refit experiment

LW

POSTPONED

Action items

Decisions