Our goal is to find the set of atom-centered parameters that provides an optimal reproduction of the EP while at the same time being transferable across a wide range of conformations. We apply two strategies to address both requirements. First, the parameters are fitted to the EP of all conformers of the same molecule simultaneously. Second, we apply a harmonic restraint towards a set of reference parameters to avoid large changes that only marginally affect the quality of the fit (see Fig. ).