The BespokeFit workflow concludes with the optimization of the torsion parameters to the QM (or alternative) reference data using an interface with the ForceBalance package . ForceBalance iteratively optimizes the FF parameters related to the torsional potential, which in the case of the class 1 additive FF employed here, is described by the following truncated Fourier series: