based approaches, which are in general intermediate in accuracy and computational expense between MM and full QM. While these approaches are currently too slow for routine molecular dynamics sampling in the condensed phase, their use for the rapid generation of reference data for parametrizing FFs is appealing, as demonstrated by the recent refinement of GAFF2 parameters against ANI2x torsional scans .