The diversity of the optimization results (see figures 5 and 8) shows the multi-modal nature of the loss function's landscape. Almost every optimization run results in a different FFParams set, depending on the optimization setup. A major disadvantage of this local optimization algorithm is that it can get stuck in a local minimum of the optimization landscape, unable to discover other minima that might be representatives of more suitable parameter sets. A consequence is that the current optimization workflow is very dependent on the initial parameter set. This might be a small issue when an already well-performing force field is re-optimized, but can be challenging for newly investigated models of molecules.