I am trying to learn best practices for force matching ab initio molecular dynamics (AIMD) data with nonlinear least-squares minimization to generate interatomic potentials for use in classical MD.
I understand the concepts of selecting appropriate reference structures (for the expected classical MD behavior) and fitting to physically-justified functional forms. However, the art of selecting the right weighting factors, the right amount of reference structures, and judging if a set of potential parameters are fit properly appears to be highly open-ended problem with no single solution.
I was wondering if anyone knows of good resources (or has any advice) on gaining mastery and judgment over this process.