By restricting the relaxation process to one step we circumvent this hurdle as the nominal energies are not of interest in this case. The second challenge that emerges is the fact that the distribution of coordination numbers of a specific particle shifts dramatically during the optimization process, from on average very low to medium/high as a more and more compact structures is formed. Our final training set has to account for this in order to be applicable.